Equations of Continuum Dynamics

It is important to emphasize that Q. itself is time-dependent since the material volume element of interest is undergoing deformation. Note that in our statement of the balance of linear momentum we have introduced notation for our description of the time derivative which differs from the conventional time derivative, indicating that we are evaluating the material time derivative. The material time derivative evaluates the time rate of change of a quantity for a given material particle. Explicitly, we write 

Though the integral form of linear momentum balance written above as eqn (2.26) serves as the starting point of many theoretical developments, it is also convenient to cast the equation in local form as a set of partial differential equations. The idea is to exploit the Reynolds transport theorem in conjunction 

The logic behind this development is related to the equation of continuity which asserts the conservation of mass and may be written as 

If we now recognize that the choice of our subregion Q was arbitrary we are left with 

This equation is the main governing equation of continuum mechanics and presides over problems ranging from the patterns formed by clouds to the evolution of dislocation structures within metals. Note that as yet no constitutive assumptions have been advanced rendering this equation of general validity. We will find later that for many problems of interest in the modeling of materials it will suffice to consider the static equilibrium version of this equation in which it is assumed that the body remains at rest. In that case, eqn (2.32) reduces to 

---

What this equation tells us is that a particular state of stress is nothing more than a linear combination (albeit perhaps a tedious one) of the entirety of components of the strain tensor. The tensor Cijn is known as the elastic modulus tensor or stiffness and for a linear elastic material provides nearly a complete description of the material properties related to deformation under mechanical loads. Eqn (2.52) is our first example of a constitutive equation and, as claimed earlier, provides an explicit statement of material response that allows for the emergence of material specificity in the equations of continuum dynamics as embodied in eqn (2.32). In particular, if we substitute the constitutive statement of eqn (2.52) into eqn (2.32) for the equilibrium case in which there are no accelerations, the resulting equilibrium equations for a linear elastic medium are given by... 

From the standpoint of the continuum simulation of processes in the mechanics of materials, modeling ultimately boils down to the solution of boundary value problems. What this means in particular is the search for solutions of the equations of continuum dynamics in conjunction with some constitutive model and boundary conditions of relevance to the problem at hand. In this section after setting down some of the key theoretical tools used in continuum modeling, we set ourselves the task of striking a balance between the analytic and numerical tools that have been set forth for solving boundary value problems. In particular, we will examine Green function techniques in the setting of linear elasticity as well as the use of the finite element method as the basis for numerical solutions. 

Cauchy Tetrahedron and Equilibrium In this problem, use the balance of linear momentum to deduce the Cauchy stress by writing down the equations of continuum dynamics for the tetrahedral volume element shown in fig. 2.7. Assign an area AS to the face of the tetrahedron with normal n and further denote the distance from the origin to this plane along the direction why h. Show that the dynamical equations may be written as... 

This chapter is organized into two main parts. To give the reader an appreciation of real fluids, and the kinds of behaviors that it is hoped can be captured by CA models, the first part provides a mostly physical discussion of continuum fluid dynamics. The basic equations of fluid dynamics, the so-called Navier-Stokes equations, are derived, the Reynolds Number is defined and the different routes to turbulence are described. Part I also includes an important discussion of the role that conservation laws play in the kinetic theory approach to fluid dynamics, a role that will be exploited by the CA models introduced in Part II. 

There are two levels, discrete particle level and continuum level, for describing and modeling of the macroscopic behaviors of dilute and condensed matters. The physics laws concerning the conservation of mass, momentum, and energy in motion, are common to both levels. For simple dilute gases, the Boltzmann equation, as shown below, provides the governing equation of gas dynamics on the discrete particle level... 

More or less as a spin-off result of the foregoing analysis determining the transport coefficients, a rigorous procedure deriving the governing equations of fluid dynamics from first principles was established. It is stressed that in classical fluid dynamics the continuum hypothesis is used so that the governing... 

The numerical solution of the equations of continuum turbulent flow with effective viscosity coefficient he dynamic coefficient of molecular viscosity)... 

At the microscale, dislocation dynamics simulations implement the equations of continuum elasticity theory to track the motion and interaction of individual dislocations under an applied stress, leading to the development of a dislocation microstructure and single-crystal plastic deformation. In our multiscale modeling... 

Continuum mechanics, one of our most successful physical theories, is readily applicable to the process industries. In continuum mechanics, the existence of molecules is ignored, and matter is treated as a continuous medium. The continuum hypothesis is valid, provided the equations of continuum mechanics are applied at sufficiently large length scales and time scales that the properties of individual molecules are not noticed. The mapping of the laws of mass, momentum, and energy conservation to the continuum results in field equations that describe the dynamics of the continuum. These field equations, variously known as the equations of motion, the equations of change, or simply the conservation equations, are nonlinear, partial differential equations that can be solved, in principle, when combined with the appropriate constitutive information and boundary conditions. 

