Coupling 3D and 1D fluid-structure interaction models for blood flow simulations

Three-dimensional (3D) simulations of blood flow in medium to large vessels are now a common practice. These models consist of the 3D Navier-Stokes equations for incompressible Newtonian fluids coupled with a model for the vessel wall structure. However, it is still computationally unaffordable to simulate very large sections, let alone the whole, of the human circulatory system with fully 3D fluid-structure interaction models. Thus truncated 3D regions have to be considered. Reduced models, one-dimensional (1D) or zero-dimensional (0D), can be used to approximate the remaining parts of the cardiovascular system at a low computational cost. These models have a lower level of accuracy, since they describe the evolution of averaged quantities, nevertheless they provide useful information which can be fed to the more complex model. More precisely, the 1D models describe the wave propagation nature of blood flow and coupled with the 3D models can act also as absorbing boundary conditions. We consider in this work the coupling of a 3D fluid-structure interaction model with a 1D hyperbolic model. We study the stability of the coupling and present some numerical results. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Coupling of the Reynolds Fluid-Film Equation with the 2D Navier-Stokes Equations

The pressure field in thin fluid films can quite precisely be calculated by Reynolds fluid-film equation. In some problems, it may be useful to couple thin fluid-films with general 2D or 3D fluid flows. In the current work, we analyze the fluid flow, pressure and temperature field in a hydrodynamic journal bearing with a rectangular oil groove. Pressure and temperature in the fluid gap are calculated by means of the Reynolds equation and the 2D energy equation. Cavitation effects are taken into account by incorporating a 2-phase cavitation approach. In order to calculate the velocity and pressure field in the oil groove, the 2D Navier-Stokes equations are used; the temperature distribution in the oil groove is computed by means of the 2D energy equation. Appropriate coupling conditions for velocity, pressure and temperature are formulated in order to couple the flow in the fluid gap with the flow in the oil groove. Thermal expansion of journal shaft and bearing housing are also taken into account, since the bearing clearance changes with increasing temperature. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Fluid-Structure Port-Hamiltonian Model for Incompressible Flows in Tubes with Time Varying Geometries

ABSTRACT

A simple and scalable finite-dimensional model based on the port-Hamiltonian framework is proposed to describe the fluid–structure interaction in tubes with time-varying geometries. For this purpose, the moving tube wall is described by a set of mass-spring-damper systems while the fluid is considered as a one-dimensional incompressible flow described by the average momentum dynamics in a set of incompressible flow sections. To couple these flow sections small compressible volumes are defined to describe the pressure between two adjacent fluid sections. The fluid-structure coupling is done through a power-preserving interconnection between velocities and forces. The resultant model includes external inputs for the fluid and inputs for external forces over the mechanical part that can be used for control or interconnection purposes. Numerical examples show the accordance of this simplified model with finite-element models reported in the literature.     

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this work, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.     

三维横向流体诱发直管振动的数值模拟

基于有限体积法和有限元法,结合动网格控制技术,建立了横向流体作用下三维弹性直管流致振动计算的数值模型,实现了计算结构动力学与计算流体力学之间的联合仿真．首先,通过对刚性管的静止绕流计算,研究了网格离散方式和不同湍流模型对圆柱类结构静止绕流流场特征的影响和预测能力,得到了适用于双向耦合分析的CFD模型;其次,利用基于双向流固耦合方法的流致振动模型,计算并分析了流体力与结构位移间的相位关系,指出流体力与位移间的相位差是由流体力引起的,同时对双向耦合和单向耦合进行了比较分析;最后通过对直管流致振动的数值计算,联合管表面压力、尾流区时均速度、分离角等时均量,分析了尾流区的流场特征.     

Convergence analysis of a projection semi-implicit coupling scheme for fluid–structure interaction problems

