Mathematical Models for Metal Matrix Composites Manufacturing

Infiltration by a pure liquid metal into a rigid preform in the presence of partial solidification is analyzed. The mathematical model which is developed takes into account the coupling between the flow through the porous medium and the thermal effects. The model has a two-scale formulation: phase change occurs at the microscopic scale, while macroscopic flow and heat transfer involve large scale quantities.     

Effect of flows between formations and capillary forces on the process of oil displacement in a stratified inhomogeneous bed

The results of a numerical investigation of the process of oil displacement in a stratified inhomogeneous formation on the basis of the two-phase flow model with account for capillary forces are presented. It is shown that in many cases the vertical inhomogeneity of oil reservoirs may not be a cause of nonuniform displacement and the non-recovery of large oil reserves by the time of water breakthrough to the extraction surface. The action of the capillary forces is an additional factor leading to equalization of the water propagation front in the inhomogeneous formation, water breakthrough delay, and intensification of the mass transfer between the layers with different permeabilities. Analysis of the contribution of the interlayer flows to the water flooding of low-permeability formation intervals calls into question the practicability of blocking high-permeability inclusions in the neighborhood of pumping wells.     

On Mathematical Models of Average Flow in Heterogeneous Formations

We consider a general model of transient flow in media of random conductivity and storativity. The flow is driven by the spatially distributed source function (x, t) and the initial head distribution h ₀(x). The function models sources and wells and can be deterministic, random or a sum of both. The deterministic source function corresponds to singularities of deterministic strength, whereas the random models the head boundary condition. In the latter case, is shown to be proportional to the hydraulic conductivity. The aim of the study is to analyze the feasibility of averaging the flow equations and of developing the mathematical model of average flow (AFM) without solving problems in detail. It is shown that the problem of averaging is reduced to deriving two constitutive equations. The first equation, the effective Darcy's law (EDL) stems from averaging Darcy's law at local scale. The second one is related to the medium ability to store a fluid and expresses the correlation between the storativity and head in terms of the mean head. Both relationships are required to be completely determined by the medium structure (conductivity and storativity statistical properties) and independent of the flow configuration (functions and h ₀). We show that if one of the constitutive equations exists, the same is true respective to the second. This reduces the problem of averaging to the classic one of deriving the EDL. For steady flows the EDL is shown to exist for flows driven by sources (wells) of either deterministic flux or head boundary conditions. No EDL can be derived if both types of sources are present in the flow domain. For unsteady flows the EDL does not exist if the initial head correlates with the medium properties. For uncorrelated initial head distribution, its random residual (due to the measurement errors and scarcity of the data) has no impact on the EDL and is immaterial. For deterministic h ₀, the only case for which the EDL exists is the flow by sources of deterministic discharge. For sources of given head boundary condition the EDL can be derived only for uniform initial head distribution. For all other cases, the EDL does not exist. The results of the study are not limited by usually adopted assumptions of weak heterogeneity and of stationarity of the formation random properties.     

Fluid penetration effects in porous media contact

The contact boundary conditions at the interface between two fluid-saturated porous bodies are derived. The general derivation is performed within the well-founded framework of the Theory of Porous Media (TPM) based on the constituent balance relations of mass, momentum, and energy accounting for finite discontinuities at the contact surface. Particular attention is drawn to the effects associated with the interstitial fluid flux across the interface. The derived contact conditions include two kinematic continuity conditions for the solid velocity and the fluid seepage velocity as well as two jump conditions for the effective solid stress and the pore-fluid pressure. As an application, the common case of biphasic porous media contact proceeding from materially incompressible constituents and inviscid fluid properties is discussed in detail.      

Mathematical modeling of methane flow in a borehole coal mining system

Safety in coal mining is greatly increased by the drainage of the methane content of coal seams through boreholes, simultaneously producing significant energy. The design of suitable drainage technology is based on the mathematical modeling of methane flow in coal seams. In the calculation of the methane pressure, the new mathematical model presented in this paper considers both the sorption phenomenon of methane depending upon the methane pressure and the fact that the variation in methane pressure can create a change in the stress condition of the rock and, as a consequence of this, a change in the permeability of the coal. The new mathematical model can be used for the numerical simulation of the flow processes in coal seams and methane drainage technology can be designed more accurately.     

Permeability Tensor for Heterogeneous Porous Medium of Fibre Type

A theoretical model which allows us to determine the permeability of a fibrous porous medium is proposed. Fibres are assumed to be parallel and nonuniform in space and material with a low volume fraction of fibres is considered. The model includes two geometric parameters: the diameter of fibres and the diameter of caverns or fissures inside the bundle of fibres. The tensor of permeability of the porous medium is determined based upon a generalized cell model. The components of permeability tensor depend on two parameters which are determined using experimental data and least-squares approximation. The influence of the geometric parameters on components of permeability tensor is discussed.     

Mathematical model of two-phase fluid nonlinear flow in low-permeability porous media with applications

