Smectic-C tilt under shear in smectic-a elastomers

Size: px
Start display at page:

Download "Smectic-C tilt under shear in smectic-a elastomers"

Transcription

1 Smectic-C tilt uner shear in smectic-a elastomers Olaf Stenull an T. C. Lubensky Department of Physics an Astronomy, University of Pennsylvania, Philaelphia, Pennsylvania 19104, USA J. M. Aams an Mark Warner Cavenish Laboratory, JJ Thomson Avenue, Cambrige CB3 0HE, Unite Kingom Receive 27 February 2008; publishe 8 August 2008 Stenull an Lubensky Phys. Rev. E 76, have argue that shear strain an tilt of the irector relative to the layer normal are couple in smectic elastomers an that the imposition of one necessarily leas to the evelopment of the other. This means, in particular, that a smectic-a elastomer subjecte to a simple shear will evelop smectic-c-like tilt of the irector. Recently, Kramer an Finkelmann e-print arxiv: ; Phys. Rev. E 78, , performe shear experiments on smectic-a elastomers using two ifferent shear geometries. One of the experiments, which implements simple shear, prouces clear evience for the evelopment of smectic-c-like tilt. Here, we generalize a moel for smectic elastomers introuce by Aams an Warner Phys. Rev. E 71, an use it to stuy the magnitue of SmC-like tilt uner shear for the two geometries investigate by Kramer an Finkelmann. Using reasonable estimates of moel parameters, we estimate the tilt angle for both geometries, an we compare our estimates to the experimental results. The other shear geometry is problematic since it introuces aitional in-plane compressions in a sheetlike sample, thus inucing instabilities that we iscuss. DOI: /PhysRevE PACS number s : Va, v, Df I. INTRODUCTION Smectic elastomers 1 are rubbery materials with the orientational properties of smectic liqui crystals 2. They possess a planelike, lamellar moulation of ensity in one irection. In the smectic-a SmA phase, the Frank irector n escribing the average orientation of constituent mesogens is parallel to the normal k of the smectic layers, whereas in the smectic-c SmC phase, there is a nonzero tilt angle between n an k. Recently, there has been some controversy about whether shear strain an tilt of the irector relative to the layer normal are couple in smectic elastomers an whether the imposition of one necessarily leas to the evelopment of the other. Beautiful experiments by Nishikawa an Finkelmann 3, where a SmA elastomer was subjecte to extensional strain along the layer normal, foun a rastic ecrease in Young s moulus at a threshol strain of about 3% accompanie by a rotation of k through an angle that sets in at the same threshol. Interpreting their x-ray ata, the authors conclue that there was no SmC-like orer over the range of strains they probe, an that the reuction in Young s moulus stems from a partial breakown of smectic layering. Recently, Aams an Warner AW 4 evelope a moel for SmA elastomers that assumes that n an k are rigily locke such that =0. This moel prouces a stress-strain curve an a curve for the rotation angle of k in full agreement with the experimental curves but without neeing to invoke a breakown of smectic layering. Eviently, because of the assumption =0, the AW moel preicts no SmC-like orer. More recently, Stenull an Lubensky 5 argue that shear strain an SmC-like orer are couple an that the imposition of one inevitably leas to the evelopment of the other. They evelope a moel base on Lagrangian elasticity that preicts, as oes the AW moel, a stress-strain curve an a curve for in full agreement with the experiment. In contrast to the interpretation of Nishikawa an Finkelmann of their ata an to the central assumption of the AW moel, Ref. 5 foun that the tilt angle n of n is not ientical to that of the layer normal,, implying that there is SmC-like orer with nonzero above the threshol strain. However, estimates of an n base on reasonable assumptions guie by the available experimental ata of the layer normal-irector coupling turn out to have the same orer of magnitue. The upshot is that Ref. 5 preicts SmC-like tilt above the threshol strain but that the angle is small. It is entirely possible that is smaller than the resolution of the experiments by Nishikawa an Finkelmann. In this case, there is no isagreement between the preictions of Ref. 5 an the experimental ata by Nishikawa an Finkelmann. The arguments of Ref. 5 imply in particular that a SmA elastomer subjecte to a shear in the plane containing k will evelop SmC-like tilt of the irector. Very recently, Kramer an Finkelmann KF 6,7 performe corresponing shear experiments using two ifferent shear geometries. One geometry, which we refer to as tilt geometry, imposes a shear that is accompanie by an effective compression of the sample along the layer normal. In this geometry, the elastomer ruptures at an impose mechanical shear angle of 13 eg 7 or 14 eg 6, an up to these values of no SmC-like tilt is etecte within the accuracy of the experiments of about 1 eg, as was the case in the earlier stretching experiments of Nishikawa an Finkelmann. The other shear geometry, which we refer to as slier geometry, imposes simple shear an can be use to probe values of exceeing 20 eg. For this geometry, the KF x-ray ata provie clear evience for the emergence of SmC-like tilt. In this paper, we generalize the moel introuce by AW 4 an use it to stuy the magnitue of SmC-like tilt uner shear for the two geometries investigate by KF. Using reasonable estimates of moel parameters, we estimate the tilt /2008/78 2 / The American Physical Society

