Law of nonlinear flow in saturated clays and radial consolidation

It was derived that micro-scale amount level of average pore radius of clay changed from 0.01 to 0.1 micron by an equivalent concept of flow in porous media.There is good agreement between the derived results and test ones.Results of experiments show that flow in micro-scale pore of saturated clays follows law of nonlinear flow.Theoretical analyses demonstrate that an interaction of solid-liquid interfaces varies inversely with permeability or porous radius.The interaction is an important reason why nonlinear flow in saturated clays occurs.An exact mathematical model was presented for nonlinear flow in micro-scaie pore of saturated clays.Dimension and physical meanings of parameters of it are definite.A new law of nonlinear flow in saturated clays was established.It can describe characteristics of flow curve of the whole process of the nonlinear flow from low hydraulic gradient to high one.Darcy law is a special case of the new law.A math- ematical model was presented for consolidation of nonlinear flow in radius direction in saturated clays with constant rate based on the new law of nonlinear flow.Equations of average mass conservation and moving boundary,and formula of excess pore pressure distribution and average degree of consolidation for nonlinear flow in saturated clay were derived by using an idea of viscous boundary layer,a method of steady state in stead of transient state and a method of integral of an equation.Laws of excess pore pressure distribution and changes of average degree of consolidation with time were obtained.Re- suits show that velocity of moving boundary decreases because of the nonlinear flow in saturated clay.The results can provide geology engineering and geotechnical engineering of saturated clay with new scientific bases.Calculations of average degree of consolidation of the Darcy flow are a special case of that of the nonlinear flow.     

Peristaltic flow of a heated Jeffrey fluid inside an elliptic duct:streamline analysis

The peristaltic flow of a heated Jeffrey fluid inside a duct with an elliptic cross-section is studied. A thorough heat transfer mechanism is interpreted by analyzing the viscous effects in the energy equation. The governing mathematical equations give dimensionless partial differential equations after simplification. The final simplified form of the mathematical equations is evaluated with respect to the relevant boundary conditions, and the exact solution is attained. The results are further i...     

Nonlocal and strain gradient effects on nonlinear forced vibration of axially moving nanobeams under internal resonance conditions

Based on the nonlocal strain gradient theory(NSGT), a model is proposed for an axially moving nanobeam with two kinds of scale effects. The internal resonanceaccompanied fundamental harmonic response of the external excitation frequency in the vicinities of the first and second natural frequencies is studied by adopting the multivariate Lindstedt-Poincaré(L-P) method. Based on the root discriminant of the frequencyamplitude equation under internal resonance conditions, theoretical analyses are p...     

Heat transfer and Helmholtz-Smoluchowski velocity in Bingham fluid flow

A mathematical study is developed for the electro-osmotic flow of a nonNewtonian fluid in a wavy microchannel in which a Bingham viscoplastic fluid model is considered. For electric potential distributions, a Poisson-Boltzmann equation is employed in the presence of an electrical double layer(EDL). The analytical solutions of dimensionless boundary value problems are obtained with the Debye-Huckel theory, the lubrication theory, and the long wavelength approximations. The effects of the Debyelen...     

Regular perturbation solution of Couette flow(non-Newtonian)between two parallel porous plates: a numerical analysis with irreversibility

The unavailability of wasted energy due to the irreversibility in the process is called the entropy generation.An irreversible process is a process in which the entropy of the system is increased.The second law of thermodynamics is used to define whether the given system is reversible or irreversible.Here,our focus is how to reduce the entropy of the system and maximize the capability of the system.There are many methods for maximizing the capacity of heat transport.The constant pressure gradient or motion of the wall can be used to increase the heat transfer rate and minimize the entropy.The objective of this study is to analyze the heat and mass transfer of an Eyring-Powell fluid in a porous channel.For this,we choose two different fluid models,namely,the plane and generalized Couette flows.The flow is generated in the channel due to a pressure gradient or with the moving of the upper lid.The present analysis shows the effects of the fluid parameters on the velocity,the temperature,the entropy generation,and the Bejan number.The nonlinear boundary value problem of the flow problem is solved with the help of the regular perturbation method.To validate the perturbation solution,a numerical solution is also obtained with the help of the built-in command NDSolve of MATHEMATICA 11.0.The velocity profile shows the shear thickening behavior via first-order Eyring-Powell parameters.It is also observed that the profile of the Bejan number has a decreasing trend against the Brinkman number.Whenηi→0(i=1,2,3),the Eyring-Powell fluid is transformed into a Newtonian fluid.     

