Fifteen formulae of wave transmission coefficient for submerged breakwaters obtained during last 3 decades are presented, compared, and analyzed in this paper. The dimensionless parameters mainly involved in this disc...Fifteen formulae of wave transmission coefficient for submerged breakwaters obtained during last 3 decades are presented, compared, and analyzed in this paper. The dimensionless parameters mainly involved in this discussion are the relative submerged depth Re/h, relative wave height Rc/Hi, relative rubble size B/D50, relative breakwater width B√ HiL0 and wave breaker index ξ. It indicates that there exist notable differences among the computed results, which mainly originate from the limited experimental conditions and different analytical methods, even though the major tendency keeps similar. It is necessary to conduct more systematic studies to obtain better understanding about the mechanism of wave transmission over submerged breakwaters.展开更多
A simple method for solving Cauchy’s problem of wave equations in higher space dimensions with initial condition of separated variables, has been given by using D’Alembert’s formula and some examples have been shown.
Laboratory experiments are performed to explore the response rule of a sandy beach profile under plunging wave on a non-uniform sediment-bed slope. The initial beach slope of combination of 1/10 and 1/20 is exposed to...Laboratory experiments are performed to explore the response rule of a sandy beach profile under plunging wave on a non-uniform sediment-bed slope. The initial beach slope of combination of 1/10 and 1/20 is exposed to regular waves and cnoidal waves respectively. The free surface elevation, process of wave propagation, wave breaking, uprush and backwash and the change of a cross-shore beach profile are measured and recorded. The beach profile under the regular waves action exhibits two parts： a sandbar profile and a beach berm profile, and only one typical profile transformation under the cnoidal waves action is obtained, which is the beach berm profile. In the laboratory experiments, it is found that the beach states under wave action related to the previous factors. In addition, they are related to the characteristic of breaking waves such as the breaking intensity of the plunging wave. A concept about the characteristic angle of the plunging wave has been put forward through the observation and analysis of the phenomenon of the laboratory experiment. A qualitative analysis about the sediment transport carrying by currents generated from the plunging wave and the state of beach profile under the wave action has been done. The quantitative analysis about the relationship between the characteristic angle and Irribarren number has been done. An available formula of equilibrium states for the sandy beach induced by the plunging wave has been established based on the relationship between Irribarren number and the beach profile. By fitting these experimental results and others＇ experimental results to three lines, the three fitting coefficients can be calculated in their formula respectively. The recommended empirical formulas can divide three states of a beach morphology profile obviously, which include a depositive beach, an erosive beach and an intermediate beach.展开更多
Air–sea exchange plays a vital role in the development and maintenance of tropical cyclones（TCs）. Although studies have suggested the dependence of air–sea fluxes on surface waves and sea spray, how these processe...Air–sea exchange plays a vital role in the development and maintenance of tropical cyclones（TCs）. Although studies have suggested the dependence of air–sea fluxes on surface waves and sea spray, how these processes modify those fluxes under TC conditions have not been sufficiently investigated based on in-situ observations.Using continuous meteorological and surface wave data from a moored buoy in the northern South China Sea,this study examines the effects of surface waves and sea spray on air–sea fluxes during the passage of Typhoon Hagupit. The mooring was within about 40 km of the center of Hagupit. Surface waves could increase momentum flux to the ocean by about 15%, and sea spray enhanced both sensible and latent heat fluxes to the atmosphere,causing Hagupit to absorb 500 W/m^2 more heat flux from the ocean. These results have powerful implications for understanding TC–ocean interaction and improving TC intensity forecasting.展开更多
Studies of the breaking criteria for solitary waves on a slope are presented in this paper. The boundary element method is used to model the processes of shoaling and breaking of solitary waves on various slopes. Empi...Studies of the breaking criteria for solitary waves on a slope are presented in this paper. The boundary element method is used to model the processes of shoaling and breaking of solitary waves on various slopes. Empirical formulae that can be used to characterize the breaking of solitary waves are presented. These include the breaking index, the wave height, the water depth, and the maximum particle velocity at the point of breaking. Comparisons with the results of other researches are given.展开更多
This papaer presents an analytic formula for the locus of shock front of sphericalblast waves with back pressure. Meanwhile, it is shown that the wave velocity is inversely proportions to the area of shock wave in the...This papaer presents an analytic formula for the locus of shock front of sphericalblast waves with back pressure. Meanwhile, it is shown that the wave velocity is inversely proportions to the area of shock wave in the Lagrangian coordinate space. If shock wave ispropagating into the gas at rest, the wave velocity is also inversely propagating into the gasat rest, the wave velocity is also inversely proportional to the area of shock front in the Eulerian coordinate space. The results obtained are generally applicable and may be suited to thewave velocity of contracting spherical shock wave.展开更多
In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models...In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.展开更多
In this study, systematic physical model tests were performed to investigate the wave forces on the twin-plate breakwater under irregular waves. Based on the experimental results, the effects of the relative plate wid...In this study, systematic physical model tests were performed to investigate the wave forces on the twin-plate breakwater under irregular waves. Based on the experimental results, the effects of the relative plate width B/L,wave height Hs/D and incident angle θ0 on the wave forces were analyzed and discussed. The results showed that:(1) The envelopes of the total wave pressure were generally symmetrical along the direction of plate width under the incident angles(θ0) being 0°, 15°, 30°, 45° and 60°. In particular, the envelopes of wave pressure atθ0=30° were larger than all other cases.(2) The synchronous pressure distribution of the breakwater under oblique wave action was more complicated comparing to the normal incident waves.(3) Based on data analysis, an empirical formula was obtained to estimate the total vertical force of the twin-plate breakwater.This empirical formula can be a good reference for the design basis of engineering applications under specified wave conditions.展开更多
The run-up on offshore structures induced by the steep regular wave is a highly nonlinear flow with a free surface. This article focuses on the investigation of the steep regular wave run-up on a single vertical cylin...The run-up on offshore structures induced by the steep regular wave is a highly nonlinear flow with a free surface. This article focuses on the investigation of the steep regular wave run-up on a single vertical cylinder by solving the Navier-Stokes equations. A numerical wave tank is established based on the open-source package to simulate the wave scattering induced by a vertical cylinder. The VOF method is applied to capture the large deformation and breaking of the free surface. The numerical model is validated by experimental results. The relative wave run-ups on the front face and the back face along the centerline of a cylinder are analyzed. The changes of the relative run-ups with the wave steepness, the relative diameter and the relative depth are studied. It is found that the relative run-ups on the front face and the back face of the cylinder depend mainly on the wave steepness and the relative diameter, while the dependence on the relative depth is weak. The empirical formulae are proposed to calculate the relative run-ups in terms of the wave steepness of incident regular waves and the relative diameter of a cylinder.展开更多
Based on the 2nd order cnoidal wave theory, the characters of shallow water standing waves and their action on vertical walls are studied in this paper. The theoretical expressions of the wave surface elevation in fro...Based on the 2nd order cnoidal wave theory, the characters of shallow water standing waves and their action on vertical walls are studied in this paper. The theoretical expressions of the wave surface elevation in front of and the wave pressure on the vertical wall are obtained. In order to verify the theoretical results, model tests were made in the State Key Laboratory of Coastal and Offshore Engineering at DUT. For the wave surface elevation in front of the wall and the wave forces on the wall at the moment when the wave surface at the wall surface goes down to the bottom of the wave trough, the calculated results coincide quite well with the experimental results. For the wave forces on the wall at the moment when the wave surface at the wall surface goes up to the top of the wave crest, the theoretical expressions are modified by the experimental results. For the convenience of practical use, calculations are made for the wave conditions which usually occur in enginering practice by use of the investigated results obtained in this paper. Empirical formulas are fitted with these calculated results for designers to use.展开更多
The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The co...The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.展开更多
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.展开更多
Accurate estimation of wave uplift force is essential to the designs of reliable coastal and marine structures.We presents a series of laboratory work here on the impact of regular waves on horizontal panels,from whic...Accurate estimation of wave uplift force is essential to the designs of reliable coastal and marine structures.We presents a series of laboratory work here on the impact of regular waves on horizontal panels,from which an empirical formula to estimate accurately the wave uplift force on panels is established.The laboratory measurements show that the wave uplift force depends mainly on the incident wave height,the wave period,the wave length,the panel width,and the clearance between the subsurface of the panel and the still water level.Among these factors,the impact of the panel width on uplift forces is relatively complicated.Result shows that the relative panel width(i.e.,the ratio of panel width to wave length)plays a more important role in estimating the wave uplift force.Based on our comprehensive laboratory measurements,we further developed an empirical formula to compute wave uplift force on horizontal panels through dimensionless analysis.Compared with other empirical formulas,this formula uses dimensionless variables of clear physical meanings,thus can describe the interaction between waves and the panels in a better way.In addition,the efficiency of the formula to estimate wave uplift force on horizontal panels is verified against existing works.Therefore,the findings in this study shall be useful for understanding the mechanism of wave uplift force on horizontal panels and numerical model validation.展开更多
In work questions of distribution of waves in a viscoelastic wedge with any corner of top is considered. The elastic cylinder with a radial crack is a wedge corner. The regional task for system of the differential equ...In work questions of distribution of waves in a viscoelastic wedge with any corner of top is considered. The elastic cylinder with a radial crack is a wedge corner. The regional task for system of the differential equations in private derivatives is decided by means of a method of straight lines that allows using a method of orthogonal prorace.展开更多
We explore the theoretical possibility of extending the usual squeezed state to those produced by nonlinear singlemode squeezing operators. We derive the wave functions of exp[-（ig/2）（（1-X2）1/2P ＋ P（1-X2）1/2）...We explore the theoretical possibility of extending the usual squeezed state to those produced by nonlinear singlemode squeezing operators. We derive the wave functions of exp[-（ig/2）（（1-X2）1/2P ＋ P（1-X2）1/2）]｜0 in the coordinate representation. A new operator＇s disentangling formula is derived as a by-product.展开更多
内孤立波广泛存在于实际海洋中,由于其诱导的水平流速大,对海洋结构物水动力性能的影响不可忽视.近年来,针对内孤立波对浮式海洋结构物作用的载荷理论逐步完善,但对于内孤立波与潜体的作用特性有待进一步研究.文中将Morison公式和Froude...内孤立波广泛存在于实际海洋中,由于其诱导的水平流速大,对海洋结构物水动力性能的影响不可忽视.近年来,针对内孤立波对浮式海洋结构物作用的载荷理论逐步完善,但对于内孤立波与潜体的作用特性有待进一步研究.文中将Morison公式和Froude-Krylov公式相结合,建立了内孤立波作用下水下潜体载荷的一种简化理论模型.基于流体体积(volume of fluid,VOF)方法,结合3类内孤立波理论模型,设计速度入口函数,开展内孤立波与水下潜体相互作用数值模拟研究,结合文献中的模型实验结果,对数值模拟方法的可靠性进行了验证.通过系列CFD模拟,结合水下潜体内孤立波载荷简化理论模型,获得了内孤立波作用下潜体模型拖曳力及惯性力系数的拟合公式.该简化理论模型对研究内孤立波作用下水下潜体的载荷特性具有一定的参考价值.展开更多
文摘Fifteen formulae of wave transmission coefficient for submerged breakwaters obtained during last 3 decades are presented, compared, and analyzed in this paper. The dimensionless parameters mainly involved in this discussion are the relative submerged depth Re/h, relative wave height Rc/Hi, relative rubble size B/D50, relative breakwater width B√ HiL0 and wave breaker index ξ. It indicates that there exist notable differences among the computed results, which mainly originate from the limited experimental conditions and different analytical methods, even though the major tendency keeps similar. It is necessary to conduct more systematic studies to obtain better understanding about the mechanism of wave transmission over submerged breakwaters.
基金Supported by the Natural Science Foundation of Hubei Province!(992P0 30 7) the National Natural Science Foun-dation of Chi
文摘A simple method for solving Cauchy’s problem of wave equations in higher space dimensions with initial condition of separated variables, has been given by using D’Alembert’s formula and some examples have been shown.
基金The National Natural Science Foundation of China under contract Nos 51239001,51179015,51409022 and 51509023the Hunan Provincial Innovation Foundation for Postgraduate under contract No.CX2015B348
文摘Laboratory experiments are performed to explore the response rule of a sandy beach profile under plunging wave on a non-uniform sediment-bed slope. The initial beach slope of combination of 1/10 and 1/20 is exposed to regular waves and cnoidal waves respectively. The free surface elevation, process of wave propagation, wave breaking, uprush and backwash and the change of a cross-shore beach profile are measured and recorded. The beach profile under the regular waves action exhibits two parts： a sandbar profile and a beach berm profile, and only one typical profile transformation under the cnoidal waves action is obtained, which is the beach berm profile. In the laboratory experiments, it is found that the beach states under wave action related to the previous factors. In addition, they are related to the characteristic of breaking waves such as the breaking intensity of the plunging wave. A concept about the characteristic angle of the plunging wave has been put forward through the observation and analysis of the phenomenon of the laboratory experiment. A qualitative analysis about the sediment transport carrying by currents generated from the plunging wave and the state of beach profile under the wave action has been done. The quantitative analysis about the relationship between the characteristic angle and Irribarren number has been done. An available formula of equilibrium states for the sandy beach induced by the plunging wave has been established based on the relationship between Irribarren number and the beach profile. By fitting these experimental results and others＇ experimental results to three lines, the three fitting coefficients can be calculated in their formula respectively. The recommended empirical formulas can divide three states of a beach morphology profile obviously, which include a depositive beach, an erosive beach and an intermediate beach.
基金Zhejiang Provincial Natural Science Foundation of China under contract No.LR15D060001the National Program on Global Change and Air-Sea Interactions under contract No.GASI-IPOVAI-04the National Natural Science Foundation of China under contract Nos 41476021,41706034 and 41321004
文摘Air–sea exchange plays a vital role in the development and maintenance of tropical cyclones（TCs）. Although studies have suggested the dependence of air–sea fluxes on surface waves and sea spray, how these processes modify those fluxes under TC conditions have not been sufficiently investigated based on in-situ observations.Using continuous meteorological and surface wave data from a moored buoy in the northern South China Sea,this study examines the effects of surface waves and sea spray on air–sea fluxes during the passage of Typhoon Hagupit. The mooring was within about 40 km of the center of Hagupit. Surface waves could increase momentum flux to the ocean by about 15%, and sea spray enhanced both sensible and latent heat fluxes to the atmosphere,causing Hagupit to absorb 500 W/m^2 more heat flux from the ocean. These results have powerful implications for understanding TC–ocean interaction and improving TC intensity forecasting.
基金The work was financially supported by the Taiwan Science Council(Project No.NSC-89-2611-E-019-064)
文摘Studies of the breaking criteria for solitary waves on a slope are presented in this paper. The boundary element method is used to model the processes of shoaling and breaking of solitary waves on various slopes. Empirical formulae that can be used to characterize the breaking of solitary waves are presented. These include the breaking index, the wave height, the water depth, and the maximum particle velocity at the point of breaking. Comparisons with the results of other researches are given.
文摘This papaer presents an analytic formula for the locus of shock front of sphericalblast waves with back pressure. Meanwhile, it is shown that the wave velocity is inversely proportions to the area of shock wave in the Lagrangian coordinate space. If shock wave ispropagating into the gas at rest, the wave velocity is also inversely propagating into the gasat rest, the wave velocity is also inversely proportional to the area of shock front in the Eulerian coordinate space. The results obtained are generally applicable and may be suited to thewave velocity of contracting spherical shock wave.
基金National Basic Research Program of China Under Grant No. 2007CB714200National Natural Science Foundation of China Under Grant No. 90715038
文摘In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order （M= 1） and the fourth order （M = 2） recursion formulas.
基金The National Natural Science Foundation of China under contract Nos 51079025 and 11272079the Research Funds from State Key Laboratory of Coastal and Offshore Engineering under contract No.LY1602
文摘In this study, systematic physical model tests were performed to investigate the wave forces on the twin-plate breakwater under irregular waves. Based on the experimental results, the effects of the relative plate width B/L,wave height Hs/D and incident angle θ0 on the wave forces were analyzed and discussed. The results showed that:(1) The envelopes of the total wave pressure were generally symmetrical along the direction of plate width under the incident angles(θ0) being 0°, 15°, 30°, 45° and 60°. In particular, the envelopes of wave pressure atθ0=30° were larger than all other cases.(2) The synchronous pressure distribution of the breakwater under oblique wave action was more complicated comparing to the normal incident waves.(3) Based on data analysis, an empirical formula was obtained to estimate the total vertical force of the twin-plate breakwater.This empirical formula can be a good reference for the design basis of engineering applications under specified wave conditions.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.11632012 and 41861144024)the National Basic Research Program of China(973 Program,Grant No.2014CB046203)
文摘The run-up on offshore structures induced by the steep regular wave is a highly nonlinear flow with a free surface. This article focuses on the investigation of the steep regular wave run-up on a single vertical cylinder by solving the Navier-Stokes equations. A numerical wave tank is established based on the open-source package to simulate the wave scattering induced by a vertical cylinder. The VOF method is applied to capture the large deformation and breaking of the free surface. The numerical model is validated by experimental results. The relative wave run-ups on the front face and the back face along the centerline of a cylinder are analyzed. The changes of the relative run-ups with the wave steepness, the relative diameter and the relative depth are studied. It is found that the relative run-ups on the front face and the back face of the cylinder depend mainly on the wave steepness and the relative diameter, while the dependence on the relative depth is weak. The empirical formulae are proposed to calculate the relative run-ups in terms of the wave steepness of incident regular waves and the relative diameter of a cylinder.
文摘Based on the 2nd order cnoidal wave theory, the characters of shallow water standing waves and their action on vertical walls are studied in this paper. The theoretical expressions of the wave surface elevation in front of and the wave pressure on the vertical wall are obtained. In order to verify the theoretical results, model tests were made in the State Key Laboratory of Coastal and Offshore Engineering at DUT. For the wave surface elevation in front of the wall and the wave forces on the wall at the moment when the wave surface at the wall surface goes down to the bottom of the wave trough, the calculated results coincide quite well with the experimental results. For the wave forces on the wall at the moment when the wave surface at the wall surface goes up to the top of the wave crest, the theoretical expressions are modified by the experimental results. For the convenience of practical use, calculations are made for the wave conditions which usually occur in enginering practice by use of the investigated results obtained in this paper. Empirical formulas are fitted with these calculated results for designers to use.
基金National Natural Science Foundation of China (50608024 and 50538050).
文摘The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.
基金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 Key R&D Program of China(No.2016YFC1402002)the National Natural Science Foundation of China(No.51579156)the Major Project of Nanjing Hydraulic Research Institute Funds(Nos.Y218005,Y218006)
文摘Accurate estimation of wave uplift force is essential to the designs of reliable coastal and marine structures.We presents a series of laboratory work here on the impact of regular waves on horizontal panels,from which an empirical formula to estimate accurately the wave uplift force on panels is established.The laboratory measurements show that the wave uplift force depends mainly on the incident wave height,the wave period,the wave length,the panel width,and the clearance between the subsurface of the panel and the still water level.Among these factors,the impact of the panel width on uplift forces is relatively complicated.Result shows that the relative panel width(i.e.,the ratio of panel width to wave length)plays a more important role in estimating the wave uplift force.Based on our comprehensive laboratory measurements,we further developed an empirical formula to compute wave uplift force on horizontal panels through dimensionless analysis.Compared with other empirical formulas,this formula uses dimensionless variables of clear physical meanings,thus can describe the interaction between waves and the panels in a better way.In addition,the efficiency of the formula to estimate wave uplift force on horizontal panels is verified against existing works.Therefore,the findings in this study shall be useful for understanding the mechanism of wave uplift force on horizontal panels and numerical model validation.
文摘In work questions of distribution of waves in a viscoelastic wedge with any corner of top is considered. The elastic cylinder with a radial crack is a wedge corner. The regional task for system of the differential equations in private derivatives is decided by means of a method of straight lines that allows using a method of orthogonal prorace.
基金Project supported by the National Natural Science Foundation of China(Grant No.11175113)
文摘We explore the theoretical possibility of extending the usual squeezed state to those produced by nonlinear singlemode squeezing operators. We derive the wave functions of exp[-（ig/2）（（1-X2）1/2P ＋ P（1-X2）1/2）]｜0 in the coordinate representation. A new operator＇s disentangling formula is derived as a by-product.
文摘内孤立波广泛存在于实际海洋中,由于其诱导的水平流速大,对海洋结构物水动力性能的影响不可忽视.近年来,针对内孤立波对浮式海洋结构物作用的载荷理论逐步完善,但对于内孤立波与潜体的作用特性有待进一步研究.文中将Morison公式和Froude-Krylov公式相结合,建立了内孤立波作用下水下潜体载荷的一种简化理论模型.基于流体体积(volume of fluid,VOF)方法,结合3类内孤立波理论模型,设计速度入口函数,开展内孤立波与水下潜体相互作用数值模拟研究,结合文献中的模型实验结果,对数值模拟方法的可靠性进行了验证.通过系列CFD模拟,结合水下潜体内孤立波载荷简化理论模型,获得了内孤立波作用下潜体模型拖曳力及惯性力系数的拟合公式.该简化理论模型对研究内孤立波作用下水下潜体的载荷特性具有一定的参考价值.
