CN110008509A - A kind of interior estimates force characteristics analysis method under consideration Background Flow Field - Google Patents
A kind of interior estimates force characteristics analysis method under consideration Background Flow Field Download PDFInfo
- Publication number
- CN110008509A CN110008509A CN201910154630.7A CN201910154630A CN110008509A CN 110008509 A CN110008509 A CN 110008509A CN 201910154630 A CN201910154630 A CN 201910154630A CN 110008509 A CN110008509 A CN 110008509A
- Authority
- CN
- China
- Prior art keywords
- wave
- interior
- ocean
- flow field
- numerical
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of interior estimates force characteristics analysis method considered under Background Flow Field, comprising the following steps: step 1: utilizing three-dimensional ocean numerical model, carries out numerical simulation to Circulation In The South China Sea；Step 2: wave number value sink in establishing；Step 3: determining effect of the interior wave to ocean structure by considering interior wave number value sink and the ocean structure network model of background process.The present invention solves the problems, such as that interior wave wave making method can not consider Background Flow Field, interior wave ocean model, interior wave Fluid Mechanics Computation method and interior wave force theoretical method are combined for the first time, form the interior wave wave making method of a complete consideration Background Flow Field and the counting system of force measuring method.This method can make interior wave to the Effect study of works closer to actual conditions, and the charge of surety of offshore structures is calculated and structure design can play directive function.
Description
Technical field
The present invention relates to technical field of image processing more particularly to a kind of interior estimates active forces considered under Background Flow Field
Characteristic analysis method.
Background technique
With the continuous development of offshore drilling and the continuous exploitation of offshore oil field, people have looked at deep-sea.China
South China Sea has oil and gas reserve abundant, it has also become the strategic resources base of China's offshore oil natural gas production and exploitation
Ground.However there is also complicated seawater stratification phenomenon, interior wave activity are very frequent at the South Sea.
Internal wave of ocean is to marine communications and transportation and marine engineering building (including ocean platform, marine riser system, sea
Bottom pipeline system, mooring and hawser, oil storage wheel etc.) there are tremendous influence and threat.On April 10th, the 1963, " long-tail in the U.S.
The accident of shark " nuclear submarine, 129 personnel are all wrecked, are the interior estimates due to large amplitude in a short time by it vertical
It draws to caused by benthos in direction.19th-century seventies, there is document announcement, by monitoring in 4 months, the petroleum of An Demanhai
Trepan is except being twisted 90 degree and elapsing to more than 30 meters after Internal soliton passes through；Nineteen ninety, the morning in stream caul-fat field
During phase extended testing system, just once occurred because interior estimates generate abrupt strong current due to cause hawser to be broken, hull collision, very
To the accident (Chen Jinghui, 1996) for breaking and racking flotation hose.Prolong early stage the oil field of the South Sea Lufeng July 14 nineteen ninety
During long test, the abrupt strong current that interior estimates generate makes the oil tanker " Ayer of semi-submersible drilling ship " South Sea VI " number with anchoring
Difficulty occurs when connecting oil pipeline for Biru "；Followed by, at August 12nd, within the time less than 5 minutes, high current makes the oil again
Wheel swings 110 degree of angles.It is numerous it turns out that, interior wave have become restrict gas and oil in sea exploitation and navigation of submarine safety weight
Factor is wanted, the harm of interior wave must be taken into consideration in the design and installation of the ocean engineering structure of China's South China Sea.
Interior estimates are one of South Sea activity interior wave phenomenons the most frequent.Luzon Strait is as the connection Pacific Ocean and the South Sea
Main thoroughfare, when astronomical tide is passed to Northern Part of South China Sea from the Pacific Ocean, since the depth of water is shallower at Luzon Strait, the max-flow of trend
Speed is more than 1m/s.When so strong astronomical tide flows through Luzon Strait variation landform, the violent signal that baroclines can be generated, is global side
The maximum interior damp generation source (Jan etc., 2008) in edge sea.Interior tide generates interior estimates by non-linear fission in communication process,
South Sea deep-sea basin is passed through, continental shelf slope area is passed to and is finally crushed in coastal waters.Its maximum horizontal flow velocity be more than 2m/s (Zhao etc.,
2012), peak swing is more than 150m (Klymak etc., 2006).
The application of remote sensing observations means can directly acquire the global information of wave characteristic physical quantity plane distribution in extra large surface,
Interior estimates group can be captured.Hsu etc. (2000) was given using 1993 to 2000 hundreds of width ERS-1/2SAR images
Wavelength-division Butut in Northern Part of South China Sea.The results show that the interior main distributed areas of wave are Luzon Strait, near the Dongsha Islands and Hainan Island
To the east of sea area.Zheng etc. (2007) using nineteen ninety-five to SAR statistics in the 2001 Northern Part of South China Sea solitary wave frequency of occurrences,
It was found that northern South China Sea throughout the year in have an interior estimates appearance, but there is significant seasonal varieties, lonely in summer
Vertical wave frequency of occurrence is most, and winter frequency of occurrence is minimum.
In the limited deep LAYER FLUID situation of bottom flat, interior estimates keep its waveform not in its communication process
Become, referred to as stationary state interior estimates, it can be using KdV (Korteweg-de Vries), eKdV (extended KdV) and MCC
(Miyata-Choi-Camassa) (Helfrich and Melville, 2006) is described in scheduling theory.KdV theory is established
Under conditions of small nonlinearity, weak dispersion and the two balance.In Mid-low latitude sea area, the effect of geotrophic effect is just displayed, and is drawn
Enter to rotate item, obtains rKdV (rotated-modified KdV).The presence of rotation item enhances frequency dispersion effect, has widened isolated
The distance between soliton in the wave wave train reduces the number (Holloway et al., 1999) of soliton in the wave train.Simultaneously
Consider high-order nonlinear and geotrophic effect, just constitutes reKdV (rotated-modified extended KdV).
Experimental study shows that in small amplitude, KdV theory is consistent with experimental result, but in large amplitude situation,
KdV theory is not obviously inconsistent with experimental result.This is because with the increase of interior estimates amplitude, the influence of high-order nonlinear item by
Cumulative big, nonlinear terms and dispersivity no longer keep balance (Walker et al., 2003).In addition, early stage KdV equation and changing
Two layers of ocean model is used into form majority.It observes and studies have shown that in the simulation effect of two layers of ocean model and practical ocean
Observed result deviation it is larger, such as simulation nonlinear velocity c it is less than normal, interior estimates amplitude is less than normal, half band-width Δ is less than normal.
To solve these problems, people start using continuous stratification ocean model, and consider the variation of Background Flow Field in a model, will
Classical KdV equation is modified to general KdV (Generalized KdV, abbreviation GKdV).Fan et al. (2008a, b) benefit
It is found during studying the damp propagation with interior estimates in Northern Part of South China Sea with GKdV, the stream of background positive pressure flood current and ebb current
To and vertical characteristics there is significant difference.They point out, in broad continental shelf sea area, interior tidal wave and interior estimates are propagated through
The effect of Background Flow Field is highly important in journey, and the different role for distinguishing background positive pressure ebb current and flood current is necessary
's.
Currently, interior wave to the Effect study of ocean structure also in the gradually exploratory stage, about interior wave mechanism of action and shadow
The understanding for ringing mechanism is still limited.Yeung etc. (1999) is used to have been solved based on the integral equation method that singular point is distributed and be restricted water supply
The radiation and diffraction problems of rectangular configuration object in deep LAYER FLUID.Cai etc. (2003,2008) is according to Northern Part of South China Sea interior estimates
Field data combination weight coefficient analysis mode obtains the upper layer flow field Bo Zhi when interior wave passes through, and proposes to utilize Morison formula
Interior wave is calculated to the active force and torque of seabed pile, it is found that interior estimates are far longer than surface wave to the active force of pile
Effect.Ye Chunsheng etc. (2005) introduces surface wave Morison formula, and interior wave acts on small ruler when having solved different water depth and frequency
The load on cylindrical body is spent, the characteristic of load vertical characteristics is given.Wei Gang etc. (2007) uses asymptotic matching process, has studied
Reflection of the interior estimates on submerged floating vertical thin and transmission rule.He Jingyi etc. (2008) has studied SPAR in LAYER FLUID
Motor imagination problem of the platform under the effect of interior wave.Sun and Huang (2012) uses flexible member Morison Equation Iterative
The interior wave wave force for calculating S-Spar has studied motor imagination of the S-Spar platform under the effect of interior wave.Guo etc. (2013) is based on
KdV equation has studied the extreme value response of the lower top tension standpipe of South Sea interior estimates effect.
With the rapid development of Fluid Mechanics Computation (CFD), interior wave is studied using CFD technology to the work of marine structure
With problem, also referred to as a kind of research hotspot.It is outstanding to learn first-class (2008) based on N-S establishing equation numerical value sink, to density stratification stream
Body wave field Small and Medium Sized Horizontal Cylinder active force carries out numerical simulation.It pays (2009) such as Dongming and is based on N-S equation, in conjunction with interior
Solitary wave KdV and mKdV are theoretical, have developed double push plate numerical wavemaking methods, coincide with KdV and mKdV Theory Solution good.You Yunxiang
Time-domain finite displacement movement equation is used Deng (2010), in conjunction in improved Morison formula study ocean interior estimates
The load and motor imagination of power leg platform.Liu Bitao etc. (2011) is counted using mKdV theoretical origin speed entrance wave making method
Value makes wave, is studied interior estimates and top tensioned riser interaction.Nishu et al., (2011) use Sesam
The interior wave of software study is to the motor imagination of floating structure and the influence of anchoring system.It is theoretical that (2012) such as plateau snow are based on MCC
Solution, by the secondary development to business software FLUENT, establishes a kind of interior estimates numerical value sink analogy method.Cheng Youliang etc.
(2013) be based on FLUENT software, using N-S equation and VOF method, to the generation of wave in cnoidal wave, develop and its with difference
Effect of the shape pier stud under different arrangement modes has carried out numerical simulation study and analysis.
As described above, in the numerically modeling of interior wave, one is the numerical model researchs under large scale, are mainly used for mould
The generation propagation of interior wave under quasi- practical marine environment and dissipation process, result is closer to the interior wave environment under physical condition
(considering the factors such as Background Flow Field), but mesh-density is often larger, can not calculate interior wave for the operative condition of marine structure.
Another kind is wave sink in Fluid Mechanics Computation under small scale, be mainly used for calculating under the conditions of small scale the feature of interior wave and
The feature of interior wave and works interaction, it is more fine compared with large scale numerical model, but in small scale in wave sink
Wave often uses theoretical method in manufacturing, no matter its waveform or velocity field all there are larger difference, Wu Faji with actual conditions
Calculate stress condition of the structure under true interior wave environment.
Summary of the invention
It is an object of the invention to solve the problems of the above-mentioned prior art, provide a kind of by South Sea three-dimensional flow number of fields
The secondary development of value simulation and CFD technology realizes the interior wave numerical wavemaking for considering Background Flow Field, and passes through interior wave numerical experiment point
Analysis obtains mechanical characteristic of the ocean structure in interior wave environment, to provide new reference for South Sea marine engineering design.
A kind of interior estimates force characteristics analysis method under consideration Background Flow Field, comprising the following steps:
Step 1: utilizing three-dimensional ocean numerical model, numerical simulation is carried out to Circulation In The South China Sea；
Step 2: wave number value sink in establishing；
Step 3: determining interior wave to sea by considering interior wave number value sink and the ocean structure network model of background process
The effect of foreign structure.
Further, analysis method as described above, the step 1 include: based on the actual temperature in the South Sea, salinity, depth
The basic datas such as degree and atmosphere forced, precipitation, heat flux count Circulation In The South China Sea using three-dimensional ocean numerical model ROMS
Value simulation, obtains the three-dimensional flow field structure of South Sea continuous stratification.
Further, analysis method as described above, the step 2 include: using Fluid Mechanics Computation (CFD) technology,
Using UDF secondary exploitation technology, two-dimensional flume numerical wavemaking technique study is carried out, establishes interior wave number value sink.In order to be flowed
The waveform of interior wave on body interface, using the variation of VOF method tracking inner boundary waveform.Simultaneously using in the progress of numerical dissipation model
Wave wave absorption, using RANS equation as governing equation, in conjunction with k- ε turbulent model, using CFD approach to interior estimates in numerical value sink
Propagate into row numerical simulation.
Further, analysis method as described above, the step 3 include: to be simulated based on three-dimensional ocean numerical model
For the three-dimensional flow field of the continuous stratification arrived as ambient field, the upper mixed layer and lower deep-sea interval velocity induced carries out vector respectively
It is average, using the Mean Horizontal Velocity in this LAYER FLUID as movement velocity, define the speed at the condition of speed entrance boundary, fortune
It is carried out making wave with speed entrance wave making method, building considers the interior wave number value sink of Background Flow Field；
According to the crucial geometric parameter of ocean structure, suitable grid is chosen, establishes ocean structure surface model.Using interior
Wave number value sink analyzes the distribution characteristics of wave velocity field in ocean structure surface；And Morison formula is combined, internal wave effect
The surface pressure distribution and its load character of lower ocean structure carry out numerical simulation and analysis, realize interior under consideration Background Flow Field
The analysis of solitary wave force characteristics.
The present invention solves the problems, such as that interior wave wave making method can not consider Background Flow Field, for the first time by interior wave ocean model, interior
Wave Fluid Mechanics Computation method and interior wave force theoretical method combine, and form a complete consideration Background Flow Field
The counting system of interior wave wave making method and force measuring method.This method can make interior wave to the Effect study of works closer to reality
Situation, calculates the charge of surety of offshore structures and structure design can play directive function.
Detailed description of the invention
Fig. 1 is analysis method flow chart of the present invention.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, the technical solution below in the present invention carries out clear
Chu is fully described by, it is clear that described embodiments are some of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Large scale is confined to for interior wave research method at present or small scale is studied, the combination of theoretical research and numerically modeling
Degree is very low.Large scale numerical model can be calculated there is set forth herein one kind with small scale numerical value and theoretical method combines
Algorithm, to effectively manufacture the interior wave environment in practical marine environment using numerical method under small scale environment, and sharp
With the theoretical method for improving Morison formula, the load that internal wave generates is calculated.
Interior estimates force characteristics analysis method under consideration Background Flow Field provided by the invention, comprising the following steps:
Step 1: utilizing three-dimensional ocean numerical model, numerical simulation is carried out to Circulation In The South China Sea
Specifically, utilizing meteorological date provides atmosphere force data for ROMS mode；Using the temperature of Climatological, salinity,
The marine atmospheres data such as precipitation, heat flux provide boundary field and initial fields for ROMS mode；Utilize CHINESE OFFSHORE tide number
The harmonic constant for 8 main partial tides that the harmonic analysis of tide result that value assimilation model is calculated provides provides tide and forces number
According to；South Sea Sea Current numerical model data are completed by Linux work station to calculate
Step 2: wave number value sink in establishing；
Specifically, the foundation of numerical value sink is carried out as CFD platform using FLUENT software.But FLUENT software itself
Do not have interior wave wave generation.Therefore we need the open source interface by the software, i.e. user's custom feature (UDF) is right
FLUENT carries out secondary development, considers that the interior wave of Background Flow Field makes wave to realize.Concrete operations are, using ROMS to large scale sea
The interior wave environmental parameter of midocean is simulated, and the interior wave velocity field with Background Flow Field is simulated.Then extract target area
Velocity field time series, and the velocity field time series extracted is compiled into the c linguistic form that can be identified by FLUENT.It
Afterwards using user's custom feature of FLUENT using the velocity field of compiling as speed entrance boundary conditional definition numerical value sink
Inlet speed, to realize that the numerical value sink of interior estimates makes wave.In addition, the speed that the numerical value sink is used in addition to inlet
Outside the condition of entrance boundary, upper and lower right boundary is all made of the boundary Gu Bi, and exit uses Outlet boundary condition setting.
Step 3: determining interior wave to sea by considering interior wave number value sink and the ocean structure network model of background process
The effect of foreign structure.
Specifically, consider to illustrate in the principle previous step of Background Flow Field.Ocean structure grid model refers to basis
The different shape of ocean structure divides various forms of grids, including regular grid and irregular grid etc..Finally according to knot
The speed of grid node substitutes into Morison formula to calculate under background stream environment interior wave to the active force of structure around structure.
Details elaboration is done to technical solution of the present invention below:
Specifically includes the following steps:
(1) it based on basic datas such as the actual temperature in the South Sea, salinity, depth and atmosphere forced, precipitation, heat fluxs, utilizes
Three-dimensional ocean numerical model carries out numerical simulation to Circulation In The South China Sea, obtains the three-dimensional flow field structure of South Sea continuous stratification.Described three
Tieing up ocean numerical models is ROMS.
(2) Fluid Mechanics Computation (CFD) technology is used, using UDF secondary exploitation technology, carries out two-dimensional flume numerical wavemaking
Technique study establishes interior wave number value sink.In order to obtain the waveform of interior wave on fluid boundary, inner boundary is tracked using VOF method
The variation of waveform.Interior wave wave absorption is carried out using numerical dissipation model simultaneously, using RANS equation as governing equation, in conjunction with k- ε turbulent flow
Mode propagates into row numerical simulation in numerical value sink to interior estimates using CFD approach.
Specifically, the foundation of numerical value sink is carried out as CFD platform using FLUENT software.But FLUENT software itself
Do not have interior wave wave generation.Therefore we need the open source interface by the software, i.e. user's custom feature (UDF) is right
FLUENT carries out secondary development, considers that the interior wave of Background Flow Field makes wave to realize.Concrete operations are, using ROMS to large scale sea
The interior wave environmental parameter of midocean is simulated, and the interior wave velocity field with Background Flow Field is simulated.Then extract target area
Velocity field time series, and the velocity field time series extracted is compiled into the c linguistic form that can be identified by FLUENT.It
Afterwards using user's custom feature of FLUENT using the velocity field of compiling as speed entrance boundary conditional definition numerical value sink
Inlet speed, to realize that the numerical value sink of interior estimates makes wave.In addition, the speed that the numerical value sink is used in addition to inlet
Outside the condition of entrance boundary, upper and lower right boundary is all made of the boundary Gu Bi, and exit uses Outlet boundary condition setting.
VOF is the abbreviation of phase volume fraction method, is referred in the discrete control domain of Fluid Mechanics Computation, not by two kinds
The interface of commaterial phase or same substance different densities interfaces is marked and is identified by mathematical measure.Concrete principle can be with
Summary are as follows: in a certain grid control area, the definition of substance phase volume fraction is carried out to the region, fractional value is 0~1,0 and 1
Respectively represent two different phases.The phase volume fraction at the interface that then two kinds of substances have a common boundary should be just some between 0~1
Value by analyzing the size of this score, and finds using interpolation method the specific location at interface.When all grids are all defined
Afterwards, a continuous interface form can be embodied by VOF method.
VOF method is the method that a kind of most common multiphase flow numerical value calculates.In this method, since interior wave results from
Density stratification water body, therefore the water body density on both sides is not identical above and below pycnocline, we utilize VOF method by upper and lower level water
Body is defined as two different fluid phases, can all track the interface between two-phase using VOF method in entire calculate.Its
Purpose first is that in order to record the waveform at interface variation, but in order to determine in wave velocity field, as we know that, interior wave upper and lower level
Speed is of different sizes contrary, if VELOCITY DISTRIBUTION will be unable to learn without the determination at interface.
(3) three-dimensional flow field for the continuous stratification simulated based on three-dimensional ocean numerical model is lured as ambient field
The upper mixed layer and lower deep-sea interval velocity led carry out vector average respectively, using the Mean Horizontal Velocity in this LAYER FLUID as fortune
Dynamic speed, defines the speed at the condition of speed entrance boundary, carries out making wave with speed entrance wave making method, building considers background
The interior wave number value sink in flow field.
(4) according to the crucial geometric parameter of ocean structure, suitable grid is chosen, establishes ocean structure surface model.Benefit
With interior wave number value sink, the distribution characteristics of wave velocity field in ocean structure surface is analyzed；And combine Morison formula, internal wave
The surface pressure distribution and its load character for acting on lower ocean structure carry out numerical simulation and analysis, realize and consider under Background Flow Field
Interior estimates force characteristics analysis.
The ocean ROMS numerical model
ROMS full name is that Regional ocean model system (Regional Ocean Modeling System) is an open source
3D region ocean model, be the original-party of baroclining on a kind of free sea that new development is got up, terrain following, three dimensional non-linear
Journey mode.It can be effective compared with horizontal pressure force gradient algorithm using new higher order horizontal barometric gradient algorithm in ROMS mode
Reduction mode computation bring error accumulation.ROMS ocean model is joined herein in connection with many new technologies including new sub- grid
Numberization scheme and many other new modules.
(1) ROMS mode main feature:
1. assumed using static(al) and Boussinesq be approximate, inside and outside mode splitting technique, enable outer mode using compared with
Short time step, the physical phenomenon of template surface smaller scale.
2. mode uses Free Surface, closing to reality.
3. the use of σ (patch bottom) coordinate system, is capable of handling significant topography variation, such as on continental shelf slope, sill
Place.This is critically important by the shallow sea physical oceanography procedure of tide driving to research.
4. parallel (MPI, OPENMP) operation is supported to obtain mode can on various supercomputers preferably
Performance executes large-scale mode computation task.
ROMS mode capabilities are powerful, have a wide range of application, which has been widely used, and to study ocean circulation, sea-gas phase mutual
Effect, interior wave, mixing etc..Mode also has the function of some other extension, such as sea ice module packet, wave module wave, big
Gas module packet, open boundaryconditions module packet etc. allow the mode to be applied to corresponding research by simply modifying.
(2) governing equation group
The equation of motion is
State equation is
ρ=ρ (T, S, P) (4)
Temperature control equation
Salinity governing equation
Continuity equation
By the parametrization of Reynolds stress and turbulent flow tracer flux come so that equation group can solve
In formula, Du, Dv, DT, DS are diffusion terms；Fu, Fv, FT, FS are forced terms；F (x, y) is Coriolis force parameter；G is weight
Power acceleration；φ (x, y, z, t) is dynamic pressure, and value is equal to P/ ρ0, P is total pressure；ρ0+ ρ (x, y, z, t) is gross density；S(x,y,
Z, t) it is salinity；T (x, y, z, t) is position temperature；ζ (x, y, z, t) is surface water level；KM is the vertical turbulent flow coefficient of viscosity；KC is vertical
Turbulent diffusivity.
When in order to save machine, ROMS uses modal cutoff technology, and the equation of motion is decomposed into the equation of motion of vertical integral
The equation of motion (three dimensional pattern) of (two-dimensional model) and three-dimensional structure.It in this way can be by the two dimensional equation of calculated free surface from whole
It is a
It is separated in equation group, using lesser time step, biggish step-length is used to the calculating of three-dimensional flow field.?
That is integrating several step two dimensional equation groups, a step three-dimensional equation group is then integrated.
(3) boundary condition
Vertical boundary condition
As z=ζ (x, y, t), mode upper boundary conditions are as follows:
As z=-h (x, y), mode downstream condition are as follows:
In formula, Ts x, Ts yIt is extra large table wind-stress；QTIt is sea surface heat flux；E-P is the difference of evaporation capacity and precipitation；TrefIt is sea
Table reference temperature.
It is zero that horizontal lateral boundary conditions, which are normal velocity component in the boundary condition of Gu Bichu,.In the boundary of south-north direction
U=0, and in east-west boundary V=0.
(4) Sigma coordinate system
In order to improve computational efficiency, waste is avoided to calculate grid, ROMS horizontal direction is fitted using orthonormal curvilinear grid
The irregular variation in coastline, vertical direction use Sigma coordinate system, will be vertical by vertical transforming function transformation function and stretch function
Coordinate scale controls in -1~0 range, and the density of adjustable interlayer, is carried out uniformly or non-uniformly according to research contents
Layering.There are following transformational relations between the vertical physical coordinates z of cartesian coordinate system and the vertical coordinate s of Sigma coordinate system:
Z=z (x, y, s) (22)
In formula, z is the height of cartesian coordinate system, and s is the vertical distance apart from the water surface, for example, 1≤s≤0, s=
0 represents at Free Surface, z=η；S=-1 represents seabed, z=-h (x, y).ROMS gives 2 transfer functions and 4 stretchings
Function selects suitable combination voluntarily collocation use for user.
Interior wave value theory
For continuous incompressible LAYER FLUID, upper and lower level fluid density can be defined as ρ respectivelyi, i=1,2, by flute
Velocity field under karr coordinate system is defined asThen the pressure in entire flow field need to meet continuity equation and move
Measure equation.
Governing equation in numerical value in wave sink can be expressed as following forms:
uix+viy+wiz=0 (23)
Subscript i=1,2 indicate upper and lower layer fluid.
Pressure and speed in above formula should meet the condition of continuity at upper and lower level fluid boundary, by fluid layering interface
It is defined as ζ, then is had:
ζt+u1ζx=w1, ζt+u2ζx=w2, p1=P2At z=ζ (27)
According to rigid lid theorem it is assumed that top free surface is assumed to be rigid surface, bottom faces are assumed to be rigid bottom,
On this boundary at two, speed flux 0:
w1(x, h1, t)=0, w2(x ,-h2, t) and=0 (28)
Sink side is defined as smooth rigid surface, and for the Glazed fence in analog physical sink, then this boundary is answered
This meets impermeable condition:
We track LAYER FLUID internal interface using Lagrangian fluid volumetric method (VOF), and this method divides fluid
Bed boundary is assumed to be a linear plane.By the way that VOF function a (x, y, z, t)=1 or 0 is arranged at subsurface gridding
Indicate different fluid items.Then the function should meet following relational expression:
at+aux+avy+awz=0 (30)
Turbulence model in calculating describes turbulent flow using recombination k- ε turbulence model (RNG k- ε), and tubulence energy k is defeated
Fortune equation and dissipative shock wave ε may be expressed as:
Wherein, pTIt is tubulence energy energy, DkAnd DεIt is dissipative term.
Active force on works is acted on for interior estimates, can by momentum conservation equation difference force and
Viscous force integrates to obtain, and therefore, in cartesian coordinate system, fluid may be expressed as: the active force of works
fx=μ ∫ [2uixx+(vx+uy)y+(wx+uz)z]dsx+∫pix dsx (32)
fy=μ ∫ [2viyy+(vx+uy)x+(wy+vz)z]dsy+∫Piy dsy (33)
fz=μ ∫ [2wizz+(wx+uz)x+(wy+vz)y]dsz+∫piz dsz (34)
Viscous force is generated by the viscous effect of fluid, and difference force is generated by the pressure of fluid, can be respectively indicated are as follows:
In formula, (fx, fy, fz) indicate three directions (x, y, z) fluid force, (sx, sy, sz) indicate structure in three directions
On front face area.
Morison formula
In practical ocean engineering, Morison etc. proposes wave force suffered by cylindrical object in seawater being decomposed into towing
Power and inertia force two algebraical sums, form be,
Wherein p (t) represents the stress load in pile unit length.Force coefficient Cd and inertia force coefficient Cm is pulled by testing
It determines, it is and closely related with series of parameters: KC (Keulegan-Carpenter) number, Reynolds number and relatively crude
Rugosity k/D, wherein k is mean roughness, and D is diameter.
For Morison formula as a kind of semi-analytic method, application range is very extensive, is applicable not only to minor diameter
Stress of the pile in regular wave field, when calculating the interior wave force that ocean structure is subject to, when column diameter and interior wave wavelength
When compared to very little, Morison formula also can be used
Script three-dimensional ocean numerical model, Fluid Mechanics Computation and Morison formula are three substantially incoherent concepts.
Research method before is also that internal wave is studied and handled respectively in three fields.But in actual conditions, it is
Effect of the interior wave to structure under practical sea situation is more accurately calculated, need in view of Background Flow Field and Nei Bo while being acted on
Situation.The present invention provides this solution to the problem, i.e., above two numerical method and a kind of theoretical method are passed through spy
Fixed algorithm and program combines, and can either consider the effect of wave in large scale, and can be in small scale situation Imitating
Effect of the interior wave to structure.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features；
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (4)
1. a kind of interior estimates force characteristics analysis method considered under Background Flow Field, which comprises the following steps:
Step 1: utilizing three-dimensional ocean numerical model, numerical simulation is carried out to Circulation In The South China Sea；
Step 2: wave number value sink in establishing；
Step 3: determining that interior wave ties ocean by considering interior wave number value sink and the ocean structure network model of background process
The effect of structure.
2. analysis method according to claim 1, which is characterized in that the step 1 includes: based on the actual temperature in the South Sea
The basic datas such as degree, salinity, depth and atmosphere forced, precipitation, heat flux, using three-dimensional ocean numerical model ROMS, to the South Sea
Circulation carries out numerical simulation, obtains the three-dimensional flow field structure of South Sea continuous stratification.
3. analysis method according to claim 1, which is characterized in that the step 2 includes: using Fluid Mechanics Computation skill
Art carries out two-dimensional flume numerical wavemaking technique study using UDF secondary exploitation technology, establishes interior wave number value sink；In order to obtain
The waveform of interior wave on fluid boundary, using the variation of VOF method tracking inner boundary waveform；It is carried out simultaneously using numerical dissipation model
Interior wave wave absorption, using RANS equation as governing equation, in conjunction with k- ε turbulent model, using CFD approach to interior estimates in numerical value sink
In propagate into row numerical simulation.
4. analysis method according to claim 1, which is characterized in that the step 3 includes: based on three-dimensional ocean Numerical-Mode
The three-dimensional flow field of the quasi- obtained continuous stratification of pattern as ambient field, distinguish by the upper mixed layer and lower deep-sea interval velocity induced
Vector average is carried out, using the Mean Horizontal Velocity in this LAYER FLUID as movement velocity, is defined at the condition of speed entrance boundary
Speed, make wave with speed entrance wave making method, building considers the interior wave number value sink of Background Flow Field；
According to the crucial geometric parameter of ocean structure, suitable grid is chosen, establishes ocean structure surface model；Utilize interior wave number
It is worth sink, analyzes the distribution characteristics of wave velocity field in ocean structure surface；And Morison formula is combined, internal wave effect is plunged into the commercial sea
The surface pressure distribution and its load character of foreign structure carry out numerical simulation and analysis, realize interior isolated under consideration Background Flow Field
Wave force specificity analysis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910154630.7A CN110008509B (en) | 2019-03-01 | 2019-03-01 | Method for analyzing internal solitary wave action force characteristics under consideration of background flow field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910154630.7A CN110008509B (en) | 2019-03-01 | 2019-03-01 | Method for analyzing internal solitary wave action force characteristics under consideration of background flow field |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110008509A true CN110008509A (en) | 2019-07-12 |
| CN110008509B CN110008509B (en) | 2020-07-07 |
Family
ID=67166212
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910154630.7A Active CN110008509B (en) | 2019-03-01 | 2019-03-01 | Method for analyzing internal solitary wave action force characteristics under consideration of background flow field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110008509B (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111079357A (en) * | 2019-10-22 | 2020-04-28 | 重庆邮电大学 | Electromagnetic interference time domain analysis method of microstrip line on PCB under action of electromagnetic wave |
| CN111275678A (en) * | 2020-01-17 | 2020-06-12 | 山东科技大学 | System and method for measuring displacement of elongated structure under action of internal waves |
| CN111353253A (en) * | 2020-03-31 | 2020-06-30 | 温州科技职业学院 | Method for realizing multi-column waves with different parameters based on Flow-3D software simulated physical wave-making means |
| CN111368429A (en) * | 2020-03-04 | 2020-07-03 | 江苏科技大学 | Method for confirming influence of terrain change on narrow slit resonance |
| CN111460724A (en) * | 2020-04-21 | 2020-07-28 | 天津大学 | Dynamic simulation method for calculating abundance of moon jellyfishes |
| CN113640800A (en) * | 2021-08-25 | 2021-11-12 | 中国人民解放军海军潜艇学院 | Inversion method for inverting isolated wave data in ocean |
| CN113742900A (en) * | 2021-08-17 | 2021-12-03 | 生态环境部华南环境科学研究所 | Device and method for simulating medium and small tide and ocean circulation |
| CN116432562A (en) * | 2023-06-13 | 2023-07-14 | 山东科技大学 | Analysis method for influence degree of density jump layer element change on internal tide |
| CN116451511A (en) * | 2023-06-16 | 2023-07-18 | 中国人民解放军国防科技大学 | Sea ice numerical simulation method, device and equipment based on ROMS model |
| CN116738136A (en) * | 2023-05-17 | 2023-09-12 | 中国海洋大学 | Calculation method for seabed methane gas release under action of internal solitary wave |
| CN117313587A (en) * | 2023-11-28 | 2023-12-29 | 西北工业大学 | Method and system for simulating interaction between internal solitary wave and background shear flow |
| CN117521556A (en) * | 2023-11-27 | 2024-02-06 | 中国民用航空飞行学院 | Submarine attitude prediction method under influence of internal solitary waves |
| CN117890080A (en) * | 2024-03-15 | 2024-04-16 | 中国海洋大学 | Method for simulating wave generation of solitary waves in shutter door rotating ocean |
| CN118153349A (en) * | 2024-05-09 | 2024-06-07 | 华南理工大学 | Attitude response solving method of anchor system structure under influence of internal solitary wave |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103063253A (en) * | 2012-12-27 | 2013-04-24 | 清华大学深圳研究生院 | Multi-sending and multi-receiving acoustic measurement method for oceanic internal waves |
| CN103454064A (en) * | 2013-09-18 | 2013-12-18 | 中国海洋大学 | First modality internal wave generator |
| CN105205200A (en) * | 2015-04-27 | 2015-12-30 | 河海大学 | Verification method for silt coast wave-induced current numerical simulation |
| CN106596046A (en) * | 2016-12-30 | 2017-04-26 | 中国海洋大学 | Stress observation device for ocean structure under internal solitary wave environment |
| CN107218931A (en) * | 2017-05-23 | 2017-09-29 | 广东贝达海洋科学有限公司 | A kind of internal wave of ocean early warning system and method based on online monitoring data |
| CN108254157A (en) * | 2018-02-06 | 2018-07-06 | 中国海洋大学 | A kind of interior wave and the experimental system of submerged body interaction |
| CN108595831A (en) * | 2018-04-22 | 2018-09-28 | 天津大学 | It is a kind of to calculate the mixed Dynamic Simulation Method of tide cause in real time |
| CN108829918A (en) * | 2018-04-28 | 2018-11-16 | 中国海洋大学 | Intelligent buoy networking simulating method and system towards oceanographic phenomena |
-
2019
- 2019-03-01 CN CN201910154630.7A patent/CN110008509B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103063253A (en) * | 2012-12-27 | 2013-04-24 | 清华大学深圳研究生院 | Multi-sending and multi-receiving acoustic measurement method for oceanic internal waves |
| CN103454064A (en) * | 2013-09-18 | 2013-12-18 | 中国海洋大学 | First modality internal wave generator |
| CN105205200A (en) * | 2015-04-27 | 2015-12-30 | 河海大学 | Verification method for silt coast wave-induced current numerical simulation |
| CN106596046A (en) * | 2016-12-30 | 2017-04-26 | 中国海洋大学 | Stress observation device for ocean structure under internal solitary wave environment |
| CN107218931A (en) * | 2017-05-23 | 2017-09-29 | 广东贝达海洋科学有限公司 | A kind of internal wave of ocean early warning system and method based on online monitoring data |
| CN108254157A (en) * | 2018-02-06 | 2018-07-06 | 中国海洋大学 | A kind of interior wave and the experimental system of submerged body interaction |
| CN108595831A (en) * | 2018-04-22 | 2018-09-28 | 天津大学 | It is a kind of to calculate the mixed Dynamic Simulation Method of tide cause in real time |
| CN108829918A (en) * | 2018-04-28 | 2018-11-16 | 中国海洋大学 | Intelligent buoy networking simulating method and system towards oceanographic phenomena |
Non-Patent Citations (5)
| Title |
|---|
| WANG Z, DONG S, DONG X, ET AL.: "Study of Surface Wave Characteristics and Extreme Parameters Estimation in a Port With Submerged Breakwater", 《ASME 2016 35TH INTERNATIONAL CONFERENCE ON OCEAN, OFFSHORE AND ARCTIC ENGINEERING》 * |
| 刘昆; 徐振华; 尹宝树: "地形和层结变化对内潮生成影响的数值模拟研究", 《海洋科学》 * |
| 王琮，刘振纹: "存在背景流的内波场中隔水管涡激振动研究", 《哈尔滨工程大学学报》 * |
| 王飞: "内孤立波作用下小尺度竖直圆柱体的水动力特性研究", 《中国博士学位论文全文数据库工程科技II辑》 * |
| 蔡树群, 何建玲, 谢皆烁: "近10年来南海孤立内波的研究进展", 《地球科学进展》 * |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111079357A (en) * | 2019-10-22 | 2020-04-28 | 重庆邮电大学 | Electromagnetic interference time domain analysis method of microstrip line on PCB under action of electromagnetic wave |
| CN111275678A (en) * | 2020-01-17 | 2020-06-12 | 山东科技大学 | System and method for measuring displacement of elongated structure under action of internal waves |
| CN111275678B (en) * | 2020-01-17 | 2023-09-29 | 山东科技大学 | System and method for measuring displacement of slender structure under action of internal wave |
| CN111368429A (en) * | 2020-03-04 | 2020-07-03 | 江苏科技大学 | Method for confirming influence of terrain change on narrow slit resonance |
| CN111353253A (en) * | 2020-03-31 | 2020-06-30 | 温州科技职业学院 | Method for realizing multi-column waves with different parameters based on Flow-3D software simulated physical wave-making means |
| CN111353253B (en) * | 2020-03-31 | 2020-12-01 | 温州科技职业学院 | Method for realizing multi-column waves with different parameters based on Flow-3D software simulated physical wave-making means |
| CN111460724A (en) * | 2020-04-21 | 2020-07-28 | 天津大学 | Dynamic simulation method for calculating abundance of moon jellyfishes |
| CN111460724B (en) * | 2020-04-21 | 2022-05-13 | 天津大学 | Dynamic simulation method for calculating abundance of moon jellyfishes |
| CN113742900A (en) * | 2021-08-17 | 2021-12-03 | 生态环境部华南环境科学研究所 | Device and method for simulating medium and small tide and ocean circulation |
| CN113742900B (en) * | 2021-08-17 | 2023-05-12 | 生态环境部华南环境科学研究所 | Medium-small tide and ocean circulation simulation device and method |
| CN113640800B (en) * | 2021-08-25 | 2023-07-28 | 中国人民解放军海军潜艇学院 | Inversion method for inverting ocean solitary wave data |
| CN113640800A (en) * | 2021-08-25 | 2021-11-12 | 中国人民解放军海军潜艇学院 | Inversion method for inverting isolated wave data in ocean |
| CN116738136B (en) * | 2023-05-17 | 2024-03-15 | 中国海洋大学 | Calculation method for seabed methane gas release under action of internal solitary wave |
| CN116738136A (en) * | 2023-05-17 | 2023-09-12 | 中国海洋大学 | Calculation method for seabed methane gas release under action of internal solitary wave |
| CN116432562A (en) * | 2023-06-13 | 2023-07-14 | 山东科技大学 | Analysis method for influence degree of density jump layer element change on internal tide |
| CN116432562B (en) * | 2023-06-13 | 2023-08-22 | 山东科技大学 | Analysis method for influence degree of density jump layer element change on internal tide |
| CN116451511B (en) * | 2023-06-16 | 2023-08-22 | 中国人民解放军国防科技大学 | Sea ice numerical simulation method, device and equipment based on ROMS model |
| CN116451511A (en) * | 2023-06-16 | 2023-07-18 | 中国人民解放军国防科技大学 | Sea ice numerical simulation method, device and equipment based on ROMS model |
| CN117521556A (en) * | 2023-11-27 | 2024-02-06 | 中国民用航空飞行学院 | Submarine attitude prediction method under influence of internal solitary waves |
| CN117313587A (en) * | 2023-11-28 | 2023-12-29 | 西北工业大学 | Method and system for simulating interaction between internal solitary wave and background shear flow |
| CN117313587B (en) * | 2023-11-28 | 2024-02-06 | 西北工业大学 | Method and system for simulating interaction between internal solitary wave and background shear flow |
| CN117890080A (en) * | 2024-03-15 | 2024-04-16 | 中国海洋大学 | Method for simulating wave generation of solitary waves in shutter door rotating ocean |
| CN117890080B (en) * | 2024-03-15 | 2024-05-28 | 中国海洋大学 | Method for simulating wave generation of solitary waves in shutter door rotating ocean |
| CN118153349A (en) * | 2024-05-09 | 2024-06-07 | 华南理工大学 | Attitude response solving method of anchor system structure under influence of internal solitary wave |
| CN118153349B (en) * | 2024-05-09 | 2024-07-16 | 华南理工大学 | Attitude response solving method of anchor system structure under influence of internal solitary wave |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110008509B (en) | 2020-07-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110008509A (en) | A kind of interior estimates force characteristics analysis method under consideration Background Flow Field | |
| Zheng et al. | A modeling study of the Satilla River estuary, Georgia. I: Flooding-drying process and water exchange over the salt marsh-estuary-shelf complex | |
| Bilskie et al. | Data and numerical analysis of astronomic tides, wind‐waves, and hurricane storm surge along the northern Gulf of Mexico | |
| Sheng et al. | Simulation of storm surge, wave, and coastal inundation in the Northeastern Gulf of Mexico region during Hurricane Ivan in 2004 | |
| Tonelli et al. | Modeling wave impact on salt marsh boundaries | |
| Sheng et al. | Simulation of storm surge, wave, currents, and inundation in the Outer Banks and Chesapeake Bay during Hurricane Isabel in 2003: The importance of waves | |
| Van der Molen et al. | Holocene tidal conditions and tide‐induced sand transport in the southern North Sea | |
| Yong et al. | Process control of the sand wave migration in Beibu Gulf of the South China Sea | |
| Johnston et al. | Energy and momentum lost to wake eddies and lee waves generated by the North Equatorial Current and tidal flows at Peleliu, Palau | |
| Soomere et al. | On dispersion properties of surface motions in the Gulf of Finland | |
| Wang et al. | Development of an integrated model system to simulate transport and fate of oil spills in seas | |
| Guo et al. | Observed deep anticyclonic cap over Caiwei Guyot | |
| Lin et al. | Responses of hydrodynamics to changes in shoreline and bathymetry in the Pearl River Estuary, China | |
| Franchello | Shoreline tracking and implicit source terms for a well balanced inundation model | |
| Liang et al. | Bottom shear stress under wave-current interaction | |
| Zhang et al. | Coupling of wave and current numerical model with unstructured quadtree grid for nearshore coastal waters | |
| Hou et al. | The morphodynamic responses to deposition-promoting projects in island and reef coasts of the Zhoushan Archipelago, China | |
| Pinheiro et al. | Boussinesq-type numerical model for wave propagation near shore and wave penetration in harbors | |
| Jiang et al. | Evaluating the effects of a symmetric instability parameterization scheme in the Xisha-Zhongsha waters, South China Sea in winter | |
| Huang et al. | Suspended sediment dynamics and influencing factors during typhoons in Hangzhou Bay, China | |
| Wang et al. | The impact of surface wave on the sediment erosion and deposition near the wellow river mouth, China. | |
| James et al. | Calibration of a 3D hydrodynamic aquaculture model | |
| Tuhaijan et al. | Wave-current interaction on offshore spar platforms | |
| Ruo-Shan et al. | Turbulent mixing of the Kuroshio waters Southeast of Taiwan | |
| Myrberg et al. | Topography, hydrography, circulation and modelling of the Baltic Sea |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |