Modeling and simulation of relaxation in viscoelastic open cell materials and structures

Open cell materials with cubic anisotropy and structures made thereof are investigated with respect to their linear viscoelastic properties, in particular their relaxation behavior. The study is concerned with the prediction of the effective behavior which results from the isotropic bulk material properties as well as the cellular architecture. Finite Element Method simulations of three-dimensional structures are employed to predict the effective response to a wide range of loading modes in the time domain.For predicting the properties of the cellular materials and structures by the Finite Element Method different modeling strategies are employed. The first approach is a periodic unit cell method modeling an infinite medium by means of periodic boundary conditions. This way the entire effective linear viscoelastic constitutive behavior can be computed. However, it is not possible to capture effects as being attributed to traction free faces and load introduction in specimens or structures. A second approach follows to account for these effects by generating finite sample models to represent situations which occur in experimental testing. Finally, an analytical constitutive material law is developed to model linear viscoelasticity for cubic anisotropy in the time domain. It is implemented into the commercial Finite Element software ABAQUS/Standard and the material parameters are gained from the unit cell investigations. This enables the simulation of structures, parts, and components which consist or contain such cellular materials.     

Solution of flow around complex‐shaped surfaces by an immersed boundary‐body conformal enrichment method

Solving the flow around objects with complex shapes may involve extensive meshing work that has to be repeated each time a change in the geometry is needed. Time consuming meshing can be avoided when the solution algorithm can tackle grids that do not fit the shape of immersed objects. This work presents applications of a recently proposed immersed boundary—body conformal enrichment method to the solution of the flow around complex shaped surfaces such as those of a metallic foam matrix. The method produces solutions of the flow satisfying accurately Dirichlet boundary conditions imposed on the immersed fluid/solid interface. The boundary of immersed objects is defined using a level‐set function, and the finite element discretization of interface elements is enriched with additional degrees of freedom, which are eliminated at element level. The method is first validated in the case of flow problems for which reference solutions on body‐conformal grids can be obtained: flow around an array of spheres and flow around periodic arrays of cylinders. Then, solutions are shown for the more complex flow inside a metallic foam matrix. A multiscale approach combining the solution at the pore level by the immersed boundary method and the macro‐scale solution with simulated permeability is used to solve actual experimental configurations. The computed pressure drop as a function of the flow rate on the macro scale configuration replicating two experimental setups is compared with the experimental data for various foam thicknesses. Copyright © 2011 National Research Council Canada     

Equivalent particle diameter and length scale for pressure drop in porous metals

The internal architecture of metal foam is significantly different from that of traditional porous media. This provides a set of challenges for understanding the fluid flow in this relatively new class of materials. This paper proposes that despite the geometrical differences between metal foam and traditional porous media, the Ergun correlation is a good fit for the linear pressure drop as a function of the Darcian velocity, provided that an appropriate equivalent particle diameter is used. The paper investigates an appropriate particle diameter considering the physics of energy dissipation, i.e. the viscous shear and the form drag. The above approach is supported by wind tunnel steady-state unidirectional pressure drop measurements for airflow through several isotropic open-cell aluminum foam samples having different porosities and pore densities. For each foam sample, the equivalent particle diameter correlated well with the surface area per unit volume of the foam. This was also very well valid for previous porous metal pressure drop data in the open literature.     

Modelling the flow of power law fluids in a packed bed using a volume-averaged equation of motion

The development of a theoretical model for the prediction of velocity and pressure drop for the flow of a viscous power law fluid through a bed packed with uniform spherical particles is presented. The model is developed by volume averaging the equation of motion. A porous microstructure model based on a cell model is used. Numerical solution of the resulting equation is effected using a penalty Galerkin finite element method. Experimental pressure drop values for dilute solutions of carboxymethylcellulose flowing in narrow tubes packed with uniformly sized spherical particles are compared to theoretical predictions over a range of operating conditions. Overall agreement between experimental and theoretical values is within 15%. The extra pressure drop due to the presence of the wall is incorporated directly into the model through the application of the no-slip boundary condition at the container wall. The extra pressure drop reaches a maximum of about 10% of the bed pressure drop without wall effect. The wall effect increases as the ratio of tube diameter to particle diameter decreases, as the Reynolds number decreases and as the power law index increases.     

