Fourth-order central compact scheme for the numerical solution of incompressible Navier–Stokes equations

This paper provides an implicit central compact scheme for the numerical solution of incompressible Navier–Stokes equations. The solution procedure is based on the artificial compressibility method that transforms the governing equations into a hyperbolic-parabolic form. A fourth-order central compact scheme with a sixth-order numerical filtering is used for the discretization of convective terms and fourth-order central compact scheme for the viscous terms. Dual-time stepping approach is applied to time discretization with backward Euler difference scheme to the pseudo-time derivative, and three point second-order backward difference scheme to the physical time derivative. An approximate factorization-based alternating direction implicit scheme is used to solve the resulting block tridiagonal system of equations. The accuracy and efficiency of the proposed numerical method is verified by simulating several two-dimensional steady and unsteady benchmark problems.     

A numerical method for the vorticity-velocity Navier-Stokes equations in two and three dimensions

This paper provides a numerical method for solving the steady-state vorticity-velocity Navier-Stokes equations in two and three dimensions. The vorticity transport equation is considered together with a Poisson equation for the velocity vector, the latter equation being parabolized in time according to the false transient approach. The two vector equations are discretized in time using the implicit Euler time stepping and the delta form of Beam and Warming. A staggered-grid spatial discretization is employed in conjunction with a deferred correction procedure. Second-order-accurate central differences are used to approximate the steady-state residuals, written in conservative form for accuracy reasons, whereas upwind differences are used for the advection terms in the implicit operator, to obtain diagonally-dominant tridiagonal systems. The discrete equations are solved sequentially by means of a robust alternating direction line-Gauss-Seidel iteration procedure combined with a simple multigrid strategy. For the model driven-cavity-flow problem in two and three dimensions, the method is found to be efficient and very accurate. For the first time, the three-dimensional discrete vorticity and velocity fields, computed using a Poisson equation for the velocity vector, are both solenoidal and satisfy their mutual relationship, exactly.     

A Fourth Order Scheme for Incompressible Boussinesq Equations

A fourth order finite difference method is presented for the 2D unsteady viscous incompressible Boussinesq equations in vorticity-stream function formulation. The method is especially suitable for moderate to large Reynolds number flows. The momentum equation is discretized by a compact fourth order scheme with the no-slip boundary condition enforced using a local vorticity boundary condition. Fourth order long-stencil discretizations are used for the temperature transport equation with one-sided extrapolation applied near the boundary. The time stepping scheme for both equations is classical fourth order Runge–Kutta. The method is highly efficient. The main computation consists of the solution of two Poisson-like equations at each Runge–Kutta time stage for which standard FFT based fast Poisson solvers are used. An example of Lorenz flow is presented, in which the full fourth order accuracy is checked. The numerical simulation of a strong shear flow induced by a temperature jump, is resolved by two perfectly matching resolutions. Additionally, we present benchmark quality simulations of a differentially-heated cavity problem. This flow was the focus of a special session at the first MIT conference on Computational Fluid and Solid Mechanics in June 2001.     

Analysis of unsteady compressible boundary layer flow via finite elements

This paper is concerned with the discrete formulation and numerical solution of unsteady compressible boundary layer flows using the Galerkin-finite element method. Linear interpolation functions for the velocity, density, temperature and pressure are used in the momentum equation and equations of continuity, energy and state. The coupled nonlinear finite element equations are approximated by a third order Taylor series expansion as temporal operator to integrate in time with Newton-Raphson type iterations performed until convergence within each time step. As an example, a boundary layer problem of a perfect gas behind a normal shock wave is solved. A comparison of the results with those by other method indicates a favorable agreement.     

Numerical solution of incompressible Navier–Stokes equations using a fractional-step approach

A fractional step method for the solution of steady and unsteady incompressible Navier–Stokes equations is outlined. The method is based on a finite-volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (third and fifth) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when fifth-order upwind differencing and a modified production term in the Baldwin–Barth one-equation turbulence model are used with adequate grid resolution.     

