CN113408179A - Dynamic simulation method for calculating real-time wave breaking-caused mixing - Google Patents

Dynamic simulation method for calculating real-time wave breaking-caused mixing Download PDF

Info

Publication number
CN113408179A
CN113408179A CN202110801898.2A CN202110801898A CN113408179A CN 113408179 A CN113408179 A CN 113408179A CN 202110801898 A CN202110801898 A CN 202110801898A CN 113408179 A CN113408179 A CN 113408179A
Authority
CN
China
Prior art keywords
wave
ocean
model
breaking
wave breaking
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110801898.2A
Other languages
Chinese (zh)
Inventor
邓增安
王梦涵
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202110801898.2A priority Critical patent/CN113408179A/en
Publication of CN113408179A publication Critical patent/CN113408179A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a dynamic simulation method for calculating real-time wave breaking-caused mixing, which comprises the following steps: (1) calculating a wave field in a certain time period through an offshore wave mode, and outputting wave parameters required by calculating wave breaking parameters; (2) establishing a wave breaking parameter calculation model, taking wave parameters as input of the wave breaking parameter calculation model, and calculating through the wave breaking parameter calculation model to obtain a wave breaking parameter field; (3) improving the Princeton ocean model to enable the improved Princeton ocean model to comprise an unsteady wave breaking parameterization calculation method; (4) and taking the wave breaking parameter field as the input of the improved Princeton ocean model, simulating the ocean dynamic process containing unsteady wave breaking in real time through the improved Princeton ocean model, and outputting the ocean hydrology element data. The invention realizes the dynamic simulation of the physical ocean model and overcomes the defect that the wave breaking parameters are set to be constant values in the original physical ocean model.

Description

