A Multiscale Theory of Swelling Porous Media: II. Dual Porosity Models for Consolidation of Clays Incorporating Physicochemical Effects

A three-scale theory of swelling clay soils is developed which incorporates physico-chemical effects and delayed adsorbed water flow during secondary consolidation. Following earlier work, at the microscale the clay platelets and adsorbed water (water between the platelets) are considered as distinct nonoverlaying continua. At the intermediate (meso) scale the clay platelets and the adsorbed water are homogenized in the spirit of hybrid mixture theory, so that, at the mesoscale they may be thought of as two overlaying continua, each having a well defined mass density. Within this framework the swelling pressure is defined thermodynamically and it is shown to govern the effect of physico-chemical forces in a modified Terzaghi's effective stress principle. A homogenization procedure is used to upscale the mesoscale mixture of clay particles and bulk water (water next to the swelling mesoscale particles) to the macroscale. The resultant model is of dual porosity type where the clay particles act as sources/sinks of water to the macroscale bulk phase flow. The dual porosity model can be reduced to a single porosity model with long term memory by using Green's functions. The resultant theory provides a rational basis for some viscoelastic models of secondary consolidation.     

Coupled Solvent and Heat Transport of a Mixture of Swelling Porous Particles and Fluids: Single Time-Scale Problem

A three-spatial scale, single time-scale model for both moisture and heat transport is developed for an unsaturated swelling porous media from first principles within a mixture theoretic framework. On the smallest (micro) scale, the system consists of macromolecules (clay particles, polymers, etc.) and a solvating liquid (vicinal fluid), each of which are viewed as individual phases or nonoverlapping continua occupying distinct regions of space and satisfying the classical field equations. These equations are homogenized forming overlaying continua on the intermediate (meso) scale via hybrid mixture theory (HMT). On the mesoscale the homogenized swelling particles consisting of the homogenized vicinal fluid and colloid are then mixed with two bulk phase fluids: the bulk solvent and its vapor. At this scale, there exists three nonoverlapping continua occupying distinct regions of space. On the largest (macro) scale the saturated homogenized particles, bulk liquid and vapor solvent, are again homogenized forming four overlaying continua: doubly homogenized vicinal fluid, doubly homogenized macromolecules, and singly homogenized bulk liquid and vapor phases. Two constitutive theories are developed, one at the mesoscale and the other at the macroscale. Both are developed via the Coleman and Noll method of exploiting the entropy inequality coupled with linearization about equilibrium. The macroscale constitutive theory does not rely upon the mesoscale theory as is common in other upscaling methods. The energy equation on either the mesoscale or macroscale generalizes de Vries classical theory of heat and moisture transport. The momentum balance allows for flow of fluid via volume fraction gradients, pressure gradients, external force fields, and temperature gradients.     

Plane Waves and Uniqueness Theorems in the Coupled Linear Theory of Elasticity for Solids with Double Porosity

This paper concerns with the coupled linear dynamical theory of elasticity for solids with double porosity. Basic properties of plane harmonic waves are established. Radiation conditions of regular vectors are given. Basic internal and external boundary value problems (BVPs) of steady vibrations are formulated, and finally, uniqueness theorems for regular (classical) solutions of these BVPs are proved.     

Interfaces of Phase Transition and Disappearance and Method of Negative Saturation for Compositional Flow with Diffusion and Capillarity in Porous Media

The specific case of interfaces separating a single-phase fluid and a two-phase continuum appears in the theory of compositional flow through porous media. They are usually called the interfaces of phase transition (PT-interfaces) or the interfaces of phase disappearing (PD-interfaces). The principle of equivalence is proved which shows that a single-phase multi-component fluid may be replaced by an equivalent fictitious two-phase fluid having specific properties. The equivalent properties are such that the extended saturation of a fictitious phase is negative. This principle enables us to develop the uniform system of two-phase equations in the overall domain in terms of the extended saturation (the NegSat model), and to apply the direct numerical simulation. In the case with diffusion, the uniform NegSat model contains a new term proportional to the gradient of saturation in the relation for flow velocity. The canonical NegSat model represents a transport equation with discontinuous nonlinearities. The qualitative analysis of this model shows that the PT-interfaces represent the shocks of the extended saturation, or, in some cases, can transform into weak shocks. The diffusion and capillarity do not destroy necessarily the shocks, but change their velocity. The analytical technique is developed which allows capturing PT-shocks. The method is illustrated by several examples of miscible gas injection in oil reservoir. In two-dimensional case, the effects of multiple shock collisions in heterogeneous media are automatically modeled. In the case of immiscible fluids and a classic interface, the suggested method converges to the VOF-method.     