SLIPPAGE SOLUTION OF GAS PRESSURE DISTRIBUTION IN PROCESS OF LANDFILL GAS SEEPAGE

A mathematical model of landfill gas migration was established under presumption of the effect of gas slippage. The slippage solutions to the nonlinear mathematical model were accomplished by the perturbation and integral transformation method. The distribution law of gas pressure in landfill site was presented under the conditions of considering and neglecting slippage effect. Sensitivity of the model input parameters was analyzed. The model solutions were compared to observation values. Results show that gas slippage effect has a large impact on gas pressure distribution. Landfill gas pressure and pressure gradient considering slippage effect is lower than that neglecting slippage effect, with reasonable agreement between model solution and measured data. It makes clear that the difference between considering and neglecting slippage effect is obvious and the effects of coupling cannot be ignored. The theoretical basis is provided for engineering design of security control and decision making of gas exploitation in landfill site. The solutions give scientific foundation to analyzing well test data in the process of low-permeability oil gas reservoir exploitation.     

Exact solutions for nonlinear transient flow model including a quadratic gradient term

The models of the nonlinear radial flow for the infinite and finite reservoirs including a quadratic gradient term were presented. The exact solution was given in real space for flow equation including quadratic gradiet term for both constant-rate and constant pressure production cases in an infinite system by using generalized Weber transform.Analytical solutions for flow equation including quadratic gradient term were also obtained by using the Hankel transform for a finite circular reservoir case. Both closed and constant pressure outer boundary conditions are considered. Moreover, both constant rate and constant pressure inner boundary conditions are considered. The difference between the nonlinear pressure solution and linear pressure solution is analyzed. The difference may be reached about 8% in the long time. The effect of the quadratic gradient term in the large time well test is considered.     

Characteristic finite difference method and application for moving boundary value problem of coupled system

The coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.A kind of characteristic finite difference schemes is put forward,from which optimal order estimates in l^2 norm are derived for the error in the approximate solutions.The research is important both theoretically and practically for the model analysis in the field,the model numerical method and software development.     

DYNAMIC ANALYSIS OF A GRADIENT ELASTIC POLYMERIC FIBER

A dynamic analysis of an elastic gradient-dependent polymeric fiber subjected to a periodic excitation is considered.A nonlinear gradient elasticity constitutive equation with strain-dependent gradient coefficients is first derived and the dispersive wave propagation properties for its linearized counterpart are briefly discussed.For the linearized problem a variational formulation is also developed to obtain related boundary conditions of both classical(standard)and non-classical(gradient)type.Analytical solutions in the form of Fourier series for the fiber’s displacement and strain fields are provided.The solutions depend on a dimensionless scale parameter(the diameter to length radio d = D/L)and,therefore,size effects are captured.     

An improved correlation of the pressure drop in stenotic vessels using Lorentz’s reciprocal theorem

A mathematical model of the human cardiovascular system in conjunction with an accurate lumped model for a stenosis can provide better insights into the pressure wave propagation at pathological conditions. In this study, a theoretical relation between pressure drop and flow rate based on Lorentz’s reciprocal theorem is derived, which offers an identity to describe the relevance of the geometry and the convective momentum transport to the drag force. A voxelbased simulator V-FLOW VOF3 D, where the vessel geometry is expressed by using volume of fluid(VOF) functions, is employed to find the flow distribution in an idealized stenosis vessel and the identity was validated numerically. It is revealed from the correlation that the pressure drop of NS flow in a stenosis vessel can be decomposed into a linear term caused by Stokes flow with the same boundary conditions, and two nonlinear terms. Furthermore, the linear term for the pressure drop of Stokes flow can be summarized as a correlation by using a modified equation of lubrication theory, which gives favorable results compared to the numerical ones. The contribution of the nonlinear terms to the pressure drop was analyzed numerically, and it is found that geometric shape and momentum transport are the primary factors for the enhancement of drag force. This work paves a way to simulate the blood flow and pressure propagation under different stenosis conditions by using 1D mathematical model.     

A meshfree-based local Galerkin method with condensation of degree of freedom for elastic dynamic analysis简

Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark＇s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.     

Three dimensional simulations of penetrative convection in a porous medium with internal heat sources

The problem of penetrative convection in a fluid saturated porous medium heated internally is analysed. The linear instability theory and nonlinear energy theory are derived and then tested using three dimensions simulation.Critical Rayleigh numbers are obtained numerically for the case of a uniform heat source in a layer with two fixed surfaces. The validity of both the linear instability and global nonlinear energy stability thresholds are tested using a three dimensional simulation. Our results show that the linear threshold accurately predicts the onset of instability in the basic steady state. However, the required time to arrive at the basic steady state increases significantly as the Rayleigh number tends to the linear threshold.     

