This study proposes a novel cyclone separator with a conical inner core to enhance particle classification efficiency in oil and gas wellhead-recovered liquids.Particle motion and force dynamics are analyzed to optimi...This study proposes a novel cyclone separator with a conical inner core to enhance particle classification efficiency in oil and gas wellhead-recovered liquids.Particle motion and force dynamics are analyzed to optimize key structural parameters,including inlet diameter(D_i),overflow pipe diameter(D_(e)),insertion depth(L_(e)),and bottom flow pipe diameter(D_(z)).Numerical simulations employ the Reynolds stress turbulence model,SIMPLEC algorithm,and discrete phase model to evaluate separation performance in a gas-liquid two-phase system.Results indicate that a smaller D_i improves fine particle separation but increases turbulence;an optimal range of D_i/D_(c)=0.35-0.4 is recommended.Larger D_(e) enhances the diversion ratio,aiding fine particle discharge(D_(e)/D_(c)=0.25-0.35).Increased Le facilitates fine particle overflow but induces vortices,whereas a smaller L_(e) stabilizes the bottom flow for larger particle separation(L_(e)/D_(c)=0.5-0.75).A reduced D_(z) enhances centrifugal force and separation efficiency but may cause turbulence;an optimal D_(z)/D_(c) of 0.6-0.65 is suggested for stability.These findings provide valuable design guidelines for improving cyclone separator performance in multiphase flow applications.展开更多
Corresponding to the sliding and the overturning failure,the elementary motion modes of caisson breakwater include the horizontal-rotational oscillation coupled motion,the horizontal sliding-rotational oscillation cou...Corresponding to the sliding and the overturning failure,the elementary motion modes of caisson breakwater include the horizontal-rotational oscillation coupled motion,the horizontal sliding-rotational oscillation coupled motion,the horizontal vibrating-uplift rocking coupled motion,and the horizontal sliding-uplift rocking coupled motion.The motion mode of a caisson will transform from one to another depending on the wave forces and the motion behaviors of the caisson.The numerical models of four motion modes of caisson are developed,and the numerical simulation procedure for joint motion process of various modes of caisson breakwater under wave excitation is presented and tested by a physical model experiment.It is concluded that the simulation procedure is reliable and can be applied to the dynamic stability analysis of caisson breakwaters.展开更多
In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform ...In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.展开更多
The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems wit...The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.展开更多
We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficienc...We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.展开更多
The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investig...The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.展开更多
Waves are important physical phenomena in an ocean,and their accurate prediction is essential for ocean engineering,maritime traffic,and marine early warning systems.This study focuses on the Qinhuangdao Sea area loca...Waves are important physical phenomena in an ocean,and their accurate prediction is essential for ocean engineering,maritime traffic,and marine early warning systems.This study focuses on the Qinhuangdao Sea area located in the Bohai Sea,China.Herein,we use on-site wind data to correct the reanalysis wind data obtained from the European Centre for Medium-Range Weather Forecasts(ECMWF),improving the accuracy of boundary conditions.Then,we use the Simulating WAves Nearshore(SWAN)model to simulate the regional wave field over time.A regional wave-parameter prediction model is then developed using a limited number of sampled data(covering only 2 years,2020–2021);the model is based on the Whale Optimization Algorithm(WOA),convolutional neural networks(CNNs),and long short-term memory(LSTM)neural networks.WOA is used to optimize the CNN and LSTM framework;in this framework,CNN extracts spatial features,and the LSTM network captures temporal features,enabling accurate short and long-term predictions of wave height,period,and direction.The experimental results showed that despite the small sample size,the model achieves a goodness of fit of 0.9957 for wave height prediction,0.9973 for period,and 0.9749 for wave direction in short-term forecasting.As the prediction step size increases,the accuracy of the model decreases.When the prediction step size reaches 9 h,the root mean square error for the prediction of wave height,period,and direction increases to 0.2060 m,0.4582 s,and32.5358°,respectively.The reliability and applicability of the model are further validated by the experimental results.Our findings highlighted the potential of the developed model in operational wave forecasting,even with a limited number of sampled data.展开更多
The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismi...The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismic wave simulation with a data-constrained finite-fault rupture model.The constraint is implemented by identifying the optimal ground motion models(GMMs)through a scoring system that selects the best-fit GMMs to mid-and far-field China Earthquake Networks Center(CENC)seismic network data;and applying the optimal GMMs to refine the rupture model parameters for near-fault intensity field simulation.The simulated near-fault seismic intensity field reproduces seismic intensities collected from Myanmar’s sparse seismic network and concentrated in≥Ⅷintensity zones within 50 km of the projected fault plane;and identifies abnormal intensity regions exhibiting≥Ⅹintensity along the Meiktila-Naypyidaw corridor and near Shwebo that are attributed to soft soil amplification effects and near-fault directivity.This framework can also be applied to post-earthquake assessments in other similar regions.展开更多
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.展开更多
The paper considers a catastrophic event-the eruption of Hunga Tonga-Hunga Ha’apai volcano on January 15,2022.The process of preparation and eruption of Hunga Tonga volcano generated tsunami waves that were observed ...The paper considers a catastrophic event-the eruption of Hunga Tonga-Hunga Ha’apai volcano on January 15,2022.The process of preparation and eruption of Hunga Tonga volcano generated tsunami waves that were observed throughout the World Ocean.This event was notable for its unprecedented global impact and the early appearance of tsunami waves at distant coastal stations.So,the first waves at tide gauge stations in Chile and Peru were recorded 4 hours earlier than the arrival time of tsunami waves to the tide gauge after the eruption of Tonga volcano.Two mechanisms are possible for the generation of early tsunami waves:acoustic Lamb waves generated by a volcanic explosion and submarine landslides that occurred on the slopes of the volcano during the preparatory phase of the eruption.In this study,numerical simulation of various pre-eruption landslide scenarios on the slope of Hunga Tonga volcano is carried out in an attempt to explain these early tsunami waves.Under computation the elastoplastic model of landslide was taken into account.Wave characteristics of a tsunami on the coast of Chile and Peru generated by a landslide process on a volcanic slope are obtained.A detailed comparison of virtual tide gauge data with observational ones is used to validate this model.The results obtained can be used to improve early warning systems.展开更多
This study employs the Smoothed Particle Hydrodynamics(SPH)method to develop a computational fluid dynamics(CFD)model for analyzing the interaction between rogue waves and mooring systems.Four floating body configurat...This study employs the Smoothed Particle Hydrodynamics(SPH)method to develop a computational fluid dynamics(CFD)model for analyzing the interaction between rogue waves and mooring systems.Four floating body configurations are investigated:(1)dual rectangular prisms,(2)rectangular prism–sphere composites,(3)sphere–rectangular prism composites,and(4)dual spheres.These configurations are systematically evaluated under varying mooring conditions to assess their hydrodynamic performance and wave attenuation capabilities.The model accurately captures the complex fluid–structure interaction dynamics between moored floating breakwaters and incident wave fields.Among the configurations,the dual rectangular prism system demonstrates superior performance in both wave dissipation and mooring force reduction.Under conditions involving dual wave makers,the influence of floating body shape and number on wave height is found to be minimal.However,dual-body arrangements consistently outperform single-body setups in terms of both energy dissipation and structural stability.From a cost-efficiency perspective,the configuration comprising two rectangular prisms connected via a single mooring system offers significant advantages in material usage and deployment feasibility.展开更多
The impact force effect on launch platform motion response represents a critical safety consideration that requires thorough investigation prior to sea-launch implementation. This paper examines a self-designed semi-s...The impact force effect on launch platform motion response represents a critical safety consideration that requires thorough investigation prior to sea-launch implementation. This paper examines a self-designed semi-submersible launch platform comprising a box-shaped deck, six columns, and two pontoons, with sufficient structural stiffness to be analyzed as a rigid body. A proprietary code based on three-dimensional linear potential theory was developed for hydrodynamic analysis of the launching process. The Cummins equation was implemented to calculate platform responses under substantial impact force. The numerical results were validated through comparison with ANSYS/Aqwa commercial software for platform motion response in both launch and non-launch cases. Additionally, two model tests were conducted in a sea-keeping wave basin at a scale ratio of λ =1:40. The numerical results demonstrated good agreement with experimental data. Both numerical and experimental findings indicate that platform motion responses result from wave-induced effects and impact force/rocket weight effects, with the latter typically predominant. Numerical simulations revealed that in moderate sea states, maximum heave and pitch motions measure 0.6 m and 1°, respectively, suggesting the viability of sea-launch operations using the designed platform under these conditions.展开更多
Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic th...Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.展开更多
The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general condit...The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.展开更多
A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational r...A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.展开更多
The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is establishe...The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.展开更多
Numerical method is popular in analysing the blast wave propagation and interaction with structures.However,because of the extremely short duration of blast wave and energy trans-mission between different grids,the nu...Numerical method is popular in analysing the blast wave propagation and interaction with structures.However,because of the extremely short duration of blast wave and energy trans-mission between different grids,the numerical results are sensitive to the finite element mesh size.Previous numerical simulations show that a mesh size acceptable to one blast scenario might not be proper for another case,even though the difference between the two scenarios is very small,indicating a simple numerical mesh size convergence test might not be enough to guarantee accu-rate numerical results.Therefore,both coarse mesh and fine mesh were used in different blast scenarios to investigate the mesh size effect on numerical results of blast wave propagation and interaction with structures.Based on the numerical results and their comparison with field test re-sults and the design charts in TM5-1300,a numerical modification method was proposed to correct the influence of the mesh size on the simulated results.It can be easily used to improve the accu-racy of the numerical results of blast wave propagation and blast loads on structures.展开更多
Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately ...Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.展开更多
The simulating waves nearshore(SWAN) model has typically been designed for wave simulations in near-shore regions. In this study, the model's applicability to the simulation of typhoon waves in the South China Sea...The simulating waves nearshore(SWAN) model has typically been designed for wave simulations in near-shore regions. In this study, the model's applicability to the simulation of typhoon waves in the South China Sea(SCS) was evaluated. A blended wind field, consisting of an interior domain based on Fujita's model and an exterior domain based on Takahashi's model, was used as the driving wind field. The waves driven by Typhoon Kai-tak over the SCS that occurred in 2012 were selected for the numerical simulation research. Sensitivity analyses of time step, grid resolution, and angle resolution were performed in order to obtain optimal model settings. Through sensitivity analyses, it can be found that the time step has a large influence on the results, while grid resolution and angle resolution have a little effect on the results.展开更多
A numerical model of flatfish cage is built based on the lumped mass method and the principle of rigid body kinematics. To validate the numerical model, a series of physical model tests are conducted in the wave flume...A numerical model of flatfish cage is built based on the lumped mass method and the principle of rigid body kinematics. To validate the numerical model, a series of physical model tests are conducted in the wave flume. The numerical results correspond well with the data sets from physical model test. The effect of weight of bottom frame, height of fish net and net shape on motion responses of fish cage and tension force on mooring lines is then analyzed. The results indicate that the vertical displacements of float collar and bottom frame decrease with the increase in the weight of bottom frame; the maximum tension force on mooring lines increases with the increasing weight of bottom frame. The inclination angles of float collar and bottom frame decrease with the increasing net height; the maximum tension force increases obviously with the increase of net height.展开更多
基金supported by the National Natural Science Foundation of China(52074341)。
文摘This study proposes a novel cyclone separator with a conical inner core to enhance particle classification efficiency in oil and gas wellhead-recovered liquids.Particle motion and force dynamics are analyzed to optimize key structural parameters,including inlet diameter(D_i),overflow pipe diameter(D_(e)),insertion depth(L_(e)),and bottom flow pipe diameter(D_(z)).Numerical simulations employ the Reynolds stress turbulence model,SIMPLEC algorithm,and discrete phase model to evaluate separation performance in a gas-liquid two-phase system.Results indicate that a smaller D_i improves fine particle separation but increases turbulence;an optimal range of D_i/D_(c)=0.35-0.4 is recommended.Larger D_(e) enhances the diversion ratio,aiding fine particle discharge(D_(e)/D_(c)=0.25-0.35).Increased Le facilitates fine particle overflow but induces vortices,whereas a smaller L_(e) stabilizes the bottom flow for larger particle separation(L_(e)/D_(c)=0.5-0.75).A reduced D_(z) enhances centrifugal force and separation efficiency but may cause turbulence;an optimal D_(z)/D_(c) of 0.6-0.65 is suggested for stability.These findings provide valuable design guidelines for improving cyclone separator performance in multiphase flow applications.
基金supported by the National Natural Science Foundation of China(Grant No.50979069)the Science and Technology Project of West China Traffic Construction(Grant No.200632800003-06)
文摘Corresponding to the sliding and the overturning failure,the elementary motion modes of caisson breakwater include the horizontal-rotational oscillation coupled motion,the horizontal sliding-rotational oscillation coupled motion,the horizontal vibrating-uplift rocking coupled motion,and the horizontal sliding-uplift rocking coupled motion.The motion mode of a caisson will transform from one to another depending on the wave forces and the motion behaviors of the caisson.The numerical models of four motion modes of caisson are developed,and the numerical simulation procedure for joint motion process of various modes of caisson breakwater under wave excitation is presented and tested by a physical model experiment.It is concluded that the simulation procedure is reliable and can be applied to the dynamic stability analysis of caisson breakwaters.
基金China Postdoctoral Science Foundation Under Grant No.20100480321National Basic Research Program of China Under Grant No. 2007CB714200
文摘In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.
基金supported by the National Natural Science Foundation of China(Grant Nos.52271278 and 52111530137)the Natural Science Found of Jiangsu Province(Grant No.BK20221389)the Newton Advanced Fellowships(Grant No.NAF\R1\180304)by the Royal Society.
文摘The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.
基金Open Fund of State Key Laboratory of Geological Processes and Mineral Resources, China University of Geo-sciences (GPMR0750)National Natural Science Foundation of China (40437018)
文摘We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.
基金Supported by the National Natural Science Foundation of China(10672052)the Advanced TalentStart-Up Foundation of Jiangsu University(08JDG018)~~
文摘The spatial growth of the disturbance in the boundary layer is directly numerically simulated, and the receptivity of the Blasius basic flow to the local two-dimensional （2-D） sustainable micro-vibration is investigated. Results show that the disturbance velocity presents the sine vibration features with the change of time, and the vibration period is identical to the vibration of the local wall. The disturbance velocity presents the fluctuation feature downstream, and the streamwise wave length approximates to the results from the Orr-Sommerfeld equation （OSE）. The growth rate from direct numerical simulation（DNS） is a little greater than that from OSE, and their trends are almost consistent. Under the condition of Re= 2 800, the disturbance amplitude gradually grows in the given computational region with the period T=30. However, it firstly increases and then decreases with the period T= 20. The disturbance harmonic of the former is obviously larger than that of the latter. The maximum streamwise and vertical disturbance velocities from DNS do not fully coincide with those from OSE at the vicinity of the local vibration wall, but coincide well with the former when they travel downstream. The 2-D disturbance induced by the local micro-vibration represents the form of Tollmien-Schlichting （T-S） wave on the boundary layer.
基金supported by the National Natural Science Foundation of China(Nos.52071057,52171247)the Liaoning Youth Elite Talent Program(No.XLYC220309)。
文摘Waves are important physical phenomena in an ocean,and their accurate prediction is essential for ocean engineering,maritime traffic,and marine early warning systems.This study focuses on the Qinhuangdao Sea area located in the Bohai Sea,China.Herein,we use on-site wind data to correct the reanalysis wind data obtained from the European Centre for Medium-Range Weather Forecasts(ECMWF),improving the accuracy of boundary conditions.Then,we use the Simulating WAves Nearshore(SWAN)model to simulate the regional wave field over time.A regional wave-parameter prediction model is then developed using a limited number of sampled data(covering only 2 years,2020–2021);the model is based on the Whale Optimization Algorithm(WOA),convolutional neural networks(CNNs),and long short-term memory(LSTM)neural networks.WOA is used to optimize the CNN and LSTM framework;in this framework,CNN extracts spatial features,and the LSTM network captures temporal features,enabling accurate short and long-term predictions of wave height,period,and direction.The experimental results showed that despite the small sample size,the model achieves a goodness of fit of 0.9957 for wave height prediction,0.9973 for period,and 0.9749 for wave direction in short-term forecasting.As the prediction step size increases,the accuracy of the model decreases.When the prediction step size reaches 9 h,the root mean square error for the prediction of wave height,period,and direction increases to 0.2060 m,0.4582 s,and32.5358°,respectively.The reliability and applicability of the model are further validated by the experimental results.Our findings highlighted the potential of the developed model in operational wave forecasting,even with a limited number of sampled data.
基金Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration under Grant No.2023C01National Natural Science Foundation of China under Grant No.52478570Distinguished Young Scholars Program of the Natural Science Foundation of Heilongjiang Province,China under Grant No.JQ2024E002。
文摘The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismic wave simulation with a data-constrained finite-fault rupture model.The constraint is implemented by identifying the optimal ground motion models(GMMs)through a scoring system that selects the best-fit GMMs to mid-and far-field China Earthquake Networks Center(CENC)seismic network data;and applying the optimal GMMs to refine the rupture model parameters for near-fault intensity field simulation.The simulated near-fault seismic intensity field reproduces seismic intensities collected from Myanmar’s sparse seismic network and concentrated in≥Ⅷintensity zones within 50 km of the projected fault plane;and identifies abnormal intensity regions exhibiting≥Ⅹintensity along the Meiktila-Naypyidaw corridor and near Shwebo that are attributed to soft soil amplification effects and near-fault directivity.This framework can also be applied to post-earthquake assessments in other similar regions.
基金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.
基金funded by the Ministry of Science and Higher Education of the Russian Federation(agreement No.075-15-2022-1127 dated July 1,2022).
文摘The paper considers a catastrophic event-the eruption of Hunga Tonga-Hunga Ha’apai volcano on January 15,2022.The process of preparation and eruption of Hunga Tonga volcano generated tsunami waves that were observed throughout the World Ocean.This event was notable for its unprecedented global impact and the early appearance of tsunami waves at distant coastal stations.So,the first waves at tide gauge stations in Chile and Peru were recorded 4 hours earlier than the arrival time of tsunami waves to the tide gauge after the eruption of Tonga volcano.Two mechanisms are possible for the generation of early tsunami waves:acoustic Lamb waves generated by a volcanic explosion and submarine landslides that occurred on the slopes of the volcano during the preparatory phase of the eruption.In this study,numerical simulation of various pre-eruption landslide scenarios on the slope of Hunga Tonga volcano is carried out in an attempt to explain these early tsunami waves.Under computation the elastoplastic model of landslide was taken into account.Wave characteristics of a tsunami on the coast of Chile and Peru generated by a landslide process on a volcanic slope are obtained.A detailed comparison of virtual tide gauge data with observational ones is used to validate this model.The results obtained can be used to improve early warning systems.
基金funding from the National Natural Science Foundation of China(No.12462028).
文摘This study employs the Smoothed Particle Hydrodynamics(SPH)method to develop a computational fluid dynamics(CFD)model for analyzing the interaction between rogue waves and mooring systems.Four floating body configurations are investigated:(1)dual rectangular prisms,(2)rectangular prism–sphere composites,(3)sphere–rectangular prism composites,and(4)dual spheres.These configurations are systematically evaluated under varying mooring conditions to assess their hydrodynamic performance and wave attenuation capabilities.The model accurately captures the complex fluid–structure interaction dynamics between moored floating breakwaters and incident wave fields.Among the configurations,the dual rectangular prism system demonstrates superior performance in both wave dissipation and mooring force reduction.Under conditions involving dual wave makers,the influence of floating body shape and number on wave height is found to be minimal.However,dual-body arrangements consistently outperform single-body setups in terms of both energy dissipation and structural stability.From a cost-efficiency perspective,the configuration comprising two rectangular prisms connected via a single mooring system offers significant advantages in material usage and deployment feasibility.
基金financially supported by the National Natural Science Foundation of China (Grant No.52401346)the Natural Science Foundation of Jiangsu Province of China (Grant No.BK20230183)the Taihu Lake Talent Plan of Wuxi。
文摘The impact force effect on launch platform motion response represents a critical safety consideration that requires thorough investigation prior to sea-launch implementation. This paper examines a self-designed semi-submersible launch platform comprising a box-shaped deck, six columns, and two pontoons, with sufficient structural stiffness to be analyzed as a rigid body. A proprietary code based on three-dimensional linear potential theory was developed for hydrodynamic analysis of the launching process. The Cummins equation was implemented to calculate platform responses under substantial impact force. The numerical results were validated through comparison with ANSYS/Aqwa commercial software for platform motion response in both launch and non-launch cases. Additionally, two model tests were conducted in a sea-keeping wave basin at a scale ratio of λ =1:40. The numerical results demonstrated good agreement with experimental data. Both numerical and experimental findings indicate that platform motion responses result from wave-induced effects and impact force/rocket weight effects, with the latter typically predominant. Numerical simulations revealed that in moderate sea states, maximum heave and pitch motions measure 0.6 m and 1°, respectively, suggesting the viability of sea-launch operations using the designed platform under these conditions.
基金The National Key R&D Program of China under contract No.2023YFE0126300the National Natural Science Foundation of China under contract Nos 42066002 and U20A2099.
文摘Swells are critical concerns regarding safety,marine transportation,and coastal engineering construction of coastal countries along the Gulf of Guinea and have been scientific problems due to the lack of systematic theoretical,numerical,and observational research.In this study,a double nesting numerical model was constructed and validated from the Atlantic Ocean to the Gulf of Guinea based on simulating waves nearshore(SWAN)to explore the swell characteristics and source tracing in the Gulf of Guinea in winter and summer seasons from 2020 to 2021.Simulation results reveal that swells are stronger and deflect more to the west in winter than summer,even though they dominate in both seasons in the Gulf of Guinea in the S-SW directional range.Simulated two-dimensional(2D)wave spectral patterns not only clarify wave composition,variation,and propagation properties from the central South Atlantic Ocean to the Gulf of Guinea,but also distinguish swell strength and directional range in winter and summer.The NW wind events induce swells which spread toward the SSE-ESE direction from the North Atlantic Ocean,big wind source generates sustained and stable S-SW swells from the South Atlantic Ocean,and corresponding swell-influenced areas are discussed.The strongest swell event in the Gulf of Guinea during the simulation was used as a case study to trace its source.A strong clockwise wind vortex within the Roaring Forties induced these large swells in the Gulf of Guinea approximately 5.5 days later,and swell propagation formed a regular isoline of peak period distribution from the South Atlantic Ocean to the Gulf of Guinea in the SSW-SW direction.
文摘The numerical simulation is based on the authors' high-order models with a dissipative term for nonlinear and dispersive wave in water of varying depth. Corresponding finite-difference equations and general conditions for open and fixed natural boundaries with an arbitrary reflection coefficient and phase shift are also given in this paper. The systematical tests of numerical simulation show that the theoretical models, the finite-difference algorithms and the boundary conditions can give good calculation results for the wave propagating in shallow and deep water with an arbitrary slope varying from gentle to steep.
基金Trans-Century Training program Fund for the Talent,Ministry of Education of China
文摘A composite model, which is the combination of Boussinesq equations and Volume of Fluid (VOF) method, has been developed for 2-D time-domain computations of nonlinear waves in a large region. The whole computational region Omega is divided into two subregions. In the near-field around a structure, Omega(2), the flow is governed by 2-D Reynolds Averaged Navier-Stokes equations with a turbulence closure model of k-epsilon equations and numerically solved by the improved VOF method; whereas in the subregion Omega(1) (Omega(1) = Omega - Omega(2)) the flow is governed by one-D Boussinesq equations and numerically solved with the predictor-corrector algorithm. The velocity and the wave surface elevation are matched on the common boundary of the two subregions. Numerical tests have been conducted for the case of wave propagation and interaction with a wave barrier. It is shown that the composite model can help perform efficient computation of nonlinear waves in a large region with the complicated flow fields near structures taken into account.
基金funded by the National Natural Science Foundation of China under NO.11202103Qing-lan Project of Jiangsu Province。
文摘The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.
基金Supported by National Natural Science Foundation of China (No.50638030, 50528808)the National Key Technologies R&D Program of China (No.2006BAJ13B02)the Australian Research Council (No.DP0774061).
文摘Numerical method is popular in analysing the blast wave propagation and interaction with structures.However,because of the extremely short duration of blast wave and energy trans-mission between different grids,the numerical results are sensitive to the finite element mesh size.Previous numerical simulations show that a mesh size acceptable to one blast scenario might not be proper for another case,even though the difference between the two scenarios is very small,indicating a simple numerical mesh size convergence test might not be enough to guarantee accu-rate numerical results.Therefore,both coarse mesh and fine mesh were used in different blast scenarios to investigate the mesh size effect on numerical results of blast wave propagation and interaction with structures.Based on the numerical results and their comparison with field test re-sults and the design charts in TM5-1300,a numerical modification method was proposed to correct the influence of the mesh size on the simulated results.It can be easily used to improve the accu-racy of the numerical results of blast wave propagation and blast loads on structures.
基金supported by National Natural Science Foundation of China(Nos.41204077,41372290,41572244,51034003,51174210,and 51304126)natural science foundation of Shandong Province(Nos.ZR2011EEZ002 and ZR2013EEQ019)State Key Research Development Program of China(No.2016YFC0600708-3)
文摘Currently, numerical simulations of seismic channel waves for the advance detection of geological structures in coal mine roadways focus mainly on modeling two- dimensional wave fields and therefore cannot accurately simulate three-dimensional （3-D） full-wave fields or seismic records in a full-space observation system. In this study, we use the first-order velocity-stress staggered-grid finite difference algorithm to simulate 3-D full-wave fields with P-wave sources in front of coal mine roadways. We determine the three components of velocity Vx, Vy, and Vz for the same node in 3-D staggered-grid finite difference models by calculating the average value of Vy, and Vz of the nodes around the same node. We ascertain the wave patterns and their propagation characteristics in both symmetrical and asymmetric coal mine roadway models. Our simulation results indicate that the Rayleigh channel wave is stronger than the Love channel wave in front of the roadway face. The reflected Rayleigh waves from the roadway face are concentrated in the coal seam, release less energy to the roof and floor, and propagate for a longer distance. There are surface waves and refraction head waves around the roadway. In the seismic records, the Rayleigh wave energy is stronger than that of the Love channel wave along coal walls of the roadway, and the interference of the head waves and surface waves with the Rayleigh channel wave is weaker than with the Love channel wave. It is thus difficult to identify the Love channel wave in the seismic records. Increasing the depth of the receivers in the coal walls can effectively weaken the interference of surface waves with the Rayleigh channel wave, but cannot weaken the interference of surface waves with the Love channel wave. Our research results also suggest that the Love channel wave, which is often used to detect geological structures in coal mine stopes, is not suitable for detecting geological structures in front of coal mine roadways. Instead, the Rayleigh channel wave can be used for the advance detection of geological structures in coal mine roadways.
基金supported by the National Natural Science Foundation of China(Grants No.51239001,51179015,and 51509023)the Open Research Foundation of the Key Laboratory of the Pearl River Estuarine Dynamics and Associated Process Regulation,the Ministry of Water Resources(Grant No.2018KJ03)+1 种基金the Key Laboratory of Water-Sediment Sciences and Water Disaster Prevention of Hunan Province(Grant No.2017SS04)the Key Laboratory of Technology for Safeguarding of Maritime Rights and Interests and Application,State Oceanic Administration(Grant No.SCS1606)
文摘The simulating waves nearshore(SWAN) model has typically been designed for wave simulations in near-shore regions. In this study, the model's applicability to the simulation of typhoon waves in the South China Sea(SCS) was evaluated. A blended wind field, consisting of an interior domain based on Fujita's model and an exterior domain based on Takahashi's model, was used as the driving wind field. The waves driven by Typhoon Kai-tak over the SCS that occurred in 2012 were selected for the numerical simulation research. Sensitivity analyses of time step, grid resolution, and angle resolution were performed in order to obtain optimal model settings. Through sensitivity analyses, it can be found that the time step has a large influence on the results, while grid resolution and angle resolution have a little effect on the results.
基金supported by the National Natural Science Foundation of China(Grant Nos.51109187,51239002 and 51221961the Project form Zhoushan Science and Technology Bureau(Grant No.2013C41002)
文摘A numerical model of flatfish cage is built based on the lumped mass method and the principle of rigid body kinematics. To validate the numerical model, a series of physical model tests are conducted in the wave flume. The numerical results correspond well with the data sets from physical model test. The effect of weight of bottom frame, height of fish net and net shape on motion responses of fish cage and tension force on mooring lines is then analyzed. The results indicate that the vertical displacements of float collar and bottom frame decrease with the increase in the weight of bottom frame; the maximum tension force on mooring lines increases with the increasing weight of bottom frame. The inclination angles of float collar and bottom frame decrease with the increasing net height; the maximum tension force increases obviously with the increase of net height.
