Exact boundary conditions for the initial value problem of convex conservation laws

The initial value problem of convex conservation laws, which includes the famous Burgers’ (inviscid) equation, plays an important rule not only in theoretical analysis for conservation laws, but also in numerical computations for various numerical methods. For example, the initial value problem of the Burgers’ equation is one of the most popular benchmarks in testing various numerical methods. But in all the numerical tests the initial data have to be assumed that they are either periodic or having a compact support, so that periodic boundary conditions at the periodic boundaries or two constant boundary conditions at two far apart spatial artificial boundaries can be used in practical computations. In this paper for the initial value problem with any initial data we propose exact boundary conditions at two spatial artificial boundaries, which contain a finite computational domain, by using the Lax’s exact formulas for the convex conservation laws. The well-posedness of the initial-boundary problem is discussed and the finite difference schemes applied to the artificial boundary problems are described. Numerical tests with the proposed artificial boundary conditions are carried out by using the Lax–Friedrichs monotone difference schemes.     

孤立波对水下山包的绕射及开路边界条件的应用比较

在无界域的数值方法中,开路边界条件(Open Boundary Condition(OBC))是一个十分重要的研究课题.对于波动问题已有许多不同形式的OBC。但对它们的有效性,许多还缺乏在直接应用中的分析比较。本文以非线性长波的Boussinesq方程为控制方程组,在方形人为边界上采用了几类较重要的OBC,计算了无界水域中孤立波对水下山包的绕射。文中既给出了令人感兴趣的三维散射图案,又对这些OBC的有效性作了分析比较,对它们的可用性提出了一些有参考价值的意见。     

双曲方程外问题的完全无反射人工边界条件

研究了数值求解双曲方程外问题的人工边界方法.在圆形的人工边界上,得到了三类等价的完全无反射的人工边界条件.给出的数值例子验证了方法的有效性.     

三维Helmholtz方程外问题的自然积分方程及其数值解

用文[2,3]提出的自然边界归化方法来处理三维Helmholtz方程的外边值问题。在简要介绍如何用球谐展开的方法得到Helmholtz问题在外球域上的自然积分方程后,给出求解该自然积分方程的一种数值方法及相应的数值算例。     

Transmission loss prediction of reactive silencers using 3-D time-domain CFD approach and plane wave decomposition technique

A predictive method is proposed to determine the transmission loss of reactive silencers using the three-dimensional (3-D) time-domain computational fluid dynamics (CFD) approach and the plane wave decomposition technique. Firstly, a steady flow computation is performed with a mass-flow-inlet boundary condition, which provides an initial condition for the following two unsteady flow computations. The first unsteady flow computation is conducted by imposing an impulse (acoustic excitation) superimposed on the constant mass flow at the inlet of the model and then adding the non-reflecting boundary condition (NRBC) when the impulse completely propagates into the silencer. The second unsteady flow computation is conducted for the case without acoustic excitation at the inlet. The time histories of pressure and velocity at the upstream monitoring point as well as history of pressure at the downstream monitoring point are recorded during the two transient computations. The differences between the two unsteady flow computational results are the corresponding acoustic quantities. Therefore, the incident sound pressure signal is obtained by using plane wave decomposition at upstream, while the transmitted sound pressure signal is just the sound pressure at downstream. Finally, those two sound pressure signals in the time-domain are transformed into the frequency-domain by Fast Fourier Transform (FFT) and then the transmission loss (TL) of silencer is determined. For the straight-through perforated tube silencers with and without flow, the numerical results agree well with the published measurements.     

Open and traction boundary conditions for the incompressible Navier–Stokes equations

We present numerical schemes for the incompressible Navier–Stokes equations (NSE) with open and traction boundary conditions. We use pressure Poisson equation (PPE) formulation and propose new boundary conditions for the pressure on the open or traction boundaries. After replacing the divergence free constraint by this pressure Poisson equation, we obtain an unconstrained NSE. For Stokes equation with open boundary condition on a simple domain, we prove unconditional stability of a first order semi-implicit scheme where the pressure is treated explicitly and hence is decoupled from the computation of velocity. Using either boundary condition, the schemes for the full NSE that treat both convection and pressure terms explicitly work well with various spatial discretizations including spectral collocation and C⁰ finite elements. Moreover, when Reynolds number is of O(1) and when the first order semi-implicit time stepping is used, time step size of O(1) is allowed in benchmark computations for the full NSE. Besides standard stability and accuracy check, various numerical results including flow over a backward facing step, flow past a cylinder and flow in a bifurcated tube are reported. Numerically we have observed that using PPE formulation enables us to use the velocity/pressure pairs that do not satisfy the standard inf–sup compatibility condition. Our results extend that of Johnston and Liu [H. Johnston, J.-G. Liu, Accurate, stable and efficient Navier–Stokes solvers based on explicit treatment of the pressure term. J. Comp. Phys. 199 (1) (2004) 221–259] which deals with no-slip boundary conditions only.     

