This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equ...This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equations. Numerical method for modelling the short waves is investigated in detail. The forces, such as Coriolis forces, wind stress, atmosphere and bottom friction, are considered. A two-dimensional implicit difference scheme of Boussinesq equations is proposed. The low-reflection outflow open boundary is suggested. By means of this model,both velocity fields of circulation current in a channel with step expansion and the wave diffraction behind a semi-infinite breakwater are computed, and the results are satisfactory.展开更多
The existing numerical models for nearshore waves are briefly introduced, and the third-generation numerical model for shallow water wave, which makes use of the most advanced productions of wave research and has been...The existing numerical models for nearshore waves are briefly introduced, and the third-generation numerical model for shallow water wave, which makes use of the most advanced productions of wave research and has been adapted well to be used in the environment of seacoast, lake and estuary area, is particularly discussed. The applied model realizes the significant wave height distribution at different wind directions. To integrate the model into the coastal area sediment, sudden deposition mechanism, the distribution of average silt content and the change of sediment sudden deposition thickness over time in the nearshore area are simulated. The academic productions can give some theoretical guidance to the applications of sediment sudden deposition mechanism for stormy waves in the coastal area. And the advancing directions of sediment sudden deposition model are prospected.展开更多
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.展开更多
In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced. The method is simple, convenient and eff ective. And the method isn 't only ...In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced. The method is simple, convenient and eff ective. And the method isn 't only limited to conservation laws unlike other methods and maybe easily extended to multi-dimension problems. Furthermore, the numerical dissipation of the method can be flexibly regulated, so that it is especially suitable for solving various discontinuity problems.The paper shows us now to use this approach to simulate deformation and breaking of a nonlinear shallow water wave on a gentle slope, and to compute two-dimensional dam failure problem.展开更多
The safety accidents caused by collapse column water diversion occur frequently, which has great hidden danger to the safety production of coal mine. Limited by the space of underground, the detection of collapse colu...The safety accidents caused by collapse column water diversion occur frequently, which has great hidden danger to the safety production of coal mine. Limited by the space of underground, the detection of collapse column on the outside of working face has been a difficult problem. Based on this, numerical simulation and imaging research were carried out in this paper. The results indicate that when a seismic source near the roadway is excited, a part of seismic wave propagates along the roadway direction, namely direct P-wave, direct S-wave and direct Love channel wave.<span style="font-family:;" "=""> </span><span style="font-family:Verdana;">When the body waves and Love channel wave propagating to the outside</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">of working face meet the interface of collapse column, the reflected Love channel wave and reflected body waves are generated.</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">Reflection </span><span style="font-family:Verdana;">body </span><span style="font-family:Verdana;">waves and direct waves are mixed in time domain, which is difficult to identify in seismic records, while reflected Love channel wave whose amplitude is relatively strong. The </span><span style="font-family:Verdana;">reflected </span><span style="font-family:Verdana;">Love channel wave which has a large interval from other wave trains in the time domain is easily recognizable in seismic record,</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">which</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">makes it suitable for advanced detection of collapse column. The signal-to-noise ratio of X component is higher than that of Y component and Z component. According to the seismic records, polarization filtering was carried out to enhance the effective wave, which removed the interference waves, and the signal was migrated to get the position parameters of collapse column interface, which was basically consistent with the model position.</span>展开更多
The Changjiang River diluted water(CDW)spreads into the East China Sea(ECS)primarily in a plume pattern,although in some years,low-salinity water lenses(LSWLs)detach from the main body of the CDW.In-situ observations ...The Changjiang River diluted water(CDW)spreads into the East China Sea(ECS)primarily in a plume pattern,although in some years,low-salinity water lenses(LSWLs)detach from the main body of the CDW.In-situ observations indicate that in August 2006,a LSWL detached from the main body of the CDW near the river mouth.In this paper,the effects of winds,tides,baroclinity and upwelling on LSWLs are explored with a threedimensional model.The results show that:(1)winds play a crucial role in these detachment events because windinduced northerly Eulerian residual currents impose an uneven force on the CDW and cut it off,thus forming a LSWL;(2)upwelling carries high-salinity water from the lower layer to the upper layer,truncating the low-salinity water tongue vertically,which is conducive to the detachment and maintenance of LSWLs;and(3)upwelling during the evolution of a LSWL is caused by the combined effects of winds and tides.The influences of windinduced upwelling are mainly near the shore,whereas the upwelling along the 30 m isobath is predominantly affected by tides,with the effect increasing from neap tide to spring tide.展开更多
When separated water droplets condense on the surface of a composite insulator,the electrical field on the insulator surface is distorted.In turn,such distortions change the trajectories of pollution particles.In this...When separated water droplets condense on the surface of a composite insulator,the electrical field on the insulator surface is distorted.In turn,such distortions change the trajectories of pollution particles.In this study,the COMSOL software is used to simulate such a process for the FXBW4-10/100 composite insulator with or without water droplets condensation under a 10 kV DC voltage.The influence of the wind speed and particles concentration on the contamination characteristics of the considered 110 kV insulator is analyzed.The results show that:1)in the presence of water droplets on the insulator surface,the ratio of electrical field force and gravity acting on the particles is large;2)the contamination on the insulator surface increases with the wind speed;3)when the wind speed is small,the relationship between the contamination amount and the pollution concentration is essentially linear.展开更多
This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are result...This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are resulted boundary and initial conditions. The method of splitting into physical processes receives system from three equations. Then we define the approximation order and investigate stability conditions of the discrete model. The sweep method was used to calculate the system of equations. This work presents surface gravity wave profiles for different propagation phases.展开更多
A reverberation model for estimating the average reverberation intensity in layered shallow water is presented.The reverberation intensity is calculated in terms of ray theory for short range and normal mode theory fo...A reverberation model for estimating the average reverberation intensity in layered shallow water is presented.The reverberation intensity is calculated in terms of ray theory for short range and normal mode theory for long range. The calculation accuracy has been improved by taking into account the effect of complex eigenvalues on the incident normal mode field. From the comparison between different scattering models it has been shown that the separable bistatic-backscattering model is acceptable. This makes it possible to calculate reverberation by using only the monostatic-backscattering coefficient and to save greatly the computing 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.展开更多
Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes...Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.展开更多
VOF (volume of fluid) method has been used to make the numerical simulation of freak wave come true. The comparisons between the numerical results and linear theoretical results corresponding to Eq.(5) have been c...VOF (volume of fluid) method has been used to make the numerical simulation of freak wave come true. The comparisons between the numerical results and linear theoretical results corresponding to Eq.(5) have been carried out to show that the numerical results have a better exhibition of nonlinear characteristics. Wavelet analysis method has been adopted to investigate the time-frequency energy spectrum of simulation freak waves and the results reveal strong nonlinear interaction enables energy to be transferred to high harmonics during the progress of its formation. Varying water depth can enhance the nonlinear interaction, making much more energy be transferred to high harmonics and freak waves with higher asymmetry be generated.展开更多
The acoustic response characteristics of shales were investigated by the acoustic transmission experiment,which is the basis of solving geological and engineering problems using the seismic or logging information duri...The acoustic response characteristics of shales were investigated by the acoustic transmission experiment,which is the basis of solving geological and engineering problems using the seismic or logging information during the process of the exploration and development of shale gas reservoirs.Based on the theory of acoustic wave and the background of acoustic transmission experiment,the initial condition,vibration source condition,boundary condition and stability condition were constructed,and the numerical simulation of acoustic transmission experiment of shales were completed through Matlab programming.The results show that under the same bedding angle,the acoustic time and attenuation coefficient of shales shown positive correlation with the bedding density;whereas under the same bedding density,the variation laws of the acoustic time and the attenuation coefficient of shales were more complex with the change of the bedding angle,that is,the acoustic time and attenuation coefficient of shales increased first,then decreased and then increased again with the increase of the bedding angle.展开更多
Wave and longshore current interaction was examined based on the numerical models.In these models,water waves in the presence of longshore currents were modeled by parabolic mild slope equation,and wave breaking induc...Wave and longshore current interaction was examined based on the numerical models.In these models,water waves in the presence of longshore currents were modeled by parabolic mild slope equation,and wave breaking induced longshore currents were modeled by shallow water equation.Water wave provided the radiation stress gradients to drive current.Wave and longshore current interactions were considered by cycling the wave and longshore current models to a steady state.The experiments for regular and irregular breaking wave induced longshore currents by Hamilton and Ebersole (2001) and Reniers and Battjes (1997) were simulated.The numerical results indicate that the present models are effective for simulating the interaction of wave and breaking wave induced longshore currents,and the numerically simulated longshore current at wave breaking point considering wave and longshore current interaction show some disagreement with those neglecting the wave-current interaction,and the breaking wave induced longshore current effect on wave transformation is not obvious.展开更多
To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a "predicted water depth" is evaluated e...To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a "predicted water depth" is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in beth typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.展开更多
The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the...The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation, A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.展开更多
A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid (CLEAR-VOF) method was developed to simulate wave and flow problems. The N...A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid (CLEAR-VOF) method was developed to simulate wave and flow problems. The Navier-Stokes equations were discretized with a three-step finite element method that has a third-order accuracy. In the CLEAR-VOF method, the VOF function F was calculated in the Lagrangian manner and allowed the complicated free surface to be accurately captured. The propagation of regular waves and solitary waves over a flat bottom, and shoaling and breaking of solitary waves on two different slopes were simulated with this model, and the numerical results agreed with experimental data and theoretical solutions. A benchmark test of dam-collapse flow was also simulated with an unstructured mesh, and the capability of the present model for wave and flow simulations with unstructured meshes, was verified. The results show that the model is effective for numerical simulation of wave and flow problems with both structured and unstructured meshes.展开更多
Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directl...Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the effciency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.展开更多
Impulse turbine, working as a typical self-rectifying turbine, is recently utilized for the oscillating water column(OWC) wave energy converters, which can rotate in the same direction under the bi-directional air f...Impulse turbine, working as a typical self-rectifying turbine, is recently utilized for the oscillating water column(OWC) wave energy converters, which can rotate in the same direction under the bi-directional air flows. A numerical model established in Fluent is validated by the corresponding experimental results. The flow fields, pressure distribution and dimensionless evaluating coefficients can be calculated and analyzed. Effects of the rotor solidity varying with the change of blade number are investigated and the suitable solidity value is recommended for different flow coefficients.展开更多
Long waves generated by a moving atmospheric pressure distribution, associated with a storm, in coastal region are investigated numerically. For simplicity the moving atmospheric pressure is assumed to be moving only ...Long waves generated by a moving atmospheric pressure distribution, associated with a storm, in coastal region are investigated numerically. For simplicity the moving atmospheric pressure is assumed to be moving only in the alongshore direction and the beach slope is assumed to be a constant in the on-offshore direction. By solving the linear shallow water equations we obtain numerical solutions for a wide range of physical parameters, including storm size (2a), storm speed (U), and beach slope (a). Based on the numerical results, it is determined that edge wave packets are generated if the storm speed is equal to or greater than the critical velocity, Ucr, which is defined as the phase speed of the fundamental edge wave mode whose wavelength is scaled by the width of the storm size. The length and the location of the positively moving edge wave packet is roughly Ut/2 〈 y 〈 Ut, where y is in the alongshore direction and t is the time. Once the edge wave packet is generated, the wavelength is the same as that of the fundamental edge wave mode corresponding to the storm speed and is independent of the storm size, which can, however, affect the wave amplitude. When the storm speed is less than the critical velocity, the primary surface signature is a depression directly correlated to the atmospheric pressure distribution.展开更多
基金Supported by the Fund of National Nature Sciences of China
文摘This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equations. Numerical method for modelling the short waves is investigated in detail. The forces, such as Coriolis forces, wind stress, atmosphere and bottom friction, are considered. A two-dimensional implicit difference scheme of Boussinesq equations is proposed. The low-reflection outflow open boundary is suggested. By means of this model,both velocity fields of circulation current in a channel with step expansion and the wave diffraction behind a semi-infinite breakwater are computed, and the results are satisfactory.
文摘The existing numerical models for nearshore waves are briefly introduced, and the third-generation numerical model for shallow water wave, which makes use of the most advanced productions of wave research and has been adapted well to be used in the environment of seacoast, lake and estuary area, is particularly discussed. The applied model realizes the significant wave height distribution at different wind directions. To integrate the model into the coastal area sediment, sudden deposition mechanism, the distribution of average silt content and the change of sediment sudden deposition thickness over time in the nearshore area are simulated. The academic productions can give some theoretical guidance to the applications of sediment sudden deposition mechanism for stormy waves in the coastal area. And the advancing directions of sediment sudden deposition model are prospected.
基金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.
文摘In this paper a new approach for designing upwind type schemes-the characterizing-integral method and its applied skills are introduced. The method is simple, convenient and eff ective. And the method isn 't only limited to conservation laws unlike other methods and maybe easily extended to multi-dimension problems. Furthermore, the numerical dissipation of the method can be flexibly regulated, so that it is especially suitable for solving various discontinuity problems.The paper shows us now to use this approach to simulate deformation and breaking of a nonlinear shallow water wave on a gentle slope, and to compute two-dimensional dam failure problem.
文摘The safety accidents caused by collapse column water diversion occur frequently, which has great hidden danger to the safety production of coal mine. Limited by the space of underground, the detection of collapse column on the outside of working face has been a difficult problem. Based on this, numerical simulation and imaging research were carried out in this paper. The results indicate that when a seismic source near the roadway is excited, a part of seismic wave propagates along the roadway direction, namely direct P-wave, direct S-wave and direct Love channel wave.<span style="font-family:;" "=""> </span><span style="font-family:Verdana;">When the body waves and Love channel wave propagating to the outside</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">of working face meet the interface of collapse column, the reflected Love channel wave and reflected body waves are generated.</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">Reflection </span><span style="font-family:Verdana;">body </span><span style="font-family:Verdana;">waves and direct waves are mixed in time domain, which is difficult to identify in seismic records, while reflected Love channel wave whose amplitude is relatively strong. The </span><span style="font-family:Verdana;">reflected </span><span style="font-family:Verdana;">Love channel wave which has a large interval from other wave trains in the time domain is easily recognizable in seismic record,</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">which</span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">makes it suitable for advanced detection of collapse column. The signal-to-noise ratio of X component is higher than that of Y component and Z component. According to the seismic records, polarization filtering was carried out to enhance the effective wave, which removed the interference waves, and the signal was migrated to get the position parameters of collapse column interface, which was basically consistent with the model position.</span>
基金The National Natural Science Foundation of China under contract No.41376012.
文摘The Changjiang River diluted water(CDW)spreads into the East China Sea(ECS)primarily in a plume pattern,although in some years,low-salinity water lenses(LSWLs)detach from the main body of the CDW.In-situ observations indicate that in August 2006,a LSWL detached from the main body of the CDW near the river mouth.In this paper,the effects of winds,tides,baroclinity and upwelling on LSWLs are explored with a threedimensional model.The results show that:(1)winds play a crucial role in these detachment events because windinduced northerly Eulerian residual currents impose an uneven force on the CDW and cut it off,thus forming a LSWL;(2)upwelling carries high-salinity water from the lower layer to the upper layer,truncating the low-salinity water tongue vertically,which is conducive to the detachment and maintenance of LSWLs;and(3)upwelling during the evolution of a LSWL is caused by the combined effects of winds and tides.The influences of windinduced upwelling are mainly near the shore,whereas the upwelling along the 30 m isobath is predominantly affected by tides,with the effect increasing from neap tide to spring tide.
基金Science and Technology Project of State Grid Corporation(GY7111053)[www.sgcc.com.cn]。
文摘When separated water droplets condense on the surface of a composite insulator,the electrical field on the insulator surface is distorted.In turn,such distortions change the trajectories of pollution particles.In this study,the COMSOL software is used to simulate such a process for the FXBW4-10/100 composite insulator with or without water droplets condensation under a 10 kV DC voltage.The influence of the wind speed and particles concentration on the contamination characteristics of the considered 110 kV insulator is analyzed.The results show that:1)in the presence of water droplets on the insulator surface,the ratio of electrical field force and gravity acting on the particles is large;2)the contamination on the insulator surface increases with the wind speed;3)when the wind speed is small,the relationship between the contamination amount and the pollution concentration is essentially linear.
文摘This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are resulted boundary and initial conditions. The method of splitting into physical processes receives system from three equations. Then we define the approximation order and investigate stability conditions of the discrete model. The sweep method was used to calculate the system of equations. This work presents surface gravity wave profiles for different propagation phases.
文摘A reverberation model for estimating the average reverberation intensity in layered shallow water is presented.The reverberation intensity is calculated in terms of ray theory for short range and normal mode theory for long range. The calculation accuracy has been improved by taking into account the effect of complex eigenvalues on the incident normal mode field. From the comparison between different scattering models it has been shown that the separable bistatic-backscattering model is acceptable. This makes it possible to calculate reverberation by using only the monostatic-backscattering coefficient and to save greatly the computing 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.
文摘Dam-break flows pose significant threats to urban areas due to their potential for causing rapid and extensive flooding. Traditional numerical methods for simulating these events struggle with complex urban landscapes. This paper presents an alternative approach using Radial Basis Functions to simulate dam-break flows and their impact on urban flood inundation. The proposed method adapts a new strategy based on Particle Swarm Optimization for variable shape parameter selection on meshfree formulation to enhance the numerical stability and convergence of the simulation. The method’s accuracy and efficiency are demonstrated through numerical experiments, including well-known partial and circular dam-break problems and an idealized city with a single building, highlighting its potential as a valuable tool for urban flood risk management.
文摘VOF (volume of fluid) method has been used to make the numerical simulation of freak wave come true. The comparisons between the numerical results and linear theoretical results corresponding to Eq.(5) have been carried out to show that the numerical results have a better exhibition of nonlinear characteristics. Wavelet analysis method has been adopted to investigate the time-frequency energy spectrum of simulation freak waves and the results reveal strong nonlinear interaction enables energy to be transferred to high harmonics during the progress of its formation. Varying water depth can enhance the nonlinear interaction, making much more energy be transferred to high harmonics and freak waves with higher asymmetry be generated.
基金This research is supported by the National Science and Technology Major Project(Grant No.2019A-3307)the National Natural Science Foundation of China(Grant No.41872167).
文摘The acoustic response characteristics of shales were investigated by the acoustic transmission experiment,which is the basis of solving geological and engineering problems using the seismic or logging information during the process of the exploration and development of shale gas reservoirs.Based on the theory of acoustic wave and the background of acoustic transmission experiment,the initial condition,vibration source condition,boundary condition and stability condition were constructed,and the numerical simulation of acoustic transmission experiment of shales were completed through Matlab programming.The results show that under the same bedding angle,the acoustic time and attenuation coefficient of shales shown positive correlation with the bedding density;whereas under the same bedding density,the variation laws of the acoustic time and the attenuation coefficient of shales were more complex with the change of the bedding angle,that is,the acoustic time and attenuation coefficient of shales increased first,then decreased and then increased again with the increase of the bedding angle.
基金The National Natural Science Foundation of China under contract Nos 50839001,51179025 and 50709004the Specialized Research Fund for the Doctoral Program of Higher Education of China under contract No.20070141032
文摘Wave and longshore current interaction was examined based on the numerical models.In these models,water waves in the presence of longshore currents were modeled by parabolic mild slope equation,and wave breaking induced longshore currents were modeled by shallow water equation.Water wave provided the radiation stress gradients to drive current.Wave and longshore current interactions were considered by cycling the wave and longshore current models to a steady state.The experiments for regular and irregular breaking wave induced longshore currents by Hamilton and Ebersole (2001) and Reniers and Battjes (1997) were simulated.The numerical results indicate that the present models are effective for simulating the interaction of wave and breaking wave induced longshore currents,and the numerically simulated longshore current at wave breaking point considering wave and longshore current interaction show some disagreement with those neglecting the wave-current interaction,and the breaking wave induced longshore current effect on wave transformation is not obvious.
基金supported by the National Natural Science Foundation of China (Grant No.10872144)the Support Plan of Science and Technology of Tianjin (Grant No.07ZCGYSH01700)
文摘To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying (WD) method is proposed. In the new method, a "predicted water depth" is evaluated explicitly based on the simplified shallow water equations and used to determine the status (wet or dry) together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in beth typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.
基金This subject was financially supported by the National Natural Science Foundation of China(Grant No. 59839330 and No.59976047) by the Visiting Scholal Foundation of State Key Hydraulic Lab.of High Speed Flows of Dalian University of Technology.
文摘The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation, A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.
基金supported by the National Natural Science Foundation of China (Grant No. 50679008)
文摘A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid (CLEAR-VOF) method was developed to simulate wave and flow problems. The Navier-Stokes equations were discretized with a three-step finite element method that has a third-order accuracy. In the CLEAR-VOF method, the VOF function F was calculated in the Lagrangian manner and allowed the complicated free surface to be accurately captured. The propagation of regular waves and solitary waves over a flat bottom, and shoaling and breaking of solitary waves on two different slopes were simulated with this model, and the numerical results agreed with experimental data and theoretical solutions. A benchmark test of dam-collapse flow was also simulated with an unstructured mesh, and the capability of the present model for wave and flow simulations with unstructured meshes, was verified. The results show that the model is effective for numerical simulation of wave and flow problems with both structured and unstructured meshes.
基金supported by National Natural Science Foundation of China (NSFC) projects (Grant Nos. 40875065 and 40805045)the research projects 2008R001 at Chinese Academy of Meteorological Sciences (CAMS) and 2008 LASWZI05 at the State Key Laboratory of Severe Weather, CAMS
文摘Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the effciency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51279190 and 51311140259)the Shandong Province Natural Science Foundation for Distinguished Young Scholars(Grant No.JQ201314)
文摘Impulse turbine, working as a typical self-rectifying turbine, is recently utilized for the oscillating water column(OWC) wave energy converters, which can rotate in the same direction under the bi-directional air flows. A numerical model established in Fluent is validated by the corresponding experimental results. The flow fields, pressure distribution and dimensionless evaluating coefficients can be calculated and analyzed. Effects of the rotor solidity varying with the change of blade number are investigated and the suitable solidity value is recommended for different flow coefficients.
基金supported by an NSF grant to Cornell University,the China Scholarship Council and a Korean government MLTMA grant Development of Korea Operational Oceanographic System (KOOS) to KORDI
文摘Long waves generated by a moving atmospheric pressure distribution, associated with a storm, in coastal region are investigated numerically. For simplicity the moving atmospheric pressure is assumed to be moving only in the alongshore direction and the beach slope is assumed to be a constant in the on-offshore direction. By solving the linear shallow water equations we obtain numerical solutions for a wide range of physical parameters, including storm size (2a), storm speed (U), and beach slope (a). Based on the numerical results, it is determined that edge wave packets are generated if the storm speed is equal to or greater than the critical velocity, Ucr, which is defined as the phase speed of the fundamental edge wave mode whose wavelength is scaled by the width of the storm size. The length and the location of the positively moving edge wave packet is roughly Ut/2 〈 y 〈 Ut, where y is in the alongshore direction and t is the time. Once the edge wave packet is generated, the wavelength is the same as that of the fundamental edge wave mode corresponding to the storm speed and is independent of the storm size, which can, however, affect the wave amplitude. When the storm speed is less than the critical velocity, the primary surface signature is a depression directly correlated to the atmospheric pressure distribution.