A square root of a matrix approach to obtain the solution to a steady-state matrix Riccati equation

By utilizing the square root of a matrix approach, a method of generating the solution to a steady-state matrix Riccati type equation (under certain restrictions) is presented. This approach not only yields a closed form expression for the Riccati solution, but also converts the original Riccati equation into other equations which may have numerical or computational advantages. An example is worked out for a second-order case.     

Steady and unsteady flow past a rotating circular cylinder at low Reynolds numbers

Results of calculations of the steady and unsteady flows past a circular cylinder which is rotating with constant angular velocity and translating with constant linear velocity are presented. The motion is assumed to be two-dimensional and to be governed by the Navier-Stokes equations for incompressible fluids. For the unsteady flow, the cylinder is started impulsively from rest and it is found that for low Reynolds numbers the flow approaches a steady state after a large enough time. Detailed results are given for the development of the flow with time for Reynolds numbers 5 and 20 based on the diameter of the cylinder. For comparison purposes the corresponding steady flow problem has been solved. The calculated values of the steady-state lift, drag and moment coefficients from the two methods are found to be in good agreement. Notable, however, are the discrepancies between these results and other recent numerical solutions to the steady-state Navier-Stokes equations. Some unsteady results are also given for the higher Reynolds numbers of 60, 100 and 200. In these cases the flow does not tend to be a steady state but develops a periodic pattern of vortex shedding.     

On the construction of K-consistent difference schemes of gas dynamics

Previously it was found by the present authors that the use of non-classical equations of state in gas dynamic computations by Eulerian difference schemes may give rise to spurious oscillations of the numerical solution in the vicinity of contact discontinuities. Two accurate methods for elimination of these oscillations are proposed below. One of these methods is based on a new definition of the K-consistency of difference schemes. Another method makes use of difference approximations of an unsteady equation for pressure which is a well-known consequence of the Euler equations. The effectiveness of the methods proposed is illustrated by computational results.     

Numerical solution of unsteady Navier–Stokes equations on curvilinear meshes

The objective of the present work is to extend our FDS-based third-order upwind compact schemes by Shah et al. (2009) [8] to numerical solutions of the unsteady incompressible Navier–Stokes equations in curvilinear coordinates, which will save much computing time and memory allocation by clustering grids in regions of high velocity gradients. The dual-time stepping approach is used for obtaining a divergence-free flow field at each physical time step. We have focused on addressing the crucial issue of implementing upwind compact schemes for the convective terms and a central compact scheme for the viscous terms on curvilinear structured grids. The method is evaluated in solving several two-dimensional unsteady benchmark flow problems.     

Structural optimization using global stress-deviation objective function via the level-set method

The paper deals with minimum stress design using a novel stress-related objective function based on the global stress-deviation measure. The shape derivative, representing the shape sensitivity analysis of the structure domain, is determined for the generalized form of the global stress-related objective function. The optimization procedure is based on the domain boundary evolution via the level-set method. The elasticity equations are, instead of using the usual ersatz material approach, solved by the extended finite element method. The Hamilton-Jacobi equation is solved using the streamline diffusion finite element method. The use of finite element based methods allows a unified numerical approach with only one numerical framework for the mechanical problem as also for the boundary evolution stage. The numerical examples for the L-beam benchmark and the notched beam are given. The results of the structural optimization problem, in terms of maximum von Mises stress corresponding to the obtained optimal shapes, are compared for the commonly used global stress measure and the novel global stress-deviation measure, used as the stress-related objective functions.     

Solution-adaptivity in modelling complex shallow flows

A solution-adaptive algorithm is presented and tested for the shallow water equations. Specifically, we focus on the two-dimensional modelling of wind-induced hydrodynamics in shallow waters, characterised by a strong influence of variable bed topography and aquatic vegetation. The numerical solution is obtained using a Godunov-type finite-volume scheme on a hierarchical Cartesian mesh, with local time stepping. A simple, but robust algorithm based on the velocity gradient is proposed to control the dynamic mesh adaptation. Simulations on a representative steady-state and unsteady benchmark problem show that solution-adaptivity is successful in reallocating cells to where they improve global accuracy more efficiently. The algorithm is then applied to model the wind-driven circulations in Lake Neusiedl, proving the algorithm’s robustness in resolving complex geometry and vegetation cover.     

