Saturated Elastic Porous Solids: Incompressible,Compressible and Hybrid Binary Models

Constitutive models for a general binary elastic-porous media are investigated by two complementary approaches. These models include both constituents treated as compressible/incompressible, a compressible solid phase with an incompressible fluid phase (hybrid model of first type), and an incompressible solid phase with a compressible fluid phase (hybrid model of second type). The macroscopic continuum mechanical approach uses evaluation of entropy inequality with the saturation condition always considered as a constraint. This constraint leads to an interface pressure acting in both constituents. Two constitutive equations for the interface pressure, one for each phase, are identified, thus closing the set of field equations. The micromechanical approach shows that the results of Didwania and de Boer can be easily extended to general binary porous media.     

On mechanics of saturated granular materials

A continuum theory of saturated granular materials is formulated. The basic balance laws for the solid phase as well as for the fluid phase are presented. The constitutive equations are derived and the basic equations of motion of the solid and fluid continua are obtained. Several cases of interest, such as incompressible granules saturated with liquids are discussed. It is shown that the theory contains, as its special cases, the Mohr-Coulomb criterion for a granular material as well as Darcy's law of flow through porous media.     

Phase equilibrium in non-fluids and non-fluid mixtures

A criterion of phase equilibrium for mixtures of materials with arbitrary symmetry (e.g. between solid and fluid or two solid mixture phases) is deduced using a rational thermodynamics approach. This criterion, known also as Maxwell relation, is expressed via the difference of chemical potential tensors (Eshelby tensors) on the singular surface dividing the bulk phases.The thermomechanical balance equations, the entropy inequality and the Maxwell relation for phase equilibrium are given first for the case of pure (one-constituent) materials of arbitrary symmetry and then for the case of mixtures (including chemically reacting ones) of arbitrary symmetry.In the special case of fluids it is shown that the chemical potential tensors reduce to the classical scalar chemical potentials and the Maxwell relations to the classical thermodynamic criterion for the phase equilibrium.     

Three Pressures in Porous Media

In a thermodynamic setting for a single phase (usually fluid), the thermodynamically defined pressure, involving the change in energy with respect to volume, is often assumed to be equal to the physically measurable pressure, related to the trace of the stress tensor. This assumption holds under certain conditions such as a small rate of deformation tensor for a fluid. For a two-phase porous medium, an additional thermodynamic pressure has been previously defined for each phase, relating the change in energy with respect to volume fraction. Within the framework of Hybrid Mixture Theory and hence the Coleman and Noll technique of exploiting the entropy inequality, we show how these three macroscopic pressures (the two thermodynamically defined pressures and the pressure relating to the trace of the stress tensor) are related and discuss the physical interpretation of each of them. In the process, we show how one can convert directly between different combinations of independent variables without re-exploiting the entropy inequality. The physical interpretation of these three pressures is investigated by examining four media: a single solid phase, a porous solid saturated with a fluid which has negligible physico-chemical interaction with the solid phase, a swelling porous medium with a non-interacting solid phase, such as well-layered clay, and a swelling porous medium with an interacting solid phase such as swelling polymers.     

A thermodynamic and dynamic subgrid closure model for two‐material cells

A general and robust subgrid closure model for two‐material cells is proposed. The conservative quantities of the entire cell are apportioned between two materials, and then, pressure and velocity are fully or partially equilibrated by modeling subgrid wave interactions. An unconditionally stable and entropy‐satisfying solution of the processes has been successfully found. The solution is valid for arbitrary level of relaxation. The model is numerically designed with care for general materials and is computationally efficient without recourse to subgrid iterations or subcycling in time. The model is implemented and tested in the Lagrange‐remap framework. Two interesting results are observed in 1D tests. First, on the basis of the closure model without any pressure and velocity relaxation, a material interface can be resolved without creating numerical oscillations and/or large nonphysical jumps in the problem of the modified Sod shock tube. Second, the overheating problem seen near the wall surface can be solved by the present entropy‐satisfying closure model. The generality, robustness, and efficiency of the model make it useful in principle in algorithms, such as ALE methods, volume of fluid methods, and even some mixture models, for compressible two‐phase flow computations. Copyright © 2013 John Wiley & Sons, Ltd.     

Hamiltonian Structure and Dynamics of a Neutrally Buoyant Rigid Sphere Interacting with Thin Vortex Rings

In a previous paper, we presented a (noncanonical) Hamiltonian model for the dynamic interaction of a neutrally buoyant rigid body of arbitrary smooth shape with N closed vortex filaments of arbitrary smooth shape, modeled as curves, in an infinite ideal fluid in \mathbbR³\mathbb{R}^3. The setting of that paper was quite general, and the model abstract enough to make explicit conclusions regarding the dynamic behavior of such systems difficult to draw. In the present paper, we examine a restricted class of such systems for which the governing equations can be realized concretely and the dynamics examined computationally. We focus, in particular, on the case in which the body is a smooth sphere. The equations of motion and Hamiltonian structure of this dynamic system, which follow from the general model, are presented. Following this, we impose the constraint of axisymmetry on the entire system and look at the case in which the rings are all circles perpendicular to a common axis of symmetry passing through the center of the sphere. This axisymmetric model, in our idealized framework, is governed by ordinary differential equations and is, relatively speaking, easily integrated numerically. Finally, we present some plots of dynamic orbits of the axisymmetric system.     

Multi‐material closure model for high‐order finite element Lagrangian hydrodynamics

We present a new closure model for single fluid, multi‐material Lagrangian hydrodynamics and its application to high‐order finite element discretizations of these equations 1 . The model is general with respect to the number of materials, dimension and space and time discretizations. Knowledge about exact material interfaces is not required. Material indicator functions are evolved by a closure computation at each quadrature point of mixed cells, which can be viewed as a high‐order variational generalization of the method of Tipton 2 . This computation is defined by the notion of partial non‐instantaneous pressure equilibration, while the full pressure equilibration is achieved by both the closure model and the hydrodynamic motion. Exchange of internal energy between materials is derived through entropy considerations, that is, every material produces positive entropy, and the total entropy production is maximized in compression and minimized in expansion. Results are presented for standard one‐dimensional two‐material problems, followed by two‐dimensional and three‐dimensional multi‐material high‐velocity impact arbitrary Lagrangian–Eulerian calculations. Published 2016. This article is a U.S. Government work and is in the public domain in the USA.     

Neglected transport equations: extended Rankine–Hugoniot conditions and <Emphasis Type="Italic">J</Emphasis> -integrals for fracture

Transport equations in integral form are well established for analysis in continuum fluid dynamics but less so for solid mechanics. Four classical continuum mechanics transport equations exist, which describe the transport of mass, momentum, energy and entropy and thus describe the behaviour of density, velocity, temperature and disorder, respectively. However, one transport equation absent from the list is particularly pertinent to solid mechanics and that is a transport equation for movement, from which displacement is described. This paper introduces the fifth transport equation along with a transport equation for mechanical energy and explores some of the corollaries resulting from the existence of these equations. The general applicability of transport equations to discontinuous physics is discussed with particular focus on fracture mechanics. It is well established that bulk properties can be determined from transport equations by application of a control volume methodology. A control volume can be selected to be moving, stationary, mass tracking, part of, or enclosing the whole system domain. The flexibility of transport equations arises from their ability to tolerate discontinuities. It is insightful thus to explore the benefits derived from the displacement and mechanical energy transport equations, which are shown to be beneficial for capturing the physics of fracture arising from a displacement discontinuity. Extended forms of the Rankine–Hugoniot conditions for fracture are established along with extended forms of J -integrals.     

ON AN INTERDISCIPLINARY TEACHING METHOD OF FLUID MECHANICS AND OTHER MECHANICS RELATED COURSES1)

流体力学是力学和机械等专业本科生的一门专业基础课。由于其涉及数学、物理的概念多、公式复杂、理论抽象而实际应用广,是公认的难教、难学的课程之一。笔者提出了将流体力学与理论力学、材料力学、电动力学等其他力学课程交叉教学的方法,充分利用各门力学课程知识间的联系,使得学生在学习流体力学的同时复习其他力学的知识,通过对比深化理解流体、固体、刚体和电磁场,最终获得对力学学科的整体视野,增强教学效果。     

Shock waves in saturated thermoelastic porous media

The objective of this paper is to develop and present the macroscopic motion equations for the fluid and the solid matrix, in the case of a saturated porous medium, in the form of coupled, nonlinear wave equations for the fluid and solid velocities. The nonlinearity in the equations enables the generation of shock waves. The complete set of equations required for determining phase velocities in the case of a thermoelastic solid matrix, includes also the energy balance equation for the porous medium as a whole, as well as mass balance equations for the two phase. In the special case of a rigid solid matrix, the wave after an abrupt change in pressure propagates only through the fluid.     

Saturated Compressible and Incompressible Porous Solids: Macro- and Micromechanical Approaches

In this paper two complementary approaches are used to describe the mechanical behavior of saturated compressible and incompressible porous solids. The macroscopic investigation is based on the mixture theory, restricted by the volume fraction concept. In the micromechanical approach, a hierarchy of conditionally ensemble averaged fluid and solid phase momentum balance equations are derived for a simple model of quasi-static liquid saturated porous media. The ensemble averaged equations for both the phases agree remarkably well with the macroscopic results. A micromechanical basis for Terzhagi's effective stress concept is presented. In addition, an expression for additional partial solid stress modifying the effective stress principle, to account for deformability of solid materials, is also derived.     

A general solution of the axisymmetric contact problem for biphasic cartilage layers

A unilateral axisymmetric contact problem for articular cartilage layers is considered. The articular cartilages bonded to subchondral bones are modeled as biphasic materials consisting of a solid phase and a fluid phase. It is assumed that the subchondral bones are rigid and shaped like bodies of revolution with arbitrary convex profiles. The obtained closed-form analytical solution is valid over time periods compared with the typical diffusion time and can be used for increasing loading.     

Plane Deformations in Incompressible Nonlinear Elasticity

The substantially general class of plane deformation fields, whose only restriction requires that the angular deformation not vary radially, is considered in the context of isotropic incompressible nonlinear elasticity. Analysis to determine the types of deformations possible, that is, solutions of the governing systems of nonlinear partial differential equations and constraint of incompressibility, is developed in general. The Mooney-Rivlin material model is then considered as an example and all possible solutions to the equations of equilibrium are determined. One of these is interpreted in the context of nonradially symmetric cavitation, i.e., deformation of an intact cylinder to one with a double-cylindrical cavity. Results for general incompressible hyperelastic materials are then discussed. The novel approach taken here requires the derivation and use of a material formulation of the governing equations; the traditional approach employing a spatial formulation in which the governing equations hold on an unknown region of space is not conducive to the study of deformation fields containing more than one independent variable. The derivation of the cylindrical polar coordinate form of the equilibrium equations for the nominal stress tensor (material formulation) for a general hyperelastic solid and a fully arbitrary cylindrical deformation field is also given. This revised version was published online in August 2006 with corrections to the Cover Date.     

功能梯度压电夹层结构中SH波的传播

研究了功能梯度压电上、下半空间和均匀压电层组成的夹层结构中SH波的传播性能,上、下功能梯度半空间的材料性能沿垂直于界面方向以指数函数形式变化。首先推导了SH传播时电弹场的解析解,然后利用界面条件得到了行列式形式的频散方程。基于推导的频散方程,通过数值算例表明了材料性能梯度变化、压电层厚度和材料组合方式对相速度的影响,结果对功能梯度压电材料在声波器件中的应用有参考价值。     

A Model for Flow and Deformation in Unsaturated Swelling Porous Media

A thermomechanical theory for multiphase transport in unsaturated swelling porous media is developed on the basis of Hybrid Mixture Theory (saturated systems can also be modeled as a special case of this general theory). The aim is to comprehensively and non-empirically describe the effect of viscoelastic deformation on fluid transport (and vice versa) for swelling porous materials. Three phases are considered in the system: the swelling solid matrix s, liquid l, and air a. The Coleman–Noll procedure is used to obtain the restrictions on the form of the constitutive equations. The form of Darcy’s law for the fluid phase, which takes into account both Fickian and non-Fickian transport, is slightly different from the forms obtained by other researchers though all the terms have been included. When the fluid phases interact with the swelling solid porous matrix, deformation occurs. Viscoelastic large deformation of the solid matrix is investigated. A simple form of differential-integral equation is obtained for the fluid transport under isothermal conditions, which can be coupled with the deformation of the solid matrix to solve for transport in an unsaturated system. The modeling theory thus developed, which involves two-way coupling of the viscoelastic solid deformation and fluid transport, can be applied to study the processing of biopolymers, for example, soaking of foodstuffs and stress-crack predictions. Moreover, extension and modification of this modeling theory can be applied to study a vast variety of problems, such as drying of gels, consolidation of clays, drug delivery, and absorption of liquids in diapers.     