2 STENULL et al. angle, an we compare our estimates to the experimental results. The outline of the remainer of this paper is as follows. In Sec. II, we evelop our moel, which generalizes the original AW moel. In Sec. III, we escribe the two experimental setups that we consier. We iscuss the respective eformation tensors for these setups an calculate for both setups the SmC-tilt angle as a function of the impose mechanical tilt angle. We aress why the tilter apparatus is problematic for experiments where shear inuces irector tilt. In Sec. IV, we iscuss our finings an we make some concluing comments. There are two Appenixes. In Appenix A, we comment on analytical calculations of for small. In Appenix B, we briefly iscuss the Euler instability in the context of the experiments by KF. II. GENERALIZING THE ADAMS-WARNER MODEL In this section, we evise a theory for smectic elastomers base on the neoclassical approach, evelope originally by Warner an Terentjev an co-workers 1,8 for nematic elastomers, an the subsequent extension of the neoclassical approach by AW to smectics. The neoclassical approach generalizes the classical theory of rubber elasticity 9 to inclue the effects of orientational anisotropy on the ranom walks of constituent polymer links. It treats large strains with the same ease as they are treate in rubbers in the absence of orientational orer, it provies irect estimates of the magnitues of elastic energies, it is characterize by a small number of parameters, an it accounts easily for incompressibility. In the neoclassical approach, one formulates elastic energy ensities in terms of the Cauchy eformation tensor =, efine by ij = R i / x j, where R x is the target space vector that measures the position in the eforme meium of a mass point that was at position x in the uneforme reference meium. In this approach, a generic moel elastic energy ensity for smectic elastomers that allows for a relative tilt between the irector an the layer normal can be written in the form f = f trace + f layer + f tilt + f semi. 2.1 Here an in the following, incompressibility of the material is assume, i.e., the eformation tensor is subject to the constraint et = =1. f trace is the usual trace formula of the neoclassical moel with the shear moulus, f trace = 1 2 Tr = = 0 = T = 1, 2.2 where = 0 = = + r 1 n 0 n 0, with n 0 escribing the uniaxial irection before eformation, is the so-calle shape tensor escribing the istribution of conformations of polymeric chains before eformation, an = 1 = = 1 r 1 nn is the inverse shape tensor after eformation. = enotes the unit matrix, an r enotes the anisotropy ratio of the uniaxial SmA state. The contribution f layer = 2 B cos 2.3 escribes changes in the spacing of smectic layers with layer normal k = = T k 0 = T k 0, 2.4 where k 0 =n 0 is the layer normal before eformation an B is the layer compression moulus an are, respectively, the layer spacing before an after eformation, which are relate via 1 = 0 = T k The tilt energy ensity, f tilt = 1 2 a t sin 2, 2.6 incorporates into the moel the preference for the irector to be parallel to the layer normal in the SmA phase. To stuy the SmC phase, we woul have to inclue a term proportional to sin 4, which, however, is inconsequential for our current purposes. Without the contribution f semi, the moel elastic energy ensity 2.1 is invariant with respect to simultaneous rotations of the smectic layers, the nematic irector, an = in the target space. To break this unphysical invariance, we inclue the semisoft term 1 f semi = 1 2 Tr = n 0n 0 = T nn =, 2.7 where is a imensionless parameter. The relation of the moel presente here to the AW moel is the following: AW assume that the layer normal an the irector are rigily locke such that the angle is constraine to zero. Moreover, the semisoft term f semi is absent in the AW moel. Essentially, we retrieve the AW moel from Eq. 2.1 by setting =0 or equivalently a t an =0. As mentione above, one of the virtues of the neoclassical approach is that it involves only a few parameters, an for most of these there exist experimental estimates, which we will now review briefly. The shear moulus for rubbery materials is typically of the orer of Pa. A typical value for the smectic layer compression moulus in smectic elastomers, as observe in experiments with small strains along the layer normal, is B 10 7 Pa, which is greater than the values in liqui smectics. In previous experiments by the Freiburg group, the anisotropy ratio was approximately r 1.1 3, an we aopt this value for our arguments here. The value of a t can be estimate from experiments by Brehmer, Zentel, Gieselmann, Germer, an Zungenmaier 11 on smectic elastomers an by Archer an Dierking 12 on liqui smectics. The former experiment inicates that a t is of the orer of 10 5 Pa at room temperature, an the latter experiment prouces a room-temperature value of the orer of 10 6 Pa. We are not aware of any experimental ata that allow us to estimate for smectic elastomers reliably. For nematic

