A stable Galerkin reduced order model for coupled fluid–structure interaction problems

A stable reduced order model (ROM) of a linear fluid–structure interaction (FSI) problem involving linearized compressible inviscid flow over a flat linear von Kármán plate is developed. Separate stable ROMs for each of the fluid and the structure equations are derived. Both ROMs are built using the ‘continuous’ Galerkin projection approach, in which the continuous governing equations are projected onto the reduced basis modes in a continuous inner product. The mode shapes for the structure ROM are the eigenmodes of the governing (linear) plate equation. The fluid ROM basis is constructed via the proper orthogonal decomposition. For the linearized compressible Euler fluid equations, a symmetry transformation is required to obtain a stable formulation of the Galerkin projection step in the model reduction procedure. Stability of the Galerkin projection of the structure model in the standard L² inner product is shown. The fluid and structure ROMs are coupled through solid wall boundary conditions at the interface (plate) boundary. An a priori energy linear stability analysis of the coupled fluid/structure system is performed. It is shown that, under some physical assumptions about the flow field, the FSI ROM is linearly stable a priori if a stabilization term is added to the fluid pressure loading on the plate. The stability of the coupled ROM is studied in the context of a test problem of inviscid, supersonic flow past a thin, square, elastic rectangular panel that will undergo flutter once the non‐dimensional pressure parameter exceeds a certain threshold. This a posteriori stability analysis reveals that the FSI ROM can be numerically stable even without the addition of the aforementioned stabilization term. Moreover, the ROM constructed for this problem properly predicts the maintenance of stability below the flutter boundary and gives a reasonable prediction for the instability growth rate above the flutter boundary. Copyright © 2013 John Wiley & Sons, Ltd.     

On the computation of steady‐state compressible flows using a discontinuous Galerkin method

Computation of compressible steady‐state flows using a high‐order discontinuous Galerkin finite element method is presented in this paper. An accurate representation of the boundary normals based on the definition of the geometries is used for imposing solid wall boundary conditions for curved geometries. Particular attention is given to the impact and importance of slope limiters on the solution accuracy for flows with strong discontinuities. A physics‐based shock detector is introduced to effectively make a distinction between a smooth extremum and a shock wave. A recently developed, fast, low‐storage p‐multigrid method is used for solving the governing compressible Euler equations to obtain steady‐state solutions. The method is applied to compute a variety of compressible flow problems on unstructured grids. Numerical experiments for a wide range of flow conditions in both 2D and 3D configurations are presented to demonstrate the accuracy of the developed discontinuous Galerkin method for computing compressible steady‐state flows. Copyright © 2007 John Wiley & Sons, Ltd.     

Some anomalies of the numerical solutions of the Euler equations

Summary The numerical solution of the unsteady Euler equations for compressible flow over a circular cylinder is obtained using standard numerical techniques. The equations, written in cylindrical coordinates, are discretized on an orthogonal grid via central differences for spatial derivatives, using a simple second order artificial viscosity form and a special treatment of the boundary conditions. Backward differences in time are employed resulting in a large system of nonlinear difference equations at each step. A direct solver (LAPACK), based on an efficient Gaussian elimination procedure for banded matrices, is used to solve the linearized system of equations. The stability of the nonunique solutions of the steady Euler equations is investigated. It is demonstrated that the symmetric solutions, with zero circulation, are not stable. For a certain Mach number range, a periodic solution is obtained where the shock oscillation persists. If a periodic circulation (within a certain frequency range) is enforced in the far field, an irregular solution emerges with unpredictable shock motions. For such a solution, the Lyapunov exponent is shown to be greater than zero, indicating the appearance of chaos.     

Discontinuous Galerkin methods for non‐linear elasticity

This paper presents the formulation and a partial analysis of a class of discontinuous Galerkin methods for quasistatic non‐linear elasticity problems. These methods are endowed with several salient features. The equations that define the numerical scheme are the Euler–Lagrange equations of a one‐field variational principle, a trait that provides an elegant and simple derivation of the method. In consonance with general discontinuous Galerkin formulations, it is possible within this framework to choose different numerical fluxes. Numerical evidence suggests the absence of locking at near‐incompressible conditions in the finite deformations regime when piecewise linear elements are adopted. Finally, a conceivable surprising characteristic is that, as demonstrated with numerical examples, these methods provide a given accuracy level for a comparable, and often lower, computational cost than conforming formulations. Stabilization is occasionally needed for discontinuous Galerkin methods in linear elliptic problems. In this paper we propose a sufficient condition for the stability of each linearized non‐linear elastic problem that naturally includes material and geometric parameters; the latter needed to account for buckling. We then prove that when a similar condition is satisfied by the discrete problem, the method provides stable linearized deformed configurations upon the addition of a standard stabilization term. We conclude by discussing the complexity of the implementation, and propose a computationally efficient approach that avoids looping over both elements and element faces. Several numerical examples are then presented in two and three dimensions that illustrate the performance of a selected discontinuous Galerkin method within the class. Copyright © 2006 John Wiley & Sons, Ltd.     

