Interpolated characteristic interface conditions for zonal grid refinement of high-order multi-block computations

When numerically calculating fluid flows around complex geometries using the high-order finite-difference method on structured grids, grid singularities are frequently observed, even if the grids are carefully generated. In multi-block computations with the generalised characteristic interface conditions (GCIC), decomposed subdomains (blocks) do not overlap but are connected by the inviscid characteristic relations. In the original theory of the GCIC, discontinuity of the metrics on the interface can be accurately treated; however, discontinuity of the grid lines on the interface is not allowed. This article proposes a theoretical extension to the GCIC by incorporating high-order interpolation methods; this extension is called GCIC with interpolation (GCIC + I). The basic concept and solution procedure of the multi-block computation with the GCIC + I are presented in detail, and two benchmark tests are conducted to validate the proposed theory.     

A study of upstream-weighted high-order differencing for approximation to flow convection

This paper is concerned with a number of upstream-weighted second- and third-order difference schemes. Also considered are the conventional upwind and central difference schemes for comparison. It commences with a general difference equation which unifies all the given first-, second- and third-order schemes. The various schemes are evaluated through the use of the general equation. The unboundedness and accuracy of the solutions by the difference schemes are assessed via various analyses: examination of the coefficients of the difference equation, Taylor series truncation error analysis, study of the upstream connection to numerical diffusion, single-cell analysis. Finally, the difference schemes are tested on one- and two-dimensional model problems. It is shown that the high-order schemes suffer less from the problem of numerical diffusion than the first-order upwind difference scheme. However, unboundedness cannot be avoided in the solutions by these schemes. Among them the linear upwind difference scheme presents the best compromise between numerical diffusion and solution unboundedness.     

Gaussian beam formulations and interface conditions for the one-dimensional linear Schrödinger equation

Gaussian beams are asymptotic solutions of linear wave-like equations in the high frequency regime. This paper is concerned with the beam formulations for the Schrödinger equation and the interface conditions while beams pass through a singular point of the potential function. The equations satisfied by Gaussian beams up to the fourth order are given explicitly. When a Gaussian beam arrives at a singular point of the potential, it typically splits into a reflected wave and a transmitted wave. Under suitable conditions, the reflected wave and/or the transmitted wave will maintain a beam profile. We study the interface conditions which specify the relations between the split waves and the incident Gaussian beam. Numerical tests are presented to validate the beam formulations and interface conditions.     

A stable high-order method for the heated cavity problem

A fourth-order method, without using extrapolation, is developed for the steady-state solution of a non-linear system of three simultaneous partial differential equations for the flow of a fluid in a heated closed cavity. The method is a finite difference method which has converged for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. The results are presented and compared with some of the accurate results available in de Vahl Davis and Jones, Shay and Schultz, and Dennis and Hudson. The method used to develop the fourth-order method presented in this paper can be used to develop high-order methods for other partial differential equations. The method was developed to be stable without using the upwinding technique.     

High-order local non-reflecting boundary conditions: a review

A common method for numerically solving wave problems in unbounded domains is based on truncating the infinite domain via an artificial boundary B, thus defining a finite computational domain, and using a special non-reflecting boundary condition (NRBC) on B. Low-order local NRBCs have been constructed and practiced since the 1970s. Exact non-local NRBCs were introduced in the 1980s. Only recently high-order local NRBCs have been devised. These NRBCs, despite being of an arbitrarily high-order, do not involve high derivatives owing to the use of specially defined auxiliary variables. This paper reviews the latter approach, explains its advantages compared to previous approaches, and discusses the different schemes which have been proposed in this context.     

一种易于并行求解Euler方程组的分区技术

在具有复杂边界的计算区域内,求解偏微分方程组时,经常需要分区和并行计算,分区方法直接关系到数值计算的并行化程度,本文在应用时间算子分裂方法求解Euler方程组的过程中,提出了一种非常容易实现并行化计算的分区技术。     

Use of characteristics method with cubic interpolation for unsteady-flow computation

