求解双曲守恒律方程的WENO型熵相容格式

通过在单元交界面处进行高阶WENO重构,得到了一种求解双曲型守恒律方程的WENO型熵相容格式。用该格式对一维Burgers方程和Euler方程进行数值模拟,结果表明,该格式具有高精度、基本无振荡性等特点。     

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems. Of critical importance, essentially non-oscillatory solutions are uniformly attained for a range of supercritical flow situations with shocks.     

Direct Riemann solvers for the time-dependent Euler equations

Approaches for finding direct, approximate solutions to the Riemann problem are presented. These result in three approximate Riemann solvers. Here we discuss the time-dependent Euler equations but the ideas are applicable to other systems. The approximate solvers are (i) assessed on local Riemann problems with exact solutions and (ii) used in conjunction with the Weighted Average Flux (WAF) method to solve the two-dimensional Euler equations numerically. The resulting numerical technique is assessed on a shock reflection problem. Comparison with experimental observation is carried out.     

一种三阶有限体积法及其在欠膨胀射流激波结构数值模拟中的应用

以喷管出口欠膨胀射流为研究对象,在Lagrange坐标系下建立欠膨胀射流二维积分形式的流动方程。通过在单元交接面处进行三阶ENO(essentially nonoscillatory)格式插值,构造得到一种适用于求解该方程的三阶ENO有限体积法。采用该格式对一维Sod激波管算例和喷管出口欠膨胀射流进行数值计算。计算结果表明,该方法具有高精度、基本无振荡的特点,能很好地捕捉包含激波、滑移线以及三波交点等复杂流场波系结构。计算得到的波系结构中马赫盘的位置与实验结果吻合很好,相对误差小于1.1%。     

Dynamics of axially accelerating beams with multiple supports

This study represents the transverse vibrations of an axially accelerating Euler–Bernoulli beam resting on multiple simple supports. This is one of the examples of a system experiencing Coriolis acceleration component that renders such systems gyroscopic. A small harmonic variation with a constant mean value for the axial velocity is assumed in the problem. The immovable supports introduce nonlinear terms to the equations of motion due to stretching of neutral axis. The method of multiple scales is directly applied to the equations of motion obtained for the general case. Natural frequency equations are presented for multiple support case. Principal parametric resonances and combination resonances are discussed. Solvability conditions are presented for different cases. Stability analysis is conducted for the solutions; approximate stable and unstable regions are identified. Some numerical examples are presented to show the effects of axial speed, number of supports, and their locations.     

A NON-OSCILLATORY NO-FREE-PARAMETER FINITE ELEMENT AND ITS APPLICATIONS IN CDF

A non-oscillatory no-free-parameter finite element method (NNFEM) is presented based on the consideration of wave propagation characteristic in different characteristic directions across a strong discontinuity through flux vector splitting in order to satisfy the increasing entropy condition. The algorithm is analysed in detail for the one-dimensional (1D) Euler equation and then extended to the 2D, axisymmetric and 3D Euler and Navier–Stokes equations. Its applications in various cases—in viscid oblique shock wave reflection, flow over a forward step, axisymmetric free jet flow, supersonic flows over 2D and 3D rectangular cavities—are given. These computational results show that the present NNFEM is efficient in practice and stable in operations and is especially capable of giving good resolution in simulating complicated separated and vortical flows interacting with shock waves. © 1997 by John Wiley & Sons, Ltd.     

Implicit,non-oscillatory containing no free parameters and dissipative (INND) scheme

An implicit,non-oscillatory.containing no free parameter and dissipative(INND)scheme solving.Navier Stokes equations is developed.This scheme is one of total variationdiminishing(TVD)algorithms.The results show that this scheme is applicable for solvingNavier Stokes equations,and that the shock-capturing capabilbility and the convergencerate are satisfactory.     

Sloshing simulation of standing wave with time-independent finite difference method for Euler equations

