A time discretization scheme of a characteristics method for a fluid–rigid system with discontinuous density

We propose a new characteristics method for the time discretization of a fluid–rigid system in the case when the densities of the fluid and the solid are different. This method is based on a global weak formulation involving only terms defined on the whole fluid–rigid domain. The main idea is to construct a characteristic function which preserves the rigidity of the solid at the discrete time levels. A convergence result for this semi-discrete scheme is then given.     

Modeling and analysis of reactive solute transport in deformable channels with wall adsorption–desorption

We show well posedness for a model of nonlinear reactive transport of chemical in a deformable channel. The channel walls deform due to fluid–structure interaction between an unsteady flow of an incompressible, viscous fluid inside the channel and elastic channel walls. Chemical solutes, which are dissolved in the viscous, incompressible fluid, satisfy a convection–diffusion equation in the bulk fluid, while on the deforming walls, the solutes undergo nonlinear adsorption–desorption physico‐chemical reactions. The problem addresses scenarios that arise, for example, in studies of drug transport in blood vessels. We show the existence of a unique weak solution with solute concentrations that are non‐negative for all times. The analysis of the problem is carried out in the context of semi‐linear parabolic PDEs on moving domains. The arbitrary Lagrangian–Eulerian approach is used to address the domain movement, and the Galerkin method with the Picard–Lindelöf theorem is used to prove existence and uniqueness of approximate solutions. Energy estimates combined with the compactness arguments based on the Aubin–Lions lemma are used to prove convergence of the approximating sequences to the unique weak solution of the problem. It is shown that the solution satisfies the positivity property, that is, that the density of the solute remains non‐negative at all times, as long as the prescribed fluid domain motion is ‘reasonable’. This is the first well‐posedness result for reactive transport problems defined on moving domains of this type. Copyright © 2015 John Wiley & Sons, Ltd.     

An exact Block–Newton algorithm for solving fluid–structure interaction problems

In this Note, we introduce a partitioned Newton based method for solving nonlinear coupled systems arising in the numerical approximation of fluid–structure interaction problems. The originality of this Schur–Newton algorithm lies in the exact Jacobians evaluation involving the fluid–structure linearized subsystems which are here fully developed. To cite this article: M.Á. Fernández, M. Moubachir, C. R. Acad. Sci. Paris, Ser. I 336 (2003).     

Numerical simulation of non-Newtonian blood flow in bypass models

The article presents the numerical investigation of non–Newtonian effects of steady blood flow in complete idealized 3–D bypass models, whose native artery is either coronary or femoral with average physiological parameters. Considering the blood to be a generalized Newtonian fluid, the shear–dependent viscosity is described by two well–known macroscopic non–Newtonian models (the Carreau–Yasuda model and the modified Cross model). The results were obtained by own developed computational software based on the pseudo–compressibility approach and on the cell–centred finite volume method defined on unstructured hexahedral grids. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Implementation of a Preconditioning–Technique in the Hybrid Flow–Solver TAU+

The present paper describes the implementation of a preconditioning method in the hybrid DLR–TAU⁺–code and its application to nearly incompressible flows. The method is designed in order to get an efficient and accurate solution even for very low Mach numbers using a time stepping scheme for the solution of the compressible Navier–Stokes equations. The algorithm is based on the work of Choi and Merkle. The numerical results obtained for inviscid and viscous flows indicate, that for Mach numbers lower than 0.1 the accuracy as well as the convergence properties are almost independent of the fluid speed, like for incompressible codes.     

Solvability,asymptotic analysis and regularity results for a mixed type interaction problem of acoustic waves and piezoelectric structures