3 SMECTIC-C TILT UNDER SHEAR IN SMECTIC-A φ Sliing Clamp Sample Fixe Clamp FIG. 1. Color online Sketch of a slier apparatus where the upper clamp slies on a horizontal bar such that the extension of the sample in the z irection remains fixe. elastomers, it has been estimate from the Freericks effect 13 an from the magnitue of the threshol to irector rotation in response to stretches applie perpenicular to the original irector 14 that 0.06 or 0.1, respectively. For our arguments here, we aopt the latter value acknowleging that may be consierably larger in smectics than in nematics. As we will see in the following, our finings o not epen sensitively on this assumption. For the eformations that we consier, appears only in the combination =a t +, that is, semisoftness simply as to the irector-layer normal coupling. Uncertainty in estimates for is expecte to stem mainly from the sprea in estimates for a t. We account for this sprea by iscussing several values of or, more precisely, for several values of /. III. EXPERIMENTAL SHEAR GEOMETRIES AND TILT ANGLES In this section, we apply the moel efine in Sec. II to stuy the behavior of SmA elastomers in shear experiments. We consier two experimental setups. In the first setup, which is perhaps the first that comes to min from a physicist s viewpoint, a simple shear is applie, i.e., a shear strain in which the externally impose isplacements all lie in a single irection. Figure 1 shows a sketch of such an experiment, an we call the apparatus sketche in it the slier apparatus. The secon setup, which is epicte in Fig. 2, is one in which opposing surfaces of the sample remain essentially parallel but in which the height of the sample its extension in the z irection ecreases upon shearing. Hence, the applie shear is, strictly speaking, not just simple shear. In the following, we will refer to this apparatus as a tilt apparatus. φ FIG. 2. Color online One of the two tilt apparatuses use by KF. The photo has been taken from Ref. 6. z z y y x x A. Slier apparatus To make our arguments as concrete an simple as possible, we now choose a specific coorinate system. That is, we choose our z irection along the uniaxial irection n 0 of the unsheare samples, n 0 = 0,0,1, an we choose our x irection as the irection along which the sample is clampe, which is perpenicular to n 0 ; cf. Fig. 1. With these coorinates, the eformation tensor for both the slier apparatus an the tilt apparatus is of the form = xx 0 xz = 0 yy zz. 3.1 As can be easily checke, for this type of eformation, the layer normal lies along the z irection for the unsheare an for the sheare samples. Thus, the tilt angle between n an k is ientical to the tilt angle between n an the z axis, an we can parametrize the irector as n = sin,0,cos. 3.2 It is worth noting that a eformation tensor of the form shown in Eq. 3.1 leas to a particularly simple expression for the semisoft contribution to f, f semi = 1/2 2 xx sin 2, which combines with the tilt energy ensity to f tilt + f semi = 1/2 xx sin 2. Thus, as inicate above, our moel epens for the experimental geometries uner consieration on an a t through a single effective parameter, viz. xx =a t + 2 xx. From the photos of the experimental samples provie in Ref. 6, see Figs. 1 an 2, it appears as if xx oes not eviate significantly from 1. In the following, we set xx =1 for simplicity 15. In this event, xx = a t +. In the slier apparatus, the extent of the sample in the z irection is fixe, an hence zz =1. The incompressibility constraint et = =1 thus manates that yy =1. The remaining nonzero component of the eformation tensor is entirely etermine by the externally impose shear, xz =tan, where is the mechanical tilt angle of the sample; see Fig. 1. The only remaining egree of freeom in the problem is, therefore, the angle. To calculate as a function of, we insert the just iscusse eformation tensor into f an then minimize f over holing fixe. For small, this can be one analytically by expaning f to harmonic orer in an by then solving the resulting linear equation of state for an equilibrium value of. This type of analysis is presente in the Appenix. To provie reliable preictions for larger shears, one has to refrain from expaning in an use numerical methos instea. To this en, we minimize f numerically assuming, base on what we iscusse at the en of Sec. II, that B/ =10 an / 0.1,1,5,20. For this minimization, we use MATHEMATICA s FinMinimum routine. Figure 3 shows the resulting curves for as a function of. The ashe horizontal line in Fig. 3 is a guie to the eye; it correspons to the angle resolution in the KF experiments, which was about 1 egree 16. For the slier geometry, the KF ata prouce clear evience for the evelopment of SmC-like tilt uner simple