Dynamic simulation method for calculating real-time wave breaking-caused mixing
Technical Field
The invention relates to a physical ocean, in particular to a dynamic simulation method for calculating real-time wave breaking mixing, which is applied to dynamic real-time simulation of a physical ocean model.
Background
Wave breaking at the surface of the ocean generates a downward input turbulent kinetic energy flux, with the momentum of the wind being transferred first to the waves and then to the flow field, primarily by way of wave breaking. During this process, the turbulence energy near the surface of the sea increases significantly, creating a sub-layer of significantly increased turbulence generation near the surface of the sea. The enhanced turbulent kinetic energy can transfer momentum to the lower ocean layer, thereby influencing the vertical distribution of near-surface circulation and simultaneously influencing the vertical distribution of the upper ocean layer temperature field. The wave breaking phenomenon happens in the sea at any moment, plays a key dynamic role in controlling the power process of the surface of the sea, can greatly influence the mixing of ocean circulation, sea surface temperature and turbulence in the whole depth, provides an energy source of the surface of the sea for the mixing in the sea, is a main physical process of exchanging substances, heat, momentum and energy in the sea, is an important factor influencing the physical and biochemical processes of the sea, and has important influence on the maintenance and change of a global climate system.
In the current wave breaking parameterization scheme, an empirical constant is generally adopted for the value estimation of the wave breaking parameters, and a fixed value is adopted as the wave breaking parameters which are not consistent with the complex wave breaking phenomenon actually occurring in the sea, so that the calculation of the turbulence energy balance of a mixing layer on the sea by wave breaking can be influenced, and the dynamic change brought by the wave breaking effect cannot be reflected. Compared with the empirical constant, the influence of wave-induced turbulence on the sea mixing process and the power process of the Bohai sea area can be more accurately expressed by utilizing the wave-ocean current coupling model to carry out coupling calculation on the wave-induced mixing caused by wave breaking.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a method for calculating real-time wave breaking-induced mixing in the ocean so as to realize dynamic simulation of a physical ocean model and improve the defect that wave breaking parameters are set to be constant values in the original physical ocean model.
The technical scheme adopted by the invention is as follows: a dynamic simulation method for calculating real-time wave breaking-caused mixing comprises the following steps:
step 1, calculating a wave field in a certain time period through an offshore wave mode, and outputting wave parameters required for calculating wave breaking parameters;
step 2, establishing a wave breaking parameter calculation model, taking the wave parameters obtained in the step 1 as the input of the wave breaking parameter calculation model, and calculating through the wave breaking parameter calculation model to obtain a wave breaking parameter field;
step 3, improving the Princeton ocean model to enable the improved Princeton ocean model to comprise an unsteady wave breaking parameterization calculation method;
and 4, taking the wave breaking parameter field obtained in the step 2 as the input of the improved Princeton ocean model, simulating an ocean dynamic process containing unsteady wave breaking in real time through the improved Princeton ocean model, outputting ocean hydrological element data serving as important parameters for analyzing an ocean physical process and an ocean biochemical process, analyzing an ocean dynamic process, and predicting and forecasting the ocean environment disaster.
Further, in step 1, the wave parameters required for calculating the wave breaking parameters include wavelength, wave spectrum peak period and sea surface effective wave height.
Further, in step 2, the wave breaking parameter field comprises a wave energy factor, a Charnock number and a water side friction speed; the wave breaking parameter calculation model comprises:
calculating the phase velocity C of the wave at the dominant frequency according to the formula (1)p(x,y):
Figure BDA0003164966140000021
Wherein, L (x, y) is wavelength, T (x, y) is wave spectrum peak period, and L (x, y) and T (x, y) are obtained by the near-shore wave mode output in the step 1;
calculation of wave energy factor alpha from equation (2)CB(x,y):
Figure BDA0003164966140000022
In the formula u*a(x, y) is the air side friction velocity, related to wind speed and roughness, u*a(x, y) is calculated according to equation (3):
Figure BDA0003164966140000023
wherein k is von Karman constant, and the value of k is 0.41; u shape10(x, y) is the wind speed ten meters above the sea surface; z is a radical of0(x, y) is the sea surface air side roughness, calculated according to equation (4):
Figure BDA0003164966140000024
in the formula, Hs(x, y) is sea surface effective wave height, and is obtained by the offshore wave mode output in the step 1;
the charnoock number β (x, y) is calculated according to equation (5):
Figure BDA0003164966140000031
calculating the water side friction velocity u according to equation (6)*w(x,y):
u*a(x,y)=(pwa)u*w(x,y) (6)
In the formula, ρwIs the sea water density, ρaFor air density, take
Figure BDA0003164966140000032
Further, the sea surface friction degree zw (x, y) is calculated according to equation (7) for calculation in the subsequent step 3:
Figure BDA0003164966140000033
wherein g is the acceleration of gravity.
Further, in step 3, the improved princeton ocean model includes:
the modified turbulence energy equation is shown in equation (8):
Figure BDA0003164966140000034
the boundary conditions of the modified turbulence energy equation are:
Figure BDA0003164966140000035
in the formula, q2(ii)/2 is the kinetic energy of turbulence; z is a vertical coordinate; t is time; s2Is the average speed of the shear to be measured,
Figure BDA0003164966140000036
u is the flow velocity in the x direction and V is the flow velocity in the y direction; n is a radical of2Is a potential density gradient of the liquid crystal,
Figure BDA0003164966140000037
ρ0is a reference to the density of the ink,
Figure BDA0003164966140000038
is the average density; b is116.6 is the model constant;
KMis a mixing coefficient of momentum, KHIs the mixing coefficient of temperature, KqIs the mixing coefficient of turbulence, KM,KH,KqThe calculation formula (2) is shown as (10):
(KM,KH,Kq)=lq(SM,SH,Sq) (10)
in the formula, SMIs KMOf a stabilization function of SHIs KHOf a stabilization function of SqIs KqOf a stabilization function of SMAnd SHIs (lN/q)2Function of (1), generally Sq=0.41SH(ii) a At the near surface (S)M,SH,Sq)=(0.30,0.49,0.20);
l is the turbulent mixing length near the sea surface, calculated from equation (11):
l=max(kzw(x,y),lz) (11)
in the formula IzIs the value of the turbulent mixing length calculated by the Mellor-Yamada2.5 order turbulence closed model in the Princeton ocean model.
Further, in step 4, the marine hydrological element data comprises marine temperature, salinity and flow velocity fields.
The invention has the beneficial effects that: the unidirectional coupling system of the invention replaces the fixed value of the wave breaking parameter with the real-time data, embodies the dynamic influence of the wave breaking effect on the ocean, improves the precision of the physical model in the process of simulating the ocean power, has important significance for more accurately analyzing the substance, heat, momentum and energy exchange in the ocean, numerous physical and biochemical processes and perfecting the ocean disaster forecasting model under extreme weather conditions, also provides powerful technical reference for correctly knowing the ocean circulation and the global climate, and lays a foundation for establishing a globalized ocean-air coupling model.
Drawings
FIG. 1: the invention discloses a flow chart of a dynamic simulation method for calculating real-time wave breaking-caused mixing.
Detailed Description
In order to further understand the contents, features and effects of the present invention, the following embodiments are illustrated and described in detail with reference to the accompanying drawings:
as shown in fig. 1, a dynamic simulation method for calculating real-time wave breaking-induced mixing includes the following steps:
step 1, a near-shore wave mode (SWAN) starts to operate, a wave field in a certain time period is calculated through the near-shore wave mode, a wave parameter data file required by calculating wave breaking parameters is output, and the operation is stopped. The wave parameters required for calculating the wave breaking parameters comprise wavelength, wave spectrum peak period and sea surface effective wave height.
SWAN is a third-generation offshore sea wave numerical mode, adopts a balance equation based on the energy conservation principle, selects a fully-implicit finite difference format, is unconditionally stable, and enables a calculation space grid and a time step length to be unlimited. In each source and sink item of the balance equation, wind input, interaction of three-phase waves and four-phase waves, bottom friction, wave breaking and white wave dissipation are considered, and the method is widely applied to wave simulation under shallow water conditions.
And 2, establishing a wave breaking parameter calculation model, starting operation after the wave breaking parameter calculation model receives the wave parameter data file, calculating to obtain a wave breaking parameter field, and stopping operation after the data file is output. Wherein, the wave breaking parameter field comprises a wave energy factor, a Charnock number and a water side friction speed.
The wave breaking parameter calculation model is a model for calculating real-time wave breaking parameters, wave parameter data such as wave height, wavelength, period and the like are required to be input, and the essence is a constant wave breaking parameterization scheme.
An unsteady wave breaking parameterization scheme:
calculating the phase velocity C of the wave at the dominant frequency according to the formula (1)p(x,y):
Figure BDA0003164966140000051
Wherein, L (x, y) is wavelength, T (x, y) is wave spectrum peak period, and L (x, y) and T (x, y) are obtained by the near-shore wave mode output in the step 1;
calculation of wave energy factor alpha from equation (2)CB(x,y):
Figure BDA0003164966140000052
In the formula u*a(x, y) is the air side friction velocity, related to wind speed and roughness, u*a(x, y) is calculated according to equation (3):
Figure BDA0003164966140000053
wherein k is von Karman constant, and the value of k is 0.41; u shape10(x, y) is the wind speed ten meters above the sea surface; z is a radical of0(x, y) is sea surface air side roughness,calculated according to equation (4):
Figure BDA0003164966140000054
in the formula, Hs(x, y) is sea surface effective wave height, and is obtained by the offshore wave mode output in the step 1;
calculating the charock number β (x, y) according to equation (5):
Figure BDA0003164966140000055
water side friction speed u*w(x, y) and air side frictional velocity u*a(x, y) is converted by equation (6):
u*a(x,y)=(ρwa)u*w(x,y) (6)
in the formula, ρwIs the sea water density, ρaFor air density, take
Figure BDA0003164966140000056
Further, the sea surface friction z is calculated from equation (7)w(x, y) for subsequent step 3 calculation:
Figure BDA0003164966140000057
wherein g is the acceleration of gravity.
And 3, improving the Princeton ocean model to enable the improved Princeton ocean model to comprise an unsteady wave breaking parameterization calculation method.
The Princeton Ocean Model (POM) is a three-dimensional oblique pressure original equation numerical value ocean mode, adopts a frog leaping finite difference format and a split operator technology, and is widely applied to the simulation of tidal current, wind-generated current, mixed layers and jump layers, hot salt circulation, ocean circulation and transportation by scholars at home and abroad. The original Princeton ocean mode is improved, and a wave breaking effect is added into a Mellor-Yamada2.5 order turbulence closed model of the Princeton ocean mode, so that the Princeton ocean mode comprises an unusual wave breaking parameterization calculation method, and the improved Mellor-Yamada2.5 order turbulence closed model is obtained.
The improved Mellor-Yamada2.5 order turbulence closed model is as follows:
on the basis of the Mellor-Yamada2.5 order turbulence closure scheme, the wave breaking effect is added by modifying the turbulence kinetic energy equation and introducing the turbulence kinetic energy input source term generated by wave breaking at the upper boundary condition of the turbulence kinetic energy equation and sea surface roughness at the upper boundary condition of the turbulence mixing length equation to determine the depth of influence of wave breaking.
The modified turbulence energy equation is shown in equation (8):
Figure BDA0003164966140000061
the boundary conditions of the modified turbulence energy equation are:
Figure BDA0003164966140000062
in the formula, q2(ii)/2 is the kinetic energy of turbulence; z is a vertical coordinate; t is time; s2Is the average speed of the shear to be measured,
Figure BDA0003164966140000063
u is the flow velocity in the x direction and V is the flow velocity in the y direction; n is a radical of2Is a potential density gradient of the liquid crystal,
Figure BDA0003164966140000064
ρ0is a reference to the density of the ink,
Figure BDA0003164966140000065
is the average density; b is116.6 is the model constant;
KMis a mixing coefficient of momentum, KHIs the mixing coefficient of temperature, KqIs the mixing coefficient of turbulence, KM,KH,KqThe calculation formula (2) is shown as (10):
(KM,KH,Kq)=lq(SM,SH,Sq) (10)
in the formula, SMIs KMOf a stabilization function of SHIs KHOf a stabilization function of SqIs KqOf a stabilization function of SMAnd SHIs (lN/q)2Function of (1), generally Sq=0.41SH(ii) a At the near surface (S)M,SH,Sq)=(0.30,0.49,0.20);
l is the turbulent mixing length near the sea surface, calculated from equation (11):
l=max(kzw(x,y),lz) (11)
in the formula IzIs the value of the turbulent mixing length calculated by the Mellor-Yamada2.5 order turbulence closed model in the Princeton ocean model.
And 4, after receiving the wave breaking parameter file, the improved Princeton ocean model starts to operate, simulates the ocean dynamic process including unsteady wave breaking in real time, outputs the ocean hydrological element data (including ocean temperature, salinity, flow velocity field and the like) file, and then stops operating.
When the improved Princeton ocean model is operated, the influence of the wave breaking effect on the ocean circulation process is reflected by the output ocean hydrological element data. Compared with the method without considering the wave breaking process or only considering the steady wave breaking process, the method provided by the invention is more suitable for the physical process occurring in the actual ocean, embodies the dynamic influence of wave breaking on the ocean, and can perfect the physical mechanism of an ocean numerical mode, thereby greatly improving the simulation and forecast precision of the model, and having important significance for researching and analyzing substances, heat, momentum, energy exchange and numerous physical and biochemical processes in the ocean.
Although the preferred embodiments of the present invention have been described above with reference to the accompanying drawings, the present invention is not limited to the above-described embodiments, which are merely illustrative and not restrictive, and those skilled in the art can make many modifications without departing from the spirit and scope of the present invention as defined in the appended claims.