In the paper, we investigate the mixed type transmission problem arising in the model of fluid–solid acoustic interaction when a piezoceramic elastic body (Ω⁺) is embedded in an unbounded fluid domain (Ω^?). The corresponding physical process is described by the boundary‐transmission problem for second‐order partial differential equations. In particular, in the bounded domain Ω⁺, we have a 4×4 dimensional matrix strongly elliptic second‐order partial differential equation, while in the unbounded complement domain Ω^?, we have a scalar Helmholtz equation describing acoustic wave propagation. The physical kinematic and dynamic relations mathematically are described by appropriate boundary and transmission conditions. With the help of the potential method and theory of pseudodifferential equations based on the Wiener–Hopf factorization method, the uniqueness and existence theorems are proved in Sobolev–Slobodetskii spaces. We derive asymptotic expansion of solutions, and on the basis of asymptotic analysis, we establish optimal Hölder smoothness results for solutions. Copyright © 2017 John Wiley & Sons, Ltd.     

A Parallel Implementation of the Algebraic Multigrid Method for Solving Problems in Dynamics of Viscous Incompressible Fluid

An algorithm for improving the scalability of the multigrid method used for solving the system of difference equations obtained by the finite volume discretization of the Navier–Stokes equations on unstructured grids with an arbitrary cell topology is proposed. It is based on the cascade assembly of the global level; the cascade procedure gradually decreases the number of processors involved in the computations. Specific features of the proposed approach are described, and the results of solving benchmark problems in the dynamics of viscous incompressible fluid are discussed; the scalability and efficiency of the proposed method are estimated. The advantages of using the global level in the parallel implementation of the multigrid method which sometimes makes it possible to speed up the computations by several fold.     

Nonlinear vibration and instability of embedded double-walled boron nitride nanotubes based on nonlocal cylindrical shell theory

This paper investigates the nonlinear vibration and instability of the embedded double-walled boron nitride nanotubes (DWBNNTs) conveying viscous fluid based on nonlocal piezoelasticity cylindrical shell theory. The elastic medium is simulated as Winkler–Pasternak foundation, and adjacent layers interactions are assumed to have been coupled by van der Walls (vdW) force evaluated based on the Lennard–Jones model. The nonlinear strain terms based on Donnell’s theory are taken into account. The Hamilton’s principle is employed to obtain coupled differential equations, containing displacement and electric potential terms. Differential quadrature method (DQM) is applied to estimate the nonlinear frequency and critical fluid velocity for clamped supported mechanical and free electric potential boundary conditions at both ends of the DWBNNTs. Results indicated that some parameters including nonlocal parameter, elastic medium’s modulus, aspect ratio and vdW force have significant influence on the vibration and instability of the DWBNNT while the fluid viscosity effect is negligible. In addition, the low aspect ratio should be taken into account for DWBNNT in optimum design of nano/micro devices.     

The finite volume method based on the Crouzeix–Raviart element for a fracture model

In this paper, we present and analyze a finite volume method based on the Crouzeix–Raviart element for the coupled fracture model, where the fluid flow is governed by Darcy's law in the one‐dimensional fracture and two‐dimensional surrounding matrix. In the numerical scheme, the pressure in the matrix and fracture is respectively approximated by the Crouzeix–Raviart elements and piecewise constant functions, and then the velocity is calculated by piecewise constant functions element by element. The existence and uniqueness of the numerical solution are discussed, and optimal order error estimates for both the pressure p and the velocity u are proved on general triangulations. We finally carry out numerical experiments, and results confirm our theoretical analysis.     

A linear,decoupled fractional time‐stepping method for the nonlinear fluid–fluid interaction

In this paper, a linear decoupled fractional time stepping method is proposed and developed for the nonlinear fluid–fluid interaction governed by the two Navier–Stokes equations. Partitioned time stepping method is applied to two‐physics problems with stiffness of the coupling terms being treated explicitly and is also unconditionally stable. As for each fluid, the velocity and pressure are respectively determined by just solving one vector‐valued quasi‐elliptic equation and the Possion equation with homogeneous Neumann boundary condition per time step. Therefore, the cost of the fluid–fluid interaction is dominant to solve four simple linear equations, which greatly reduces the computational cost of the whole system. The method exploits properties of the fluid–fluid system to establish its stability and convergence with the same results as the standard scheme. Finally, numerical experiments are presented to show the performance of the proposed method.     