4 STENULL et al. Θ (eg) (i) Sliing (ii) (iii) (iv) L Lsinχ Lcosχ φ χ Lcosχ-2 FIG. 4. Schematic of shear an compression arising in a nonieal tilter φ (eg) FIG. 3. Color online The tilt angle between the layer normal an the irector as a function of the mechanical tilt angle in the slier apparatus for i / =0.1, ii / =1, iii / =5, an iv / =20. The ashe line correspons to the resolution for in the experiments by KF. shear, an our theoretical estimates agree well with the experimental ata. However, given that thus far only two ata points are available for 0, it cannot be juge reliably whether our theoretical curve coul be fitte to an experimental curve with more ata points. It is encouraging, though, that our curve for / =0.1 agrees with the available experimental points within their errors. B. Tilt apparatus Now we turn to the tilt geometry. The essential ifference between the slier an the tilt geometry is that in the former, the height of the sample L z remains constant, L z =L z 0 times the number of smectic layers, whereas in an ieal tilter, L z is not constant but rather ecreases as increases, L z =L z cos. This ifference has far-reaching consequences. In the slier, the sample height remains larger than the moifie natural height of the sample create by the irector tilt, L z cos. In other wors, the sample is uner effective tension. In the tilter, on the other han, the sampleheight L z cos can be smaller than L z cos, i.e., the sample can be uner effective compression. For a sample uner effective zz compression, one has to worry, experimentally an theoretically, about all sorts of complications. Most notable are perhaps buckling an wrinkling. The theory of buckling of elastic sheets is well establishe 17, an we now briefly comment on buckling in the context of the tilt geometry. As mentione above, a sample clampe into the tilt apparatus can evelop a buckling instability if the height of the sample impose by tilt, L z cos, is smaller than the natural height of the sample create by the irector tilt, L z cos. As we will see below, our results for as a function of imply that cos cos, i.e., buckling is possible. Alternatively, this can be seen by calculating the engineering stress eng zz = f / zz, which can be one using the numerical approach outline above. eng zz is positive for 0 in the slier apparatus, whereas it turns out to be negative for 0 in the tilt apparatus. Thus, the sample is effectively uner tension in the slier apparatus, whereas it is effectively uner compression in the tilt apparatus, opening the possibility for buckling in the y irection in the latter. The angle c at which buckling sets in is expecte to be comparable to that for the well-known Euler-Strut instability 17, c Euler = arccos 1 2 L y 3L z 2, 3.3 where L y is the thickness of the sample in the y irection, an where clampe rather than hinge bounary conitions are assume. A brief erivation of Eq. 3.3 is given in Appenix B. In the experiments of KF, L z =5.0 mm an L y =0.45 mm 16, which leas to Euler 15 eg. The experimental samples buckle immeiately before they rupture 16 at =13 or 14 eg, which is very close to our estimate for c. The observation that buckling occurs immeiately before rupturing might suggest that the former actually triggers the latter. A etaile analysis of sample wrinkling in the tilt geometry is beyon the scope of the present paper. However, a few comments about wrinkling are in orer. The wavelength of wrinkling is expecte to be much shorter than that of buckling. Thus, wrinkling can be harer to etect by visual inspection of an experimental sample than buckling, an there is a risk that it remains unnotice. The main problem is, however, that wrinkling can act to effectively reuce the mechanical shear in the sample. When esigning a tilt experiment, one thus has to be very careful to avoi wrinkling. If not, the effects of wrinkling can bias the ata, an there is the risk of unerestimating the SmC-tilt significantly. Our calculation of the SmC-tilt will be for an ieal tilter. Figure 4 shows the tilter of Fig. 2 schematically in a nontilte an a tilte configuration to make clear that shear an compression are complex in the nonieal case because the frame axles allowing angle change are offset from the corners of the sample. Using elementary trigonometry, one can euce for the shears an the angle L xz = sin, L 2 zz = L cos 2, 3.4 L 2 tan tan = 1 2/L sec. 3.5 In particular, the relation between the shear angle an the eformation components zz an xz is not simple. Returning to an ieal tilter, =0,, we calculate the SmC-tilt as a function of the applie mechanical shear, suppressing buckling an wrinkling. Then, the eformation tensor is of the same form 3.1 as for the slier apparatus, but with the essential ifference that here xz =sin an zz =cos. As we i for the slier apparatus, we assume xx =1. The incompressibility constraint then implies that yy