The numerical solutions of standing waves for Euler equations with the nonlinear free surface boundary condition in a two-dimensional(2D) tank are studied.The irregular tank is mapped onto a fixed square domain through proper mapping functions. A staggered mesh system is employed in a 2D tank to calculate the elevation of the transient fluid.A time-independent finite difference method,which is developed by Bang-fuh Chen,is used to solve the Euler equations for incompressible and inviscid fluids.The numer...     

New Steady and Self-Similar Solutions of the Euler Equations

Exact steady and self-similar solutions of the Euler equations are considered, which possess the property of partial invariance with respect to a certain six-parameter Lie group. New examples of vortex motion of a swirled liquid in curved channels are presented. A classification is given for self-similar solutions of the reduced system with two independent variables, which admits a three-parameter group of extensions, whereas the initial system of the Euler equations possesses a two-parameter group.     

Numerical simulation of vortical ideal fluid flow through curved channel

A numerical algorithm to study the boundary‐value problem in which the governing equations are the steady Euler equations and the vorticity is given on the inflow parts of the domain boundary is developed. The Euler equations are implemented in terms of the stream function and vorticity. An irregular physical domain is transformed into a rectangle in the computational domain and the Euler equations are rewritten with respect to a curvilinear co‐ordinate system. The convergence of the finite‐difference equations to the exact solution is shown experimentally for the test problems by comparing the computational results with the exact solutions on the sequence of grids. To find the pressure from the known vorticity and stream function, the Euler equations are utilized in the Gromeka–Lamb form. The numerical algorithm is illustrated with several examples of steady flow through a two‐dimensional channel with curved walls. The analysis of calculations shows strong dependence of the pressure field on the vorticity given at the inflow parts of the boundary. Plots of the flow structure and isobars, for different geometries of channel and for different values of vorticity on entrance, are also presented. Copyright © 2003 John Wiley & Sons, Ltd.     

WELL-POSEDNESS OF INITIAL VALUE PROBLEM FOR EULER EQUATIONS OF INVISCID COMPRESSIBLE ADIABATIC FLUID ＊

Introduction Itisreasonabletoconsiderarealfluidasanidealoneinfluidmechanicsundermany conditions.Forinstance,ingeneral,forthedistributionofthefluidfieldaroundtheaerocraft,mostpartofthefluidfieldmayberegardedasidealfluidexceptforasmallpartwhereeffects ofviscousandheatconductioninthethinlayernearthesurfacemustbeconsidered.Evenif thefluidiscompletelysupposedasidealonethroughoutthefluidfield,thequitereasonable resultsarealsogained,thereforestudyingidealfluidhasnotonlytheoreticalsignificancebut also…     

Nonstationary cylindrical vortex in an ideal fluid

Equations of rotationally symmetric motion of an ideal incompressible fluid are considered. A class of solutions to these equations, described by a hyperbolic equation of the fourth order with one space variable, for which an initial boundary-value problem is formulated, is distinguished. The new class of exact solutions of the Euler equations was used to describe the a nonstationary cylindrical vortex in an ideal fluid.     

Assessment of shock capturing schemes for discontinuous Galerkin method

This paper carries out systematical investigations on the performance of several typical shock-capturing schemes for the discontinuous Galerkin （DG） method, including the total variation bounded （TVB） limiter and three artificial diffusivity schemes （the basis function-based （BF） scheme, the face residual-based （FR） scheme, and the element residual-based （ER） scheme）. Shock-dominated flows （the Sod problem, the Shu- Osher problem, the double Mach reflection problem, and the transonic NACA0012 flow） are considered, addressing the issues of accuracy, non-oscillatory property, dependence on user-specified constants, resolution of discontinuities, and capability for steady solutions. Numerical results indicate that the TVB limiter is more efficient and robust, while the artificial diffusivity schemes are able to preserve small-scale flow structures better. In high order cases, the artificial diffusivity schemes have demonstrated superior performance over the TVB limiter.     

Non-oscillation Criteria for Hypoelastic Models under Simple Shear Deformation

