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.展开更多
In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibi...In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.展开更多
3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces o...The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.展开更多
Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of t...Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference (FD) method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.展开更多
Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time st...Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference (symplectic FFD) method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.展开更多
Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulatio...Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulation of the acoustic wave equation to marine seismic modeling using the staggered-grid finite-difference method. The scheme is developed using a fourth-order spatial and a second-order temporal operator. Then, we define a stability coefficient (SC) and calculate its maximum value under the stability condition. Based on the dispersion relationship, we conduct a detailed dispersion analysis for submarine sediments in terms of the phase and group velocity over a range of angles, stability coefficients, and orders. We also compare the numerical solution with the exact solution for a P-wave line source in a homogeneous submarine model. Additionally, the numerical results determined by a Marmousi2 model with a rugged seafloor indicate that this method is sufficient for modeling complex submarine structures.展开更多
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.展开更多
Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled me...Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.展开更多
A steady state finite difference method is used to calculate EM fields generated by an MW broadcasting antenna. The effects of buildings on the wave propagation is studied based on the numerical results and field meas...A steady state finite difference method is used to calculate EM fields generated by an MW broadcasting antenna. The effects of buildings on the wave propagation is studied based on the numerical results and field measurements. Both the algorithm and results are useful in the design of broadcasting antennas, as well as in the selection of transmitting sites.展开更多
The elastic Helmholtz equation is capable of readily simulating attenuation and dispersion behaviors of the elastic wave and performing full-wavefield modeling in wave-equation-based elastic inversions and migrations....The elastic Helmholtz equation is capable of readily simulating attenuation and dispersion behaviors of the elastic wave and performing full-wavefield modeling in wave-equation-based elastic inversions and migrations.However,solving the elastic Helmholtz equation using a finite-difference frequency-domain(FDFD)method is computationally prohibitive especially in heterogeneous media with fine-scale heterogeneities.The FDFD method usually leads to a large discrete linear system of the elastic Helmholtz equation.We develop a multiscale method of FDFD to solve the elastic Helmholtz equation in isotropic media based on the general framework of heterogeneous multiscale method(HMM).The HMM framework decomposes the elastic Helmholtz problem into a series of microscale problems and a macroscale problem.The idea of multiscale basis functions is introduced to decouple the coupled microscale and macroscale problems and to capture fine-scale heterogeneity in medium properties.A reconstruction-based downscaling coupling and a flux-based upscaling coupling are used to convey the fine-scale medium heterogeneity to a coarse scale.The dimension of the resulting linear system is much smaller than those of linear systems generated with the conventional FDFD methods.We use a homogeneous model and a heterogeneous model to investigate the effects of the size of local sampling domains and the coarse-element number per S-wave wavelength on the accuracy of our new method,and employ two highly heterogeneous models to demonstrate the superiority in terms of the efficiency and memory consumption of our method based on the optimal local sampling-domain size and stable coarse-mesh discretization.The results demonstrate that our new method can approximate the finescale reference FDFD solutions with a significant decrease in computational complexity.展开更多
The relativistic Hartree-Bogoliubov(RHB)theory is a powerful tool for describing exotic nuclei near drip lines.The key technique is to solve the RHB equation in the coordinate space to obtain the quasi-particle states...The relativistic Hartree-Bogoliubov(RHB)theory is a powerful tool for describing exotic nuclei near drip lines.The key technique is to solve the RHB equation in the coordinate space to obtain the quasi-particle states.In this paper,we solve the RHB equation with the Woods-Saxon-type mean-field and Delta-type pairing-field potentials by using the finite-difference method(FDM).We inevitably obtain spurious states when using the common symmetric central difference formula(CDF)to construct the Hamiltonian matrix,which is similar to the problem resulting from solving the Dirac equation with the same method.This problem is solved by using the asymmetric difference formula(ADF).In addition,we show that a large enough box is necessary to describe the continuum quasi-particle states.The canonical states obtained by diagonalizing the density matrix constructed by the quasi-particle states are not particularly sensitive to the box size.Part of the asymptotic wave functions can be improved by applying the ADF in the FDM compared to the shooting method with the same box boundary condition.展开更多
With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving larg...With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.展开更多
This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the...This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo's amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^(3+):YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp(α'L),which is previously employed to describe the relationship between echo's efficiency and thickness,should be modified as 1.3 · 0.09 exp(2.4 ·α'L) with the propagation of echo considered.展开更多
The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the ...The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell's equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.展开更多
This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SE...This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.展开更多
Resistant starch(RS)comprises starch fractions that resist digestion in the small intestine and reach the colon,where they are fermented by the microbiota.Resistant starch harbors functional properties and healthpromo...Resistant starch(RS)comprises starch fractions that resist digestion in the small intestine and reach the colon,where they are fermented by the microbiota.Resistant starch harbors functional properties and healthpromoting ingredients that can regulate blood glucose and lipid levels,prevent cancer,and enhance the quality of life.Consequently,new technologies for the preparation of RS are continually being developed to support its industrial production.This review describes the structural and nutritional properties of RS and examines recent advancements in RS preparation methods.Emphasis is placed on how RS structure influences its properties and the physiological mechanisms in vivo.This review aims to stimulate further research into the preparation methods,functional characteristics,and utilization of RS,thereby supporting ongoing developments in the food industry.展开更多
Traditional targeted analyses often overlook unknown or emerging contaminants,highlighting the significance of nontarget and suspect screening approaches.A novel and high-sensitivity methodology for nontarget analysis...Traditional targeted analyses often overlook unknown or emerging contaminants,highlighting the significance of nontarget and suspect screening approaches.A novel and high-sensitivity methodology for nontarget analysis of organic pollutants in human serum was newly-developed based on gas chromatography coupled with quadrupole time-of-flight high-resolution mass spectrometry.The extraction protocol employing an acetonitrile-ethyl acetate(9:1,V:V)mixture significantly improved the extraction efficiency while minimizing matrix effect.A hybridized analytical strategy integrating nontarget and suspect screening was developed to achieve comprehensive identification and classification of pollutants,employing the National Institute of Standards and Technology(NIST)20 library and Agilent Technologies Personal Compound Database and Library(PCDL).This approach successfully characterized 273 organic contaminants spanning 12 categories,including polycyclic aromatic hydrocarbons(PAHs)and their derivatives,esters,and phenolic compounds in human serum,with a significant increase in detection specificity compared to conventional workflows.The methodology used serum samples of the workers from coking industry,revealing widespread contamination dominated by PAHs and PAH derivatives.Among the target analytes,three were identified solely by NIST and six solely by PCDL,indicating the complementary benefits of combining these different databases.Notably,this work reported the first confirmed detection of 2-naphthalenamine in human serum.This optimized approach demonstrates enhanced sensitivity and reliability in serum analysis,advancing biomonitoring capabilities and providing a deep understanding of human exposure to environmental pollutants.展开更多
Advancements in tumor immunotherapy highlight the significant potential of antibody drugs,a key category of biological agents,for treating cancer and autoimmune diseases.This paper begins by defining and classifying k...Advancements in tumor immunotherapy highlight the significant potential of antibody drugs,a key category of biological agents,for treating cancer and autoimmune diseases.This paper begins by defining and classifying key targets in tumor immunity,as well as discussing their structural and functional characteristics.Subsequently,it elaborates on innovative technologies for antibody drug screening,which,when integrated with contemporary molecular biology,biotechnology,and computational biology,have substantially enhanced the efficiency and accuracy of target identification and antibody drug screening processes.Despite the promising prospects of tumor immunotherapy,certain limitations persist in its practical implementation.In conclusion,this paper offers a comprehensive examination of the cutting-edge developments in tumor immunotherapy,focusing on the aspects of tumor immunotherapy itself,critical targets for immunotherapy,and novel technologies and methodologies for antibody screening.This analysis is crucial for advancing the field of tumor immunotherapy and for enhancing both therapeutic efficacy and safety.Furthermore,research and development(R&D)of antibody drugs in other domains,such as autoimmune and inflammatory diseases,can benefit from it.展开更多
基金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.
基金the National Natural Science Foundation of China(No.40774056)Program of Excellent Team in Harbin Institute of Technology
摘要In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
摘要3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
摘要In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations: the first-order (velocity-stress) and the second- order (displacement-stress) separate elastic wave equa- tions, with the first-order (velocity-stress) and the second- order (displacement-stress) full (or mixed) elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
基金Taif University Researchers Supporting Project Number(TURSP-2020/164),Taif University,Taif,Saudi Arabia.
摘要The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.
基金This research was supported by the National Nature Science Foundation of China (No. 41074100) and the Program for NewCentury Excellent Talents in the University of the Ministry of Education of China (No. NCET- 10-0812).
摘要Prestack reverse time migration (RTM) is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference (FD) method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.
基金supported by National Natural Science Foundation of China(41504109,41404099)the Natural Science Foundation of Shandong Province(BS2015HZ008)the project of "Distinguished Professor of Jiangsu Province"
摘要Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference (symplectic FFD) method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.
基金Supported by the National Natural Science Foundation of China(Nos. 41206043, 40930845)the Open Foundation of Key Laboratory of Marine Geology and Environment of Chinese Academy of Sciences(No. MGE2011KG07)+1 种基金the Knowledge Innovation Program of Chinese Academy of Sciences (No. KZCX2-YW-229)the National Basic Research Program of China (973 Program) (No. 2009CB219505)
摘要Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulation of the acoustic wave equation to marine seismic modeling using the staggered-grid finite-difference method. The scheme is developed using a fourth-order spatial and a second-order temporal operator. Then, we define a stability coefficient (SC) and calculate its maximum value under the stability condition. Based on the dispersion relationship, we conduct a detailed dispersion analysis for submarine sediments in terms of the phase and group velocity over a range of angles, stability coefficients, and orders. We also compare the numerical solution with the exact solution for a P-wave line source in a homogeneous submarine model. Additionally, the numerical results determined by a Marmousi2 model with a rugged seafloor indicate that this method is sufficient for modeling complex submarine structures.
摘要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.
基金Supported by Major Project of National Science and Technology of China(No.2016ZX05026-002-003)National Natural Science Foundation of China(No.41374108)
摘要Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.
摘要A steady state finite difference method is used to calculate EM fields generated by an MW broadcasting antenna. The effects of buildings on the wave propagation is studied based on the numerical results and field measurements. Both the algorithm and results are useful in the design of broadcasting antennas, as well as in the selection of transmitting sites.
基金the financial support received from the National Natural Foundation of China(Grant No.42374163)the Key Program of Natural Science Foundation of Sichuan Province(Grant No.2023NSFSC0019)the Open Fund of Key Laboratory of Earth Exploration and Information Techniques(Chengdu University of Technology),Ministry of Education(Grant No.EEME202505)。
摘要The elastic Helmholtz equation is capable of readily simulating attenuation and dispersion behaviors of the elastic wave and performing full-wavefield modeling in wave-equation-based elastic inversions and migrations.However,solving the elastic Helmholtz equation using a finite-difference frequency-domain(FDFD)method is computationally prohibitive especially in heterogeneous media with fine-scale heterogeneities.The FDFD method usually leads to a large discrete linear system of the elastic Helmholtz equation.We develop a multiscale method of FDFD to solve the elastic Helmholtz equation in isotropic media based on the general framework of heterogeneous multiscale method(HMM).The HMM framework decomposes the elastic Helmholtz problem into a series of microscale problems and a macroscale problem.The idea of multiscale basis functions is introduced to decouple the coupled microscale and macroscale problems and to capture fine-scale heterogeneity in medium properties.A reconstruction-based downscaling coupling and a flux-based upscaling coupling are used to convey the fine-scale medium heterogeneity to a coarse scale.The dimension of the resulting linear system is much smaller than those of linear systems generated with the conventional FDFD methods.We use a homogeneous model and a heterogeneous model to investigate the effects of the size of local sampling domains and the coarse-element number per S-wave wavelength on the accuracy of our new method,and employ two highly heterogeneous models to demonstrate the superiority in terms of the efficiency and memory consumption of our method based on the optimal local sampling-domain size and stable coarse-mesh discretization.The results demonstrate that our new method can approximate the finescale reference FDFD solutions with a significant decrease in computational complexity.
基金Supported by the National Natural Science Foundation of China(11775119,2175109)the Natural Science Foundation of Tianjin,China(19JCYBJC30800)。
摘要The relativistic Hartree-Bogoliubov(RHB)theory is a powerful tool for describing exotic nuclei near drip lines.The key technique is to solve the RHB equation in the coordinate space to obtain the quasi-particle states.In this paper,we solve the RHB equation with the Woods-Saxon-type mean-field and Delta-type pairing-field potentials by using the finite-difference method(FDM).We inevitably obtain spurious states when using the common symmetric central difference formula(CDF)to construct the Hamiltonian matrix,which is similar to the problem resulting from solving the Dirac equation with the same method.This problem is solved by using the asymmetric difference formula(ADF).In addition,we show that a large enough box is necessary to describe the continuum quasi-particle states.The canonical states obtained by diagonalizing the density matrix constructed by the quasi-particle states are not particularly sensitive to the box size.Part of the asymptotic wave functions can be improved by applying the ADF in the FDM compared to the shooting method with the same box boundary condition.
基金Project(2025ZD1009704)supported by the National Science and Technology Major Project of ChinaProjects(2023JJ30659,2022JJ30706)supported by Hunan Provincial Natural Science Foundation,China。
摘要With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.
基金Project supported by Tianjin Research Program Application Foundation and Advanced Technology,China(Grant No.15JCQNJC01100)
摘要This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo's amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^(3+):YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp(α'L),which is previously employed to describe the relationship between echo's efficiency and thickness,should be modified as 1.3 · 0.09 exp(2.4 ·α'L) with the propagation of echo considered.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11304074,61475042,and 11274088)the Natural Science Foundation of Hebei Province,China(Grant Nos.A2015202320 and GCC2014048)the Key Subject Construction Project of Hebei Province University,China
摘要The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell's equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.
基金National Natural Science Foundation of China under Grant Nos.U2139208 and 52278516Key Laboratory of Earthquake Engineering and Engineering Vibration,China Earthquake Administration under Grant No.2024D15Key Laboratory of Soft Soil Characteristic and Engineering Environment,Tianjin Chengjian University under Grant No.2022SCEEKL003。
摘要This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.
基金financially supported by the National Key Research and Development Program of China(2023YFD2100803)the National Natural Science Foundation of China(32372387)+2 种基金the Science and Technology Major Project of Heilongjiang China(2021ZX12B07)Collaborative Innovation Achievement Project of“Double First-class”Disciplines in Heilongjiang Province(LJGXCG202080LJGXCG202083)。
摘要Resistant starch(RS)comprises starch fractions that resist digestion in the small intestine and reach the colon,where they are fermented by the microbiota.Resistant starch harbors functional properties and healthpromoting ingredients that can regulate blood glucose and lipid levels,prevent cancer,and enhance the quality of life.Consequently,new technologies for the preparation of RS are continually being developed to support its industrial production.This review describes the structural and nutritional properties of RS and examines recent advancements in RS preparation methods.Emphasis is placed on how RS structure influences its properties and the physiological mechanisms in vivo.This review aims to stimulate further research into the preparation methods,functional characteristics,and utilization of RS,thereby supporting ongoing developments in the food industry.
基金supported by the National Key Research and Development Project(Nos.2023YFC3905102 and 2024YFC3713201)the National Natural Science Foundation of China(Nos.42207485 and 42407567).
摘要Traditional targeted analyses often overlook unknown or emerging contaminants,highlighting the significance of nontarget and suspect screening approaches.A novel and high-sensitivity methodology for nontarget analysis of organic pollutants in human serum was newly-developed based on gas chromatography coupled with quadrupole time-of-flight high-resolution mass spectrometry.The extraction protocol employing an acetonitrile-ethyl acetate(9:1,V:V)mixture significantly improved the extraction efficiency while minimizing matrix effect.A hybridized analytical strategy integrating nontarget and suspect screening was developed to achieve comprehensive identification and classification of pollutants,employing the National Institute of Standards and Technology(NIST)20 library and Agilent Technologies Personal Compound Database and Library(PCDL).This approach successfully characterized 273 organic contaminants spanning 12 categories,including polycyclic aromatic hydrocarbons(PAHs)and their derivatives,esters,and phenolic compounds in human serum,with a significant increase in detection specificity compared to conventional workflows.The methodology used serum samples of the workers from coking industry,revealing widespread contamination dominated by PAHs and PAH derivatives.Among the target analytes,three were identified solely by NIST and six solely by PCDL,indicating the complementary benefits of combining these different databases.Notably,this work reported the first confirmed detection of 2-naphthalenamine in human serum.This optimized approach demonstrates enhanced sensitivity and reliability in serum analysis,advancing biomonitoring capabilities and providing a deep understanding of human exposure to environmental pollutants.
基金supported by the National Natural Science Foundation of China National(Grant Nos:32470999,31970882,81773261,81903140,82041012,82322055,82421005,82473278,92169115)the Shanghai Rising-Star Program(Grant No.:23QA1405800)+3 种基金The Shanghai Outstanding Academic Leader Program(Grant No.:23XD1424800)the Shanghai Key Laboratory of Cell Engineering(Grant No.:14DZ2272300)Yizhang Outstanding Academic Leader Program(Grant No.:JCYZRC-B-008)Cross-disciplinary research fund project of the Ninth People's Hospital affiliated to Shanghai Jiao Tong University School of Medicine(Grant No.:JCJC202410).
摘要Advancements in tumor immunotherapy highlight the significant potential of antibody drugs,a key category of biological agents,for treating cancer and autoimmune diseases.This paper begins by defining and classifying key targets in tumor immunity,as well as discussing their structural and functional characteristics.Subsequently,it elaborates on innovative technologies for antibody drug screening,which,when integrated with contemporary molecular biology,biotechnology,and computational biology,have substantially enhanced the efficiency and accuracy of target identification and antibody drug screening processes.Despite the promising prospects of tumor immunotherapy,certain limitations persist in its practical implementation.In conclusion,this paper offers a comprehensive examination of the cutting-edge developments in tumor immunotherapy,focusing on the aspects of tumor immunotherapy itself,critical targets for immunotherapy,and novel technologies and methodologies for antibody screening.This analysis is crucial for advancing the field of tumor immunotherapy and for enhancing both therapeutic efficacy and safety.Furthermore,research and development(R&D)of antibody drugs in other domains,such as autoimmune and inflammatory diseases,can benefit from it.