5 SMECTIC-C TILT UNDER SHEAR IN SMECTIC-A Θ (eg) (i) Tilting (ii) φ (eg) FIG. 5. Color online The tilt angle between the layer normal an the irector as a function of the mechanical tilt angle in the tilt apparatus for i / =0.1, ii / =1, iii / =5, an iv / =20. The ashe line correspons to the resolution for in the experiments by KF. =1/cos, leaving the angle as the only egree of freeom. We calculate as a function of in exactly the same way as above, i.e., we minimize f numerically assuming the values of the moel parameters iscusse in Sec. II. The resulting curves are shown in Fig. 5. The ip in Fig. 5 occurs because of a competition between, which woul prefer =0, an B, which woul prefer =. If B, then it is possible to get curves where stays close to 0. The ashe horizontal line in Fig. 5 inicates the angle resolution of about 1 egree of the KF experiments 16. As mentione in Sec. III A, the ata points for the slier are compatible with / 0.1. The elastomers use by KF in the slier an the tilter are ientical, an, therefore, curve i of Fig. 5 shoul escribe the SmC-tilt if the tilter use by KF were ieal or nearly so. Note from Fig. 5, however, that this implies that the SmC-tilt for an ieal tilter at =13or14eg shoul be significantly larger than the experimental resolution, which is incompatible with the finings of Refs. 6,7. Since any effective compression of the sample along the layer normal, as in a tilter, promotes rather than hampers irector rotation, we obtain a SmC-tilt at a given in a tilter if it occurs at the same angle in the slier. It is likely that the failure of KF to observe irector rotation in their tilter is ue to mechanical instability that effectively reuces the mechanical shear of the sample an thus leas to a systematic suppression of the SmC-tilt in the tilt geometry. IV. DISCUSSION AND CONCLUDING REMARKS In summary, we have evelope a neoclassical moel for smectic elastomers, an we use this moel to stuy SmA elastomers uner shear in the plane containing the irector an the layer normal. In particular, we investigate the tilt angle between the layer normal an the irector for two ifferent experimental setups as a function of the mechanical tilt angle measuring the impose shear. Our moel buils upon the neoclassical moel for smectic elastomers by AW. In their original work, AW chose not to consier the possibility of SmC orer in their theory: they force the layer normal an the nematic irector to be parallel. In our present theory, the funamental tenet is ifferent in (iii) (iv) that the nematic irector n an the eformation tensor = are inepenent quantities that must be allowe to seek their equilibrium in the presence of impose strains or stresses. The layer normal k is etermine entirely by =. Thus, n an k are inepenent variables that rotate to minimize the free energy they are not locke together. In the present paper, we have focuse for simplicity on iealize monoomain samples. Bounary conitions in real experiments prevent this simple scenario an prouce a microstructure structure of the type that AW iscuss in Ref. 4 base on their original moel that prohibits SmC orering. Our current theory, which amits SmC orering, preicts the same type of microstructure; stretching will prouce a polyomain layer structure just as the AW theory oes. In particular, the existence of SmC orering oes not contraict the appearance of optical clouiness. To avoi unue repetition, we refrain from iscussing this here in more etail an refer the reaer to Ref. 4. The main fining of our present work is that the tilt angle between n an k is nonzero if a shear in the plane containing k is impose. This angle epens on material parameters, the experimental setup, an the magnitue of the impose shear. Figures 3 an 5 epict our results for the slier apparatus an the tilt apparatus use by KF. As expecte, in both setups the angle ecreases as the value of the parameter / increases, at fixe. This is because the irector is becoming more strongly anchore to the layer normal. If one keeps / fixe an varies B/ instea, then for B/, approaches because in this limit f layer locks to =. AsB/ is ecrease, then ecreases at fixe, because the f tilt term starts to compete with the f layer term. In the opposite limit of very small B/, is locke to =0 by virtue of the tilt an semisoft contributions to f. We refrain to show the corresponing curves, which are similar to those epicte in Figs. 3 an 5, to save space. The main ifference between our results for the two setups is that for B an large, say 40 egrees or so, is roughly one orer of magnitue smaller than in the slier apparatus, whereas it is of the orer of in the tilt apparatus. In closing, we woul like to stress once more that the experiment using the slier geometry prouces clear evience of the evelopment of SmC-like tilt uner shear as preicte originally in Ref. 5. Our theory evelope here allowe us to unerstan in some etail how the magnitue of this tilt epens on the shear geometry an on moel an material parameters, such as, which can be reasonably estimate from nematic analogues. Given these estimates, the resulting estimates for the tilt angle vary over a range from being small enough to being essentially unobservable with the x-ray equipment use by the Freiburg group to exceeing the experimental resolution. While the latter is consistent with the shear experiments in the slier geometry, the former is consistent with the stretching experiments by Nishikawa an Finkelmann in which no SmC-like tilt was etecte. We believe that the reason for the iscrepancy between the rotations in the slier an tilter apparatus lies in a possible wrinkling of the samples in the tilt apparatus use by KF. Mechanical instability leas to an effective reuction of the mechanical shear of the sample causing a systematic unerestimation of the SmC-like tilt in the tilt geometry