The Flow Problem of Fluids Flow in a Fractal Reservoir with Double Porosity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｆｌｏｗｔｈｅｏｒｙａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｏｆｆｌｕｉｄｓｆｌｏｗｉｎａｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒｈａｖｅｃｏｎｔｉｎｕａｌｌｙｇｏｎｅｄｅｅｐｉｎｔｏｓｔｕｄｙｓｉｎｃｅＣｈａｎｇａｎｄＹｏｒｔｓｏｓ[1]ｂｕｉｌｔｔｈｅｆｌｏｗｍｏｄｅｌｏｆｆｌｕｉｄｔｈｒｏｕｇｈａｆｒａｃｔａｌｒｅｓｅｒｖｏｉｒ.ＴＯＮＧＤｅｎｇ＿ｋｅ[2 ]ｐｒｅｓｅｎｔｅｄｔｈｅｅｘａｃｔｓｏｌｕｔｉｏｎａｎｄｉｔｓｐｒｅｓｓｕｒｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓｆｏｒｔｈｅｖａ…     

Dynamic Model of a Periodic Medium with Double Porosity

The paper presents a unified mathematical approach for describing the dynamic stressstrain state of mechanical structures from heterogeneous materials possessing a double coupled system of pore channels filled with fluid. New dynamic equations describing the oscillations of poroelastic systems based on the developed model of a continuous medium with additional degrees of freedom in the form of various pressures of the components constituting the liquid phase of the material are obtained. The equations and the method of obtaining them have a greater degree of generalization than those encountered in the literature. Theoretical results can be used to study the propagation of vibrations in fractured geological rocks saturated with liquid, to develop technical systems of new structural materials with a porous structure, for the analysis of micro streams of fluid in the hierarchical system of microporous bone tissue.     

Swelling and instability of a gel annulus

We study the swelling of a gel annulus attached to a rigid core when it is immersed in a solvent.For equilibrium states,the free-energy function of the gel can be converted into a strain energy function,and as a result the gel can be treated as a compressible hyperelastic material.Asymptotic methods are used to study the inhomogeneous swelling in order to obtain the leading-order solution.Some analytical insights are then deduced.Because of the compressive hoop stress in this state,at some stage instability can occur,leading to wrinkles in the gel.An incremental deformation theory in nonlinear elasticity is used to conduct a linear bifurcation analysis for understanding such instability.More specifically,the critical loading for the onset of a wrinkled state is obtained.Detailed discussions on the behaviors of various physical quantities in this critical state are given.It is found that the critical mode number,while insensitive to the material parameters,is greatly influenced by the ratio of the inner and outer radii of the gel.Also,an interesting finding is that the critical swelling ratio is an increasing function of this geometrical parameter,which implies that a thin annulus is more likely to be unstable than a thick one.     

Drag Force Calculations in Polydisperse DEM Simulations with the Coarse-Grid Method: Influence of the Weighting Method and Improved Predictions Through Artificial Neural Networks

Transport in Porous Media - Several methods are employed for drag force calculations with the discrete element method depending on the desired accuracy and the number of particles involved. For...     

Torsional Swelling of a Hyperelastic Tube with Helically Wound Reinforcement

This paper examines the combination of radial deformation with torsion for a circular cylindrical tube composed of a transversely isotropic hyperelastic material subject to finite deformation swelling. The stored energy function involves separate matrix and fiber contributions such that the fiber contribution is minimized when the fiber direction is at a natural length. This natural length is not affected by the swelling. Hence swelling preferentially expands directions that are orthogonal to the fiber. The swelling itself is described via a swelling field that prescribes the local free volume at each location in the body. Such a treatment is a relatively simple generalization of the conventional incompressible theory. The direction of transverse isotropy associated with the fiber reinforcement is described by a helical winding about the tube axis. The swelling induced preferential expansion orthogonal to this direction induces the torsional aspect of the deformation. For a specific class of strain energy functions we find that the twist increases with swelling and approaches a limiting asymptotic value as the swelling becomes large. The fibers reorient such that fibers at the inner portion of the tube assume a more circumferential orientation whereas, at least for small and moderate swelling, the fibers in the outer portion of the tube assume a more axial orientation. For large swelling the fibers in the outer portion of the tube reorient beyond the axial orientation, and so are described by helices with orientation in the opposite sense to that in the reference configuration.      

Numerical Analysis of Thermally Coupled Flow Problems with Interfaces and Phase-change Effects

In this work a fixed mesh finite element approach is presented to solve thermally coupled flow problems including moving interfaces between immiscible fluids and phase-change effects. The weak form of the full incompressible Navier-Stokes equations is obtained using a generalized streamline operator (GSO) technique that enables the use of equal order interpolation of the primitive variables of the problem: velocity, pressure and temperature. The interfaces are defined with a mesh of marker points whose motion is obtained applying a Lagrangian scheme. Moreover, a temperature-based formulation is considered to describe the phase-change phenomena. The proposed methodology is used in the analysis of a filling of a step mould and a gravity-driven flow of an aluminium alloy in an obstructed vertical channel.     

Effect of sample size on the response of DEM samples with a realistic grading

