The size and shape of the effective test area are crucial to consider when short-crested waves are created by segmented wavemakers. The range of the effective test area of short-crested waves simulated by two-sided se...The size and shape of the effective test area are crucial to consider when short-crested waves are created by segmented wavemakers. The range of the effective test area of short-crested waves simulated by two-sided segmented wavemakers is analyzed in this paper. The experimental investigation on the wave field distribution of short-crested waves generated by two-sided segmented wavemakers is conducted by using an array of wave gauges. Wave spectra and directional spreading function are analyzed and the results show that when the main direction is at a certain angle with the normal line of wave generators, the wave field of 3D short-crested waves generated by two-sided segmented wavemakers has good spatial uniformity within the model test area. The effective test area can provide good wave environments for seakeeping model tests of various ocean engineering structures in the deep ocean engineering basin.展开更多
The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric an...The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.展开更多
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.展开更多
Considering the serious electromagnetic wave(EMW)pollution problems and complex application condition,there is a pressing need to amalgamate multiple functionalities within a single substance.However,the effective int...Considering the serious electromagnetic wave(EMW)pollution problems and complex application condition,there is a pressing need to amalgamate multiple functionalities within a single substance.However,the effective integration of diverse functions into designed EMW absorption materials still faces the huge challenges.Herein,reduced graphene oxide/carbon foams(RGO/CFs)with two-dimensional/three-dimensional(2D/3D)van der Waals(vdWs)heterostructures were meticulously engineered and synthesized utilizing an efficient methodology involving freeze-drying,immersing absorption,secondary freeze-drying,followed by carbonization treatment.Thanks to their excellent linkage effect of amplified dielectric loss and optimized impedance matching,the designed 2D/3D RGO/CFs vdWs heterostructures demonstrated commendable EMW absorption performances,achieving a broad absorption bandwidth of 6.2 GHz and a reflection loss of-50.58 dB with the low matching thicknesses.Furthermore,the obtained 2D/3D RGO/CFs vdWs heterostructures also displayed the significant radar stealth properties,good corrosion resistance performances as well as outstanding thermal insulation capabilities,displaying the great potential in complex and variable environments.Accordingly,this work not only demonstrated a straightforward method for fabricating 2D/3D vdWs heterostructures,but also outlined a powerful mixeddimensional assembly strategy for engineering multifunctional foams for electromagnetic protection,aerospace and other complex conditions.展开更多
Combining 3D printing with precursor-derived ceramic for fabricating electromagnetic(EM) wave-absorbing metamaterials has attracted great attention. This study presents a novel ultraviolet-curable polysiloxane precurs...Combining 3D printing with precursor-derived ceramic for fabricating electromagnetic(EM) wave-absorbing metamaterials has attracted great attention. This study presents a novel ultraviolet-curable polysiloxane precursor for digital light processing(DLP) 3D printing to fabricate ceramic parts with complex geometry, no cracks and linear shrinkage. Guiding with the principles of impedance matching, attenuation, and effective-medium theory, we design a crosshelix-array metamaterial model based on the complex permittivity constant of precursor-derived ceramics. The corresponding ceramic metamaterials can be successfully prepared by DLP printing and subsequent pyrolysis process, achieving a low reflection coefficient and a wide effective absorption bandwidth in the X-band even under high temperature. This is a general method that can be extended to other bands, which can be realized by merely adjusting the unit structure of meta-materials. This strategy provides a novel and effective avenue to achieve “target-design-fabricating” ceramic metamaterials, and it exposes the downstream applications of highly efficient and broad EM wave-absorbing materials and structures with great potential applications.展开更多
Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coeff...Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.展开更多
A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves...A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves produced by an explosive source. These numerical results indicate that the code MMIC3D has the ability in computing cases such as 3D shock waves produced by air explosion, vortex region of the shock wave, the Mach wave, and reflected waves behind rigid boundaries.展开更多
In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numer...In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases: The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver(used for solving the discretized governing equations in the computational domain) upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.展开更多
Computations for air gap response of a semisubmersible platform based on a 3D numerical wave tank approach are presented.The developed method is in time domain and can consider nonlinearities associated with incident ...Computations for air gap response of a semisubmersible platform based on a 3D numerical wave tank approach are presented.The developed method is in time domain and can consider nonlinearities associated with incident wave and hydrostatic forces exactly in determining the body response, but the interaction hydrodynamics of radiation and diffraction are based on simplified linearization assumptions. The incident wave can be defined by any suitable wave theory and here defined by a fully nonlinear numerical wave model. After verifying the present computations results in its degenerated linearized version against the usual linear 3D Green function–based frequency-domain results for air gap predictions, systematic comparative studies are undertaken between linear and the approximate nonlinear solutions. It is found that nonlinear computations can yield considerably conservative predictions as compared to fully linear calculations, amounting to a difference of up to 30%–40% in the minimum air gap in steep ambient incident waves at high and moderate frequencies.展开更多
In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform ...In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order (M=1) recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.展开更多
A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity m...A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity model and earthquake relocation for this region are obtained using Pg and Sg phase readings of 9 988 earthquakes from 1992 to 1999 recorded at 193 seismic stations within central-western China by SPHYPIT90 and SPHREL3D90 programs. A lateral inhomogeneous structure of P-wave velocity in this region was obtained. Ob-vious contrast of P-wave velocities was revealed on both sides of active fault zones. Relocated epicenters of 6 459 events show clear lineation along active faults, which indicated a close correlation between seismicity and the active faults in this region. Focal depths of 82% relocated events ranged from 0 to 20 km, which is in good agreement with that from double-difference earthquake location algorithm.展开更多
The phenomenon of electrical attraction and repulsion between charged particles is well known, and described mathematically by Coulomb’s Law, yet until now there has been no explanation for why this occurs. There has...The phenomenon of electrical attraction and repulsion between charged particles is well known, and described mathematically by Coulomb’s Law, yet until now there has been no explanation for why this occurs. There has been no mechanistic explanation that reveals what causes the charged particles to accelerate, either towards or away from each other. This paper gives a detailed explanation of the phenomena of electrical attraction and repulsion based on my previous work that determined the exact wave-function solutions for both the Electron and the Positron. It is revealed that the effects are caused by wave interactions between the wave functions that result in Electromagnetic reflections of parts of the particle’s wave functions, causing a change in their momenta.展开更多
In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain...In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.展开更多
Textile-reinforced composites,due to their excellent highstrength-to-low-mass ratio, provide promising alternatives to conventional structural materials in many high-tech sectors. 3D braided composites are a kind of a...Textile-reinforced composites,due to their excellent highstrength-to-low-mass ratio, provide promising alternatives to conventional structural materials in many high-tech sectors. 3D braided composites are a kind of advanced composites reinforced with 3D braided fabrics; the complex nature of 3D braided composites makes the evaluation of the quality of the product very difficult. In this investigation,a defect recognition platform for 3D braided composites evaluation was constructed based on dual-tree complex wavelet packet transform( DT-CWPT) and backpropagation( BP) neural networks. The defects in 3D braided composite materials were probed and detected by an ultrasonic sensing system. DT-CWPT method was used to analyze the ultrasonic scanning pulse signals,and the feature vectors of these signals were extracted into the BP neural networks as samples. The type of defects was identified and recognized with the characteristic ultrasonic wave spectra. The position of defects for the test samples can be determined at the same time. This method would have great potential to evaluate the quality of 3D braided composites.展开更多
In this paper, we firstly derive the stability conditions of high-order staggered-grid schemes for the three-dimensional (3D) elastic wave equation in heterogeneous media based on the energy method. Moreover, the plan...In this paper, we firstly derive the stability conditions of high-order staggered-grid schemes for the three-dimensional (3D) elastic wave equation in heterogeneous media based on the energy method. Moreover, the plane wave analysis yields a sufficient and necessary stability condition by the von Neumann criterion in homogeneous case. Numerical computations for 3D wave simulation with point source excitation are given.展开更多
We present a theoretical investigation of the scattering of high frequency S0 Lamb mode from a circular blind hole defect in a plate based on the 3D theory. The SO wave is incident at the frequency above the A1 mode c...We present a theoretical investigation of the scattering of high frequency S0 Lamb mode from a circular blind hole defect in a plate based on the 3D theory. The SO wave is incident at the frequency above the A1 mode cut-off frequency, in which the popular approximate plate theories are inapplicable. Due to the non-symmetric blind hole defect, the scattered fields will contain higher order converted modes in addition to the fundamental SO and AO modes. The far-field scattering amplitudes of various propagating Lamb modes for different hole sizes are inspected. The results are compared with those of lower frequencies and some different phenomena are found. Two-dimensional Fourier transform (2DFT) results of transient scattered Lamb and SH wave signals agree well with the analytical dispersion curves, which check the validity of the solutions from another point of view.展开更多
Generally, FD coefficients can be obtained by using Taylor series expansion (TE) or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a...Generally, FD coefficients can be obtained by using Taylor series expansion (TE) or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares (LS) can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional (2D) forward modeling to three-dimensional (3D) and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.展开更多
Using over 3 500 first P arrival times recorded by nine digital seismic stations from Hainan Digital Seismic Net-work during 1999~2005,a 3-D P-wave velocity model of the crust under Hainan Island and adjacent regions...Using over 3 500 first P arrival times recorded by nine digital seismic stations from Hainan Digital Seismic Net-work during 1999~2005,a 3-D P-wave velocity model of the crust under Hainan Island and adjacent regions has been determined. The results show that the pattern of velocity anomalies in the shallower upper crust is somewhat associated with the surface geological tectonics in the region. A relative low-velocity anomaly appears north of the Wangwu-Wenjiao fault zone and a relative high-velocity anomaly appears south of the Wangwu-Wenjiao fault zone,corresponding to the depressed areas in north Hainan Island,where many volcanoes are frequently active and geothermal values are relatively higher,and the uplifted and stable regions in central and south of the Hainan Is-land. In the middle and lower crust velocities are relatively lower in east Hainan than those in west Hainan,possi-bly suggesting the existence of the upwelling of hot materials from the mantle in east Hainan. The pattern of veloc-ity anomalies also indicates that NW faults,i.e.,the Puqian-Qinglan fault,may be shallower,while the E-W Wangwu-Wenjiao fault may be deeper,which perhaps extends down to Moho depth or deeper.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant No.51239007)
文摘The size and shape of the effective test area are crucial to consider when short-crested waves are created by segmented wavemakers. The range of the effective test area of short-crested waves simulated by two-sided segmented wavemakers is analyzed in this paper. The experimental investigation on the wave field distribution of short-crested waves generated by two-sided segmented wavemakers is conducted by using an array of wave gauges. Wave spectra and directional spreading function are analyzed and the results show that when the main direction is at a certain angle with the normal line of wave generators, the wave field of 3D short-crested waves generated by two-sided segmented wavemakers has good spatial uniformity within the model test area. The effective test area can provide good wave environments for seakeeping model tests of various ocean engineering structures in the deep ocean engineering basin.
文摘The wave/particle duality of particles in Physics is well known. Particles have properties that uniquely characterize them from one another, such as mass, charge and spin. Charged particles have associated Electric and Magnetic fields. Also, every moving particle has a De Broglie wavelength determined by its mass and velocity. This paper shows that all of these properties of a particle can be derived from a single wave function equation for that particle. Wave functions for the Electron and the Positron are presented and principles are provided that can be used to calculate the wave functions of all the fundamental particles in Physics. Fundamental particles such as electrons and positrons are considered to be point particles in the Standard Model of Physics and are not considered to have a structure. This paper demonstrates that they do indeed have structure and that this structure extends into the space around the particle’s center (in fact, they have infinite extent), but with rapidly diminishing energy density with the distance from that center. The particles are formed from Electromagnetic standing waves, which are stable solutions to the Schrödinger and Classical wave equations. This stable structure therefore accounts for both the wave and particle nature of these particles. In fact, all of their properties such as mass, spin and electric charge, can be accounted for from this structure. These particle properties appear to originate from a single point at the center of the wave function structure, in the same sort of way that the Shell theorem of gravity causes the gravity of a body to appear to all originate from a central point. This paper represents the first two fully characterized fundamental particles, with a complete description of their structure and properties, built up from the underlying Electromagnetic waves that comprise these and all fundamental particles.
基金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.
基金Key Project Process Mechanism and Prediction of Geological Hazards (2001CB711005-1-3) and State Key Basic Research Project Mechanism and Prediction of Continental Earthquakes (G1998040702). sponsored by the Ministry of Science and Techno
基金provided by Guizhou Provincial Science and Technology Projects for Platform and Talent Team Plan(GCC[2023]007)Fok Ying Tung Education Foundation(171095)National Natural Science Foundation of China(11964006).
文摘Considering the serious electromagnetic wave(EMW)pollution problems and complex application condition,there is a pressing need to amalgamate multiple functionalities within a single substance.However,the effective integration of diverse functions into designed EMW absorption materials still faces the huge challenges.Herein,reduced graphene oxide/carbon foams(RGO/CFs)with two-dimensional/three-dimensional(2D/3D)van der Waals(vdWs)heterostructures were meticulously engineered and synthesized utilizing an efficient methodology involving freeze-drying,immersing absorption,secondary freeze-drying,followed by carbonization treatment.Thanks to their excellent linkage effect of amplified dielectric loss and optimized impedance matching,the designed 2D/3D RGO/CFs vdWs heterostructures demonstrated commendable EMW absorption performances,achieving a broad absorption bandwidth of 6.2 GHz and a reflection loss of-50.58 dB with the low matching thicknesses.Furthermore,the obtained 2D/3D RGO/CFs vdWs heterostructures also displayed the significant radar stealth properties,good corrosion resistance performances as well as outstanding thermal insulation capabilities,displaying the great potential in complex and variable environments.Accordingly,this work not only demonstrated a straightforward method for fabricating 2D/3D vdWs heterostructures,but also outlined a powerful mixeddimensional assembly strategy for engineering multifunctional foams for electromagnetic protection,aerospace and other complex conditions.
基金supported by the National Science Fund for Distinguished Young Scholars(52025034)National Natural Science Foundation of China(21975204)Innovation Team of Shaanxi Sanqin Scholars。
文摘Combining 3D printing with precursor-derived ceramic for fabricating electromagnetic(EM) wave-absorbing metamaterials has attracted great attention. This study presents a novel ultraviolet-curable polysiloxane precursor for digital light processing(DLP) 3D printing to fabricate ceramic parts with complex geometry, no cracks and linear shrinkage. Guiding with the principles of impedance matching, attenuation, and effective-medium theory, we design a crosshelix-array metamaterial model based on the complex permittivity constant of precursor-derived ceramics. The corresponding ceramic metamaterials can be successfully prepared by DLP printing and subsequent pyrolysis process, achieving a low reflection coefficient and a wide effective absorption bandwidth in the X-band even under high temperature. This is a general method that can be extended to other bands, which can be realized by merely adjusting the unit structure of meta-materials. This strategy provides a novel and effective avenue to achieve “target-design-fabricating” ceramic metamaterials, and it exposes the downstream applications of highly efficient and broad EM wave-absorbing materials and structures with great potential applications.
文摘Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.
文摘A numerical method is presented that simulates 3D explosive field problems. A code MMIC3D using this method can be used to simulate the propagation and reflected effects of all kinds of rigid boundaries to shock waves produced by an explosive source. These numerical results indicate that the code MMIC3D has the ability in computing cases such as 3D shock waves produced by air explosion, vortex region of the shock wave, the Mach wave, and reflected waves behind rigid boundaries.
文摘In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases: The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver(used for solving the discretized governing equations in the computational domain) upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.
文摘Computations for air gap response of a semisubmersible platform based on a 3D numerical wave tank approach are presented.The developed method is in time domain and can consider nonlinearities associated with incident wave and hydrostatic forces exactly in determining the body response, but the interaction hydrodynamics of radiation and diffraction are based on simplified linearization assumptions. The incident wave can be defined by any suitable wave theory and here defined by a fully nonlinear numerical wave model. After verifying the present computations results in its degenerated linearized version against the usual linear 3D Green function–based frequency-domain results for air gap predictions, systematic comparative studies are undertaken between linear and the approximate nonlinear solutions. It is found that nonlinear computations can yield considerably conservative predictions as compared to fully linear calculations, amounting to a difference of up to 30%–40% in the minimum air gap in steep ambient incident waves at high and moderate frequencies.
基金China Postdoctoral Science Foundation Under Grant No.20100480321National Basic Research Program of China Under Grant No. 2007CB714200
文摘In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order (M=1) recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.
文摘A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity model and earthquake relocation for this region are obtained using Pg and Sg phase readings of 9 988 earthquakes from 1992 to 1999 recorded at 193 seismic stations within central-western China by SPHYPIT90 and SPHREL3D90 programs. A lateral inhomogeneous structure of P-wave velocity in this region was obtained. Ob-vious contrast of P-wave velocities was revealed on both sides of active fault zones. Relocated epicenters of 6 459 events show clear lineation along active faults, which indicated a close correlation between seismicity and the active faults in this region. Focal depths of 82% relocated events ranged from 0 to 20 km, which is in good agreement with that from double-difference earthquake location algorithm.
文摘The phenomenon of electrical attraction and repulsion between charged particles is well known, and described mathematically by Coulomb’s Law, yet until now there has been no explanation for why this occurs. There has been no mechanistic explanation that reveals what causes the charged particles to accelerate, either towards or away from each other. This paper gives a detailed explanation of the phenomena of electrical attraction and repulsion based on my previous work that determined the exact wave-function solutions for both the Electron and the Positron. It is revealed that the effects are caused by wave interactions between the wave functions that result in Electromagnetic reflections of parts of the particle’s wave functions, causing a change in their momenta.
基金supported by the National Natural Science Foundation of China(Project U1901602&41790465)Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory(Guangzhou)(GML2019ZD0203)+2 种基金Shenzhen Key Laboratory of Deep Offshore Oil and Gas Exploration Technology(Grant No.ZDSYS20190902093007855)Shenzhen Science and Technology Program(Grant No.KQTD20170810111725321)the leading talents of Guangdong province program(Grant No.2016LJ06N652).
文摘In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.
基金National Natural Science Foundation of China(No.51303131)
文摘Textile-reinforced composites,due to their excellent highstrength-to-low-mass ratio, provide promising alternatives to conventional structural materials in many high-tech sectors. 3D braided composites are a kind of advanced composites reinforced with 3D braided fabrics; the complex nature of 3D braided composites makes the evaluation of the quality of the product very difficult. In this investigation,a defect recognition platform for 3D braided composites evaluation was constructed based on dual-tree complex wavelet packet transform( DT-CWPT) and backpropagation( BP) neural networks. The defects in 3D braided composite materials were probed and detected by an ultrasonic sensing system. DT-CWPT method was used to analyze the ultrasonic scanning pulse signals,and the feature vectors of these signals were extracted into the BP neural networks as samples. The type of defects was identified and recognized with the characteristic ultrasonic wave spectra. The position of defects for the test samples can be determined at the same time. This method would have great potential to evaluate the quality of 3D braided composites.
文摘In this paper, we firstly derive the stability conditions of high-order staggered-grid schemes for the three-dimensional (3D) elastic wave equation in heterogeneous media based on the energy method. Moreover, the plane wave analysis yields a sufficient and necessary stability condition by the von Neumann criterion in homogeneous case. Numerical computations for 3D wave simulation with point source excitation are given.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11474195,11274226 and 61171145
文摘We present a theoretical investigation of the scattering of high frequency S0 Lamb mode from a circular blind hole defect in a plate based on the 3D theory. The SO wave is incident at the frequency above the A1 mode cut-off frequency, in which the popular approximate plate theories are inapplicable. Due to the non-symmetric blind hole defect, the scattered fields will contain higher order converted modes in addition to the fundamental SO and AO modes. The far-field scattering amplitudes of various propagating Lamb modes for different hole sizes are inspected. The results are compared with those of lower frequencies and some different phenomena are found. Two-dimensional Fourier transform (2DFT) results of transient scattered Lamb and SH wave signals agree well with the analytical dispersion curves, which check the validity of the solutions from another point of view.
基金supported by the National Natural Science Foundation of China(No.41474110)Shell Ph.D. Scholarship to support excellence in geophysical research
文摘Generally, FD coefficients can be obtained by using Taylor series expansion (TE) or optimization methods to minimize the dispersion error. However, the TE-based FD method only achieves high modeling precision over a limited range of wavenumbers, and produces large numerical dispersion beyond this range. The optimal FD scheme based on least squares (LS) can guarantee high precision over a larger range of wavenumbers and obtain the best optimization solution at small computational cost. We extend the LS-based optimal FD scheme from two-dimensional (2D) forward modeling to three-dimensional (3D) and develop a 3D acoustic optimal FD method with high efficiency, wide range of high accuracy and adaptability to parallel computing. Dispersion analysis and forward modeling demonstrate that the developed FD method suppresses numerical dispersion. Finally, we use the developed FD method to source wavefield extrapolation and receiver wavefield extrapolation in 3D RTM. To decrease the computation time and storage requirements, the 3D RTM is implemented by combining the efficient boundary storage with checkpointing strategies on GPU. 3D RTM imaging results suggest that the 3D optimal FD method has higher precision than conventional methods.
基金The special project of Detection of Haikou City Earthquake Active Faults from the Tenth Five-year Plan of China Earthquake Administration (0106512)Joint Seismological Science Foundation of China (105086)CAS Key Laboratory of Marginal Sea Geology (MSGL0503).
文摘Using over 3 500 first P arrival times recorded by nine digital seismic stations from Hainan Digital Seismic Net-work during 1999~2005,a 3-D P-wave velocity model of the crust under Hainan Island and adjacent regions has been determined. The results show that the pattern of velocity anomalies in the shallower upper crust is somewhat associated with the surface geological tectonics in the region. A relative low-velocity anomaly appears north of the Wangwu-Wenjiao fault zone and a relative high-velocity anomaly appears south of the Wangwu-Wenjiao fault zone,corresponding to the depressed areas in north Hainan Island,where many volcanoes are frequently active and geothermal values are relatively higher,and the uplifted and stable regions in central and south of the Hainan Is-land. In the middle and lower crust velocities are relatively lower in east Hainan than those in west Hainan,possi-bly suggesting the existence of the upwelling of hot materials from the mantle in east Hainan. The pattern of veloc-ity anomalies also indicates that NW faults,i.e.,the Puqian-Qinglan fault,may be shallower,while the E-W Wangwu-Wenjiao fault may be deeper,which perhaps extends down to Moho depth or deeper.