Claims (5)

1. A dynamic simulation method for calculating real-time wave breaking-induced mixing is characterized by comprising the following steps:
step 1, calculating a wave field in a certain time period through an offshore wave mode, and outputting wave parameters required for calculating wave breaking parameters;
step 2, establishing a wave breaking parameter calculation model, taking the wave parameters obtained in the step 1 as the input of the wave breaking parameter calculation model, and calculating through the wave breaking parameter calculation model to obtain a wave breaking parameter field;
step 3, improving the Princeton ocean model to enable the improved Princeton ocean model to comprise an unsteady wave breaking parameterization calculation method;
and 4, taking the wave breaking parameter field obtained in the step 2 as the input of the improved Princeton ocean model, simulating an ocean dynamic process containing unsteady wave breaking in real time through the improved Princeton ocean model, outputting ocean hydrological element data serving as important parameters for analyzing an ocean physical process and an ocean biochemical process, analyzing an ocean dynamic process, and predicting and forecasting the ocean environment disaster.
2. The dynamic simulation method for calculating real-time wave-breaking-induced mixing according to claim 1, wherein in step 1, the wave parameters required for calculating the wave-breaking parameters include wavelength, wave spectrum peak period and sea surface effective wave height.
3. The method of claim 1, wherein in step 2, the wave breaking parameter field comprises wave energy factor, Charnock number and water side friction speed; the wave breaking parameter calculation model comprises:
calculating the phase velocity C of the wave at the dominant frequency according to the formula (1)p(x,y):
Figure FDA0003164966130000011
Wherein, L (x, y) is wavelength, T (x, y) is wave spectrum peak period, and L (x, y) and T (x, y) are obtained by the near-shore wave mode output in the step 1;
calculation of wave energy factor alpha from equation (2)CB(x,y):
Figure FDA0003164966130000012
In the formula u*a(x, y) is the air side friction velocity, related to wind speed and roughness, u*a(x, y) is calculated according to equation (3):
Figure FDA0003164966130000021
wherein k is von Karman constant, and the value of k is 0.41; u shape10(x, y) is the wind speed ten meters above the sea surface; z is a radical of0(x, y) is the sea surface air side roughness, calculated according to equation (4):
Figure FDA0003164966130000022
in the formula, Hs(x, y) is sea surface effective wave height, and is obtained by the offshore wave mode output in the step 1;
the charnoock number β (x, y) is calculated according to equation (5):
Figure FDA0003164966130000023
calculating the water side friction velocity u according to equation (6)*w(x,y):
u*a(x,y)=(ρwa)u*w(x,y) (6)
In the formula, ρwIs the sea water density, ρaFor air density, take
Figure FDA0003164966130000024
Further, the sea surface friction z is calculated from equation (7)w(x, y) for subsequent step 3 calculation:
Figure FDA0003164966130000025
wherein g is the acceleration of gravity.
4. The method of claim 1, wherein in step 3, the improved Princeton ocean model comprises:
the modified turbulence energy equation is shown in equation (8):
Figure FDA0003164966130000026
the boundary conditions of the modified turbulence energy equation are:
Figure FDA0003164966130000027
in the formula, q2(ii)/2 is the kinetic energy of turbulence; z is a vertical coordinate; t is time; s2Is the average speed of the shear to be measured,
Figure FDA0003164966130000028
u is the flow velocity in the x direction and V is the flow in the y directionSpeed; n is a radical of2Is a potential density gradient of the liquid crystal,
Figure FDA0003164966130000029
ρ0is a reference to the density of the ink,
Figure FDA00031649661300000210
is the average density; b is116.6 is the model constant;
KMis a mixing coefficient of momentum, KHIs the mixing coefficient of temperature, KqIs the mixing coefficient of turbulence, KM,KH,KqThe calculation formula (2) is shown as (10):
(KM,KH,Kq)=lq(SM,SH,Sq) (10)
in the formula, SMIs KMOf a stabilization function of SHIs KHOf a stabilization function of SqIs KqOf a stabilization function of SMAnd SHIs (lN/q)2Function of (1), generally Sq=0.41SH(ii) a At the near surface (S)M,SH,Sq)=(0.30,0.49,0.20);
l is the turbulent mixing length near the sea surface, calculated from equation (11):
l=max(kzw(x,y),lz) (11)
in the formula IzIs the value of the turbulent mixing length calculated by the Mellor-Yamada2.5 order turbulence closed model in the Princeton ocean model.
5. The dynamic simulation method for calculating real-time wave-breaking-induced mixing of claim 1, wherein in step 4, the marine hydrological element data includes marine temperature, salinity and flow velocity field.
CN202110801898.2A 2021-07-15 2021-07-15 Dynamic simulation method for calculating real-time wave breaking-caused mixing Pending CN113408179A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110801898.2A CN113408179A (en) 2021-07-15 2021-07-15 Dynamic simulation method for calculating real-time wave breaking-caused mixing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110801898.2A CN113408179A (en) 2021-07-15 2021-07-15 Dynamic simulation method for calculating real-time wave breaking-caused mixing