The specified-time-interval (STI) scheme has been used commonly in applying the method of characteristics (MOC) to unsteady open-channel flow problems. However, with the use of STI scheme, the numerical error for the simulation results can always be induced due to the interpolation used to approximate the characteristics trajectory. Hence, in order to remedy the numerical errors caused by the interpolation, one needs to seek some kind of interpolation technique with higher-order accuracy. Instead of the linear interpolation technique, which has been used very commonly and can induce serious numerical diffusion, the Holly--Preissmann two-point, method, which is a cubic interpolation technique with fourth-order of accuracy, is proposed here to integrate with the method of characteristics for the computation of one-dimensional unsteady flow in open channel. The concept of reachback and reachout in space and time directions for the characteristics is also introduced to assure the model stability. The computed results from this new model are compared with those computed by using the Preissmann four-point scheme and the multimode method of characteristics with linear interpolation.     

Efficient high-order immersed interface methods for heat equations with interfaces

An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.     

Generalized stroh formalism for anisotropic elasticity for general boundary conditions

The Stroh formalism is most elegant when the boundary conditions are simple, namely, they are prescribed in terms of traction or displacement. For mixed boundary conditions such as there for a slippery boundary, the concise matrix expressions of the Stroh formalism are destroyed. We present a generalized Stroh formalism which is applicable to a class of general boundary conditions. The general boundary conditions include the simple and slippery boundary conditions as special cases. For Green's functions for the half space, the general solution is applicable to the case when the surface of the half-space is a fixed, a free, a slippery, or other more general boundary. For the Griffith crack in the infinite space, the crack can be a slit-like crack with free surfaces, a rigid line inclusion (which is sometimes called an anticrack), or a rigid line with slippery surface or with other general surface conditions. It is worth mention that the modifications required on the Stroh formalism are minor. The generalized formalism and the final solutions look very similar to those of unmodified version. Yet the results are applicable to a rather wide range of boundary conditions.     

Radiation boundary conditions for finite element solutions of generalized wave equations

On the basis of the dispersion relation of the generalized linear wave equation we derive a radiation boundary condition (RBC) that explicitly incorporates the physical parameters of the governing equation into the form of the boundary condition. Using finite element techniques we investigate the properties of the generalized RBC by examining forced and unforced solutions to the telegraph and Klein-Gordon equations in one dimension. The results show that within the limits of the physical parameters of the problem the generalized RBC is an improvement over the Sommerfeld RBC when the governing equation contains additional terms that influence the propagation. These gains are achieved without introducing any computational overhead. A two-dimensional example suggests that the 1D findings can generalize to higher dimensions.     

Holberg's optimisation for high-order compact finite difference staggered schemes

Two optimised high-order compact finite difference (FD) staggered schemes are presented in this communication. Following Holberg's optimisation strategy, the least squares problem to minimising the group velocity (MGV) error, for the fourth- and sixth-order pentadiagonal schemes, is formulated. For a fixed level of group velocity accuracy, the optimised spectrum of wave number and the optimised coefficients for the schemes, are analytically evaluated. The spectral accuracy of these schemes has been verified by several comparisons with the FD staggered schemes obtained following Kim and Lee's (1996) optimisation procedure. Fewer group and phase velocity errors, greater resolution in terms of absolute error and resolving efficiency have been achieved by the optimised schemes proposed. High-order accuracy in time is obtained by marching the solution with an optimised Runge–Kutta scheme. Next, the comparison in terms of the number of grid points per wavelength required to achieve a standard accuracy for distances expressed in terms of the number of wavelengths travelled is presented. Numerical results from benchmark tests for the one-dimensional shallow water equations are presented.     

Supercritical equilibrium solutions of axially moving beams with hybrid boundary conditions

In this paper supercritical equilibria and critical speeds of axially moving beams constrained by sleeves with torsion springs are deduced. Transverse vibration of the beams is governed by a nonlinear integro-partial-differential equation. In the supercritical regime, the corresponding static equilibrium equation for the hybrid boundary conditions is analytically solved for the equilibria and the critical speeds. In the view of the non-trivial equilibrium, comparisons are made among the integro-partial-differential equation, a nonlinear partial-differential equation for transverse vibration, and coupled equations for planar motion under the hybrid boundary conditions.     

A finite difference technique for solving the Newtonian jet swell problem

A finite difference technique that incorporates a numerical mapping has been successfully applied to analyse both planar and axisymmetric Newtonian jets. A pressure gradient equation and a free-surface slope equation have been derived for free-surface iteration. The computation of pressure inside the jet surface using the pressure gradient equation is stable and accurate at high Reynolds numbers. The free-surface slope equation is needed for updating the free surface and is applicable for jets with strong surface tension effects. The present development can simulate the Newtonian jets for Reynolds numbers as high as 2000 and capillary number as low as 10^?5. Numerical predictions by the present technique are close to the results of previous finite element simulations.     