基于特征正交分解的三维壁板非线性颤振分析

基于von Karman 大变形理论及活塞理论建立超音速流中壁板的气动弹性方程。采用特征正交分解法 (POD)结合向伽辽金法(Galerkin)的映射这样一种半解析法建立降阶模型(ROM)求解三维壁板的非线性气动弹性问题,并与传统的Galerkin法对比。发现并证明了POD数值模态与伽辽金法简谐基函数之间转换矩阵的正交性,从而简化了POD降阶模型的建立过程。通过数值算例考察了POD法的准确性、收敛性及高效性。结果表明POD降阶模型能够以更少的模态,更高的计算效率达到与Galerkin法同样的精度。以长宽比4为例,POD法以2个模态,3s的时间计算了壁板的振动响应;而Galerkin法需要16个模态,900s的时间。     

Global and local POD models for the prediction of compressible flows with DG methods

Proper orthogonal decomposition (POD) allows to compress information by identifying the most energetic modes obtained from a database of snapshots. In this work, POD is used to predict the behavior of compressible flows by means of global and local approaches, which exploit some features of a discontinuous Galerkin spatial discretization. The presented global approach requires the definition of high‐order and low‐order POD bases, which are built from a database of high‐fidelity simulations. Predictions are obtained by performing a cheap low‐order simulation whose solution is projected on the low‐order basis. The projection coefficients are then used for the reconstruction with the high‐order basis. However, the nonlinear behavior related to the advection term of the governing equations makes the use of global POD bases quite problematic. For this reason, a second approach is presented in which an empirical POD basis is defined in each element of the mesh. This local approach is more intrusive with respect to the global approach but it is able to capture better the nonlinearities related to advection. The two approaches are tested and compared on the inviscid compressible flow around a gas‐turbine cascade and on the compressible turbulent flow around a wind turbine airfoil.     

A non-linearly stable implicit finite element algorithm for hypersonic aerodynamics

A generalized curvilinear co-ordinate Taylor weak statement implicit finite element algorithm is developed for the two-dimensional and axisymmetric compressible Navier-Stokes equations for ideal and reacting gases. For accurate hypersonic simulation, air is modelled as a mixture of five perfect gases, i.e. molecular and atomic oxygen and nirogen as well as nitric oxide. The associated pressure is then determined via Newton solution of the classical chemical equilibrium equation system. The directional semi-discretization is achieved using an optimal metric data Galerkin finite element weak statement, on a developed ‘companion conservation law system’, permitting classical test and trial space definitions. Utilizing an implicit Runge-Kutta scheme, the terminal algorithm is then non-linearly stable, and second-order accurate in space and time on arbitrary curvilinear co-ordinates. Subsequently, a matrix tensor product factorization procedure permits an efficient numerical linear algebra handling for large Courant numbers. For ideal- and real-gas hypersonic flows, the algorithm generates essentially non-oscillatory numerical solutions in the presence of strong detached shocks and boundary layer inviscid flow interactions.     

A hybrid stabilization approach for reduced-order models of compressible flows with shock-vortex interaction

Model order reduction approaches, such as proper orthogonal decomposition (POD)-Galerkin projection, provide a systematic manner to construct Reduced-Order Models (ROM) from pregenerated high-fidelity datasets. The current study focuses on the stabilization of ROMs built from high-fidelity simulation data of a supersonic flow passing a circular cylinder, which features strong interactions between shockwaves and vortices. As shown in previous literatures and the current study, an implicit subspace correction (ISC) method is efficient in the stabilization of similar problems, but its accuracy is not consistent when applied on different ROMs; on the other hand, an eigenvalue reassignment (ER) method delivers superb accuracy when the mode number is small, but becomes too expensive and less robust as the number increases. A Hybrid method is proposed here to balance the computational cost while improving the overall robustness/accuracy in ROM stabilization. The Hybrid method first handles the majority of the modes using the ISC method and then applies the ER method to fine tune a smaller number of modes under a constraint for accuracy. Furthermore, when the typical L2 inner product is changed to a symmetry inner product in both POD computation and Galerkin projection, the performance of the stabilized ROMs is substantially improved for all methods.     

