Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be o...Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be one of highly accurate methods. In this paper, we propose a new non-reflecting recursive algorithm for reverse-time migration by introducing the wave impedance function into the acoustic wave equation and the algorithm for the surface relief case is derived from the coordinate transformation principle. Using the exploding reflector principle and the zero-time imaging condition of poststack reverse- time migration, poststack numerical simulation and reverse-time migration with complex conditions can be realized. The results of synthetic and real data calculations show that the method effectively suppresses unwanted internal reflections and also deals with the seismic imaging problems resulting from surface relief. So, we prove that this method has strong adaptability and practicality.展开更多
The responses of the pavement in service are the basis for the design of the semi-rigid base course asphalt pavement. Due to the dynamic characteristics of wheel loacis and the temperature loads, the dynamic response ...The responses of the pavement in service are the basis for the design of the semi-rigid base course asphalt pavement. Due to the dynamic characteristics of wheel loacis and the temperature loads, the dynamic response analysis is very significant. In this article, the dynamic analysis of asphalt pavement under moving wheel loads is carried out using finite dement method canpled with non-reflective boundary method. The influences of the base modulus, thickness, the vehicle velocity, the tire pressure, and the contact condition at the interface are studied using parametric analysis. The results of numerical analysis show that it is not appropriate to simply increase the base modulus or thickness in the design. It would be beneficial if the base design is optimized synthetically. The increase of damping is also beneficial to the pavements because of the surface deflection and the stresses declination. Furthermore, the good contact condition at the interface results in good performance because it combines every layer of the pavement to work together. As overload aggravates the working condition of the pavement, it is not allowed.展开更多
Computational aeroacoustics （CAA） is an interdiscipline of aeroacoustics and computational fluid dynamics （CFD） for the investigation of sound generation and propagation from various aeroacoustics problems. In thi...Computational aeroacoustics （CAA） is an interdiscipline of aeroacoustics and computational fluid dynamics （CFD） for the investigation of sound generation and propagation from various aeroacoustics problems. In this review, the foundation and research scope of CAA are introduced firstly. A review of the early advances and applications of CAA is then briefly surveyed, focusing on two key issues, namely, high order finite difference scheme and non-reflecting boundary condition. Furthermore, the advances of CAA during the past five years are highlighted. Finally, the future prospective of CAA is briefly discussed.展开更多
In this paper, a new kind of discrete non-reflecting boundary conditions is developed. It can be used for a variety of wave equations such as the acoustic wave equation, the isotropic and anisotropic elastic wave equa...In this paper, a new kind of discrete non-reflecting boundary conditions is developed. It can be used for a variety of wave equations such as the acoustic wave equation, the isotropic and anisotropic elastic wave equations and the equations for wave propagation in multi-phase media and so on. In this kind of boundary conditions, the composition of all artificial reflected waves, but not the individual reflected ones, is considered and eliminated. Thus, it has a uniform formula for different wave equations. The velocity CA of the composed reflected wave is determined in the way to make the reflection coefficients minimal, the value of which depends on equations. In this paper, the construction of the boundary conditions is illustrated and CA is found, numerical results are presented to illustrate the effectiveness of the boundary conditions.展开更多
The hydrodynamic problem of a two-dimensional wedge entering water is studied based on Smoothed Particle Hydrodynamics （SPH） model. A non-reflection boundary treatment for SPH method is proposed to reduce the reflec...The hydrodynamic problem of a two-dimensional wedge entering water is studied based on Smoothed Particle Hydrodynamics （SPH） model. A non-reflection boundary treatment for SPH method is proposed to reduce the reflection of sound waves. The boundary pressure is obtained using an improved coupling boundary treatment approach, which is validated by comparing the simulation results with experimental and analytical results in literature. A series of cases with different initial entering velocities are simulated. The maximum force on the wedge and the corresponding time required to reach it for the different cases of initial entering velocities of the wedge are obained and fitted into formulas against the initial entering velocity of the wedge. The maximum drag coefficients of the wedge for the different cases with Froude number greater than 2 are all near the value of 0.91.展开更多
The multi-transmitting formula (MTF) governed by a single artificial speed is analytically developed into a generalized MTF governed by a few artificial speeds to improve its capacity in simultaneous simulation of sev...The multi-transmitting formula (MTF) governed by a single artificial speed is analytically developed into a generalized MTF governed by a few artificial speeds to improve its capacity in simultaneous simulation of several one-way waves propagating at different speeds.The generalized MTF is then discretized and further generalized using the space extrapolation to improve its accuracies in numerical simulation of transient waves at large angles of incidence.The above two successive generalizitions of MTF based on the notion of normal transmission lead to a compact formula of local non-reflecting boundary condition.The formula not only provides a general representation of the major schemes of existing local boundary conditions but can be used to generate new schemes,which combine advantages of different schemes.展开更多
We present the development of a non-reflecting boundary condition,based on the Local One-Dimensional Inviscid(LODI)approach,for Lattice Boltzmann Models working with multi-speed stencils.We test and evaluate the LODI ...We present the development of a non-reflecting boundary condition,based on the Local One-Dimensional Inviscid(LODI)approach,for Lattice Boltzmann Models working with multi-speed stencils.We test and evaluate the LODI implementation with numerical benchmarks,showing significant accuracy gains with respect to the results produced by a simple zerogradient condition.We also implement a simplified approach,which allows handling the unknown distribution functions spanning several layers of nodes in a unified way,still preserving a comparable level of accuracy with respect to the standard formulation.展开更多
A new weak boundary procedure for hyperbolic problems is presented.We consider high order finite difference operators of summation-by-parts form with weak boundary conditions and generalize that technique.The new boun...A new weak boundary procedure for hyperbolic problems is presented.We consider high order finite difference operators of summation-by-parts form with weak boundary conditions and generalize that technique.The new boundary procedure is applied near boundaries in an extended domain where data is known.We show how to raise the order of accuracy of the scheme,how to modify the spectrum of the resulting operator and how to construct non-reflecting properties at the boundaries.The new boundary procedure is cheap,easy to implement and suitable for all numerical methods,not only finite difference methods,that employ weak boundary conditions.Numerical results that corroborate the analysis are presented.展开更多
We present a model for simulating wave propagation in stratified magnetoatmospheres. The model is based on equations of ideal MHD together with gravitational source terms. In addition, we present suitable boundary dat...We present a model for simulating wave propagation in stratified magnetoatmospheres. The model is based on equations of ideal MHD together with gravitational source terms. In addition, we present suitable boundary data and steady statesto model wave propagation. A finite volume framework is developed to simulate thewaves. The framework is based on HLL and Roe type approximate Riemann solversfor numerical fluxes, a positivity preserving fractional steps method for discretizingthe source and modified characteristic and Neumann type numerical boundary conditions. Second-order spatial and temporal accuracy is obtained by using an ENO piecewise linear reconstruction and a stability preserving Runge-Kutta method respectively.The boundary closures are suitably modified to ensure mass balance. The numericalframework is tested on a variety of test problems both for hydrodynamic as well asmagnetohydrodynamic configurations. It is observed that only suitable choices of HLLsolvers for the numerical fluxes and balanced Neumann type boundary closures yieldstable results for numerical wave propagation in the presence of complex magneticfields.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 40974073)the National 863 Program (Grant No.2007AA060504)the National 973 Program (Grant No. 2007CB209605) and CNPC Geophysical Laboratories
文摘Presently the research based on the accurate seismic imaging methods for surface relief, complex structure, and complicated velocity distributions is of great significance. Reverse-time migration is considered to be one of highly accurate methods. In this paper, we propose a new non-reflecting recursive algorithm for reverse-time migration by introducing the wave impedance function into the acoustic wave equation and the algorithm for the surface relief case is derived from the coordinate transformation principle. Using the exploding reflector principle and the zero-time imaging condition of poststack reverse- time migration, poststack numerical simulation and reverse-time migration with complex conditions can be realized. The results of synthetic and real data calculations show that the method effectively suppresses unwanted internal reflections and also deals with the seismic imaging problems resulting from surface relief. So, we prove that this method has strong adaptability and practicality.
基金This paper is supported by the Programfor New Century Excellent Talents in University(2004)and Highway Depart ment of Henan Province
文摘The responses of the pavement in service are the basis for the design of the semi-rigid base course asphalt pavement. Due to the dynamic characteristics of wheel loacis and the temperature loads, the dynamic response analysis is very significant. In this article, the dynamic analysis of asphalt pavement under moving wheel loads is carried out using finite dement method canpled with non-reflective boundary method. The influences of the base modulus, thickness, the vehicle velocity, the tire pressure, and the contact condition at the interface are studied using parametric analysis. The results of numerical analysis show that it is not appropriate to simply increase the base modulus or thickness in the design. It would be beneficial if the base design is optimized synthetically. The increase of damping is also beneficial to the pavements because of the surface deflection and the stresses declination. Furthermore, the good contact condition at the interface results in good performance because it combines every layer of the pavement to work together. As overload aggravates the working condition of the pavement, it is not allowed.
基金Project supported by the National Basic Research Program of China(No.2012CB720202)the National Natural Science Foundation of China(No.51476005)the 111 Project of China(No.B07009)
文摘Computational aeroacoustics （CAA） is an interdiscipline of aeroacoustics and computational fluid dynamics （CFD） for the investigation of sound generation and propagation from various aeroacoustics problems. In this review, the foundation and research scope of CAA are introduced firstly. A review of the early advances and applications of CAA is then briefly surveyed, focusing on two key issues, namely, high order finite difference scheme and non-reflecting boundary condition. Furthermore, the advances of CAA during the past five years are highlighted. Finally, the future prospective of CAA is briefly discussed.
基金Supported by the National Natural Sciences Foundation of China and State Key Laboratory of Scientific and Engineering Computing.
文摘In this paper, a new kind of discrete non-reflecting boundary conditions is developed. It can be used for a variety of wave equations such as the acoustic wave equation, the isotropic and anisotropic elastic wave equations and the equations for wave propagation in multi-phase media and so on. In this kind of boundary conditions, the composition of all artificial reflected waves, but not the individual reflected ones, is considered and eliminated. Thus, it has a uniform formula for different wave equations. The velocity CA of the composed reflected wave is determined in the way to make the reflection coefficients minimal, the value of which depends on equations. In this paper, the construction of the boundary conditions is illustrated and CA is found, numerical results are presented to illustrate the effectiveness of the boundary conditions.
基金supported by the National Natural Science Foundation of China (Grant No.10802049)the Shanghai Leading Academic Discipline Project (Grant No.B206)the Research Found of State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University (Grant No.GKZD010806)
文摘The hydrodynamic problem of a two-dimensional wedge entering water is studied based on Smoothed Particle Hydrodynamics （SPH） model. A non-reflection boundary treatment for SPH method is proposed to reduce the reflection of sound waves. The boundary pressure is obtained using an improved coupling boundary treatment approach, which is validated by comparing the simulation results with experimental and analytical results in literature. A series of cases with different initial entering velocities are simulated. The maximum force on the wedge and the corresponding time required to reach it for the different cases of initial entering velocities of the wedge are obained and fitted into formulas against the initial entering velocity of the wedge. The maximum drag coefficients of the wedge for the different cases with Froude number greater than 2 are all near the value of 0.91.
基金Project supported by the Joint Seismological Science Foundation.
文摘The multi-transmitting formula (MTF) governed by a single artificial speed is analytically developed into a generalized MTF governed by a few artificial speeds to improve its capacity in simultaneous simulation of several one-way waves propagating at different speeds.The generalized MTF is then discretized and further generalized using the space extrapolation to improve its accuracies in numerical simulation of transient waves at large angles of incidence.The above two successive generalizitions of MTF based on the notion of normal transmission lead to a compact formula of local non-reflecting boundary condition.The formula not only provides a general representation of the major schemes of existing local boundary conditions but can be used to generate new schemes,which combine advantages of different schemes.
文摘We present the development of a non-reflecting boundary condition,based on the Local One-Dimensional Inviscid(LODI)approach,for Lattice Boltzmann Models working with multi-speed stencils.We test and evaluate the LODI implementation with numerical benchmarks,showing significant accuracy gains with respect to the results produced by a simple zerogradient condition.We also implement a simplified approach,which allows handling the unknown distribution functions spanning several layers of nodes in a unified way,still preserving a comparable level of accuracy with respect to the standard formulation.
基金supported by the National Science Foundation under Award No.0948304 and by the Southern California Earthquake Center.SCEC is funded by NSF Cooperative Agreement EAR-0529922 and USGS Cooperative Agreement 07HQAG0008(SCEC contribution number 1806).The work by the last author was carried out within the Swedish e-science Research Centre(SeRC)and supported by the Swedish Research Council(VR).
文摘A new weak boundary procedure for hyperbolic problems is presented.We consider high order finite difference operators of summation-by-parts form with weak boundary conditions and generalize that technique.The new boundary procedure is applied near boundaries in an extended domain where data is known.We show how to raise the order of accuracy of the scheme,how to modify the spectrum of the resulting operator and how to construct non-reflecting properties at the boundaries.The new boundary procedure is cheap,easy to implement and suitable for all numerical methods,not only finite difference methods,that employ weak boundary conditions.Numerical results that corroborate the analysis are presented.
文摘We present a model for simulating wave propagation in stratified magnetoatmospheres. The model is based on equations of ideal MHD together with gravitational source terms. In addition, we present suitable boundary data and steady statesto model wave propagation. A finite volume framework is developed to simulate thewaves. The framework is based on HLL and Roe type approximate Riemann solversfor numerical fluxes, a positivity preserving fractional steps method for discretizingthe source and modified characteristic and Neumann type numerical boundary conditions. Second-order spatial and temporal accuracy is obtained by using an ENO piecewise linear reconstruction and a stability preserving Runge-Kutta method respectively.The boundary closures are suitably modified to ensure mass balance. The numericalframework is tested on a variety of test problems both for hydrodynamic as well asmagnetohydrodynamic configurations. It is observed that only suitable choices of HLLsolvers for the numerical fluxes and balanced Neumann type boundary closures yieldstable results for numerical wave propagation in the presence of complex magneticfields.