An enhanced treatment of boundary conditions in implicit smoothed particle hydrodynamics简

In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme （Implicit SPH）. Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares （MLS） interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.     

Effect of electric boundary conditions on crack propagation in ferroelectric ceramics

In this paper, the effect of electric boundary conditions on Mode I crack propagation in ferroelectric ceramics is studied by using both linear and nonlinear piezoelectric fracture mechanics. In linear analysis, impermeable cracks under open circuit and short circuit are analyzed using the Stroh formalism and a rescaling method. It is shown that the energy release rate in short circuit is larger than that in open circuit. In nonlinear analysis, permeable crack conditions are used and the nonlinear effect of domain switching near a crack tip is considered using an energy-based switching criterion proposed by Hwang et al.（Acta Metal. Mater.,1995）. In open circuit, a large depolarization field induced by domain switching makes switching much more diffcult than that in short circuit. Analysis shows that the energy release rate in short circuit is still larger than that in open circuit, and is also larger than the linear result. Consequently,whether using linear or nonlinear fracture analysis, a crack is found easier to propagate in short circuit than in open circuit, which is consistent with the experimental observations of Kounga Njiwa et al.（Eng. Fract. Mech., 2006）.     

Effect of quadratic pressure gradient term on a one-dimensional moving boundary problem based on modified Darcy’s law

A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equations with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be stably numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy’s flow problem, the exact analytical solution for the case of Darcy’s flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclusion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equations; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pressure gradient increasing for the one-dimensional problem.     

Interval finite difference method for steady-state temperature field prediction with interval parameters

A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.     

A cellular scale numerical study of the effect of mechanical properties of erythrocytes on the near-wall motion of platelets

The effect of mechanical properties of erythrocytes on the near-wall motion of platelets was numerically studied with the immersed boundary method. Cells were modeled as viscous-fluid-filled capsules surrounded by hyper-elastic membranes with negligible thickness. The numerical results show that with the increase of hematocrit, the near-wall approaching of platelets is enhanced, with which platelets exhibit larger deformation and orientation angle of its near-wall tank-treading motion, and the lateral force pushing platelets to the wall is increased with larger fluctuation amplitude. Meanwhile the near-wall approaching is reduced by increasing the stiffness of erythrocytes.     

The application of three-dimensional hydroelastic analysis of ship structures in Pekeris hydro-acoustic waveguide environment简

The hydroelastic analysis and sonoelastic analysis methods are incorporated with the Green＇s function of the Pekeris ocean hydro-acoustic waveguide model to produce a three-dimensional sonoelastic analysis method for ships in the ocean hydro-acoustic environment. The seabed condition is represented by a penetrable boundary of prescribed density and sound speed. This method is employed in this paper to predict the vibration and acoustic radiation of a 1 500 t Small Water Area Twin Hull （SWATH） ship in shallow sea acoustic environment. The wet resonant frequencies and radiation sound source levels are predicted and compared with the measured results of the ship in trial.     

Stability analysis of a thin-walled cylinder in turning operation using the semi-discretization method

In this work, the stability of a flexible thin cylindrical workpiece in turning is analyzed. A process model is derived based on a finite element representation of the workpiece flexibility and a nonlinear cutting force law. Repeated cutting of the same surface due to overlapping cuts is modeled with the help of a time delay. The stability of the so obtained system of periodic delay differential equations is then determined using an approximation as a time-discrete system and Floquet theory. The time-discrete system is obtained using the semi-discretization method. The method is implemented to analyze the stability of two different workpiece models of different thicknesses for different tool positions with respect to the jaw end. It is shown that the stability chart depends on the tool position as well as on the thickness.     

Frequency domain fundamental solutions for a poroelastic half-space

In frequency domain, the fundamental solutions for a poroelastic half-space are re-derived in the context of Biot＇s theory. Based on Biot＇s theory, the governing field equations for the dynamic poroelasicity are established in terms of solid displacement and pore pressure. A method of potentials in cylindrical coordinate system is proposed to decouple the homogeneous Biot＇s wave equations into four scalar Helmholtz equations, and the general solutions to these scalar wave equations are obtained. After that, spectral Green＇s functions for a poroelastic full-space are found through a decomposition of solid displacement, pore pressure, and body force fields. Mirror-image technique is then applied to construct the half-space fundamental solutions.Finally, transient responses of the half-space to buried point forces are examined.