IntroductionItisasuccessfulexampleinadevelopmentstoryofscienceandtechnologyformechanicsoffluidsinporousmediatocombinewithengineeringtechnology .Fieldsinfluencedbythemechanicsinvolveddevelopmentofoil_gasandgroundwaterresources,controlonseawaterintrusionandsubsidenceandgeologichazards,geotechnicalengineeringandbioengineering ,andairlineindustry[1～ 7].Aproblemonnonlinearflowinlow_permeabilityporousmediaisbutonlyabasiconeindifferentkindsofengineeringfields,butalsooneoffrontlineresearchfieldsofmod…     

Review on Scale Dependent Characterization of the Microstructure of Porous Media

The paper discusses local porosity theory and its relation with other geometric characterization methods for porous media such as correlation functions and contact distributions. Special emphasis is placed on the charcterization of geometric observables through Hadwigers theorem in stochastic geometry. The four basic Minkowski functionals are introduced into local porosity theory, and for the first time a relationship is established between the Euler characteristic and the local percolation probabilities. Local porosity distributions and local percolation probabilities provide a scale dependent characterization of the microstructure of porous media that can be used in an effective medium approach to predict transport.     

A hyperbolic mathematical modeling for describing the transition saturated/unsaturated in a rigid porous medium

This work proposes a mathematical model to study the filling up of an unsaturated porous medium by a liquid identifying the transition from unsaturated to saturated flow and allowing a small super saturation. As a consequence the problem remains hyperbolic even when saturation is reached. This important feature enables obtaining numerical solution for any initial value problem and allows employing Glimm’s scheme associated with an operator splitting technique for treating drag and viscous effects. A mixture theory approach is used to build the mechanical model, considering a mixture of three overlapping continuous constituents: a solid (porous medium), a liquid (Newtonian fluid) and a very low-density gas (to account for the mixture compressibility). The constitutive assumption proposed for the pressure gives rise to a continuous function of the fluid fraction. The complete solution of the Riemann problem associated with the system of conservation laws, as well as four examples, considering all the four possible connections, namely, 1-shock/2-shock, 1-rarefaction/2-rarefaction, 1-rarefaction/2-shock and 1-shock/2-rarefaction are presented.     

Mathematical modeling of magnetorheological fluids

Magnetorheological (MR) fluids are a class of smart materials whose rheological properties may be rapidly modified by the application of a magnetic field. These materials typically consist of micron-sized ferrous particles dispersed in a fluid. In the present paper, we consider the full system of equations as well as the Clausius-Duhem inequality for moving isotropic MR fluids in an electro-magnetic field. We present the material constitutive relations for a non-Newtonian incompressible MR fluid. To illustrate the validity of the constitutive relations, the flow of a MR fluid between two parallel fixed plates under the influence of a constant magnetic field perpendicular to the flow direction is considered.Received: 14 July 2003, Accepted: 18 May 2004, Published online: 22 February 2005 Correspondence to: A. Dorfmann     

一种计入横向剪切变形的各向异性板弯曲分析新方法

本文提出一种各向异性板弯曲分析新方法,既考虑了横向剪切变形影响,又具有经典理论的简明性,易被工程人员接受。算例表明,本文方法分析精度较高。     

Linear and Nonlinear Aerodynamic Theory of Interaction between Flexible Long Structure and Wind

ＩｎｔｒｏｄｕｃｔｉｏｎＡｌｏｎｇｗｉｔｈｔｈｅｉｎｃｒｅａｓｅｏｆｔｈｅｆｌｅｘｉｂｉｌｉｔｙｏｆｓｌｅｎｄｅｒｓｔｒｕｃｔｕｒｅｓｅｔｉｎｔｈｅｗｉｎｄｆｉｅｌｄ ,ｓｕｃｈａｓｔｈｅｌｏｎｇ＿ｓｐａｎｂｒｉｄｇｅａｎｄｈｉｇｈ＿ｒｉｓｅｂｕｉｌｄｉｎｇ ,ｔｈｅａｃｔｉｏｎｓｏｆｕｎｓｔｅａｄｙａｅｒｏｄｙｎａｍｉｃｆｏｒｃｅｓｗｉｌｌａｃ ｃｅｌｅｒａｔｅ,ａｎｄｈｏｗｔｏｓｏｌｖｅｔｈｅａｅｒｏｅｌａｓｔｉｃｐｒｏｂｌｅｍｓｉｓｇｅｔｔｉｎｇｍｏｒｅｉｍｐｏｒｔａｎｔ.Ｔｈｅａｅｒｏｅｌａｓ…     

Mathematical modeling of the motion of a portion of helium under pulsed injection over a fixed bed of cenospheres

A mathematical model is constructed and an analytical solution is obtained for the problem of a one-dimensional steady flow of a mixture of different gases with hollow permeable particles. The case of a one-dimensional unsteady flow of such a mixture is analyzed numerically. The numerical solutions are compared with experimental data on the motion of the peak concentration of helium in a fixed bed filled with cenospheres (solid hollow permeable spherical particles). The permeability of cenosphere walls and the drag coeficient of cenospheres in the gas flow are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 92–102, May–June, 2007.     