Continuum mechanics is the mechanical analog of classical electrodynamics, in which a set of field equations (Maxwell s equations) describe the dynamics of the relevant variables of the electrical and magnetic fields. Whereas Maxwell s equations are linear unless the constitutive behavior is nonlinear, the equations of continuum mechanics are nonlinear, regardless of the constitutive behavior of the materials of interest. The inherent nonlinearity of the conservation equations, which is due to convective transport of momentum, energy, and... 

A mathematical description of continuum fluid dynamics can proceed from two fundamentally different points of view (1) A macroscopic, or top-down ([hass87], [hass88]), approach, in which the equations of motions are derived as the most... 

The first equation gives the diserete version of Newton s equation the second equation gives energy c onservation. We make two comments (1) Notice that while energy eouseivation is a natural consequence of Newton s equation in continuum mechanics, it becomes an independent property of the system in Lee s discrete mechanics (2) If time is treated as a conventional parameter and not as a dynamical variable, the discretized system is not tiine-translationally invariant and energy is not conserved. Making both and t , dynamical variables is therefore one way to sidestep this problem. 

Steady state models of the automobile catalytic converter have been reported in the literature 138), but only a dynamic model can do justice to the demands of an urban car. The central importance of the transient thermal behavior of the reactor was pointed out by Vardi and Biller, who made a model of the pellet bed without chemical reactions as a onedimensional continuum 139). The gas and the solid are assumed to have different temperatures, with heat transfer between the phases. The equations of heat balance are ... 

In literature, some researchers regarded that the continuum mechanic ceases to be valid to describe the lubrication behavior when clearance decreases down to such a limit. Reasons cited for the inadequacy of continuum methods applied to the lubrication confined between two solid walls in relative motion are that the problem is so complex that any theoretical approach is doomed to failure, and that the film is so thin, being inherently of molecular scale, that modeling the material as a continuum ceases to be valid. Due to the molecular orientation, the lubricant has an underlying microstructure. They turned to molecular dynamic simulation for help, from which macroscopic flow equations are drawn. This is also validated through molecular dynamic simulation by Hu et al. [6,7] and Mark et al. [8]. To date, experimental research had "got a little too far forward on its skis however, theoretical approaches have not had such rosy prospects as the experimental ones have. Theoretical modeling of the lubrication features associated with TFL is then urgently necessary. 

This section provides an alternative measurement for a material parameter the one in the ensemble averaged sense to pave the way for usage of continuum theory from a hope that useful engineering predictions can be made. More details can be found in Ref. [15]. In fact, macroscopic flow equations developed from molecular dynamics simulations agree well with the continuum mechanics prediction (for instance. Ref. [16]). 

Crystals lack some of the dynamic complexity of solutions, but are still a challenging subject for theoretical modeling. Long-range order and forces in crystals cause their spectrum of vibrational frequencies to appear more like a continuum than a series of discrete modes. Reduced partition function ratios for a continuous vibrational spectrum can be calculated using an integral, rather than the hnite product used in Equation (3) (Kieffer 1982),... 

When considering boundary conditions, a useful dimensionless hydrodynamic number is the Knudsen number, Kn = X/L, the ratio of the mean free path length to the characteristic dimension of the flow. In the case of a small Knudsen number, continuum mechanics will apply, and the no-slip boundary condition assumption is valid. In this formulation of classical fluid dynamics, the fluid velocity vanishes at the wall, so fluid particles directly adjacent to the wall are stationary, with respect to the wall. This also ensures that there is a continuity of stress across the boundary (i.e., the stress at the lower surface—the wall—is equal to the stress in the surface-adjacent liquid). Although this is an approximation, it is valid in many cases, and greatly simplifies the solution of the equations of motion. Additionally, it eliminates the need to include an extra parameter, which must be determined on a theoretical or experimental basis. 

Solids undergoing martensitic phase transformations are currently a subject of intense interest in mechanics. In spite of recent progress in understanding the absolute stability of elastic phases under applied loads, the presence of metastable configurations remains a major puzzle. In this overview we presented the simplest possible discussion of nucleation and growth phenomena in the framework of the dynamical theory of elastic rods. We argue that the resolution of an apparent nonuniqueness at the continuum level requires "dehomogenization" of the main system of equations and the detailed description of the processes at micro scale. 

