Multiple 3D flaws in coupled electro-magneto-thermo-elastic multiphase composites by extended hypersingular integral equation method

This work presents extended hypersingular integral equation (E-HIE) method to analyze the multiple 3D mixed-mode flaws problem in fully coupled electro-magneto-thermo-elastic multiphase composites under extended electro-magneto-thermo-elastic coupled loads through intricate theoretical analysis and numerical simulations. First, the problem is reduced to solving a set of E-HIEs. Analytical solutions for the extended singular stresses, the extended stress intensity factors (E-SIFs), the extended energy release rate and the extended strain energy density factors (E-SEDFs) near the flaws front are obtained. Then, the numerical method for the E-HIEs for two 3D flaws subjected to extended coupled loads is proposed. Finally, numerical solutions of E-SIFs and E-SEDFs of some examples are given, and the effect of flaws orientation, interaction and shielding is discussed.     

Mixed-mode stress intensity factors of 3D interface crack in fully coupled electromagnetothermoelastic multiphase composites

This contribution presents an extended hypersingular intergro-differential equation (E-HIDE) method for modeling the 3D interface crack problem in fully coupled electromagnetothermoelastic anisotropic multiphase composites under extended electro-magneto-thermo-elastic coupled loads through theoretical analysis and numerical simulations. First, based on the extended boundary element method, the 3D interface crack problem is reduced to solving a set of E-HIDEs coupled with extended boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended singular stress indices around the interface crack front terminating at the interface is analyzed by the extended main-part analysis. The extended stress intensity factors near the crack front are defined. In addition, a numerical method for a 3D interface crack problem subjected to extended loads is proposed, in which the extended displacement discontinuities are approximated by the product of basic density functions and polynomials. Finally, the radiation distribution of extended stress intensity factors at the interface crack surface are calculated, and the results are presented toward demonstrating the applicability of the proposed method.     

Modelling Biogrout: A New Ground Improvement Method Based on Microbial-Induced Carbonate Precipitation

Biogrout is a new soil reinforcement method based on microbial-induced carbonate precipitation. Bacteria are placed and reactants are flushed through the soil, resulting in calcium carbonate precipitation, causing an increase in strength and stiffness of the soil. Due to this precipitation, the porosity of the soil decreases. The decreasing porosity influences the permeability and therefore the flow. To analyse the Biogrout process, a model was created that describes the process. The model contains the concentrations of the dissolved species that are present in the biochemical reaction. These concentrations can be solved from a advection?Cdispersion?Creaction equation with a variable porosity. Other model equations involve the bacteria, the solid calcium carbonate concentration, the (decreasing) porosity, the flow and the density of the fluid. The density of the fluid changes due to the biochemical reactions, which results in density driven flow. The partial differential equations are solved by the Standard Galerkin finite-element method. Simulations are done for some 1D and 2D configurations. A 1D configuration can be used to model a column experiment and a 2D configuration may correspond to a sheet or a cross section of a 3D configuration.     

Construction of hybrid 1D-0D networks for efficient and accurate blood flow simulations

The one-dimensional (1D) modeling of blood flow in complex networks of vessels and cardiovascular models can result in computationally expensive simulations. The complexity of such networks has significantly increased in the last years, in terms of both enhanced anatomical detail and modeling of physiological mechanisms and mechanical characteristics. To address such issue, the main goal of this work is to present a novel methodology to construct hybrid networks of coupled 1D and 0D vessels and to perform computationally efficient and accurate blood flow simulations in such networks. Departing from both the 1D and lumped-parameter (0D) nonlinear models for blood flow, we propose high-order numerical coupling strategies to solve the 1D, 0D, and hybrid coupling of vessels at junctions. To effectively construct hybrid networks, we explore different a-priori model selection criteria focusing in obtaining the best possible trade-off between computational cost of the simulations and accuracy of the computed solutions for the hybrid network with respect to the 1D network. The achievement of the expected order of accuracy is verified in several test cases. The novel methodology is applied to two different arterial networks, the 37-artery network and the reduced ADAN56 model, where, in order to identify the best performing a-priori model selection criteria, the quantitative assessment of CPU times and errors and the qualitative comparison between results are carried out and discussed.     

A computationally efficient 3D finite‐volume scheme for violent liquid–gas sloshing