Particle Flow Code Method-Based Meso-scale Identification for Seepage Failure of Soil Levee

The development and validation of a grid-based pore-scale numerical modelling methodology applied to five different commercial metal foam samples is described. The 3-D digital representation of the foam geometry was obtained by the use of X-ray microcomputer tomography scans, and macroscopic properties such as porosity, specific surface and pore size distribution are directly calculated from tomographic data. Pressure drop measurements were performed on all the samples under a wide range of flow velocities, with focus on the turbulent flow regime. Airflow pore-scale simulations were carried out solving the continuity and Navier–Stokes equations using a commercial finite volume code. The feasibility of using Reynolds-averaged Navier–Stokes models to account for the turbulence within the pore space was evaluated. Macroscopic transport quantities are calculated from the pore-scale simulations by averaging. Permeability and Forchheimer coefficient values are obtained from the pressure gradient data for both experiments and simulations and used for validation. Results have shown that viscous losses are practically negligible under the conditions investigated and pressure losses are dominated by inertial effects. Simulations performed on samples with varying thickness in the flow direction showed the pressure gradient to be affected by the sample thickness. However, as the thickness increased, the pressure gradient tended towards an asymptotic value.     

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围. 此外,引入一种混合式边界处理方法来模拟多孔介质惯性流, 结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点; 而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力-压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.      

Foam flow in vertical gas wells under liquid loading: Critical velocity and pressure drop prediction

Foam lift is one of the most cost effective methodologies for unloading gas wells. The surfactants are either injected intermittently or continuously to lift the liquid to the surface. By reducing the gravitational gradient and increasing the frictional gradient, the critical velocity at which liquid loading occurs is shifted to lower gas velocities. Currently, we do not have a methodology to predict the critical velocity (at the transition boundary of annular and intermittent flow) and the pressure drop under foam flow conditions.To address this, we measured several foam flow characteristics in both small scale and large scale facilities. Small scale facility involved measurement of foam carryover capacity as a function of time and surfactant concentration. Large scale facility involved measurement of liquid holdup, pressure drop, fraction of gas trapped in foam and foam holdup in 40-ft 2-in. and 4-in. tubing.We developed closure relationships for liquid hold up, foam holdup, fraction of gas trapped in the foam and interfacial friction factor by combining the small scale data with the data collected in the large scale experiments. These closure relationships are applicable to four different surfactants tested. A new transition criterion was developed and successfully used to predict onset of liquid loading under foam flow. Using a force balance over the gas core in annular flow, we developed a new procedure to calculate the pressure drop under foam flow conditions. We compared our model results with actual measurements in the large scale facility. Our model was reasonably able to predict the pressure drop within ±30%. The reason for such a large variance is that the small scale facility was not able to capture all the characteristics of the foam which were observed in the large scale facility. It is very difficult to reproduce the foam characteristics exactly in two different experiments. This is discussed further in this paper.The procedure developed is the only one currently available to calculate the pressure drop under the foam flow conditions using the small scale data. It is superior to conventional annular flow pressure drop prediction models which are currently available in the literature.     

Simulation of Entry-Region Flow in Open-Cell Metal Foam and Experimental Validation

Open-cell metal foam is distinguished from traditional porous media by its very high porosities (often greater than 90 %), and its web-like open structure and good permeability. As such, the foam is a very attractive core for many engineered systems, e.g., heat exchangers, filtration devices, catalysts, and reactors. The flow field inside the foam is rather complex due to flow reversal and vigorous mixing. This complexity is increased by the possible presence of an entry region. The entrance region in metal foam is usually underestimated and ignored, just like its counterpart in traditional porous media. In this paper, the actual entry length is determined by simulation and direct experiment on commercial open-cell aluminum foam. It is shown to be dependent on flow velocity and to reach a constant value for higher velocities. The complex and intrinsically random architecture of the foam is idealized using a unit geometrical model, in order to numerically investigate the flow field and pressure drop inside the foam. The Navier–Stokes equations are solved directly, and velocity and pressure fields are obtained for various approach velocities using a commercial numerical package. The entry length is ascertained from the behavior of the velocity field close to the entrance. Comparisons to experimental data were also carried out. The commercial foam that was used in the experiment had 10 ppi and porosity of 91.2 %. Air was forced to flow inside the foam using an open-loop wind tunnel. Good qualitative agreement between the modeling and experimental results are obtained. The agreement lends confidence to the modeling approach and the determined entry length.     