We now modify the 2D continuum equations of step motion, Eqs. (7) and (8), in order to study some aspects of the dynamics of faceting. We assume the system is in the nucleation regime where the critical width Wc is much larger than the average step spacing In the simplest approximation discussed here, we incorporate the physics of the two state critical width model into the definition of the effective interaction term V(w) in Eq. (2), which in turn modifies the step chemical potential terms in Eqs. (7) and (8). Again we set V(w) = w/ l/w) as in Eq. (4) but now we use the /from Eq. (10) that takes account of reconstruction if a terrace is sufficiently wide. Note that this use of the two state model to describe an individual terrace with width w is more accurate than is the use of Eq. (10) to describe the properties of a macroscopic surface with average slope s = Mw. 

The proposed mechanisms of models to explain the drag reduction phenomenon are based on either a molecular approach or fluid dynamical continuum considerations, but these models are mainly empirical or semi-empirical in nature. Models constructed from the equations of motion (or energy) and from the constitutive equations of the dilute polymer solutions are generally not suitable for use in engineering applications due to the difficulty of placing numerical values on all the parameters. In the absence of a more generally accurate model, semi-empirical ones remain the most useful for applications. 

Continuum Dynamics. In this appruach, fluid properties, such as velocity, density, pressure, temperature, viscosity, and conductivity, among others, arc assumed to be physically meaningful functions of three spatial variables t. . 1 . and. n. and lime i. Nonlinear partial differential equations are set up to relate these variables. Such equations have nil general solutions even for the most restrictive boundary conditions. Bui solutions are carried out for very idealized flows. Couetle flow is one of these. See Fig. I. 

When motion of the fluid consists of only small fluctuations about a state of near-rest, Lhe continuum equations are linearized by neglecting nonlinear terms and they become lhe equalions of acoustics. A large variety of fluid motions are described as sound waves when the small-motion or acoustic description can be used, the principle of superposition is valid. This powerful principle allows addition of simple simultaneous motions to represent a more complex motion, such as the sound reaching lhe audience from the instruments of a symphony orchestra. The superposition principle does not apply to large-scale (nonacoustical) motions, and the subject of fluid dynamics (in distinction from acoustics) treats nonlinear flows. i.e.. those that cannot be described as superpositions of other flows. 

The sets of equations are solved by the assumption of periodic waves and, by expansion in powers of the wave number, a relation is found for the limiting case of long waves so that the elements of the dynamical matrix elastic constants of the continuum. It is also possible to derive the Raman frequencies from the lattice dynamics analysis but this does not seem to have been done for polymer crystals, though they have been derived for example, for NaCl and for diamond. 

The Eulerian continuum approach, based on a continuum assumption of phases, provides a field description of the dynamics of each phase. The Lagrangian trajectory approach, from the study of motions of individual particles, is able to yield historical trajectories of the particles. The kinetic theory modeling for interparticle collisions, extended from the kinetic theory of gases, can be applied to dense suspension systems where the transport in the particle phase is dominated by interparticle collisions. The Ergun equation provides important flow relationships, which are useful not only for packed bed systems, but also for some situations in fluidized bed systems. 

Now that we have settled on a model, one needs to choose the appropriate algorithm. Three methods have been used to study polymers in the continuum Monte Carlo, molecular dynamics, and Brownian dynamics. Because the distance between beads is not fixed in the bead-spring model, one can use a very simple set of moves in a Monte Carlo simulation, namely choose a monomer at random and attempt to displace it a random amount in a random direction. The move is then accepted or rejected based on a Boltzmann weight. Although this method works very well for static and dynamic properties in equilibrium, it is not appropriate for studying polymers in a shear flow. This is because the method is purely stochastic and the velocity of a mer is undefined. In a molecular dynamics simulation one can follow the dynamics of each mer since one simply solves Newton s equations of motion for mer i,... 

The equations of motion (6) or (10) do not per se describe irreversible dynamics. This is only possible if the number N of the bath oscillators tends toward infinity, with their angular frequencies forming a continuum in this limit. 

Note that for the asymptotic equations of Eqns. (2) and (3) to be valid, r

characteristic length, and is normally the crack length or the remaining ligament, whichever is the smaller, of a fracture specimen. Also, the above asymptotic equations are not valid for an orthotropic elastic continuum, such as a ceramic fiber/ceramic matrix composite. While the static crack tip state for an orthotropic elastic continuum has been derived, to the author s knowledge, no dynamic counterpart is available to date. Nevertheless, the above crack tip state should be applicable to particulate/whisker-filled ceramic matrix composites which macroscopically behave like an isotropic homogeneous continuum. 

---