Application of the generalized multipole technique in band structure calculation of two‐dimensional solid/fluid phononic crystals

A multiple monopole method based on the generalized multipole technique is presented for the calculation of band structures of two‐dimensional mixed solid/fluid phononic crystals. In this method, the fields are expanded by using the fundamental solutions with multiple origins. Besides the sources used to expand the wave fields, an extra monopole source is introduced as the external excitation. By varying the frequency of the excitation, the eigenvalues can be localized as the extreme points of an appropriately chosen function. By sweeping the frequency range of interest and sweeping the boundary of the irreducible first Brillouin zone, the band structure of the phononic crystals can be obtained. The method can consider the fluid–solid interface conditions and the transverse wave mode in the solid component strictly. Some typical examples are illustrated to discuss the accuracy of the present method. Copyright © 2014 John Wiley & Sons, Ltd.     

Coupling of mixed finite element and stabilized boundary element methods for a fluid–solid interaction problem in 3D

We introduce and analyze the coupling of a mixed finite element and a boundary element for a three‐dimensional time‐harmonic fluid–solid interaction problem. We consider a formulation in which the Cauchy stress tensor and the rotation are the main variables in the elastic structure and use the usual pressure formulation in the acoustic fluid. The mixed variational formulation in the solid is completed with boundary integral equations relating the Cauchy data of the acoustic problem on the coupling interface. A crucial point in our formulation is the stabilization technique introduced by Hiptmair and coworkers to avoid the well‐known instability issue appearing in the boundary element method treatment of the exterior Helmholtz problem. The main novelty of this formulation, with respect to a previous approach, consists in reducing the computational domain to the solid media and providing a more accurate treatment of the far field effect. We show that the continuous problem is well‐posed and propose a conforming Galerkin method based on the lowest‐order Arnold–Falk–Winther mixed finite element. Finally, we prove that the numerical scheme is convergent with optimal order.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1211–1233, 2014     

Stability analysis of a single-walled carbon nanotube conveying pulsating and viscous fluid with nonlocal effect

For a single-walled carbon nanotube (CNT) conveying fluid, the internal flow is considered to be pulsating and viscous, and the resulting instability and parametric resonance of the CNT are investigated by the method of averaging. The partial differential equation of motion based on the nonlocal elasticity theory is discretized by the Galerkin method and the averaging equations for the first two modes are derived. The stability regions in frequency–amplitude plane are obtained and the influences of nonlocal effect, viscosity and some system parameters on the stability of CNT are discussed in detail. The results show that decrease of nonlocal parameter and increase of viscous parameter both increase the fundamental frequency and critical flow velocity; the dynamic stability of CNT can be enhanced due to nonlocal effect; the contributions of the fluid viscosity on the stability of CNT depend on flow velocity; the axial tensile force can only influence natural frequencies of the system however the viscoelastic characteristic of materials can enhance the dynamic stability of CNT. The conclusions drawn in this paper are thought to be helpful for the vibration analysis and structural design of nanofluidic devices.     

Multi-Fracture Reservoir Simulations for Long Time Physical Processes

Reservoirs with multi-fracture techniques are developed and frequently used for oil and gas industry. Recently, they are also used for deep geothermal reservoirs especially for Hot Dry Rock (HDR). The analysis of the reservoir is generally interested in long time physical properties (10–100 years), e.g. fluid flow, heat transport etc. Typical CFD simulations are limited in this context. Here we developed a fluid flow and heat transport modeling in a multi-fracture reservoir based on the so-called Mixed Dimensional Model (MDM), which describes the different characteristic flows and the heat transport in different dimensions. In the mathematical point of view, these models are discretized based on the Cellular Automaton (CA) method combined with other necessary numerical techniques. The different cases of fluid flow and heat transport in multi-fracture reservoirs have been simulated and shown physical results very reasonably with less computational time. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Divergence-free finite elements for the numerical solution of a hydroelastic vibration problem