An analysis of the time integration algorithms for the finite element solutions of incompressible Navier–Stokes equations based on a stabilised formulation

This work is concerned with the analysis of time integration procedures for the stabilised finite element formulation of unsteady incompressible fluid flows governed by the Navier–Stokes equations. The stabilisation technique is combined with several different implicit time integration procedures including both finite difference and finite element schemes. Particular attention is given to the generalised-α method and the linear discontinuous in time finite element scheme. The time integration schemes are first applied to two model problems, represented by a first order differential equation in time and the one dimensional advection–diffusion equation, and subjected to a detailed mathematical analysis based on the Fourier series expansion. In order to establish the accuracy and efficiency of the time integration schemes for the Navier–Stokes equations, a detailed computational study is performed of two standard numerical examples: unsteady flow around a cylinder and flow across a backward facing step. It is concluded that the semi-discrete generalised-α method provides a viable alternative to the more sophisticated and expensive space–time methods for simulations of unsteady flows of incompressible fluids governed by the Navier–Stokes equations.     

Method of steady-state simulation of a gas network

A method is described for the steady-state simulation of an arbitrary gas network. Basic equations for a steady-state analysis are given. For controllable elements of the network, such as sources, compressors, regulators, and valves, the equation linking the inlet pressure, outlet pressure, and the flow through each unit is formulated. The set of equations for the steady-state simulation of a gas network and numerical methods for their solution are also described. Two examples are given to illustrate the properties of the method.     

On the application of the minimum polynomial extrapolation method to incompressible flows with heat transfer

In this paper, the Minimum Polynomial Extrapolation method (MPE) is used to accelerate the convergence of the Characteristic–Based–Split (CBS) scheme for the numerical solution of steady state incompressible flows with heat transfer. The CBS scheme is a fractional step method for the solution of the Navier–Stokes equations while the MPE method is a vector extrapolation method which transforms the original sequence into another sequence converging to the same limit faster then the original one without the explicit knowledge of the sequence generator. The developed algorithm is tested on a two-dimensional benchmark problem (buoyancy–driven convection problem) where the Navier–Stokes equations are coupled with the temperature equation. The obtained results show the feature of the extrapolation procedure to the CBS scheme and the reduction of the computational time of the simulation.     

Nonlinear circuit analysis using the method of harmonic balance—A review of the art. Part I. Introductory concepts

The harmonic balance method is a technique for the numerical solution of nonlinear analog circuits operating in a periodic, or quasi-periodic, steady-state regime. The method can be used to efficiently derive the continuous-wave response of numerous nonlinear microwave components including amplifiers, mixers, and oscillators. Its efficiency derives from imposing a predetermined steady-state form for the circuit response onto the nonlinear equations representing the network, and solving for the set of unknown coefficients in the response equation. Its attractiveness for nonlinear microwave applications results from its speed and ability to simply represent the dispersive, distributed elements that are common at high frequencies. The last decade has seen the development and application of harmonic balance techniques to model analog circuits, particularly microwave circuits. The first part of this paper reviews the fundamental achievements made during this time. The second part covers the extension of the method to quasi-periodic regimes, optimization analysis, and practical application. A critical assessment of the various types of harmonic balance techniques is given. The different sampling and Fourier transform methods are compared, and numerical speed and precision results are given enabling a quantitative analysis of the merits of the major variants of the harmonic balance technique. Examples of designs which have been modeled using the harmonic balance technique and built both in hybrid and MMIC form are presented.     

相关观测融合稳态Kalman滤波器及其最优性

