Abstract:Accurate three-dimensional(3D)velocity models are essential for fitting high-frequency seismic waveform records.This process usually requires regional-scale 3D numerical simulations that are computationally e...Abstract:Accurate three-dimensional(3D)velocity models are essential for fitting high-frequency seismic waveform records.This process usually requires regional-scale 3D numerical simulations that are computationally expensive,especially with sparse seismic networks.Because of the significance of source domain modeling,we propose a hybrid waveform simulation approach that combines the 3D spectral-element method(SEM)with the displacement representation theorem.By separating near-source wavefield excitation from long-distance wave propagation to stations,only the source domain wavefield needs to be recomputed when the local velocity and source models change.We apply the method to the 2019 M_(w)5.0 Changning shallow earthquake to verify its flexibility and effectiveness.We compare high-frequency waveforms computed with different regional velocity models against observations.Results show that the hybrid method achieves accuracy comparable to full SEM 3D simulations while reducing computation costs by more than two orders of magnitude when the structure of the source region updates.Our results further indicate that high-frequency waveforms are highly sensitive to shallow structures.Introducing low-velocity shallow layers into the source region improves near-field waveform fits,indicating pronounced low-velocity sediments in the Changning area.Large surface-wave time delays suggest that shallow velocities within the Sichuan Basin are lower than those in existing published models.In addition,an Interferometric Synthetic Aperture Radar(InSAR)-derived finite-fault model outperforms the point-source model in near-field waveform fitting and better reproduces rupture directivity.The proposed method is practical for high-frequency waveform modeling in areas with complex subsurface structures and rupture processes.展开更多
With the depletion of fossil fuels and increasing environmental concerns,the development of renewable energy,such as wave energy,has become a critical component of global energy strategies.However,challenges persist i...With the depletion of fossil fuels and increasing environmental concerns,the development of renewable energy,such as wave energy,has become a critical component of global energy strategies.However,challenges persist in the field testing methodologies for wave energy converters(WECs).In this paper,a numerical wave field of the Dawanshan Island Sea Area in Zhuhai City is constructed based on the MIKE21 SW wave model and by using an NCEP wind field driving model.In conjunction with the IEC-62600-100 standard,by taking site testing of the“Wanshan”wave energy converter on which a sea trial has been conducted in Dawanshan Island of Zhuhai city as an example,research on-site testing method for a wave energy converter has been carried out.The wave measurement position for the“Wanshan”converter was determined by combining statistically analyzed field data with a validated numerical wave model.By comparing a valid wave height at the position where a wave rider is located with a valid wave height at the position where the“Wanshan”wave energy converter is situated,the correlation coefficient between simulation and observed data reached 0.90,with a root-mean-square error of 0.19.The representativeness of wave measurement data during site testing is verified and can be used as a basis for calculating the input energy of the“Wanshan”wave energy converter.展开更多
In this paper,we investigate the periodic traveling wave solutions problem for a single population model with advection and distributed delay.By the bifurcation analysis method,we can obtain periodic traveling wave so...In this paper,we investigate the periodic traveling wave solutions problem for a single population model with advection and distributed delay.By the bifurcation analysis method,we can obtain periodic traveling wave solutions for this model under the influence of advection term and distributed delay.The obtained results indicate that weak kernel and strong kernel can both deduce the existence of periodic traveling wave solutions.Finally,we apply the main results in this paper to Logistic model and Nicholson’s blowflies model.展开更多
In multi-component oil and gas exploration using ocean bottom nodes,converted wave data is rich in lithological and fracture information.One of the urgent problems to be solved is how to construct an accurate shear wa...In multi-component oil and gas exploration using ocean bottom nodes,converted wave data is rich in lithological and fracture information.One of the urgent problems to be solved is how to construct an accurate shear wave velocity model of the shallow sea bottom by leveraging the seismic wave information at the fluid-solid interface in the ocean,and improve the lateral resolution of marine converted wave data.Given that the dispersion characteristics of surface waves are sensitive to the S-wave velocity of subsurface media,and that Scholte surface waves,which propagate at the interface between liquid and solid media,exist in the data of marine oil and gas exploration,this paper proposes a Scholte wave inversion and modeling method based on oil and gas exploration using ocean bottom nodes.By using the method for calculating the Scholte wave dispersion spectrum based on the Bessel kernel function,the accuracy of dispersion spectrum analysis is improved,and more accurate dispersion curves are picked up.Through the adaptive weighted least squares Scholte wave dispersion inversion algorithm,the Scholte wave dispersion equation for liquid-solid media is solved,and the shear wave velocity model of the shallow sea bottom is calculated.Theoretical tests and applications of realdata have proven that this method can significantly improve the lateral resolution of converted wave data,provide high-quality data for subsequent inversion of marine multi-component oil and gas exploration data and reservoir reflection information,and contribute to the development of marine oil and gas exploration technology.展开更多
This study applies Ensemble Optimal Interpolation(EnOI)to assimilate individual spectral components derived from National Data Buoy Center(NDBC)buoy directional spectra into the WAVEWATCHⅢ(WW3)wave model during tropi...This study applies Ensemble Optimal Interpolation(EnOI)to assimilate individual spectral components derived from National Data Buoy Center(NDBC)buoy directional spectra into the WAVEWATCHⅢ(WW3)wave model during tropical cyclone(TC)Isaias(2020).The analysis provides a comprehensive evaluation of the assimilation’s impact on wave parameters,frequency spectra,and directional spectra.Two series of assimilation experiments—one based on spectral components(Exp^(*)-DaSpec)and the other on significant wave height(Exp^(*)-DaSWH)—are evaluated against a non-assimilated control run(Exp-NoDa).Particular focus is placed on six key parameters:SWH,mean wave period(MWP),mean wave direction(MWD),mean wave directional spread(MWS),dominant wave period(DWP),and dominant wave direction(DWD).Sensitivity analyses suggest 400 km and 0.30 as appropriate values for the localization radius and observation error variance,respectively,though no single setting is optimal across all wave parameters.Exp4-DaSWH and Exp4-DaSpec are therefore selected as representative experiments.In non-independent validation,Exp4-DaSpec generally outperforms Exp-NoDa across MWP,MWD,MWS,DWP,and DWD,demonstrating closer agreement with observed frequency spectra and directional patterns.In independent validation,Exp4-DaSpec maintains superior overall performance,whereas Exp4-DaSWH shows only limited improvement.Exp4-DaSWH may capture the spectral peak near 0.1 Hz,although its directional characteristics remain largely similar to those of Exp-NoDa.Importantly,the assimilation experiments significantly improve SWH in both non-independent and independent validations,with Exp4-DaSWH slightly outperforming Exp4-DaSpec overall.展开更多
This paper investigates the impact of the model top and damping layer on the numerical simulation of tropical cyclones(TCs)and reveals the significant role of stratospheric gravity waves(SGWs).TCs can generate SGWs,wh...This paper investigates the impact of the model top and damping layer on the numerical simulation of tropical cyclones(TCs)and reveals the significant role of stratospheric gravity waves(SGWs).TCs can generate SGWs,which propagate upward and outward into the stratosphere.These SGWs can reach the damping layer,which is a consequence of the numerical scheme employed,where they can affect the tangential circulation through the dragging and forcing processes.In models with a higher top boundary,this tangential circulation develops far from the TC and has minimal direct impact on TC intensity.By comparison,in models with a lower top(e.g.,20 km),the damping layer is located just above the top of the TC.The SGW dragging in the damping layer and the consequent tangential force can thus induce ascent outside the eyewall,promote latent heat release,tilt the eyewall,and enlarge the inner-core radius.This process will reduce inner-core vorticity advection within the boundary layer,and eventually inhibits the intensification of the TC.This suggests that when the thickness of the damping layer is 5 km,the TC numerical model top height should be at least higher than 20 km to generate more accurate simulations.展开更多
Rock discontinuities such as joints widely exist in natural rock masses,and wave attenuation through rock masses is mainly caused by discontinuities.The displacement discontinuity model(DDM)has been widely used in the...Rock discontinuities such as joints widely exist in natural rock masses,and wave attenuation through rock masses is mainly caused by discontinuities.The displacement discontinuity model(DDM)has been widely used in theoretical and numerical analysis of wave propagation across rock discontinuity.However,the circumstance under which the DDM is applicable to predict wave propagation across rock discontinuity remains poorly understood.In this study,theoretical analysis and ultrasonic laboratory tests were carried out to examine the theoretical applicability of the DDM for wave propagation,where specimens with rough joints comprising regular rectangular asperities of different spacings and heights were prepared by 3D printing technology.It is found that the theoretical applicability of the DDM to predict wave propagation across rock discontinuity is determined by three joint parameters,i.e.the dimensionless asperity spacing(L),the dimensionless asperity height(H)and the groove density(D).Through theoretical analysis and laboratory tests,the conditions under which the DDM is applicable are derived as follows:and,.With increase in the groove density,the thresholds of the dimensionless asperity spacing and the dimensionless asperity height show a decreasing trend.In addition,the transmission coefficient in the frequency domain decreases with increasing groove density,dimensionless asperity spacing or dimensionless asperity height.The findings can facilitate our understanding of DDM for predicting wave propagation across rock discontinuity.展开更多
The internal solitary wave(ISW)represents a frequent and severe oceanic dynamic phenomenon observed in the South China Sea,exposing marine structures to sudden loads.This paper examines the prediction model of interac...The internal solitary wave(ISW)represents a frequent and severe oceanic dynamic phenomenon observed in the South China Sea,exposing marine structures to sudden loads.This paper examines the prediction model of interaction loads between ISW and FPSO,accounting for varying attack angles and incorporating ISW theories.The research demonstrates that the horizontal and transverse forces on FPSO under internal solitary waves(ISWs)comprise wave pressure difference force and viscous force,while the vertical force primarily consists of vertical wave pressure difference force.The wave pressure difference force is determined using the Froude-Krylov equation.The viscous force is derived from the tangential particle velocity induced by ISW and the viscous coefficient.The viscous coefficient formula is obtained through regression analysis of experimental data with different ISW attack angles.The research reveals that the horizontal viscous coefficient C_(vx)decreases as Reynolds number(R_(e))increases,while the transverse viscous coefficient C_(vy)initially increases and subsequently decreases with the growth of the Keulegan-Carpenter number(KC).Moreover,changes in wave propagation direction significantly affect the extreme magnitudes of both horizontal and transverse forces,and simultaneously modify the transverse force orientation,while having minimal impact on the vertical force.Additionally,the forces increase with the ISW’s amplitude.For horizontal and transverse forces,a thinner upper fluid layer generates larger forces.Comparative analysis of experimental,numerical,and theoretical results indicates strong agreement between theoretical predictions and experimental and numerical outcomes.展开更多
We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground...We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.展开更多
Studies of wave-current interactions are vital for the safe design of structures.Regular waves in the presence of uniform,linear shear,and quadratic shear currents are explored by the High-Level Green-Naghdi model in ...Studies of wave-current interactions are vital for the safe design of structures.Regular waves in the presence of uniform,linear shear,and quadratic shear currents are explored by the High-Level Green-Naghdi model in this paper.The five-point central difference method is used for spatial discretization,and the fourth-order Adams predictor-corrector scheme is employed for marching in time.The domain-decomposition method is applied for the wave-current generation and absorption.The effects of currents on the wave profile and velocity field are examined under two conditions:the same velocity of currents at the still-water level and the constant flow volume of currents.Wave profiles and velocity fields demonstrate substantial differences in three types of currents owing to the diverse vertical distribution of current velocity and vorticity.Then,loads on small-scale vertical cylinders subjected to regular waves and three types of background currents with the same flow volume are investigated.The maximum load intensity and load fluctuation amplitude in uniform,linear shear,and quadratic shear currents increase sequentially.The stretched superposition method overestimates the maximum load intensity and load fluctuation amplitude in opposing currents and underestimates these values in following currents.The stretched superposition method obtains a poor approximation for strong nonlinear waves,particularly in the case of the opposing quadratic shear current.展开更多
It is well known that aluminum and copper exhibit structural phase transformations in quasi-static and dynamic measurements,including shock wave loading.However,the dependence of phase transformations in a wide range ...It is well known that aluminum and copper exhibit structural phase transformations in quasi-static and dynamic measurements,including shock wave loading.However,the dependence of phase transformations in a wide range of crystallographic directions of shock loading has not been revealed.In this work,we calculated the shock Hugoniot for aluminum and copper in different crystallographic directions([100],[110],[111],[112],[102],[114],[123],[134],[221]and[401])of shock compression using molecular dynamics(MD)simulations.The results showed a high pressure(>160 GPa for Cu and>40 GPa for Al)of the FCC-to-BCC transition.In copper,different characteristics of the phase transition are observed depending on the loading direction with the[100]compression direction being the weakest.The FCC-to-BCC transition for copper is in the range of 150–220 GPa,which is consistent with the existing experimental data.Due to the high transition pressure,the BCC phase transition in copper competes with melting.In aluminum,the FCC-to-BCC transition is observed for all studied directions at pressures between 40 and 50 GPa far beyond the melting.In all considered cases we observe the coexistence of HCP and BCC phases during the FCC-to-BCC transition,which is consistent with the experimental data and atomistic calculations;this HCP phase forms in the course of accompanying plastic deformation with dislocation activity in the parent FCC phase.The plasticity incipience is also anisotropic in bothmetals,which is due to the difference in the projections of stress on the slip plane for different orientations of the FCC crystal.MD modeling results demonstrate a strong dependence of the FCC-to-BCC transition on the crystallographic direction,in which the material is loaded in the copper crystals.However,MD simulations data can only be obtained for specific points in the stereographic direction space;therefore,for more comprehensive understanding of the phase transition process,a feed-forward neural network was trained using MD modeling data.The trained machine learning model allowed us to construct continuous stereographic maps of phase transitions as a function of stress in the shock-compressed state of metal.Due to appearance and growth of multiple centers of new phase,the FCC-to-BCC transition leads to formation of a polycrystalline structure from the parent single crystal.展开更多
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.展开更多
Based on the Longuet-Higgins wave model theory, the previews studies have shown that freak waves can be generated in finite space and time successfully. However, as to generating high nonlinear freak waves, the simula...Based on the Longuet-Higgins wave model theory, the previews studies have shown that freak waves can be generated in finite space and time successfully. However, as to generating high nonlinear freak waves, the simulation results will be unrealistic. Therefore, a modified phase modulation method for simulating high nonlinear freak waves was developed. The surface elevations of some wave components at certain time and place are positive by modulating the corresponding random initial phases, then the total surface elevation at the focused point is enhanced and furthermore a freak wave event is generated. The new method can not only make the freak wave occur at certain time and place, but also make the simulated wave surface time series satisfy statistical properties of the realistic sea state and keep identical with the target wave spectrum. This numerical approach is of good precision and high efficiency by the comparisons of the simulated freak waves and the recorded freak waves.展开更多
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
Four focusing models for generation of freak waves are presented. An extreme wave focusing model is presented on the basis of the enhanced High-Order Spectral （HOS） method and the importance of the nonlinear wave-wa...Four focusing models for generation of freak waves are presented. An extreme wave focusing model is presented on the basis of the enhanced High-Order Spectral （HOS） method and the importance of the nonlinear wave-wave interaction is evaluated by comparison of the calculated results with experimental and theoretical data. Based on the modification of the Longuet-Higgins model, four wave models for generation of freak waves （a. extreme wave model ＋ random wave model; b. extreme wave model ＋ regular wave model; e. phase interval modulation wave focusing model; d. number modulation wave focusing model with the same phase） are proposed. By use of different energy distribution techniques in the four models, freak wave events are obtained with different Hmax/Hs in finite space and time.展开更多
High-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth with an arbitrary sloping bottom are presented in this article. First, the formal derivations to any high order of ...High-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth with an arbitrary sloping bottom are presented in this article. First, the formal derivations to any high order of mu(= h/lambda, depth to deep-water wave length ratio) and epsilon(= a/h, wave amplitude to depth ratio) for velocity potential, particle velocity vector, pressure and the Boussinesq-type equations for surface elevation eta and horizontal velocity vector (U) over right arrow at any given level in water are given. Then, the exact explicit expressions to the fourth order of mu are derived. Finally, the linear solutions of eta, (U) over right arrow, C (phase-celerity) and C-g (group velocity) for a constant water depth are obtained. Compared with the Airy theory, excellent results can be found even for a water depth as large as the wave legnth. The present high-order models are applicable to nonlinear regular and irregular waves in water of any varying depth (from shallow to deep) and bottom slope (from mild to steep).展开更多
The LAGFD-WAM wave model is a third generation wave model. In the present paper the physical aspect of the model was shown in great detail including energy spectrum balance equation, complicated characteristics equati...The LAGFD-WAM wave model is a third generation wave model. In the present paper the physical aspect of the model was shown in great detail including energy spectrum balance equation, complicated characteristics equations and source functions.展开更多
Typhoon-generated waves are simulated with two numerical wave models, the SWAN model for the coastal and Yangtze Estuary domain, nested within the WAVEWATCHIII （WW3） for the basin-scale East China Sea domain. Typhoo...Typhoon-generated waves are simulated with two numerical wave models, the SWAN model for the coastal and Yangtze Estuary domain, nested within the WAVEWATCHIII （WW3） for the basin-scale East China Sea domain. Typhoon No. 8114 is chosen because it was very strong, and generated high waves in the Estuary. WW3 was implemented for the East China Sea coarse-resolution computational domain, to simulate the waves over a large spatial scale and provide boundary conditions for SWAN model simulations, implemented on a fine-resolution nested domain for the Yangtze Estuary area. The Takahashi wind model is applied to the simulation of the East China Sea scale （3-hourly） and Yangtze Estuary scale （1-hourly） winds. Simulations of significant wave heights in the East China Sea show that the highest waves are on the right side of the storm track, and maxima tend to occur at the eastern deep-water open boundary of the Yangtze Estuary. In the Yangtze Estuary, incoming swell is dominant over locally generated waves before the typhoon approaches the Estuary. As the typhoon approaches the Estuary, wind waves and swell coexist, and the wave direction is mainly influenced by the swell direction and the complex topography.展开更多
In this study, the combined actions of waves and tidal currents in estuarine and coastal areas are considered and a 2D mathematical model for sediment transport by waves and tidal currents has been established in orth...In this study, the combined actions of waves and tidal currents in estuarine and coastal areas are considered and a 2D mathematical model for sediment transport by waves and tidal currents has been established in orthogonal curvilinear coordinates. Non-equilibrium transport equations of suspended load and bed load are used in the model. The concept of background concentration is introduced, and the formula of sediment transport capacity of tidal currents for the Oujiang River estuary is obtained. The Dou Guoren formula is employed for the sediment transport capacity of waves. Sediment transport capacity in the form of mud and the intensity of back silting are calculated by use of Luo Zaosen＇ s formula. The calculated tidal stages are in good agreement with the field data, and the calculated velocities and flow directions of 46 vertical lines for 8 cross sections are also in good agreement with the measured data. On such a basis, simulations of back silting after excavation of the waterway with a sand bar under complicated boundary conditions in the navigation channel induced by suspended load, bed load and mud by waves and tidal currents are discussed.展开更多
基金supported by the National Key Research and Development Program of China(Grant No.2020YFA0710600)the National Natural Science Foundation of China(Grant Nos.42325401 and U1939202).
文摘Abstract:Accurate three-dimensional(3D)velocity models are essential for fitting high-frequency seismic waveform records.This process usually requires regional-scale 3D numerical simulations that are computationally expensive,especially with sparse seismic networks.Because of the significance of source domain modeling,we propose a hybrid waveform simulation approach that combines the 3D spectral-element method(SEM)with the displacement representation theorem.By separating near-source wavefield excitation from long-distance wave propagation to stations,only the source domain wavefield needs to be recomputed when the local velocity and source models change.We apply the method to the 2019 M_(w)5.0 Changning shallow earthquake to verify its flexibility and effectiveness.We compare high-frequency waveforms computed with different regional velocity models against observations.Results show that the hybrid method achieves accuracy comparable to full SEM 3D simulations while reducing computation costs by more than two orders of magnitude when the structure of the source region updates.Our results further indicate that high-frequency waveforms are highly sensitive to shallow structures.Introducing low-velocity shallow layers into the source region improves near-field waveform fits,indicating pronounced low-velocity sediments in the Changning area.Large surface-wave time delays suggest that shallow velocities within the Sichuan Basin are lower than those in existing published models.In addition,an Interferometric Synthetic Aperture Radar(InSAR)-derived finite-fault model outperforms the point-source model in near-field waveform fitting and better reproduces rupture directivity.The proposed method is practical for high-frequency waveform modeling in areas with complex subsurface structures and rupture processes.
基金supported by the“National Ocean Technology Center Innovation Fund”under Project No.N3220Z002,led by Ning Jia.The official website of the National Ocean Technology Center is accessible at:http://www.notcsoa.org.cn/.
文摘With the depletion of fossil fuels and increasing environmental concerns,the development of renewable energy,such as wave energy,has become a critical component of global energy strategies.However,challenges persist in the field testing methodologies for wave energy converters(WECs).In this paper,a numerical wave field of the Dawanshan Island Sea Area in Zhuhai City is constructed based on the MIKE21 SW wave model and by using an NCEP wind field driving model.In conjunction with the IEC-62600-100 standard,by taking site testing of the“Wanshan”wave energy converter on which a sea trial has been conducted in Dawanshan Island of Zhuhai city as an example,research on-site testing method for a wave energy converter has been carried out.The wave measurement position for the“Wanshan”converter was determined by combining statistically analyzed field data with a validated numerical wave model.By comparing a valid wave height at the position where a wave rider is located with a valid wave height at the position where the“Wanshan”wave energy converter is situated,the correlation coefficient between simulation and observed data reached 0.90,with a root-mean-square error of 0.19.The representativeness of wave measurement data during site testing is verified and can be used as a basis for calculating the input energy of the“Wanshan”wave energy converter.
基金Supported by the National Natural Science Foundation of China(12261050)Science and Technology Project of Department of Education of Jiangxi Province(GJJ2201612 and GJJ211027)Natural Science Foundation of Jiangxi Province of China(20212BAB202021)。
文摘In this paper,we investigate the periodic traveling wave solutions problem for a single population model with advection and distributed delay.By the bifurcation analysis method,we can obtain periodic traveling wave solutions for this model under the influence of advection term and distributed delay.The obtained results indicate that weak kernel and strong kernel can both deduce the existence of periodic traveling wave solutions.Finally,we apply the main results in this paper to Logistic model and Nicholson’s blowflies model.
基金financially supported by the Scientific Research and Technology Development Project of China National Petroleum Corporation(No.2021ZG02)titled"Development of Seismic Data Processing Software for Ocean Nodes(OBN)"。
文摘In multi-component oil and gas exploration using ocean bottom nodes,converted wave data is rich in lithological and fracture information.One of the urgent problems to be solved is how to construct an accurate shear wave velocity model of the shallow sea bottom by leveraging the seismic wave information at the fluid-solid interface in the ocean,and improve the lateral resolution of marine converted wave data.Given that the dispersion characteristics of surface waves are sensitive to the S-wave velocity of subsurface media,and that Scholte surface waves,which propagate at the interface between liquid and solid media,exist in the data of marine oil and gas exploration,this paper proposes a Scholte wave inversion and modeling method based on oil and gas exploration using ocean bottom nodes.By using the method for calculating the Scholte wave dispersion spectrum based on the Bessel kernel function,the accuracy of dispersion spectrum analysis is improved,and more accurate dispersion curves are picked up.Through the adaptive weighted least squares Scholte wave dispersion inversion algorithm,the Scholte wave dispersion equation for liquid-solid media is solved,and the shear wave velocity model of the shallow sea bottom is calculated.Theoretical tests and applications of realdata have proven that this method can significantly improve the lateral resolution of converted wave data,provide high-quality data for subsequent inversion of marine multi-component oil and gas exploration data and reservoir reflection information,and contribute to the development of marine oil and gas exploration technology.
基金The Youth Talent Support Program for Doctoral Students of the China Association for Science and Technologythe Key Laboratory of Space Ocean Remote Sensing and Application at Ministry of Natural Resources under contract No.2023CFO005+1 种基金the National Natural Science Foundation of China under contract No.42176011the Fundamental Research Funds for the Central Universities under contract No.24CX03001A.
文摘This study applies Ensemble Optimal Interpolation(EnOI)to assimilate individual spectral components derived from National Data Buoy Center(NDBC)buoy directional spectra into the WAVEWATCHⅢ(WW3)wave model during tropical cyclone(TC)Isaias(2020).The analysis provides a comprehensive evaluation of the assimilation’s impact on wave parameters,frequency spectra,and directional spectra.Two series of assimilation experiments—one based on spectral components(Exp^(*)-DaSpec)and the other on significant wave height(Exp^(*)-DaSWH)—are evaluated against a non-assimilated control run(Exp-NoDa).Particular focus is placed on six key parameters:SWH,mean wave period(MWP),mean wave direction(MWD),mean wave directional spread(MWS),dominant wave period(DWP),and dominant wave direction(DWD).Sensitivity analyses suggest 400 km and 0.30 as appropriate values for the localization radius and observation error variance,respectively,though no single setting is optimal across all wave parameters.Exp4-DaSWH and Exp4-DaSpec are therefore selected as representative experiments.In non-independent validation,Exp4-DaSpec generally outperforms Exp-NoDa across MWP,MWD,MWS,DWP,and DWD,demonstrating closer agreement with observed frequency spectra and directional patterns.In independent validation,Exp4-DaSpec maintains superior overall performance,whereas Exp4-DaSWH shows only limited improvement.Exp4-DaSWH may capture the spectral peak near 0.1 Hz,although its directional characteristics remain largely similar to those of Exp-NoDa.Importantly,the assimilation experiments significantly improve SWH in both non-independent and independent validations,with Exp4-DaSWH slightly outperforming Exp4-DaSpec overall.
基金supported by the National Natural Science Foundation of China(Grant Nos.42475016,42192555 and 42305085)the China Postdoctoral Science Foundation(Grant No.2023M741615)the 2023 Graduate Research Innovation Project of Hunan Province(Grant No.CX20230011)。
文摘This paper investigates the impact of the model top and damping layer on the numerical simulation of tropical cyclones(TCs)and reveals the significant role of stratospheric gravity waves(SGWs).TCs can generate SGWs,which propagate upward and outward into the stratosphere.These SGWs can reach the damping layer,which is a consequence of the numerical scheme employed,where they can affect the tangential circulation through the dragging and forcing processes.In models with a higher top boundary,this tangential circulation develops far from the TC and has minimal direct impact on TC intensity.By comparison,in models with a lower top(e.g.,20 km),the damping layer is located just above the top of the TC.The SGW dragging in the damping layer and the consequent tangential force can thus induce ascent outside the eyewall,promote latent heat release,tilt the eyewall,and enlarge the inner-core radius.This process will reduce inner-core vorticity advection within the boundary layer,and eventually inhibits the intensification of the TC.This suggests that when the thickness of the damping layer is 5 km,the TC numerical model top height should be at least higher than 20 km to generate more accurate simulations.
基金supported by the National Key R&D Program of China (No.2022YFC3004602)the National Natural Science Foundation of China (No.52325404)the Shenzhen Science and Technology Program (No.JCYJ20220818095605012).
文摘Rock discontinuities such as joints widely exist in natural rock masses,and wave attenuation through rock masses is mainly caused by discontinuities.The displacement discontinuity model(DDM)has been widely used in theoretical and numerical analysis of wave propagation across rock discontinuity.However,the circumstance under which the DDM is applicable to predict wave propagation across rock discontinuity remains poorly understood.In this study,theoretical analysis and ultrasonic laboratory tests were carried out to examine the theoretical applicability of the DDM for wave propagation,where specimens with rough joints comprising regular rectangular asperities of different spacings and heights were prepared by 3D printing technology.It is found that the theoretical applicability of the DDM to predict wave propagation across rock discontinuity is determined by three joint parameters,i.e.the dimensionless asperity spacing(L),the dimensionless asperity height(H)and the groove density(D).Through theoretical analysis and laboratory tests,the conditions under which the DDM is applicable are derived as follows:and,.With increase in the groove density,the thresholds of the dimensionless asperity spacing and the dimensionless asperity height show a decreasing trend.In addition,the transmission coefficient in the frequency domain decreases with increasing groove density,dimensionless asperity spacing or dimensionless asperity height.The findings can facilitate our understanding of DDM for predicting wave propagation across rock discontinuity.
基金supported by JUST Start-up Fund for Science Research,the Jiangsu Natural Science Foundation(Grant No.BK20210885).
文摘The internal solitary wave(ISW)represents a frequent and severe oceanic dynamic phenomenon observed in the South China Sea,exposing marine structures to sudden loads.This paper examines the prediction model of interaction loads between ISW and FPSO,accounting for varying attack angles and incorporating ISW theories.The research demonstrates that the horizontal and transverse forces on FPSO under internal solitary waves(ISWs)comprise wave pressure difference force and viscous force,while the vertical force primarily consists of vertical wave pressure difference force.The wave pressure difference force is determined using the Froude-Krylov equation.The viscous force is derived from the tangential particle velocity induced by ISW and the viscous coefficient.The viscous coefficient formula is obtained through regression analysis of experimental data with different ISW attack angles.The research reveals that the horizontal viscous coefficient C_(vx)decreases as Reynolds number(R_(e))increases,while the transverse viscous coefficient C_(vy)initially increases and subsequently decreases with the growth of the Keulegan-Carpenter number(KC).Moreover,changes in wave propagation direction significantly affect the extreme magnitudes of both horizontal and transverse forces,and simultaneously modify the transverse force orientation,while having minimal impact on the vertical force.Additionally,the forces increase with the ISW’s amplitude.For horizontal and transverse forces,a thinner upper fluid layer generates larger forces.Comparative analysis of experimental,numerical,and theoretical results indicates strong agreement between theoretical predictions and experimental and numerical outcomes.
基金supported by the National Key R&D Program of China(No.2023YFA1606701)the National Natural Science Foundation of China(Nos.12175042,11890710,11890714,12047514,12147101,and 12347106)+1 种基金Guangdong Major Project of Basic and Applied Basic Research(No.2020B0301030008)China National Key R&D Program(No.2022YFA1602402).
文摘We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.
基金Supported by the Development and Application Project of Ship CAE Software.
文摘Studies of wave-current interactions are vital for the safe design of structures.Regular waves in the presence of uniform,linear shear,and quadratic shear currents are explored by the High-Level Green-Naghdi model in this paper.The five-point central difference method is used for spatial discretization,and the fourth-order Adams predictor-corrector scheme is employed for marching in time.The domain-decomposition method is applied for the wave-current generation and absorption.The effects of currents on the wave profile and velocity field are examined under two conditions:the same velocity of currents at the still-water level and the constant flow volume of currents.Wave profiles and velocity fields demonstrate substantial differences in three types of currents owing to the diverse vertical distribution of current velocity and vorticity.Then,loads on small-scale vertical cylinders subjected to regular waves and three types of background currents with the same flow volume are investigated.The maximum load intensity and load fluctuation amplitude in uniform,linear shear,and quadratic shear currents increase sequentially.The stretched superposition method overestimates the maximum load intensity and load fluctuation amplitude in opposing currents and underestimates these values in following currents.The stretched superposition method obtains a poor approximation for strong nonlinear waves,particularly in the case of the opposing quadratic shear current.
基金founded by the Ministry of Science and Higher Education of the Russian Federation,State assignments for research,registration No.1024032600084-8-1.3.2Study of the grain growth and the formation of polycrystalline structure as a result of phase transition(Section 6)was founded by the Russian Science Foundation,Project No.24-71-00078+3 种基金https://rscf.ru/en/project/24-71-00078/(accessed on 01 December 2025).Study of the orientation dependence of the phase transition of aluminum in Section 3 was founded by the Russian Science Foundation,Project No.24-19-00684https://rscf.ru/en/project/24-19-00684/(accessed on 01 December 2025).
文摘It is well known that aluminum and copper exhibit structural phase transformations in quasi-static and dynamic measurements,including shock wave loading.However,the dependence of phase transformations in a wide range of crystallographic directions of shock loading has not been revealed.In this work,we calculated the shock Hugoniot for aluminum and copper in different crystallographic directions([100],[110],[111],[112],[102],[114],[123],[134],[221]and[401])of shock compression using molecular dynamics(MD)simulations.The results showed a high pressure(>160 GPa for Cu and>40 GPa for Al)of the FCC-to-BCC transition.In copper,different characteristics of the phase transition are observed depending on the loading direction with the[100]compression direction being the weakest.The FCC-to-BCC transition for copper is in the range of 150–220 GPa,which is consistent with the existing experimental data.Due to the high transition pressure,the BCC phase transition in copper competes with melting.In aluminum,the FCC-to-BCC transition is observed for all studied directions at pressures between 40 and 50 GPa far beyond the melting.In all considered cases we observe the coexistence of HCP and BCC phases during the FCC-to-BCC transition,which is consistent with the experimental data and atomistic calculations;this HCP phase forms in the course of accompanying plastic deformation with dislocation activity in the parent FCC phase.The plasticity incipience is also anisotropic in bothmetals,which is due to the difference in the projections of stress on the slip plane for different orientations of the FCC crystal.MD modeling results demonstrate a strong dependence of the FCC-to-BCC transition on the crystallographic direction,in which the material is loaded in the copper crystals.However,MD simulations data can only be obtained for specific points in the stereographic direction space;therefore,for more comprehensive understanding of the phase transition process,a feed-forward neural network was trained using MD modeling data.The trained machine learning model allowed us to construct continuous stereographic maps of phase transitions as a function of stress in the shock-compressed state of metal.Due to appearance and growth of multiple centers of new phase,the FCC-to-BCC transition leads to formation of a polycrystalline structure from the parent single crystal.
基金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 Key Technology Program,the Ministry of Education of China under contract No.104061
文摘Based on the Longuet-Higgins wave model theory, the previews studies have shown that freak waves can be generated in finite space and time successfully. However, as to generating high nonlinear freak waves, the simulation results will be unrealistic. Therefore, a modified phase modulation method for simulating high nonlinear freak waves was developed. The surface elevations of some wave components at certain time and place are positive by modulating the corresponding random initial phases, then the total surface elevation at the focused point is enhanced and furthermore a freak wave event is generated. The new method can not only make the freak wave occur at certain time and place, but also make the simulated wave surface time series satisfy statistical properties of the realistic sea state and keep identical with the target wave spectrum. This numerical approach is of good precision and high efficiency by the comparisons of the simulated freak waves and the recorded freak waves.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
基金supported by the National Natural Science Foundation of China (Grant No.50779004)
文摘Four focusing models for generation of freak waves are presented. An extreme wave focusing model is presented on the basis of the enhanced High-Order Spectral （HOS） method and the importance of the nonlinear wave-wave interaction is evaluated by comparison of the calculated results with experimental and theoretical data. Based on the modification of the Longuet-Higgins model, four wave models for generation of freak waves （a. extreme wave model ＋ random wave model; b. extreme wave model ＋ regular wave model; e. phase interval modulation wave focusing model; d. number modulation wave focusing model with the same phase） are proposed. By use of different energy distribution techniques in the four models, freak wave events are obtained with different Hmax/Hs in finite space and time.
文摘High-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth with an arbitrary sloping bottom are presented in this article. First, the formal derivations to any high order of mu(= h/lambda, depth to deep-water wave length ratio) and epsilon(= a/h, wave amplitude to depth ratio) for velocity potential, particle velocity vector, pressure and the Boussinesq-type equations for surface elevation eta and horizontal velocity vector (U) over right arrow at any given level in water are given. Then, the exact explicit expressions to the fourth order of mu are derived. Finally, the linear solutions of eta, (U) over right arrow, C (phase-celerity) and C-g (group velocity) for a constant water depth are obtained. Compared with the Airy theory, excellent results can be found even for a water depth as large as the wave legnth. The present high-order models are applicable to nonlinear regular and irregular waves in water of any varying depth (from shallow to deep) and bottom slope (from mild to steep).
文摘The LAGFD-WAM wave model is a third generation wave model. In the present paper the physical aspect of the model was shown in great detail including energy spectrum balance equation, complicated characteristics equations and source functions.
基金This project is supported bythe Canadian Panel on Energy Research and Development (Offshore Environmental Fac-tors Program) , ONR (US Office of Naval Research) via GoMOOS-the Gulf of Maine Ocean Observing System,Petroleum Research Atlantic Canada (PRAC) ,and the CFCAS (Canada Foundation for Climate and AtmosphericStudies) ,Canadian Panel on Energy Research and Development (Offshore Environmental Factors Program) .It is al-so supported bythe Advanced Doctoral Fund of the Ministry of Education of China (Grant No.20030294010)
文摘Typhoon-generated waves are simulated with two numerical wave models, the SWAN model for the coastal and Yangtze Estuary domain, nested within the WAVEWATCHIII （WW3） for the basin-scale East China Sea domain. Typhoon No. 8114 is chosen because it was very strong, and generated high waves in the Estuary. WW3 was implemented for the East China Sea coarse-resolution computational domain, to simulate the waves over a large spatial scale and provide boundary conditions for SWAN model simulations, implemented on a fine-resolution nested domain for the Yangtze Estuary area. The Takahashi wind model is applied to the simulation of the East China Sea scale （3-hourly） and Yangtze Estuary scale （1-hourly） winds. Simulations of significant wave heights in the East China Sea show that the highest waves are on the right side of the storm track, and maxima tend to occur at the eastern deep-water open boundary of the Yangtze Estuary. In the Yangtze Estuary, incoming swell is dominant over locally generated waves before the typhoon approaches the Estuary. As the typhoon approaches the Estuary, wind waves and swell coexist, and the wave direction is mainly influenced by the swell direction and the complex topography.
基金This work was supported bythe National Basic Research Program(973) of China (Grant No.2003CB415206) andthe National Natural Science Foundation of China (Grant No.50379027 and No.50479004)
文摘In this study, the combined actions of waves and tidal currents in estuarine and coastal areas are considered and a 2D mathematical model for sediment transport by waves and tidal currents has been established in orthogonal curvilinear coordinates. Non-equilibrium transport equations of suspended load and bed load are used in the model. The concept of background concentration is introduced, and the formula of sediment transport capacity of tidal currents for the Oujiang River estuary is obtained. The Dou Guoren formula is employed for the sediment transport capacity of waves. Sediment transport capacity in the form of mud and the intensity of back silting are calculated by use of Luo Zaosen＇ s formula. The calculated tidal stages are in good agreement with the field data, and the calculated velocities and flow directions of 46 vertical lines for 8 cross sections are also in good agreement with the measured data. On such a basis, simulations of back silting after excavation of the waterway with a sand bar under complicated boundary conditions in the navigation channel induced by suspended load, bed load and mud by waves and tidal currents are discussed.