Exact wave field simulation for finite-volume scattering problems

An exact boundary condition is presented for scattering problems involving spatially limited perturbations of arbitrary magnitude to a background model in generally inhomogeneous acoustic media. The boundary condition decouples the wave propagation on a perturbed domain while maintaining all interactions with the background model, thus eliminating the need to regenerate the wave field response on the full model. The method, which is explicit, relies on a Kirchhoff-type integral extrapolation to update the boundary condition at every time step of the simulation. The Green's functions required for extrapolation through the background model are computed efficiently using wave field interferometry.     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法,该方程求解声辐射和散射问题具有相同的影响矩阵,能够同时求解薄壳体气动和振动噪声的辐射和散射现象,以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上,可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

三维边界层内诱导横流失稳模态的感受性机理

边界层感受性问题是层流向湍流转捩的初始阶段,在转捩过程中起关键性作用,尤其是三维边界层流动.因此,研究三维边界层感受性问题对进一步理解层流向湍流转捩机理以及湍流成因具有重要的理论意义.采用数值方法研究自由来流湍流与三维壁面局部粗糙相互作用下三维边界层的感受性问题,确定是否能在三维边界层内寻找一种新的横流失稳模态;确定在何种条件下三维边界层内能诱导出定常、非定常的横流失稳模态;探索自由来流湍流的强度、展向波数和法向波数以及三维壁面局部粗糙的大小和结构类型等因素在自由来流湍流与三维壁面局部粗糙作用下三维边界层内被激发出的感受性过程中有何影响,并确定何种横流失稳模态在三维边界层感受性过程中占据何种地位.对自由来流湍流与三维壁面局部粗糙作用激发三维边界层内感受性问题的深入研究,将有助于完善流动稳定性与湍流理论,为层流向湍流转捩过程的预测与控制提供合理的理论依据.     

Boundary element-free method for elastodynamics

The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.     

Boundary element-free method for elastodynamics

1 Introduction In recent years, more and more attention has been paid to researches on the meshless (or meshfree) method, which makes it a hot direction of computational mechanics[1,2]. The meshless method is the approximation based on nodes, then the large deformation and crack growth problems can be simulated with the method without the re-meshing technique. And the meshless method has some advantages over the traditional computa- tional methods, such as finite element method (FEM) and boun…     

基于表面阻抗边界条件的等离子体薄涂层电磁散射的时域有限差分分析

基于表面阻抗边界条件时域有限差分(FDTD)方法研究了一维斜入射情况下非磁化等离子体薄涂层涂敷金属材料的电磁散射特性, 该方法忽略对薄层背景材料进行网格剖分, 大大减少了计算量. 首先推导了理想导体涂敷等离子体薄涂层的表面阻抗频域表达式, 然后代入边界条件并变换到时域, 再用分段线性递推卷积方法将时域表达式离散得到FDTD迭代式. 编程计算了垂直及斜入射情形下的平行极化和垂直极化反射系数, 通过验证算例与解析解对比, 结果表明该方法的准确性和有效性. 最后利用该方法分析了不同入射角对反射系数的影响.     

An Efficient Algorithm for Truncating Spatial Domain in Modeling Light Scattering by Finite-Difference Technique

The finite-difference time domain technique is one of the most robust and accurate numerical methods for the solution of light scattering by small particles with arbitrary composition and geometry. In practice, this method requires that the spatial domain for the computation of near-field be truncated. An absorbing boundary condition must be imposed in conjunction with this truncation. The performance of this boundary condition is essential to the stability of numerical computations and the reliability of results. In the present study, a new boundary condition, referred to as the mixed T algorithm, has been developed, which is a generalization of the transmitting boundary condition originally developed by Liao and co-workers. The present algorithm does not require spatial interpolation for wave values at interior grid points. In addition, it produces two minima of spurious reflections at small and large incident angles, allowing efficient absorption of the scattered waves at the boundary for large incident angles. When the third-order mixed T algorithm is used, the reflection coefficient of the boundary is less than 1% for incident angles from 0° to about 70°. We find that the numerical instability associated with the transmitting boundary condition is caused by the location-dependent amplitude of outgoing waves in the vicinity of the boundary. For this reason, the mixed T algorithm is stabilized by consistently introducing diffusive coefficients into the boundary equation. When the stabilized algorithm is applied, the near-field within the truncated domain can be computed by using single-precision arithmetic without overflows for more than 10⁵steps in the time-marching iteration. Finally, the new absorbing boundary condition is validated by carrying out numerical experiments involving the propagation of a TM wave excited by a sinusoidal point source, simultaneous simulation of the wave propagation in small and large domains, and the scattering of a TM wave by an infinite circular cylinder.     