对于带相关观测噪声和带不同观测阵的多传感器系统, 用加权最小二乘 (Weighted least squares, WLS) 法提出了两种相关观测融合稳态Kalman滤波算法. 其原理是用加权局部观测方程得到一个融合观测方程, 它伴随状态方程实现观测融合稳态Kalman滤波. 用信息滤波器证明了它们功能等价于集中式融合稳态Kalman滤波算法, 因而具有渐近全局最优性, 且可减少计算负担. 它们可应用于多通道自回归滑动平均 (Autoregressive moving average, ARMA) 信号观测融合滤波和反卷积. 两个数值仿真例子验证了它们的功能等价性.     

A Newton method for the resolution of steady stochastic Navier-Stokes equations

We present a Newton method to compute the stochastic solution of the steady incompressible Navier-Stokes equations with random data (boundary conditions, forcing term, fluid properties). The method assumes a spectral discretization at the stochastic level involving a orthogonal basis of random functionals (such as Polynomial Chaos or stochastic multi-wavelets bases). The Newton method uses the unsteady equations to derive a linear equation for the stochastic Newton increments. This linear equation is subsequently solved following a matrix-free strategy, where the iterations consist in performing integrations of the linearized unsteady Navier-Stokes equations, with an appropriate time scheme to allow for a decoupled integration of the stochastic modes. Various examples are provided to demonstrate the efficiency of the method in determining stochastic steady solution, even for regimes where it is likely unstable.     

Development of a parallel implicit solver of fluid modeling equations for gas discharges

A parallel fully implicit PETSc-based fluid modeling equations solver for simulating gas discharges is developed. Fluid modeling equations include: the neutral species continuity equation, the charged species continuity equation with drift-diffusion approximation for mass fluxes, the electron energy density equation, and Poisson's equation for electrostatic potential. Except for Poisson's equation, all model equations are discretized by the fully implicit backward Euler method as a time integrator, and finite differences with the Scharfetter–Gummel scheme for mass fluxes on the spatial domain. At each time step, the resulting large sparse algebraic nonlinear system is solved by the Newton–Krylov–Schwarz algorithm. A 2D-GEC RF discharge is used as a benchmark to validate our solver by comparing the numerical results with both the published experimental data and the theoretical prediction. The parallel performance of the solver is investigated.     

Nonlinear circuit analysis using the method of harmonic balance—a review of the art. II. Advanced concepts

The harmonic balance method is a technique for the numerical solution of nonlinear analog circuits operating in a periodic, or quasi-periodic, steady-state regime. The method can be used to efficiently derive the continuous-wave response of numerous nonlinear microwave components including amplifiers, mixers, and oscillators. Its efficiency derives from imposing a predetermined steady-state form for the circuit response onto the nonlinear equations representing the network, and solving for the set of unknown coefficients in the response equation. Its attractiveness for nonlinear microwave applications results from its speed and ability to simply represent the dispersive, distributed elements that are common at high frequencies. The last decade has seen the development and application of harmonic balance techniques to model analog circuits, particularly microwave circuits. The first part of this article reviewed the fundamental achievements made during this time. In this part, the extension of the method to quasi-periodic regimes, optimization analysis, oscillator analysis, studies of various convergence strategies, and practical applications are discussed. A critical assessment of the various types of harmonic balance techniques is given. Examples of designs which have been modeled using the harmonic balance technique and built both in hybrid and MMIC form are presented.     

Pade type model order reduction for multivariable systems using routh approximation

A method of obtaining reduced order models for multivariable systems is described. It is shown that the method has several advantages, e.g. the reduced order models retain the steady-state value and stability of the original system. Irrespective of whether the original multivariable system is described in state space form or in the transfer matrix form, the proposed method yields the reduced order models in state space form. In this method the Routh approximation is used to formulate the common denominator polynomial of a reduced order model. This is used to describe the structure of A_r matrix. The matrices B_r and C_r are chosen appropriately and some of the elements of B_r/C_r matrices are specified in such a way that after matching time moments/Markov parameters, the resulting equations are linear in the unknown elements of B_r and C_r matrices. The procedure is illustrated via a numerical example.