CN108563867A - A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM - Google Patents
A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM Download PDFInfo
- Publication number
- CN108563867A CN108563867A CN201810330662.3A CN201810330662A CN108563867A CN 108563867 A CN108563867 A CN 108563867A CN 201810330662 A CN201810330662 A CN 201810330662A CN 108563867 A CN108563867 A CN 108563867A
- Authority
- CN
- China
- Prior art keywords
- grid
- wrf
- wind
- cfd
- openfoam
- 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
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
Abstract
The invention belongs to computation wind engineering fields, and disclose a kind of method for realizing WRF and CFD coupled simulation wind fields based on OpenFOAM, and the method includes the steps of：1) complicated landform CFD model is established, fining grid is generated；2) mesoscale WRF modeling complex topography wind fields are established, the wind field data of WRF models are handled；3) it realizes that WRF models and complicated landform CFD model are coupled in borderline data, obtains the fining grid wind field data changed over time on complicated landform CFD model 4 sides interface and top surface；4) the complicated landform CFD model after coupling is calculated, to obtain complex topography wind field simulated behavior result.The present invention can use the program language of computer supported, realize that WRF models and complicated landform CFD model are coupled in borderline data, be conducive to the precision for promoting complicated landform wind-resources assessment and Numerical Simulation of Wind.
Description
Technical field
The invention belongs to computation wind engineering fields, WRF and CFD couplings are realized based on OpenFOAM more particularly, to one kind
The method for closing modeling wind field.
Background technology
When studying the wind field of complicated landform, primary problem is how to give rational entrance boundary condition.Entrance boundary
Reasonability be whether wind tunnel test or numerical simulation an important factor for influencing result of calculation, in given entrance boundary
When all done a large amount of simplification, result and actual conditions deviation may be caused larger.
In recent years, with the fast development of weather forecast pattern, Meso-scale meteorology Forecast Mode is widely applied,
Middle WRF patterns are most popular.WRF is American National Center for Atmospheric Research, National Environmental forecasting centre and multiple universities, research aircraft
The meso-scale model system of new generation that structure is researched and developed jointly.WRF patterns use Arakawa C staggered-meshes in the horizontal direction
Point, vertical direction then use terrain following mass coordinate, have flexible, expansible, easy care and use computer platform extensive
The characteristics of.WRF modular systems such as improve from cloud scale to synoptic scale at different scales significant weather feature forecast precision, the emphasis
The horizontal grid for considering 1~10 kilometer, in conjunction with advanced numerical method and Data Assimilation technology, using improved physics mistake
Journey scheme, while there is multinest and be easily positioned in the ability of diverse geographic location, it can be very good to adapt to from idealization
Research to the applications such as operational forecast needs, and with convenient for further strengthening perfect flexibility.
Currently in the world about the simulation of complex topography wind field mostly using Mesoscale Meteorology, pass through nested net
Lattice are then wide by the method and atmospheric boundary layer wind of power NO emissions reduction or statistics NO emissions reduction to a certain minimum internal layer scale
Line distribution theory or semiempirical formula obtain high-resolution wind speed profile data.Then with observation website field data into
Row comparative analysis finally obtains the numerical simulation result and wind for meeting certain required precision to improved mode and analogy method
It can resource map.But Meso-scale meteorology Forecast Mode generally use landform follows coordinate, passes through the coordinate transform of equation group
Complicated landform is described, on numerical computation method based on difference scheme, and is needed to landform in the pretreatment stage of simulation
In various degree smooth is carried out, with the stability calculated, when encountering extremely precipitous landform, it is possible to integrate
The case where spilling, thus simple mesoscale model be difficult to complicated mima type microrelief wind field carry out detailed analysis, to wind energy resources into
Row accurate evaluation.
Numerical model based on CFD approach is suitable for the simulation in the fine flow field in boundary layer under MODEL OVER COMPLEX TOPOGRAPHY, has
Required expense is low, the calculating time is short, computational efficiency is high, is convenient for the features such as simulation true wind field characteristic.But limitation is, a side
Face, costly, closure is strong for CFD business softwares, and user does not know detailed calculating process, and custom feature is carried out to software
Or discrete based on governing equation be from source when programming, difficulty is high, heavy workload, very inconvenient.In this regard, OpenFOAM is
A set of complete CFD class libraries write by C++ freely increased income, OpenFOAM is not only by program code external disclosure, but also its is soft
The structure of part program and the architecture design of software also increasing income completely.Therefore user can farthest expand program to realize
The self-defined or extended software of user function is functional.In addition to this, OpenFOAM also supports a variety of grids, can handle
Complicated geometry, and parallel computation etc. is decomposed in support area.
On the other hand, it is typically only capable to hypothetical boundary inflow condition and uses mean wind speed or power exponent Wind outline state, it can not
Ensure that in numerical simulation can include correctly true mesoscale circulation information.Therefore the coupled mode of WRF and CFD is intricately
The important development direction of shape numerical simulation is conducive to the precision for promoting complicated landform wind-resources assessment and Numerical Simulation of Wind.
A kind of wind of coupling WRF and OpenFOAM patterns is disclosed in Chinese invention patent specification CN106326625A
Field stimulation method, technical solution provided by the invention realizes WRF patterns and the coupling of OpenFOAM patterns calculates, by WRF's
The data NO emissions reduction of several kilometers of horizontal resolutions is the 30m resolution datas of OpenFOAM, can improve WRF under complicated landform
Simulation of wind effect.But when the point of each series of vertical on offside boundary carries out fitting of a polynomial, fitting function is wind speed
Component u, v, w function between height z, i.e. u (z), v (z) and w (z) respectively do not account for u, v, w and plane coordinates x, y
Relationship, and have in same lateral boundaries the fitting function of different level position, each fitting function acts only on one section of level
Range causes not accounting for similar influence of multiple fitting functions to CFD Grid Velocities, difference is generated with actual conditions.And
And when OpenFOAM patterns being driven to be calculated, only with the Reynolds average stress model of stable state, this model, which can not parse, to be had
The microvortex of stronger pulsatile characteristics.
Shen refining etc. is on the basis of Mesoscale Meteorology (WRF) using multi-scale Coupling to mountain area valley bridge site wind field
Fining numerical simulation is carried out, the analog result of WRF is handled using piecemeal differential polynomial method, difference has been carried out to inlet
Quantity piecemeal, it is more in massif complex region piecemeal near the ground, and it is more sparse far from ground area dividing, it is right respectively after piecemeal
The wind field in every piece of region carries out fitting of a polynomial.Result of calculation shows that piecemeal is more, and result is coincide more with meteorological data
It is good.But piecemeal difference is rule of thumb carried out, workload is huge, and efficiency is low；And subjective experience piecemeal difference result is uncertain
Property is big.
Above-mentioned two research when being calculated, boundary condition setting in addition to speed entrance, also pressure export or other
Boundary condition is simplified, and has the case where speed to have larger difference with practical each boundary, all directions, without abundant
Play WRF and advantage of the CFD coupled mode in complicated landform numerical simulation.
Invention content
For the disadvantages described above or Improvement requirement of the prior art, the present invention provides realize WRF and CFD based on OpenFOAM
The method of coupled simulation wind field, this method can use the program language of computer supported, parametrization to establish complicated landform CFD moulds
Type establishes WRF modeling complex topography wind fields, realizes WRF models and complicated landform CFD model in borderline data coupling
It closes, is conducive to the precision for promoting complicated landform wind-resources assessment and Numerical Simulation of Wind.
To achieve the above object, it is proposed, according to the invention, provide a kind of based on OpenFOAM realization WRF and CFD coupled simulations
The method of wind field, which is characterized in that the method includes the steps of：
1) complicated landform CFD model is established based on GIS data in OpenFOAM, generates fining grid；
2) mesoscale WRF modeling complex topography wind fields are established, the wind field data of WRF models are handled；
3) it realizes that WRF models and complicated landform CFD model are coupled in borderline data, obtains complicated landform CFD model 4
The fining grid wind field data changed over time on a side interface and top surface；
4) OpenFOAM is used to calculate the complicated landform CFD model after coupling, to obtain complex topography wind field characteristic mould
Quasi- result.
Preferably, step 1) specifically includes following sub-step：
1.1) center longitude for determining target obtains the ASTER of target region from GIS GIS-Geographic Information System
Terrain data under GDEM30 meters of resolution ratio geodetic coordinates handles GIS data, to obtain WGS_1984_UTM using ArcMap
Terrain data grid node under coordinate system；
1.2) parametrization establishes three-dimensional level terrain CFD model and generates fining grid：The meter of complicated landform CFD model
It calculates domain and is set as rectangular body region, complicated landform influential on computational domain, computational domain packet are found according to actual topography and geomorphology
Containing these complicated landforms, and use unstructured prismatic mesh to be fully fitted complex geometry boundary, wherein in computational domain
Different zones using the different grid of resolution ratio, complicated landform CFD model uses adjacent mesh size ratio in the vertical direction
Value is the σ grids of definite value；
1.3) number between the terrain data grid node of computational domain and three-dimensional level terrain CFD model grid node is realized
According to conversion：According to the coordinate of each soffit grid back end in three-dimensional level terrain CFD model, it is found in terrain data
Then corresponding position in grid node uses the elevation of 3 terrain data grid nodes around the position, passes through triangle
Coplanar linear interpolation obtains the elevation of soffit grid node, meanwhile, in three-dimensional level terrain CFD model, on vertical direction its
The elevation of remaining grid node can generate respective change with the variation of soffit grid node elevation, after finally obtaining modification elevation
Complicated landform CFD model；
1.4) OpenFOAM is imported by converter gambitToFoam, complicated landform CFD model is converted to
The identifiable formats of OpenFOAM.
Preferably, step 2) specifically includes following sub-step：
2.1) mesoscale WRF modeling complex topography wind fields are established, WRF models are cuboid, and horizontal plane direction uses
Multilayer nest grid, the wherein size of the size Yu outer layer coarse grids of WRF models inner-layer thin fine grid at odd number than relationship, with
Realize two-way nested-grid ocean；
2.2) the wind field data for the different time that processing WRF models calculate, respectively obtain 4 around complicated landform CFD model
A series of coarse resolution wind field data (t, x, y, z, u, v, w) on a side interface and top surface on grid nodes, wherein t tables
Show the time, (x, y, z) indicates the three dimensional space coordinate of each grid node locations, three under (u, v, w) representation space rectangular coordinate system
The wind speed component in a direction, concrete processing procedure are as follows：
For each side interface of complicated landform CFD model or top surface, the WRF model meshes lists that will intersect first with it
Member centre coordinate be projected on corresponding side interface or top surface, obtain each grid node locations three dimensional space coordinate (x,
Y, z), then corresponding position of each grid node in WRF models is found, pass through the WRF model wind field numbers around each grid node
According to the wind speed (u, v, w) of each grid node being obtained using distance-reverse weighting function interpolation, in this way, it is possible to obtain intricately
Around shape CFD model on 4 sides interface and top surface on a series of grid nodes coarse resolution wind field data (t, x, y, z, u,
V, w).
Preferably, step 3) specifically includes following sub-step：
3.1) in OpenFOAM read boundary primary condition file, obtain complicated landform CFD model 4 sides interface and
The number of grid is refined on top surface；
3.2) for each side interface of complicated landform CFD model or top surface, first according to the fining net on the face
The number of lattice determines the rectangular space coordinate of each fining grid, then finds each fining grid node and is obtained in step 2.2)
Coarse resolution wind field back end in corresponding position, pass through it is each fining grid node around coarse resolution wind field data
Node obtains the wind speed (u of each fining grid node using distance-reverse weighting function interpolation*, v*, w*), wherein (u*, v*, w*)
The wind speed component in lower three directions of representation space rectangular coordinate system obtains 4 sides of complicated landform CFD model by this method
The fining grid wind field data changed over time on interface and top surface；
3.3) consider as a result flux conservation is attached to the wind of top surface using the total flux of computational domain divided by the area of top surface
Fast w1In, obtain final top surface wind speed；
3.4) the fining grid wind field data of adjacent time are carried out according to the time interval of update boundary condition linear
Fining of each time grid wind field data that linear interpolation obtains are stored in different files by interpolation.
Preferably, step 4) specifically includes following sub-step：
4.1) for the boundary condition of complicated landform CFD model, other than bottom surface is set as wall surface, 4 lateral boundaries of surrounding
Face and top surface are disposed as speed entrance；
4.2) first the fining grid wind field data of initial time are written in the boundary primary condition file of OpenFOAM,
And OpenFOAM is driven to calculate until wind field result is stablized, then using previous result of calculation as initial fields, according to the time
Interval after updating the boundary condition file of OpenFOAM successively and solving control file, then drives OpenFOAM to be calculated,
In, tri- kinds of large eddy simulation, k- ω or k- ε turbulence models may be used in when calculating, calculate and use PIMPLE algorithms, final to obtain more
Accurate complex topography wind field simulated behavior result.
In general, through the invention it is contemplated above technical scheme is compared with the prior art, can obtain down and show
Beneficial effect：
(1) present invention can be based on GIS data parametrization and establish complicated landform CFD model, and it is complicated to establish WRF modelings
Landform wind field realizes that WRF models and complicated landform CFD model are coupled in borderline data, is calculated and is obtained using OpenFOAM
More accurate complex topography wind field simulated behavior result.Above-mentioned steps avoid lengthy and jumbled manual operation flow, use computer branch
The program language held, it is so efficient that implement the technical solution in OpenFOAM.
(2) in the wind field data and realization WRF models and complicated landform CFD model of processing WRF models in borderline number
When according to coupling, use distance-reverse weighting function into row interpolation, Er Qiechong more efficient, simple and direct compared to curve matching
Point considers the influence of the stronger data of correlation, avoids the feelings that curve matching occurs relatively large deviation when fluctuations in wind speed is larger
Condition, precision higher；In addition, directly using the grid of complicated landform CFD model as sub- boundary, avoids and rule of thumb divided
The subjective uncertainty of block.
(3) the 4 sides interface and top surface of complicated landform CFD model are disposed as speed entrance, consider flux conservation into
Row adjustment, is more in line with actual conditions.
Description of the drawings
Fig. 1 is that the present invention is based on the flow charts of the method for OpenFOAM realization WRF and CFD coupled simulation wind fields；
Fig. 2 a and Fig. 2 b be the complicated landform CFD model established based on GIS data in OpenFOAM of the present invention respectively with
Refine grid distribution schematic diagram；
Fig. 3 is the wind field data that present invention processing WRF models calculate, and obtains 4 lateral boundaries around complicated landform CFD model
A series of schematic diagram of coarse resolution wind field data on face and top surface on grid nodes；
Fig. 4 is that the present invention realizes that WRF models and complicated landform CFD model are coupled in borderline data, is obtained intricately
The schematic diagram of the fining grid wind field data changed over time on shape CFD model 4 sides interface and top surface；
Fig. 5 is to realize WRF and CFD coupled simulation complex topography wind fields in the embodiment of the present invention based on OpenFOAM, obtain
12 hours wind speed profile change schematic diagrams.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.As long as in addition, technical characteristic involved in the various embodiments of the present invention described below
It does not constitute a conflict with each other and can be combined with each other.
Referring to Fig.1~Fig. 5, a kind of side realizing WRF and CFD coupled simulation wind fields based on OpenFOAM provided by the invention
Method, flow chart are as shown in Figure 1.
1, complicated landform CFD model is established based on GIS data in OpenFOAM, generates fining grid；
1.1) center longitude for determining target obtains the ASTER of target region from GIS GIS-Geographic Information System
Terrain data under GDEM30 meters of resolution ratio geodetic coordinates handles GIS data, to obtain WGS_1984_UTM using ArcMap
Terrain data grid node under coordinate system；
Embodiment 1 choose target center longitude coordinate be：113.328778 ° of longitude, 22.055561 ° of latitude, from
The 30 meters of resolution ratio terrain data ranges obtained in GIS GIS-Geographic Information System are：21~23 ° of north latitude, 113~114 ° of east longitude.
1.2) parametrization establishes three-dimensional level terrain CFD model and generates fining grid：Three-dimensional level terrain CFD model
Computational domain be set as rectangular body region, complicated landform influential on computational domain is found according to actual topography and geomorphology, calculate
Domain needs to include these complicated landforms, and uses unstructured prismatic mesh to be fully fitted complex geometry boundary, wherein
In the grid that the different zones of computational domain use resolution ratio different, three-dimensional level terrain CFD model uses phase in the vertical direction
Adjacent size of mesh opening ratio is the σ grids of definite value；
In embodiment 1, computational domain is set as rectangular body region, and length, width and height take 100km, 100km and 20km, use non-respectively
Structuring prismatic mesh, unstructured grid refer to that the internal point in net region does not have identical adjacent cells, parameter
Change control grid density distribution, fine and closely woven size of mesh opening is 20m, and maximum mesh size is 100m, and three-dimensional level terrain CFD model exists
Vertical direction carry out mesh refinement near the ground, minimum grid size are 1m, take equal ratio to divide upwards, and grid is most in computational domain
Big growth rate is 1.2；
1.3) number between the terrain data grid node of computational domain and three-dimensional level terrain CFD model grid node is realized
According to conversion：According to the coordinate of each soffit grid back end in three-dimensional level terrain CFD model, it is found in terrain data
Then corresponding position in grid node uses the elevation of 3 terrain data grid nodes around the position, passes through triangle
Coplanar linear interpolation obtains the elevation of soffit grid node, meanwhile, in three-dimensional level terrain CFD model, on vertical direction its
The elevation of remaining grid node can generate respective change with the variation of soffit grid node elevation, after finally obtaining modification elevation
Complicated landform CFD model, as shown in Fig. 2 a, Fig. 2 b；
1.4) OpenFOAM is imported by converter gambitToFoam, complicated landform CFD model is converted to
The format that OpenFOAM can be identified.
2, the wind field data calculated using the coarse resolution wind field in WRF simulated targets region, processing WRF；
2.1) mesoscale WRF modeling complex topography wind fields are established, WRF models are cuboid, and horizontal plane direction uses
Multilayer nest grid, wherein the size of WRF model inner-layer thin fine grids and the size of outer layer coarse grids at odd number than relationship,
Realize two-way nested-grid ocean；
In embodiment 1, WRF models use three layers of two-way nested-grid ocean grid, ecto-entad plane grid size is respectively 13.5,
4.5,1.5km, innermost layer grid cover the complicated landform CFD model of embodiment 1, and WRF models divide 50 in the vertical direction
Layer, it is close under being dredged on grid, mesh refinement is carried out in 1km or less, divides 15 layers, Parameterization Scheme is set as：The simple ice of 3 classes of WSM
Scheme microphysical processes scheme, rrtm long-wave radiation schemes, Dudhia shortwave radiation schemes, Monin-Obukhov ground layers
Scheme, heat diffusion scheme, YSU PBL schemes and shallow convection Kain-Fritsch cumulus parameterization schemes, the time of integration
From 00 point to 12 points of October 22 in 2010, totally 12 hours, the output in each hour of WRF the model calculations was primary；
2.2) the wind field data for the different time that processing WRF models calculate, respectively obtain 4 around complicated landform CFD model
A series of coarse resolution wind field data (t, x, y, z, u, v, w) on a side interface and top surface on grid nodes, wherein t tables
Show the time, (x, y, z) indicates the three dimensional space coordinate of each grid node locations, three under (u, v, w) representation space rectangular coordinate system
The wind speed component in a direction, it is specific as follows：
For each side interface of complicated landform CFD model or top surface, the WRF model meshes lists that will intersect first with it
Member centre coordinate be projected on corresponding side interface or top surface, obtain each grid node locations three dimensional space coordinate (x,
Y, z), then corresponding position of each grid node in WRF models is found, pass through the WRF model wind field numbers around each grid node
According to the wind speed (u, v, w) of each grid node being obtained using distance-reverse weighting function interpolation, in this way, it is possible to obtain intricately
Around shape CFD model on 4 sides interface and top surface on a series of grid nodes coarse resolution wind field data (t, x, y, z, u,
V, w).
In embodiment 1, the wind field data that processing WRF models calculate obtain 4 lateral boundaries around complicated landform CFD model
The schematic diagrames of coarse resolution wind field data on face and top surface on grid nodes a series of is as shown in figure 3, Dark grey face in figure
In some boundary face of complicated landform CFD model, the centre coordinate of the WRF model meshes units intersected with it is projected in
On the face, the three dimensional space coordinate (x, y, z) for obtaining grid node locations passes through wind when using distance-reverse weighting function interpolation
Fast u1And u2Interpolation obtains the wind speed u of grid node, and by v1, v2, v3 and v4 interpolation obtains the wind speed v of grid node, passes through
w1, w2, w3And w4Interpolation obtains the wind speed w of grid node.
3, it realizes that WRF models and complicated landform CFD model are coupled in borderline data, obtains complicated landform CFD model 4
The schematic diagram of the fining grid wind field data changed over time on a side interface and top surface is as shown in Figure 4；
3.1) in OpenFOAM read boundary primary condition file, obtain complicated landform CFD model 4 sides interface and
The number of grid is refined on top surface；
3.2) for each side interface of complicated landform CFD model or top surface, first according to the fining net on the face
Lattice number determines the rectangular space coordinate of each fining grid, then finds what each fining grid node was obtained in step 2.2)
Corresponding position in coarse resolution wind field back end passes through the coarse resolution wind field data section around each fining grid node
Point, interpolation obtain the wind speed (u of each fining grid node*, v*, w*), (u*, v*, w*) lower three of representation space rectangular coordinate system
The wind speed component in direction, for side interface, interpolation method chosen distance inverse weight method, for top surface, interpolation method selection
The coplanar linear interpolation method of triangle obtains on the 4 sides interface and top surface of complicated landform CFD model at any time by this method
The fining grid wind field data of variation；
3.3) total flux of computational domain divided by the area of top surface are used, the wind speed w of top surface is as a result attached to1In, it obtains most
Whole top surface wind speed；
3.4) the fining grid wind field data of adjacent time are carried out according to the time interval of update boundary condition linear
Fining of each time grid wind field data that linear interpolation obtains are stored in different files by interpolation.
4, the complicated landform CFD model after coupling is calculated using OpenFOAM, to obtain more accurate complicated landform wind
Field characteristic analog result；
4.1) for the boundary condition of complicated landform CFD model, other than bottom surface is set as wall surface, 4 lateral boundaries of surrounding
Face and top surface are disposed as speed entrance；
4.2) first the fining grid wind field data of initial time are written in the boundary primary condition file of OpenFOAM,
And OpenFOAM is driven to calculate until wind field result is stablized, then using previous result of calculation as initial fields, according to the time
Interval after updating the boundary condition file of OpenFOAM successively and solving control file, then drives OpenFOAM to be calculated,
In, tri- kinds of large eddy simulation, k- ω or k- ε turbulence models may be used in when calculating, calculate and use PIMPLE algorithms, final to obtain more
Accurate complex topography wind field simulated behavior result.
Fig. 5 is to realize WRF and CFD coupled simulation complex topography wind fields in the embodiment of the present invention based on OpenFOAM, obtain
12 hours wind speed profile change schematic diagrams.
Coupled mode of the technical solution provided by the invention based on WRF and CFD simulates complex topography wind field, is conducive to be promoted
The precision of complicated landform wind-resources assessment and Numerical Simulation of Wind.
As it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, not to
The limitation present invention, all within the spirits and principles of the present invention made by all any modification, equivalent and improvement etc., should all include
Within protection scope of the present invention.
Claims (5)
1. a method of based on OpenFOAM realize WRF and CFD coupled simulation wind fields, which is characterized in that this method include with
Lower step：
1) complicated landform CFD model is established based on GIS data in OpenFOAM, generates fining grid；
2) mesoscale WRF modeling complex topography wind fields are established, the wind field data of WRF models are handled；
3) it realizes that WRF models and complicated landform CFD model are coupled in borderline data, obtains 4 sides of complicated landform CFD model
The fining grid wind field data changed over time in boundary face and top surface；
4) OpenFOAM is used to calculate the complicated landform CFD model after coupling, to obtain complex topography wind field simulated behavior knot
Fruit.
2. a kind of method for realizing WRF and CFD coupled simulation wind fields based on OpenFOAM according to claim 1, special
Sign is that step 1) specifically includes following sub-step：
1.1) center longitude for determining target obtains the ASTER of target region from GIS GIS-Geographic Information System
Terrain data under GDEM30 meters of resolution ratio geodetic coordinates handles GIS data, to obtain WGS_1984_UTM using ArcMap
Terrain data grid node under coordinate system；
1.2) parametrization establishes three-dimensional level terrain CFD model and generates fining grid：The computational domain of complicated landform CFD model
It is set as rectangular body region, complicated landform influential on computational domain is found according to actual topography and geomorphology, computational domain includes this
A little complicated landforms, and use unstructured prismatic mesh to be fully fitted complex geometry boundary, wherein computational domain not
With region using the different grid of resolution ratio, complicated landform CFD model use in the vertical direction adjacent mesh dimension ratio for
The σ grids of definite value；
1.3) realize that the data between the terrain data grid node of computational domain and three-dimensional level terrain CFD model grid node turn
It changes：According to the coordinate of each soffit grid back end in three-dimensional level terrain CFD model, it is found in terrain data grid
Then corresponding position in node uses the elevation of 3 terrain data grid nodes around the position, coplanar by triangle
Linear interpolation obtains the elevation of soffit grid node, meanwhile, in three-dimensional level terrain CFD model, its co-net on vertical direction
The elevation of lattice node can generate respective change with the variation of soffit grid node elevation, finally obtain the complexity after modification elevation
Landform CFD model；
1.4) OpenFOAM is imported by converter gambitToFoam, complicated landform CFD model is converted to OpenFOAM can
The format of identification.
3. a kind of method for realizing WRF and CFD coupled simulation wind fields based on OpenFOAM according to claim 1, special
Sign is that step 2) specifically includes following sub-step：
2.1) mesoscale WRF modeling complex topography wind fields are established, WRF models are cuboid, and horizontal plane direction uses multilayer
Nested grid, wherein the size of WRF model inner-layer thin fine grids and the size of outer layer coarse grids at odd number than relationship, with reality
Existing two-way nested-grid ocean；
2.2) the wind field data for the different time that processing WRF models calculate, respectively obtain 4 sides around complicated landform CFD model
A series of coarse resolution wind field data (t, x, y, z, u, v, w) in boundary face and top surface on grid nodes, wherein when t is indicated
Between, (x, y, z) indicates the three dimensional space coordinate of each grid node locations, lower three sides of (u, v, w) representation space rectangular coordinate system
To wind speed component, concrete processing procedure is as follows：
For each side interface of complicated landform CFD model or top surface, first by the WRF model meshes units intersected with it
Centre coordinate is projected on corresponding side interface or top surface, obtains the three dimensional space coordinate (x, y, z) of each grid node locations,
Corresponding position of each grid node in WRF models is found again, by the WRF model wind field data around each grid node, is made
The wind speed (u, v, w) of each grid node is obtained with distance-reverse weighting function interpolation, in this way, it is possible to obtain complicated landform CFD
A series of coarse resolution wind field data (t, x, y, z, u, v, w) around model on 4 sides interface and top surface on grid nodes.
4. a kind of method for realizing WRF and CFD coupled simulation wind fields based on OpenFOAM according to claim 3, special
Sign is that step 3) specifically includes following sub-step：
3.1) boundary primary condition file is read in OpenFOAM, obtains complicated landform CFD model 4 sides interface and top surface
The number of upper fining grid；
3.2) for each side interface of complicated landform CFD model or top surface, first according to the fining grid on the face
Number determine it is each fining grid rectangular space coordinate, then find each fining grid node obtained in step 2.2) it is thick
Corresponding position in resolution ratio wind field back end passes through the coarse resolution wind field data section around each fining grid node
Point obtains the wind speed (u of each fining grid node using distance-reverse weighting function interpolation*, v*, w*), wherein (u*, v*, w*) table
Show the wind speed component in lower three directions of rectangular coordinate system in space, by this method, obtains 4 lateral boundaries of complicated landform CFD model
The fining grid wind field data changed over time on face and top surface；
3.3) consider as a result flux conservation is attached to the wind speed w of top surface using the total flux of computational domain divided by the area of top surface1
In, obtain final top surface wind speed；
3.4) linear interpolation is carried out to the fining grid wind field data of adjacent time according to the time interval of update boundary condition,
Fining of each time grid wind field data that linear interpolation obtains are stored in different files.
5. a kind of method for realizing WRF and CFD coupled simulation wind fields based on OpenFOAM according to claim 1, special
Sign is that step 4) specifically includes following sub-step：
4.1) for the boundary condition of complicated landform CFD model, other than bottom surface is set as wall surface, surrounding 4 sides interface and
Top surface is disposed as speed entrance；
4.2) first the fining grid wind field data of initial time are written in the boundary primary condition file of OpenFOAM, and driven
Dynamic OpenFOAM is calculated until wind field result is stablized, then using previous result of calculation as initial fields, according to time interval,
After the boundary condition file of update OpenFOAM and solution control file successively, then OpenFOAM is driven to be calculated, wherein meter
Tri- kinds of large eddy simulation, k- ω or k- ε turbulence models may be used when calculation, calculate and use PIMPLE algorithms, final acquisition is more accurate
Complex topography wind field simulated behavior result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810330662.3A CN108563867A (en) | 2018-04-13 | 2018-04-13 | A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810330662.3A CN108563867A (en) | 2018-04-13 | 2018-04-13 | A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108563867A true CN108563867A (en) | 2018-09-21 |
Family
ID=63534954
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810330662.3A Pending CN108563867A (en) | 2018-04-13 | 2018-04-13 | A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108563867A (en) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109165476A (en) * | 2018-10-16 | 2019-01-08 | 广东电网有限责任公司 | A kind of modeling method and simulation of wind method of modularization wind-field model |
| CN109583096A (en) * | 2018-12-03 | 2019-04-05 | 华润电力技术研究院有限公司 | A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling |
| CN109614638A (en) * | 2018-10-30 | 2019-04-12 | 华中科技大学 | A kind of city wind environment CFD analogy method of indirect modeling |
| CN109635317A (en) * | 2018-10-30 | 2019-04-16 | 华中科技大学 | A kind of high altitude localities CFD artificially generated terrain topological method |
| CN109871568A (en) * | 2018-12-24 | 2019-06-11 | 中国电力科学研究院有限公司 | Consider the pipe gallery heat extraction Ventilating Design Method of wall surface heat transfer |
| CN110268864A (en) * | 2019-08-01 | 2019-09-24 | 南京林业大学 | A kind of berry picking method and power source rated output pressure value based on negative pressure absorbing principle determines method |
| CN110321642A (en) * | 2019-07-09 | 2019-10-11 | 中国气象局上海台风研究所 | A kind of typhoon simulation of wind method near the ground for improving kinematics effect |
| CN111339714A (en) * | 2020-02-17 | 2020-06-26 | 珠江水利委员会珠江水利科学研究院 | Multi-scale hydrodynamic coupling method based on FVCOM and OpenFOAM models |
| CN112149365A (en) * | 2020-09-29 | 2020-12-29 | 华能新能源股份有限公司 | Microscale wind model system and method |
| CN112434361A (en) * | 2020-11-16 | 2021-03-02 | 东南大学 | Intelligent ventilation system suitable for office buildings and data processing method |
| CN113095012A (en) * | 2021-05-17 | 2021-07-09 | 华能陇东能源有限责任公司 | Splicing and fusing method for numerical simulation calculation results of wind power plant flow field partitions |
| CN113160409A (en) * | 2021-02-26 | 2021-07-23 | 中国华能集团清洁能源技术研究院有限公司 | Complex terrain modeling method and system suitable for refined wind field simulation |
| CN113420458A (en) * | 2021-07-14 | 2021-09-21 | 广东电网有限责任公司广州供电局 | Wind field data refining processing method and device, terminal and medium |
| CN113435010A (en) * | 2021-05-26 | 2021-09-24 | 中国再保险(集团)股份有限公司 | Large-scale refined terrain digital simulation method and device |
| CN113868970A (en) * | 2021-08-13 | 2021-12-31 | 中国人民解放军国防科技大学 | Airport area multi-source wind field fusion method based on numerical simulation model and mesoscale meteorological model |
| CN113868971A (en) * | 2021-08-13 | 2021-12-31 | 中国人民解放军国防科技大学 | Airport area three-dimensional refined wind field reconstruction method based on numerical simulation model and historical wind field characteristics |
| CN114424268A (en) * | 2019-09-20 | 2022-04-29 | 亚洲航测株式会社 | Super-resolution stereoscopic vision processing system and program therefor |
| CN114492220A (en) * | 2021-11-23 | 2022-05-13 | 华北电力大学 | OpenFOAM platform-based method and system for predicting flow characteristics of ocean atmospheric boundary layer |
| CN114692355A (en) * | 2022-04-29 | 2022-07-01 | 交通运输部规划研究院 | Motor vehicle pollution diffusion simulation method and system by coupling WRF (write driver software) and OpenFOAM (open-world automatic learning machine) models |
| CN115408962A (en) * | 2022-11-02 | 2022-11-29 | 南京信息工程大学 | Wind field reconstruction method and system based on CFD simulation and wind lidar |
| CN115908739A (en) * | 2023-03-08 | 2023-04-04 | 成都流体动力创新中心 | Method and system for rapidly generating high-fidelity structural grid of complex terrain surface |
| CN116341410A (en) * | 2023-02-07 | 2023-06-27 | 成都流体动力创新中心 | Typical terrain wind field simulation method and system considering different stabilities |
| CN117197383A (en) * | 2023-11-03 | 2023-12-08 | 成都流体动力创新中心 | Terrain extension method, equipment and medium based on characteristic dimension of complex terrain |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106326625A (en) * | 2015-07-09 | 2017-01-11 | 中国电力科学研究院 | Method for simulating wind field by coupling WRF (weather research and forecasting) and OpenFOAM modes |
| CN107833282A (en) * | 2017-11-16 | 2018-03-23 | 广东电网有限责任公司电力科学研究院 | A kind of terrain modeling and mess generation method and device |
-
2018
- 2018-04-13 CN CN201810330662.3A patent/CN108563867A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106326625A (en) * | 2015-07-09 | 2017-01-11 | 中国电力科学研究院 | Method for simulating wind field by coupling WRF (weather research and forecasting) and OpenFOAM modes |
| CN107833282A (en) * | 2017-11-16 | 2018-03-23 | 广东电网有限责任公司电力科学研究院 | A kind of terrain modeling and mess generation method and device |
Non-Patent Citations (6)
| Title |
|---|
| ORKUN TEMEL ET AL.: "Coupled WRF-OpenFOAM study of wind flow over complex terrain", 《JOURNAL OF WIND ENGINEERING AND INDUSTRIAL AERODYNAMICS》 * |
| 刘炜: "大型客机座舱气流组织的模拟与逆向设计", 《中国博士学位论文全文数据库工程科技Ⅱ辑（月刊）》 * |
| 刘震卿: "桥梁墩台局部冲刷与CFD仿真研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑（月刊）》 * |
| 周鸥杰: "西部山区某悬索桥桥址处风环境实测及CFD数值模拟", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑（月刊）》 * |
| 杨剑: "WRF与CFD嵌套的局地台风风场数值模拟研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑（月刊）》 * |
| 詹永照等: "物体仿真中的纹理表示", 《计算机辅助设计与图形学学报》 * |
Cited By (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109165476A (en) * | 2018-10-16 | 2019-01-08 | 广东电网有限责任公司 | A kind of modeling method and simulation of wind method of modularization wind-field model |
| CN109614638A (en) * | 2018-10-30 | 2019-04-12 | 华中科技大学 | A kind of city wind environment CFD analogy method of indirect modeling |
| CN109635317A (en) * | 2018-10-30 | 2019-04-16 | 华中科技大学 | A kind of high altitude localities CFD artificially generated terrain topological method |
| CN109614638B (en) * | 2018-10-30 | 2023-04-07 | 华中科技大学 | Indirect modeling urban wind environment CFD simulation method |
| CN109635317B (en) * | 2018-10-30 | 2023-04-07 | 华中科技大学 | CFD (computational fluid dynamics) simulation terrain topology method for high-altitude area |
| CN109583096A (en) * | 2018-12-03 | 2019-04-05 | 华润电力技术研究院有限公司 | A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling |
| CN109871568A (en) * | 2018-12-24 | 2019-06-11 | 中国电力科学研究院有限公司 | Consider the pipe gallery heat extraction Ventilating Design Method of wall surface heat transfer |
| CN109871568B (en) * | 2018-12-24 | 2023-07-14 | 中国电力科学研究院有限公司 | Heat exhausting and ventilating design method for comprehensive pipe rack considering wall heat transfer |
| CN110321642A (en) * | 2019-07-09 | 2019-10-11 | 中国气象局上海台风研究所 | A kind of typhoon simulation of wind method near the ground for improving kinematics effect |
| CN110321642B (en) * | 2019-07-09 | 2022-04-26 | 中国气象局上海台风研究所 | Typhoon near-ground wind field simulation method for improving kinematic effect |
| CN110268864B (en) * | 2019-08-01 | 2021-07-06 | 南京林业大学 | Berry picking method based on negative pressure absorption principle and power source rated output pressure value determining method |
| CN110268864A (en) * | 2019-08-01 | 2019-09-24 | 南京林业大学 | A kind of berry picking method and power source rated output pressure value based on negative pressure absorbing principle determines method |
| CN114424268B (en) * | 2019-09-20 | 2023-09-29 | 亚洲航测株式会社 | Super-resolution stereoscopic vision processing system and storage medium |
| CN114424268A (en) * | 2019-09-20 | 2022-04-29 | 亚洲航测株式会社 | Super-resolution stereoscopic vision processing system and program therefor |
| CN111339714A (en) * | 2020-02-17 | 2020-06-26 | 珠江水利委员会珠江水利科学研究院 | Multi-scale hydrodynamic coupling method based on FVCOM and OpenFOAM models |
| CN112149365A (en) * | 2020-09-29 | 2020-12-29 | 华能新能源股份有限公司 | Microscale wind model system and method |
| CN112149365B (en) * | 2020-09-29 | 2023-06-30 | 华能新能源股份有限公司 | Micro-scale wind model system and method |
| CN112434361A (en) * | 2020-11-16 | 2021-03-02 | 东南大学 | Intelligent ventilation system suitable for office buildings and data processing method |
| CN113160409A (en) * | 2021-02-26 | 2021-07-23 | 中国华能集团清洁能源技术研究院有限公司 | Complex terrain modeling method and system suitable for refined wind field simulation |
| CN113095012A (en) * | 2021-05-17 | 2021-07-09 | 华能陇东能源有限责任公司 | Splicing and fusing method for numerical simulation calculation results of wind power plant flow field partitions |
| CN113435010A (en) * | 2021-05-26 | 2021-09-24 | 中国再保险(集团)股份有限公司 | Large-scale refined terrain digital simulation method and device |
| CN113435010B (en) * | 2021-05-26 | 2023-09-12 | 中国再保险(集团)股份有限公司 | Digital simulation method and device for large-scale fine terrain |
| CN113420458A (en) * | 2021-07-14 | 2021-09-21 | 广东电网有限责任公司广州供电局 | Wind field data refining processing method and device, terminal and medium |
| CN113868971A (en) * | 2021-08-13 | 2021-12-31 | 中国人民解放军国防科技大学 | Airport area three-dimensional refined wind field reconstruction method based on numerical simulation model and historical wind field characteristics |
| CN113868970A (en) * | 2021-08-13 | 2021-12-31 | 中国人民解放军国防科技大学 | Airport area multi-source wind field fusion method based on numerical simulation model and mesoscale meteorological model |
| CN114492220A (en) * | 2021-11-23 | 2022-05-13 | 华北电力大学 | OpenFOAM platform-based method and system for predicting flow characteristics of ocean atmospheric boundary layer |
| CN114692355B (en) * | 2022-04-29 | 2022-10-18 | 交通运输部规划研究院 | Motor vehicle pollution diffusion simulation method and system by coupling WRF (write driver software) and OpenFOAM (open-world automatic learning machine) models |
| CN114692355A (en) * | 2022-04-29 | 2022-07-01 | 交通运输部规划研究院 | Motor vehicle pollution diffusion simulation method and system by coupling WRF (write driver software) and OpenFOAM (open-world automatic learning machine) models |
| CN115408962A (en) * | 2022-11-02 | 2022-11-29 | 南京信息工程大学 | Wind field reconstruction method and system based on CFD simulation and wind lidar |
| CN116341410A (en) * | 2023-02-07 | 2023-06-27 | 成都流体动力创新中心 | Typical terrain wind field simulation method and system considering different stabilities |
| CN116341410B (en) * | 2023-02-07 | 2024-01-30 | 成都流体动力创新中心 | Typical terrain wind field simulation method and system considering different stabilities |
| CN115908739A (en) * | 2023-03-08 | 2023-04-04 | 成都流体动力创新中心 | Method and system for rapidly generating high-fidelity structural grid of complex terrain surface |
| CN115908739B (en) * | 2023-03-08 | 2023-05-05 | 成都流体动力创新中心 | Method and system for rapidly generating high-fidelity structural grid on surface of complex terrain |
| CN117197383A (en) * | 2023-11-03 | 2023-12-08 | 成都流体动力创新中心 | Terrain extension method, equipment and medium based on characteristic dimension of complex terrain |
| CN117197383B (en) * | 2023-11-03 | 2024-02-09 | 成都流体动力创新中心 | Terrain extension method, equipment and medium based on characteristic dimension of complex terrain |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108563867A (en) | A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM | |
| CN106326625B (en) | A kind of simulation of wind method coupling WRF and OpenFOAM mode | |
| CN112100720B (en) | Artificial intelligence city design form layout method for improving wind environment | |
| Liu et al. | LES study on the turbulent flow fields over complex terrain covered by vegetation canopy | |
| Piroozmand et al. | Coupled CFD framework with mesoscale urban climate model: Application to microscale urban flows with weak synoptic forcing | |
| CN108664705B (en) | OpenFOAM-based method for simulating surface roughness of complex terrain | |
| CN102722909B (en) | Assembly line topology network dynamic simulation method based on adaptive-dimensional DEM (dynamic effect model) | |
| CN108763825B (en) | Numerical simulation method for simulating wind field of complex terrain | |
| CN111768502A (en) | Non-structural grid two-dimensional flood simulation system based on GPU acceleration technology | |
| CN115994496B (en) | Urban park high resolution atmospheric CO 2 Numerical simulation method of concentration three-dimensional field | |
| Ferragut et al. | Comparison between 2.5-D and 3-D realistic models for wind field adjustment | |
| CN115329691B (en) | CFD and GIS-based ultra-large city wind environment simulation method | |
| CN114580310A (en) | Method for realizing scale reduction processing of WRF (hand-wrenching simulation) wind field based on PALM (PALM fiber laser) | |
| CN110335355A (en) | A kind of high automatic calculating method of large-scale shallow water lake | |
| CN114357571A (en) | Inversion method and system for atmospheric boundary layer wind field characteristics in constructed building environment | |
| CN109948214B (en) | Urban multi-scale wind environment numerical simulation method | |
| CN110175370A (en) | A kind of REGION OF WATER INJECTION OILFIELD recognition methods of city charge for remittance | |
| CN108376183A (en) | A kind of city CA model construction method based on principle of maximum entropy | |
| CN109583096A (en) | A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling | |
| CN115470731A (en) | Method and system for predicting wind power of wind field based on microclimate | |
| Liu et al. | Large-eddy simulations of the flow over an isolated three-dimensional hill | |
| CN109635317B (en) | CFD (computational fluid dynamics) simulation terrain topology method for high-altitude area | |
| CN104299262A (en) | Three-dimensional cloud simulating method based on speed field flow line | |
| CN111291520A (en) | Intelligent ocean anchoring buoy station location site selection method for optimizing space efficiency | |
| Houda et al. | A CFD model for simulating urban flow in complex morphological street network |
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: 20180921 |
|
| RJ01 | Rejection of invention patent application after publication |