Domain Decomposition Methods for Sensitivity Analysis of a Nonlinear Aeroelasticity Problem

Abstract

We consider the nonlinear aeroelasticity problem of the interaction between a viscous, incompressible fluid and Lin elastic solid undergoing large displacement. The non-linearities of the problem formulation include the solid and fluid governing equations. as well as thc dependence of the How geometry on the solid deformation. The resulting coupling is thus two-way. We develop domain-decomposition methods for solution and sensitivity analysis of the coupled problem. The domain decomposition is in the form of a block-Gauss-Seidel-like prcconditioncr that decomposes ihc coupled-domain problem into distinct nonovcrlapping fluid and solid subdotnain problems. The preconditioner thus enables exploitation or single-domain algorithms for solid and fluid mechanics discretization and solution. On the other hand, two-way fluid-solid coupling is retained within the residuals, which is essential for correct sensitivities. Sensitivities of field quantities can be found with little additional work beyond that required for solving the coupled fluid-solid system. The methodology developed here is illustrated by the solution of a problem of viscous incompressible flow about an infinite clastic cylinder. Sensitivities of the resulting velocity and displacement fields with respect to elastic modulus and fluid viscosity are computed.     

Finite element computations for unsteady fluid and elastic membrane interaction problems

The interaction between the hydrodynamic forces of a flow field and the elastic forces of adjacent deformable boundaries is described by elastohydrodynamics, a coupled fluid–elastic membrane problem. Direct numerical solution of the unsteady, highly non-linear equations requires that the dynamic evolution of both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. This paper describes a numerical algorithm based on the deformable spatial domain space–time (DSD/ST) finite element method for the unsteady motion of an incompressible, viscous fluid with elastic membrane interaction. The unsteady Navier–Stoke and elastic membrane equations are solved separately using an iterative procedure by the GMRES technique with an incomplete lower-upper (ILU) decomposition at every time instant. One-dimensional, two-dimensional and deformable domain model problems are used to demonstrate the capabilities and accuracy of the present algorithm. Both steady state and transient problems are studied. © 1997 John Wiley & Sons, Ltd.     

Fluid–solid interaction problems with thermal convection using the immersed element‐free Galerkin method

In this work, the immersed element‐free Galerkin method (IEFGM) is proposed for the solution of fluid–structure interaction (FSI) problems. In this technique, the FSI is represented as a volumetric force in the momentum equations. In IEFGM, a Lagrangian solid domain moves on top of an Eulerian fluid domain that spans over the entire computational region. The fluid domain is modeled using the finite element method and the solid domain is modeled using the element‐free Galerkin method. The continuity between the solid and fluid domains is satisfied by means of a local approximation, in the vicinity of the solid domain, of the velocity field and the FSI force. Such an approximation is achieved using the moving least‐squares technique. The method was applied to simulate the motion of a deformable disk moving in a viscous fluid due to the action of the gravitational force and the thermal convection of the fluid. An analysis of the main factors affecting the shape and trajectory of the solid body is presented. The method shows a distinct advantage for simulating FSI problems with highly deformable solids. Copyright © 2009 John Wiley & Sons, Ltd.     

基于特征分裂有限元准隐格式的共轭传热整体耦合数值模拟方法

共轭传热现象在科学和工程领域中大量存在.随着计算能力的发展,对共轭传热现象进行准确有效的数值模拟,成为科学研究和工程设计上的重要挑战.共轭传热数值模拟的方法可以分为两大类:分区耦合和整体耦合.本文采用有限元法对共轭传热问题进行整体耦合模拟.固体传热求解采用标准的伽辽金有限元方法.流动求解采用基于特征分裂的有限元方法.该...     

Reduced order modelling method via proper orthogonal decomposition (POD) for flow around an oscillating cylinder

This paper presents reduced order modelling (ROM) in fluid–structure interaction (FSI). The ROM via the proper orthogonal decomposition (POD) method has been chosen, due to its efficiency in the domain of fluid mechanics. POD-ROM is based on a low-order dynamical system obtained by projecting the nonlinear Navier–Stokes equations on a smaller number of POD modes. These POD modes are spatial and temporally independent. In FSI, the fluid and structure domains are moving, owing to which the POD method cannot be applied directly to reduce the equations of each domain. This article proposes to compute the POD modes for a global velocity field (fluid and solid), and then to construct a low-order dynamical system. The structure domain can be decomposed as a rigid domain, with a finite number of degrees of freedom. This low-order dynamical system is obtained by using a multiphase method similar to the fictitious domain method. This multiphase method extends the Navier–Stokes equations to the solid domain by using a penalisation method and a Lagrangian multiplier. By projecting these equations on the POD modes obtained for the global velocity field, a nonlinear low-order dynamical system is obtained and tested on a case of high Reynolds number.     

A new approach for computing the steady state fluid–structure interaction response of periodic problems

