Lubrication Forces. Ranga University of Edinburgh, UK (Dated: June 21, 2017) I. COMPUTING FORCES AND TORQUES
|
|
- Donna Harmon
- 5 years ago
- Views:
Transcription
1 Lubrication Forces Ranga University of Edinburgh UK Dated: June I. COMPUTING FORCES AND TORQUES The lubrication force between a pair of particles say p = 1 q = are given as a product between the Gr Resistance Matrix which only depends on the positions of the two particles a vector containing particle background fluid velocities in Refs. [1 3]. Our task below is to simplify the relation between the lubrication forces torques given in Refs. [1 3] for spheres of different sizes into a form similar to the one given by Ball&Melrose for spheres of the same size [4]. 1 a FIG. 1. Sketch of the two spheres with the normal vector indicated by the arrow. We have sphere p = 1 of radius a at position x 1 separated by a distance h from sphere q = of radius b at position x as shown in Fig. 1. The superscripts denote the particle number later the subscript will denote the component index. The total center to center distance is given by R = x 1 x = a + b + h. The separation between the two spheres is non-dimensionalised as ξ = R a b a + b the normal vector connecting the two spheres is given by = n = x x 1 x x 1 = x x 1 R h a + b/ 1. The velocity vectors of the spheres are given by U 1 U the angular velocities by ω 1 ω. The background fluid velocity is assumed to be linear with uniform velocity given by v angular velocity given by Ω rate of strain tensor given by E. According to Ref. [] the lubrication force acting on sphere 1 is given by h b F 1 /µ =A v x 1 U 1 + A 1 v x U + B Ω ω 1 + B 1 Ω ω + G E + G 1 E 3 the torque is given by T 1 /µ =B v x 1 U 1 + B 1 v x U + C Ω ω 1 + C 1 Ω ω + H E + H 1 E 4 where the bold symbols denote resistance functions which operate on a vector say χ = χ j as specified below A pq ij = Xpq A n in j + Y pq A δ ij n i n j B pq ij = Y pq B ɛ ijkn k C pq ij = Xpq C n in j + Y pq C δ ij n i n j G pq ijk = Xpq G n in j δ ij /3n k + Y pq G n iδ jk + n j δ ik n i n j n k H pq ijk = Y pq H ɛ ikln l n j + ɛ jkl n l n i. 5
2 Further Bpq ij = B qp ji Gpq ijk = Gqp formulation. Therefore we also have jki Hpq ijk = Hqp jki are required due to symmetry reasons of the resistance matrix B ij = Bji = YB ɛ ijk n k B ij 1 = Bji 1 = YB 1 ɛ ijk n k G ijk = G jki = XG n j n k δ jk /3n i + YG n j δ ki + n k δ ji n j n k n i G 1 ijk = G 1 jki = XG 1 n j n k δ jk /3n i + YG 1 n j δ ki + n k δ ji n j n k n i H ijk = Hjki = YH ɛ ijl n l n k + ɛ ikl n l n j H ijk 1 = Hjki 1 = YH 1 ɛ ijl n l n k + ɛ ikl n l n j. 6 The different components of the resistance functions specified by symbols such as X pq A Y pq A... are given in Sec. VI the relations between them are specified in Eq. 54. It follows that the lubrication force given in Eq. 3 can be written in component form as F 1 i /µ = X A n i n j + A δ ij n i n j U j U 1 j ɛ ijk Y B Ω j ωj 1 + YB 1 Ω j ωj n k + XA n i n j + YA δ ij n i n j vj x 1 vj x + XG + XG 1 n j n k δ jk /3n i Ejk + YG + YG 1 n j δ ki + n k δ ji n j n k n i Ejk 7 the torque given by Eq. 4 in component form is T 1 i /µ =ɛ ijk B U j U 1 j n k ɛ ijk B v j x v j x 1 n k + δ ij n i n j C Ω j ω 1 j + Y 1 C Ω j ω j H + Y 1 H ɛ ijl n l n k + ɛ ikl n l n j E jk. 8 Our objective is to make use of the relations between the different terms in Eq. 54 the simplifications in Sec. V to convert the forces torques above into simpler forms.
3 3 II. COMPUTING THE FORCE ACTING ON THE SPHERES A. Grouping terms in Eq. 7 Making use of the relations in Eq. 54 Sec. V we can get X A n i n j v j x 1 v j x + X G + X 1 G n j n k δ jk /3n i E jk =X A E jkn i n j n k R + a + b = X A E jkn i n j n k a + bξ/ 0 9 A δ ij n i n j vj x 1 vj x + YG + YG 1 n j δ ki + n k δ ji n j n k n i Ejk = YA δij n i n j Ejkn k + ɛ ikl Ω k n l R + Y A δ ij n i n j n k Ejka + b = YA ɛ ikl Ω k n l R YA δ ij n i n j Ejkn k a + bξ/ YA ɛ ikl Ω k n l R 10 A ɛ ijk Ω j n k R B + Y 1 B ɛ ijk Ω j n k = A ɛ ijk Ω j n j R a + b = A ɛ ijk Ω j n j a + bξ/ 0 Using the results from Sec. II A Eq. 7 simplifies to B. Simplified Lubrication force acting on sphere 1 F 1 i /µ = X A n i n j + A δ ij n i n j U j U 1 j + ɛ ijk Y B ω 1 j + Y 1 B ω j nk. 1 Eq.1 can be written as F 1 i /µ = X A n i n j + A δ ij n i n j Uj Uj 1 1 a + b ɛ ijk A ω 1 j + ω j n k + ba + 4b a4a + b B ɛ ijk ω 1 j ω j n k /. 13 C. Force acting on sphere According to the form given by Ref. [] F /µ =A 1 v x 1 U 1 + A v x U + B 1 Ω ω 1 + B Ω ω 14 + G 1 E + G E And similar to Eq. 7 we can write Eq. 14 after using the relations given in Eq. 54 as F i /µ = X A n i n j + A δ ij n i n j U j U 1 j + ɛ ijk Y B Ω j ωj 1 + YB 1 Ω j ωj n k XA n i n j + YA δ ij n i n j vj x 1 vj x XG + XG 1 n j n k δ jk /3n i Ejk YG + YG 1 n j δ ki + n k δ ji n j n k n i Ejk. 15 Comparing Eq. 7 Eq. 15 we get F i = F 1 i.
4 4 III. LUBRICATION TORQUES A. Grouping terms in Eq. 8 Making use of the relations in Eq. 54 Sec. V we can get which implies that H + Y 1 H ɛ ijl n l n k + ɛ ikl n l n j E jk =a + b B ɛ ijl n l n k E jk YC + YC 1 δ ij n i n j Ω j = a + byb δ ij n i n j Ω j B ɛ ijl v j x v j x 1 n l = a + b1 + ξ/ B ɛijl n k n l E jk δ ij n i n j Ω j 16 YC + YC 1 δ ij n i n j Ω j ɛ ijk YB vj x vj x 1 n k YH + YH 1 ɛ ijl n l n k + ɛ ikl n l n j Ejk = ξa + byb ɛijl n k n l Ejk + δ ij n i n j Ω 17 j / = Oξ 0. Using the above equations we get B. Torque acting on sphere 1 T 1 i /µ = B ɛ ijk U j U 1 j n k δ ij n i n j C ω 1 j + Y 1 C ω j. 18 We can try further simplification by noting Therefore Eq. 18 becomes simply T 1 i /µ = B ɛ ijk U j U 1 j n k + a + b B δ ij n i n j ω 1 j + ω j / + 1 4a/bY 1 C δ ij n i n j ω 1 j ω j / 19 C. Torque on sphere For the torque acting on sphere we have which turns out to be T /µ =B 1 v x 1 U 1 + B v x U + C 1 Ω ω 1 + C Ω ω + H 1 E + H E T i /µ =ɛ ijk Y 1 B U j U 1 j n k ɛ ijk Y 1 B v j x v j x 1 n k + δ ij n i n j Y 1 C Ω j ω 1 j + Y C Ω j ω j Y H + Y 1 H ɛ ijl n l n k + ɛ ikl n l n j E jk. Analogous to the previous calculations this makes Eq. 1 to be Therefore Eq. becomes simply T i /µ = Y 1 B ɛ ijk U j U 1 j n k δ ij n i n j Y 1 C ω 1 j + Y C ω j. 0 1 T i /µ =Y 1 B ɛ ijk U j U 1 j n k + a + by 1 B δ ij n i n j ω 1 j + ω j / 1 4b/aδ ij n i n j Y 1 C ω 1 j ω j / 3 Equally Ti ba + 4b /µ = a4a + b Y B ɛ ijk Uj Uj 1 n k + 1 4b/aY 1 C δ ij n i n j ω 1 j ω j / ba + ba + 4b YB δ ij n i n j ωj 1 + ωj ab + 4a 4
5 5 IV. TESTING THE SIMPLIFIED EXPRESSIONS A. Comparison of the forces torques to Ref. [4] For equal spheres b = a Eq. 13 for the force becomes F 1 i /µ = X A n i n j + A δ ij n i n j U j U 1 j ɛ ijk A a ω 1 j + ω j nk = X A n i n j + A δ ij n i n j U 1 j + a ω 1 n j U j a ω n j. 5 is equal to the form given in Ref. [4]. For the torques we have T 1 i /µ = a A ɛ ijk U j U 1 j n k a A δ ij n i n j ω 1 j + ω j 3Y 1 C δ ij n i n j ω 1 j ω j / T i /µ = a A ɛ ijk U j U 1 j n k a A δ ij n i n j ω 1 j + ω j + 3Y 1 C δ ij n i n j ω 1 j ω j /. The first two terms are the same as Ref. [4]. However the third term is which is 4 times the value given in Ref. [4]. 3Y 1 C / = 3 5 πa3 ln1/ξ B. Shearing motion of two rigid surfaces This is an example problem in Chap. 9 of Ref. []. Sphere 1 is moving past a fixed sphere along the x-axis at a constant velocity U. We have U 1 = U ˆx U = 0 ω 1 = 0 ω = 0 v = 0. The force on particle 1 given by Eq. 1 after substituting for the different values of velocities becomes after substituting for A F 1 x = µ A U 6 from Sec. VI we get Fx 1 6πµaU = 4β + β + β β 3 ln ξ 1. 7 The torque on particle 1 given by Eq. 18 after substituting for the different values of velocities becomes T 1 y = µ B U 8 1 z a b U x FIG.. Sketch of the two spheres with the normal vector indicated by the arrow. after substituting for B from Sec. VI we get T 1 y 8πµa U These results are the same as that Eq. 9.4 Eq. 9.5 in Ref. []. = β4 + β β. 9
6 6 C. Shearing motion of two rigid surfaces due to rotation This is an example problem in Chap. 9 of Ref. []. Sphere 1 is rotating near a fixed sphere about the y-axis at a constant angular velocity ω. We have U 1 = 0 U = 0 ω 1 = ωŷ ω = 0 v = 0. The force on particle 1 given by Eq. 1 after substituting for the different values of velocities becomes after substituting for B F 1 x = µ B ω 30 from Sec. VI we get Fx 1 β4 + β 8πµa = ω β. 31 The torque on particle 1 given by Eq. 18 after substituting for the different values of velocities becomes after substituting for C T 1 y = µ C ω 3 from Sec. VI we get These results are the same as that Eq. 9.6 Eq. 9.7 in Ref. []. T 1 y 1 z a b x FIG. 3. Sketch of the two spheres with the normal vector indicated by the arrow. 8πµa 3 ω = β β. 33 D. Comparison of the expressions with lubricate/poly pair style in LAMMPS In the LAMMPS implementation of lubricate/poly we have where the velocities V 1 j Fi 1 /µ = XA n i n j + YA δ ij n i n j Vj Vj 1 34 are defined as The previous equations imply that Vj 1 = Uj 1 vj x 1 + aɛ jlm ωl 1 Ω l n m aejmn m 35 Vj = Uj vj x bɛ jlm ωl Ω l n m + bejmn m. 36 V j V 1 j = U j U 1 j ɛ jlm aω 1 l + bω l n m v j x v j x 1 + a + be jm + ɛ jlm Ω l n m U j U 1 j ɛ jlm aω 1 l + bω l n m = U j bɛ jlm ω l n m U 1 j + aɛ jlm ω 1 l n m 37 by using the definition for the infinite velocities given in Sec. V noting that R a + b. Using the above Eq. 34 becomes Fi 1 /µ = XA n i n j + YA δ ij n i n j Uj Uj 1 ɛ ilm YA aωl 1 + bωl n m 38 which is equal to Eq. 5 when both particles are of the same size a=b if X A = X A Y A = Y A. 1. Difference between the terms in lubricate/poly Eq. 1 If we set XA = X A Y A = Y A in Eq. 38 we compare it with Eq. 1 we obtain that the pre-factors to the tangential forces due to rotation of the spheres are different between them while the force contribution due to linear velocity components are the same. The difference in the rotational contributions can be found out to be B + a A = 1πa b b a 5a + b 3 log 1/ξ Y 1 B + b A = 1πa b a b 5a + b 3 log 1/ξ 39 for ω 1 ω contributions to the rotational forces respectively. For a typical case of particles with the aspect ratio β = b/a = 1.4 we can calculate that lubricate/poly will under-estimate force due to rotation on particle 1.% over-estimate the force due to rotation on particle by 8.% even if other terms in the force class were corrected.
7 7 V. SIMPLIFICATIONS By definition of the rate of strain tensor E angular velocity vector ω we know that v j x = E jk x k + ɛ jkl Ω k x l. It implies Taking the dot product of Eq. 40 by n i n j we have v j x v j x 1 = E jkx k x 1 k + ɛ jlk Ω l x k x 1 k = E jkn k + ɛ jlk Ω l n k R. 40 n i n j vj x vj x 1 = Ejkn k n i n j + ɛ jlk Ω l n k n i n j R = Ejkn i n j n k R by dot product with the unit tensor δ ij vj x vj x 1 = Ejkn k δ ij + ɛ jlk Ω l n k δ ij R = Eik n k + ɛ ilk Ω l n k R. From the above two equations we have 41 4 Similarly we can also get δ ij n i n j v j x v j x 1 = E jkδ ij n i n j n k + ɛ ilk Ω l n k R. 43 ɛ ijk vj x vj x 1 n k = Rɛ ijk Ejmn m + ɛ jlm Ω l n m n k = R ɛ ijk Ejmn m n k δ kl δ im δ km δ il Ω m n l n k = R ɛ ijk Ejmn m n k δ im n m n i Ω m = 1 + ξ/a + b ɛ ijk E jmn m n k δ ij n i n j Ω j. 44 VI. NEAR FIELD FORMS OF THE RESISTANCE FUNCTIONS For the sake of completeness the resistance functions as specified in Chap. of Ref. [] up to Oln ξ 1 are given in this section. Ref. [3] only includes the leading terms of the following expressions its normal vector n points in the opposite direction to ours. Note that β = b/a for all the functions below the symmetry relations for the resistance matrix Chap. 7 Ref. [] are satisfied by the expressions below. The components of the A resistance functions are β XA = 6πa 1 + β 3 ξ 1 + β1 + 7β + β 51 + β 3 ln ξ 1 β XA = 6πb 1 + β 1 3 ξ 1 + β β 1 + β 51 + β 1 3 ln ξ 1 4β + β + β YA 45 = 6πa β 3 ln ξ 1 4β YA 1 + β 1 + β = 6πb β 1 3 ln ξ 1 The components of the B resistance functions are XA = XA = XA 1 = XA 1 YA = YA = YA 1 = YA 1. β4 + β B = 4πa 51 + β Y β B = 4πb β β
8 8 YB = YB 1 YB 1 = YB. 48 The components of the C resistance functions are β YC = 8πa 3 YC 1 = 8πa 3 β 51 + β β β YC = 8πb β 1 YC 1 = 8πb 3 β β 1 49 X pq C = O1. 50 The components of the G resistance functions are 3β XG = 4πa 1 + β 3 ξ 1 + 3β1 + 1β 4β β 3 ln ξ 1 3β XG = 4πb 1 + β 1 3 ξ 1 + 3β β 1 4β β 1 3 ln ξ 1 β4 β + 7β YG = 4πa β 3 ln ξ 1 Y G = 4πb β 1 4 β 1 + 7β β 1 3 ln ξ 1 51 XG = XG 1 XG = XG 1 YG = YG 1 YG = YG 1. 5 The components of the H resistance functions are Y β β H = 8πa β β H = 8πb 3 1 β β 1 Y 1 β H = 8πa β 01 + β Y β H = 8πb β β 1. 53
9 9 A. Relations between different components of the resistance functions To make the analysis simple we can note that the following relations based on the formulae given in Ref. [] for the leading terms Oξ 1 Olnξ 1 of the lubrication force C = 4 a b Y 1 C Y C = 4 b a Y 1 C C = a b Y C Y 1 C = Y 1 C X G + X 1 G = a + bx A G + Y 1 G = a + by A B + Y 1 B = a + b A YC + YC 1 = a + byb YC 1 + YC = a + byb 1 YC + YC 1 = YH + YH 1 YC 1 + YC = YH 1 + YH. Other simplifications can be done based on the above relations 54 YB ωj 1 + YB 1 ωj = YB + YB 1 ωj 1 + ωj / + YB YB 1 ωj 1 ωj / 55 = a + bya ωj 1 + ωj ba + 4b / + 1 YB ωj 1 ωj / a4a + b 56 C ω 1 j + Y 1 C ω j = C + Y 1 C ω 1 j + ω j / + C Y 1 C ω 1 j ω j / = a + b B / + 4a/b 1Y 1 C / Y 1 C ω 1 j + Y C ω j = Y 1 C + Y C ω 1 j + ω j / + Y 1 C Y C ω 1 j ω j / = a + by 1 B / + 1 4b/aY 1 C / VII. STRESSLETS AND STRESS According to Ref. [] the stresslet acting on particle 1 is given by S 1 /µ =G v x 1 U 1 + G 1 v x U + H Ω ω 1 + H 1 Ω ω + M E + M 1 E 59 where the resistance function M is given by M pq ijkl = Xpq M d0 ijkl + Y pq M d1 ijkl + Zpq M d ijkl 60 d 0 ijkl = 3 n in j δ ij /3n k n l δ kl /3 d 1 ijkl = 1 n in k δ jl + n j n k δ il + n i n l δ jk + n j n l δ ik 4n i n j n k n l 61
10 10 we can neglect Z pq M = O1 contributions. We can also find the following relationships to simplify our calculations including terms only upto Olog 1/ξ XM + XM 1 = a + bxg /3 YM + YM 1 = a + byg YH + YH 1 = a + byg. 6 Eq. 59 can be written in index notation as S 1 ij/µ = X G n i n j δ ij /3n k + G n i δ jk + n j δ ik n i n j n k v k x 1 U 1 k + XG 1 n i n j δ ij /3n k + YG 1 n i δ jk + n j δ ik n i n j n k vk x Uk + YH ɛ ikl n l n j + ɛ jkl n l n i Ω k ωk 1 + YH 1 ɛ ikl n l n j + ɛ jkl n l n i Ω k ωk + + XM d 0 M d 1 ijkl ijkl + XM 1 d 0 ijkl + Y M 1 d 1 ijkl Ekl Ekl. 63 The above can be simplified to be S 1 ij/µ = X G n i n j δ ij /3n k + G n i δ jk + n j δ ik n i n j n k U k U 1 k ɛ ikl n l n j + ɛ jkl n l n i YH ωk 1 + YH 1 ωk XG n i n j δ ij /3n k + YG n i δ jk + n j δ ik n i n j n k Ekl n l + ɛ kql Ω q n l R + a + byg ɛ ikl n l n j + ɛ jkl n l n i Ω k + a + bxg d 0 ijkl /3 + a + byg d 1 ijkl Ekl. 64 Grouping G terms in the above equation we have 1 = n i δ jk + n j δ ik n i n j n k E kl n l + ɛ kql Ω q n l R + d 1 ijkl E kl + ɛ ikl n l n j + ɛ jkl n l n i Ω k which we can simplify using the following identitities a + b 65 d 1 ijkl E kl = 1 n in k δ jl + n j n k δ il + n i n l δ jk + n j n l δ ik 4n i n j n k n l E kl = n i n k δ jl n j n l E kl n i δ jk + n j δ ik n i n j n k E kl n l = n i n k δ jl n j n l E kl ɛ ikl n l n j + ɛ jkl n l n i Ω k = n i ɛ jkl Ω k n l n i δ jk + n j δ ik n i n j n k ɛ kql Ω q n l = n i ɛ jkl Ω k n l 66 to Grouping some of the X G terms in Eq. 64 we have 1 = Oξ. 67 = n i n j n k n k δ ij /3 E kl n l + ɛ kql Ω q n l R + /3d 0 ijkl E kl a + b 68 which can again be simplified based on the following relation /3d 0 ijkl E kl = n i n j δ ij /3n k n l δ kl /3d 0 ijkl E kl = n i n j n k n l E kl 69 to be = Oξ. 70
11 From 1 Eq. 64 becomes S 1 ik/µ = X G n i n k δ ik /3n j + n k G δ ji n j n i U j U 1 j n i ɛ kjl n l H ω 1 j + Y 1 H ω j. 71 Similarly we can get S ik/µ = X 1 G n i n k δ ik /3n j + n k Y 1 G δ ji n j n i U j U 1 j n i ɛ kjl n l Y 1 H ω 1 j + Y H ω j. 7 Note that the stresslets on particles 1 are not equal in magnitude i.e. S 1 ik S ik. The stress acting between the two particles can be obtained as a sum of the individual stresslets since the total stress is a sum of all the stresslets acting on all the particles. So we have S 1 ik + S ik/µ = X G + X 1 G n i n k δ ik /3n j + n k G + Y 1 G δ ji n j n i U j U 1 j n i ɛ kjl n l H + Y 1 H ω 1 j + Y 1 H + Y H ω j =a + b X A n i n j + A δ ij n i n j U j U 1 j n k + a + bɛ kjl B ω 1 j + Y 1 B ω j n l n i a + bx A U j U 1 j n j δ ik /3 73 The anti-symmetric part of the stress acting between two particles is given by A 1 ij/µ = X G + X 1 G n i n k δ ik /3n j + n k G + Y 1 G δ ji n j n i U j U 1 j n i ɛ kjl n l H + Y 1 H ω 1 j + Y 1 H + Y H ω j =a + b X A n i n j + A δ ij n i n j U j U 1 j n k + a + bɛ kjl B ω 1 j + Y 1 B ω j n l n i a + bx A U j U 1 j n j δ ik /3 74 Let us now compare the stress tensor obtained by dyadic of force of particle 1 to the distance vector which is F 1 i n k R/µ = X A n i n j + A δ ij n i n j U j U 1 j n k R + Rɛ kjl Y B ω 1 j + Y 1 B ω j nl n i. 75 The difference between the two equations apart from the Oξ differences is only in the isotropic part δ ik which should vanish if the stress satisfies continuity equation σ = 0. [1] D. J. Jeffrey Physics of Fluids A: Fluid Dynamics [] S. Kim S. J. Karrila Microhydrodynamics Butterworth-Heinemann [3] R. Mari R. Seto J. F. Morris M. M. Denn Journal of Rheology [4] R. Ball J. Melrose Physica A: Statistical Mechanics its Applications
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationFaculty of Engineering, Mathematics and Science School of Mathematics
Faculty of Engineering, Mathematics and Science School of Mathematics JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Module MA3429: Differential Geometry I Trinity Term 2018???, May??? SPORTS
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationEXERCISES FOR CHAPTER 2
Agronomy 615 Statistical Genetics W. E. Nyquist March 2004 EXERCISES FOR CHAPTER 2 Exercise 2.1. Define the word population as you think it should be defined and discuss your understanding of the use of
More informationIndex Notation for Vector Calculus
Index Notation for Vector Calculus by Ilan Ben-Yaacov and Francesc Roig Copyright c 2006 Index notation, also commonly known as subscript notation or tensor notation, is an extremely useful tool for performing
More informationA.1 Appendix on Cartesian tensors
1 Lecture Notes on Fluid Dynamics (1.63J/2.21J) by Chiang C. Mei, February 6, 2007 A.1 Appendix on Cartesian tensors [Ref 1] : H Jeffreys, Cartesian Tensors; [Ref 2] : Y. C. Fung, Foundations of Solid
More information1. Basic Operations Consider two vectors a (1, 4, 6) and b (2, 0, 4), where the components have been expressed in a given orthonormal basis.
Questions on Vectors and Tensors 1. Basic Operations Consider two vectors a (1, 4, 6) and b (2, 0, 4), where the components have been expressed in a given orthonormal basis. Compute 1. a. 2. The angle
More informationUNIVERSITY OF DUBLIN
UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationBy drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.
Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,
More informationMATH45061: SOLUTION SHEET 1 V
1 MATH4561: SOLUTION SHEET 1 V 1.) a.) The faces of the cube remain aligned with the same coordinate planes. We assign Cartesian coordinates aligned with the original cube (x, y, z), where x, y, z 1. The
More informationVECTORS, TENSORS AND INDEX NOTATION
VECTORS, TENSORS AND INDEX NOTATION Enrico Nobile Dipartimento di Ingegneria e Architettura Università degli Studi di Trieste, 34127 TRIESTE March 5, 2018 Vectors & Tensors, E. Nobile March 5, 2018 1 /
More informationCartesian Tensors. e 2. e 1. General vector (formal definition to follow) denoted by components
Cartesian Tensors Reference: Jeffreys Cartesian Tensors 1 Coordinates and Vectors z x 3 e 3 y x 2 e 2 e 1 x x 1 Coordinates x i, i 123,, Unit vectors: e i, i 123,, General vector (formal definition to
More informationGenerating, from scratch, the near-field asymptotic forms of scalar resistance functions for two unequal rigid spheres in low-reynolds-number flow
Generating, from scratch, the near-field asymptotic forms of scalar resistance functions for two unequal rigid spheres in low-reynolds-number flow Adam K. Townsend arxiv:8.86v3 [physics.flu-dyn] 5 Jan
More informationSolutions. 3. (a) The easiest way to solve this problem is the following. Consider the vector: A i = ǫ ijk j k φ (16) Since ǫ ijk = ǫ ikj,
Solutions 7 Vector and tensor analysis: 1. 2. 3. a) The easiest way to solve this problem is the following. Consider the vector: A i = ǫ ijk j k φ 16) Since ǫ ijk = ǫ ikj, A i = ǫ ikj j k φ = ǫ ikj k j
More informationVector analysis and vector identities by means of cartesian tensors
Vector analysis and vector identities by means of cartesian tensors Kenneth H. Carpenter August 29, 2001 1 The cartesian tensor concept 1.1 Introduction The cartesian tensor approach to vector analysis
More informationHomework 7-8 Solutions. Problems
Homework 7-8 Solutions Problems 26 A rhombus is a parallelogram with opposite sides of equal length Let us form a rhombus using vectors v 1 and v 2 as two adjacent sides, with v 1 = v 2 The diagonals of
More informationThe Navier-Stokes equations for an incompressible fluid in indicial notation are. i u i = 0 (9.1)
Chapter 9 Viscous flow The Navier-Stokes equations for an incompressible fluid in indicial notation are i u i = 0 9.1 ρ t + u j j u i = i p + µ 2 j u i 9.2 When these are applied to practical problems,
More informationFormalism of the Tersoff potential
Originally written in December 000 Translated to English in June 014 Formalism of the Tersoff potential 1 The original version (PRB 38 p.990, PRB 37 p.6991) Potential energy Φ = 1 u ij i (1) u ij = f ij
More informationLECTURE 9: MOVING FRAMES IN THE NONHOMOGENOUS CASE: FRAME BUNDLES. 1. Introduction
LECTURE 9: MOVING FRAMES IN THE NONHOMOGENOUS CASE: FRAME BUNDLES 1. Introduction Until now we have been considering homogenous spaces G/H where G is a Lie group and H is a closed subgroup. The natural
More informationNIELINIOWA OPTYKA MOLEKULARNA
NIELINIOWA OPTYKA MOLEKULARNA chapter 1 by Stanisław Kielich translated by:tadeusz Bancewicz http://zon8.physd.amu.edu.pl/~tbancewi Poznan,luty 2008 ELEMENTS OF THE VECTOR AND TENSOR ANALYSIS Reference
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More information1.3 LECTURE 3. Vector Product
12 CHAPTER 1. VECTOR ALGEBRA Example. Let L be a line x x 1 a = y y 1 b = z z 1 c passing through a point P 1 and parallel to the vector N = a i +b j +c k. The equation of the plane passing through the
More informationDiscontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method
Discontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method Pradipto and Hisao Hayakawa Yukawa Insitute for Theoretical Physics Kyoto University Rheology of disordered
More informationOn Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers
Continuum Mech. Thermodyn. (1996) 8: 247 256 On Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers I-Shih Liu Instituto de Matemática Universidade do rio de Janeiro, Caixa Postal
More informationCharacterization of 2 n -Periodic Binary Sequences with Fixed 2-error or 3-error Linear Complexity
Characterization of n -Periodic Binary Sequences with Fixed -error or 3-error Linear Complexity Ramakanth Kavuluru Department of Computer Science, University of Kentucky, Lexington, KY 40506, USA. Abstract
More informationShear stress of a monolayer of rough spheres
J. Fluid Mech. (2002), vol. 452, pp. 425 441. c 2002 Cambridge University Press DOI: 10.1017/S002211200100688 Printed in the United Kingdom 425 Shear stress of a monolayer of rough spheres By HELEN J.
More informationNew Algorithms for Solving Interval Valued Fuzzy Relational Equations
Intern. J. Fuzzy Mathematical Archive Vol. 1, 2013, 92-98 ISSN: 2320 3242 (P), 2320 3250 (online) Published on 24 April 2013 www.researchmathsci.org International Journal of New Algorithms for Solving
More information1. Tensor of Rank 2 If Φ ij (x, y) satisfies: (a) having four components (9 for 3-D). (b) when the coordinate system is changed from x i to x i,
1. Tensor of Rank 2 If Φ ij (x, y satisfies: (a having four components (9 for 3-D. Φ i j (x 1, x 2 = β i iβ j jφ ij (x 1, x 2. Example 1: ( 1 0 0 1 Φ i j = ( 1 0 0 1 To prove whether this is a tensor or
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
More informationTensor Calculus. arxiv: v1 [math.ho] 14 Oct Taha Sochi. October 17, 2016
Tensor Calculus arxiv:1610.04347v1 [math.ho] 14 Oct 2016 Taha Sochi October 17, 2016 Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT. Email: t.sochi@ucl.ac.uk.
More informationThermodynamics for fluid flow in porous structures
Communications to SIMAI Congress, ISSN 1827-9015, Vol. 1 (2006) DOI: 10.1685/CSC06105 Thermodynamics for fluid flow in porous structures M.E. Malaspina University of Messina, Department of Mathematics
More informationVectors. (Dated: August ) I. PROPERTIES OF UNIT ANSTISYMMETRIC TENSOR
Vectors Dated: August 25 2016) I. PROPERTIE OF UNIT ANTIYMMETRIC TENOR ɛijkɛ klm = δ il δ jm δ im δ jl 1) Here index k is dummy index summation index), which can be denoted by any symbol. For two repeating
More informationModule 4M12: Partial Differential Equations and Variational Methods IndexNotationandVariationalMethods
Module 4M12: Partial Differential Equations and ariational Methods IndexNotationandariationalMethods IndexNotation(2h) ariational calculus vs differential calculus(5h) ExamplePaper(1h) Fullinformationat
More informationChapter 7. Kinematics. 7.1 Tensor fields
Chapter 7 Kinematics 7.1 Tensor fields In fluid mechanics, the fluid flow is described in terms of vector fields or tensor fields such as velocity, stress, pressure, etc. It is important, at the outset,
More informationElastic Fields of Dislocations in Anisotropic Media
Elastic Fields of Dislocations in Anisotropic Media a talk given at the group meeting Jie Yin, David M. Barnett and Wei Cai November 13, 2008 1 Why I want to give this talk Show interesting features on
More informationElastic Wave Theory. LeRoy Dorman Room 436 Ritter Hall Tel: Based on notes by C. R. Bentley. Version 1.
Elastic Wave Theory LeRoy Dorman Room 436 Ritter Hall Tel: 4-2406 email: ldorman@ucsd.edu Based on notes by C. R. Bentley. Version 1.1, 20090323 1 Chapter 1 Tensors Elastic wave theory is concerned with
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationThe two loop D0-brane metric in N = 2 sigma models
The two loop D0-brane metric in N = sigma models Thomas Wynter SPhT, CEA-Saclay, 91011 Gif-sur Yvette, France Abstract We investigate, up to two loop order, the physical metric seen by a D0-brane probe
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The Koiter shell equation in a coordinate free description CSC/1-0 Chemnitz Scientific Computing Preprints ISSN 1864-0087 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific
More informationMATH45061: SOLUTION SHEET 1 II
MATH456: SOLUTION SHEET II. The deformation gradient tensor has Cartesian components given by F IJ R I / r J ; and so F R e x, F R, F 3 R, r r r 3 F R r, F R r, F 3 R r 3, F 3 R 3 r, F 3 R 3 r, F 33 R
More informationInvariant multipole theory of induced macroscopic fields in homogeneous dielectrics
Invariant multipole theory of induced macroscopic fields in homogeneous dielectrics by Allard Welter Submitted in fulfilment of the academic requirements for the degree of PhD in Physics in the School
More informationPaolo Lorenzoni Based on a joint work with Jenya Ferapontov and Andrea Savoldi
Hamiltonian operators of Dubrovin-Novikov type in 2D Paolo Lorenzoni Based on a joint work with Jenya Ferapontov and Andrea Savoldi June 14, 2015 Paolo Lorenzoni (Milano-Bicocca) Hamiltonian operators
More informationMicrostructure and Particle-Phase Stress in a Dense Suspension
Microstructure and Particle-Phase Stress in a Dense Suspension Jim Jenkins, Cornell University Luigi La Ragione, Politecnico di Bari Predict microstructure and particle stresses in slow, steady, pure straining
More informationOn Expected Gaussian Random Determinants
On Expected Gaussian Random Determinants Moo K. Chung 1 Department of Statistics University of Wisconsin-Madison 1210 West Dayton St. Madison, WI 53706 Abstract The expectation of random determinants whose
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationAdd-on unidirectional elastic metamaterial plate cloak
Add-on unidirectional elastic metamaterial plate cloak Min Kyung Lee *a and Yoon Young Kim **a,b a Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro, Gwanak-gu, Seoul,
More informationMathematical Preliminaries
Mathematical Preliminaries Introductory Course on Multiphysics Modelling TOMAZ G. ZIELIŃKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Vectors, tensors, and index notation. Generalization of the concept
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationAppendix A Spin-Weighted Spherical Harmonic Function
Appendix A Spin-Weighted Spherical Harmonic Function Here, we review the properties of the spin-weighted spherical harmonic function. In the past, this was mainly applied to the analysis of the gravitational
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationOn the bulk viscosity of suspensions
J. Fluid Mech. (2006), vol. 554, pp. 109 12. c 2006 Cambridge University Press doi:10.1017/s002211200600948 Printed in the United Kingdom 109 On the bulk viscosity of suspensions By JOHN F. BRADY, ADITYA
More informationA Brief Introduction to Tensor
Gonit Sora Article 1-13 Downloaded from Gonit Sora Gonit Sora: http://gonitsora.com Introductory Article Rupam Haloi Department of Mathematical Sciences, Tezpur University, Napaam - 784028, India. Abstract:
More informationClassical Mechanics Solutions 1.
Classical Mechanics Solutions 1. HS 2015 Prof. C. Anastasiou Prologue Given an orthonormal basis in a vector space with n dimensions, any vector can be represented by its components1 ~v = n X vi e i. 1)
More informationFree rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012
Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate
More informationPhysics 110. Electricity and Magnetism. Professor Dine. Spring, Handout: Vectors and Tensors: Everything You Need to Know
Physics 110. Electricity and Magnetism. Professor Dine Spring, 2008. Handout: Vectors and Tensors: Everything You Need to Know What makes E&M hard, more than anything else, is the problem that the electric
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationClassical Mechanics in Hamiltonian Form
Classical Mechanics in Hamiltonian Form We consider a point particle of mass m, position x moving in a potential V (x). It moves according to Newton s law, mẍ + V (x) = 0 (1) This is the usual and simplest
More informationLectures notes. Rheology and Fluid Dynamics
ÉC O L E P O L Y T E C H N IQ U E FÉ DÉR A L E D E L A U S A N N E Christophe Ancey Laboratoire hydraulique environnementale (LHE) École Polytechnique Fédérale de Lausanne Écublens CH-05 Lausanne Lectures
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Engineering Acoustics Session 1aEA: Thermoacoustics I 1aEA7. On discontinuity
More informationTensors, and differential forms - Lecture 2
Tensors, and differential forms - Lecture 2 1 Introduction The concept of a tensor is derived from considering the properties of a function under a transformation of the coordinate system. A description
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More information1 Curvature of submanifolds of Euclidean space
Curvature of submanifolds of Euclidean space by Min Ru, University of Houston 1 Curvature of submanifolds of Euclidean space Submanifold in R N : A C k submanifold M of dimension n in R N means that for
More informationQMI PRELIM Problem 1. All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work.
QMI PRELIM 013 All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work. Problem 1 L = r p, p = i h ( ) (a) Show that L z = i h y x ; (cyclic
More informationDynamics of the Coarse-Grained Vorticity
(B) Dynamics of the Coarse-Grained Vorticity See T & L, Section 3.3 Since our interest is mainly in inertial-range turbulence dynamics, we shall consider first the equations of the coarse-grained vorticity.
More informationExam in Fluid Mechanics 5C1214
Eam in Fluid Mechanics 5C1214 Final eam in course 5C1214 13/01 2004 09-13 in Q24 Eaminer: Prof. Dan Henningson The point value of each question is given in parenthesis and you need more than 20 points
More informationHomogenization in Elasticity
uliya Gorb November 9 13, 2015 uliya Gorb Homogenization of Linear Elasticity problem We now study Elasticity problem in Ω ε R 3, where the domain Ω ε possess a periodic structure with period 0 < ε 1 and
More informationUseful Formulae ( )
Appendix A Useful Formulae (985-989-993-) 34 Jeremić et al. A.. CHAPTER SUMMARY AND HIGHLIGHTS page: 35 of 536 A. Chapter Summary and Highlights A. Stress and Strain This section reviews small deformation
More informationTwo-variable Wilson polynomials and the generic superintegrable system on the 3-sphere
Two-variable Wilson polynomials and the generic superintegrable system on the 3-sphere Willard Miller, [Joint with E.G. Kalnins (Waikato) and Sarah Post (CRM)] University of Minnesota Special Functions
More informationDynamic simulation of flexible fibers composed of linked rigid bodies
Dynamic simulation of flexible fibers composed of linked rigid bodies Russell F. Ross Department of Engineering Mechanics and Astronautics, University of Wisconsin Madison, Madison, Wisconsin 53706 Daniel
More information7.1 Creation and annihilation operators
Chapter 7 Second Quantization Creation and annihilation operators. Occupation number. Anticommutation relations. Normal product. Wick s theorem. One-body operator in second quantization. Hartree- Fock
More information2 Tensor Notation. 2.1 Cartesian Tensors
2 Tensor Notation It will be convenient in this monograph to use the compact notation often referred to as indicial or index notation. It allows a strong reduction in the number of terms in an equation
More informationDiscontinuous Shear Thickening
Discontinuous Shear Thickening dynamic jamming transition Ryohei Seto, Romain Mari, Jeffrey F. Morris, Morton M. Denn Levich Institute, City College of New York First experimental data Williamson and Hecker
More informationHÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling
HÅLLFASTHETSLÄRA, LTH Examination in computational materials modeling TID: 2016-28-10, kl 14.00-19.00 Maximalt 60 poäng kan erhållas på denna tenta. För godkänt krävs 30 poäng. Tillåtet hjälpmedel: räknare
More informationLagrangian. µ = 0 0 E x E y E z 1 E x 0 B z B y 2 E y B z 0 B x 3 E z B y B x 0. field tensor. ν =
Lagrangian L = 1 4 F µνf µν j µ A µ where F µν = µ A ν ν A µ = F νµ. F µν = ν = 0 1 2 3 µ = 0 0 E x E y E z 1 E x 0 B z B y 2 E y B z 0 B x 3 E z B y B x 0 field tensor. Note that F µν = g µρ F ρσ g σν
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationMicrostructural studies for rheology. Chapter 7: Microstructural studies for rheology. Separation of length scales. Micro & macro views
Chapter 7: Microstructural studies for rheology Microstructural studies for rheology To calculate the flow of complex fluids, need governing equations, in particular, the constitutive equation relating
More informationA Brief Introduction to Tensors
A Brief Introduction to Tensors Jay R Walton Fall 2013 1 Preliminaries In general, a tensor is a multilinear transformation defined over an underlying finite dimensional vector space In this brief introduction,
More informationCONTINUUM MECHANICS. lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern
CONTINUUM MECHANICS lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern Contents Tensor calculus. Tensor algebra.................................... Vector algebra.................................
More informationChapter 6: Momentum Analysis
6-1 Introduction 6-2Newton s Law and Conservation of Momentum 6-3 Choosing a Control Volume 6-4 Forces Acting on a Control Volume 6-5Linear Momentum Equation 6-6 Angular Momentum 6-7 The Second Law of
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More information1 Index Gymnastics and Einstein Convention
Index Gymnastics and Einstein Convention Part (i) Let g δ µνe µ e ν be the metric. Then: Part (ii) x x g (x, x) δ µνe µ (x σ e σ) e ν (x ρ e ρ) δ µνx σ x ρ e µ (e σ) e ν (e ρ) δ µνx σ x ρ δ µ σ δ ν ρ δ
More informationStructure relations for the symmetry algebras of classical and quantum superintegrable systems
UNAM talk p. 1/4 Structure relations for the symmetry algebras of classical and quantum superintegrable systems Willard Miller miller@ima.umn.edu University of Minnesota UNAM talk p. 2/4 Abstract 1 A quantum
More informationContents. Physical Properties. Scalar, Vector. Second Rank Tensor. Transformation. 5 Representation Quadric. 6 Neumann s Principle
Physical Properties Contents 1 Physical Properties 2 Scalar, Vector 3 Second Rank Tensor 4 Transformation 5 Representation Quadric 6 Neumann s Principle Physical Properties of Crystals - crystalline- translational
More informationPHYS 705: Classical Mechanics. Rigid Body Motion Introduction + Math Review
1 PHYS 705: Classical Mechanics Rigid Body Motion Introduction + Math Review 2 How to describe a rigid body? Rigid Body - a system of point particles fixed in space i r ij j subject to a holonomic constraint:
More informationON THE GEOMETRY OF TANGENT BUNDLE OF A (PSEUDO-) RIEMANNIAN MANIFOLD
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I.CUZA IAŞI Tomul XLIV, s.i.a, Matematică, 1998, f1 ON THE GEOMETRY OF TANGENT BUNDLE OF A (PSEUDO-) RIEMANNIAN MANIFOLD BY V. OPROIU and N. PAPAGHIUC 0. Introduction.
More informationLecture XXXIV: Hypersurfaces and the 3+1 formulation of geometry
Lecture XXXIV: Hypersurfaces and the 3+1 formulation of geometry Christopher M. Hirata Caltech M/C 350-17, Pasadena CA 91125, USA (Dated: May 4, 2012) I. OVERVIEW Having covered the Lagrangian formulation
More information(a) Show that electric field and magnetic field have units (force)/area or energy/volume.
Problem. Units (a) how that electric field and magnetic field have units (force)/area or energy/volume. (b) A rule of thumb that you may need in the lab is that coaxial cable has a capcitance of 2 pf/foot.
More informationMACROSCOPIC DYNAMIC EQUATIONS FOR NEMATIC LIQUID CRYSTALLINE SIDE-CHAIN POLYMERS
Mol.Cryst.Liq.Cryst.199, 407-418, (1991) MACROSCOPIC DYNAMIC EQUATIONS FOR NEMATIC LIQUID CRYSTALLINE SIDE-CHAIN POLYMERS HARALD PLEINER and HELMUT R. BRAND,+ FB 7 Physik, Universität Essen, D 4300 Essen
More informationAnomalous effect of turning off long-range mobility interactions in Stokesian Dynamics
Anomalous effect of turning off long-range mobility interactions in Stokesian Dynamics Adam K. Townsend, a), b) and Helen J. Wilson Department of Mathematics, University College London, Gower Street, London
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. # 02 Conservation of Mass and Momentum: Continuity and
More informationProblem Set 4. Eshelby s Inclusion
ME34B Elasticity of Microscopic Structures Wei Cai Stanford University Winter 24 Problem Set 4. Eshelby s Inclusion Chris Weinberger and Wei Cai c All rights reserved Problem 4. (5 Spherical inclusion.
More informationPhysical Properties. Reading Assignment: 1. J. F. Nye, Physical Properties of Crystals -chapter 1
Physical Properties Reading Assignment: 1. J. F. Nye, Physical Properties of Crystals -chapter 1 Contents 1 Physical Properties 2 Scalar, Vector 3 Second Rank Tensor 4 Transformation 5 Representation Quadric
More informationIntroduction to Tensor Notation
MCEN 5021: Introduction to Fluid Dynamics Fall 2015, T.S. Lund Introduction to Tensor Notation Tensor notation provides a convenient and unified system for describing physical quantities. Scalars, vectors,
More informationCosmology. April 13, 2015
Cosmology April 3, 205 The cosmological principle Cosmology is based on the principle that on large scales, space (not spacetime) is homogeneous and isotropic that there is no preferred location or direction
More informationResearch Article A Note on Pseudo-Umbilical Submanifolds of Hessian Manifolds with Constant Hessian Sectional Curvature
International Scolarly Researc Network ISRN Geometry Volume 011, Article ID 374584, 1 pages doi:10.540/011/374584 Researc Article A Note on Pseudo-Umbilical Submanifolds of Hessian Manifolds wit Constant
More informationChapter 3 Stress, Strain, Virtual Power and Conservation Principles
Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress
More information