We describe a semi‐implicit volume‐of‐fluid free‐surface‐modelling methodology for flow problems involving violent free‐surface motion. For efficient computation, a hybrid‐unstructured edge‐based vertex‐centred finite volume discretisation is employed, while the solution methodology is entirely matrix free. Pressures are solved using a matrix‐free preconditioned generalised minimum residual algorithm and explicit time‐stepping is employed for the momentum and interface‐tracking equations. The high resolution artificial compressive (HiRAC) volume‐of‐fluid method is used for accurate capturing of the free surface in violent flow regimes while allowing natural applicability to hybrid‐unstructured meshes. The code is parallelised for solution on distributed‐memory architectures and evaluated against 2D and 3D benchmark problems. Good parallel scaling is demonstrated, with almost linear speed‐up down to 6000 cells per core. Finally, the code is applied to an industrial‐type problem involving resonant excitation of a fuel tank, and a comparison with experimental results is made in this violent sloshing regime. Copyright © 2015 John Wiley & Sons, Ltd.     

Lattice Boltzmann modeling of pore‐scale fluid flow through idealized porous media

In order to capture the hydro‐mechanical impacts on the solid skeleton imposed by the fluid flowing through porous media at the pore‐scale, the flow in the pore space has to be modeled at a resolution finer than the pores, and the no‐slip condition needs to be enforced at the grain–fluid interface. In this paper, the lattice Boltzmann method (LBM), a mesoscopic Navier–Stokes solver, is shown to be an appropriate pore‐scale fluid flow model. The accuracy and lattice sensitivity of LBM as a fluid dynamics solver is demonstrated in the Poiseuille channel flow problem (2‐D) and duct flow problem (3‐D). Well‐studied problems of fluid creeping through idealized 2‐D and 3‐D porous media (J. Fluid Mech. 1959; 5 (2):317–328, J. Fluid Mech. 1982; 115 :13–26, Int. J. Multiphase Flow 1982; 8 (4):343–360, Phys. Fluids A 1989; 1 (1):38–46, Int. J. Numer. Anal. Meth. Geomech. 1999; 23 :881–904, Int. J. Numer. Anal. Meth. Geomech. 2010; DOI: 10.1002/nag.898, Int. J. Multiphase Flow 1982; 8 (3):193–206) are then simulated using LBM to measure the friction coefficient for various pore throats. The simulation results agree well with the data reported in the literature. The lattice sensitivity of the frictional coefficient is also investigated. Copyright © 2010 John Wiley & Sons, Ltd.     

人工黏性增强三维移动粒子半隐算法稳定性

通过综合比较移动粒子半隐方法(moving-particle semi-implicitmethod, MPS)各种稳定性计算方法, 考虑了多种可能导致计算不稳定的因素, 并首次提出加入人工黏性来抑制非物理压力振荡, 得到了较为稳定的三维MPS算法.采用各种稳定性方案对比模拟了三维立方流体旋转状态下的变形及三维静水压问题, 并进一步与商业软件的流体体积函数方法(volumeof fluid, VOF)模型计算结果对比, 验证了其正确性.结果表明:对三维情形, 单纯应用已有的稳定性算法仍难以满足MPS计算稳定性的要求, 而进一步辅以该文提出的人工黏性方法则可以在准确性及稳定性方面均取得较好的效果, 且计算结果显示, 通过该文提出的开关模式增加人工黏性并不影响压强泊松方程的求解精度, 因而不会对流场的求解产生负面作用.     

A finite difference method for 3D flows about bodies of complex geometry in rectangular co-ordinate systems

A finite difference simulation method is developed for 3D flow about a body of complex geometry. The Navier–Stokes equation is approximated by a high-order-accurate difference scheme in the framework of rectangular co-ordinate systems. The configuration of the 3D body is represented by use of both surface porosity and volume porosity and the no-slip body boundary conditions are approximately implemented on the boundary cells. The validity of the method is demonstrated by a numerical test of flow past a sphere at a Reynolds number of 1000. The complicated structure of separated vortices is well revealed by this test computation. The versatility of the method is shown by application to an ocean-engineering problem of flow about a bay with an island.     

Analysis of velocity equation of steady flow of a viscous incompressible fluid in channel with porous walls