In this paper, we provide a convergence analysis of a projection semi-implicit scheme for the simulation of fluid–structure systems involving an incompressible viscous fluid. The error analysis is performed on a fully discretized linear coupled problem: a finite element approximation and a semi-implicit time-stepping strategy are respectively used for space and time discretization. The fluid is described by the Stokes equations, the structure by the classical linear elastodynamics equations (linearized elasticity, plate or shell models) and all changes of geometry are neglected. We derive an error estimate in finite time and we prove that the time discretization error for the coupling scheme is at least ${\sqrt{\delta t}}In this paper, we provide a convergence analysis of a projection semi-implicit scheme for the simulation of fluid–structure systems involving an incompressible viscous fluid. The error analysis is performed on a fully discretized linear coupled problem: a finite element approximation and a semi-implicit time-stepping strategy are respectively used for space and time discretization. The fluid is described by the Stokes equations, the structure by the classical linear elastodynamics equations (linearized elasticity, plate or shell models) and all changes of geometry are neglected. We derive an error estimate in finite time and we prove that the time discretization error for the coupling scheme is at least ?{dt}{\sqrt{\delta t}}. Finally, some numerical experiments that confirm the theoretical analysis are presented.     

A 3D non-Newtonian fluid-structure interaction model for blood flow in arteries

The mathematical modelling and numerical simulation of the human cardiovascular system is playing nowadays an important role in the comprehension of the genesis and development of cardiovascular diseases. In this paper we deal with two problems of 3D modelling and simulation in this field, which are very often neglected in the literature. On the one hand blood flow in arteries is characterized by travelling pressure waves due to the interaction of blood with the vessel wall. On the other hand, blood exhibits non-Newtonian properties, like shear-thinning, viscoelasticity and thixotropy. The present work is concerned with the coupling of a generalized Newtonian fluid, accounting for the shear-thinning behaviour of blood, with an elastic structure describing the vessel wall, to capture the pulse wave due to the interaction between blood and the vessel wall. We provide an energy estimate for the coupling and compare the numerical results with those obtained with an equivalent fluid-structure interaction model using a Newtonian fluid.     

The coupling system of Navier–Stokes equations and elastic Navier–Lame equations in a blood vessel

In this paper, the blood flow problem is considered in a blood vessel, and a coupling system of Navier–Stokes equations and linear elastic equations, Navier–Lame equations, in a cylinder with cylindrical elastic shell is given as the governing equations of the problem. We provide two finite element models to simulating the three-dimensional Navier–Stokes equations in the cylinder while the asymptotic expansion method is used to solving the linearly elastic shell equations. Specifically, in order to discrete the Navier–Stokes equations, the dimensional splitting strategy is constructed under the cylinder coordinate system. The spectral method is adopted along the rotation direction while the finite element method is used along the other directions. By using the above strategy, we get a series of two-dimensional-three-components (2D-3C) fluid problems. By introduce the S-coordinate system in E³ and employ the thickness of blood vessel wall as the expanding parameter, the asymptotic expansion method can be established to approximate the solution of the 3D elastic problem. The interface contact conditions can be treated exactly based on the knowledge of tensor analysis. Finally, numerical test shows that our method is reasonable.     

A Lagrangian-based SPH-DEM model for fluid–solid interaction with free surface flow in two dimensions

A Lagrangian-based SPH-DEM coupling model is proposed to study fluid–solid interaction (FSI) problems with free-surface flow. In this model, SPH uses an incompressible divergence-free scheme for simulating complex flow problems. Based on the Mohr–Coulomb criterion with tension cut, the DEM describes the characteristics of solid deformation and failure by means of contact models between particles. The coupling mechanism between SPH and DEM is realised by the decoupling of the force field during the process of fluid–solid interaction. That is, the motions of fluid and solid particles are reflected by the Navier–Stokes equations and interactions among solid particles are determined by Newton's second law in the DEM. To demonstrate the applicability of the SPH-DEM model, three case studies are used to verify the different fluid interaction situations with rigid bodies, deformable objects, and granular assemblies, respectively. The results of the proposed model shows good agreement with experimental data and indicates that it is capable of capturing the features of solid movement, deformation and failure under complex flow conditions with convincing accuracy and high efficiency.     

A beam segment element for dynamic analysis of large aqueducts