Pseudo‐characteristic formulation and dynamic boundary conditions for computational aeroacoustics

The present paper deals with the use of the pseudo‐characteristic formulation of the Navier–Stokes and Euler equations recently introduced by Sesterhenn (Comput. Fluid. 2001; 30 :37–67) for the simulation of acoustic wave propagation. The emphasis is put on the formulation of an efficient method on structured curvilinear grids, along with the definition and implementation of efficient boundary conditions. The cases of inflow, outflow, rigid/compliant walls and walls with prescribed impedance are addressed. The proposed boundary conditions are assessed on generic cases. The pseudo‐characteristic formulation enables a straightforward and optimal use of high‐order upwind dispersion‐relation‐preserving schemes, yielding an efficient method. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical schemes with high order of accuracy for the computation of shock waves

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｏｍｐｕｔａｔｉｏｎｏｆｆｌｏｗｆｉｅｌｄｗｉｔｈｓｈｏｃｋｗａｖｅｓｈａｓｂｅｅｎｔｈｅｓｕｂｊｅｃｔｏｆｒｅｓｅａｒｃｈｆｏｒｍａｎｙｙｅａｒｓ.Ｔｈｅｒｅａｒｅｂａｓｉｃａｌｌｙｔｗｏｃａｔｅｇｏｒｉｅｓｏｆｍｅｔｈｏｄｓ,ｎａｍｅｌｙ,ｓｈｏｃｋｆｉｔｔｉｎｇｍｅｔｈｏｄａｎｄｓｈｏｃｋｃａｐｔｕｒｉｎｇｍｅｔｈｏｄ.Ｔｈｅｆｏｒｍｅｒｄｉｖｉｄｅｓｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｄｏｍａｉｎｉｎｔｏｓｕｂ＿ｄｏｍａｉｎｓｂｙｔｈｅｓｈｏｃｋｗａｖｅｓ.Ｉｎｅ…     

Analytical solutions for anisotropic bimaterials under several boundary conditions on the interface

Summary Analytical solutions are proposed for the stress and displacement fields induced by in-plane loading of a bimaterial under several boundary conditions. The two joined orthotropic layers forming the bimaterial are allowed to possess different types of anisotropy with their planes of elastic symmetry arbitrarily inclined with respect to the horizontal.     

实现有限元并行分布计算的一种新策略

在几种典型的计算网络上，给出了实现有限元并行分布计算的一种全新策略。它对子结构的划分方式没有任何限制，使结构划分方式对通讯不产生任何影响，并利用所谓的Σ通讯完成有关迭代计算。这种策略广泛适用于多项式加速法的并行分布迭代计算，使有限元并行分布计算的算法及程序与具体的计算网络有很好的分离性，同时也很大程度地保留了已有串行有限元算法及程序的优点。以预处理的共轭斜量法为例，在ＩｎｍｏｓＴ８００Ｔｒａｎｓｐｕｔｅｒ系统上实现了有限元并行分布计算。通过数值算例，验证了本文方法的可行性与有效性。     

Generalized impedance boundary conditions for thin dielectric coatings with variable thickness

We derive the so-called Generalized Impedance Boundary Conditions (GIBC) that model thin dielectric coatings with variable width located on a perfect electric or magnetic conducting surface. We treat the 2-D electromagnetic problem for both TM and TE polarizations. The expressions of the GIBCs are obtained up to the third order (with respect to the coating width). The order of convergence is numerically validated through various numerical examples using a finite element type method. Particular attention is given to the cases where the inner boundary has corner singularities.     

Differential-equation-based representation of truncation errors for accurate numerical simulation

High-order compact finite difference schemes for two-dimensional convection-diffusion-type differential equations with constant and variable convection coefficients are derived. The governing equations are employed to represent leading truncation terms, including cross-derivatives, making the overall O(h⁴) schemes conform to a 3 × 3 stencil. We show that the two-dimensional constant coefficient scheme collapses to the optimal scheme for the one-dimensional case wherein the finite difference equation yields nodally exact results. The two-dimensional schemes are tested against standard model problems, including a Navier-Stokes application. Results show that the two schemes are generally more accurate, on comparable grids, than O(h²) centred differencing and commonly used O(h) and O(h³) upwinding schemes.