Continuity Of Normal Stress Research Articles

Bending and free vibration analysis of symmetric and unsymmetric functionally graded CNT reinforced sandwich beams containing softcore

In the present work, bending and free vibration analyses of functionally graded carbon nanotube-reinforced (FG-CNTR) sandwich beams are carried out using finite element-based higher-order zigzag theory. Face sheets are assumed to be made up of FG-CNTR composite, and the core is assumed to be made up of balsa wood (softcore). The present formulation also takes into account transverse normal stresses. The computational model incorporates transverse shear stress and transverse normal stress continuity condition at interfaces. Zero transverse shear stress condition at the bottom and top surfaces of the beam is also satisfied. The principle of minimum potential energy is employed for carrying out bending analysis, while Hamilton’s principle is adopted for free vibration analysis. The investigation is carried out for different gradation laws which govern the distribution of CNTs across the thickness of face sheets. The influence of the core’s thickness on stresses and displacements is also critically analyzed in the present work. It has been observed that the thickness of the core and CNT gradation law significantly affect the mechanical behavior of the sandwich FG-CNTRC beam.

Domain decomposition and partitioning methods for mixed finite element discretizations of the Biot system of poroelasticity

We develop non-overlapping domain decomposition methods for the Biot system of poroelasticity in a mixed form. The solid deformation is modeled with a mixed three-field formulation with weak stress symmetry. The fluid flow is modeled with a mixed Darcy formulation. We introduce displacement and pressure Lagrange multipliers on the subdomain interfaces to impose weakly continuity of normal stress and normal velocity, respectively. The global problem is reduced to an interface problem for the Lagrange multipliers, which is solved by a Krylov space iterative method. We study both monolithic and split methods. In the monolithic method, a coupled displacement-pressure interface problem is solved, with each iteration requiring the solution of local Biot problems. We show that the resulting interface operator is positive definite and analyze the convergence of the iteration. We further study drained split and fixed stress Biot splittings, in which case we solve separate interface problems requiring elasticity and Darcy solves. We analyze the stability of the split formulations. Numerical experiments are presented to illustrate the convergence of the domain decomposition methods and compare their accuracy and efficiency.

Axisymmetric Slow Motion of a Porous Spherical Particle in a Viscous Fluid Using Time Fractional Navier–Stokes Equation

In this paper, we present the unsteady translational motion of a porous spherical particle in an incompressible viscous fluid. In this case, the modified Navier–Stokes equation with fractional order time derivative is used for conservation of momentum external to the particle whereas modified Brinkman equation with fractional order time derivative is used internal to the particle to govern the fluid flow. Stress jump condition for the tangential stress along with continuity of normal stress and continuity of velocity vectors is used at the porous–liquid interface. The integral Laplace transform technique is employed to solve the governing equations in fluid and porous regions. Numerical inversion code in MATLAB is used to obtain the solution of the problem in the physical domain. Drag force experienced by the particle is obtained. The numerical results have been discussed with the aid of graphs for some specific flows, namely damping oscillation, sine oscillation and sudden motion. Our result shows a significant contribution of the jump coefficient and the fractional order parameter to the drag force.

A hybrid mimetic spectral element method for three-dimensional linear elasticity problems

We introduce a domain decomposition structure-preserving method based on a hybrid mimetic spectral element method for three-dimensional linear elasticity problems in curvilinear conforming structured meshes. The method is an equilibrium method which satisfies pointwise equilibrium of forces. The domain decomposition is established through hybridization which first allows for an inter-element normal stress discontinuity and then enforces the normal stress continuity using a Lagrange multiplier which turns out to be the displacement in the trace space. Dual basis functions are employed to simplify the discretization and to obtain a higher sparsity. Numerical tests supporting the method are presented.

Effect of magnetic field on the slow motion of a porous spheroid: Brinkman’s model