Steady flow of a viscous incompressible fluid in a channel, driven by suction or injection of the fluid through the channel walls, is investigated. The velocity equation of this problem is reduced to nonlinear ordinary differential equation with two boundary conditions by appropriate transformation and convert the two‐point boundary‐value problem for the similarity function into an initial‐value problem in which the position of the upper channel. Then obtained differential equation is solved analytically using differential transformation method and compare with He's variational iteration method and numerical solution. These methods can be easily extended to other linear and nonlinear equations and so can be found widely applicable in engineering and sciences. Copyright © 2009 John Wiley & Sons, Ltd.     

Onset of instability of a flag in uniform flow

This paper numerically and analytically studies the onset of instability of a flag in uniform flow. The three-dimensional (3D) simulation is performed by using an immersed-boundary method coupled with a nonlinear finite element method. The global stability, bistability and instability are identified in the 3D simulations. The Squire's theorem is extended to analyze the stability of the fluid-flag system with 3D initial perturbations. It is found that if a parallel flow around the flag admits an unstable 3D disturbance for a certain value of the flutter speed, then a two-dimensional (2D) disturbance at a lower flutter speed is also admitted. In addition, the growth rate of 2D disturbance is larger than that of the 3D disturbance.     

Homogenized modeling for vascularized poroelastic materials

A new mathematical model for the macroscopic behavior of a material composed of a poroelastic solid embedding a Newtonian fluid network phase (also referred to as vascularized poroelastic material), with fluid transport between them, is derived via asymptotic homogenization. The typical distance between the vessels/channels (microscale) is much smaller than the average size of a whole domain (macroscale). The homogeneous and isotropic Biot’s equation (in the quasi-static case and in absence of volume forces) for the poroelastic phase and the Stokes’ problem for the fluid network are coupled through a fluid-structure interaction problem which accounts for fluid transport between the two phases; the latter is driven by the pressure difference between the two compartments. The averaging process results in a new system of partial differential equations that formally reads as a double poroelastic, globally mass conserving, model, together with a new constitutive relationship for the whole material which encodes the role of both pore and fluid network pressures. The mathematical model describes the mutual interplay among fluid filling the pores, flow in the network, transport between compartments, and linear elastic deformation of the (potentially compressible) elastic matrix comprising the poroelastic phase. Assuming periodicity at the microscale level, the model is computationally feasible, as it holds on the macroscale only (where the microstructure is smoothed out), and encodes geometrical information on the microvessels in its coefficients, which are to be computed solving classical periodic cell problems. Recently developed double porosity models are recovered when deformations of the elastic matrix are neglected. The new model is relevant to a wide range of applications, such as fluid in porous, fractured rocks, blood transport in vascularized, deformable tumors, and interactions across different hierarchical levels of porosity in the bone.     

3D phase field modeling of electrohydrodynamic multiphase flows

A 3D phase field model is developed to investigate the electrohydrodynamic (EHD) two phase flows. The explicit finite difference method, enhanced by parallel computing, is employed to solve the coupled nonlinear governing equations for the electric field, the fluid flow field and free surface deformation. Numerical tests indicate that an appropriate interpolation of densities within the interface is critical in ensuring numerical stability for highly stratified flows. The 3D phase field model compares well with the Taylor theory for the deformation of a single dielectric droplet in an electric field. Computed results show that the deformation of a leaky dielectric droplet in an electric field undergoes various stages before it reaches the final oblate shape. This is caused by the free charge relaxation near the fluid–fluid interface. The coalescence of four droplets in an electric field illustrates a truly 3D deformation behavior and a complex evolving fluid flow field associated with the participating droplets. The coalescence is a result of combined actions produced by the global electric force, the circulatory flows generated by the local electrohydrodynamic stress and the electrically-induced deformation. The 3D phase field model is also applied in modeling of an electrohydrodynamic patterning process for manufacturing nanoscaled structures, in which complex 3D flow structures develop as the electrically-induced deformation evolves.     

3-D inviscid self-excited vibrations of a blade row in the last stage turbine