Non‐intrusive reduced‐order modeling for multiphase porous media flows using Smolyak sparse grids

In this article, we describe a non‐intrusive reduction method for porous media multiphase flows using Smolyak sparse grids. This is the first attempt at applying such an non‐intrusive reduced‐order modelling (NIROM) based on Smolyak sparse grids to porous media multiphase flows. The advantage of this NIROM for porous media multiphase flows resides in that its non‐intrusiveness, which means it does not require modifications to the source code of full model. Another novelty is that it uses Smolyak sparse grids to construct a set of hypersurfaces representing the reduced‐porous media multiphase problem. This NIROM is implemented under the framework of an unstructured mesh control volume finite element multiphase model. Numerical examples show that the NIROM accuracy relative to the high‐fidelity model is maintained, whilst the computational cost is reduced by several orders of magnitude. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamic and quasi-static bending of saturated poroelastic Timoshenko cantilever beam

Based on the three-dimensional Gurtin-type variational principle of the incompressible saturated porous media, a one-dimensional mathematical model for dynamics of the saturated poroelastic Timoshenko cantilever beam is established with two assumptions, i.e., the deformation satisfies the classical single phase Timoshenko beam and the movement of the pore fluid is only in the axial direction of the saturated poroelastic beam. Under some special cases, this mathematical model can be degenerated into the Euler-Bernoulli model, the Rayleigh model, and the shear model of the saturated poroelastic beam, respectively. The dynamic and quasi-static behaviors of a saturated poroelastic Timoshenko cantilever beam with an impermeable fixed end and a permeable free end subjected to a step load at its free end are analyzed by the Laplace transform. The variations of the deflections at the beam free end against time are shown in figures. The influences of the interaction coefficient between the pore fluid and the solid skeleton as well as the slenderness ratio of the beam on the dynamic/quasi-static performances of the beam are examined. It is shown that the quasi-static deflections of the saturated poroelastic beam possess a creep behavior similar to that of viscoelastic beams. In dynamic responses, with the increase of the slenderness ratio, the vibration periods and amplitudes of the deflections at the free end increase, and the time needed for deflections approaching to their stationary values also increases. Moreover, with the increase of the interaction coefficient, the vibrations of the beam deflections decay more strongly, and, eventually, the deflections of the saturated poroelastic beam converge to the static deflections of the classic single phase Timoshenko beam.     

Contact Interaction Between a Laminated Shell of Revolution and a Rigid or Elastic Foundation

A method is proposed to solve the contact problem for laminated anisotropic shells of revolution. The method is based on a two-dimensional model that accounts for transverse shears and reduction. Also the method is based on the method of successive approximations, the generalized pseudo-force method, and a numerical-analytical method of solving boundary-value problems. The results obtained for a cylindrical shell of complex thickness structure are compared with those obtained in three-dimensional formulation__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 5, pp. 68–75, May 2005.     

Laboratory observation of nonlocal dispersion

This work presents the results of a one-dimensional experimental investigation of contaminant transport in heterogeneous porous media. The usual transport equations fail to adequately predict dispersion in such systems, and new theories to account for the distinctions have not yet been examined experimentally. We use a one-dimensional porous media which is heterogeneous on the scale of observation to determine if the phenomena predicted by the new theories are observable.The experimental media are constructed from distinct layers of spherical glass beads packed into cylindrical columns of Lucite. Flow was in the direction perpendicular to the layers. Dispersion was measured by recording the concentration of a chloride tracer as a function of time and position. The scale of measurement was finer than the scale of the heterogeneity. The results show that the mixing between miscible fluids was affected by transitions in the system parameters, before the transitions were encountered by the mixing zone. This newly observed phenomenon has been interpreted as a nonlocal effect, and it begins to verify the new predictive theories.     

轴向扩散不可压饱和多孔弹性梁的拟静态弯曲

建立了横观各向同性不可压饱和多孔弹性梁拟静态弯曲的数学模型,并给出了一般的求解方法.作为例子,研究了端部不同渗透条件对梁中点承受突加常集中载荷作用的饱和多孔悬臂梁拟静态弯曲的影响,给出了挠度和孔隙流体压力等效力偶沿梁轴线的分布以及随时间的响应曲线.结果表面:端部渗透条件对饱和多孔弹性梁的弯曲行为有显著的影响,梁的弯曲挠度既可随时间单调递增、亦可单调递减,其性态依赖于梁端部的渗透条件.同时发现不同于经典单相弹性梁,由于孔隙流体压力的作用,不承受载荷作用的梁段亦发生弯曲,并且Mandel-Cryer效应亦存在于不可压饱和多孔弹性梁的拟静态响应中,这些结果有助于揭示传热管道、植物根茎等力学行为的机理.