Generalized finite element method using proper orthogonal decomposition

A methodology is presented for generating enrichment functions in generalized finite element methods (GFEM) using experimental and/or simulated data. The approach is based on the proper orthogonal decomposition (POD) technique, which is used to generate low‐order representations of data that contain general information about the solution of partial differential equations. One of the main challenges in such enriched finite element methods is knowing how to choose, a priori, enrichment functions that capture the nature of the solution of the governing equations. POD produces low‐order subspaces, that are optimal in some norm, for approximating a given data set. For most problems, since the solution error in Galerkin methods is bounded by the error in the best approximation, it is expected that the optimal approximation properties of POD can be exploited to construct efficient enrichment functions. We demonstrate the potential of this approach through three numerical examples. Best‐approximation studies are conducted that reveal the advantages of using POD modes as enrichment functions in GFEM over a conventional POD basis. Copyright © 2009 John Wiley & Sons, Ltd.     

Efficient solution of time‐domain boundary integral equations arising in sound‐hard scattering

We consider the efficient numerical solution of the three‐dimensional wave equation with Neumann boundary conditions via time‐domain boundary integral equations. A space‐time Galerkin method with C^∞‐smooth, compactly supported basis functions in time and piecewise polynomial basis functions in space is employed. We discuss the structure of the system matrix and its efficient parallel assembly. Different preconditioning strategies for the solution of the arising systems with block Hessenberg matrices are proposed and investigated numerically. Furthermore, a C++ implementation parallelized by OpenMP and MPI in shared and distributed memory, respectively, is presented. The code is part of the boundary element library BEM4I. Results of numerical experiments including convergence and scalability tests up to a thousand cores on a cluster are provided. The presented implementation shows good parallel scalability of the system matrix assembly. Moreover, the proposed algebraic preconditioner in combination with the FGMRES solver leads to a significant reduction of the computational time. Copyright © 2015 John Wiley & Sons, Ltd.     

Robin‐Neumann transmission conditions for fluid‐structure coupling: Embedded boundary implementation and parameter analysis

Partitioned procedures are appealing for solving complex fluid‐structure interaction (FSI) problems, as they allow existing computational fluid dynamics (CFD) and computational structural dynamics algorithms and solvers to be combined and reused. However, for problems involving incompressible flow and strong added‐mass effect (eg, heavy fluid and slender structure), partitioned procedures suffer from numerical instability, which typically requires additional subiterations between the fluid and structural solvers, hence significantly increasing the computational cost. This paper investigates the use of Robin‐Neumann transmission conditions to mitigate the above instability issue. Firstly, an embedded Robin boundary method is presented in the context of projection‐based incompressible CFD and finite element–based computational structural dynamics. The method utilizes operator splitting and a modified ghost fluid method to enforce the Robin transmission condition on fluid‐structure interfaces embedded in structured non–body‐conforming CFD grids. The method is demonstrated and verified using the Turek and Hron benchmark problem, which involves a slender beam undergoing large transient deformation in an unsteady vortex‐dominated channel flow. Secondly, this paper investigates the effect of the combination parameter in the Robin transmission condition, ie, α_f, on numerical stability and solution accuracy. This paper presents a numerical study using the Turek and Hron benchmark problem and an analytical study using a simplified FSI model featuring an Euler‐Bernoulli beam interacting with a two‐dimensional incompressible inviscid flow. Both studies reveal a trade‐off between stability and accuracy: smaller values of α_f tend to improve numerical stability, yet deteriorate the accuracy of the partitioned solution. Using the simplified FSI model, the critical value of α_f that optimizes this trade‐off is derived and discussed.     

On the a priori and a posteriori error analysis of a two‐fold saddle‐point approach for nonlinear incompressible elasticity