Presented here is a three-dimensional (3-D) nonlinear time-marching method for the aeroelastic behaviour of an oscillating turbine blade row. The approach has been based in the solution of a coupled fluid–structure problem where the aerodynamic and structural dynamic equations are integrated simultaneously in time. This provides the correct formulation of a coupled problem, as the interblade phase angle (IBPA) at which stability (instability) would occur is also a part of solution. The ideal gas flow around multiple interblade passages (with periodicity in the entire annulus) is described by the unsteady Euler equations in conservative form, which are integrated by using the explicit monotonic second-order accurate Godunov–Kolgan finite-volume scheme and a moving hybrid H–O (or H–H) grid. The fluid and the structural equations are solved using the modal superposition method. An aeroelasticity prediction of a turbine blade of 0.765 m is presented. The natural frequencies and modal shapes of the blade were calculated by using 3-D finite element models. The instability regions for five mode shapes and the distribution of the aerodamping coefficient along the blade length were shown for harmonic oscillations with an assumed IBPA. The coupled fluid–structure oscillations in which the IBPA is part of the solution are shown.     

Implementation of some higher-order convection schemes on non-uniform grids

A generalized formulation is applied to implement the quadratic upstream interpolation (QUICK) scheme, the second-order upwind (SOU) scheme and the second-order hybrid scheme (SHYBRID) on non-uniform grids. The implementation method is simple. The accuracy and efficiency of these higher-order schemes on non-uniform grids are assessed. Three well-known bench mark convection-diffusion problems and a fluid flow problem are revisited using non-uniform grids. These are: (1) transport of a scalar tracer by a uniform velocity field; (2) heat transport in a recirculating flow; (3) two-dimensional non-linear Burgers equations; and (4) a two-dimensional incompressible Navier-Stokes flow which is similar to the classical lid-driven cavity flow. The known exact solutions of the last three problems make it possible to thoroughly evaluate accuracies of various uniform and non-uniform grids. Higher accuracy is obtained for fewer grid points on non-uniform grids. The order of accuracy of the examined schemes is maintained for some tested problems if the distribution of non-uniform grid points is properly chosen.     

Numerical simulations of viscous flows using a meshless method

This paper uses the element‐free Galerkin (EFG) method to simulate 2D, viscous, incompressible flows. The control equations are discretized with the standard Galerkin method in space and a fractional step finite element scheme in time. Regular background cells are used for the quadrature. Several classical fluid mechanics problems were analyzed including flow in a pipe, flow past a step and flow in a driven cavity. The flow field computed with the EFG method compared well with those calculated using the finite element method (FEM) and finite difference method. The simulations show that although EFG is more expensive computationally than FEM, it is capable of dealing with cases where the nodes are poorly distributed or even overlap with each other; hence, it may be used to resolve remeshing problems in direct numerical simulations. Flows around a cylinder for different Reynolds numbers are also simulated to study the flow patterns for various conditions and the drag and lift forces exerted by the fluid on the cylinder. These forces are calculated by integrating the pressure and shear forces over the cylinder surface. The results show how the drag and lift forces oscillate for high Reynolds numbers. The calculated Strouhal number agrees well with previous results. Copyright © 2008 John Wiley & Sons, Ltd.     

A least‐square weighted residual method for level set formulation

In this paper, a least‐square weighted residual method (LSWRM) for level set (LS) formulation is introduced to achieve interface capturing in two‐dimensional (2D) and three‐dimensional (3D) problems. An LSWRM was adopted for two semi‐discretized advection and reinitialization equations of the LS formulation. The present LSWRM provided good mathematical properties such as natural numerical diffusion and the symmetry of the resulting algebraic systems for the advection and reinitialization equations. The proposed method was validated by solving some 2D and 3D benchmark problems such as those involving a rotating slotted disk, the rotation of a slotted sphere, and a time‐reversed single‐vortex flow and a deformation problem of a spherical fluid. The numerical results were compared with those obtained from essentially non‐oscillatory type formulations and particle LS methods. Further, the proposed LSWRM for the LS formulation was coupled with a splitting finite element method code to solve the incompressible Navier–Stokes equations, and then, the collapse of a 3D broken dam flow was well simulated; in the simulation, the entrapping of air and the splashing of the surge front of water were reproduced. The mass conservation of the present method was found to be satisfactory during the entire simulation. Copyright © 2011 John Wiley & Sons, Ltd.     