Two-layer debris mixture flows on arbitrary terrain with mass exchange at the base and the interface

Granular-fluid gravity driven flow down arbitrary topographic terrain is modelled as a two-layer system of a solid-fluid mixture layer overlain by a slurry, consisting of a particle-laden fluid. The lower layer is dynamically treated as a two-phase flow with two constituent mass and momentum balance laws. By contrast, the slurry is described by mass and momentum balances for the mixture as a whole and a diffusive mass balance for the suspended particle phase. At the base, the mixture interacts with the stagnant base by solid-fluid deposition or erosion. At the mixture-slurry interface, solid and fluid mass exchanges are equally taken into account, but the free surface is treated as material and tractionless. The dynamical equations are formulated in three-dimensional form as general balance laws of mass and momentum in each layer. Intrinsic expressions of the jump conditions of mass and momentum are given for the basal and interface surfaces. The field equations are put into dimensionless form and presented relative to topography adjusted coordinates. These equations are further simplified and approximated by a depth-averaging procedure using an order parameter ${\varepsilon = H/L}$ , where H and L are typical thickness and length scales of the gravity current. Detailed proposals are worked out for the parameterizations of the solid and fluid mass flows across the basal surface and layer interface.     

On the Equations of Capillarity

Some conceptual ambiguities in the derivation of the equations of capillarity on the basis of the principle of virtual work are addressed, and hypotheses are proposed toward obtaining a physically correct characterization in general circumstances. It is shown that under the hypotheses, the classical equations of capillarity for an interface of an incompressible fluid with a fluid of negligible density can be obtained on the basis of global phenomenological reasoning, without recourse to consideration of intermolecular attractions. More generally, the procedure is applied to derive the specific equations arising from a compressible fluid configuration with idealized pressure-density relationship in a capillary tube, and a general necessary condition for existence of a solution is established. It is shown that for symmetric domains, the condition is also sufficient for existence of a unique symmetric solution.     

Homogenized modeling for vascularized poroelastic materials

A new mathematical model for the macroscopic behavior of a material composed of a poroelastic solid embedding a Newtonian fluid network phase (also referred to as vascularized poroelastic material), with fluid transport between them, is derived via asymptotic homogenization. The typical distance between the vessels/channels (microscale) is much smaller than the average size of a whole domain (macroscale). The homogeneous and isotropic Biot’s equation (in the quasi-static case and in absence of volume forces) for the poroelastic phase and the Stokes’ problem for the fluid network are coupled through a fluid-structure interaction problem which accounts for fluid transport between the two phases; the latter is driven by the pressure difference between the two compartments. The averaging process results in a new system of partial differential equations that formally reads as a double poroelastic, globally mass conserving, model, together with a new constitutive relationship for the whole material which encodes the role of both pore and fluid network pressures. The mathematical model describes the mutual interplay among fluid filling the pores, flow in the network, transport between compartments, and linear elastic deformation of the (potentially compressible) elastic matrix comprising the poroelastic phase. Assuming periodicity at the microscale level, the model is computationally feasible, as it holds on the macroscale only (where the microstructure is smoothed out), and encodes geometrical information on the microvessels in its coefficients, which are to be computed solving classical periodic cell problems. Recently developed double porosity models are recovered when deformations of the elastic matrix are neglected. The new model is relevant to a wide range of applications, such as fluid in porous, fractured rocks, blood transport in vascularized, deformable tumors, and interactions across different hierarchical levels of porosity in the bone.     

An elastoplasticity model in porous media theories

In this article, porous media theories are referred to as mixture theories extended by the well-known concept of volume fractions. This approach implies the diverse field functions of both the porous solid matrix and the pore fluid to be represented by average functions of the macroscale.The present investigations are based on a binary model of incompressible constituents, solid skeleton, and pore liquid, where, in the constitutive range, use is made of the second-grade character of general heterogeneous media. Within the framework of geometrically finite theories, the paper offers a set of constitutive equations for the solid matrix, the viscous pore liquid and the different interactions between the constituents. The constitutive model applies to saturated as well as to empty solid materials, taking into account the physical nonlinearities based on elasto-plastic solid deformations. In particular, the constitutive model concentrates on granular materials like soil or concrete, where the elastic deformations are usually small and the plastic range is governed by kinematically hardening properties.