In this paper, we reconsider the a priori and a posteriori error analysis of a new mixed finite element method for nonlinear incompressible elasticity with mixed boundary conditions. The approach, being based only on the fact that the resulting variational formulation becomes a two‐fold saddle‐point operator equation, simplifies the analysis and improves the results provided recently in a previous work. Thus, a well‐known generalization of the classical Babu?ka–Brezzi theory is applied to show the well‐posedness of the continuous and discrete formulations, and to derive the corresponding a priori error estimate. In particular, enriched PEERS subspaces are required for the solvability and stability of the associated Galerkin scheme. In addition, we use the Ritz projection operator to obtain a new reliable and quasi‐efficient a posteriori error estimate. Finally, several numerical results illustrating the good performance of the associated adaptive algorithm are presented. Copyright © 2006 John Wiley & Sons, Ltd.     

Parallel computing of high‐speed compressible flows using a node‐based finite‐element method

An efficient parallel computing method for high‐speed compressible flows is presented. The numerical analysis of flows with shocks requires very fine computational grids and grid generation requires a great deal of time. In the proposed method, all computational procedures, from the mesh generation to the solution of a system of equations, can be performed seamlessly in parallel in terms of nodes. Local finite‐element mesh is generated robustly around each node, even for severe boundary shapes such as cracks. The algorithm and the data structure of finite‐element calculation are based on nodes, and parallel computing is realized by dividing a system of equations by the row of the global coefficient matrix. The inter‐processor communication is minimized by renumbering the nodal identification number using ParMETIS. The numerical scheme for high‐speed compressible flows is based on the two‐step Taylor–Galerkin method. The proposed method is implemented on distributed memory systems, such as an Alpha PC cluster, and a parallel supercomputer, Hitachi SR8000. The performance of the method is illustrated by the computation of supersonic flows over a forward facing step. The numerical examples show that crisp shocks are effectively computed on multiprocessors at high efficiency. Copyright © 2003 John Wiley & Sons, Ltd.     

Weight‐adjusted discontinuous Galerkin methods: Matrix‐valued weights and elastic wave propagation in heterogeneous media

Weight‐adjusted inner products are easily invertible approximations to weighted L² inner products. These approximations can be paired with a discontinuous Galerkin (DG) discretization to produce a time‐domain method for wave propagation which is low storage, energy stable, and high‐order accurate for arbitrary heterogeneous media and curvilinear meshes. In this work, we extend weight‐adjusted DG methods to the case of matrix‐valued weights, with the linear elastic wave equation as an application. We present a DG formulation of the symmetric form of the linear elastic wave equation, with upwind‐like dissipation incorporated through simple penalty fluxes. A semidiscrete convergence analysis is given, and numerical results confirm the stability and high‐order accuracy of weight‐adjusted DG for several problems in elastic wave propagation.     

Bending Analysis of Curved Sandwich Beams with Functionally Graded Core

In this article, bending analysis of curved sandwich beams with transversely and functionally graded (FG) core is studied. The Euler–Bernoulli beam theory is used to model the thin face-sheets and high-order shear theory is used to analyze the core. Equilibrium/field equations, compatibility and boundary conditions are used to derive the set of governing equations. The numerical solution of the governing nonlinear differential equations is based on the series Fourier–Galerkin method. Finally, the effect of geometric properties on radial deflection of core and the effect of core radius and Young's modulus on radial deflection, circumferential displacement, and stresses are investigated.     

Moving‐least‐squares‐particle hydrodynamics—I. Consistency and stability

The Smooth‐Particle‐Hydrodynamics (SPH) method is derived in a novel manner by means of a Galerkin approximation applied to the Lagrangian equations of continuum mechanics as in the finite‐element method. This derivation is modified to replace the SPH interpolant with the Moving‐Least‐Squares (MLS) interpolant of Lancaster and Saulkaskas, and define a new particle volume which ensures thermodynamic compatibility. A variable‐rank modification of the MLS interpolants which retains their desirable summation properties is introduced to remove the singularities that occur when divergent flow reduces the number of neighbours of a particle to less than the minimum required. A surprise benefit of the Galerkin SPH derivation is a theoretical justification of a common ad hoc technique for variable‐h SPH. The new MLSPH method is conservative if an anti‐symmetric quadrature rule for the stiffness matrix elements can be supplied. In this paper, a simple one‐point collocation rule is used to retain similarity with SPH, leading to a non‐conservative method. Several examples document how MLSPH renders dramatic improvements due to the linear consistency of its gradients on three canonical difficulties of the SPH method: spurious boundary effects, erroneous rates of strain and rotation and tension instability. Two of these examples are non‐linear Lagrangian patch tests with analytic solutions with which MLSPH agrees almost exactly. The examples also show that MLSPH is not absolutely stable if the problems are run to very long times. A linear stability analysis explains both why it is more stable than SPH and not yet absolutely stable and an argument is made that for realistic dynamic problems MLSPH is stable enough. The notion of coherent particles, for which the numerical stability is identical to the physical stability, is introduced. The new method is easily retrofitted into a generic SPH code and some observations on performance are made. Copyright © 1999 John Wiley & Sons, Ltd.     