Formulation of thermo-hydro-mechanical coupling behavior of unsaturated soils based on hybrid mixture theory

Thermo-Hydro-Mechanical （THM） coupling pro- cesses in unsaturated soils are very important in both theoretical researches and engineering applications. A coupled formulation based on hybrid mixture theory is derived to model the THM coupling behavior of unsaturated soils. The free-energy and dissipative functions for different phases are derived from Taylor＇s series expansions. Constitutive relations for THM coupled behaviors of unsaturated soils, which include deformation, entropy change, fluid flow, heat conduction, and dynamic compatibility conditions on the interfaces, are then established. The number of field equations is shown to be equal to the number of unknown variables; thus, a closure of this coupling problem is established. In addition to modifications of the physical conservation equations with coupling effect terms, the constitutive equations, which consider the coupling between elastoplastic deformation of the soil skeleton, fluid flow, and heat transfer, are also derived.     

Finite Difference Method for the Acoustic Radiation of an Elastic Plate Excited by a Turbulent Boundary Layer: A Spectral Domain Solution

A finite difference method is developed to study, on a two-dimensional model, the acoustic pressure radiated when a thin elastic plate, clamped at its boundaries, is excited by a turbulent boundary layer. Consider a homogeneous thin elastic plate clamped at its boundaries and extended to infinity by a plane, perfectly rigid, baffle. This plate closes a rectangular cavity. Both the cavity and the outside domain contain a perfect fluid. The fluid in the cavity is at rest. The fluid in the outside domain moves in the direction parallel to the system plate/baffle with a constant speed. A turbulent boundary layer develops at the interface baffle/plate. The wall pressure fluctuations in this boundary layer generates a vibration of the plate and an acoustic radiation in the two fluid domains. Modeling the wall pressure fluctuations spectrum in a turbulent boundary layer developed over a vibrating surface is a very complex and unresolved task. Ducan and Sirkis [1] proposed a model for the two-way interactions between a membrane and a turbulent flow of fluid. The excitation of the membrane is modeled by a potential flow randomly perturbed. This potential flow is modified by the displacement of the membrane. Howe [2] proposed a model for the turbulent wall pressure fluctuations power spectrum over an elastomeric material. The model presented in this article is based on a hypothesis of one-way interaction between the flow and the structure: the flow generates wall pressure fluctuations which are at the origin of the vibration of the plate, but the vibration of the plate does not modify the characteristics of the flow. A finite difference scheme that incorporates the vibration of the plate and the acoustic pressure inside the fluid cavity has been developed and coupled with a boundary element method that ensures the outside domain coupling. In this paper, we focus on the resolution of the coupled vibration/interior acoustic problem. We compare the results obtained with three numerical methods: (a) a finite difference representation for both the plate displacement and the acoustic pressure inside the cavity; (b) a coupled method involving a finite difference representation for the displacement of the plate and a boundary element method for the interior acoustic pressure; (c) a boundary element method for both the vibration of the plate and the interior acoustic pressure. A comparison of the numerical results obtained with two models of turbulent wall pressure fluctuations spectrums - the Corcos model [3] and the Chase model [4] - is proposed. A difference of 20 dB is found in the vibro-acoustic response of the structure. In [3], this difference is explained by calculating a wavenumber transfer function of the plate. In [6], coupled beam-cavity modes for similar geometry are calculated by the finite difference method. This revised version was published online in July 2006 with corrections to the Cover Date.     

Development of a three-dimensional solver based on the local domain-free discretization and immersed boundary method and its application for incompressible flow problems

Recently, the author and two other coauthors have proposed a two-dimensional hybrid local domain-free discretization and immersed boundary method (LDFD-IBM), which can be used to solve the flow problem with complex geometries. In this paper, the LDFD-IBM is extended to solve a three-dimensional unsteady incompressible flow with the complex computational domain. The technical issues related to the implementation of the LDFD-IBM in three-dimensional problems are discussed in detail, particularly for the discretization of Navier-Stokes equations, mesh strategies for a three-dimensional flow, and the fast algorithm on the identification of the status of mesh nodes (ie, to identify if the mesh node is located in the solid domain, in the fluid domain, or near the immersed boundary). Numerical tests show that the LDFD-IBM can accurately solve three-dimensional incompressible problems with ease.