Hydrodynamics of viscous fluid through porous slit with linear absorption

The exact solutions for the viscous fluid through a porous slit with linear ab-sorption are obtained. The Stokes equation with non-homogeneous boundary conditions is solved to get the expressions for the velocity components, pressure distribution, wall shear stress, fractional absorption, and leakage flux. The volume flow rate and mean flow rate are found to be useful in obtaining a convenient form of the longitudinal velocity component and pressure difference. The points of the maximum velocity components for a fixed axial distance are identified. The value of the linear absorption parameter is ran-domly chosen, and the rest available data of the rat kidney to the tabulate pressure drop and fractional absorption are incorporated. The effects of the linear absorption, uniform absorption, and flow rate parameters on the flow properties are discussed by graphs. It is found that forward flow occurs only if the volume flux per unit width is greater than the absorption velocity throughout the length of the slit, otherwise back flow may occur. The leakage flux increases with the increase in the linear absorption parameter. Streamlines are drawn to help the analysis of the flow behaviors during the absorption of the fluid flow through the renal tubule and purification of blood through an artificial kidney.     

Transport modeling of sedimenting particles in a turbulent pipe flow using Euler–Lagrange large eddy simulation

A volume-filtered Euler–Lagrange large eddy simulation methodology is used to predict the physics of turbulent liquid–solid slurry flow through a horizontal periodic pipe. A dynamic Smagorinsky model based on Lagrangian averaging is employed to account for the sub-filter scale effects in the liquid phase. A fully conservative immersed boundary method is used to account for the pipe geometry on a uniform cartesian grid. The liquid and solid phases are coupled through volume fraction and momentum exchange terms. Particle–particle and particle–wall collisions are modeled using a soft-sphere approach. Three simulations are performed by varying the superficial liquid velocity to be consistent with the experimental data by Dahl et al. (2003). Depending on the liquid flow rate, a particle bed can form and develop different patterns, which are discussed in light of regime diagrams proposed in the literature. The fluctuation in the height of the liquid-bed interface is characterized to understand the space and time evolution of these patterns. Statistics of engineering interest such as mean velocity, mean concentration, and mean streamwise pressure gradient driving the flow are extracted from the numerical simulations and presented. Sand hold-up calculated from the simulation results suggest that this computational strategy is capable of predicting critical deposition velocity.     

Reynolds Stress Budgets in Couette and Boundary Layer Flows

Reynolds stress budgets for both Couette and boundary layer flows are evaluated and presented. Data are taken from direct numerical simulations of rotating and non-rotating plane turbulent Couette flow and turbulent boundary layer with and without adverse pressure gradient. Comparison of the total shear stress for the two types of flows suggests that the Couette case may be regarded as the high Reynolds number limit for the boundary layer flow close to the wall. The limit values of turbulence statistics close to the wall for the boundary layer for increasing Reynolds number approach the corresponding Couette flow values. The direction of rotation is chosen so that it has a stabilizing effect, whereas the adverse pressure gradient is destabilizing. The pressure-strain rate tensor in the Couette flow case is presented for a split into slow, rapid and Stokes terms. Most of the influence from rotation is located to the region close to the wall, and both the slow and rapid parts are affected. The anisotropy for the boundary layer decreases for higher Reynolds number, reflecting the larger separation of scales, and becomes close to that for Couette flow. The adverse pressure gradient has a strong weakening effect on the anisotropy. All of the data presented here are available on the web [36]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Counter-current three-phase fluidization in a turbulent contact absorber: A CFD simulation