The major goal of this work is to analyze the magnetic effect on the creeping viscous flow past a porous spheroidal particle, a particle of slightly deformed spherical shape. Brinkman’s model is proposed to govern the flow in the porous media. Boundary value problem considers the conditions of continuity of velocity components, continuity of normal stresses, and stress jump boundary condition for tangential stress. A transverse magnetic field of uniform nature is applied to the flow. An expression for the drag force acting on the spheroidal particle is derived analytically. The effects of the physical parameters involved in the flow like permeability, deformation, Hartmann number’s, viscosity ratio, and stress jump coefficient parameters are visualized through graphs and tables. The applied magnetic field seems to suppress the flow of fluid that leads to the increase in the drag experienced on the porous spheroid. It is also observed that the increase in the deformation, stress jump, and permeability decreases the drag coefficient. Our results without magnetic effect match with the results reported earlier in the literature.

Flow past composite cylindrical shell of porous layer with a liquid core: magnetic effect

The creeping flow of a magnetic fluid perpendicular to a porous cylindrical shell is investigated, employing the unit cell model. The viscous fluid is assumed to be flowing in three zones divided as fluid, annular porous, and cavity regions, respectively. We apply a uniform magnetic field in a transverse direction of flow and then emphasize the influence of the Hartmann layers which are developed in their vicinity. Modified Stokes and modified Brinkman’s equations are employed in the liquid and porous regions, respectively. Happel and Kuwabara cell models are used as the interface conditions at the cell surface, and at the fluid–porous interface, continuity of velocity components, continuity of normal stresses, and stress jump condition for tangential stresses are applied. An expression for Kozeny constant for the cylindrical shell is presented. Representation of Kozeny constant under the influences of the pertinent parameters such as Hartmann numbers, stress jump coefficient, fractional void volume, viscosity ratios, permeability, and separation is displayed through graphs and a table. The results are compared with the cases which do not involve magnetic effect. They reveal the strong impact of Hartmann’s numbers on the resisting force experienced on the cylindrical shell. The results agree well with previous available works.

INFLUENCE OF MAGNETIC FIELD ON THE STOKES FLOW THROUGH POROUS SPHEROID: HYDRODYNAMIC PERMEABILITY OF A MEMBRANE USING CELL MODEL TECHNIQUE