6 STENULL et al. Another experimental approach that coul be use, at least in principle, to investigate is to measure changes in the layer spacing, as was one by Nishikawa an Finkelmann 3. It shoul be note, however, that in this approach one measures cos rather than irectly. For small values of SmC tilt, say =1 egree, cos , which is practically inistinguishable from the cos =1 pertaining to SmA orer. For a refine comparison between experiment an theory, it woul be useful to critically evaluate an potentially improve the esign of the tilt apparatuses regaring samplewrinkling. Also, it woul be interesting to perform shear an stretching experiments with an angle resolution better than 1 eg. We hope that our work stimulates interest in such experiments. ACKNOWLEDGMENTS This work was supporte in part by the National Science Founation uner Grants No. DMR an No. DMR NSF MRSEC T.C.L.. We are grateful to D. Kramer an H. Finkelmann for iscussing with us some of the etails of their experiments. APPENDIX A: ANALYTICAL CONSIDERATION FOR SMALL ANGLES For arbitrary SmC-tilt angle, the total elastic energy ensity f is, for the eformations iscusse in Sec. III, a fairly complicate conglomerate of trigonometric functions an powers thereof of an. Thus, we resorte in Sec. III to a numerical approach to etermine the equilibrium value of over a wie range of the impose mechanical tilt angle. Here, we focus on the regime of where is small, such that it is justifie to expan f to harmonic orer in. In this case, the equations of state are linear in, of course, an are therefore reaily solve. For the slier apparatus, we obtain r 1 r tan = r + 1 r 1 r + r tan 2, an for the tilt apparatus we get A1 2 1 r 1 r sin 2 = r rb sin r 1 r cos 2. A2 Equations A1 an A2 imply that the SmC-tilt in both setups is ientical for small, r 1 r = r + 1 r 2 + O 3. A3 This is consistent with Figs. 3 an 5, where the initial slopes are ientical for ientical values of though this is perhaps somewhat har to see because the orinate scales are ifferent in the two figures. An interesting relate question, which we have not aresse in the main text because there are currently no experimental ata available, is that of the stress that is cause by the impose shear strains. From the above, it is straightforwar to calculate this stress for small. Inserting Eq. A3 into the aforementione harmonic in total elastic energy ensity provies us with an effective f in terms of. From this effective f, we reaily extract the engineering or first Piola-Kirchhoff stress as eng xz = f xz = r 2 r + 1 r 2 + O 3. A4 Equation A4 highlights a problem that arises when the tilt an the semisoft contributions f tilt an f semi to the total elastic energy ensity are missing, an one allows n an = to be inepenent quantities. Setting =0 in Eq. A4 leas to zero stress for nonzero, i.e., this truncate moel preicts soft elasticity. This soft elasticity, however, is not compatible with SmA elastomers crosslinke in the SmA phase, such as the experimental samples of Refs. 3,6, where the anisotropy irection is permanently frozen into the system. APPENDIX B: EULER INSTABILITY IN THE TILT APPARATUS In this appenix, we give a brief erivation of Eq We are intereste primarily in a rough estimate for c, an therefore, for simplicity, we assume that we can ignore the effects of smectic layering. In the following, we employ, for convenience, the Lagrangian formulation of elasticity theory. In the Lagrangian formulation, the elastic energy ensity of a thin elastomeric film with thickness L y an height L z that is compresse along the z irection can be written as f film = zu y Y 2u 2 zz, B1 where the choice of coorinates is the same as epicte in Fig. 2. u y is the y component of the elastic isplacement u =R x, an u 2 zz is the zz component of the Cauchy-Saint- Venant strain tensor, u ij = 1 2 iu j + j u i + i u k j u k. an Y 2 are, respectively, the bening moulus an Young s moulus of the film, which are given in the incompressible limit by = /3 L 3 y an Y 2 =3 L y 17. At leaing orer, z u z = L z /L z, where we have use that the height change L z is negative when the sample is effectively compresse as in the tilt apparatus. This leas to u zz = L z /L z zu y 2, B2 when we concentrate on the parts of u zz that are most important with respect to buckling in the y irection. Next, we substitute the strain B2 into Eq. B1 an switch to Fourier space. To leaing orer in the elastic isplacement, this prouces f film = 2 q 1 2 z Y 2 L z L z z q 2 ũ y q ũ y q, B3 where q is the wave vector conjugate to x, ũ y is the Fourier transform of u y, an so on. Equation B3 makes it transparent that buckling occurs for

7 SMECTIC-C TILT UNDER SHEAR IN SMECTIC-A L z L z = q 2 Y z = L 2 y 2 9 q 2 z. B4 The smallest value of q z is etermine by the specifics of the bounary conitions. From Fig. 2 it appears as if the sample of KF is clampe such that it prefers to stay parallel to the clamps in their immeiate vicinity. In this case, the smallest value of q z is q z =2 /L z. For hinge bounary conitions, in comparison, its smallest value woul be q z = /L z. Using the former value an exploiting that, approximately, L z /L z =1 cos, we obtain the estimate Eq. 3.3 for the onset of the Euler instability in the tilt apparatus. 1 For a review on liqui crystal elastomers, see M. Warner an E. M. Terentjev, Liqui Crystal Elastomers Clarenon, Oxfor, For a review, see, P. G. e Gennes an J. Prost, The Physics of Liqui Crystals Clarenon, Oxfor, 1993 ; S. Chanrasekhar, Liqui Crystals Cambrige University Press, Cambrige, E. Nishikawa an H. Finkelmann, Macromol. Chem. Phys. 200, J. M. Aams an M. Warner, Phys. Rev. E 71, O. Stenull an T. C. Lubensky, Phys. Rev. E 76, D. Kramer an H. Finkelmann, e-print arxiv: D. Kramer an H. Finkelmann, preceing paper, Phys. Rev. E 78, P. Blaon, E. Terentjev, an M. Warner, J. Phys. II 4, L. R. G. Treloar, The Physics of Rubber Elasticity Clarenon, Oxfor, Note that the layer compression term formulate by AW is of the form / 0 1 2, whereas setting =0 in Eq. 2.3 prouces / This form is chosen for consistency with the Chen-Lubensky moel of smectics. This term can be factorize into / / / for 0. Thus the ifference is merely a factor of 4 in the efinition of the smectic moulus B in the two representations. 11 M. Brehmer, R. Zentel, F. Gieselmann, R. Germer, an P. Zungenmaier, Liq. Cryst. 21, P. Archer an I. Dierking, Phys. Rev. E 72, E. M. Terentjev, M. Warner, R. B. Meyer, an J. Yamamoto, Phys. Rev. E 60, H. Finkelmann, I. Kunler, E. M. Terentjev, an M. Warner, J. Phys. II 7, We checke the valiity of this assumption by allowing xx to relax to its equilibrium value, i.e., we repeate our analysis using the xx that minimizes f rather than xx =1. The results obtaine by using this more involve approach o not iffer significantly from those obtaine using our simplifie approach with xx =1. 16 D. Kramer an H. Finkelmann private communication. 17 L. D. Lanau an E. M. Lifshitz, Theory of Elasticity, 3r e. Pergamon, New York,

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects

3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Homework 7 Due 18 November at 6:00 pm

Homework 7 Due 18 November at 6:00 pm Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Generalized hydrodynamic equations for nematic liquid crystals

Generalized hydrodynamic equations for nematic liquid crystals PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalize hyroynamic equations for nematic liqui crystals Tiezheng Qian 1,2 an Ping Sheng 1 1 Department of Physics, The Hong Kong University of Science

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Entanglement is not very useful for estimating multiple phases

Entanglement is not very useful for estimating multiple phases PHYSICAL REVIEW A 70, 032310 (2004) Entanglement is not very useful for estimating multiple phases Manuel A. Ballester* Department of Mathematics, University of Utrecht, Box 80010, 3508 TA Utrecht, The

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 4. Electrostatics of Macroscopic Media Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Model for Dopant and Impurity Segregation During Vapor Phase Growth

Model for Dopant and Impurity Segregation During Vapor Phase Growth Mat. Res. Soc. Symp. Proc. Vol. 648, P3.11.1-7 2001 Materials Research Society Moel for Dopant an Impurity Segregation During Vapor Phase Growth Craig B. Arnol an Michael J. Aziz Division of Engineering

More information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role of parameters in the stochastic dynamics of a stick-slip oscillator Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

12.11 Laplace s Equation in Cylindrical and

12.11 Laplace s Equation in Cylindrical and SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

G j dq i + G j. q i. = a jt. and

G j dq i + G j. q i. = a jt. and Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints

OF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints EUROSTEEL 2, August 3 September 2, 2, Buapest, Hungary A NUMERICAL EVALUATION OF CHS T JOINTS UNDER AXIAL LOADS Raphael S. a Silva a, Luciano R. O. e Lima b, Pero C. G. a S. Vellasco b, José G. S. a Silva

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Math 1B, lecture 8: Integration by parts

Math 1B, lecture 8: Integration by parts Math B, lecture 8: Integration by parts Nathan Pflueger 23 September 2 Introuction Integration by parts, similarly to integration by substitution, reverses a well-known technique of ifferentiation an explores

More information

Optimization of Geometries by Energy Minimization

Optimization of Geometries by Energy Minimization Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.

More information

Math 1271 Solutions for Fall 2005 Final Exam

Math 1271 Solutions for Fall 2005 Final Exam Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek

More information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear Dielectric Response of Periodic Composite Materials onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics

More information

arxiv:physics/ v4 [physics.class-ph] 9 Jul 1999

arxiv:physics/ v4 [physics.class-ph] 9 Jul 1999 AIAA-99-2144 PROPULSION THROUGH ELECTROMAGNETIC SELF-SUSTAINED ACCELERATION arxiv:physics/9906059v4 [physics.class-ph] 9 Jul 1999 Abstract As is known the repulsion of the volume elements of an uniformly

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Integration Review. May 11, 2013

Integration Review. May 11, 2013 Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In

More information

Code_Aster. Detection of the singularities and calculation of a map of size of elements

Code_Aster. Detection of the singularities and calculation of a map of size of elements Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & Transducers 2015 by IFSA Publishing, S. L. Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Code_Aster. Detection of the singularities and computation of a card of size of elements

Code_Aster. Detection of the singularities and computation of a card of size of elements Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size

More information

Asymptotics of a Small Liquid Drop on a Cone and Plate Rheometer

Asymptotics of a Small Liquid Drop on a Cone and Plate Rheometer Asymptotics of a Small Liqui Drop on a Cone an Plate Rheometer Vincent Cregan, Stephen B.G. O Brien, an Sean McKee Abstract A cone an a plate rheometer is a laboratory apparatus use to measure the viscosity

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Ductility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2

Ductility and Failure Modes of Single Reinforced Concrete Columns. Hiromichi Yoshikawa 1 and Toshiaki Miyagi 2 Ductility an Failure Moes of Single Reinforce Concrete Columns Hiromichi Yoshikawa 1 an Toshiaki Miyagi 2 Key Wors: seismic capacity esign, reinforce concrete column, failure moes, eformational uctility,

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

Magnetic helicity evolution in a periodic domain with imposed field

Magnetic helicity evolution in a periodic domain with imposed field PHYSICAL REVIEW E 69, 056407 (2004) Magnetic helicity evolution in a perioic omain with impose fiel Axel Branenburg* Norita, Blegamsvej 17, DK-2100 Copenhagen Ø, Denmark William H. Matthaeus University

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS SEQUENTIAL CIRCUITS CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents

More information

The Second Order Contribution to Wave Crest Amplitude Random Simulations and NewWave