A computational fluid dynamics study of three-phase counter-current fluidization occurring in a turbulent contact absorber was performed. A two-dimensional, transient Eulerian multi-fluid model was used, in which the dispersed solid phase was modeled employing a kinetic theory of granular flow. The grid independence of the model, the effect of wall boundary conditions, the choice of granular temperature model, the effects of order of discretization scheme and drag models were studied for a base case setting. The results of simulations were validated against experimental results obtained from the literature. Once the model settings were finalized, simulations were performed for different gas and liquid velocities to predict the hydrodynamics of the absorber. Computed bed expansion and pressure drop were compared with experimental data. Good agreement between the two was found for low velocities of gas and liquid.     

The prediction of stratified two-phase flow with a two-equation model of turbulence

A two-equation model is applied to a stratified two-phase flow system to predict turbulent transport mechanisms in both phases.In the present analysis, the effects of interfacial waves on the flow field are formulated in terms of boundary conditions for the gas-liquid interface. For the gas phase, the wavy interface has such flow separation effects as a rough surface in a single-phase flow. While for the liquid phase, the waves generate turbulant energy which is transported progressively toward a lower wall region. The analytical results are in good agreement with available data regarding pressure drop, holdup and velocity profiles.     

A hybrid numerical approach to predict the vibrational responses of panels excited by a turbulent boundary layer

In this work, a hybrid numerical approach to predict the vibrational responses of planar structures excited by a turbulent boundary layer is presented. The approach combines an uncorrelated wall plane wave technique with the finite element method. The wall pressure field induced by a turbulent boundary layer is obtained as a set of uncorrelated wall pressure plane waves. The amplitude of these plane waves are determined from the cross spectrum density function of the wall pressure field given either by empirical models from literature or from experimental data. The response of the planar structure subject to a turbulent boundary layer excitation is then obtained from an ensemble average of the different realizations. The numerical technique is computationally efficient as it rapidly converges using a small number of realizations. To demonstrate the method, the vibrational responses of two panels with simply supported or clamped boundary conditions and excited by a turbulent flow are considered. In the case study comprising a plate with simply supported boundary conditions, an analytical solution is employed for verification of the method. For both cases studies, numerical results from the hybrid approach are compared with experimental data measured in two different anechoic wind tunnels.     

Low-Prandtl-number Bénard–Marangoni convection in a vertical magnetic field

The effect of a homogeneous magnetic field on surface-tension-driven Bénard convection is studied by means of direct numerical simulations. The flow is computed in a rectangular domain with periodic horizontal boundary conditions and the free-slip condition on the bottom wall using a pseudospectral Fourier–Chebyshev discretization. Deformations of the free surface are neglected. Two- and three-dimensional flows are computed for either vanishing or small Prandtl number, which are typical of liquid metals. The main focus of the paper is on a qualitative comparison of the flow states with the non-magnetic case, and on the effects associated with the possible near-cancellation of the nonlinear and pressure terms in the momentum equations for two-dimensional rolls. In the three-dimensional case, the transition from a stationary hexagonal pattern at the onset of convection to three-dimensional time-dependent convection is explored by a series of simulations at zero Prandtl number.     

Analysis of Differential Equations Modelling the Reactive Flow through a Deformable System of Cells

A system of model equations coupling fluid flow, deformation of solid structure and chemical reactions is formulated starting from processes in biological tissue. The main aim of this paper is to analyse this non-standard system, where the elasticity modules are functionals of a concentration and the diffusion coefficients of the chemical substances are functions of their concentrations. A new approach and new methods are required and adapted to these nonlinearities and the transmission conditions on the interface solid–fluid. Strong solutions for the initial and boundary value problem are constructed under suitable regularity assumptions on the data, and stability estimates of the solutions with respect to the initial and boundary values are proved. These estimates imply uniqueness directly. The approach of the paper can be used in more general problems modeling reactive flow and transport and its interaction with elastic cell structures. In a forthcoming paper the approach of this paper is used for getting the upscaled system modeling reactive flow through biological tissue on the macroscopic scale, starting from a system on the cell level.     

Onset of Mobilization and the Fraction of Trapped Foam in Porous Media