Ten objective rates, spinning or non-spinning, are critically examined from the viewpoint of Sturm's theorems in ordinary differential equations. Upon developing implication relations of oscillatory, non-oscillatory, and disconjugate behavior, we establish oscillation and non-oscillation criteria which pick out the objective stress rates that lead to oscillatory and non-oscillatory responses in simple shear deformation, respectively. Among the hypoelastic equations associated with the spinning objective rates examined, the Jaumann equation is an oscillatory minorant, the homogeneous Xiao–Bruhns–Meyers equation is a non-oscillatory majorant, and the homogeneous Green–Naghdi equation is a disconjugate majorant. If (Sturm comparable) non-spinning objective rates are also taken into consideration, then the Durban–Baruch equation becomes an oscillatory minorant, but the other two equations remain to play the same roles. The Jaumann equation is a Sturm majorant for all the other nine homogeneous hypoelastic equations, and the homogeneous Szabó–Balla-2 equation is a Sturm minorant for all the other nine homogeneous hypoelastic equations. Most of the solutions of the zeroth-grade hypoelastic equations at simple shear have already been published, except for those of Szabó and Balla, to which the closed-form solutions are derived here. Moreover, all solutions are extended to include the effect of initial stresses. This revised version was published online in August 2006 with corrections to the Cover Date.     

A General Formulation and Solution Scheme for Fractional Optimal Control Problems

Accurate modeling of many dynamic systems leads to a set of Fractional Differential Equations (FDEs). This paper presents a general formulation and a solution scheme for a class of Fractional Optimal Control Problems (FOCPs) for those systems. The fractional derivative is described in the Riemann–Liouville sense. The performance index of a FOCP is considered as a function of both the state and the control variables, and the dynamic constraints are expressed by a set of FDEs. The Calculus of Variations, the Lagrange multiplier, and the formula for fractional integration by parts are used to obtain Euler–Lagrange equations for the FOCP. The formulation presented and the resulting equations are very similar to those that appear in the classical optimal control theory. Thus, the present formulation essentially extends the classical control theory to fractional dynamic system. The formulation is used to derive the control equations for a quadratic linear fractional control problem. An approach similar to a variational virtual work coupled with the Lagrange multiplier technique is presented to find the approximate numerical solution of the resulting equations. Numerical solutions for two fractional systems, a time-invariant and a time-varying, are presented to demonstrate the feasibility of the method. It is shown that (1) the solutions converge as the number of approximating terms increase, and (2) the solutions approach to classical solutions as the order of the fractional derivatives approach to 1. The formulation presented is simple and can be extended to other FOCPs. It is hoped that the simplicity of this formulation will initiate a new interest in the area of optimal control of fractional systems.     

Upwinded finite difference schemes for dense snow avalanche modeling

Avalanche dynamics models are used by engineers and land‐use planners to predict the reach and destructive force of snow avalanches. These models compute the motion of the flowing granular core of dense snow avalanches from initiation to runout. The governing differential equations for the flow height and velocity can be approximated by a hyperbolic system of equations of first‐order with respect to time, formally equivalent to the Euler equations of a one‐dimensional isentropic gas. In avalanche practice these equations are presently solved analytically by making restrictive assumptions regarding mountain topography and avalanche flow behaviour. In this article the one‐dimensional dense snow avalanche equations are numerically solved using the conservative variables and stable upwinded and total variation diminishing finite difference schemes. The numerical model is applied to simulate avalanche motion in general terrain. The proposed discretization schemes do not use artificial damping, an important requirement for the application of numerical models in practice. In addition, non‐physical M‐wave solutions are not encountered as in previous attempts to solve this problem using Eulerian finite difference methods and non‐conservative variables. The simulation of both laboratory experiments and a field case study are presented to demonstrate the newly developed discretization schemes. Copyright © 2000 John Wiley & Sons, Ltd.     

Kinetic theory based multi-level adaptive finite difference WENO schemes for compressible Euler equations