Publications (1)

Publication Number Publication Date
CN113408179A true CN113408179A (en) 2021-09-17

Family

ID=77686502

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110801898.2A Pending CN113408179A (en) 2021-07-15 2021-07-15 Dynamic simulation method for calculating real-time wave breaking-caused mixing

Country Status (1)

Country Link
CN (1) CN113408179A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116702659A (en) * 2023-08-08 2023-09-05 自然资源部第一海洋研究所 Rapid forecasting method and system for sea wave crushing parameters

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106021629A (en) * 2015-07-07 2016-10-12 国家海洋局海洋减灾中心 Method and device for establishing security risk evaluation system in ocean and coastal environment
CN110096792A (en) * 2019-04-28 2019-08-06 天津大学 A kind of Dynamic Simulation Method calculating unsteady Langmuir circulation
CN111241756A (en) * 2020-01-08 2020-06-05 大连理工大学 Numerical wave water tank wave making method based on Fluent UDF
CN111597506A (en) * 2020-04-26 2020-08-28 河海大学 Prediction method for near-shore wave breaking parameters and wave height

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106021629A (en) * 2015-07-07 2016-10-12 国家海洋局海洋减灾中心 Method and device for establishing security risk evaluation system in ocean and coastal environment
CN110096792A (en) * 2019-04-28 2019-08-06 天津大学 A kind of Dynamic Simulation Method calculating unsteady Langmuir circulation
CN111241756A (en) * 2020-01-08 2020-06-05 大连理工大学 Numerical wave water tank wave making method based on Fluent UDF
CN111597506A (en) * 2020-04-26 2020-08-28 河海大学 Prediction method for near-shore wave breaking parameters and wave height

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GEORGE MELLOR等: "Wave Breaking and Ocean Surface Layer Thermal Response", 《PHYSICAL OCEANOGRAPHY》 *
M. BAKHODAY PASKYABI等: "The influence of surface gravity waves on the injection of turbulence in the upper ocean", 《NONLINEAR PROCESS IN GEOPHYSICS》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116702659A (en) * 2023-08-08 2023-09-05 自然资源部第一海洋研究所 Rapid forecasting method and system for sea wave crushing parameters
CN116702659B (en) * 2023-08-08 2023-11-10 自然资源部第一海洋研究所 Rapid forecasting method and system for sea wave crushing parameters

