The Mathematics of Thermal Diffusion Shocks
|
|
- Melanie Francis
- 6 years ago
- Views:
Transcription
1 The Mathematics of Thermal Diffusion Shocks Vitalyi Gusev,WalterCraig 2, Roberto LiVoti 3, Sorasak Danworahong 4, and Gerald J. Diebold 5 Université du Maine, av. Messiaen, 7285 LeMans, Cedex 9 France, 2 Deartment of Mathematics McMaster University, Hamilton, Ontario, Canada L8S 4K, 3 Universita di Roma La Saienza, Diartimento di Energetica, Via Scara, 4 6 Roma,Italy 4 School of Science, Walailak University, 222 Thaiburi, Thasala District, Nakorn Si Thammarat,86, Thailand, 5 Deartment Chemistry, Brown University, Providence, RI, USA, 292 Thermal diffusion, also known as the Ludwig-Soret effect, refers to the searation of mixtures in a temerature gradient. For a binary mixture the time deendence of the change in concentration of each secies is governed by a nonlinear artial differential equation in sace and time. Here, an exact solution of the Ludwig-Soret equation without mass diffusion for a sinusoidal temerature field is given. The solution shows that counter-roagating shock waves are roduced which slow and eventually come to a halt. Exressions are found for the shock time for two limiting values of the starting density fraction. The effects of diffusion are shown on the develoment of the concentration rofile in time and sace are found by numerical integration of the artial integration of he nonlinear differential equation. I. Introduction The searation of liquid mixtures in a thermal gradient was discovered by Ludwig[] in 856, and was described theoretically somewhat later by Soret[2]. The searation of mixtures by a thermal gradient takes lace not only in liquids, but also in gases and solids as is described by nonequilibrium thermodynamics[3 6]. Recently, a new method of roducing the Ludwig-Soret effect has been introduced where crossed laser beams roduce thermal gradients[7 ]. The method relies on interference in the electric fields of two hase coherent light beams to write an otical grating in a weakly absorbing liquid mixture. At the nodes of the electric field, there is no absortion, while at the antinodes heat is deosited resulting in a temerature field that is sinusoidal in sace[]. A robe laser is directed at the grating at the Bragg angle, and roduces diffracted beams in resonse to either aeriodic change in absortion or by a eriodic change in index of refraction that forms in the mixture as a result of thermal diffusion. A recording of the intensity of the diffracted beam gives the time develoment of searation of the comonents of the mixture. Here we discuss solution of the Ludwig-Soret equation in a one-dimensional geometry for a sinusoidal temerature distribution. A derivation of an exact solution of the Ludwig-Soret roblem with negligible mass diffusion is given in Section II. The time develoment of the shock fronts[2] are given in Section III. Section IV gives an exression for the shock velocity and exressions for the shock formation time for two limiting cases of the starting dentsity fraction. The effects of diffusion are shown through numerical solutions to the Ludwig-Soret equation. Derivation of the concentration rofile in sace at longtimesisgiven. II. Solution to the Thermal Diffusion Equation Thermal diffusion in a binary mixture is governed by the couled differential equations [4] T ρc P = (λ T + κ c ) () t c = (c c 2 D T + D c ), (2) t where t is the time, ρ is the solution density, c P is the secific heat caacity, λ is the thermal conductivity, T is the temerature above ambient, κ is a constant, D is the thermal diffusion coefficient, D is the mass diffusion coefficient, and c and c 2 are density fractions of the secies that make u the solution. Since the total mass of the solution is the sum of the masses of the two comonents of the mixture, the densities of the two secies satisfy the relation ρ + ρ 2 = ρ, whereρ and ρ 2 are the densities of the two comonents and ρ is the density of the mixture. If the exression relating the densities of the two secies is divided by ρ, the result is an exression for the density fractions c and c 2 of each of the comonents which obey the relation c + c 2 =. The first term on the right hand side of Eq. is the Fourier law for heat conduction; the second describes the Dufour effect, which refers to roduction of a temerature change as a result of imosition of a concentration gradient. The first term on the right hand side of Eq. 2 describes thermal diffusion, while the seconddescribesmassdiffusion. Consider a temerature field of the form T = T [ + sin(kx)], where 2T is the eak temerature, K is a wavenumber determined by the otical fringe sacing in the grating, and x is the coordinate. If the Dufour effect is taken as negligible [4], then Eq.2 becomes c(z,t) τ = α z {c(z, t)[ c(z,t)] cos z} + 2 c(z,t) z 2, (3) Tyeset by REVTEX
2 2 where c 2 has been written in terms of c, and the subscrit droed from the density fraction c, where the thermal diffusion factor α, isdefined as α = D T /D, and the dimensionless quantities τ and z are defined by τ = K 2 Dt, and z = Kx. If the last term in Eq. 3 describing diffusion is ignored, then the differential equation of motion for c(z,t) can be written c τ = f z, (4) where a flux f(c, z) is defined as f(c, z) = αc( c)cosz and τ is given by τ = K 2 D. Equation 4 is the differential form of a conservation equation that exresses the buildu of c in a volume as a consequence of a flux change in sace. Since for a eriodic temerature field the density fraction must also be eriodic in z, it follows that c(2 π, τ) = c(,τ). Integration of Eq. 4 over one eriod of the temerature field, i.e. from z =to z =2π, gives the integral form of Eq. 4 as 2π c(z, t) dz =2πc, (5) where c is the density fraction at time t =,assumed to be a constant throughout the cell. The integral for the density fraction over z exresses simle mass conservation for the Ludwig-Soret effect. The Eulerian descrition of the rofile c = c(z, τ) by Eq. 3 can be transformed[3] into a Lagrangian descrition yielding the couled air of ordinary differential equations, dz f(c, β) = = α(2c ) cos z (6) dτ c β) = f(c, = αc( c)sinz, (7) dτ β that determine the motion of oints with coordinates z = z(τ,c, z ) and c = c(τ,c, z ) on the zc lane for a oint initially at (c, z ) at time τ =. As discussed in Ref. [3], Eqs. 6 and 7 form a Hamiltonian system of equations as found in classical mechanics, with the flux function in the resent roblem taking on the role of the Hamiltonian function in classical mechanics. Numerical integration of Eqs. 6 and 7 can be used to give the z coordinate and the value of c as a function of time for any oint initially at (z,c ), as shown in Fig.. Since the temerature field has a eriod of 2π, the solutions to Eq. 3 must also have the same eriodicity. As can be seen from Fig., the oints π/2 and the oint 3π/2 reresent cold oints of the grating so that the solutions to Eq. 3 are identical at these oints, and for any oints dislaced from these oints by a multile of 2π. From FIG. : Density fraction c versus dimensionless distance z for τ =,.5, 3.5, and 6., for a starting density fraction c =.3. The curves can be determined by numerical integration of Eqs. 6 and 7 or through solution of Eqs. 2 and 3. The vertical lines indicates the ositions of the shock fronts. The shock velocity is found by determining the line that makes the areas of the darkened regions adjacent to each of the vertical lines equal. Eqs. 6 and 7 the sloe of the curves on a hase ortrait is found to be c( c)sinz = dz (2c ) cos z, which can be written as d [c( c)cosz] =. (8) dz It follows then that the quantity in brackets in Eq. 8 is a constant, taken here to be k, so that c( c)cosz = k. (9) The hase ortrait for the Ludwig-Soret equation, as shown in Fig. 2, gives the trajectories of oints with secified starting values of c and z.by differentiation of Eq. 9 with resect to c, it can be seen that the maximum value of the left hand side occurs when cos z =,an =/2; hence the values of k are restricted to k /4. By solving Eq. 7 and Eq. 9 for sin z and cos z, resectively, squaring and adding the resulting exressions, it is easily shown that k c( c) 2 which can be written + c( c) c2 ( c) 2 k 2 2 µ 2 dτ =, = ±dτ. ()
3 3 FIG. 2: Phase ortrait for the density fraction c versus dimensionless distance along the grating z from Eq. 9. Points on the zc lanemovealongcurvesofconstantk as time rogresses. It is convenient to define the dearture of the density fraction from a value of /2 through the variable c,and to exress the quantity in the denominator of Eq. in terms of two arameters a and b, defined by a = q c = c k b = q 4 k. Owing to the restriction on the magnitude of k, the inequality a b is automatically satisfied. Hence, Eq. then can be integrated from a starting density fraction c at time τ =along a curve of constant k to a density fraction c at time τ giving c 2 c 2 τ (a2 c 2 )(b 2 c 2 ) = ± dτ, () where the integration is taken over a trajectory in hase sace with a constant value of k. The left hand side of Eq. can be exressed as an ellitic integral[4] of the first kind F, so that the solution to the Ludwig-Soret roblem for a sinusoidal temerature field, neglecting diffusion is F arcsin c 2, 4 k 4 k 4 +k q (2) F arcsin c 2, 4 k 4 k = ±τ 4 + k. 4 +k It can be seen that Eq. 2 gives τ exlicitly for any starting oint (z,c ) that arrives at the newoint c; the new coordinate z is determined through c ( c )cosz = c( c)cosz = k. (3) FIG. 3: Density fraction c versus dimensionless distance z from Eq. 6. with α =for τ =,, 2,and 3.. The oints lotted as dots ( ) are those that originally lie in the region <z<π/2 at τ =. The oints lotted as squares ( ) are those that lie in the region π/2 <z< at time τ =. From the oint of view of calculating density fraction rofiles as shown in Fig., it is convenient to secify values of z, c and τ and to solve Eq. 2 for the new values of c through a numerical root search. Alternately, the oints (c, z) can be substituted directly into Eq. 2 throughuseofeq. 3todeterminek. Note that for ositive values of α, theminussignisusedineq.2for oints moving in the hot regions of the grating <z<π, and the lus sign is used for motion in the cold regions π<z<2π, or π <z<. Inversion of Eq. 2 to give c exlicitly as a function of z for any given time τ, does not aear to be ossible; however it is straightforward to determine the rofiles of the density fraction versus coordinate through numerical methods. III. Motion of the Points in Phase Sace Consider the motion of a oint located between and π/2 with a value of c < /2. It can be seen from Fig. that oints in this region move downwards and towards the c axis. For such oints the motion is determined first by Eq. 2 with a negative sign is used, and after assing the c axis at z =, by Eq. 2 with the ositive sign. Any oint in this region starting at (c,z ) reaches its minimum value at z =, which from Eq. 3, is given by = 2 r 4 c ( c )cosz. (4)
4 4 The time required to arrive at is therefore τ min = c min c x2 ( x) 2 [c ( c )cosz ] 2, (5) and the total time τ, from the starting oint to the oint (c, z) is τ = τ min+ c x2 ( x) 2 [c ( c )cosz ] 2. (6) which can be written in terms of ellitic integrals as in Eq. 2. Since the shocks in the region π/2 <z< originate from oints in the region >z<π/2, Eq. 2is exression that determines the time of shock formation. Similar reasoning gives the maximum density fraction c max for any oint starting at (c,z ) located between and π/2 as c max = 2 + r 4 c ( c )cosz, (7) where τ max, thetimeittakestomovefromc to to c max, is given by τ max = τ min + c max x2 ( x) 2 [c ( c )cosz ] 2. (8) The curves in Fig. 3 show how oints originally in the region between and π/2 move downwards to their minimum values at z =, ass to the left of the c axis and then rise uwards as time rogresses. IV Shock Formation A. Determination of the shock velocity It is clear from examination of Fig. 2 that the solutions to Eq. 2 exhibit nonfunctional behavior after a certain time. Similar solutions for different values of c show that nonfunctional behavior is found for any value of c if τ is allowed to become large. When the sloe of any curve becomes unbounded, i.e.. immediately before the density fraction rofile becomes nonfunctional, shocks are formed, and the mathematics must be treated in a different manner. Consider the vertical line to the right of the density fraction maximum at 3π/2 indicating the osition of the shock front on the curve in Fig. for τ =3.5. For the shock front, the density fraction immediately to the left of the front is denoted c l and that to the right of the shock by c r. The area of the curve bisected by the vertical line, can be equated to the roduct of the velocity of the shock dz sh /dτ, the difference in c and the time interval dτ,that is, dz cl sh dτ (c l c r )dτ = dz, c r which, together with Eq. 6, gives dz cl sh dτ (c l c r )dτ = α cos zdτ (2c ). (9) c r After evaluation of the integral in Eq. 9 the shock velocity becomes dz sh = α(c l + c r ) cos z. (2) dτ It can be seen that two shocks, one on either side of the oint z =3π/2 are formed, each of which travels with the same seed towards the cold region of the grating, the different velocities being determined by the change in sign of the cosine function on either side of z =3π/2. The left-going shock comes to a halt when c l =and c r =, that is, where a comlete searation of the mixture has been attained, at which oint the rofile of c in sace is asquarewavecenteredatz =3π/2. The shock velocity can be obtained equally by writingeq. 6asdz sh /dτ = [f(c l,z) f(c r,z)]/(c l c r ), which, after substitution of the exression for the flux function, gives Eq. 2. As noted in Ref. [2], Eq. 2, which exresses the thermal diffusion shock velocity in terms of density fractions on either side of the shock, is an exact analog of the Rankine-Hugoniot relations for one-dimensional fluid shocks: Eq. 2 gives the thermal diffusion shock seed in terms of the "state variables" c l and c r, the Rankine-Hugoniot relations exress the shock seed in terms of ratios of the state variables of the fluid on either side of the shock. The time of formation of the shock can be found from numerical integration of the Hamiltonian system of equations and finding the oint in sace and time where the rofile reverses direction. Figure 4 shows a lot of the dimensionless time τ required for shock formation versus the initial density fraction c found from Eqs. 6 and 7. Determination of a general exression for the time of shock formation in closed form is difficult; however, solutions for two limiting cases, c. c,an = can be found. B. Shock formation time: c. /2 Consider the region between π/2 and π/2. When c is aroximately equal to /2, the shock must form at a oint to the left of the c axis, but near the oint z =. If the cosine function is aroximated as cos z = z 2 /2 and substituted into Eq. 9 then the rofiles of shown in
5 for oints in the region π/2 z. The modified density fraction c and coordinate z then become z = c cos(φ + τ ) cos 2 φ sin 2 τ sin φ c = c cos(φ + τ (26) ) cos 2 φ sin 2 cos φ. τ The condition for the aearance of a shock is that dz / should aroach zero, which from Eqs. 26 can be found as dz = 5 FIG. 4: Time for formation of a shock τ versus the starting density fraction c from numerical integration of the Hamiltonian system of equations. The oint (N) at c = /2, is the limining value of τ from Eq. 28. The oints ( ) are from Eq.3for c Fig. reduce to circles described by µ 2 z 2 + c 2 = k (2) on the zc lane for k 4. Substitution of z from Eq. 2 into Eq. 7 with the sine function aroximated as sin z = z gives, which can be integrated to give arcsin q c arcsin 4 + k q = ±dτ, (22) 4 + k c 2 c q 4 + k Further maniulation of Eq. 23 gives for z< = τ. (23) 2 c c 2 c 2 + z 2 + c z c 2 + z 2 =sinτ, (24) where z has been defined as z = z/2 2 and τ = τ / 2. It is convenient at this oint to exress the coordinates c and z in Eq. 24 in cylindrical coordinates, c = r cos φ and z = r sin φ, where φ is an angle measured clockwise with resect to the negative c axis Equation 24 gives, after some algebraic maniulation, r = c cos(φ + τ ) cos 2 φ sin 2 (25) τ cos(2φ + τ )[cos 2 φ sin 2 τ ]+sin(2φ)cos(φ + τ )sinφ sin(2φ)cos(φ + τ )sinφ sin(2φ + τ )[cos 2 φ sin 2 τ ] (27) The shock for c aroximately equal to zero is exected to aear at φ = π so that Eq. 27 reduces to dz / = cot τ =;hence, it follows that τ = π/2, for the time of shock generation, and the time for shock formation is τ ' π 2. (28) As can be seen in Fig. 4 the shock time given by Eq. 28 is in excellent agreement with the numerical results. C. Shock formation time: c = A second case where the time of shock formation time can be found in closed form is when the initial density fraction is small. Consider oints in the region π/2 < z<π/2. Starting from Eq. 6, the time for the density fraction to reach a value c starting from the oint (z,c ) canbeshowntobegivenby τ = c c x 2 ( x) 2 c 2 ( c)2 cos 2 z +2 c. x 2 ( x) 2 c 2 ( c ) 2 cos 2 z (29) It is clear from lots found from the solution to Eq. 6 or through numerical integration of Eqs. 6 and 7, that the oints first forming the shock lie originally in the region <z<π/2 and which move downward and to the left along the curves shown in Fig. 2. For small c,theshock formsatvaluesofc ' /2 near π/2 from oints z on the ositive z axis where z '. Equation 29 can then be aroximated as τ ' /2 c x2 c 2 min +2 c x2 c 2 min
6 6 giving τ as τ ' ln c +ln à c +! c 2 c2 min where ' c cos z. for small c is thus Thetimefortheshocktoform τ ' ln c (3) The lot in Fig. 4 gives two oints from Eq. 3. Evaluation of the shock time for values smaller than those given in Fig. 4 shows that the aroximation given by Eq. 3 gets successively better for small values of c. For examle, for c =5 5, numerical integration gives τ =.6 whereas Eq. 3 gives τ =9.9 FIG. 5: Density fraction c versus dimensionless distance z for α=4. The times for the lots are τ =,.2,.4, and.8 for lots with successively higher eak values of c V. Effects of Mass Diffusion A. Results from numerical integration The effects of mass diffusion were determined by numerical integration of Eq. 3 using the finite difference method, which gives a solution by considering small changes in the density fraction resulting from a change in time τ and sace z. A solution for c is obtained by selecting successively smaller values of τ and z until the solution converges. The first order derivative of the density fraction with resect to z at a time t k at the oint z k is given, for instance, by c z = c(t k,z i+ ) c(t k,z i ). (3) 2 z Since derivatives must be comuted at the end oints of the integration range in sace, the eriodicity of the solution is used to give values for the required derivatives. The modulus function, mod(a, b), isintroduced to reresents the remaining integer of division of a by b. For examle, mod(3, 2) is equal to one and mod(4, 2) is equal to zero. If a is less than b, the value of the modulus function is a, where a is a ositive integer, including zero. The exressions for the first and second sace derivatives become c z = c(t k,z mod(i+,imax)) c(t k,z mod(i,imax)) 2 z 2 c z 2 = c(t k,z mod(i+,imax)) 2c(t k,z i )+c(t k,z mod(i,imax)) z 2 To simlify the notation, z i and τ k are now relaced by their indices, giving Eq. 3 as c(k +,i)=c(k, i)+α T τ cos(2π i )( 2c(k, i)) i max (32) c(k, mod(i +,imax )) c(k, mod(i,i max )) 2 z α T τ sin(2π i µ τ )(c(k, i) c 2 (k, i)) + i max z 2 [c(k, mod(i +,i max )) 2c(k, i)+c(k, mod(i,i max ))]. The lots shown in Fig. 5 indicate that the effect of diffusion is to remove the high satial frequencies in the distribution. Since the satial rofile of c deends on α and τ only, it is ossible to calculate the satial Fourier comonents of the distribution after obtaining c from numerical integration, and hence to determine the intensity of light diffracted from an absortion grating, as was done in Ref. [2]. By minimizing the error in a fit to the intensities of the diffracted light beams, the thermal diffusion factor can be determined from exerimental data. B. Density fraction distribution for long times When τ becomes long, thermal diffusion is exactly balanced by mass diffusion and the satial rofile of the density fraction must aroach a limiting form. It can be seen that when /dτ =, Eq. 3 reduces to an ordinary differential equation in z, a first integral of which is given by ˆT + αc( c)d dz dz = C, (33)
7 7 where C is a constant, and a dimensionless satial derivative of the temerature field d ˆT/dz has been substituted for the cosine function. For temerature fields that are eriodic in z it is exected that in the region <z<2π, both /dz and d ˆT/dz will be zero at the some value of z. For the resent roblem, both of these quantities are zero at the oint z =3π/2; hence the constant C must be zero and Eq. 33 can be integrated as c(z) c(z ) c( c) = α The general solution to 34 is thus, z z d ˆT dz. (34) dz c(z) = +Fe αt (z) (35) where F is a constant given by [ c(z )] ex[ α ˆT (z )]/c(z ) For the resent roblem with a sinusoidal temerature field, the long time distribution becomes c(z) =. +Feαsin z (36) An equivalent exression for the linearized Eq. 3 has been given in Ref. [3]. Since the value of c(z ) is not known without solution to Eq. 3, F (α, c ) can be determined for a given value of α and c through use of the mass conservation law, Eq. 5. Deending on the value of α, curves of c versus α from Eq. 36 can resemble a sinusoidal function for small α, or a nearly square wave for large α. VI. Discussion The exressions in Eqs. 2 and 6 rovide an exact solution to the thermal diffusion roblem without mass diffusion. If Eq. 3 is divided by α so that the time derivative is exressed with resect to τ, it is clear that the effects of diffusion are small when α is large or when the gradient of the density fraction is small, such as at short times. The conditions for the time of aearance of the shocks is straightforward, however, the comlexity of an exlicit exression for the shock time in terms of α and c aears to be great. The formation of a shocks follows from the nonlinearity of the deendent variable c in Eq. 3. Perhas the unique feature of shocks roduced by thermal diffusion in sinusoidal temerature fields is that they come as counter-roagating airs of waves, the direction of which is governed by the sign of the cosine function in Eq. 2 to the right and left of z =3π/2.. The left-going wave near z =3π/2 has the roerty that it slows to a seed of zero when c l aroaches one, and c r aroaches zero. At this oint in time, a comlete searation of the mixtures, ossible only in the absence of diffusion, takes lace. The arallel between thermal diffusion shocks and fluid shocks has been given in Ref. [2]. In both cases, the mathematics are treated by ignoring the dissiative term in the equations of motion to find the roerties of the shock waves. In fluid shocks, viscosity causes a broadening of the shock front, just as mass diffusion broadens the thermal diffusion shock fronts as shown by the results of numerical integration given here. The extent to which thermal diffusion shocks are visible in the laboratory is clearly a function of how large the thermal diffusion factor can be made. ACKNOWLEDGEMENT SD and GD acknowledge the suort of the US Deartment of Energy, Office of Basic Energy Studies under grant ER WC acknowledges artial suort for his research from the NSF under grant DMS-728, the NSERC under grant # , and the Canada Research Chairs Program [] C. Ludwig, Sitzber. Akad. Wiss. Vien Math.-naturw. 2, 539 (856). [2]C.Soret,Arch.Sci.Phys.Nat.Geneve3, 48 (879). [3] D.D.Fitts,Nonequilibrium Thermodynamics, McGraw- Hill, New York, 962. [4] S. R. degroot and P. Mazur, Non-Equilibrium Thermodynamics, North Holland, Amsterdam, 962, As noted in this reference the Dufort effect is generally small and can be ignored. [5]E.A.Mason,R.J.Munn,andR.J.Smigh, Thermal diffusion in gases, in Advances in Atomic and Molecular Physics, edited by D. R. Bates, age 33, New York, 966, Academic Press. [6]R.D.Present, Kinetic Theory of Gases, McGraw Hill, New York, 958. [7] K. Thyagarajan and P. Lallemand, Ot. Commun 36, 54 (978). [8] F. Bloisi, Ot. Commun. 68, 87 (988). [9] W. Köhler, J. Chem. Phys. 98, 66 (993). [] W. Köhler and S. Wiegand, Measurement of transort coefficients by an otical grating technique, in Thermal Nonequilibrium Phenomena in Fluid Mixtures, editedby W. Köhler and S. Wiegand, age 9, Sringer Verlag, Berlin, 2. [] H. J. Eichler, P. Gunter, and K. W. Pohl, Laser-Induced Dynamic Gratings, Sringer, Berlin, 985. [2] S. Danworahong, W. Craig, V. Gusev, and G. J. Diebold, submitted for ublication. [3] W. L. Craig, S. Danworahong, and G. J. Diebold, Physical Review Letters 92, 259 (24). [4] L. M. Milne-Thomson, Handbook of Mathematical Functions with Formulas, Grahs, and Mathematical Tables, volume 55, age 596, National Bureau of Standards, Washington D.C., 964, Some authors (e.g. I. Gradshyteyn and Ryzhik, Table of Integrals, Series, and Products, Academic Press, New York, 98) give the ellitic
8 8 integral with the second argument as a/b instead of a 2 /b 2, which we believe to be correct.
SELF-SIMILAR FLOW UNDER THE ACTION OF MONOCHROMATIC RADIATION BEHIND A STRONG CYLINDRICAL SHOCK WAVE IN A NON-IDEAL GAS
SELF-SIMILAR FLOW UNDER THE ACTION OF MONOCHROMATIC RADIATION BEHIND A STRONG CYLINDRICAL SHOCK WAVE IN A NON-IDEAL GAS *J. P. Vishwakarma and Vijay Kumar Pandey Deartment of Mathematics & Statistics,
More informationdn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential
Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system
More informationANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE
ANALYTICAL MODEL FOR THE BYPASS ALE IN A LOOP HEAT PIPE Michel Seetjens & Camilo Rindt Laboratory for Energy Technology Mechanical Engineering Deartment Eindhoven University of Technology The Netherlands
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.
More informationUniversity of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009
University of North Carolina-Charlotte Deartment of Electrical and Comuter Engineering ECG 4143/5195 Electrical Machinery Fall 9 Problem Set 5 Part Due: Friday October 3 Problem 3: Modeling the exerimental
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationChapter 6: Sound Wave Equation
Lecture notes on OPAC0- ntroduction to Acoustics Dr. Eser OLĞAR, 08 Chater 6: Sound Wave Equation. Sound Waves in a medium the wave equation Just like the eriodic motion of the simle harmonic oscillator,
More information16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE
16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE H. Yamasaki, M. Abe and Y. Okuno Graduate School at Nagatsuta, Tokyo Institute of Technology 459, Nagatsuta, Midori-ku, Yokohama,
More informationThickness and refractive index measurements using multiple beam interference fringes (FECO)
Journal of Colloid and Interface Science 264 2003 548 553 Note www.elsevier.com/locate/jcis Thickness and refractive index measurements using multile beam interference fringes FECO Rafael Tadmor, 1 Nianhuan
More informationATM The thermal wind Fall, 2016 Fovell
ATM 316 - The thermal wind Fall, 2016 Fovell Reca and isobaric coordinates We have seen that for the synotic time and sace scales, the three leading terms in the horizontal equations of motion are du dt
More informationFactors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier
Australian Journal of Basic and Alied Sciences, 5(12): 2010-2020, 2011 ISSN 1991-8178 Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doed Fiber Amlifier
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More informationSpin Diffusion and Relaxation in a Nonuniform Magnetic Field.
Sin Diffusion and Relaxation in a Nonuniform Magnetic Field. G.P. Berman, B. M. Chernobrod, V.N. Gorshkov, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 V.I. Tsifrinovich
More informationLow field mobility in Si and GaAs
EE30 - Solid State Electronics Low field mobility in Si and GaAs In doed samles, at low T, ionized imurity scattering dominates: τ( E) ------ -------------- m N D πe 4 ln( + γ ) ------------- + γ γ E 3
More informationPulse Propagation in Optical Fibers using the Moment Method
Pulse Proagation in Otical Fibers using the Moment Method Bruno Miguel Viçoso Gonçalves das Mercês, Instituto Suerior Técnico Abstract The scoe of this aer is to use the semianalytic technique of the Moment
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationPhase transition. Asaf Pe er Background
Phase transition Asaf Pe er 1 November 18, 2013 1. Background A hase is a region of sace, throughout which all hysical roerties (density, magnetization, etc.) of a material (or thermodynamic system) are
More informationAll-fiber Optical Parametric Oscillator
All-fiber Otical Parametric Oscillator Chengao Wang Otical Science and Engineering, Deartment of Physics & Astronomy, University of New Mexico Albuquerque, NM 87131-0001, USA Abstract All-fiber otical
More informationSpeed of sound measurements in liquid Methane at cryogenic temperature and for pressure up to 10 MPa
LNGII - raining Day Delft, August 07 Seed of sound measurements in liquid Methane at cryogenic temerature and for ressure u to 0 MPa Simona Lago*, P. Alberto Giuliano Albo INRiM Istituto Nazionale di Ricerca
More informationPHYS 301 HOMEWORK #9-- SOLUTIONS
PHYS 0 HOMEWORK #9-- SOLUTIONS. We are asked to use Dirichlet' s theorem to determine the value of f (x) as defined below at x = 0, ± /, ± f(x) = 0, - < x
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationChapter 5. Transient Conduction. Islamic Azad University
Chater 5 Transient Conduction Islamic Azad University Karaj Branch 1 Transient Conduction Many heat transfer roblems are time deendent Changes in oerating conditions in a system cause temerature variation
More informationMODULE 2: DIFFUSION LECTURE NO. 2
PTEL Chemical Mass Transfer Oeration MODULE : DIFFUSIO LECTURE O.. STEDY STTE MOLECULR DIFFUSIO I FLUIDS UDER STGT D LMIR FLOW CODITIOS.. Steady state diffusion through a constant area Steady state diffusion
More informationHighly improved convergence of the coupled-wave method for TM polarization
. Lalanne and G. M. Morris Vol. 13, No. 4/Aril 1996/J. Ot. Soc. Am. A 779 Highly imroved convergence of the couled-wave method for TM olarization hilie Lalanne Institut d Otique Théorique et Aliquée, Centre
More informationAdiabatic Shear Bands in Simple and Dipolar Plastic Materials
Adiabatic Shear Bands in Simle and Diolar Plastic Materials T W \-1RIGHT us Army Ballistic Research Laboratory Aberdeen Proving Ground, MD 215 R C BATRA University of Missouri-Rolla Rolla, Missouri 6541
More informationOptical Fibres - Dispersion Part 1
ECE 455 Lecture 05 1 Otical Fibres - Disersion Part 1 Stavros Iezekiel Deartment of Electrical and Comuter Engineering University of Cyrus HMY 445 Lecture 05 Fall Semester 016 ECE 455 Lecture 05 Otical
More informationONE. The Earth-atmosphere system CHAPTER
CHAPTER ONE The Earth-atmoshere system 1.1 INTRODUCTION The Earth s atmoshere is the gaseous enveloe surrounding the lanet. Like other lanetary atmosheres, it figures centrally in transfers of energy between
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
More informationThe Second Law: The Machinery
The Second Law: The Machinery Chater 5 of Atkins: The Second Law: The Concets Sections 3.7-3.9 8th Ed, 3.3 9th Ed; 3.4 10 Ed.; 3E 11th Ed. Combining First and Second Laws Proerties of the Internal Energy
More informationarxiv: v1 [nucl-th] 26 Aug 2011
The Viscosity of Quark-Gluon Plasma at RHIC and the LHC Ulrich Heinz, Chun Shen and Huichao Song Deartment of Physics, The Ohio State University, Columbus, Ohio 436, USA Lawrence Berkeley National Laboratory,
More informationA review on Lamb's atmospheric oscillations using initial value problem approach
Journal of Physics: Conference Series OPEN ACCESS A review on Lamb's atmosheric oscillations using initial value roblem aroach To cite this article: Ángel De Andrea González 04 J. Phys.: Conf. Ser. 56
More informationWolfgang POESSNECKER and Ulrich GROSS*
Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND
More informationMET 4302 Midterm Study Guide 19FEB18
The exam will be 4% short answer and the remainder (6%) longer (1- aragrahs) answer roblems and mathematical derivations. The second section will consists of 6 questions worth 15 oints each. Answer 4.
More informationSetting up the Mathematical Model Review of Heat & Material Balances
Setting u the Mathematical Model Review of Heat & Material Balances Toic Summary... Introduction... Conservation Equations... 3 Use of Intrinsic Variables... 4 Well-Mixed Systems... 4 Conservation of Total
More informationJohn Weatherwax. Analysis of Parallel Depth First Search Algorithms
Sulementary Discussions and Solutions to Selected Problems in: Introduction to Parallel Comuting by Viin Kumar, Ananth Grama, Anshul Guta, & George Karyis John Weatherwax Chater 8 Analysis of Parallel
More informationε(ω,k) =1 ω = ω'+kv (5) ω'= e2 n 2 < 0, where f is the particle distribution function and v p f v p = 0 then f v = 0. For a real f (v) v ω (kv T
High High Power Power Laser Laser Programme Programme Theory Theory and Comutation and Asects of electron acoustic wave hysics in laser backscatter N J Sircombe, T D Arber Deartment of Physics, University
More informationInternational Journal of Mathematics Trends and Technology- Volume3 Issue4-2012
Effect of Hall current on Unsteady Flow of a Dusty Conducting Fluid through orous medium between Parallel Porous Plates with Temerature Deendent Viscosity and Thermal Radiation Harshbardhan Singh and Dr.
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationCentral Force Motion Challenge Problems
Central Force Motion Challenge Problems Problem 1: Ellitic Orbit A satellite of mass m s is in an ellitical orbit around a lanet of mass m which is located at one focus of the ellise. The satellite has
More informationChapter 2 Introductory Concepts of Wave Propagation Analysis in Structures
Chater 2 Introductory Concets of Wave Proagation Analysis in Structures Wave roagation is a transient dynamic henomenon resulting from short duration loading. Such transient loadings have high frequency
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationLight at a Standstill Tim Kuis June 13, 2008
Light at a Standstill Tim Kuis June 13, 008 1. Introduction There is something curious about the seed of light. It is the highest obtainable seed. Nothing can travel faster. But how slow can light go?
More informationOptical self-energy of superconducting Pb in the terahertz region
Otical self-energy of suerconducting Pb in the terahertz region T. Mori, 1 E. J. Nicol, 2, * S. Shiizuka, 1 K. Kuniyasu, 3 T. Nojima, 3 N. Toyota, 1 and J. P. Carbotte 4 1 Physics Deartment, Graduate School
More informationDomain Dynamics in a Ferroelastic Epilayer on a Paraelastic Substrate
Y. F. Gao Z. Suo Mechanical and Aerosace Engineering Deartment and Princeton Materials Institute, Princeton University, Princeton, NJ 08544 Domain Dynamics in a Ferroelastic Eilayer on a Paraelastic Substrate
More informationAnalysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia
roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.
More informationNONRELATIVISTIC STRONG-FIELD APPROXIMATION (SFA)
NONRELATIVISTIC STRONG-FIELD APPROXIMATION (SFA) Note: SFA will automatically be taken to mean Coulomb gauge (relativistic or non-diole) or VG (nonrelativistic, diole-aroximation). If LG is intended (rarely),
More informationFUGACITY. It is simply a measure of molar Gibbs energy of a real gas.
FUGACITY It is simly a measure of molar Gibbs energy of a real gas. Modifying the simle equation for the chemical otential of an ideal gas by introducing the concet of a fugacity (f). The fugacity is an
More informationREFLECTION AND TRANSMISSION BAND STRUCTURES OF A ONE-DIMENSIONAL PERIODIC SYSTEM IN THE PRESENCE OF ABSORPTION
Armenian Journal of Physics, 0, vol. 4, issue,. 90-0 REFLECTIO AD TRASMISSIO BAD STRUCTURES OF A OE-DIMESIOAL PERIODIC SYSTEM I THE PRESECE OF ABSORPTIO A. Zh. Khachatrian State Engineering University
More informationHypersonic flow: introduction
Hyersonic flow: introduction Van Dyke: Hyersonic flow is flow ast a body at high ach number, where nonlinearity is an essential feature of the flow. Also understood, for thin bodies, that if is the thickness-to-chord
More informationThe individual electric and magnetic waves are in phase. The fields peak at the same position at the same time.
1 Part 3: Otics 3.1: Electromagnetic Waves An electromagnetic wave (light wave) consists of oscillating electric and magnetic fields. The directions of the electric and magnetic fields are erendicular.
More informationR g. o p2. Lecture 2: Buoyancy, stability, convection and gravity waves
Lecture : Clarifications of lecture 1: Hydrostatic balance: Under static conditions, only gravity will work on the fluid. Why doesn't all the fluid contract to the ground? Pressure builds u and resists
More informationBrownian Motion and Random Prime Factorization
Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationSimplifications to Conservation Equations
Chater 5 Simlifications to Conservation Equations 5.1 Steady Flow If fluid roerties at a oint in a field do not change with time, then they are a function of sace only. They are reresented by: ϕ = ϕq 1,
More informationOn the relationship between sound intensity and wave impedance
Buenos Aires 5 to 9 Setember, 16 Acoustics for the 1 st Century PROCEEDINGS of the nd International Congress on Acoustics Sound Intensity and Inverse Methods in Acoustics: Paer ICA16-198 On the relationshi
More informationSection 4: Electromagnetic Waves 2
Frequency deendence and dielectric constant Section 4: Electromagnetic Waves We now consider frequency deendence of electromagnetic waves roagating in a dielectric medium. As efore we suose that the medium
More informationLiquid water static energy page 1/8
Liquid water static energy age 1/8 1) Thermodynamics It s a good idea to work with thermodynamic variables that are conserved under a known set of conditions, since they can act as assive tracers and rovide
More informationKeywords: pile, liquefaction, lateral spreading, analysis ABSTRACT
Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationMath 121: Calculus 1 - Fall 2012/2013 Review of Precalculus Concepts
Introduction Math : Calculus - Fall 0/0 Review of Precalculus Concets Welcome to Math - Calculus, Fall 0/0! This roblems in this acket are designed to hel you review the toics from Algebra and Precalculus
More informationTemperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures
Indian Journal of Pure & Alied Physics Vol. 44, December 2006,. 953-958 Temerature, current and doing deendence of non-ideality factor for n and nn unch-through structures Khurshed Ahmad Shah & S S Islam
More informationBulletin of the Seismological Society of America
Bulletin of the Seismological Society of America Vol. 69 August 1979 No. 4 ELEMENTARY SOLUTIONS TO LAMB'S PROBLEM FOR A POINT SOURCE AND THEIR RELEVANCE TO THREE-DIMENSIONAL STUDIES OF SPONTANEOUS CRACK
More informationRANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES
RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationI Poles & zeros. I First-order systems. I Second-order systems. I E ect of additional poles. I E ect of zeros. I E ect of nonlinearities
EE C28 / ME C34 Lecture Chater 4 Time Resonse Alexandre Bayen Deartment of Electrical Engineering & Comuter Science University of California Berkeley Lecture abstract Toics covered in this resentation
More informationMODELING OF CONSTANT VOLUME COMBUSTION OF PROPELLANTS FOR ARTILLERY WEAPONS
MODELING OF CONSTANT VOLUME COMBUSTION OF PROPELLANTS FOR ARTILLERY WEAPONS Lt.col. lect. dr.ing. Daniel ANTONIE, Col.(r). rof.as. dr.ing. Sorin GHEORGHIAN, Col.(r). C.S.II dr.ing. Nicolae MĂRUNȚELU MILITARY
More informationCompressible Flow Introduction. Afshin J. Ghajar
36 Comressible Flow Afshin J. Ghajar Oklahoma State University 36. Introduction...36-36. he Mach Number and Flow Regimes...36-36.3 Ideal Gas Relations...36-36.4 Isentroic Flow Relations...36-4 36.5 Stagnation
More informationNanoscale fluctuations and surface tension measurements in droplets using phase-resolved low-coherence interferometry
Nanoscale fluctuations and surface tension measurements in drolets using hase-resolved low-coherence interferometry Ru Wang, 1 Taewoo Kim, 2 Mustafa Mir, 2 and Gabriel Poescu 2, * 1 Quantitative Light
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationSimple geometric interpretation of signal evolution in phase-sensitive fibre optic parametric amplifier
Simle geometric interretation of signal evolution in hase-sensitive fibre otic arametric amlifier A.A. REDYUK,,,* A.E. BEDNYAKOVA,, S.B. MEDVEDEV, M.P. FEDORUK,, AND S.K. TURITSYN,3 Novosibirsk State University,
More informationOn Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada
On Using FASTEM2 for the Secial Sensor Microwave Imager (SSM/I) March 15, 2001 Godelieve Deblonde Meteorological Service of Canada 1 1. Introduction Fastem2 is a fast model (multile-linear regression model)
More informationAutomatic Generation and Integration of Equations of Motion for Linked Mechanical Systems
Automatic Generation and Integration of Equations of Motion for Linked Mechanical Systems D. Todd Griffith a, John L. Junkins a, and James D. Turner b a Deartment of Aerosace Engineering, Texas A&M University,
More information3.4 Design Methods for Fractional Delay Allpass Filters
Chater 3. Fractional Delay Filters 15 3.4 Design Methods for Fractional Delay Allass Filters Above we have studied the design of FIR filters for fractional delay aroximation. ow we show how recursive or
More informationTHE FIRST LAW OF THERMODYNAMICS
THE FIRST LA OF THERMODYNAMIS 9 9 (a) IDENTIFY and SET UP: The ressure is constant and the volume increases (b) = d Figure 9 Since is constant, = d = ( ) The -diagram is sketched in Figure 9 The roblem
More informationChemical Kinetics and Equilibrium - An Overview - Key
Chemical Kinetics and Equilibrium - An Overview - Key The following questions are designed to give you an overview of the toics of chemical kinetics and chemical equilibrium. Although not comrehensive,
More informationNumerical Modeling of Powder Flow during Coaxial Laser Direct Metal Deposition Comparison between Ti-6Al-4V Alloy and Stainless Steel 316L
Numerical Modeling of Powder Flow during Coaxial Laser Direct Metal Deosition Comarison between Ti-6Al-4V Alloy and Stainless Steel 316L S. Morville 1, M. Carin *1, D. Carron 1, P. Le Masson 1, M. Gharbi,
More informationSUMS OF TWO SQUARES PAIR CORRELATION & DISTRIBUTION IN SHORT INTERVALS
SUMS OF TWO SQUARES PAIR CORRELATION & DISTRIBUTION IN SHORT INTERVALS YOTAM SMILANSKY Abstract. In this work we show that based on a conjecture for the air correlation of integers reresentable as sums
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More informationAtmosphere, Ocean and Climate Dynamics Answers to Chapter 4
Atmoshere, Ocean and Climate Dynamics Answers to Chater 4 1. Show that the buoyancy frequency, Eq.(4.22), may be written in terms of the environmental temerature rofile thus N 2 = g µ dte T E dz + Γ d
More informationA Closed-Form Solution to the Minimum V 2
Celestial Mechanics and Dynamical Astronomy manuscrit No. (will be inserted by the editor) Martín Avendaño Daniele Mortari A Closed-Form Solution to the Minimum V tot Lambert s Problem Received: Month
More informationON POLYNOMIAL SELECTION FOR THE GENERAL NUMBER FIELD SIEVE
MATHEMATICS OF COMPUTATIO Volume 75, umber 256, October 26, Pages 237 247 S 25-5718(6)187-9 Article electronically ublished on June 28, 26 O POLYOMIAL SELECTIO FOR THE GEERAL UMBER FIELD SIEVE THORSTE
More informationKinetics of Protein Adsorption and Desorption on Surfaces with Grafted Polymers
1516 Biohysical Journal Volume 89 Setember 2005 1516 1533 Kinetics of Protein Adsortion and Desortion on Surfaces with Grafted Polymers Fang Fang,* Javier Satulovsky, y and Igal Szleifer* *Deartment of
More informationWave Drift Force in a Two-Layer Fluid of Finite Depth
Wave Drift Force in a Two-Layer Fluid of Finite Deth Masashi Kashiwagi Research Institute for Alied Mechanics, Kyushu University, Jaan Abstract Based on the momentum and energy conservation rinciles, a
More information5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:
5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum
More informationwhether a process will be spontaneous, it is necessary to know the entropy change in both the
93 Lecture 16 he entroy is a lovely function because it is all we need to know in order to redict whether a rocess will be sontaneous. However, it is often inconvenient to use, because to redict whether
More informationDufour-Driven Thermosolutal Convection of the Veronis Type
JOURNAL OF MATEMATICAL ANALYI AND APPLICATION 8, 56958 998 ARTICLE AY975753 Dufour-Driven Thermosolutal Convection of the Veronis Tye ari Mohan Deartment of Mathematics, imachal Pradesh Uniersity, ummer
More informationOn Optimization of Power Coefficient of HAWT
Journal of Power and Energy Engineering, 14,, 198- Published Online Aril 14 in Scies htt://wwwscirorg/journal/jee htt://dxdoiorg/1436/jee1448 On Otimization of Power Coefficient of HAWT Marat Z Dosaev
More informationModeling Volume Changes in Porous Electrodes
Journal of The Electrochemical Society, 53 A79-A86 2006 003-465/2005/53/A79/8/$20.00 The Electrochemical Society, Inc. Modeling olume Changes in Porous Electrodes Parthasarathy M. Gomadam*,a,z John W.
More informationDimensional perturbation theory for Regge poles
Dimensional erturbation theory for Regge oles Timothy C. Germann Deartment of Chemistry, University of California, Berkeley, California 94720 Sabre Kais Deartment of Chemistry, Purdue University, West
More informationON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS
#A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,
More informationA Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands
A Secial Case Solution to the Persective -Point Problem William J. Wolfe California State University Channel Islands william.wolfe@csuci.edu Abstract In this aer we address a secial case of the ersective
More informationUnsteady Flow of a Dusty Conducting Fluid through porous medium between Parallel Porous Plates with Temperature Dependent Viscosity and Heat Source
Volume Issue3 3- June www.ijsret.org ISSN 78-88 Unsteady Flow of a Dusty Conducting Fluid through orous medium between Parallel Porous Plates with Temerature Deendent Viscosity and Heat Source Shalini
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationThe University of the State of New York REGENTS HIGH SCHOOL EXAMINATION COURSE III. Wednesday, August 16, :30 to 11:30 a.m.
The University of the State of New York REGENTS HIGH SCHOOL EXAMINATION THREE-YEAR SEQUENCE FOR HIGH SCHOOL MATHEMATICS COURSE III Wednesday, August 6, 000 8:0 to :0 a.m., only Notice... Scientific calculators
More informationMontgomery self-imaging effect using computer-generated diffractive optical elements
Otics Communications 225 (2003) 13 17 www.elsevier.com/locate/otcom Montgomery self-imaging effect using comuter-generated diffractive otical elements J urgen Jahns a, *, Hans Knuertz a, Adolf W. Lohmann
More informationLECTURE 3 BASIC QUANTUM THEORY
LECTURE 3 BASIC QUANTUM THEORY Matter waves and the wave function In 194 De Broglie roosed that all matter has a wavelength and exhibits wave like behavior. He roosed that the wavelength of a article of
More informationEmittance Growth Caused by Surface Roughness
Emittance Growth Caused by Surface Roughness he hang, Chuanxiang Tang Tsinghua University, Beijing Oct. 17th, 2016 Motivation What causes the emittance growth Dowell s equations of QE & emittance for bulk
More informationREFINED STRAIN ENERGY OF THE SHELL
REFINED STRAIN ENERGY OF THE SHELL Ryszard A. Walentyński Deartment of Building Structures Theory, Silesian University of Technology, Gliwice, PL44-11, Poland ABSTRACT The aer rovides information on evaluation
More information