Weakly coupled analysis of a blade in multiphase mixing vessel

Two or more physical systems frequently interact with each other, where the independent solution of one system is impossible without a simultaneous solution of the others. An obvious coupled system is that of a dynamic fluid-structure interaction. [8] In this paper a computational analysis of the fluid-structure interaction in a mixing vessel is presented. In mixing vessels the fluid can have a significant influence on the deformation of blades during mixing, depending on speed of mixing blades and fluid viscosity. For this purpose a computational weakly coupled analysis has been performed to determine the multiphase fluid influences on the mixing vessel structure. The multiphase fluid field in the mixing vessel was first analyzed with the computational fluid dynamics (CFD) code CFX. The results in the form of pressure were then applied to the blade model, which was the analysed with the structural code MSC.visualNastran forWindows, which is based on the finite element method (FEM). (© 2004 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

基于GMRES算法的弹性结构强耦合小变形流激振动分析

使用混合广义变分原理,将基于Lagrange表述的小位移变形结构振动问题与基于Euler描述的不可压缩粘性流动问题,统一在功率平衡的框架下建立流固系统的耦合控制方程．用有限元格式做空间离散后,再用广义梯形法将有限元控制方程转化为增量型的线性方程组,该方程组的系数矩阵具有非对称性,其中元素含对流效应和时间因子．将GMRES算法与振动分析的Newmark法和流动分析的Hughes预测多修正法结合,发展成一种基于GMRES-Hughes-Newmark的稳定算法,用于计算具有复杂几何边界的强耦合流激振动问题．以混流式水轮机叶道为数值算例的计算表明,模拟结果与试验实测结果吻合较好．     

Construction of a high order fluid-structure interaction solver

Accuracy is critical if we are to trust simulation predictions. In settings such as fluid-structure interaction, it is all the more important to obtain reliable results to understand, for example, the impact of pathologies on blood flows in the cardiovascular system. In this paper, we propose a computational strategy for simulating fluid structure interaction using high order methods in space and time.First, we present the mathematical and computational core framework, Life, underlying our multi-physics solvers. Life is a versatile library allowing for 1D, 2D and 3D partial differential solves using h/p type Galerkin methods. Then, we briefly describe the handling of high order geometry and the structure solver. Next we outline the high-order space-time approximation of the incompressible Navier-Stokes equations and comment on the algebraic system and the preconditioning strategy. Finally, we present the high-order Arbitrary Lagrangian Eulerian (ALE) framework in which we solve the fluid-structure interaction problem as well as some initial results.     

Modeling concept and numerical simulation of ultrasonic wave propagation in a moving fluid-structure domain based on a monolithic approach

In the present study, we propose a novel multiphysics model that merges two time-dependent problems – the Fluid-Structure Interaction (FSI) and the ultrasonic wave propagation in a fluid-structure domain with a one directional coupling from the FSI problem to the ultrasonic wave propagation problem. This model is referred to as the “eXtended fluid-structure interaction (eXFSI)” problem. This model comprises isothermal, incompressible Navier–Stokes equations with nonlinear elastodynamics using the Saint-Venant Kirchhoff solid model. The ultrasonic wave propagation problem comprises monolithically coupled acoustic and elastic wave equations. To ensure that the fluid and structure domains are conforming, we use the ALE technique. The solution principle for the coupled problem is to first solve the FSI problem and then to solve the wave propagation problem. Accordingly, the boundary conditions for the wave propagation problem are automatically adopted from the FSI problem at each time step. The overall problem is highly nonlinear, which is tackled via a Newton-like method. The model is verified using several alternative domain configurations. To ensure the credibility of the modeling approach, the numerical solution is contrasted against experimental data.     

Numerical modeling on inlet aperture effects on flow pattern in primary settling tanks

Inlets should be designed to dissipate the kinetic energy or velocity head of the mixed liquor and to prevent short-circuiting, mitigate the effects of density currents, and minimize blanket disturbances. Flow in primary settling tank is simulated by means of computational fluid dynamics. The fluid is assumed incompressible and non-buoyant. A two-dimensional computational and one phase fluid dynamics model was built to simulate the flow properties in the settling tank including the velocity profiles, the flow separation area and kinetic energy. In this study, the RNG turbulent model was solved with the Navier–Stokes equations. In order to evaluate hydraulic influences on the velocity profile, separation length and kinetic energy, three different of opening positions and two and three aperture in inlets were simulated. The flow model uses to apply a fixed-grid of cells that are all rectangular faces; the fluid moves through the grid and free surfaces are tracked with the volume-of-fluid (VOF) technique. Effects of numbers and locations of inlet apertures on the flow field are presented and the results show the positions of inlet apertures are affected on the flow pattern in the settling basin and increasing the numbers of slots can reduce kinetic energy in the inlet zone and produce uniform flow.     

Pure Lagrangian and semi-Lagrangian finite element methods for the numerical solution of Navier–Stokes equations

In this paper we propose a unified formulation to introduce Lagrangian and semi-Lagrangian velocity and displacement methods for solving the Navier–Stokes equations. This formulation allows us to state classical and new numerical methods. Several examples are given. We combine them with finite element methods for spatial discretization. In particular, we propose two new second-order characteristics methods in terms of the displacement, one semi-Lagrangian and the other one pure Lagrangian. The pure Lagrangian displacement methods are useful for solving free surface problems and fluid-structure interaction problems because the computational domain is independent of the time and fluid–solid coupling at the interphase is straightforward. However, for moderate to high-Reynolds number flows, they can lead to high distortion in the mesh elements. When this happens it is necessary to remesh and reinitialize the transformation to the identity. In order to assess the performance of the obtained numerical methods, we solve different problems in two space dimensions. In particular, numerical results for a sloshing problem in a rectangular tank and the flow in a driven cavity are presented.     

Geometric multigrid for an implicit-time immersed boundary method

The immersed boundary (IB) method is an approach to fluid-structure interaction that uses Lagrangian variables to describe the deformations and resulting forces of the structure and Eulerian variables to describe the motion and forces of the fluid. Explicit time stepping schemes for the IB method require solvers only for Eulerian equations, for which fast Cartesian grid solution methods are available. Such methods are relatively straightforward to develop and are widely used in practice but often require very small time steps to maintain stability. Implicit-time IB methods permit the stable use of large time steps, but efficient implementations of such methods require significantly more complex solvers that effectively treat both Lagrangian and Eulerian variables simultaneously. Several different approaches to solving the coupled Lagrangian-Eulerian equations have been proposed, but a complete understanding of this problem is still emerging. This paper presents a geometric multigrid method for an implicit-time discretization of the IB equations. This multigrid scheme uses a generalization of box relaxation that is shown to handle problems in which the physical stiffness of the structure is very large. Numerical examples are provided to illustrate the effectiveness and efficiency of the algorithms described herein. These tests show that using multigrid as a preconditioner for a Krylov method yields improvements in both robustness and efficiency as compared to using multigrid as a solver. They also demonstrate that with a time step 100–1000 times larger than that permitted by an explicit IB method, the multigrid-preconditioned implicit IB method is approximately 50–200 times more efficient than the explicit method.     

Numerical analysis of two partitioned methods for uncoupling evolutionary MHD flows

Magnetohydrodynamics (MHD) studies the dynamics of electrically conducting fluids, involving Navier–Stokes (NSE) equations in fluid dynamics and Maxwell equations in eletromagnetism. The physical processes of fluid flows and electricity and magnetism are quite different and numerical simulations of each subprocess can require different meshes, time steps, and methods. In most terrestrial applications, MHD flows occur at low‐magnetic Reynold numbers. We introduce two partitioned methods to solve evolutionary MHD equations in such cases. The methods we study allow us at each time step to call NSE and Maxwell codes separately, each possibly optimized for the subproblem's respective physics. Complete error analysis and computational tests supporting the theory are given.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1083–1102, 2014     

Advanced FEM for free surface flow with application to landslides

An advanced space-time finite element method is presented to investigate movements of landslides and their interaction with flexible structures. The mechanics of liquefied soil is described by Navier-Stokes-equations for visco-plastic non-newtonian fluid. Likewise the fluid the kinematics of the structure is described by velocities, taking large rotations into account. This leads to a monolithic fluid-structure interaction approach considering the multi-field problem as a whole. The discretized model equations are assembled in a single set of algebraic equations, which are solved by applying Newton-Raphson scheme. Free surface motion of landslide is described by the level-set method. To reduce computational effort the fragmented finite element method is used, where only active finite elements are evaluated. A pde-based extrapolation of the velocity-field is applied to ensure an accurate transport of distance function, which defines the profile in space and time of the free surfaces. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems

The governing dynamics of fluid flow is stated as a system of partial differential equations referred to as the Navier-Stokes system. In industrial and scientific applications, fluid flow control becomes an optimization problem where the governing partial differential equations of the fluid flow are stated as constraints. When discretized, the optimal control of the Navier-Stokes equations leads to large sparse saddle point systems in two levels. In this paper, we consider distributed optimal control for the Stokes system and test the particular case when the arising linear system can be compressed after eliminating the control function. In that case, a system arises in a form which enables the application of an efficient block matrix preconditioner that previously has been applied to solve complex-valued systems in real arithmetic. Under certain conditions, the condition number of the so preconditioned matrix is bounded by 2. The numerical and computational efficiency of the method in terms of number of iterations and execution time is favorably compared with other published methods.     

Optimal boundary control with critical penalization for a PDE model of fluid–solid interactions

We study the finite-horizon optimal control problem with quadratic functionals for an established fluid-structure interaction model. The coupled PDE system under investigation comprises a parabolic (the fluid) and a hyperbolic (the solid) dynamics; the coupling occurs at the interface between the regions occupied by the fluid and the solid. We establish several trace regularity results for the fluid component of the system, which are then applied to show well-posedness of the Differential Riccati Equations arising in the optimization problem. This yields the feedback synthesis of the unique optimal control, under a very weak constraint on the observation operator; in particular, the present analysis allows general functionals, such as the integral of the natural energy of the physical system. Furthermore, this work confirms that the theory developed in Acquistapace et al. (Adv Diff Eq, [2])—crucially utilized here—encompasses widely differing PDE problems, from thermoelastic systems to models of acoustic-structure and, now, fluid-structure interactions.     

Pressure Boundary Conditions for Blood Flows

Simulations of blood flows in arteries require numerical solutions of fluid-structure interactions involving Navier-Stokes equations coupled with large displacement visco-elasticity for the vessels. Among the various simplifications which have been proposed, the surface pressure model leads to a hierarchy of simpler models including one that involves only the pressure. The model exhibits fundamental frequencies which can be computed and compared with the pulse. Yet unconditionally stable time discretizations can be constructed by combining implicit time schemes with Galerkin-characteristic discretization of the convection terms in the Navier-Stokes equations. Such problems with prescribed pressure on the walls will be shown to be efficient and accurate as an approximation of the full fluid structure interaction problem.     

Applications of level set methods in computational biophysics

We describe in this paper two applications of Eulerian level set methods to fluid-structure problems arising in biophysics. The first one is concerned with three-dimensional equilibrium shapes of phospholipidic vesicles. This is a complex problem, which can be recast as the minimization of the curvature energy of an immersed elastic membrane, under a constant area constraint. The second deals with isolated cardiomyocyte contraction. This problem corresponds to a generic incompressible fluid-structure coupling between an elastic body and a fluid. By the choice of these two quite different situations, we aim to bring evidence that Eulerian methods provide efficient and flexible computational tools in biophysics applications.     

Convergence of a Lagrangian scheme for a compressible Naviers–Stokes model defined on a domain depending on time

This paper deals with the numerical resolution of the Navier–Stokes equations defined on a time dependent domain. We give an existence result for a fluid-structure interaction problem in which the boundary is governed by a thin plate operator. We propose to solve the fluid equations with the characteristics method. We approach the total derivative with a “regularized” finite difference scheme and we study the convergence of the discrete problem towards the continuous one.     

The Strong Stability and Instability of a Fluid-Structure Semigroup

The strong stability problem for a fluid-structure interactive partial differential equation (PDE) is considered. The PDE comprises a coupling of the linearized Stokes equations to the classical system of elasticity, with the coupling occurring on the boundary interface between the fluid and solid media. Because of the nature of the unbounded coupling between fluid and structure, the resolvent of the associated semigroup generator will not be a compact operator. In consequence, the classical solution to the stability problem, by means of the Nagy-Foias decomposition, will not avail here. Moreover, it is not practicable to write down explicitly the resolvent of the fluid-structure generator; this situation thus makes it problematic to use the well-known semigroup stability result of Arendt-Batty and Lyubich-Phong. When a locally supported boundary dissipative mechanism is in place, we derive here a result of strong decay for this fluid-structure PDE. In the absence of said dissipative mechanism, we show the lack of asymptotic decay for solutions corresponding to arbitary initial data of finite energy.     

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.     

Variational Formulations for Viscous Flow

For physical systems, the dynamics of which is formulated within the framework of Lagrange formalism the dynamics is completely defined by only one function, namely the Lagrangian. As well-known the whole conservative Newtonian mechanics has been successfully embedded into this methodical concept. Different from this, in continuum theories many open questions remain up to date, especially when considering dissipative processes. The viscous flow of a fluid, given by the Navier-Stokes equations is a typical example for this. In this contribution a special approach for finding a Lagrangian for viscous flow is suggested and discussed. The equations of motion resulting from the respective Lagrangian are compared to the Navier-Stokes equations and differences are discussed. For a simple flow example their solution is compared to the one resulting from Navier-Stokes equations. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Fluid dynamics and thermodynamics as a unified field theory

We study the problem of consistency of equations of continuum dynamics (using the Euler equations and the continuity equation as examples) and thermodynamic equations of state (for the specific free energy, entropy, and volume). We propose a variant of the Hamiltonian formulation of a model that combines the fluid dynamics of a potential flow of a compressible fluid or gas and local equilibrium thermodynamics into a unified field theory. Thermodynamic equations of state appear in this model as second-class constraint equations. As a consistency condition, there arises another second-class constraint requiring that the product of density and temperature should be independent of time. The model provides an in-principle possibility of finding the time dependence of the specific entropy of the arising dynamical system.     

Multiscale Considerations for Sound Generation in Low Mach Number Flows

Classical approaches in aeroacoustics are mainly based on analytical solutions of the linear wave equations which are valid in the far field. The sound generation is approximated by source terms obtained from a flow simulation. This procedure designated as the ‘acoustic analogy’ was initiated in the classical work of Lighthill. In order to tackle the sound generation problem at low Mach numbers, we consider a multiple scale asymptotic analysis. As we deal with a fluid flow generating the sound itself, the asymptotic ansatz uses one time scale given by the flow convection, but two space scales due to the difference in fluid and sound velocity. The insight given by this analysis is used to obtain source terms describing the sound generation and perturbation equations for the sound propagation. Numerical results are shown for the example of a co‐rotating vortex pair.     

Partitioned Algorithms for Fluid-Structure Interaction Problems in Haemodynamics

We consider the fluid-structure interaction problem arising in haemodynamic applications. The finite elasticity equations for the vessel are written in Lagrangian form, while the Navier-Stokes equations for the blood in Arbitrary Lagrangian Eulerian form. The resulting three fields problem (fluid/ structure/ fluid domain) is formalized via the introduction of three Lagrange multipliers and consistently discretized by p-th order backward differentiation formulae (BDFp). We focus on partitioned algorithms for its numerical solution, which consist in the successive solution of the three subproblems. We review several strategies that all rely on the exchange of Robin interface conditions and review their performances reported recently in the literature. We also analyze the stability of explicit partitioned procedures and convergence of iterative implicit partitioned procedures on a simple linear FSI problem for a general BDFp temporal discretizations.