Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground...We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.展开更多
A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a n...A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a new mitigation method is urgently expected to be explored.This study proposed a novel asymptotic linear method(ALM)for micro pressure wave(MPW)mitigation to achieve a constant gradient of initial c ompression waves(ICWs),via a study with various open ratios on hoods.The properties of ICWs and MPWs under various open ratios of hoods were analyzed.The results show that as the open ratio increases,the MPW amplitude at the tunnel exit initially decreases before rising.At the open ratio of 2.28%,the slope of the ICW curve is linearly coincident with a supposed straight line in the ALM,which further reduces the MPW amplitude by 26.9%at 20 m and 20.0%at 50 m from the exit,as compared to the unvented hood.Therefore,the proposed method effectively mitigates MPW and quickly determines the upper limit of alleviation for the MPW amplitude at a fixed train-tunnel operation condition.All achievements provide a ne w potential measure for the adaptive design of tunnel hoods.展开更多
Offshore structures are constantly subjected to the complex forces of the marine environment,including wind,sea waves,currents,and seismic loadings.Among these,wind and sea wave forces persist throughout the structure...Offshore structures are constantly subjected to the complex forces of the marine environment,including wind,sea waves,currents,and seismic loadings.Among these,wind and sea wave forces persist throughout the structure’s lifetime.This study proposes a dynamic analysis approach that incorporates both time and frequency domain methods to investigate the structural responses of offshore structures under the combined effects of wind and wave forces.A wind wave-pier coupling dynamic model is first developed using a small-scale single pier,with corresponding dynamic equilibrium equations established.Fluctuating wind and sea waves are simulated using the weighted amplitude wave superposition(WAWS)method and linear superposition,respectively.Wind and wave load histories are then derived via Fourier transforms.The structural dynamic responses under different loading scenarios(wind only,wave only,and combined wind and wave)are analyzed using the Newmarkβmethod.Additionally,the effects of varying wind and wave parameters on structural responses are evaluated.The simulation results demonstrate that the structural responses to wind-wave coupling are smaller than the superimposed effects of wind and wave forces acting independently.When wind speeds are relatively low,wave forces dominate structural displacement and serve as the primary source of vibration.展开更多
Using a recognition model of atmospheric gravity waves(AGWs),we identified 519 AGW events from the OH airglow images observed at the Dandong and Lhasa stations from 2015 to 2017.The 317 AGW events detected at the Dand...Using a recognition model of atmospheric gravity waves(AGWs),we identified 519 AGW events from the OH airglow images observed at the Dandong and Lhasa stations from 2015 to 2017.The 317 AGW events detected at the Dandong station have wavelengths ranging from 30 to 60 km,periods from 14 to 20 min,horizontal speeds from 30 to 60 m/s,and relative intensities from 0.4%to 0.6%,respectively.The parameters of 202 events recorded at the Lhasa station mainly vary within 15-35 km in horizontal wavelength,4-6 min in period,40-100 m/s in horizontal velocity,and 0.1%-0.3%in relative intensity.The occurrence rate peaks in winter and summer at Dandong and the peak in summer are absent at Lhasa because of the lack of convective weather.The seasonal propagation directions of the waves are influenced by both the wind field-filtering effect and the distribution of wave sources.In spring,because of the southeastward background wind field,fewer southeastward events are observed at the Dandong station.The situation at the Lhasa station is similar.In summer,both the Lhasa and Dandong stations are dominated by northeastward AGWs,which can be attributed to the southwestward wind.In autumn,ray-tracing results show that the events at Dandong mainly originate from wind shear,whereas the events at the Lhasa station are triggered by convective weather.The location of the wave sources determines the trend of the propagation directions at the Dandong and Lhasa stations in autumn.In winter,because of the eastward wind,more events are propagating to the southwest at the Dandong station.展开更多
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be effi...3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.展开更多
Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based ...Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.展开更多
A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations....A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.展开更多
A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element...A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.展开更多
The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6×6 Luttinger Kohn model.The effect of the number and period of pla...The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6×6 Luttinger Kohn model.The effect of the number and period of plane waves used for expansion on the stability of energy eigenvalues is examined.For practical calculation,it should choose the period large sufficiently to ensure the envelope functions vanish at the boundary and the number of plane waves large enough to ensure the energy eigenvalues keep unchanged within a prescribed range.展开更多
The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data...The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.展开更多
The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident...The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.展开更多
A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the l...A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.展开更多
A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated wate...A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.展开更多
In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident bou...In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.展开更多
A modified domain reduction method(MDRM) that introduces damping terms to the original DRM is presented in this paper. To verify the proposed MDRM and compare the computational accuracy of these two methods, a numeric...A modified domain reduction method(MDRM) that introduces damping terms to the original DRM is presented in this paper. To verify the proposed MDRM and compare the computational accuracy of these two methods, a numerical test is designed. The numerical results of the MDRM and DRM are compared using an extended meshed model. The results show that the MDRM significantly improved the computational accuracy of the DRM. Then, the MDRM is compared with two existing conventional methods, namely Liao's transmitting boundary and viscous-spring boundary with Liu's method. The MDRM shows its great advancement in computational accuracy, stability and range of applications. This paper also discusses the influence of boundary location on computational accuracy. It can be concluded that smaller models tend to have larger errors. By introducing two dimensionless parameters, φ_1 and φ_2, the rational distance between the observation point and the MDRM boundary is suggested. When φ_1 >2 or φ_2>13, the relative PGA error can be limited to 5%. In practice, the appropriate model size can be chosen based on these two parameters to achieve desired computational accuracy.展开更多
In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for th...In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for the corresponding problem is obtained.展开更多
A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedur...A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.展开更多
Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orien...Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orientations, etc. of the three array arrangements, i. e., linear array, T-type array and pentagon array, on the resolution of the directional spreading of waves, are investigated experimentally. This study can be used as a reference in the experimental study and the field measurement of directional waves.展开更多
A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid （VOF） method. An active ＇absorbing wave-maker based on linear wave theory is set on the left boun...A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid （VOF） method. An active ＇absorbing wave-maker based on linear wave theory is set on the left boundary of the wave flume. The progressive waves and the absorbing waves are generated simultaneously at the active wave generating-absorbing boundary. The absorbing waves are generated to eliminate the waves coming back to the generating boundary due to reflection from the outflow boundary and the structures. SIRW method proposed by Frigaard and Brorsen （1995） is used to separate the incident waves and reflected waves. The digital filters are designed based on the surface elevation signals of the two wave gauges. The corrected velocity of the wave-maker paddle is the output from the digital filter in real time. The numerical results of regular and irregular waves by the active absorbing-generating boundary are compared with the numerical results by the ordinary generating boundary to verify the performance of the active absorbing-generator boundary. The differences between the initial incident waves and the estimated incident waves are analyzed.展开更多
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
基金supported by the National Key R&D Program of China(No.2023YFA1606701)the National Natural Science Foundation of China(Nos.12175042,11890710,11890714,12047514,12147101,and 12347106)+1 种基金Guangdong Major Project of Basic and Applied Basic Research(No.2020B0301030008)China National Key R&D Program(No.2022YFA1602402).
文摘We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.
基金Project(24A0006)supported by the Key Project of Scientific Research Fund of Hunan Provincial Department of Education,ChinaProject(2024JJ5430)supported by the Natural Science Foundation of Hunan Province,ChinaProjects(2024JK2045,2023RC3061)supported by the Science and Technology Innovation Program of Hunan Province,China。
文摘A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a new mitigation method is urgently expected to be explored.This study proposed a novel asymptotic linear method(ALM)for micro pressure wave(MPW)mitigation to achieve a constant gradient of initial c ompression waves(ICWs),via a study with various open ratios on hoods.The properties of ICWs and MPWs under various open ratios of hoods were analyzed.The results show that as the open ratio increases,the MPW amplitude at the tunnel exit initially decreases before rising.At the open ratio of 2.28%,the slope of the ICW curve is linearly coincident with a supposed straight line in the ALM,which further reduces the MPW amplitude by 26.9%at 20 m and 20.0%at 50 m from the exit,as compared to the unvented hood.Therefore,the proposed method effectively mitigates MPW and quickly determines the upper limit of alleviation for the MPW amplitude at a fixed train-tunnel operation condition.All achievements provide a ne w potential measure for the adaptive design of tunnel hoods.
基金Project(2022YFB2302700)supported by the National Key Research and Development Program of China。
文摘Offshore structures are constantly subjected to the complex forces of the marine environment,including wind,sea waves,currents,and seismic loadings.Among these,wind and sea wave forces persist throughout the structure’s lifetime.This study proposes a dynamic analysis approach that incorporates both time and frequency domain methods to investigate the structural responses of offshore structures under the combined effects of wind and wave forces.A wind wave-pier coupling dynamic model is first developed using a small-scale single pier,with corresponding dynamic equilibrium equations established.Fluctuating wind and sea waves are simulated using the weighted amplitude wave superposition(WAWS)method and linear superposition,respectively.Wind and wave load histories are then derived via Fourier transforms.The structural dynamic responses under different loading scenarios(wind only,wave only,and combined wind and wave)are analyzed using the Newmarkβmethod.Additionally,the effects of varying wind and wave parameters on structural responses are evaluated.The simulation results demonstrate that the structural responses to wind-wave coupling are smaller than the superimposed effects of wind and wave forces acting independently.When wind speeds are relatively low,wave forces dominate structural displacement and serve as the primary source of vibration.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFF0711402)the Specialized Research Fund for State Key Laboratories。
文摘Using a recognition model of atmospheric gravity waves(AGWs),we identified 519 AGW events from the OH airglow images observed at the Dandong and Lhasa stations from 2015 to 2017.The 317 AGW events detected at the Dandong station have wavelengths ranging from 30 to 60 km,periods from 14 to 20 min,horizontal speeds from 30 to 60 m/s,and relative intensities from 0.4%to 0.6%,respectively.The parameters of 202 events recorded at the Lhasa station mainly vary within 15-35 km in horizontal wavelength,4-6 min in period,40-100 m/s in horizontal velocity,and 0.1%-0.3%in relative intensity.The occurrence rate peaks in winter and summer at Dandong and the peak in summer are absent at Lhasa because of the lack of convective weather.The seasonal propagation directions of the waves are influenced by both the wind field-filtering effect and the distribution of wave sources.In spring,because of the southeastward background wind field,fewer southeastward events are observed at the Dandong station.The situation at the Lhasa station is similar.In summer,both the Lhasa and Dandong stations are dominated by northeastward AGWs,which can be attributed to the southwestward wind.In autumn,ray-tracing results show that the events at Dandong mainly originate from wind shear,whereas the events at the Lhasa station are triggered by convective weather.The location of the wave sources determines the trend of the propagation directions at the Dandong and Lhasa stations in autumn.In winter,because of the eastward wind,more events are propagating to the southwest at the Dandong station.
基金supported by NSFC(Nos.41274120,41404085,and 41504084)
文摘3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.
基金supported by the China State Key Science and Technology Project on Marine Carbonate Reservoir Characterization (No. 2011ZX05004-003)the Basic Research Programs of CNPC during the 12th Five-Year Plan Period (NO.2011A-3603)+1 种基金the Natural Science Foundation of China (No.41104066)the RIPED Young Professional Innovation Fund (NO.2010-13-16-02, 2010-A-26-02)
文摘Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.
基金Supported by the National Natural Science Foundation of China(50976072,51106099,10902070)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)the Science Foundation for the Excellent Youth Scholar of Higher Education of Shanghai(slg09003)~~
文摘A numerical simulation of the toroidal shock wave focusing in a co-axial cylindrical shock tube is inves- tigated by using discontinuous Galerkin （DG） finite element method to solve the axisymmetric Euler equations. For validating the numerical method, the shock-tube problem with exact solution is computed, and the computed results agree well with the exact cases. Then, several cases with higher incident Mach numbers varying from 2.0 to 5.0 are simulated. Simulation results show that complicated flow-field structures of toroidal shock wave diffraction, reflection, and focusing in a co-axial cylindrical shock tube can be obtained at different incident Mach numbers and the numerical solutions appear steep gradients near the focusing point, which illustrates the DG method has higher accuracy and better resolution near the discontinuous point. Moreover, the focusing peak pres- sure with different grid scales is compared.
基金Project supported by the National Basic Research Program of China (973Project) (No.2002CB412709) and the National Natural Science Foundation of China (Nos.50278012,10272027,19832010)
文摘A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.
文摘The valence subband energies and wave functions of a tensile strained quantum well are calculated by the plane wave expansion method within the 6×6 Luttinger Kohn model.The effect of the number and period of plane waves used for expansion on the stability of energy eigenvalues is examined.For practical calculation,it should choose the period large sufficiently to ensure the envelope functions vanish at the boundary and the number of plane waves large enough to ensure the energy eigenvalues keep unchanged within a prescribed range.
基金supported by the National Science and Technology Major Project of China(Grant No. 2011ZX05004-003,2011ZX05014-006-006)the National Key Basic Research Program of China(Grant No. 2013CB228602)the Natural Science Foundation of China(Grant No. 40974066)
文摘The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.
基金National Natural Science Foundation of China under Grant No.51478247National Key Research and Development Program of China under Grant No.2016YFC1402800
文摘The method of inputting the seismic wave determines the accuracy of the simulation of soil-structure dynamic interaction. The wave method is a commonly used approach for seismic wave input, which converts the incident wave into equivalent loads on the cutoff boundaries. The wave method has high precision, but the implementation is complicated, especially for three-dimensional models. By deducing another form of equivalent input seismic loads in the fi nite element model, a new seismic wave input method is proposed. In the new method, by imposing the displacements of the free wave fi eld on the nodes of the substructure composed of elements that contain artifi cial boundaries, the equivalent input seismic loads are obtained through dynamic analysis of the substructure. Subsequently, the equivalent input seismic loads are imposed on the artifi cial boundary nodes to complete the seismic wave input and perform seismic analysis of the soil-structure dynamic interaction model. Compared with the wave method, the new method is simplifi ed by avoiding the complex processes of calculating the equivalent input seismic loads. The validity of the new method is verifi ed by the dynamic analysis numerical examples of the homogeneous and layered half space under vertical and oblique incident seismic waves.
基金the National Natural Science Foundation of China(50478014)the National 973 Program(2007CB714200)the Beijing Natural Science Foundation(8061003).
文摘A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.
基金The National Natural Science Foundation of China under contract No.51279023the Public Science and Technology Research Funds Projects of Ocean under contract No.201205023+1 种基金the Special Funds for Postdoctoral Innovative Projects of Liaoning Province of China under contract No.2011921018the Special Funds for Talent Projects of Dalian Ocean University under contract No.SYYJ2011004
文摘A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51309050 and 51221961)the National Basic Research Program of China(973 Program,Grant Nos.2013CB036101 and 2011CB013703)
文摘In this paper, a numerical model is developed based on the High Order Spectral （HOS） method with a non-periodic boundary. A wave maker boundary condition is introduced to simulate wave generation at the incident boundary in the HOS method. Based on the numerical model, the effects of wave parameters, such as the assumed focused amplitude, the central frequency, the frequency bandwidth, the wave amplitude distribution and the directional spreading on the surface elevation of the focused wave, the maximum generated wave crest, and the shifting of the focusing point, are numerically investigated. Especially, the effects of the wave directionality on the focused wave properties are emphasized. The numerical results show that the shifting of the focusing point and the maximum crest of the wave group are dependent on the amplitude of the focused wave, the central frequency, and the wave amplitude distribution type. The wave directionality has a definite effect on multidirectional focused waves. Generally, it can even out the difference between the simulated wave amplitude and the amplitude expected from theory and reduce the shifting of the focusing points, implying that the higher order interaction has an influence on wave focusing, especially for 2D wave. In 3D wave groups, a broader directional spreading weakens the higher nonlinear interactions.
基金National Natural Science Foundation of China under Grant Nos.91315301,51478279the State Key Laboratory Basic Theory Foundation of the Ministry of Science and Technology of China under the Grant SLDRCE08-A-07
文摘A modified domain reduction method(MDRM) that introduces damping terms to the original DRM is presented in this paper. To verify the proposed MDRM and compare the computational accuracy of these two methods, a numerical test is designed. The numerical results of the MDRM and DRM are compared using an extended meshed model. The results show that the MDRM significantly improved the computational accuracy of the DRM. Then, the MDRM is compared with two existing conventional methods, namely Liao's transmitting boundary and viscous-spring boundary with Liu's method. The MDRM shows its great advancement in computational accuracy, stability and range of applications. This paper also discusses the influence of boundary location on computational accuracy. It can be concluded that smaller models tend to have larger errors. By introducing two dimensionless parameters, φ_1 and φ_2, the rational distance between the observation point and the MDRM boundary is suggested. When φ_1 >2 or φ_2>13, the relative PGA error can be limited to 5%. In practice, the appropriate model size can be chosen based on these two parameters to achieve desired computational accuracy.
基金Supported by the National Natural Sciences Foundation of China(40676016 and 10471039)the National Key Project for Basic Research(2003CB415101-03 and 2004CB418304)+2 种基金the Key Project of the Chinese Academy of Sciences(KZCX3-SW-221)in part by E-Institutes of Shanghai Municipal Education Commission(N.E03004)the Natural Science Foundation of Zeijiang,China(Y606268).
文摘In this article the travelling wave solution for a class of nonlinear reaction diffusion problems are considered. Using the homotopic method and the theory of travelling wave transform, the approximate solution for the corresponding problem is obtained.
基金The project was financially supported by the National Natural Science Foundation of China under the Grant No. 19732004 the National Science Fund for Distinguished Young Scholars under the Grant No. 50029002
文摘A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.
文摘Measuring multi-directional waves with the wave gauge array is one of the fundamental and easily realised methods. In this paper, the wave gauge array is described and the effects of the gauge spacing, the array orientations, etc. of the three array arrangements, i. e., linear array, T-type array and pentagon array, on the resolution of the directional spreading of waves, are investigated experimentally. This study can be used as a reference in the experimental study and the field measurement of directional waves.
基金supported by the‘New Century Outstanding Talent’Scheme of the Ministry of Education of China(Grant No.NCET-07-0135)
文摘A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid （VOF） method. An active ＇absorbing wave-maker based on linear wave theory is set on the left boundary of the wave flume. The progressive waves and the absorbing waves are generated simultaneously at the active wave generating-absorbing boundary. The absorbing waves are generated to eliminate the waves coming back to the generating boundary due to reflection from the outflow boundary and the structures. SIRW method proposed by Frigaard and Brorsen （1995） is used to separate the incident waves and reflected waves. The digital filters are designed based on the surface elevation signals of the two wave gauges. The corrected velocity of the wave-maker paddle is the output from the digital filter in real time. The numerical results of regular and irregular waves by the active absorbing-generating boundary are compared with the numerical results by the ordinary generating boundary to verify the performance of the active absorbing-generator boundary. The differences between the initial incident waves and the estimated incident waves are analyzed.