Transparent boundary conditions for the Helmholtz equation in some ramified domains with a fractal boundary

The paper addresses a class of boundary value problems in some self-similar ramified domains, with the Laplace or Helmholtz equations. Much stress is placed on transparent boundary conditions which allow the solutions to be computed in subdomains. A self similar finite element method is proposed and tested. It can be used for numerically computing the spectrum of the Laplace operator with Neumann boundary conditions, as well as the eigenmodes. The eigenmodes are normalized by means of a perturbation method and the spectral decomposition of a compactly supported function is carried out. Finally, a numerical method for the wave equation is addressed.     

Modified periodic-shell boundary condition for dissipative particle dynamics simulations

We have developed the modified periodic-shell boundary condition (BC) for dissipative particle dynamics (DPD) simulations, which enables us to simulate an outer flow problem around an obstacle using a small simulation region. In order to clarify the validity of this BC, we have treated a uniform flow past a circular cylinder. The present BC has been compared with the ordinary BC such as the uniform flow condition. Also, the present results have been compared with those of the numerical results of the Navier–Stokes equation. The ordinary uniform BC is seen to give rise to significantly distorted flow fields and also to significant disappearance of dissipative particles from the simulation region. In contrast, for the present modified periodic-shell BC, the number density of dissipative particles is kept almost constant during a simulation run, and the flow field is in reasonable agreement with the result, which has been obtained by numerical simulations of the Navier–Stokes equation.     

The evaluation of non-reflecting boundary conditions for duct acoustic computation

A simple method is used to quantify the performance of non-reflecting boundary conditions for duct acoustic applications. The method uses a two-dimensional wavesplitting technique to decompose the linearized Euler equations into its constitutive modes, allowing the magnitude of reflected waves from outflow boundaries to be accurately determined. Realistic conditions are simulated by conducting the boundary condition analysis using acoustic waves with characteristics often found in duct and turbomachinery problems. For this paper, the method is used to investigate the performance of three different buffer zone implementations and the Perfectly Matched Layer as non-reflecting boundary conditions. The effect of the damping properties of these boundary conditions and the incident acoustic wave characteristics on performance are considered. Results indicate that the buffer zone boundary condition using explicit damping of the solution vector after each timestep produces the best non-reflecting performance. A deterioration in performance was observed for incident waves at high angles relative to the boundary normal for all implementations.     

Direct assessment of lattice Boltzmann hydrodynamics and boundary conditions for recirculating flows

A hydrodynamic boundary condition is developed for lattice Boltzmann hydrodynamics using a square, orthogonal grid. A constraint based on energy considerations is developed to provide closure for the equations which govern the particle distribution at the boundaries. This boundary condition is applied to the two-dimensional, steady flow of an incompressible fluid behind a grid, known as Kovasznay flow. The results are compared to those using alternate boundary conditions using the known exact solution. The hydrodynamic boundary condition produces quadratic spatial convergence, while alternate techniques fail to maintain this second-order accuracy.     

一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法

基于伸展坐标完全匹配层方程和辅助微分方程 方法, 给出了一种在时域有限差分(FDTD) 计算中适用于常见色散介质的通用边界条件算法. 该算法适用于任意阶的FDTD空间差分, 并且由于所采用的D-H方程独立于计算区域, 所以可以被用于截断任意电介质. 数值试验结果表明, 与卷积完全匹配层 算法相比较, 所提出的吸收边界条件算法不仅通用性强、 计算复杂度低、 计算时间短, 并且吸收效果有明显的提高.     

Probability of boundary conditions in quantum cosmology

One of the main interest in quantum cosmology is to determine boundary conditions for the wave function of the universe which can predict observational data of our universe. For this purpose, we solve the Wheeler–DeWitt equation for a closed universe with a scalar field numerically and evaluate probabilities for boundary conditions of the wave function of the universe. To impose boundary conditions of the wave function, we use exact solutions of the Wheeler–DeWitt equation with a constant scalar field potential. These exact solutions include wave functions with well known boundary condition proposals, the no-boundary proposal and the tunneling proposal. We specify the exact solutions by introducing two real parameters to discriminate boundary conditions, and obtain the probability for these parameters under the requirement of sufficient e-foldings of the inflation. The probability distribution of boundary conditions prefers the tunneling boundary condition to the no-boundary boundary condition. Furthermore, for large values of a model parameter related to the inflaton mass and the cosmological constant, the probability of boundary conditions selects an unique boundary condition different from the tunneling type.