The present work concerns an analysis of the creeping flow of steady, axisymmetric Stokes flow of an electrically conducting, viscous, incompressible fluid through a swarm of porous spheroidal particles in the presence of a uniform magnetic field. Cell model technique has been used to model the problem. Four known boundary conditions, Happel's, Kuwabara's, Kvashnin's, and Cunningham's (Mehta-Morse's), are used to find the hydrodynamic permeability of the membrane built up by porous spheroidal particles by using the perturbation method. The stress jump boundary condition for tangential stresses, along with the continuity of normal stress and velocity components, is used at the fluid-porous interface. The variation of the dimensionless hydrodynamic permeability of the membrane with the stress jump coefficient, the Hartmann number, and the dimensionless permeability of the porous region and particle volume fraction are discussed. The presented model can be used for the evaluation of changing hydrodynamic permeability of a membrane under the influence of a uniform magnetic field in pressure-driven processes (nano-, reverse osmosis, and microfiltration).

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.

Viscous flow past a porous sphere within a nonconcentric fictitious spherical cell

The quasi steady axisymmetrical flow of an incompressible viscous fluid past an assemblage of porous sphere situated at an arbitrary position within a virtual spherical cell along the line connecting their centers is analyzed using a combined analytical–numerical technique. At the fluid–porous interface, the stress jump boundary condition for the tangential stresses along with continuity of normal stress and velocity components are employed. The Brinkman model governs the flow inside the porous particle and the flow in the fictitious envelope medium is governed by Stokes equations. A general solution is constructed from the superposition of the fundamental solutions in the two spherical coordinate systems based on both the porous particle and virtual spherical cell. Boundary conditions on the particle’s surface and the fictitious spherical envelope for four known cell models are satisfied by a collocation method for truncated series. Numerical solutions for the drag force exerted on the porous sphere in the presence of the cell are obtained for various cases of the effective distance between the center of the porous particle and the fictitious envelope, the volume ratio of the porous particle over the surrounding sphere, the viscosity ratio, the stress jump coefficient, and a coefficient that is proportional to the permeability. Streamlines in and around porous sphere are presented for the unit cell models at different values of relevant physical parameters. In the limits of the motions of the porous particle in the concentric position with the cell surface and near the cell surface with a small curvature, the numerical values of the normalized drag force are in good agreement with the available values in the literature.

Modified equations of finite-size layered plates made of orthotropic material. Comparison of the results of numerical calculations with analytical solutions

This paper describes the modified bending equations of layered orthotropic plates in the first approximation. The approximation of the solution of the equation of the three-dimensional theory of elasticity by the Legendre polynomial segments is used to obtain differential equations of the elastic layer. For the approximation of equilibrium equations and boundary conditions of three-dimensional theory of elasticity, several approximations of each desired function (stresses and displacements) are used. The stresses at the internal points of the plate are determined from the defining equations for the orthotropic material, averaged with respect to the plate thickness. The construction of the bending equations of layered plates for each layer is carried out with the help of the elastic layer equations and the conjugation conditions on the boundaries between layers, which are conditions for the continuity of normal stresses and displacements. The numerical solution of the problem of bending of the rectangular layered plate obtained with the help of modified equations is compared with an analytical solution. It is determined that the maximum error in determining the stresses does not exceed 3 %.

Modeling elastic wave propagation in fluid-filled boreholes drilled in nonhomogeneous media: BEM-MLPG versus BEM-FEM coupling

The efficiency of two coupling formulations, the boundary element method (BEM)-meshless local Petrov–Galerkin (MLPG) versus the BEM-finite element method (FEM), used to simulate the elastic wave propagation in fluid-filled boreholes generated by a blast load, is compared. The longitudinal geometry is assumed to be invariant in the axial direction (2.5D formulation). The material properties in the vicinity of the borehole are assumed to be nonhomogeneous as a result of the construction process and the ageing of the material. In both models, the BEM is used to tackle the propagation within the fluid domain inside the borehole and the unbounded homogeneous domain. The MLPG and the FEM are used to simulate the confined, damaged, nonhomogeneous, surrounding borehole, thus utilizing the advantages of these methods in modeling nonhomogeneous bounded media. In both numerical techniques the coupling is accomplished directly at the nodal points located at the common interfaces. Continuity of stresses and displacements is imposed at the solid–solid interface, while continuity of normal stresses and displacements and null shear stress are prescribed at the fluid–solid interface. The performance of each coupled BEM-MLPG and BEM-FEM approach is determined using referenced results provided by an analytical solution developed for a circular multi-layered subdomain. The comparison of the coupled techniques is evaluated for different excitation frequencies, axial wavenumbers and degrees of freedom (nodal points).

Effect of magnetic field on the hydrodynamic permeability of a membrane built up by porous spherical particles

In this paper, an analysis of steady, axi-symmetric Stokes flow of an electrically conducting viscous incompressible fluid through spherical particle covered by porous shell in presence of uniform magnetic field is presented. To model flow through the swarm of spherical particles, cell model technique has been used, i.e. porous spherical shell is assumed to be confined within a hypothetical cell of the same geometry. At the fluid-porous interface, the stress jump boundary condition for tangential stresses along with continuity of normal stress and velocity components are used. Four known boundary conditions on the hypothetical surface were considered and compared: Happel’s, Kuwabara’s, Kvashnin’s and Cunningham’s (Mehta−Morse’s) condition. The effect of stress jump coefficient, Hartmann number, and dimensionless permeability of the porous region as well as particle volume fraction on the hydrodynamic permeability and streamlines were discussed. The patterns of streamlines were also obtained.

Multiscale coupling of compliant and rigid walls blood flow models

SummaryNumerical methods based on geometrical multiscale models of blood flows solve for averaged flow statistics on a network of vessels while providing more detailed information about fluid dynamics in a specific region of interest. In such an approach, a 3D model based on the Navier–Stokes equations posed in a domain with rigid walls is often used to describe blood flow dynamics in the refined region. While ignoring elasticity effects in 3D models is plausible in certain applications and saves computational time significantly, coupling such models with 1D flow models may result in non‐physiological phenomena in the computed solutions. Thus, the immediate coupling conditions based on continuity of normal stresses, flow rate, pressure, or a combination of thereof do not account for the inconsistency between elasticity effects in the 1D model and the non‐compliance of the 3D model. In this paper, we introduce and study an auxiliary absorbing 0D model, which is placed at the interface between 1D and 3D models. A virtual device mimics the effect of the 3D model compliance and hence reduces pressure wave reflection and instabilities caused by the inconsistency. The absorbing model is developed from basic mechanical principles. As a result, parameters of the 0D model can be designed based on hemodynamic data. We analyze the stability of the geometrical multiscale model and perform several numerical experiments to assess its computational efficiency. Copyright © 2016 John Wiley & Sons, Ltd.

Fluid-structure interaction between an incompressible, viscous 3D fluid and an elastic shell with nonlinear Koiter membrane energy

The long-time existence of a weak solution is proved for a nonlinear, fluid-structure interaction (FSI) problem between an incompressible, viscous fluid and a semilinear cylindrical Koiter membrane shell with inertia. No axial symmetry is assumed in the problem. The fluid flow is driven by the time dependent dynamic pressure data prescribed at the inlet and outlet boundaries of the 3D cylindrical fluid domain. The fluid and the elastic structure are fully coupled via continuity of velocity and continuity of normal stresses. Global existence of a weak solution is proved as long as the lateral walls of the cylinder do not touch each other. The main novelty of the work is the nonlinearity in the structure model: the model accounts for the fully nonlinear Koiter membrane energy, supplemented with a small linear fourth-order derivative term modeling the bending rigidity of shells. The existence proof is constructive, and it is based on an operator splitting scheme. A version of this scheme can be implemented for the numerical simulation of the underlying FSI problem by extending the FSI solver, developed by the authors in [5], to include the nonlinearity in the structure model discussed in this manuscript.

Hydrodynamic permeability of a membrane composed of porous spherical particles in the presence of uniform magnetic field

This work concerns the flow of an incompressible viscous fluid past a porous sphere in presence of transverse applied uniform magnetic field, using particle-in-cell method. The Brinkman equations are used in porous region and the Stokes equations for non-porous region. At the fluid-porous interface, the stress jump boundary condition for tangential stresses along with continuity of normal stress and velocity components are used. Four known boundary conditions on the hypothetical surface are considered and compared: Happel’s, Kuwabara’s, Kvashnin’s and Cunningham’s (Mehta-Morse’s condition). The hydrodynamic drag force experienced by a porous spherical particle in a cell and hydrodynamic permeability of membrane built up by porous spherical particles are evaluated. The patterns of streamlines are also obtained and discussed. The effect of stress jump coefficient, Hartmann number, dimensionless specific permeability of the porous particle and particle volume fraction on the hydrodynamic permeability and streamlines are discussed. Some previous results for hydrodynamic drag force and dimensionless hydrodynamic permeability have been verified.

Slow Motion of a Porous Sphere Translating Along the Axis of a Circular Cylindrical Pore Subject to a Stress Jump Condition

The coupled flow problem of an incompressible axisymmetrical quasisteady motion of a porous sphere translating in a viscous fluid along the axis of a circular cylindrical pore is discussed using a combined analytical–numerical technique. At the fluid–porous interface, the stress jump boundary condition for the tangential stress along with continuity of normal stress and velocity components are employed. The flow through the porous particle is governed by the Brinkman model and the flow in the outside porous region is governed by Stokes equations. A general solution for the field equations in the clear region is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are satisfied first at the cylindrical pore wall by the Fourier transforms and then on the surface of the porous particle by a collocation method. The collocation solutions for the normalized hydrodynamic drag force exerted by the clear fluid on the porous particle is calculated with good convergence for various values of the ratio of radii of the porous sphere and pore, the stress jump coefficient, and a coefficient that is proportional to the permeability. The shape effect of the cylindrical pore on the axial translation of the porous sphere is compared with that of the particle in a spherical cavity; it found that the porous particle in a circular cylindrical pore in general attains a lower hydrodynamic drag than in a spherical envelope.

Stokes flow past an assemblage of axisymmetric porous spherical shell-in-cell models: effect of stress jump condition

The quasisteady axisymmetrical flow of an incompressible viscous fluid past an assemblage of porous concentric spherical shell-in-cell model is studied. Boundary conditions on the cell surface that correspond to the Happel, Kuwabara, Kvashnin and Cunningham/Mehta-Morse models are considered. At the fluid-porous interfaces, the stress jump boundary condition for the tangential stresses along with continuity of normal stress and velocity components are employed. The Brinkman’s equation in the porous region and the Stokes equation for clear fluid are used. The hydrodynamic drag force acting on the porous shell by the external fluid in each of the four boundary conditions on the cell surface is evaluated. It is found that the normalized mobility of the particles (the hydrodynamic interaction among the porous shell particles) depends not only on the permeability of the porous shells and volume fraction of the porous shell particles, but also on the stress jump coefficient. As a limiting case, the drag force or mobility for a suspension of porous spherical shells reduces to those for suspensions of impermeable solid spheres and of porous spheres with jump.

ON THE INTERFACE LAW BETWEEN A DEFORMABLE POROUS MEDIUM CONTAINING A VISCOUS FLUID AND AN ELASTIC BODY

In this paper we undertake a rigorous derivation of the interface conditions between a poroelastic medium (the pay zone) and an elastic body (the non-pay zone). We assume that the poroelastic medium contains a pore structure of the characteristic size ε and that the fluid/structure interaction regime corresponds to diphasic Biot's law. The question is challenging because the Biot's equations for the poroelastic part contain one unknown more than the Navier equations for the non-pay zone. The solid part of the pay zone (the matrix) is elastic and the pores contain a viscous fluid. The fluid is assumed viscous and slightly compressible. We study the case when the contrast of property is of order ε2, i.e. the normal stress of the elastic matrix is of the same order as the fluid pressure. We assume a periodic matrix and obtain the a priori estimates. Then we let the characteristic size of the inhomogeneities tend to zero and pass to the limit in the sense of the two-scale convergence. The obtained effective equations represent a two-scale system for three displacements and two pressures, coupled through the interface conditions with the Navier equations for the elastic displacement in the non-pay zone. We prove that the appropriate interface conditions at the interface between an elastic and a poroelastic medium are: (i) the effective displacement continuity, (ii) the effective normal stress continuity and (iii) the normal Darcy velocity zero from the poroelastic side. In addition we determine the effective boundary conditions for the contact between a poroelastic body and a rigid obstacle, giving us the effective outer boundary conditions.

On the Slow Viscous Flow through a Swarm of Solid Spherical Particles Covered by Porous Shells

This paper concerns the slow viscous flow through a swarm of concentric clusters of porous spherical particles. An aggregate of clusters of porous spherical particles is considered as a hydro-dynamically equivalent to a porous spherical shell enclosing a solid spherical core. The Brinkman equation inside and the Stokes equation outside the porous spherical shell in their stream function formulations are used. As boundary conditions, continuity of velocity, continuity of normal stress and stress-jump condition at the porous and fluid interface, the continuity of velocity components on the solid spherical core are employed. On the hypothetical surface, uniform velocity and Happel boundary conditions are used. The drag force experienced by each porous spherical shell in a cell is evaluated. As a particular case, the drag force experienced by a porous sphere in a cell with jump is also investigated. The earlier results reported for the drag force by Davis and Stone[5] for the drag force experienced by a porous sphere in a cell without jump, Happel[2] for a solid sphere in a cell and Qin and Kaloni[4] for a porous sphere in an unbounded medium have been then deduced. Representative results are presented in graphical form and discussed.

Mass Transport in a Microchannel Bioreactor with a Porous Wall

A two-dimensional flow model, incorporating mass transport, has been developed to simulate flow in a microchannel bioreactor with a porous wall. A two-domain method was implemented which was based on finite volume method. For the porous-fluid interface, a stress jump condition was used with continuity of normal stress; and the mass interfacial conditions were continuities of mass and mass flux. Two parameters are defined to characterize the mass transports in the fluid and porous regions. The porous Damkohler number is the ratio of consumption to diffusion of the substrates in the porous medium. The fluid Damkohler number is the ratio of substrate consumption in the porous medium to substrate convection in the fluid region. The concentration results are found to be well correlated by the use of a reaction-convection distance parameter which incorporates the effects of axial distance, substrate consumption and convection.