A rapid simulation of nano‐particle transport in a two‐dimensional human airway using POD/Galerkin reduced‐order models

An approach is proposed for the rapid prediction of nano‐particle transport and deposition in the human airway, which requires the solution of both the Navier–Stokes and advection–diffusion equations and for which computational efficiency is a challenge. The proposed method builds low‐order models that are representative of the fully coupled equations by means of the Galerkin projection and proper orthogonal decomposition technique. The obtained reduced‐order models (ROMs) are a set of ordinary differential equations for the temporal coefficients of the basis functions. The numerical results indicate that the ROMs are highly efficient for the computation (the speedup factor is approximately 3 × 10³) and have reasonable accuracy compared with the full model (relative error of ≈7 × 10^?3). Using ROMs, the deposition of particles is studied for 1≤d_n≤100 nm, where d_n is the diameter of a nano‐particle. The effectiveness of this approach is promising for applications of health risk assessment. Copyright © 2015 John Wiley & Sons, Ltd.     

Weighted residual solutions in the time domain

The Galerkin and subdomain forms of the weighted residual method are used to generate recursive equations in time for the numerical solution of a system of ordinary differential equations. The single-step methods that result from a linear interpolation equation match currently available methods whose stability and oscillation properties are known. A three-level scheme developed by combining two linear elements is shown to be unconditionally unstable. Two of the three schemes obtained using a quadratic interpolation equation and quadratic weighting functions are also shown to be unconditionally unstable. The third scheme is unconditionally stable, but the calculated values for a numerical solution of u? + u = 0, u(0) = 1 are not as accurate as the values obtained using the single-step central difference method.     

Trefftz boundary elements—multi‐region formulations

This paper is concerned with an effective numerical implementation of the Trefftz boundary element method, for the analysis of two‐dimensional potential problems, defined in arbitrarily shaped domains. The domain is first discretized into multiple subdomains or regions. Each region is treated as a single domain, either finite or infinite, for which a complete set of solutions of the problem is known in the form of an expansion with unknown coefficients. Through the use of weighted residuals, this solution expansion is then forced to satisfy the boundary conditions of the actual domain of the problem, leading thus to a system of equations, from which the unknowns can be readily determined. When this basic procedure is adopted, in the analysis of multiple‐region problems, proper boundary integral equations must be used, along common region interfaces, in order to couple to each other the unknowns of the solution expansions relative to the neighbouring regions. These boundary integrals are obtained from weighted residuals of the coupling conditions which allow the implementation of any order of continuity of the potential field, across the interface boundary, between neighbouring regions. The technique used in the formulation of the region‐coupling conditions drives the performance of the Trefftz boundary element method. While both of the collocation and Galerkin techniques do not generate new unknowns in the problem, the technique of Galerkin presents an additional and unique feature: the size of the matrix of the final algebraic system of equations which is always square and symmetric, does not depend on the number of boundary elements used in the discretization of both the actual and region‐interface boundaries. This feature which is not shared by other numerical methods, allows the Galerkin technique of the Trefftz boundary element method to be effectively applied to problems with multiple regions, as a simple, economic and accurate solution technique. A very difficult example is analysed with this procedure. The accuracy and efficiency of the implementations described herein make the Trefftz boundary element method ideal for the study of potential problems in general arbitrarily‐shaped two‐dimensional domains. Copyright © 1999 John Wiley & Sons, Ltd.     

A Galerkin boundary element formulation with dual reciprocity for elastodynamics

In this paper, a Galerkin boundary integral equation method for two‐dimensional elastodynamic problems is presented. The formulation makes use of a static fundamental solution to weight the dynamic equilibrium equations. Following the Galerkin approach, the equations are weighted again with the interpolation functions used in the discretization of the unknowns. For the numerical integration, a regularization process is followed to deal with the integrands containing strong singularities. The implementation of the dual reciprocity method to transfer the domain integrals to the boundary is also presented in the context of the Galerkin formulation. Finally, the Hubolt integration scheme was used for the time‐marching process. Several numerical examples are presented to demonstrate the accuracy of the method. Copyright © 2000 John Wiley & Sons, Ltd.