Large aqueduct structure is a complex structure that is commonly used in hydraulic engineering, especially in large-scale water conveying projects. The analysis of dynamic response for an aqueduct structure is necessary if the aqueduct is built in an earthquake area. Traditional 3D finite element analysis is time consuming and the existing simplified response method cannot take into account all the effects, such as the bending-torsion coupling effect and the constrained torsion, of the deformations of the thin wall structure of the aqueduct body. For this special structure, a simple and yet accurate model for dynamic analysis is needed. In this paper, a beam segment element is developed and used for the calculation of dynamic response for aqueduct structures. With the frame of the aqueduct being modeled using beam element, the proposed model can calculate the dynamic response of the whole aqueduct structures. Results are compared with that of a general purpose finite element analysis software using 3D finite element model. Good agreement is achieved between the two models. However, the proposed model needs less elements and much less computing time.     

Phase-Field Modeling of Hydraulic Fracture

Hydraulically driven fracture has gained more and more research activity in the last few years, especially due to the growing interest of the petroleum industry. Key challenge for a powerful simulation of this scenario is an effective modeling and numerical implementation of the behavior of the solid skeleton and the fluid phase, the mechanical coupling between the two phases as well as the incorporation of the fracture process. Existing models for hydraulic fracturing can be found for example in [1], where the crack path is predetermined, or in [2] who use a phase field fracture model in an elastic framework, however without incorporating the fluid flow. In this work we propose a new compact model structure for the Biot-type fluid transport in porous media at finite strains based on only two constitutive functions, that is the free energy function ψ and a dissipation potential ϕ that includes the incorporation of an additional Poiseuille-type fluid flow in cracks. This formulation is coupled to a phase field approach for fracture and is fully variational in nature, as shown in [3]. In contrast to formulations with a sharp-crack discontinuity, the proposed regularized approach has the main advantage of a straight-forward modeling of complex crack patterns including branching. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Strongly coupling of partitioned fluid–solid interaction solvers using reduced-order models

In this work a powerful technique is described which allows the implicit coupling of partitioned solvers in fluid–structure interaction (FSI) problems. The flow under consideration is governed by the Navier–Stokes equations for incompressible viscous fluids and modeled with the finite volume method. The structure is represented by a finite element formulation. The method allows the use of a black box fluid and structural solver because it builds up a reduced order model of the fluid and structural problem during the coupling process. Each solution of the fluid/structural solver in the coupling process can be seen as a sensitivity response of an applied displacement/pressure mode. The applied modes and their responses are used to build up a reduced-order model. The proposed model is used to predict the unsteady flow fields of a particular flow-induced vibrational phenomenon – a fixed cubic rigid body is submerged in an incompressible fluid flow (water), an elastic plate is attached to the rigid body in the centre of the downstream face, and the vortices, which separate from the corners of the rigid body upstream, generate lift forces which excite continuous oscillations of the elastic plate downstream. The computational results show that a fairly good convergence solution is achieved by using the reduced-order model that is based on only a few displacement and stress modes, which largely reduces the computational cost, compared with traditional approaches. At the same time, comparison of the numerical results of the model with available experimental data validates the methodology and assesses its accuracy.     

Evaluating flood extent mapping of two hydraulic models, 1D HEC‐RAS and 2D LISFLOOD‐FP in comparison with aerial imagery observations in Gorgan flood plain,Iran

We compared flood mapping techniques using a one‐dimensional (1D) hydraulic model HEC‐RAS and two‐dimensional (2D) LISFLOOD‐FP for a 10‐km reach of Gorgan River in Iran. Both models were run using the same hydrologic input data. The input into the models was a steady discharge of 90 cm, corresponds to a flood peak occurred on March 25, 2012. Flood maps generated using these two models were compared with an observed flood inundation map, using F‐statistic. The roughness coefficients of the models were calibrated by maximizing the value of the F‐statistic. Based on the F‐statistic, LISFLOOD‐FP gives a slightly better result (F = 0.69) than HEC‐RAS (F = 0.67). Visual comparison of the flood extents generated by the two models showed reasonably good agreement. Validation was done using a flood event occurred on May 31, 2014. The LISFLOOD‐FP model gave a better result for validation as well. The 2D model showed more consistency in comparison with the 1D model.     

Global existence of classical solutions for the 3D generalized compressible Oldroyd-B model