The Second Order Contribution to Wave Crest Amplitude Random Simulations and NewWave The Secon Orer Contribution to Wave Crest Amplitue Ranom Simulations an NewWave Thomas A.A. Acock Department of Engineering Science University of Oxfor Parks Roa Oxfor Unite Kingom Scott Draper School

More information

PARALLEL-PLATE CAPACITATOR

PARALLEL-PLATE CAPACITATOR Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.

More information

An Anisotropic Hardening Model for Springback Prediction

An Anisotropic Hardening Model for Springback Prediction An Anisotropic Harening Moel for Springback Preiction Danielle Zeng an Z. Ceric Xia Scientific Research Laboratories For Motor Company Dearborn, MI 48 Abstract. As more Avance High-Strength Steels (AHSS

More information

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies

Chapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation

More information

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of

More information

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Bair, faculty.uml.eu/cbair University of Massachusetts Lowell 1. Pre-Einstein Relativity - Einstein i not invent the concept of relativity,

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In the usual geometric derivation of Bragg s Law one assumes that crystalline Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A

More information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

Fall 2016: Calculus I Final

Fall 2016: Calculus I Final Answer the questions in the spaces provie on the question sheets. If you run out of room for an answer, continue on the back of the page. NO calculators or other electronic evices, books or notes are allowe

More information

Situation awareness of power system based on static voltage security region

Situation awareness of power system based on static voltage security region The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

Research Article When Inflation Causes No Increase in Claim Amounts

Research Article When Inflation Causes No Increase in Claim Amounts Probability an Statistics Volume 2009, Article ID 943926, 10 pages oi:10.1155/2009/943926 Research Article When Inflation Causes No Increase in Claim Amounts Vytaras Brazauskas, 1 Bruce L. Jones, 2 an

More information

Vertical shear plus horizontal stretching as a route to mixing

Vertical shear plus horizontal stretching as a route to mixing Vertical shear plus horizontal stretching as a route to mixing Peter H. Haynes Department of Applie Mathematics an Theoretical Physics, University of Cambrige, UK Abstract. The combine effect of vertical

More information

LeChatelier Dynamics

LeChatelier Dynamics LeChatelier Dynamics Robert Gilmore Physics Department, Drexel University, Philaelphia, Pennsylvania 1914, USA (Date: June 12, 28, Levine Birthay Party: To be submitte.) Dynamics of the relaxation of a

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

Impurities in inelastic Maxwell models

Impurities in inelastic Maxwell models Impurities in inelastic Maxwell moels Vicente Garzó Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain Abstract. Transport properties of impurities immerse in a granular gas unergoing homogenous

More information

A Modification of the Jarque-Bera Test. for Normality

A Modification of the Jarque-Bera Test. for Normality Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam

More information

Non-Equilibrium Continuum Physics TA session #10 TA: Yohai Bar Sinai Dislocations

Non-Equilibrium Continuum Physics TA session #10 TA: Yohai Bar Sinai Dislocations Non-Equilibrium Continuum Physics TA session #0 TA: Yohai Bar Sinai 0.06.206 Dislocations References There are countless books about islocations. The ones that I recommen are Theory of islocations, Hirth

More information

Damage identification based on incomplete modal data and constrained nonlinear multivariable function

Damage identification based on incomplete modal data and constrained nonlinear multivariable function Journal of Physics: Conference Series PAPER OPEN ACCESS Damage ientification base on incomplete moal ata an constraine nonlinear multivariable function To cite this article: S S Kourehli 215 J. Phys.:

More information

3.2 Shot peening - modeling 3 PROCEEDINGS

3.2 Shot peening - modeling 3 PROCEEDINGS 3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:

More information

Linearized Theory: Sound Waves

Linearized Theory: Sound Waves Linearized Theory: Sound Waves In the linearized limit, Λ iα becomes δ iα, and the distinction between the reference and target spaces effectively vanishes. K ij (q): Rigidity matrix Note c L = c T in

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945

More information

PHYS 414 Problem Set 2: Turtles all the way down

PHYS 414 Problem Set 2: Turtles all the way down PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian

More information

Placement and tuning of resonance dampers on footbridges

Placement and tuning of resonance dampers on footbridges Downloae from orbit.tu.k on: Jan 17, 19 Placement an tuning of resonance ampers on footbriges Krenk, Steen; Brønen, Aners; Kristensen, Aners Publishe in: Footbrige 5 Publication ate: 5 Document Version

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N.

Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N. Submitte to the Journal of Rheology Moeling the effects of polyispersity on the viscosity of noncolloial har sphere suspensions Paul M. Mwasame, Norman J. Wagner, Antony N. Beris a) epartment of Chemical

More information

arxiv:hep-th/ v1 3 Feb 1993

arxiv:hep-th/ v1 3 Feb 1993 NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Adhesive Wear Theory of Micromechanical Surface Contact

Adhesive Wear Theory of Micromechanical Surface Contact International Journal Of Computational Engineering esearch ijceronline.com Vol. Issue. hesive Wear Theory of Micromechanical Surface Contact iswajit era Department of Mechanical Engineering National Institute

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information