A Thermomechanical Model of Gels
|
|
- Adam Mason
- 5 years ago
- Views:
Transcription
1 A Thermomechanical Model of Gels A DISSERTATION SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Minsu Kim IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY Yoichiro Mori, Advisor September, 04
2 c Minsu Kim 04 ALL RIGHTS RESERVED
3 Acknowledgements First of all, I would like to thank Professor Yoichiro Mori for the direction to complete this thesis as well as for his unwavering support. His sincere guidance has influenced me in many aspects: from mutual trust to the research commitment and to teaching philosophy. I feel fortunate to have studied with him during my graduate course. Also, I would like to thank Professor Maria-Carme Calderer for her generosity and for the invaluable support in the beginning of my research. Broadly-based references she provided helped me have a variety of views. I also appreciate Summer Research Assistantships she offered consistently. In particular, I want to thank Professor Richard McGehee for his encouragement when I experienced frustration. It would be impossible to have come this far without his motivation and understanding. I am extremely grateful to my family and husband for their unconditional and emotional support throughout this long journey. Particularly, I am much obliged to my sister, Yeon-Kyoung for her steadfast support. Finally, I deeply thank God for His guidance and protection always. i
4 Dedication To my husband, Jaebeom and my sister, Yeon-Kyoung. ii
5 Abstract In this thesis, a model of polymer gels on thermal effects is developed and analyzed, and numerical simulations for a one-dimensional geometry is performed. Gels are useful materials in many manufacturing as well as biomedicine, by use of which a variety of apparatus such as drug delivery devices are designed. It is important in physics to observe how gels change under different temperature, and in particular, we examine how gels respond to the gradient of temperature. Even though temperature is a key factor to be considered, no physical theory with thermal gradients has been presented. That becomes a motivation for this work, and we construct a thermodynamic theory of gels and apply it to investigate several phenomena in the one-dimensional geometry, such as the evolution of the interface between the gel and the environmental fluid. iii
6 Contents Acknowledgements Dedication Abstract List of Tables List of Figures i ii iii vi vii Introduction Preliminaries 6. Basic Principles of Continuum Mechanics Basic Laws of Thermodynamics Mechanics of Gels 3. A Mechanical Model of Gels A Theory of Thermal Effects in Gels 5 4. Energy Balance Boundary Conditions Deduced by Thermodynamics Second Law Evolution Equations for Temperature Entropy Inequality and the Resulting Relations Galilean Invariance iv
7 5 Gels in One-dimensional Geometry Governing Equations in One-Dimensional Space Referential Form for the Equations of Gels Numerical Algorithmn Discretization Algorithm Convergence Studies Computation Results 5 6. Nondimensionalization Analytic Estimates and Simulations Conclusion 59 Bibliography 6 v
8 List of Tables 5. Triples of h, t, N Divided Differences with M 0 = vi
9 List of Figures 3. System Setup Positions of Two Interfaces of the Moving Gel to the Right Positions of Two Interfaces of the Moving Gel to the Left vii
10 Chapter Introduction Gels are cross-linked polymer networks immersed but not dissolved in a solvent. The cross-links in polymer networks cause solid-like behavior of gels, but addition of solvent helps gels to become softer and swollen by degrees. Gels present itself in nature in animal and plant tissues. Gels are also found in a wide range of applications from foods to many manufacturing devices like contact lenses, hair gel and paints, to medical devices such as drug delivery devices or implants. In this thesis, we develop and analyze a model of polymer gels driven by thermal effects. We have incorporated thermal effects into the equations of gel dynamics. This is particularly relevant for implanted biomedical devices. After the device is implanted into the body, it should adapt to the body environment to work stably in the body. One of the factors that influence the device is temperature, and it is important to consider how the device responds to temperature change. Mathematical modeling and analysis of thermodynamics of gels is thus significant to the biomedicine. Nonetheless, no physical theory for such system was available. That becomes a motivation for this work, and we take this system into account and develop thermal effects on it. We assume that the gel is non-ionic, immiscible, and incompressible. Due to the assumption of non-ionic gel, electrical effects can be neglected. The immiscibility means that the responses in the mixture are depending on volume fractions. On the contrary to this, we call a mixture miscible if corresponding constitutive equations are independent of the volume fractions. Similarly, the incompressibility characterizes a mixture as follows. For a single body, an incompressible one is refered to a body in which its
11 divergence of velocity is free. Mass balance yields zero material time derivative of density, that is, the mass density is constant. In the case of a mixture, we distinguish two kinds of densities, the bulk density and the intrinsic or true density for each component. ρ i denotes the bulk density of the ith component which stands for the mass of the ith constituent per unit volume of mixture. On the other hand, the intrinsic density of the ith component denoted by γ i is its mass per unit true volume which the ith constituent occupies by itself. These two densities are intimately related by the volume fraction φ i as ρ i = γ i φ i, i =,. The mixture is said to be incompressible if the intrinsic density is independent of time for each constituent in the mixture. For convenience sake, in this paper, we identically set the intrinsic densities to be one, and deal with volume fractions for polymer and solvent in the gel instead of bulk densities. In addition, we assume that inertial effects are neglected in the system to characterize its slow motion. The incompressibility assumption on the mixture is mathematically set as each of the components having constant mass density when found in the single phase. This, together with the saturation assumption, that is φ + φ =, gives the relation v φ + v φ = 0. We point out that, in some gel literature, this is known as the incompressibility condition of the mixture. Here, v i represents the velocity field of the individual components. We postulate laws of balance of mass and linear momentum for each of the components, and the law of balance of energy for the whole system. Under constant temperature, the free energy density ψ of the gel, which characterizes the mechanical behavior, is often assumed to consist of two energy functions, the stored elastic energy function of the polymer network and the Flory-Huggins energy of mixing. However, inclusion of thermal energy is necessary when temperature is not constant. Since thermal effects are present, the internal energy density ε is also taken into consideration in order to account for the non-recoverable part of the energy due to heat and temperature as well as for the mechanical contribution. The thermodynamic relation ε = ψ + θη, where θ > 0 denotes absolute temperature, defines the entropy densityη of the system. The remaining constitutive equations include the stress tensor of each component, consisting of the sum of elastic and viscous parts, the pressure, the heat flux vector, and the drag forces between fluid and solid in the gel. Pressure is introduced as a Lagrange
12 3 multiplier associated with the incompressibility assumption on the mixture. We assume that the system, comprising the gel and the surrounding fluid, is thermally isolated. The Second Law of Thermodynamics states that the total entropy is non-decreasing in such a thermally insulated system. We follow the procedure developed in [] and found as well in [] and [3] that uses the Second Law of Thermodynamics to derive restrictions on the constitutive equations of the material. We point out that the approach to derive the constitutive equations is not a trivial generalization of the case with a single phase. In particular, it takes special care to define the total internal energy density of the mixture. The approach presented here is different from previous works on multi-component mixtures with thermal effects [5], [6] which define the total energy as the sum of energies of components in the mixture. Indeed, measuring the energy values of each constituent separately is almost impossible. Instead, the total energy of mixture is given as a whole in most cases. For this reason, in this thesis, we take a different point of view for the description of energy. It is, in fact, the other way around. Given representation of free energy density per unit volume for the whole mixture, we derive both the internal energy density and entropy density for the gel from thermodynamic relations, and at the same time those quantities associated with the individual components are derived by the thermodynamic property of energy that the energy function be extensive. In the system concerned here, the gel is entirely submerged in its own fluid without touching the boundary of the system and fluid moves through the gel freely. This is a special feature of the gel model presented here; the boundary conditions follow from postulating balance laws across the interface between the gel and its surrounding fluid. This approach was first developed in [3], for the purely mechanical case as well as the case when families of ions are present in the system. In all such situations, the problem becomes a system of partial differential equations with a free boundary that being the interface between gel and fluid, assumed to move with the normal component of the polymer velocity at that boundary. Before giving the outline of this dissertation, we summarize some related results concerning mixture theory. The classical book [4] deals extensibly with mathematical models of mixture theories, in abstract physical frameworks. Although such a reference is an outstanding resource to assist on the studies of mixture theories, none of the
13 4 models there was directly relevant to our problems. The paper [3], on the dynamics of polyelectrolyte gels, without including thermal effects, starts with postulating mechanical equations for each component as well as the balance of energy of the system. Proper boundary conditions at the interface between the gel and the outer fluid are postulated, as well as the diffusion equations of ion species. Then dissipation laws of the energy that include the Flory-Huggins mixing energy, the polymer elastic energy and the ionic and electrostatic energies are shown. These give rise to a governing system of equations that it is subsequently studied in the one dimensional case. A mathematical model of dynamics of gels is studied by B. Chabaud and M. C. Calderer [8], with focus on the immiscible and incompressible gel mechanics, neglecting ionic and thermal effects. The article analyzes well-posedness of the coupled inertialess elasticity-stokes problem, with boundary conditions ranging from impermeable to fully permeable. In this article, residual stresses are taken into account, which poses a challenge to the linear stability analysis. Further, they consider two different kinds of solvent type: inviscid and viscous. In [6], Dunwoody states basic balance equations for a mixture that consists of constituents more than two with different temperature postulating a general entropy inequality. Such governing equations are subsequently applied to a binary solid-fluid mixture at a single temperature under the assumption that there are no chemical reactions between components. However, a significant difference with the present approach is in that they postulate the equation of balance of energy for the gel, assuming that its boundary moves with the average velocity of the components. However, this does not seem to be a correct assumption for the present model, in part because the gel boundary is defined to be exactly the polymer boundary, and so it does move with the polymer velocity. Kenyon, in [7], concerns a mixture of N-fluid like constituents and a porous, deformable solid N having the same absolute temperature. In the article, chemical reactions are excluded, which results in zero growth of mass of each species. Balance equations and constitutive theory are also obtained by following the approach by Coleman & Noll. Finally, the equilibrium results are discussed, as usual. In [], an application of the mixture theory to a two-phase flow is presented. The speed of boundary evolution is assumed to be that of either one of the two fluids, and so, it is also different from the assumption in the present work. This thesis has the following organization. Chapter provides a brief review of
14 5 fundamental principles of Continuum Mechanics and basic laws of Thermodynamics. In section., equations of mass conservation and balance of linear momentum for a single component elastic matter are presented. Section. discusses the equation of balance of energy and illustrates the procedure of derivation of constitutive equations and dissipation inequalities from the Second Law of Thermodynamics. Chapter 3 describes the system under consideration and presents a mechanical model of the gel as well as conditions for the interface between the gel and the surrounding fluid. In Chapter 4, we study the thermodynamics of gels. Based upon the first and the second laws of thermodynamics, in particular, section 4. and section 4. establish energy conservation equations and additional boundary conditions in the system. Section 4.3 presents a splitting of the energy of the whole gel into energies of the individual elements, for which the extensive property of energy is used, and in consequence, derives time-dependent equations for temperatures. Thermodynamics constitutive restrictions are obtained in section 4.4. Lastly, section 4.5 proves that the resulting equations for our system are invariant in Galilean relativity. The outcome of this chapter is the formulation of the initial boundary value problem for the evolution of the gel in an arbitrary domain submerged into the fluid. Chapter 5 contains the governing equations of the gel in one dimensional geometry and a numerical scheme for solving the system. Section 5. provides the equations in one spatial dimension, and they are expressed in the referential forms in section 5.. Subsequently, a numerical algorithm is proposed and convergence studies are carried out in section 5.3. Chapter 6 demonstrates computational results. In section 6., we determine relatively crucial parameters in the motion of the gel on thermal effects by nondimensionalization of the governing equations. In section 6., we investigate conditions under which boundaries of the gel moves in the same direction by analytical computations together with numerical simulations. Chapter 7, the final chapter of this thesis, presents conclusion as well as a summary of the work.
15 Chapter Preliminaries This chapter gives a review of basic laws in continuum mechanics and thermodynamics.. Basic Principles of Continuum Mechanics We introduce some of the concepts used in this dissertation for the case of a single elastic material, occupying a bounded reference domain Ω R 3, with smooth boundary Ω. Let Ω t be the image of Ω under the time-dependent deformation map x = ϕx, t, X Ω. The deformation gradient is denoted by F = ϕx, t, and we assume that F is continuous and it satisfies the invertibility condition, det F > 0. Remark. Given the invertibility assumption on the deformation map, either material or spatial representations of a function, gx, t = gϕx, t, t = GX, t, is considered as needed. Let B Ω denote any sub-body of the body, and B t the deformed body. We write ρ = ρx, t for the mass density at the spatial point x in the current configuration. The 6
16 balance equations of mass and linear momentum are then the following: d ρx, t dx = ρv n da, dt B t B t d ρx, tvx, t dx = ρvv n da + T n da + b dx, dt B t B t B t B t 7 where n denotes the outward unit normal vector from the boundary, v the velocity field, T the Cauchy stress tensor, b the body force, and B t represents the boundary of B t. Since B t is arbitrary, we can formally express the previous balance equations in the local form ρ + ρv = 0,. t v ρ + v v = T + b.. t. Basic Laws of Thermodynamics This section provides the first two laws of thermodynamics and constitutive restrictions. The first law of thermodynamics is a statement of energy conservation and the second law states nonnegative production of entropy in a thermodynamic system. Here, we follow the approaches developed in [3] and [] on the original work on reduction of constitutive equations by means of the Second Law of Thermodynamics []. We associate with each dynamical process the specific internal energy ε = εx, t, the specific entropy η = ηx, t and the specific free energy ψ = ψx, t. The net internal energy and entropy of B t, respectively, are given by EB t = ρεx, t dx, B t SB t = ρηx, t dx. B t Note. The term specific indicates the one measured per unit mass. The heat entering the part B t per unit time is defined as Q = q n da + B t ρr dx, B t
17 where q denotes the heat flux vector and r the heat supply density heat per unit mass entering the body per unit time. The equation of energy balance is then d ρε dx = ερv + q n da + T v n da + ρr dx, dt B t B t B t B t and the arbitrariness of B t together with absence of heat supply leads to the local form ε ρ + v ε = T : v q..3 t 8 Now let us denote the material time derivative of a function with a superimposed dot. That is, fx, t := f t + f dx x dt. The production of entropy per unit time in the part B t is defined by Γ := Ṡ + q B t θ n da ρr dx..4 B t θ The Second Law of Thermodynamics states that Γ 0, for all B t Ω t. The local form of this inequality, combined with equations.-.3 and Reynolds transport theorem, becomes where D = v + vt. ρθ η ε + T : D q θ 0,.5 θ The following constitutive equations define a thermoelastic material: where g = θ. T = T θ, g, F, F, η = ηθ, g, F, F, ε = εθ, g, F, F, q = qθ, g, F, F,
18 By using the chain rule, we can write the material time derivative for η: and, likewise, for ε. Here, η = η θ θ + η g ġ + η F : F + η F : F, F denotes the material time derivative of F. Substituting the resulting expressions into.5, and gathering analogous terms yield ρ θ η θ ε θ + ρ θ η θ g ε ġ + T : D g + ρ θ η F ε : F + ρ θ η F F ε F : F q g 0..6 θ It can be verified by appropriate choice of body force and heat radiation term r that at any given place and time, it is possible to find a deformation, stress, and heat flux satisfying the balance laws of momentum and energy and giving to the quantities θ, g, F, θ, ġ, F, F arbitrarily assigned values. According to the Clausius-Duhem inequality, the left hand side of.6 must be non-negative, identically in these six quantities. This implies θ η θ ε θ = 0, θ η g ε g = 0, Using the thermodynamic relation the first equation in.7 gives θ η F ε F = 0,.7 ρ θ η F ε : F + T : D 0, F q θ g 0. ε = ψ + θη, η = ψ θ, which together with the second and third equations yields ψ = ψθ, F, η = ηθ, F, ε = εθ, F..8 9
19 We now assume that the Cauchy stress tensor can be written as the sum of an elastic and a dissipative contribution T = T elas + T visc. The fourth relation in.7 together with the chain rule expression and equation.8 gives Further, we get the residual inequality F = vf, T elas = ρ ψ F F T. T visc : D 0. In the case that Newtonian dissipation is assumed, we get T visc = νd, where ν > 0 denotes the Newtonian dissipation. 0 Likewise, a sufficient condition to satisfy the last residual inequality in.7 is known as the Fourier law of heat, that is, where λ is the thermal conductivity. q = λ θ, λ > 0,
20 Chapter 3 Mechanics of Gels 3. A Mechanical Model of Gels The system in our consideration consists of a gel and solvent. In the system, the gel is fully submerged under the solvent and assumed to be an incompressible and immiscible mixture of two components, elastic polymer and the same solvent as surroundings. Let the system be fixed and we denote it by U. The domain, denoted by Ω, where polymer networks occupy in the reference configuration is mapped to a region, represented by Ω t, after time t by a motion x = ϕ t X with det X ϕ t X > 0, which is a smooth one to one function that assigns to each point X Ω a spatial point x at time t. Accordingly, we denote Γ t = Ω t for the boundary of the gel and R t = U \ Ω t Γ t for the surrounding fluid region. We assume that Γ t does not touch U for all time. Let v,,f be the velocities of the polymer and fluid in Ω t, and the fluid in R t, respectively. For simplicity of notation, we use subscripts i =,, f in order to refer to polymer, solvent in Ω t, and the outer fluid, respectively. The velocity of polymer v is defined by v x, t ϕ tx. t
21 Figure 3.: System Setup The system consists of the gelω t and the surrounding fluidr t. φ, v, p, F, θ represent volume fraction, velocity, pressure, the deformation gradient, and temperature, respectively. Subscripts,, and f stand for polymer networks, fluid in and out of the gel in order. Let F = X ϕ t X be the deformation gradient. Then, for the deformation gradient evaluated in Ω t, F = F ϕ t, it follows that F t + v F = v F. 3. We regard that each current position x Ω t holds both polymer and fluid as their portion. The volume fractions of polymer networks and solvent in the gel are denoted by φ and φ, respectively, and they satisfy φ + φ =. 3. It is assumed that the two components in the gel have common temperature θ x, t = θ x, t = θx, t > 0 at a point x in Ω t, and it has the same value on the boundary of the gel as θ f x, t, temperature at x in R t. Also, we assume that the temperature on U keeps constant for all time. By the incompressibility, the balance of mass for each constituent in the gel can be
22 written in terms of the volume fraction 3 φ i t + v iφ i = 0, i =,, 3.3 with which the saturation assumptionφ + φ = implies φ v + φ v = The above 3.4 is sometimes called the incompressibility condition of mixture. Since we neglect inertial effects, the balance of linear momentum for ith component in Ω t reduces to T i φ i p + f i = 0, f + f = 0, 3.5 where T i is the Cauchy stress exerting on ith constituent, p is the pressure that enforces the incompressibility condition 3.4 and f i is the momentum supply on the component i due to the local interaction with the other component. We have assumed that there is no external body force. We note that the local mass balance 3.3 can be referentially written as where φ ref stands for the volume fraction φ at time t = 0. φ detf = φ ref, 3.6 Similarly, under the assumptions that the fluid is also incompressible and that no body force occurs in R t, we have T f p f = 0, 3.7 v f = We suppose that stresses T, T, T f are given by the following: T = T elas + T visc, T,f = T,f visc, Ti visc = µ i v i + v i T = µ i S v i, i =,, f where S v i indicates the symmetric part of the velocity gradient and µ,,f represent shear viscosities.,
23 4 The boundary conditions are given as follows. Let n be the unit normal vector which points outward from Ω t to R t. The consevation of mass of fluid gives rise to v f v n = φ v v n w, 3.9 and moreover no-slip along the surface Γ t is imposed: v f v = v v r. 3.0 The force balance across the interface Γ t is T f n T n T n + [p]n = 0, [p] = p p f. 3. Finally, on U, we set v f = 0. 3.
24 Chapter 4 A Theory of Thermal Effects in Gels In this chapter, we concern two laws of thermodynamics, balance of energy and entropy inequality. In particular, a control volume approach is used. Let ω 0 be a small fixed volume in time inside Ω t. We first set up the equations of balance of energy and additional boundary conditions, and from these we derive the evolution equations for temperatures θ and θ f as well as constitutive restrictions. Then Galilean relativity is discussed. The independent variables are given by φ, φ, θ, θ f, p, p f, T elas, v i, i =,, f. Densities of the internal energy, the free energy, and the entropy of the gel and surrounding fluid, and heat flux vectors into Ω t and R t across Γ t are represented by ε, ε f, ψ, ψ f, η, η f, q, q f, repectively, in which first six are density functions per unit volume and the last two the ones per unit area. We postulate constitutive equations of the form T elas = T elas F, φ, φ, θ, ε = εf, φ, φ, θ, ψ = ψf, φ, φ, θ, η = ηf, φ, φ, θ, q = qf, φ, φ, θ, θ, 5
25 and 6 ε f = ε f θ f, ψ f = ψ f θ f, η f = η f θ f, q f = q f θ f, θ f. We assume that the system U is adiabatic. Forms of constitutive equations are constraints by the Second Law of Thermodynamics, which requires that entropy in an isolated system is not decreasing in time. 4. Energy Balance Let ε i, i =,, f, be the internal energy density functions for each component i. Then the internal energy density ε for the whole gel may be written by ε = φ ε + φ ε. 4. Under the assumption that external heat is not supplied, the first law of thermodynamics in the control volume ω 0 is stated that the amount of change of internal energy in a fixed volume ω 0 in time is equal to the sum of energy flux across the boundary of ω 0, power expended on the boundary, and the heat flux through the boundary, and it has the form d ε dx = φ ε v + φ ε v n da dt ω 0 ω 0 [ + T φ piv n + T φ piv n ] da q n da. ω 0 ω 0 Since ω 0 is arbitrary, this leads to the local form of balance of energy in Ω t : ε t + φ ε v + φ ε v = T v + T v pφ v + φ v q. 4. Similarly, the local version of energy balance in R t is ε f t + ε f v f = T f v f p f v f q f. 4.3
26 The energy conservation in the system allows an equation on the interface Γ t : 7 ε f ε [p]w + n T n φ w n T f n w T n r + q f q n = This boundary condition is obtained as follows. Consider a tiny volume ω which Γ t cuts through and moves along the interface. Let U be the total internal energy in the volume ω. The rate of change of U is expressed as du dt = d ε dx + d ε f dx. 4.5 dt ω Ω t dt ω R t By applying the Reynolds transport theorem, the first term becomes d ε ε dx = dt ω Ω t ω Ω t t + v ε dx ε = ω Ω t t dx + εv n da. ω Ω t 4.6 By use of equations 4.-4., 4.6 can be written as φ ε v + φ ε v + T v + T v pφ v + φ v q dx ω Ω t + φ ε + φ ε v n da, ω Ω t which yields ε φ v v n + T n v + T n v pφ v n pφ v n q n da. ω Ω t Likewise, the second term of 4.5 is computed as d ε f ε f dx = dt ω R t ω R t t + v ε f dx ε f = ω R t t dx ε f v n da ω R t = ε f v f v n T f n v f + p f v f n + q f n da. ω R t Since ω Ω t = ω Ω t ω Γ t and ω R t = ω R t ω Γ t, the first
27 law of thermodynamics implies 0 = ε φ v v n + T n v + T n v pφ v n pφ v n q n ω Γ t = ω Γ t + ε f v f v n T f n v f + p f v f n + q f n da ε f ε p w + T n v + T n v T f n v f + p f v f pv n + q f q n da, where 3.9 is used, and by the arbitrariness of ω, we obtain ε f ε pw + T n v + T n v T f n v f + p f v f pv + q f q n = 0. Let us simplify this equation. By use of 3., we see T n v + T n v T f n v f = T f n v f v + T n v v + [p]n v = T f n nw + r + T n n w + r + [p]n v φ 8 and also = n T f n w T f n r + n T n φ w + T n r + [p]n v, therefore we have [p f v f pv ] n = [p f v f v + p f v pv ] n = p f w [p]v n, ε f ε pw + n T n φ w n T f n w T f n r + T n r + p f w + q f q n = 0. Using 3. again results in 4.4: ε f ε [p]w + n T n φ w n T f n w T n r + q f q n = Boundary Conditions Deduced by Thermodynamics Second Law Some boundary conditions are deduced by considering variations in temperature. We take the entropy into consideration. Let us denote the entropy density functions by
28 η i, i =,, f, for each constituent i and by η = φ η + φ η for the entire gel, Ω t. For any fixed ω 0, as in the preceding section, the change of entropy in time can be written as d q η dx = φ η v + φ η v n da dt ω 0 ω 0 ω 0 θ n da + Q dx, ω 0 where the first and second integrals on the right hand side represent the rates of entropy flow across ω 0 conveyed by the material coming in due to convection and caused by heat, respectively, and the third integral represents the net entropy production. Since ω 0 was arbitrary, we have the local version of the second law of thermodynamics: η q t + φ η v + φ η v + = Q 0, in Ω t. 4.7 θ In the same way, the local form in R t is given by η f t + η f v f + qf θ f 9 = Q f Here, Q and Q f mean the respective net entropy production densities in Ω t and R t. Now we claim that d dt U η dx = Q dx + Ω t Q f dx + R t Q Γ da, Γ t 4.9 where qf Q Γ = η f η w + q n, 4.0 θ f θ and the second law of thermodynamics implies that Q Γ The following is the proof of 4.9. First we use the Reynolds transport theorem,
29 4.7 and 3.9 to get d η dx dt Ω t η = Ω t t dx + ηv n da Γ t q = φ η v + φ η v + Q dx + φ η + φ η v n da Ω t θ Γ t = Q dx + φ η v + φ η v n q Ω t Γ t θ n + φ η + φ η v n da = Q dx + η φ v v n q Ω t Γ t θ n da = Q dx + η w q Ω t Γ t θ n da. A similar computation provides d η f dx = Q f dx + η f w + q f n da, dt R t R t Γ t θ f and sum of these two results brings the equation 4.9, since d η dx = d η dx + d η f dx. dt U dt Ω t dt R t 0 Using 4.4 and thermodynamic relations ε i = ψ i + θ i η i, i =,, f, in conjunction with the assumption that θ = θ f on Γ t, 4.0 takes the form qf Q Γ = η f η w + q n θ f θ = η f η w + θ q f q n = T n T f n n n + [p] + ψ ψ f w + θ φ θ T n r. From 4. we can impose the boundary conditions on Γ t : T µ w = n T f n n n + [p] + ψ ψ f, 4. φ µ r = T n, 4.3 where µ and µ are non-negative constants.
30 4.3 Evolution Equations for Temperature In this section we derive equations for temperature development in time of the system. In order to obtain them, first of all, we must specify T i, f i, ε, ε i and q, q f. We shall take as our starting point an expression for the Helmholtz free energy density ψ in terms of φ, F and θ. We recall that F satisfies the equations t F + v F = v F, 4.4 φ detf = φ ref. 4.5 Given ψφ, F, θ, from the usual thermodynamic relations, we have ε = ψ + θη, η = ψ θ. 4.6 We need expressions for ε i. The crucial insight here is that it is possible to determine the expressions for ε i given only the expression for ε. Note that 4. is hardly sufficient to obtain an expression for ε i given ε. Consider a spatially uniform gel of volume V, so that φ, F and θ take the same values everywhere. The total energy E of this gel is given by Define E = V εφ, F, θ. 4.7 V i = V φ i, G = detf /3 F = φ φ ref /3 F. 4.8 View E as a function of V i, G and θ instead of V, φ i, F and θ, and call this function Ê V /3 ÊV, V, G, θ = V + V ε, φ /3 V V + V ref G, θ. 4.9 V + V The energy per unit volume of polymer and fluid should be given as the partial derivative of Ê with respect to V and V respectively ε i = Ê V i. 4.0
31 We have used the symbol ε i to emphasize the fact that ε i is a function of V, V, G and θ instead of φ, F and θ. The function Ê is a one-homogeneous function in V and V ÊsV, sv, G, θ = sêv, V, G, θ for s > By taking the derivative of the above in s and letting s =, we have Euler s formula for homogeneous functions Dividing by V, we have Ê = V Ê V + V Ê V = V ε + V ε. 4. ε = Ê V = V V ε + V V ε = φ ε + φ ε. 4.3 This is nothing other than 4. except that the functions are in terms of V, V, G and θ. For the above expressions to be useful, we must check that they are independent of the arbitrary volume parameter V. This can be seen as follows. Dividing 4.9 by V, we see that ε may be viewed as a function of φ and φ instead of V and V φ /3 εφ, φ, G, θ = φ + φ ε, φ /3 φ φ + φ ref G, θ. 4.4 φ + φ It is easily checked that 4.0 is equivalent to the following ε i φ, φ, G, θ = ε φ i. 4.5 We shall henceforth view ε, ε i as functions of φ, φ, G and θ instead of V, V, G and θ. By restricting ε i to φ + φ =, we obtain the expressions for ε i in terms of φ, F and θ /3 F, θ ε i φ, F, θ = ε i φ, φ, φ φ ref. 4.6 Similarly to ε, we define ψ, ψ i, η and η i as follows. φ /3 ψφ, φ, G, θ = φ + φ ψ, φ /3 φ φ + φ ref G, θ, 4.7 φ + φ φ /3 ηφ, φ, G, θ, = φ + φ η, φ /3 φ φ + φ ref G, θ, 4.8 φ + φ ψ i = ψ φ i, η i = η φ i. 4.9
32 Note that, similarly to ψ, ε and η, the above quantities are also related to one another through Given expressions for ε i in terms of ε, we now simplify 4.. The left hand side of 4. yields We have ε t + φ ε v + φ ε v = ε φ φ t + ε φ φ t + ε G : G t + ε θ ε φ φ t + ε φ φ t + φ ε v + φ ε v = ε φ v ε φ v + φ φ ε ε =φ v + φ v, φ φ θ t + φ ε v + φ ε v. φ ε φ v + φ ε φ v where we used 3.3 and 4.5 in the first equality. The third term in the last line of 4.30 can be written as ε G : G t = φ/3 ref φ + φ 4/3 φ /3 ε F : = φ ε 3 φ t F : F + ε F : F t = 3 φ φ v ε F t φ /3 ref φ /3 F : F + ε F : v F + v F ε = v 3 φ φ F : F + ε F : F ε + v : F F T ε 3 F : F I, ε F : F ε F ij, x k k = i,j F ij 4.3 where I is the 3 3 identity matrix, k denotes the k-th component of a vector, x k are the Cartesian spatial coordinates and F ij denotes the ij element of the deformation gradient. We used 4.4, 4.8 in the first equality, 3. in the second equality and 3.3 and 4.4 in the third equality.
33 The right hand side of 4. yields T v + T v pφ v + φ v q =T : v + T : v f v f v q, where we used 3.5. Combining 4.30, 4.3, 4.3 and 4.33, we find that 4. can be rewritten as follows. θ ε C v t = v φ φ + v : T ε ε v φ F F T + 3 φ + 3 φ φ ε F : F ε F : F I + f + T : v q, where C v = ε/ θ is the heat capacity at constant volume in Ω t. + ε F : F f Because ε f depends only on temperature, 4.3 can be expressed as 4.34 ε f θ f θ f t + v f ε f = T f : v f + v f T f v f p f q f, 4.35 where 3.8 has been used. By 3.7, 4.35 is reduced to C fv θ f t = T f : v f v f ε f q f, 4.36 as above C fv = ε f / θ f denotes the heat capacity at constant volume in R t. 4.4 Entropy Inequality and the Resulting Relations The purpose of this section is to derive thermodynamic restrictions of the system from the entropy inequality. The same point of view is taken for the energy and entropy densities as in Section 4.3. We have seen in Section 4. that the local form of the net entropy production inside the gel was given by 4.7: Q = η q t + φ η v + φ η v +, θ
34 which can be written by the chain rule as Q = η φ φ t + η φ φ t + η G : G t + η θ q θ t + φ η v + φ η v +. θ 4.37 Some calculus similar to and 5 η θ θ t = ε θ θ θ t = θ C θ v t 4.38 lead 4.37 to η Q = φ v φ + v : η F F T 3 η + φ v φ η F : F I v 3 φ φ + θ C θ q v t + θ η F : F. + η F : F Using 4.34, then by collecting like terms, we have Q = v g f + v g f + v : T T elas + v : T q θ, θ θ where ε η g = φ θ + ψ φ φ φ 3φ F : F + ψ F : F, ε η g = φ θ, φ φ T elas = ψ F F T ψ 3 F : F I We now claim that g + g = By 4.3 and 4.6, we have ε = φ ε φ + φ ε φ. 4.4 Take the gradient on both sides of the above ε ε ε ε : G + G θ θ = φ + φ. 4.4 φ φ
35 Proceeding as in the first line of 4.3, we have ε 3 φ φ F : F + ε F In exactly the same fashion, we find that η 3 φ φ F : F + η F ε ε : F + θ θ = φ φ η η : F + θ θ = φ φ 6 ε + φ φ η + φ φ Multiply 4.44 by θ and subtract this from Using the definition of g and g in 4.39, we have The right hand side vanishes by ε g + g = θ θ η θ θ by Thus, the net entropy production of a small volume in the gel can be simply written Q = θ v v g f + v : T T elas + v : T q θ. θ θ In a completely similar way, we get Q f = η f t + η f v f + qf = η f θ f θ f t + η f v f + = θ f ε θ f θ f θ f qf t + η f v f + = θ f C fv θ f t + η f v f + θ f qf θ f qf = θ f T f : v f v f ε f q f + η f v f + θ f qf θ f, 4.46 where 4.36 is used in the last equality. Note that v f ε f = εf ηf v f θ f = v f θ f = v f η f, θ f θ f θ f θ f
36 because of 4.38 and the fact that ε f is a function of θ f only. By use of this identity and 3.8, 4.46 can be written as 7 Q f = θ f T f : v f θ f q f θ f. Therefore, we have the following relations so that the non-negativity of the net entropy production density, referred to as the local version of the second law of thermodynamics, is satisfied for any choices {φ, v, v, v f, θ, θ f }: f = g κv v, f = g κv v = f, q = λ θ, q f = λ f θ f, 4.47 in which κ is the friction coefficient, and λ, λ f are the heat conductance. Those coefficients are positive and κ, λ may depend on the volume fraction φ of polymer networks. On applying all these restrictions 4.47 to 4.34 and 4.36, we are led to the evolution equations of a thermoelastic gel and surrounding fluid of the form: θ C v t = κ v v + µ S v + µ S v + λ θ η η θ φ v + φ v φ φ η + θv φ 3φ F : F + η F : F η + θ v : 3 F : F I η F F T 4.48 and C fv θ f t = µ f S v f v f ε f + λ f θ f We end this section with the following proposition that summarizes the results.
37 8 Proposition Suppose that the balance laws of mass 3.3, linear momentum , and energy and the constraint 3. hold. Given the thermodynamic relations 4.6, suppose that the entropy inequality holds for all processes {φ, v, v, v f, θ, θ f }. Then the constitutive restrictions 4.39 and 4.47 hold. 4.5 Galilean Invariance In this section, we prove the Galilean invariance of equations in the motion of the gel. Let x, t be a pair of position and time in another inertial frame, say S, whose motion is supposed to be uniform with a constant velocity v c relative to the present frame of reference. Let the velocities for polymer networks and fluid in the gel in S be denoted by v i = v i x, t, i =,, respectively. By assuming that t = t, we have then x = x v c t, v x, t = v x, t v c, 4.50 v x, t = v x, t v c. Easy calculus shows that for any scalar function f, f = f, f t = v c f + f t, and accordingly, v i = v i, v i = v i, from which we can easily see that equations 3.,3.4 and 3.5 hold in S. Here has been used to represent the gradient notation in S. As for the evolution equation for temperature 4.48, since C v θ t = C v θ t + C vv c θ
38 and η κ v v + µ S v + µ S v + λ θ θ φ v φ η + θ v φ : F + η η : F + θ v : : F 3φ F F 3 F = κ v v + µ S v + µ S v + λ θ η η θ φ v v c + φ v v c φ φ η + θv v c φ 3φ F : F + η F : F + θ v : 9 η + φ v φ I η T F F η 3 F : F I η F F T, in order for 4.48 to be invariant under Galilean transformation, it should be satisfied that η η C v θ = θ φ + φ η φ φ φ 3φ F : F η F : F. From C v = ε θ = θ η θ, we need only show that η η η θ θ = φ + φ η φ φ φ 3φ F : F η F It is clearly true due to : F. 4.5 Same argument applies to equations in R t, therefore, we conclude that the governing equations in both the gel and its environment fluid are the same in any intertial frame of reference.
39 Chapter 5 Gels in One-dimensional Geometry We restrict all the equations that we have derived in the preceding chapter to one spatial dimension. When elaborating a numerical scheme, the use of spatial description can involve difficulties because of moving interfaces of the gel. Thus, we carry out the spatial to referential transformation, and build a numerical algorithm. 5. Governing Equations in One-Dimensional Space We consider the thermoelastic motion of a one-dimensional gel with Dirichlet boundary conditions. Let U = [0, ] R be the region containing the gel and fluid, and we express the domain of polymer networks at time t as Ω t = L l t, L r t, 0 < L l t < L r t <, and Ω 0 = L l 0, L r 0 for the reference domain. The boundaries L i t move with velocities v L i, t, i = l, r. We assume that the velocity of surrounding fluid is zero at the ends x = 0,, that is, v f = 0 at x = 0,. Furthermore, we specify temperature values at both ends: θ f = θ 0 at x = 0, θ f = θ at x =. 30
40 The deformation gradient is a positive number: F = ϕ x := α, α > 0, and by the relation φ det F = φ ref, it is proportional to the inverse of φ : 3 α = φ ref φ, φ ref = φ x, We let the total free energy densities at each point in Ω t and R t, respectively, be of the form ψ = k B aθφ ln φ + bθφ ln φ + χθφ φ cθ ln θ + φ W elas F, 5. ψ f = c f θ f ln θ f, 5.3 where k B is the Boltzmann constant, aθ, bθ and the Flory-Huggins interaction parameter, χθ, are linear functions of θ, W elas F represents the elastic stored energy function, and c, c f are positive constants. Because of 5., we can view the elastic energy, W elas, as a function of φ, and thus we can write 5. in the form ψ = k B aθφ ln φ + bθφ ln φ + χθφ φ cθ ln θ + φ W elas φ. 5.4 Also, by 4.7, we have ψφ, φ, θ =k B [aθφ ln φ φ + φ φ + φ c θ ln θ + φ W elas φ + bθφ ln φ φ + φ φ + φ. + χθ φ ] φ φ + φ From the thermodynamic relation η = ψ/ θ, entropy density function is given by η = k B [a φ θφ ln + b φ θφ ln + χ θ φ ] φ φ + φ φ + φ φ + φ + φ + φ c ln θ +, where we have used the prime notation to denote the derivatives over θ. In a similar fashion, we get η f = c f ln θ f +.
41 The constant c has the same value as heat capacity C v : 3 C v = c, since C v = ε θ = θ ψ + θ η = ψ ψ θ = θ ψ θ θ θ = θ η θ = c. By the same reasoning, we obtain also C fv = c f. Further, the associated values are given by T elas = 0, T visc i = µ i v ix, i =,, f, f =f = g κv v, ε η g = g = φ φ θ x φ x = φ ψ φ + θ η x. x φ For simple notation, we have used subscripts to indicate a partial differentiation. Therefore, from and 4.48 it follows that in Ω t = L l t, L r t, φ + φ =, 5.5 φ t + φ v x = 0, 5.6 φ v + φ v x = 0, 5.7 µ v xx φ p x + f = 0, 5.8 µ v xx φ p x + f = 0, 5.9 θ C v t = κv v + µ v x + µ v x + λθ xx 5.0 η η θ φ v φ + φ v x φ x = κv v + µ v x + µ v x + λθ xx +θk B φ x v v a θφ + b θφ χ θφ φ c θ x φ v + φ v,
42 and from 3.8, 3.7, 4.49, we obtain that in R t = 0, L l t L r t,, 33 v f x = 0, 5. µ f v f xx p f x = 0, 5. C fv θ f t = µ f v f x c f θ fx v f + λ f θ f xx, 5.3 where we have assumed λ and λ f are constants. The equations for boundaries ,4.4,4. yield at x = L l t, at x = L r t, v v f = φ v v = w, 5.4 v f v = v v r = 0, 5.5 µ f v f x µ v x µ v x + [p] = 0, 5.6 ε f ε [p]w + µ v x φ w µ f v f x w + λ f θ f x λθ x = 0, 5.7 µ w = µ f v f x φ µ v x + [p] + ψ ψ f, 5.8 v f v = φ v v = w, 5.9 v f v = v v r = 0, 5.0 µ f v f x µ v x µ v x + [p] = 0, 5. ε f ε [p] w + µ v x φ w µ f v f x w + λθ x λ f θ f x = 0, 5. µ w = µ f v f x φ µ v x + [p] + ψ ψ f. 5.3 We note that the interfaces, x = L i t, move with velocities v L i, t, i = l, r, respectively. In addition, and θ = θ f, at x = L l t, L r t, 5.4 v f = 0, θ f = θ 0 at x = 0, 5.5 v f = 0, θ f = θ at x =. 5.6
43 34 From 5. and , we find v f 0 in R t, so it follows from 5., that p f is constant. Now, we set p f = 0 in L r t, for simplicity. Moreover, along with 5.5 and 5.7, both equations 5.4 and 5.9 imply that in Ω t, Adding 5.8 and 5.9, we have Integrating on x gives φ v + φ v = 0 or v = φ v φ. 5.7 p x = µ v xx + µ v xx. p = µ v x + µ v x + kt. By comparing with 5.6 and 5., kt p f = constant in 0, L l t and kt 0 in L r t,. Then, 5.8 can be rewritten by φ µ v xx φ µ v xx = κv v g. Therefore, by use of 5.7, we have the following governing equations for {φ, v, θ, θ f } : in Ω t = L l t, L r t, φ t + φ v x = 0, 5.8 φ v φ µ v xx + φ µ = κ v g, 5.9 φ xx φ [ ] θ C v t = κ φ v + µ v x φ v + µ + λθ xx 5.30 φ x +k B θv φ x a θ + φ b θ φ χ θ, φ where [ g =φ x k B aθ φ bθφ + χθφ φ φ φ W elas φ + φ W elas φ ].
44 Equations for the interfaces and surrounding fluid are given by at x = L l t, φ φ v µ v = µ v x + µ + ψ ψ f, 5.3 φ φ x µ ε ε f + p p f v + φ v φ v + λθ x λ f θ f φ x = 0, 5.3 x 35 θ = θ f, 5.33 at L r t, φ φ v µ v = µ v x + µ + ψ ψ f, 5.34 φ φ x µ ε ε f + pv + φ v φ v + λθ x λ f θ f φ x = 0, 5.35 x θ = θ f, 5.36 where } ψ ψ f =k B {bθ aθ φ + bθ ln φ + φ χθ + c f cθ ln θ φ W elas φ, } ε ε f =k B {bθ aθ φ + bθ ln φ + φ χθ c f cθ φ W elas φ } θk B {b θ a θ φ + b θ ln φ + φ χ θ, in R t = 0, L l t L r t,, and at each end of the system, C fv θ f t = λ f θ f xx, 5.37 θ f = θ 0 at x = 0, θ f = θ at x = Referential Form for the Equations of Gels In developing a numerical method for the spatially one-dimensional model, moving boundaries of the gel can be a major obstacle to using the Eulerian description. To
45 36 cope with the problem, we adopt the Lagrangian approach in which each material particles are identified and tracked. We discretize the gel using a Lagrangian mesh so that the Lagrangian nodes move at the polymer network velocity v. The grid is extended to the surrounding fluid domain as will be described below. For that reason, we now rewrite the governing equations in the Lagrangian frame. Let the material description for the functions be expressed as x = ϕx, t, X Ω 0 = L l 0, L r 0, F = ϕ X = α > 0, V X, t = v ϕx, t, t, α = αx, t, ΦX, t = φ ϕx, t, t, Φ 0 X = φ ref ϕx, 0, 0, ΘX, t = θϕx, t, t, HΦ, t = W elas ΦX, t, and ΘX, t = θ f X, t, for X R 0. We make a Eulerian to the Lagrangian transformation, in particular, v x = v x = v X X x = v X F = V α X, v xx = v = V X x x X α X x = V α X α X θ t x, t = θx, t θ x x, t v x, t [ ] = t θ ϕx, t, t X X=ϕ x,t, X θϕx, t, t α v ϕx, t, t = Θ t X, t Θ X, tv X, t. α X Here the notation t represents the partial derivative with respect to t holding X X fixed.
46 Letting linear functions aθ, bθ, and χθ be the following form 37 aθ = a θ + a 0, bθ = b θ + b 0, χθ = z θ + z 0, we arrive at the following equations in the referential coordinates in Ω t : where ΦV + Φ V = 0, 5.39 µ V Φ P α X α X α X + F = 0, 5.40 µ V Φ P α X α X α X + F = 0, 5.4 Θ C v t = κv V + µ V α + µ V X α + λ Θ X α X α X Θ [ ] Φ +k B α V V a Φ + b Φ z Φ Φ 5.4 X + C v Θ α X ΦV V, Φ = Φ, F X, t = F X, t, and [ F X, t = κv V Φ k B a Φ + b Φ z Φ Φ Θ α X + k B a 0 Φ + b 0 Φ z 0 Φ Φ + Φ Φ H Φ + Φ ] H. Φ Φ Next, for the temperature Θ in R t, we use a moving mesh frame. Let x l = ϕ t L l be the position of the left boundary, L l, at time t. Then, there follows and θ f t = θ f t θ f x = x x = ϕl l, t L l X, for X 0, L l, + x θ f X t x = Θ X t + Θ V L l L l X, θf = L l Θ = x x ϕl l, t X Ll ϕl l, t Θ X.
47 By transformation, the evolution equation of temperature for the outer fluid on the left side takes the form C Θ fv t + V L l X Θ Ll Θ = λ f L l X ϕl l, t X. 38 In the exactly similar manner, we let x r = ϕ t L r be the position of the right boundary, L r, at time t. For the fluid on the right side x r,, we have ϕlr, t x = X, for X 0, L r, L r and hence C Θ X fv t + V Θ Lr Θ L r = λ f L r X ϕl r, t X. The boundary conditions and can be expressed in the referential setting, respectively, as and µ α µ V R + P P f V V X + µ α µ α Φ V X P + P f = 0, V X + E + P P f = 0, µ V α Φ X V + λ α Θ X λ f L l Θ ϕ t L l X = 0, R + P V µ α µ V V X + µ α µ α Φ µ V α Φ X V + λ α V X P = 0, V X + E + P = 0, Θ X λ f L r Θ ϕ t L r X = 0, where E = ψ ψ f, R = ε ε f, and those are computed by EX, t =k B Θ b a Φ + b ln Φ + z Φ + c f cθ ln Θ + k B b 0 a 0 Φ + b 0 ln Φ + z 0 Φ Φ H Φ, RX, t =k B b 0 a 0 Φ + b 0 ln Φ + z 0 Φ c f cθ Φ H Φ.
48 In addition, Θ0, t = θ 0, Θ, t = θ Numerical Algorithmn In this section, we construct an algorithm to compute the solution of the governing equations in Lagrangian coordinates. It is a time-marching algorithm which includes Newton s method and the backward Euler method at each time step Discretization We discretize the space [0, ] to have M + grid points: X 0 = 0, X, X,, X }{{} m,, X m+n,, X }{{} M =, =L l =L r for some integers m, n, M such that 0 < m < m + n < M. The function values by the deformation mapping and temperatures at the grid points are denoted by ϕx i, t = ϕ i, ΘX i, t = Θ i, i = m, m +,, m + n, ΘX i, t = Θ i, i {,,, m } {m +,, M }. We assume that the pressure in the gel is defined at midpoints, P j/, j = m +, m + 3,..., m + n, and its spatial derivative at the kth point is calculated as P P k+/ P k /, k = m +,, m + n. X X k For each point X i or X i+/ i = m,, m + n, we let the values of the deformation gradient α, volume fraction Φ, velocities V, W of the polymer and fluid,
49 and temperature Θ be of the form, respectively: 40 α i+ = ϕ i+ ϕ i X, Φ i+ V i = ϕ t X i, t := V i, V i := W i, V i+ Θ i+ = V i+ + V i, W i+ = Θ i+ + Θ i. = Φ 0, α i+ i+ = W i+ + W i, Then Φ i, the volume fraction of polymer networks at a node point, is computed by taking the average of the values defined at midpoints, and α at the node is determined by the relation 5.: Φ i+ = Φ i+ + Φ i++, α i+ = Φ0 i+, i = m,, m + n. Φ i+ From direct calculation Φ X X i, t Φ i+ Φ i, X [ V X i, t V X i + X V X α X X α X α X [ = V X X α i+ X i+ α i [ ] = V i+ V i V i V i, X ϕ i+ ϕ i ϕ i ϕ i [ ] Θ X i, t Θ i+ Θ i Θ i Θ i, X α X X ϕ i+ ϕ i ϕ i ϕ i X i X V X we obtain the discretized equations as follows: for i = m,, m + n, X i ] ] Φ i V i + Φ i W i = 0, 5.43
50 for i = m +,, m + n, [ ] µ V i+ V i V i V i Φ i α i X ϕ i+ ϕ i ϕ i ϕ i α i P i+ P i X 4 + F X i = 0, 5.44 [ ] µ W i+ W i W i W i Φ P i i+ α i X ϕ i+ ϕ i ϕ i ϕ i α i P i X + F X i = 0, 5.45 C v Θ t X i = κv i W i + µ + λ α i X Θ i+ α i + C v α i Vi+ V i + µ X αi ] [ Θ i+ Θ i ϕ i+ ϕ i Θ i Θ i ϕ i ϕ i Θ i X + k B Θ i α i V i W i Φ i V i W i [ a Φ i + b Φ i z Φ i Φ i Wi+ W i X ]Φ i+ Φ i X, 5.46 where F X i, t = κv i W i [ Φ i+ Φ i k B a Φ i + b Φ i z Φ i Φ i Θ i α i X + k B a 0 Φ i + b 0 Φ i z 0 Φ i Φ i + Φ i Φ i Hi+ H i H Φ i+ Φ i+ H i i + Φ i H ] i H i. Φ i+ Φ i Φ i Φ i Since moving interfaces influence the lengths of the left and right regions that the outside fluid occupy at each time, we repartition those areas. Let x m = ϕ t X m be the position of the left boundary, X m, at time t. We divide 0, x m into m equal intervals with the length X l. Similarly, the fluid region on the right, x m+n,, is rezoned to have M m n intervals of the uniform length X r. Then application of the above rearrangement gives rise to the equations for the surrounding
Nonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More information4 Constitutive Theory
ME338A CONTINUUM MECHANICS lecture notes 13 Tuesday, May 13, 2008 4.1 Closure Problem In the preceding chapter, we derived the fundamental balance equations: Balance of Spatial Material Mass ρ t + ρ t
More informationAnalysis and numerics of the mechanics of gels
Analysis and numerics of the mechanics of gels A DISSERTATION SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Brandon Michael Chabaud IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More informationLevel Set Tumor Growth Model
Level Set Tumor Growth Model Andrew Nordquist and Rakesh Ranjan, PhD University of Texas, San Antonio July 29, 2013 Andrew Nordquist and Rakesh Ranjan, PhD (University Level Set of Texas, TumorSan Growth
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationCONSTITUTIVE THEORIES: BASIC PRINCIPLES
1 6.161.6 CONSTITUTIVE THEORIES: BASIC PRINCIPLES I-Shih Liu Instituto de Matemática, Universidade Federal do Rio de Janeiro, Brasil Keywords: Frame of reference, change of observer, Euclidean objectivity,
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationA dynamic model of polyelectrolyte gels. A dissertation submitted to the faculty of the graduate school of the university of minnesota by.
A dynamic model of polyelectrolyte gels A dissertation submitted to the faculty of the graduate school of the university of minnesota by Haoran Chen In partial fulfillment of the requirements for the degree
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationarxiv: v2 [math.ap] 16 Oct 2012
EFFECTS OF PERMEABILITY AND VISCOSITY IN LINEAR POLYMERIC GELS B. CHABAUD AND M. C. CALDERER arxiv:121.3813v2 [math.ap] 16 Oct 212 Abstract. We propose and analyze a mathematical model of the mechanics
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationOn pore fluid pressure and effective solid stress in the mixture theory of porous media
On pore fluid pressure and effective solid stress in the mixture theory of porous media I-Shih Liu Abstract In this paper we briefly review a typical example of a mixture of elastic materials, in particular,
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationA CONTINUUM MECHANICS PRIMER
A CONTINUUM MECHANICS PRIMER On Constitutive Theories of Materials I-SHIH LIU Rio de Janeiro Preface In this note, we concern only fundamental concepts of continuum mechanics for the formulation of basic
More informationIn this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case
Section.. Thermoelasticity In this section, thermoelasticity is considered. By definition, the constitutive relations for F, θ, Gradθ. This general case such a material depend only on the set of field
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationChapter 2: Basic Governing Equations
-1 Reynolds Transport Theorem (RTT) - Continuity Equation -3 The Linear Momentum Equation -4 The First Law of Thermodynamics -5 General Equation in Conservative Form -6 General Equation in Non-Conservative
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationCONSTITUTIVE THEORIES FOR THERMOELASTIC SOLIDS IN LAGRANGIAN DESCRIPTION USING GIBBS POTENTIAL YUSSHY MENDOZA
CONSTITUTIVE THEORIES FOR THERMOELASTIC SOLIDS IN LAGRANGIAN DESCRIPTION USING GIBBS POTENTIAL BY YUSSHY MENDOZA B.S. in Mechanical Engineering, Saint Louis University, 2008 Submitted to the graduate degree
More informationComputational Modelling of Mechanics and Transport in Growing Tissue
Computational Modelling of Mechanics and Transport in Growing Tissue H. Narayanan, K. Garikipati, E. M. Arruda & K. Grosh University of Michigan Eighth U.S. National Congress on Computational Mechanics
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationLONG-TIME EXISTENCE OF CLASSICAL SOLUTIONS TO A 1-D SWELLING GEL
LONG-IME EXISENCE OF CLASSICAL SOLUIONS O A -D SWELLING GEL M. CARME CALDERER AND ROBIN MING CHEN Abstract. In this paper we derived a model which describes the swelling dynamics of a gel and study the
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationPhysical Conservation and Balance Laws & Thermodynamics
Chapter 4 Physical Conservation and Balance Laws & Thermodynamics The laws of classical physics are, for the most part, expressions of conservation or balances of certain quantities: e.g. mass, momentum,
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationRole of thermodynamics in modeling the behavior of complex solids
IWNET Summer School 2015 Role of thermodynamics in modeling the behavior of complex solids Bob Svendsen Material Mechanics RWTH Aachen University Microstructure Physics and Alloy Design Max-Planck-Institut
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationInterpreting Differential Equations of Transport Phenomena
Interpreting Differential Equations of Transport Phenomena There are a number of techniques generally useful in interpreting and simplifying the mathematical description of physical problems. Here we introduce
More informationChapter 3. Forces, Momentum & Stress. 3.1 Newtonian mechanics: a very brief résumé
Chapter 3 Forces, Momentum & Stress 3.1 Newtonian mechanics: a very brief résumé In classical Newtonian particle mechanics, particles (lumps of matter) only experience acceleration when acted on by external
More informationIntroduction to Continuum Mechanics
Introduction to Continuum Mechanics I-Shih Liu Instituto de Matemática Universidade Federal do Rio de Janeiro 2018 Contents 1 Notations and tensor algebra 1 1.1 Vector space, inner product........................
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationIII.3 Momentum balance: Euler and Navier Stokes equations
32 Fundamental equations of non-relativistic fluid dynamics.3 Momentum balance: Euler and Navier tokes equations For a closed system with total linear momentum P ~ with respect to a given reference frame
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationTHE STOKES SYSTEM R.E. SHOWALTER
THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes
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 informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationThermodynamic form of the equation of motion for perfect fluids of grade n
Thermodynamic form of the equation of motion for perfect fluids of grade n Henri Gouin To cite this version: Henri Gouin. Thermodynamic form of the equation of motion for perfect fluids of grade n. Comptes
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 informationFluid Mechanics Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 15 Conservation Equations in Fluid Flow Part III Good afternoon. I welcome you all
More informationFrom the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)
Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationBulletin of the Transilvania University of Braşov Vol 10(59), No Series III: Mathematics, Informatics, Physics, 83-90
Bulletin of the Transilvania University of Braşov Vol 10(59), No. 1-2017 Series III: Mathematics, Informatics, Physics, 83-90 GENERALIZED MICROPOLAR THERMOELASTICITY WITH FRACTIONAL ORDER STRAIN Adina
More informationAA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations
AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationApplication of Chimera Grids in Rotational Flow
CES Seminar Write-up Application of Chimera Grids in Rotational Flow Marc Schwalbach 292414 marc.schwalbach@rwth-aachen.de Supervisors: Dr. Anil Nemili & Emre Özkaya, M.Sc. MATHCCES RWTH Aachen University
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
More informationNotes on Entropy Production in Multicomponent Fluids
Notes on Entropy Production in Multicomponent Fluids Robert F. Sekerka Updated January 2, 2001 from July 1993 Version Introduction We calculate the entropy production in a multicomponent fluid, allowing
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationFundamental equations of relativistic fluid dynamics
CHAPTER VI Fundamental equations of relativistic fluid dynamics When the energy density becomes large as may happen for instance in compact astrophysical objects, in the early Universe, or in high-energy
More informationNonlinear stability of steady flow of Giesekus viscoelastic fluid
Nonlinear stability of steady flow of Giesekus viscoelastic fluid Mark Dostalík, V. Průša, K. Tůma August 9, 2018 Faculty of Mathematics and Physics, Charles University Table of contents 1. Introduction
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationCoordinate systems and vectors in three spatial dimensions
PHYS2796 Introduction to Modern Physics (Spring 2015) Notes on Mathematics Prerequisites Jim Napolitano, Department of Physics, Temple University January 7, 2015 This is a brief summary of material on
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationSupplement to Molecular Gas TitleBirkhäuser, Boston, 007 Dynamic Version Authors Sone, Yoshio Citation Yoshio Sone. 008 Issue Date 008-09-0 URL http://hdl.handle.net/433/66098 Right Type Book Textversion
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationLecture 3: 1. Lecture 3.
Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector
More information6.2 Formulation of thermomechanical constitutive equations
6 CONSTITUTIVE EQUATIONS 61 The need for constitutive equations Basic principles of continuum mechanics, namely, conservation of mass, balance of momenta, and conservation of energy, discussed in Chapter
More informationMath Tune-Up Louisiana State University August, Lectures on Partial Differential Equations and Hilbert Space
Math Tune-Up Louisiana State University August, 2008 Lectures on Partial Differential Equations and Hilbert Space 1. A linear partial differential equation of physics We begin by considering the simplest
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
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 informationChapter 1. Viscosity and the stress (momentum flux) tensor
Chapter 1. Viscosity and the stress (momentum flux) tensor Viscosity and the Mechanisms of Momentum Transport 1.1 Newton s law of viscosity ( molecular momentum transport) 1.2 Generalization of Newton
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationCourse Syllabus: Continuum Mechanics - ME 212A
Course Syllabus: Continuum Mechanics - ME 212A Division Course Number Course Title Academic Semester Physical Science and Engineering Division ME 212A Continuum Mechanics Fall Academic Year 2017/2018 Semester
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationA NONLOCAL VECTOR CALCULUS, NONLOCAL VOLUME-CONSTRAINED PROBLEMS, AND NONLOCAL BALANCE LAWS
Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company A NONLOCAL VECTOR CALCULUS, NONLOCAL VOLUME-CONSTRAINED PROBLEMS, AND NONLOCAL BALANCE LAWS QIANG DU Department
More information