An adaptive open boundary condition

A significant improvement of the open boundary condition which was originally studied by Engquist and Majda is given. The method given here is applicable without limitation of the angle between the direction of the incident wave and the normal to the boundary. It also improves the precision of the method from order one to two. The test examples show that this method is much better than the method mentioned above. Inst. of Software, Academia Sinica     

Finite element method for shallow water equation including open boundary condition

A method to deal with an open boundary condition in the analysis of water surface waves, the tide, etc. by means of the finite element method is proposed in this paper. The present method has two important features relating to the treatment of the open boundary condition. The first feature is to consider the non-reflective virtual boundary condition which has been developed in the numerical wave analysis method. The incident wave conditions without spurious reflected waves can be imposed at the open boundary. The second feature is to identify the amplitude of the components of incident waves in terms of observed water elevations in the field of standing waves. This can be done as a parameter identification based on an optimization problem by applying the conjugate gradient method. The applicability of this method to wave propagation problems is verified by several numerical computations.     

Energy flux equalization: An open boundary condition for non-linear free surface flows

In the simulation of non-linear free surface flows in a finite domain a major concern is with the radiation condition to be applied at the ‘open’ boundary. No general theoretical radiation conditions are known to exist. In this paper a new open boundary condition is formulated based on energy flux equalization between the non-linear inner domain and a linear outer domain. The non-linear flow in the inner domain is solved numerically using a semi-Lagrangian procedure. The energy flux arriving at the open boundary is removed using a new open boundary condition which acts as a linear wave absorber. For the cases studied the performance of this boundary condition is found to be quite good.     

基于吸收边界条件的有限元方法和离散切比雪夫展开的二维缺陷识别

结合考虑了吸收边界条件的有限元正演方法,通过离散切比雪夫多项式在特定离散点的展开逼近待识别函数,实现了无限区域内二维缺陷的识别。正问题计算中引入的吸收边界条件虽然不能严格地模拟无限域对近场波动的影响,但能够以满足实际需要的精度模拟这一影响,且具有解耦特性,极大地简化了数值计算,在反演计算中,离散点取在切比雪夫多项式的零点处,使逼近时最大偏差最小,数值算例表明此方法的有效性。     

Non-reflecting boundary conditions applicable to general purpose CFD simulators

In simulations of propagating blast waves the effects of artificial reflections at open boundaries can seriously degrade the accuracy of the computations. In this paper, a boundary condition based on a local approximation by a plane traveling wave is presented. The method yields small artificial reflections at open boundaries. The derivation and the theory behind these so-called plane-wave boundary conditions are presented. The method is conceptually simple and is easy to implement in two and three dimensions. These non-reflecting boundary conditions are employed in the three-dimensional computational fluid dynamics (CFD) solver FLACS, capable of simulating gas explosions and blast-wave propagation in complex geometries. Several examples involving propagating waves in one and two dimensions, shock tube and an example of a simulation of a propagating blast wave generated by an explosion in a compressor module are shown. The numerical simulations show that artificial reflections due to the boundary conditions employed are negligible. © 1998 John Wiley & Sons, Ltd.     

Discrete non‐local absorbing boundary condition for exterior problems governed by Helmholtz equation

The finite element method is employed to approximate the solutions of the Helmholtz equation for water wave radiation and scattering in an unbounded domain. A discrete, non‐local and non‐reflecting boundary condition is specified at an artificial external boundary by the DNL method, yielding an equivalent problem that is solved in a bounded domain. This procedure formulates a boundary value problem in a bounded region by imposing a relation in the discrete medium between the nodal values at the two last layers. For plane geometry, this relation can be found by straightforward eigenvalue decomposition. For circular geometry, the plane condition is applied at the external layer and this condition is condensed through a structured annular region, resulting in a condition at an inner radius. Exterior problems with a bounded internal physical obstacle are considered. It is well‐known that these kind of problems are well‐posed, and have a unique solution. Numerical studies based on standard Galerkin methodology examine the dependence of the DNL condition with respect to the circular annular region width. The DNL condition is compared with local boundary conditions of several orders. Numerical examples confirm the important improvement in accuracy obtained by the DNL method over standard conditions. Copyright © 1999 John Wiley & Sons, Ltd.     

An immersed boundary method for simulation of inviscid compressible flows

In this paper, an immersed boundary method for simulating inviscid compressible flows governed by Euler equations is presented. All the mesh points are classified as interior computed points, immersed boundary points (interior points closest to the solid boundary), and exterior points that are blanked out of computation. The flow variables at an immersed boundary point are determined via the approximate form of solution in the direction normal to the wall boundary. The normal velocity is evaluated by applying the no‐penetration boundary condition, and therefore, the influence of solid wall in the inviscid flow is taken into account. The pressure is computed with the local simplified momentum equation, and the density and the tangential velocity are evaluated by using the constant‐entropy relation and the constant‐total‐enthalpy relation, respectively. With a local coordinate system, the present method has been extended easily to the three‐dimensional case. The present work is the first endeavor to extend the idea of hybrid Cartesian/immersed boundary approach to compressible inviscid flows. The tedious task of handling multi‐valued points can be eliminated, and the overshoot resulting from the extrapolation for the evaluation of flow variables at exterior points can also be avoided. In order to validate the present method, inviscid compressible flows over fixed and moving bodies have been simulated. All the obtained numerical results show good agreement with available data in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

Interface damage modeled by spring boundary conditions for in-plane elastic waves

In-plane elastic wave propagation in the presence of a damaged interface is investigated. The damage is modeled as a distribution of small cracks and this is transformed into a spring boundary condition. First the scattering by a single interface crack is determined explicitly in the low frequency limit for the case of a plane wave normally incident to the interface. The transmission at an interface with a random distribution of small cracks is then determined and is compared to periodically distributed cracks. The cracked interface is then described by a distributed spring boundary condition. As an illustration the dispersion relation of the first modes in a thick plate with a damaged interface in the middle is given.     

Inverse scattering problem from an impedance crack via a composite method

In this paper we consider an inverse scattering problem from an open arc with impedance boundary conditions on both sides of the crack. Our aim is to recover both the impedance function and the unknown crack simultaneously from the far-field pattern with only one incident wave. Making the most out of the direct problem, a straightforward method of iterative nature is developed for the inverse problem. The ill-posedness of this problem is considered by incorporating the Tikhonov regularization. Numerical examples are provided at the end of the paper to show the feasibility of our method.     

Transmission boundary condition for mixed variable finite-difference method

All kinds of numerical methods based on the discretization of first-order velocity-stress hyperbolic system of elastic wave equations are named mixed variable finite-difference method in this paper, and the transmission boundary condition of artificial boundaries is studied here by the mixed variable finite-difference method. The transmission condition of complex geometrical boundaries of a transversely isotropic medium is presented based on the characteristic variables of wavefields propagating in the normal direction of the boundary. The boundary condition proposed in this paper is a local artificial boundary condition, with which the computation cost is very low. Elastic wave propagations in transversely isotropic medium are modelled by the staggered grid finite-difference method incorporated with the boundary condition presented. Numerical results and analysis of reflection coefficients show that the reflections of incoming waves by artificial boundaries are efficiently reduced.Supported by National Natural Science Foundation and Liaoning Province Science Foundation.