This paper shows that for DEM simulations of triaxial tests using samples with a grading that is repre- sentative of a real soil, the sample size significantly influences the observed material response. Four DEM samples with identical initial states were produced： three cylindrical samples bounded by rigid wails and one bounded by a cubical periodic cell, When subjected to triaxial loading, the samples with rigid boundaries were more dilative, stiffer and reached a higher peak stress ratio than the sample enclosed by periodic boundaries. For the rigid-wall samples, dilatancy increased and stiffness decreased with increasing sample size, The periodic sample was effectively homogeneous, The void ratio increased and the contact density decreased close to the rigid walls, This heterogeneity reduced with increasing sample size. The positions of the critical state lines （CSLs） of the overall response in e-log p＇ space were sensitive to the sample size, although no difference was observed between their slopes. The critical states of the interior regions of the rigid-wall-bounded samples approached that of the homogeneous periodic sample with increasing sample size. The ultimate strength of the material at the critical state is independent of sample size.     

Coupled Diffusion-Dissolution Around a Fracture Channel: The Solute Congestion Phenomenon

The paper studies the coupled diffusion-dissolution process in reactive porous media, separated by a fracture channel. An aggressive solute, corresponding for e.g., to a complete demineralization that dissolves the solid skeleton of the surrounding porous material, is prescribed at the inlet of the fracture. By means of asymptotic dimensional analysis it is shown that for large times the diffusion length in the fracture develops with the quadratic root of time. In comparison with the 1D-Stefan Problem, in which the dissolution front evolves with the square root of time, this indicates that the overall solute evacuation through the structure slows down in time. This phenomenon is referred to as a diffusive solute congestion in the fracture. This asymptotic behavior is confirmed by means of model-based simulation, and the relevant material parameters, related to only the chemical equilibrium condition, are identified. The obtained results suggest that the presence of a small crack does not significantly increase the propagation of the dissolution front in the porous bulk, and hence the overall chemical degradation of the porous material. The same applies to other diffusion induced demineralization, mineralization, sorption and melting processes, provided that the convective transport of the solute in the crack is small in comparison with the solute diffusion. The result is relevant for several problems in durability mechanics: calcium leaching of concrete in nuclear waste containment, mineralization and demineralization in bone remodeling, chloride penetration, etc.     

Studies on Bifurcation and Chaos of a String-Beam Coupled System with Two Degrees-of-Freedom

In this paper, research on nonlinear dynamic behavior of a string-beam coupled system subjected to parametric and external excitations is presented. The governing equations of motion are obtained for the nonlinear transverse vibrations of the string-beam coupled system. The Galerkin's method is employed to simplify the governing equations to a set of ordinary differential equations with two degrees-of-freedom. The case of 1:2 internal resonance between the modes of the beam and string, principal parametric resonance for the beam, and primary resonance for the string is considered. The method of multiple scales is utilized to analyze the nonlinear responses of the string-beam coupled system. Based on the averaged equation obtained here, the techniques of phase portrait, waveform, and Poincare map are applied to analyze the periodic and chaotic motions. It is found from numerical simulations that there are obvious jumping phenomena in the resonant response–frequency curves. It is indicated from the phase portrait and Poincare map that period-4, period-2, and periodic solutions and chaotic motions occur in the transverse nonlinear vibrations of the string-beam coupled system under certain conditions. An erratum to this article is available at .     

Application of a screw conveyor with axial tilt blades on a shearer drum and investigation of conveying performance based on DEM

Screw conveyors are widely employed in industrial fields for conveying bulk materials. The shearer drum which uses the screw conveying principle is responsible for excavation and conveying coal particles onto the chain conveyor. Screw conveyor performance is affected by potential factors, such as the blade axial tilt angle and style, core shaft form and diameter. The effect of blade axial tilt angle on the conveying performance was investigated with the help of DEM. In the case of the screw conveyor, the mass flow rate, and particle axial velocity increased with increasing positive axial tilt angle, and declined with increasing negative axial tilt angle. In the case of the drum, the mass flow rate, particle axial velocity, and loading rate first increased and then decreased with increasing positive axial tilt angle, and decreased with increasing negative axial tilt angle. These results can be considered as a benchmark for screw conveyor and drum structural designs with axial tilt screw blades.     

Application of a screw conveyor with axial tilt blades on a shearer drum and investigation of conveying performance based on DEM

Screw conveyors are widely employed in industrial fields for conveying bulk materials.The shearer drum which uses the screw conveying principle is responsible f...     

多重网格法与Newton-Raphson法耦合求线弹流数值解

利用多重网格算法,在最粗网格上,采用Newton-Raphson迭代法模式精解线接触弹流非线性方程组,充分利用了多重网格法与Newton-Raphson法各自的优点,计算实践表明,求解弹流问题的数值解过程在收敛与稳定性方面均有较大改善,且有相当宽的载荷参数适用范围。