In this paper we proposed the kinetic framework based fifth-order adaptive finite difference WENO schemes abbreviated as WENO-AO-K schemes to solve the compressible Euler equations, which are quasi-linear hyperbolic equations that can admit discontinuous solutions like shock and contact waves. The formulation of the proposed schemes is based on the kinetic theory where one can recover the Euler equations by applying a suitable moment method strategy to the Boltzmann equation. The kinetic flux vector splitting strategy is used in WENO-AO framework, which produces the computationally expensive error and exponential functions. Thus, to reduce the computational cost, a physically more relevant peculiar velocity based splitting strategy is used, which is more efficient than the kinetic flux vector splitting. High order of accuracy in time is achieved using the third-order total variation diminishing Runge–Kutta (TVD-RK) scheme. Several one- and two-dimensional test cases are solved for the compressible Euler equations using the proposed fifth-order WENO-AO-K schemes and the results are compared with conventional WENO-AO scheme. Proposed schemes capture the complex flow features in a smooth region accurately, and discontinuity is also well resolved. Error analysis of the proposed schemes shows optimal convergence rates in various norms.     

A new efficient method to evaluate exact stiffness and mass matrices of non-uniform beams resting on an elastic foundation

In this paper, a new efficient method to evaluate the exact stiffness and mass matrices of a non-uniform Bernoulli–Euler beam resting on an elastic Winkler foundation is presented. The non-uniformity may result from variable cross-section and/or from inhomogeneous linearly elastic material. It is assumed that there is no abrupt variation in the cross-section of the beam so that the Euler–Bernoulli theory is valid. The method is based on the integration of the exact shape functions which are derived from the solution of the axial deformation problem of a non-uniform bar and the bending problem of a non-uniform beam which are both formulated in terms of the two displacement components. The governing differential equations are uncoupled with variable coefficients and are solved within the framework of the analog equation concept. According to this, the two differential equations with variable coefficients are replaced by two linear ones pertaining to the axial and transverse deformation of a substitute beam with unit axial and bending stiffness, respectively, under ideal load distributions. The key point of the method is the evaluation of the two ideal loads which in this work is achieved by approximating them by two polynomials. More specifically, the axial ideal load is approximated by a linear polynomial while the transverse one by a cubic polynomial. The numerical implementation of the method is simple, and the results are compared favorably to those obtained by exact solutions available in literature.     

Application of Lie transformation group methods to classical linear theories of rods and plates

In the present paper, a class of partial differential equations governing various rod and plate theories of Bernoulli–Euler and Poisson–Kirchhoff type is studied by Lie transformation group methods. A system of equations determining the generators of the admitted point Lie groups (symmetries) is derived and the general statement of the associated group-classification problem is given. A simple relation is deduced allowing to recognize easily the variational symmetries among the “ordinary” symmetries of a self-adjoint equation of the class examined. Explicit formulae for the conserved currents of the corresponding (via Bessel-Hagen’s extension of Noether’s theorem) conservation laws are suggested. Solutions of group-classification problems are given for subclasses of equations of the foregoing type governing stability and vibration of rods, fluid conveying pipes and plates resting on variable elastic foundations. The obtained group-classification results are used to derive conservation laws and group-invariant solutions readily applicable in rod dynamics and plate statics and dynamics. New generalized symmetries and conservation laws for the theories of Timoshenko beams, Reissner–Mindlin plates and three-dimensional elastostatics are presented.     

Modelling gas dynamics in 1D ducts with abrupt area change

Most gas dynamic computations in industrial ducts are done in one dimension with cross-section-averaged Euler equations. This poses a fundamental difficulty as soon as geometrical discontinuities are present. The momentum equation contains a non-conservative term involving a surface pressure integral, responsible for momentum loss. Definition of this integral is very difficult from a mathematical standpoint as the flow may contain other discontinuities (shocks, contact discontinuities). From a physical standpoint, geometrical discontinuities induce multidimensional vortices that modify the surface pressure integral. In the present paper, an improved 1D flow model is proposed. An extra energy (or entropy) equation is added to the Euler equations expressing the energy and turbulent pressure stored in the vortices generated by the abrupt area variation. The turbulent energy created by the flow–area change interaction is determined by a specific estimate of the surface pressure integral. Model’s predictions are compared with 2D-averaged results from numerical solution of the Euler equations. Comparison with shock tube experiments is also presented. The new 1D-averaged model improves the conventional cross-section-averaged Euler equations and is able to reproduce the main flow features.