In this paper, we prove the global existence of small classical solutions to the 3D generalized compressible Oldroyd-B system. It can be seen as compressible Euler equations coupling the evolution of stress tensor τ. The result mainly shows that singularity of solutions to compressible Euler equations can be prevented by the coupling of viscoelastic stress tensor. Moreover, unlike most complex fluids containing compressible Euler equations, the irrotational condition ∇×u=0 would not be proposed here to achieve the global well-posedness.     

Numerical solution of 2D and 3D backward facing step flows

The work deals with numerical modelling of flow through 2-dimensional (2D) and 3-dimensional (3D) backward facing step. In laminar case, we apply several higher order upwind and central discretizations and compare numerical results with measurements. The turbulent regime is considered in 2D as well as in 3D and influence of secondary flow is observed. Different modifications of low-Re two equation turbulence models and an explicit algebraic Reynolds stress model (EARSM) are considered. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Solution Properties of a 3D Stochastic Euler Fluid Equation

We prove local well-posedness in regular spaces and a Beale–Kato–Majda blow-up criterion for a recently derived stochastic model of the 3D Euler fluid equation for incompressible flow. This model describes incompressible fluid motions whose Lagrangian particle paths follow a stochastic process with cylindrical noise and also satisfy Newton’s second law in every Lagrangian domain.

     

Numerical study of forebody wake effect on axisymmetric parachute opening shock and drag reduction

ABSTRACT

Parachute–forebody distance is a parameter which is amongst the most critical factors to be considered in forebody wake effect. In this study, a new axisymmetric parachute–forebody coupling model is developed. Axisymmetric wrinkling membrane element is built to assess the dynamic response of the parachute canopy membrane under fluid pressure. Besides, fluid model and its further implementation on the fluid structure analysis are discussed. With the proposed method, the wake effect on both the opening shock during inflation state and the drag reduction during steady state can be obtained efficiently. Finally, numerical model is validated with published experimental result and further employed to investigate the influence of distance parameters on fluid–parachute coupling behaviour. On the basis of numerical results, failure distance during the inflation process and critical forebody–parachute distance are determined. The results show that forebody–parachute distance has a strong influence on flow behaviour around the parachute in both inflation state and steady descent state.     

Development of POD reduced-order model and its closure scheme for 2D Rayleigh–Bénard convection

Reduced-order model (ROM) based on proper orthogonal decomposition (POD) is a fast computational fluid dynamics (CFD) method and has been widely applied to pure flow or heat conduction problems in the past. In this paper, the typical 2D Rayleigh–Bénard convection (RBC) in a square cavity was set as a research target. Firstly, the POD-ROM of 2D RBC problem was constructed at Ra = 10⁷, Pr = 0.71. Combining with direct numerical simulation (DNS) databases, a closure model (CM) was then proposed to correct the evolution process of POD-ROM. Based on the proposed CM, we realized the prediction of flow evolution for a new flow case under the parameters different from that used to get its POD eigenmodes. It showed that the proposed POD-ROM with CM could be able to predict the dynamics of new flow cases. Moreover, the corresponding method proposed in the present study can be also easily extended to other types of flow-heat coupling problems, such as natural heat convection, etc.     

Stability of Blowup for a 1D Model of Axisymmetric 3D Euler Equation

The question of the global regularity versus finite- time blowup in solutions of the 3D incompressible Euler equation is a major open problem of modern applied analysis. In this paper, we study a class of one-dimensional models of the axisymmetric hyperbolic boundary blow-up scenario for the 3D Euler equation proposed by Hou and Luo (Multiscale Model Simul 12:1722–1776, 2014) based on extensive numerical simulations. These models generalize the 1D Hou–Luo model suggested in Hou and Luo Luo and Hou (2014), for which finite-time blowup has been established in Choi et al. (arXiv preprint. arXiv:1407.4776, 2014). The main new aspects of this work are twofold. First, we establish finite-time blowup for a model that is a closer approximation of the three-dimensional case than the original Hou–Luo model, in the sense that it contains relevant lower-order terms in the Biot–Savart law that have been discarded in Hou and Luo Choi et al. (2014). Secondly, we show that the blow-up mechanism is quite robust, by considering a broader family of models with the same main term as in the Hou–Luo model. Such blow-up stability result may be useful in further work on understanding the 3D hyperbolic blow-up scenario.