Similar Documents

Publication Publication Date Title
Kantha et al. On the effect of surface gravity waves on mixing in the oceanic mixed layer
Hodges et al. Modeling basin‐scale internal waves in a stratified lake
Vermeer et al. Wind turbine wake aerodynamics
Castro et al. A mechanistic model of bubble entrainment in turbulent free surface flows
CN110096792B (en) Dynamic simulation method for calculating unsteady Langmuir circulation
Sullivan et al. Turbulent flow over water waves in the presence of stratification
Chen et al. Enhancement of wind stress evaluation method under storm conditions
Jenkins et al. Modelling the effect of ocean waves on the atmospheric and ocean boundary layers
St. Laurent et al. AN INTROUDUCTION TO THE SPECIAL ISSUE ON INTERNAL WAVES
CN113408179A (en) Dynamic simulation method for calculating real-time wave breaking-caused mixing
Ruan et al. The evolution and arrest of a turbulent stratified oceanic bottom boundary layer over a slope: Downslope regime
Skyllingstad et al. A numerical modeling study of warm offshore flow over cool water
Moghimi et al. On the modeling of wave-enhanced turbulence nearshore
Zhang et al. Numerical simulation of turbulent natural convection heat transfer in an MW-class offshore wind turbine nacelle based on a multi-feature acquisition meshing technique
Meirink et al. Modelling low-Reynolds-number effects in the turbulent air flow over water waves
CN110490354B (en) Optimization method for calculating tidal front position based on ROMS simulation result
CN115564150A (en) Oil spill track prediction method based on stokes drift and tide mixing mechanism
Ruan et al. Bottom boundary potential vorticity injection from an oscillating flow: A PV pump
Nam et al. Modeling undertow due to random waves
Meirink et al. Note on the growth rate of water waves propagating at an arbitrary angle to the wind
Ben Ali et al. Mean wind flow reconstruction of a high-rise building based on variational data assimilation using sparse pressure measurements
Zhang et al. Experimental and computational aerodynamic investigations of very thick wind turbine airfoils
Liu et al. Global estimates of wind energy input to subinertial motions in the Ekman-Stokes layer
Zhou et al. Experimental study on the performance of a wave pump for artificial upwelling in irregular waves
Curto et al. Large scale circulations in shallow lakes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20210917

RJ01 Rejection of invention patent application after publication