Usually, foam in a porous medium flows through a small and spatially varying fraction of available pores, while the bulk of it remains trapped. The trapped foam is under a pressure gradient corresponding to the pressure gradient imposed by the flowing foam and continuous wetting liquid. The imposed pressure gradient and coalescence of the stationary foam lamellae periodically open flow channels in the trapped foam region. Foam lamellae in each of these channels flow briefly, but channels are eventually plugged by smaller bubbles entering into the trapped region. The result is a cycling of flow channels that open and close throughout the trapped foam, leading to intermittent pulsing of foam flow in that region.The dynamic behavior of foam trapped in porous media is modeled here with a pore network simulator. We predict the magnitude of the pressure drop leading to the onset of flow of foam lamellae in the region containing trapped foam. This mobilization pressure drop depends only on the number of lamellae in the flow path and on the geometry of the pores that make up this path.The principles learned in this study allow us to predict the fraction of foam that is trapped in a porous medium under given flow conditions. We present here the first analytic expression for the trapped foam fraction as a function of the pressure gradient, and of the mean and standard deviation of the pore size distribution. This expression provides a missing piece for the continuum foam flow models based on the moments of the volume-averaged population balance of foam bubbles.     

Multiscale Modeling of Colloid and Fluid Dynamics in Porous Media Including an Evolving Microstructure

We consider colloidal dynamics and single-phase fluid flow within a saturated porous medium in two space dimensions. A new approach in modeling pore clogging and porosity changes on the macroscopic scale is presented. Starting from the pore scale, transport of colloids is modeled by the Nernst?CPlanck equations. Here, interaction with the porous matrix due to (non-)DLVO forces is included as an additional transport mechanism. Fluid flow is described by incompressible Stokes equations with interaction energy as forcing term. Attachment and detachment processes are modeled by a surface reaction rate. The evolution of the underlying microstructure is captured by a level set function. The crucial point in completing this model is to set up appropriate boundary conditions on the evolving solid?Cliquid interface. Their derivation is based on mass conservation. As a result of an averaging procedure by periodic homogenization in a level set framework, on the macroscale we obtain Darcy??s law and a modified averaged convection?Cdiffusion equation with effective coefficients due to the evolving microstructure. These equations are supplemented by microscopic cell problems. Time- and space-dependent averaged coefficient functions explicitly contain information of the underlying geometry and also information of the interaction potential. The theoretical results are complemented by numerical computations of the averaged coefficients and simulations of a heterogeneous multiscale scenario. Here, we consider a radially symmetric setting, i.e., in particular we assume a locally periodic geometry consisting of circular grains. We focus on the interplay between attachment and detachment reaction, colloidal interaction forces, and the evolving microstructure. Our model contributes to the understanding of the effects and processes leading to porosity changes and pore clogging from a theoretical point of view.     

Direct modeling of flow of FENE fluids

Direct simulations of macromolecular fluids are carried out for flows between parallel plates and in expanding and contracting channels. The macromolecules are modeled as FENE dumbbells with soft disks or Lennard-Jones dumbbell-dumbbell interactions. The results are presented in terms of profiles and contour plots of velocity, pressure, temperature, density, and flow fields. In addition the data for potential energy, shear stress, and the normal components of the stress tensor are collected. In general, an excellent agreement is found between the simulated profiles and the well-known flow structures, such as flow separation and formation of viscous eddies, indicating that micro-hydrodynamics is a viable tool in linking macroscopic phenomena with the underlying physical mechanisms. The simulations are performed in the Newtonian regime, for medium-size systems comprising up to 3888 dumbbells. This number is sufficiently large to control boundary and particle number effects. The flow is induced by gravity. The traditional stochastic (thermal) and periodic boundary conditions are employed. Also, diffusive boundary conditions, which could include a stagnant fluid layer and repulsive potential walls, are developed. The scaling problems, which are related to the application of a large external force in a microscopic system (of the size of the order 100 Å), result in extreme pressure and temperature gradients. In addition, the viscosity and thermal conductivity coefficients obtained from velocity and temperature profiles of the channel flow are presented. These results are confirmed independently from modeling of Couette flow by the SLLOD equations of motion and from the Evans algorithm for thermal conductivity.