A special type of fluid–structure interaction (FSI) problems are problems with periodic boundary conditions like in turbomachinery. The steady state FSI response of these problems is usually calculated with similar techniques as used for transient FSI analyses. This means that, when the fluid and structure problem are not simultaneously solved with a monolithic approach, the problem is partitioned into a fluid and structural part and that each time step coupling iterations are performed to account for strong interactions between the two sub-domains. This paper shows that a time-partitioned FSI computation can be very inefficient to compute the steady state FSI response of periodic problems. A new approach is introduced in which coupling iterations are performed on periodic level instead of per time step. The convergence behaviour can be significantly improved by implementing existing partitioned solution methods as used for time step coupling (TSC) algorithms in the time periodic coupling (TPC) framework. The new algorithm has been evaluated by comparing the convergence behaviour to TSC algorithms. It is shown that the number of fluid–structure evaluations can be considerably reduced when a TPC algorithm is applied instead of a TSC. One of the most appealing advantages of the TPC approach is that the structural problem can be solved in the frequency domain resulting in a very efficient algorithm for computing steady state FSI responses.     

A finite-element approach in stream-tube method for solving fluid and solid mechanics problems

This paper presents a theoretical formulation in which the stream-tube method (STM) is examined through a variational approach for solving solid strain and fluid flow problems with finite elements. The analysis considers a reference domain, used as computational domain, related to the physical domain by an unknown transformation function to be determined numerically. Mass conservation is automatically verified by STM. The variational approach leads to eliminate the pressure in fluid problems and avoids to set up a mixed displacement–pressure procedure in the case of incompressible solids. Examples are given for fluid flows, applications and comparisons are also provided in the bending problem in elasticity.     

Nonlinear acceleration of coupled fluid–structure transient thermal problems by Anderson mixing

Conjugate heat‐transfer problems are typically solved using partitioned methods where fluid and solid subdomains are evaluated separately by dedicated solvers coupled through a common boundary. Strongly coupled schemes for transient analysis require fluid and solid problems to be solved many times each time step until convergence to a steady state. In many practical situations, a fairly simple and frequently employed fixed‐point iteration process is rather ineffective; it leads to a large number of iterations per time step and consequently to long simulation times. In this article, Anderson mixing is proposed as a fixed‐point convergence acceleration technique to reduce computational cost of thermal coupled fluid–solid problems. A number of other recently published methods with applications to similar fluid–structure interaction problems are also reviewed and analyzed. Numerical experiments are presented to illustrate relative performance of these methods on a test problem of rotating pre‐swirl cavity air flow interacting with a turbine disk. It is observed that performance of Anderson mixing method is superior to that of other algorithms in terms of total iteration counts. Additional computational savings are demonstrated by reusing information from previously solved time steps. Copyright © All rights reserved 2012 Rolls‐Royce plc.     

An Initial and Boundary Value Problem Modeling of Fish-like Swimming

In this paper we consider an initial and boundary value problem that models the self-propelled motion of solids in a bidimensional viscous incompressible fluid. The self-propelling mechanism, consisting of appropriate deformations of the solids, is a simplified model of the propulsion mechanism of fish-like swimmers. The governing equations consist of the Navier–Stokes equations for the fluid, coupled to Newton’s laws for the solids. Since we consider the case in which the fluid–solid system fills a bounded domain we have to tackle a free boundary value problem. The main theoretical result in the paper asserts the global existence and uniqueness (up to possible contacts) of strong solutions of this problem. The second novel result of this work is the provision of a numerical method for the fluid–solid system. This method provides a simulation of the simultaneous displacement of several swimmers and is tested on several examples.     

The effect of local thermal non-equilibrium on conduction in porous channels with a uniform heat source

We examine the effect of local thermal non-equilibrium on the steady state heat conduction in a porous layer in the presence of internal heat generation. A uniform source of heat is present in either the fluid or the solid phase. A two-temperature model is assumed and analytical solutions are presented for the resulting steady-state temperature profiles in a uniform porous slab. Attention is then focussed on deriving simple conditions which guarantee local thermal equilibrium.     

Asymptotic analysis for micropolar fluidsAnalyse asymptotique des fluides micropolaires

The flow of a micropolar fluid through a wavy constricted channel which depends on a small parameter ε?1 is considered. The asymptotic solution is built and justified thanks to a study of the boundary layers terms. The Stokes and Navier–Stokes problems set in a tube structure were previously considered. The method of partial asymptotic decomposition of domain (MAPPD) is also applied and justified for the micropolar flow problem. This method reduces the initial problem to the problem set in the boundary layers domain. To cite this article: D. Dupuy et al., C. R. Mecanique 332 (2004).     

A modular,operator‐splitting scheme for fluid–structure interaction problems with thick structures