In this paper, we analyze a divergence-free finite element method to solve a fluid–structure interaction spectral problem in the three-dimensional case. The unknowns of the resulting formulation are the fluid and solid displacements and the fluid pressure on the interface separating both media. The resulting mixed eigenvalue problem is approximated by using appropriate basis of the divergence-free lowest order Raviart–Thomas elements for the fluid, piecewise linear elements for the solid and piecewise constant elements for the interface pressure. It is proved that eigenvalues and eigenfunctions are correctly approximated and some numerical results are reported in order to assess the performance of the method.     

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.     

A residual-based a posteriori error estimator for a two-dimensional fluid–solid interaction problem

In this paper, we develop an a posteriori error analysis of a mixed finite element method for a fluid–solid interaction problem posed in the plane. The media are governed by the acoustic and elastodynamic equations in time-harmonic regime, respectively, and the transmission conditions are given by the equilibrium of forces and the equality of the normal displacements of the solid and the fluid. The coupling of primal and dual-mixed finite element methods is applied to compute both the pressure of the scattered wave in the linearized fluid and the elastic vibrations that take place in the elastic body. The finite element subspaces consider continuous piecewise linear elements for the pressure and a Lagrange multiplier defined on the interface, and PEERS for the stress and rotation in the solid domain. We derive a reliable and efficient residual-based a posteriori error estimator for this coupled problem. Suitable auxiliary problems, the continuous inf-sup conditions satisfied by the bilinear forms involved, a discrete Helmholtz decomposition, and the local approximation properties of the Clément interpolant and Raviart–Thomas operator are the main tools for proving the reliability of the estimator. Then, Helmholtz decomposition, inverse inequalities, and the localization technique based on triangle-bubble and edge-bubble functions are employed to show the efficiency. Finally, some numerical results confirming the reliability and efficiency of the estimator are reported.     

Laguerre spectral approximation of Stokes’ first problem for third-grade fluid

A Laguerre–Galerkin method is proposed and analyzed for Quasilinear parabolic differential equation which arises from Stokes’ first problem for a third-grade fluid on a semi-infinite interval. By reformulating this equation with suitable functional transforms, it is shown that the Laguerre–Galerkin approximations are convergent on a semi-infinite interval with spectral accuracy. An efficient and accurate algorithm based on the Laguerre–Galerkin approximations to the transformed equations is developed and implemented. Effects of non-Newtonian parameters on the flow phenomena are analyzed and documented.     

Stokes' first problem for a Newtonian fluid in a non‐Darcian porous half‐space using a Laguerre–Galerkin method

A Laguerre–Galerkin method is proposed and analysed for the Stokes' first problem of a Newtonian fluid in a non‐Darcian porous half‐space on a semi‐infinite interval. It is well known that Stokes' first problem has a jump discontinuity on boundary which is the main obstacle in numerical methods. By reformulating this equation with suitable functional transforms, it is shown that the Laguerre–Galerkin approximations are convergent on a semi‐infinite interval with spectral accuracy. An efficient and accurate algorithm based on the Laguerre–Galerkin approximations of the transformed equations is developed and implemented. Numerical results indicating the high accuracy and effectiveness of this algorithm are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

Dynamics of Strongly Nonlinear Kinks and Solitons in a Two‐Layer Fluid

Perturbation theory is developed for interaction of strongly nonlinear solitary waves close to the limiting, tabletop solitons (Π‐solitons). The method is based on representing each soliton as a compound of two kinks so that the interaction of N solitons is treated as the interaction of 2N kinks. As an example the Miyata–Choi–Camassa equations for a two‐layer fluid is considered. Equations for kink coordinates are obtained and analyzed. Some nontrivial features of two‐soliton interaction characteristic of the strongly nonlinear case are established.