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.展开更多
In this paper the parameterizational approach of nonlinear source function and the implicit scheme of the model are discussed in detail. The matching problem is solved between time and space steps using the characteri...In this paper the parameterizational approach of nonlinear source function and the implicit scheme of the model are discussed in detail. The matching problem is solved between time and space steps using the characteristics inlaid scheme with very strong physical meaning. The computational comparison in typical winds shows some improvements to the WAM model. That the hindcast results of the model for typhoon cases are in good agreement with real data illustrates its applicability to wave forecast and engineering study.展开更多
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.展开更多
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.展开更多
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.展开更多
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).展开更多
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.展开更多
In a probe and drogue aerial refueling system, the bow wave of the receiver aircraft will produce a strong aerodynamic effect on the drogue once the receiver follows the drogue at a close distance. It is a major diffi...In a probe and drogue aerial refueling system, the bow wave of the receiver aircraft will produce a strong aerodynamic effect on the drogue once the receiver follows the drogue at a close distance. It is a major difficulty of docking control in the probe and drogue refueling. This paper analyses the bow wave effect and presents a simple method to model it. Firstly, the inviscid flow around the receiver is modeled based on the stream function defined by basic stream singularities. Secondly, a correction function is developed to eliminate the error caused by the absence of air vis- cosity. Then, the aerodynamic coefficients are used to calculate the induced aerodynamic force on the drogue. The obtained model is in an analytical form that can be easily applied to the controller design and the real-time simulations. In the verification part, computational fluid dynamics （CFD） simulation tests are conducted to validate the obtained flow fields and aerodynamic forces. Finally, the modeling method is applied to an F-16 receiver aircraft in a previously developed autonomous aerial refueling simulation system. The simulations results are analyzed and compared with the NASA flight-test data, which demonstrates the effectiveness of the proposed method.展开更多
This paper pnesents a third gneration shallow Whter disode spedtal wave nbotal medeIYE-WAM based on the spedtal action balance equation. The mode accounts for all edevan effectsof currents on waves, incuding tmpotally...This paper pnesents a third gneration shallow Whter disode spedtal wave nbotal medeIYE-WAM based on the spedtal action balance equation. The mode accounts for all edevan effectsof currents on waves, incuding tmpotally and spatialy varying depth and current inded refraction,sttalning and fequency shift and also explidtly takeS into aanunt all source terms, speclally adePth-limited breaking dheipation. In addition, an energy forcing scheme is propond and applied to themode’s open boundaries to areUn for the propagution of sedIs into the study spstem The upwinddiffeIenng scheme and a standard hybrid diffdrencing scheme for the propagaion terrn and a simpleEuler method for the source teme are employed.展开更多
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.展开更多
The history of forecasting wind waves by wave energy conservation equation Is briefly described. Several currently used wave numerical models for shallow water based on different wave theories are discussed. Wave ener...The history of forecasting wind waves by wave energy conservation equation Is briefly described. Several currently used wave numerical models for shallow water based on different wave theories are discussed. Wave energy conservation models for the simulation of shallow water waves are introduced, with emphasis placed on the SWAN model, which takes use of the most advanced wave research achievements and has been applied to several theoretical and field conditions. The characteristics and applicability of the model, the finite difference numerical scheme of the action balance equation and its source terms computing methods are described in detail. The model has been verified with the propagation refraction numerical experiments for waves propagating in following and opposing currents; finally, the model is applied to the Haian Gulf area to simulate the wave height and wave period field there, and the results are compared with observed data.展开更多
文摘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.
文摘In this paper the parameterizational approach of nonlinear source function and the implicit scheme of the model are discussed in detail. The matching problem is solved between time and space steps using the characteristics inlaid scheme with very strong physical meaning. The computational comparison in typical winds shows some improvements to the WAM model. That the hindcast results of the model for typhoon cases are in good agreement with real data illustrates its applicability to wave forecast and engineering study.
基金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.
基金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.
基金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.
文摘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).
基金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.
基金supported by the National Natural Science Foundation of China(Nos.61473012 and 51375462)
文摘In a probe and drogue aerial refueling system, the bow wave of the receiver aircraft will produce a strong aerodynamic effect on the drogue once the receiver follows the drogue at a close distance. It is a major difficulty of docking control in the probe and drogue refueling. This paper analyses the bow wave effect and presents a simple method to model it. Firstly, the inviscid flow around the receiver is modeled based on the stream function defined by basic stream singularities. Secondly, a correction function is developed to eliminate the error caused by the absence of air vis- cosity. Then, the aerodynamic coefficients are used to calculate the induced aerodynamic force on the drogue. The obtained model is in an analytical form that can be easily applied to the controller design and the real-time simulations. In the verification part, computational fluid dynamics （CFD） simulation tests are conducted to validate the obtained flow fields and aerodynamic forces. Finally, the modeling method is applied to an F-16 receiver aircraft in a previously developed autonomous aerial refueling simulation system. The simulations results are analyzed and compared with the NASA flight-test data, which demonstrates the effectiveness of the proposed method.
基金Supported by the National Eighty-Five-Year Project D09920109 and Chinese Academy of Sciences and State Education Commission
文摘This paper pnesents a third gneration shallow Whter disode spedtal wave nbotal medeIYE-WAM based on the spedtal action balance equation. The mode accounts for all edevan effectsof currents on waves, incuding tmpotally and spatialy varying depth and current inded refraction,sttalning and fequency shift and also explidtly takeS into aanunt all source terms, speclally adePth-limited breaking dheipation. In addition, an energy forcing scheme is propond and applied to themode’s open boundaries to areUn for the propagution of sedIs into the study spstem The upwinddiffeIenng scheme and a standard hybrid diffdrencing scheme for the propagaion terrn and a simpleEuler method for the source teme are employed.
基金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.
基金"333"Project Scientific Research Foundation of Jiangsu ProvinceScience Fundation of Hohai University(3853)
文摘The history of forecasting wind waves by wave energy conservation equation Is briefly described. Several currently used wave numerical models for shallow water based on different wave theories are discussed. Wave energy conservation models for the simulation of shallow water waves are introduced, with emphasis placed on the SWAN model, which takes use of the most advanced wave research achievements and has been applied to several theoretical and field conditions. The characteristics and applicability of the model, the finite difference numerical scheme of the action balance equation and its source terms computing methods are described in detail. The model has been verified with the propagation refraction numerical experiments for waves propagating in following and opposing currents; finally, the model is applied to the Haian Gulf area to simulate the wave height and wave period field there, and the results are compared with observed data.