We present an operator‐splitting scheme for fluid–structure interaction (FSI) problems in hemodynamics, where the thickness of the structural wall is comparable to the radius of the cylindrical fluid domain. The equations of linear elasticity are used to model the structure, while the Navier–Stokes equations for an incompressible viscous fluid are used to model the fluid. The operator‐splitting scheme, based on the Lie splitting, separates the elastodynamics structure problem from a fluid problem in which structure inertia is included to achieve unconditional stability. We prove energy estimates associated with unconditional stability of this modular scheme for the full nonlinear FSI problem defined on a moving domain, without requiring any sub‐iterations within time steps. Two numerical examples are presented, showing excellent agreement with the results of monolithic schemes. First‐order convergence in time is shown numerically. Modularity, unconditional stability without temporal sub‐iterations, and simple implementation are the features that make this operator‐splitting scheme particularly appealing for multi‐physics problems involving FSI. Copyright © 2013 John Wiley & Sons, Ltd.     

A quasi‐dynamic procedure for coupled thermal simulations

This paper outlines the development and adaptation of a coupling strategy for transient temperature analysis in a solid via a conjugate heat transfer method. This study proposes a quasi‐dynamic coupling procedure to bridge the temporal disparities between the fluid and the solid. In this approach, dynamic thermal modeling in the solid is coupled with a sequence of steady states in the fluid. This quasi‐dynamic algorithm has been applied to the problem of convective heat transfer over, and transient conduction heat transfer within, a flat plate using the severe thermal conditions of a solid propellant rocket. Two different coupled thermal computations have been performed. In the first one—referred to as the reference computation—the coupling period is equal to the smallest solid time constant. In the second one, a very large coupling period is used. The results show that the procedure can predict accurate transient temperature fields at a reasonable computational cost. The simulation CPU time is approximately reduced by up to 90%, while maintaining a very good accuracy. All the details of the numerical test case are given in the paper. This application illustrates the capabilities and the overall efficiency of this coupled approach in a solid transient problem using long term simulations of time dependent flows. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical approximation of the heat transfer between domains separated by thin walls

In this paper, we analyse the numerical approximation of the heat transfer problem between two subdomains that we will consider filled with a fluid and separated by a thin solid wall. First of all, we state the problem in the whole domain with discontinuous physical properties. As an alternative and under certain assumptions on the separating walls, a classical Robin boundary condition between the fluid domains is obtained, thus eliminating the solid wall, and according to which the heat flux is proportional to the temperature difference between the two subdomains. Apart from discussing the relation between both approaches, we consider their numerical approximation, considering different alternatives for the first case, that is, the case in which temperatures are also computed in the solid wall. Copyright © 2006 John Wiley & Sons, Ltd.     

Exact solutions of linearized equations of convection of a weakly compressible fluid

A mathematical model of fluid convection under microgravity conditions is considered. The equation of state is used in a form that allows considering the fluid as a weakly compressible medium. Based on the previously proposed mathematical model of convection of a weakly compressible fluid, unsteady convective motion in a vertical band, with a heat flux periodic in time set on the solid boundaries of this band, is considered. This model of convection allows one to study the problem with the boundary thermal model oscillating in an antiphase rather than in-phase mode, while the latter was required for the model of microconvection of an isothermally incompressible fluid. Exact solutions for velocity components and temperature are derived, and the trajectories of fluid particles are constructed. For comparison, the trajectories predicted by the classical Oberbeck-Boussinesq model of convection and by the microconvection model are presented.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 2, pp. 52–63, March–April, 2005.     

Nonlinear Stability of Convection in a Porous Layer with Solid Partitions

We show that for many classes of convection problem involving a porous layer, or layers, interleaved with finite but non-deformable solid layers, the global nonlinear stability threshold is exactly the same as the linear instability one. The layer(s) of porous material may be of Darcy type, Brinkman type, possess an anisotropic permeability, or even be such that they are of local thermal non-equilibrium type where the fluid and solid matrix constituting the porous material may have different temperatures. The key to the global stability result lies in proving the linear operator attached to the convection problem is a symmetric operator while the nonlinear terms must satisfy appropriate conditions.     

Onset of convection in a horizontal porous channel with uniform heat generation using a thermal nonequilibrium model

This paper considers the onset of free convection in a horizontal fluid-saturated porous layer with uniform heat generation. Attention is focused on cases where the fluid and solid phases are not in local thermal equilibrium, and where two energy equations describe the evolution of the temperature of each phase. Standard linearized stability theory is used to determine how the criterion for the onset of convection varies with the inter-phase heat transfer coefficient, H, and the porosity-modified thermal conductivity ratio, γ. We also present asymptotic solutions for small values of H. Excellent agreement is obtained between the asymptotic and the numerical results.     

固体颗粒在热对流下沉降的直接数值模拟研究

在任意拉格朗日-欧拉算法模拟等温惰性颗粒两相流的基础上,增加对能量方程的联立求解,在热对流条件下对固体颗粒在不同雷诺数下的沉降规律进行了直接数值模拟.结果表明:热对流引起了流场流动的变化和不对称,颗粒在热流体中沉降,热对流产生的力加速了冷颗粒的运动,尾部形成了涡脱落;颗粒在冷流体中沉降,热对流产生的力阻碍了冷颗粒的运动,尾部形成了羽流. 随着雷诺数的增大,颗粒经历了稳定沉降、周期性摆动到不稳定、无规则运动3个阶段.