CN109948214A - The multiple dimensioned wind environment method for numerical simulation in city - Google Patents
The multiple dimensioned wind environment method for numerical simulation in city Download PDFInfo
- Publication number
- CN109948214A CN109948214A CN201910179920.7A CN201910179920A CN109948214A CN 109948214 A CN109948214 A CN 109948214A CN 201910179920 A CN201910179920 A CN 201910179920A CN 109948214 A CN109948214 A CN 109948214A
- Authority
- CN
- China
- Prior art keywords
- city
- model
- numerical simulation
- boundary
- wide area
- 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
-
- 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
Abstract
Present disclose provides a kind of multiple dimensioned wind environment method for numerical simulation in city, comprising: generates the circular City Terrain model centered on city, and is inlayed by model by architecture ensemble model insertion City Terrain model, establishes wide area city overall model；Wide area city overall model is extracted using Boolean subtraction calculation, wide area city model is divided by round outside but spuare inside form and is simulated；Local block model is inlayed in extraction, and carries out interpolation to boundary, is provided accurate wind environment boundary condition for local block model and is simulated.The disclosure can avoid providing reasonable flow velocity boundary it is assumed that directly after the macroscopic view simulation of wind of city using the unreasonable boundary that becomes a mandarin for local block model using nested model method in the city of large scale is integrally simulated；Simultaneously as using city overall model uses circular boundary regime, same nested grid, the workload that greatly mitigation models and grid repeatedly divides can be used to different wind speed directions.
Description
Technical field
This disclosure relates to the multiple dimensioned wind environment numerical simulation side in urban wind environmental analysis field more particularly to a kind of city
Method.
Background technique
Currently, China is in the boost phase of Development of Urbanization, and since populous, soil is rare, per capita resources phase
To scarcity, urban development is caused to be had to using the high-intensitive mode developed.Under the influence of this development mode, urban heat island effect
Should constantly it enhance, atmosphere pollution is serious, and haze weather takes place frequently.The reason of leading to air quality severe exacerbation, is big in addition to plant emissions
It is outer caused by gas pollution sources, due also to it causes pollutant difficult caused by city layout the problem of the wind environments such as quiet wind and vortex
In diffusion, promote the formation of high density pollution air, affects the air quality in city some areas.Therefore, by controlling city
City's air-flow improves city drafting efficiency, optimizes moist heat, it has also become the important means solved these problems.
Three kinds of main methods of present analysis city wind environment are field measurement, wind tunnel experiment and computer numerical mould respectively
It is quasi-.Field survey method is highly prone to the limitation of long-term observation data acquisition, in the urban planning to face the future, is difficult to apply
Obtain the Wind Data of actual measurement.Wind tunnel experiment shrinkage scale model research test method is with modelling is at high cost, the period is long, is difficult to
The features such as studying different schemes simultaneously, time and effort consuming are generally only used in and study special building concentration or Minor Construction group
In the cities physical parameter researchs such as wind pressure, wind speed and air velocity distribution variation.
It is often that the groups of building for studying area are single from the surrounding enviroment of reality in practical study in order to adapt to this limitation
It solely separates, drops it off in the computational domain of numerical simulation, the initial data such as wind direction, wind speed that the meteorological observatory of use provides are simultaneously
Carry out analog study.Even the stronger work station of computing capability, in a wide range of wind environment research of City-level, in order to subtract
The landform of periphery or the building on periphery are often removed, city scope to be studied are individually put into calculating by the workload calculated less
Domain range is studied.
It there is problems in the prior art based on above-mentioned: first, the simulation of local environment and the perimeter strip inputted
Part and actual wind field environment be not corresponding.Under the influence of being in neighboring buildings environment due to the local environment of block, such as directly adopt
The original wind field numerical value conditions setting of meteorological observatory, will not meet practical wind field due to boundary condition, then tie simulation
The not accurate enough defect of fruit.Second, it is huge to will lead to calculation amount, generally as simulated using the method for urban global refined net
General purpose computer even work station is all difficult to meet huge calculation amount requirement.Third, study it is windy under the conditions of, simulation
Model must repeat to build, modeling work amount is significantly greatly increased, and then limit the multi-scheme digital-to-analogue optimization to many practical problems
Research.
Summary of the invention
(1) technical problems to be solved
It is set forth above at least partly to solve present disclose provides a kind of multiple dimensioned wind environment method for numerical simulation in city
The technical issues of.
(2) technical solution
According to one aspect of the disclosure, a kind of multiple dimensioned wind environment method for numerical simulation in city is provided, comprising: step
S100: the circular City Terrain model centered on city is generated, and is inlayed by model by architecture ensemble model insertion
In City Terrain model, wide area city overall model to be processed is established；Step S200: fortune is subtracted by boolean from step S100
It calculates and extracts wide area city overall model, wide area city model is divided by round outside but spuare inside form and simulated；Step S300:
It is extracted from step S200 and inlays local block model, and interpolation is carried out to boundary, provide accurate wind for local block model
Ambient boundary condition is simultaneously simulated.
In some embodiments of the present disclosure, step S100 includes: step S110: generating the circle centered on target cities
The City Terrain model of shape；Step S120: being modeled for architecture ensemble, generates urban architecture group model；Step S130: will
The urban architecture group model that the circular City Terrain model that step S110 is generated is generated with step S120 carries out Boolean calculation and adds
Operation obtains wide area city overall model.
In some embodiments of the present disclosure, step S200 further include: step S210: establishing height is 5-6H rotary table,
In, H is urban architecture maximum height；Step S220: it carries out Boolean calculation to subtract operation including: to be subtracted with the rotary table that step S210 is established
Remove City Terrain model and urban architecture group model；Step S230: rectangular area delimited, city scope is confined；Step S240: will
The result that step S220 and step S230 are obtained carries out Boolean calculation, by whole with the wide area city of round outside but spuare inside form segmentation
Body Model；Step S250: the wide area city overall model after step S240 segmentation is simulated.
In some embodiments of the present disclosure, step S300 includes: step S310: defining the calculating side of local block model
Boundary；Step S320: the interpolation in the overall model of city is calculated, and interpolation is mounted in the model boundary of local block；Step
S330: using the good boundary condition of interpolation, simulated conditions are set, local area block model carries out numerical simulation.
In some embodiments of the present disclosure, step S200 further include: step S260: the wide area city that step S250 is obtained
City's overall model carries out grid dividing.
In some embodiments of the present disclosure, step S300 further include: step S340: the local street that step S330 is obtained
Section model carries out grid dividing.
The circular City Terrain centered on target cities is generated in some embodiments of the present disclosure, in step S110
Model, further includes: step S111: the smallest circumference that target cities can be included is established；Step S112: by step S111
The area of obtained smallest circumference is reciprocal multiplied by maximum obstruction rate, obtains the periphery concentric with the smallest circumference in step S111
Product, as City Terrain model Calculation Plane range；Wherein, maximum obstruction rate is less than 5%.
In some embodiments of the present disclosure, step S250 further include: step S260: the wide area city that step S250 is obtained
City's overall model is arranged using common flow rate boundary.
In some embodiments of the present disclosure, the interpolation in step S320 is one or more in velocity field, pressure field and K E
It is a.
In some embodiments of the present disclosure, in step S250 or step S330, using K- ε two-equation model, initialization
The field flow of computational domain, and carry out numerical simulation calculation；K- ε two-equation model includes:
Continuity equation
The equation of momentum
K equation
E equation
Wherein,uiFor the fast component under cartesian coordinate system；xiFor cartesian coordinate system
Under coordinate components；T is time term；Re is reynolds number Re=UD/ ν；U a and D are the characteristic length of speed and simulation；P is pressure
Power item；giComponent for acceleration of gravity in the direction i；ρ is density；ν molecular viscosity coefficient is (μ=ρ ν)；νt(or μt=ρ νt)
For the turbulent flow viscosity of Turbulent Model；K is turbulent fluctuation kinetic energy；ε is turbulent fluctuation kinetic energy consumption rate；C1ε=1.44, C2ε=
1.92,Cμ=0.09, σk=1.0, Cμ=0.09 and σε=1.3 be model constants.
(3) beneficial effect
It can be seen from the above technical proposal that the multiple dimensioned wind environment method for numerical simulation in disclosure city at least have it is following
One of beneficial effect or in which a part:
(1) general boundary condition is cooperated to be arranged using round computational domain, using the city integral net disposably generated
Lattice, as long as changing the component of the wind speed to become a mandarin, do not need to draw again under the premise of only change becomes a mandarin wind speed size and direction
Subnetting lattice can greatly improve multi-state, the model Efficiency of multi-angle wind speed.
(2) method with interpolation boundary is inlayed using model, avoided in traditional analog method since the boundary that becomes a mandarin is by week
Side groups of building influence, and are unsatisfactory for index or logarithm vertical velocity profile problem.
(3) disclosure directly obtains interpolating on sides from the overall model of city, provides standard for one or more Local Model
True boundary condition goes out to flow than conventional method more scientific using symmetrical border hypothesis using the hypothesis of no gradient or lateral boundaries.
(4) analogy method of the disclosure stage by stage, under conditions of avoiding taking peripheral environment, the part that need to accurately simulate
Environment can individually delimit computational domain, and range size is fully controllable, thus make not limited by grid sum when partial simulation,
The details of building or block is looked after when grid dividing as far as possible, and then the accuracy of simulation is greatly improved, to effectively solve
Overall model modeling in city is complicated in conventional method, low efficiency, and since boundary condition setting can not when Local Model is simulated
Consider influence of the neighboring buildings to wind field, and then the problem of not being inconsistent with practical wind field occurs in analog result.
Detailed description of the invention
Fig. 1 is the schematic diagram of the multiple dimensioned wind environment method for numerical simulation in embodiment of the present disclosure city.
Fig. 2 is the multiple dimensioned wind environment method for numerical simulation flow diagram in embodiment of the present disclosure certain city of northern China.
Fig. 3 is the schematic diagram that the gradient wind formed in Different Ground roughness becomes.
Fig. 4 is certain city of northern China groups of building three-dimensional modeling scene figure in the embodiment of the present disclosure.
Fig. 5 is that certain city of northern China overall model computational domain divides schematic diagram in the embodiment of the present disclosure.
Fig. 6 a is certain city of northern China overall model grid dividing overall effect figure in the embodiment of the present disclosure.
Fig. 6 b is certain city of northern China overall model grid dividing part detailed map in the embodiment of the present disclosure.
Fig. 7 a is that wind speed field plan view is integrally simulated in certain city of northern China in the embodiment of the present disclosure.
Fig. 7 b is corresponding high levels of three-dimensional distribution map in Fig. 7 a.
Fig. 8 a is certain city of northern China entirety simulated pressure field plan view in the embodiment of the present disclosure.
Fig. 8 b is corresponding high levels of three-dimensional distribution map in Fig. 8 a.
Fig. 9 is embodiment of the present disclosure local block model location mark figure in the overall model of wide area city.
Figure 10 is that the Local grid of embodiment of the present disclosure local block model divides schematic diagram.
Figure 11 is near-earth/nearly building surface wind speed field schematic diagram of local block model in the embodiment of the present disclosure.
Figure 12 is local block model in the embodiment of the present disclosure in 5 meters of elevation horizontal cross-section wind speed field schematic diagrames.
Figure 13 is local block model in the embodiment of the present disclosure in 5 meters of elevation horizontal cross-section pressure field schematic diagrames.
Specific embodiment
Present disclose provides a kind of multiple dimensioned wind environment method for numerical simulation in city, comprising: step S100: generating with city
Centered on circular City Terrain model, and inlayed by model by architecture ensemble model insertion City Terrain model,
Establish wide area city overall model；Step S200: extracting wide area city overall model by Boolean subtraction calculation from step S100,
Divide wide area city model by round outside but spuare inside form and is simulated；Step S300: it is extracted from step S200 and inlays office
Domain block model, and interpolation is carried out to boundary, accurate wind environment boundary condition is provided for local block model and is simulated.
The disclosure can be avoided in the city of large scale is integrally simulated using the unreasonable boundary vacation that becomes a mandarin using nested model method
It is fixed, directly after the macroscopic view simulation of wind of city, reasonable flow velocity boundary is provided for local block model；Simultaneously as using
City overall model uses circular boundary regime, and same nested grid can be used to different wind speed directions, greatly mitigates modeling
And the workload that grid repeatedly divides.
It needs to illustrate in advance, numerical simulation is to utilize Fluid Mechanics Computation theory CFD software (such as ANSYS-
FLUEUNT, STAR-CCM and FLOW3D etc.) carry out wind environment simulation method, it have fast and convenient, accurate and effective, to people
Power material resources expend the advantages that relatively small, visual result and are widely used.CFD method for numerical simulation can be by urban morphology
Wind environment simulation, dissect urban morphology and wind environment coupled relation, the overall process of visualization, image conversion is provided, be system,
It the understanding of science and understands city wind environment rule in depth, provides strong support.
But since CFD is in city or building simulation process, computational domain need to expand the 10- of be selected in institute's survey region
15 times of ranges are limited by general Computing ability, and current CFD numerical simulation software is generally only applicable to Small and Medium Sized
(several hectares are arrived tens of sq-kms).In order to adapt to this limitation, in practical study, the groups of building in area will be often studied from reality
Surrounding enviroment in be separately separated out, drop it off in the computational domain of numerical simulation, and use meteorological observatory provide wind direction,
The initial data such as wind speed simultaneously carry out analog study.Even the stronger work station of computing capability, in a wide range of wind of City-level
In Environmental Studies, in order to reduce the workload of calculating, often the landform of periphery or the building on periphery are removed, it individually will be to be studied
City scope is put into computational domain range and is studied.
For the purposes, technical schemes and advantages of the disclosure are more clearly understood, below in conjunction with specific embodiment, and reference
The disclosure is further described in attached drawing.
Disclosure some embodiments will be done referring to appended attached drawing in rear and more comprehensively describe to property, some of but not complete
The embodiment in portion will be shown.In fact, the various embodiments of the disclosure can be realized in many different forms, and should not be construed
To be limited to this several illustrated embodiment；Relatively, these embodiments are provided so that the disclosure meets applicable legal requirement.
In first exemplary embodiment of the disclosure, a kind of multiple dimensioned wind environment numerical simulation side in city is provided
Method.As shown in Figure 1, the multiple dimensioned wind environment method for numerical simulation in disclosure city includes:
Step S100: the circular City Terrain model centered on city is generated, and is inlayed by model and builds city
It builds in group model insertion City Terrain model, establishes wide area city overall model to be processed.
Step S200: wide area city overall model is extracted by Boolean subtraction calculation from step S100, by round outside but spuare inside
Form segmentation wide area city model and simulated.
Step S300: extracting from step S200 and inlay local block model, and carries out interpolation to boundary, is local block
Model provides accurate wind environment boundary condition and is simulated.
It is further described below in conjunction with specific implementation method of the Fig. 2 to the disclosure, the ANSYS- that the present embodiment uses
FLUENT software is simulated.It is since it can be with SKETCHUP or other CAD modeling softwares etc. using ANSYS-FLUENT
Seamless link.Slightly adjust other equally applicable CFD softwares of disclosed technique scheme.The disclosure can be using nested model method
During the city of large scale is integrally simulated, avoid assuming such as Fig. 3 using the unreasonable boundary that becomes a mandarin, directly in city macroscopic view wind field mould
After quasi-, reasonable flow velocity boundary is provided for local block model；Simultaneously as using city overall model uses circular side
Boundary's form, therefore same nested grid, the work that greatly mitigation models and grid repeatedly divides can be used to different wind speed directions
Amount.Below by taking the simulation of certain city overall model local environment as an example, it is specifically described:
Step S100 can be refined further
Step S110: the circular City Terrain model centered on target cities is generated.It is simulated and is calculated due to wind environment
The size in region is directly closely related with the authenticity of analog result, and the small flow field of simulated domain can be distorted, but excessively increases and calculate
The problems such as region causes grid number more, calculation amount and added cost.Therefore, round computational domain size is meeting maximum obstruction rate
The requirement of index less than 5%, and be minimized under this condition in the satisfaction.It should be noted that blocking rate refers to simulation wind direction
The ratio between the building projected area in direction and the free cross-section of computational domain.Using it is circular be one kind be different from traditional calculations domain divide
Method, since wind environment simulation will usually carry out multi-scheme multi-angle inflow angle wind and compare, common modeling and grid are raw
At considerable time and efforts can be occupied.Using round computational domain it is possible to prevente effectively from using traditional analog method because of incoming flow side
To difference, the drawbacks of different numerical models need to be built repeatedly, thus significant increase computational efficiency.It is mainly adopted in specific implementation
Three-dimensional modeling is carried out to certain City Terrain and groups of building with SKETCH software.According to Google-earth architectural plane position and E
The depth of building data of Urban Architecture group, as shown in Figure 4.
Step S120: being modeled for architecture ensemble, generates urban architecture group model.
Step S130: the architecture ensemble that the circular City Terrain model and step S120 that step S110 is generated generate
Model carries out Boolean calculation and adds operation, obtains wide area city overall model.
Specifically, step S200 includes:
Step S210: establishing height is 5-6H rotary table, wherein H is urban architecture maximum height.
Step S220: it carries out Boolean calculation to subtract operation including: that the step S210 rotary table established is subtracted City Terrain and is built
Build 3D solid, obtain include the 3D solid of city space wind environment simulation model.
Step S230: rectangular area delimited, city scope is confined.It should be noted that the rectangular extent actually confined can
It is bigger than city actual range by 10 or so.
Step S240: the result that step S220 and step S230 are obtained carries out Boolean calculation, by with round outside but spuare inside
The wide area city overall model of form segmentation obtains including the rectangle computational domain of the 3D solid of city space and the outer circle of periphery
The wind environment simulation model of interior side.
As shown in figure 5, grid dividing is carried out to the 3D solid computational domain of city space to after step S240, it is main to use
Maneuverable tetrahedral grid, and utilize the grid number of Boundary Layer Method control pedestrian level.In order to accurately simulate near-earth
The change of gradient of wind field is using body fitted anisotropic mesh technology close within the scope of about 2 meters of ground.Further, to step S240 net
Lattice divide, and mainly use body fitted anisotropic mesh partitioning technology, control the size magnification ratio of first grid height and adjacent mesh
(being not more than 1.2) ensure that wind speed prolongs vertical distributing and can be distributed by digit rate or index percent close to ground underlying surface in this way.Cause
This, the main purpose of this step is to obtain macroscopical flow field and provide boundary condition for local high-precision analog, such as Fig. 6 a, Fig. 6 b institute
Show.
Step S250: the wide area city overall model after step S240 segmentation is simulated.Specifically use K- ε both sides journey
Model initializes the field flow of computational domain, and carries out numerical simulation calculation；K- ε two-equation model includes:
Continuity equation
The equation of momentum
K equation
E equation
Wherein,uiFor the fast component under cartesian coordinate system；xiFor cartesian coordinate system
Under coordinate components；T is time term；Re is reynolds number Re=UD/ ν；U a and D are the characteristic length of speed and simulation；P is pressure
Power item；giComponent for acceleration of gravity in the direction i；ρ is density；ν molecular viscosity coefficient is (μ=ρ ν)；νt(or μt=ρ νt)
For the turbulent flow viscosity of Turbulent Model；K is turbulent fluctuation kinetic energy；ε is turbulent fluctuation kinetic energy consumption rate；C1ε=1.44, C2ε=
1.92,Cμ=0.09, σk=1.0, Cμ=0.09 and σε=1.3 be model constants.
Step S300: it is extracted from the 3D solid space of the city space step S200 and inlays local block model (local
Block model can also be obtained by individually carrying out Geometric Modeling in the range of the overall model of wide area city), as shown in figure 9,
And interpolation is carried out to boundary, accurate wind environment boundary condition is provided for local block model and is simulated.Further include:
Step S310: the calculating boundary of local block model is defined.It should be noted that Ying Jinliang when design calculates boundary
The boundary name easily identified using wide area city overall model difference.Such as: with all directions boundary or corresponding block name.
Step S320: calculating velocity field in the three-dimensional in the overall model of city, the interpolation such as pressure field and K E, and by interpolation
It is mounted in the model boundary of local block.
Step S330: opening ANSYS-FLUENT calls in local block model meshes and utilizes the good boundary condition of interpolation,
And other simulated conditions are set, local area block model carries out numerical simulation.Physical simulation method can refer to step S250
The K- ε two-equation model of middle offer, details are not described herein.Analog result, will be preliminary as shown in Fig. 7 a, Fig. 7 b, Fig. 8 a and Fig. 8 b
Analog result is exported in the quasi- external file accurately simulated with FLUENT interpolative mode.
Step S340: grid dividing is carried out to the local block model that step S330 is obtained.
For identified goal in research targetedly subdivided meshes, as shown in Figure 10.For example pedestrian is carried out to wind
Boundary layer can be added in building surface and earth's surface by comfort level research.Influence to contour of building to vortex can then be built
Beyond the region of objective existence profile subrange grid is built to be encrypted.
Define the calculating boundary of local block model, the side that when setting uses city overall model interpolation to easily identify as far as possible
Boundary's name.
The simulation three-dimensional result that is ready for for calling in city overall model be interpolated into and calls in local block model side
In boundary.
Other simulated conditions are set, local area block model carries out numerical simulation, and local area block model result carries out
Reprocessing analysis and report.As shown in figures 11-13.
So far, attached drawing is had been combined the embodiment of the present disclosure is described in detail.It should be noted that in attached drawing or saying
In bright book text, the implementation for not being painted or describing is form known to a person of ordinary skill in the art in technical field, and
It is not described in detail.In addition, the above-mentioned definition to each element and method be not limited in mentioning in embodiment it is various specific
Structure, shape or mode, those of ordinary skill in the art simply can be changed or be replaced to it.
According to above description, those skilled in the art should have the multiple dimensioned wind environment method for numerical simulation in disclosure city
Clear understanding.
In conclusion the disclosure provides a kind of multiple dimensioned wind environment method for numerical simulation in city, using nested model method
It can avoid using the unreasonable boundary that becomes a mandarin in the city of large scale is integrally simulated it is assumed that directly in city macroscopic view wind field mould
After quasi-, reasonable flow velocity boundary is provided for local block model；Simultaneously as using city overall model uses circular side
Same nested grid, the workload that greatly mitigation models and grid repeatedly divides can be used to different wind speed directions in boundary's form.
It anticipates unless there are specified to be mutually on the contrary, the numerical parameter in this specification and appended claims is approximation, energy
Enough bases pass through the resulting required characteristic changing of content of this disclosure.Specifically, all be used in specification and claim
The middle content for indicating composition, the number of reaction condition etc., it is thus understood that repaired by the term of " about " in all situations
Decorations.Under normal circumstances, the meaning expressed refers to include by specific quantity ± 10% variation in some embodiments, some
± 5% variation in embodiment, ± 1% variation in some embodiments, in some embodiments ± 0.5% variation.
Furthermore word "comprising" does not exclude the presence of element or step not listed in the claims.It is located in front of the element
Word "a" or "an" does not exclude the presence of multiple such elements.
In addition, unless specifically described or the step of must sequentially occur, there is no restriction in the above institute for the sequence of above-mentioned steps
Column, and can change or rearrange according to required design.And above-described embodiment can be based on the considerations of design and reliability, that
This mix and match is used using or with other embodiments mix and match, i.e., the technical characteristic in different embodiments can be freely combined
Form more embodiments.
Algorithm and display are not inherently related to any particular computer, virtual system, or other device provided herein.
Various general-purpose systems can also be used together with enlightenment based on this.As described above, it constructs required by this kind of system
Structure be obvious.In addition, the disclosure is also not for any particular programming language.It should be understood that can use various
Programming language realizes content of this disclosure described herein, and the description done above to language-specific is to disclose this public affairs
The preferred forms opened.
The disclosure can by means of include several different elements hardware and by means of properly programmed computer come
It realizes.The various component embodiments of the disclosure can be implemented in hardware, or to run on one or more processors
Software module is realized, or is implemented in a combination thereof.It will be understood by those of skill in the art that can be used in practice micro-
Processor or digital signal processor (DSP) are some or all in the relevant device according to the embodiment of the present disclosure to realize
The some or all functions of component.The disclosure be also implemented as a part for executing method as described herein or
Whole device or device programs (for example, computer program and computer program product).Such journey for realizing the disclosure
Sequence can store on a computer-readable medium, or may be in the form of one or more signals.Such signal can
To download from internet website, perhaps it is provided on the carrier signal or is provided in any other form.
Similarly, it should be understood that in order to simplify the disclosure and help to understand one or more of each open aspect,
Above in the description of the exemplary embodiment of the disclosure, each feature of the disclosure is grouped together into single implementation sometimes
In example, figure or descriptions thereof.However, the disclosed method should not be interpreted as reflecting the following intention: i.e. required to protect
The disclosure of shield requires features more more than feature expressly recited in each claim.More precisely, as following
Claims reflect as, open aspect is all features less than single embodiment disclosed above.Therefore,
Thus the claims for following specific embodiment are expressly incorporated in the specific embodiment, wherein each claim itself
All as the separate embodiments of the disclosure.
Particular embodiments described above has carried out further in detail the purpose of the disclosure, technical scheme and beneficial effects
Describe in detail it is bright, it is all it should be understood that be not limited to the disclosure the foregoing is merely the specific embodiment of the disclosure
Within the spirit and principle of the disclosure, any modification, equivalent substitution, improvement and etc. done should be included in the guarantor of the disclosure
Within the scope of shield.
Claims (10)
1. a kind of multiple dimensioned wind environment method for numerical simulation in city, comprising:
Step S100: the circular City Terrain model centered on city is generated, and is inlayed by model by architecture ensemble
In model insertion City Terrain model, wide area city overall model to be processed is established；
Step S200: wide area city overall model is extracted by Boolean subtraction calculation from step S100, passes through round outside but spuare inside shape
Formula segmentation wide area city model is simultaneously simulated；
Step S300: extracting from step S200 and inlay local block model, and carries out interpolation to boundary, is local block model
Accurate wind environment boundary condition is provided and is simulated.
2. the multiple dimensioned wind environment method for numerical simulation in city according to claim 1, wherein the step S100 includes:
Step S110: the circular City Terrain model centered on target cities is generated；
Step S120: being modeled for architecture ensemble, generates urban architecture group model；
Step S130: the urban architecture group model that the circular City Terrain model and step S120 that step S110 is generated generate
It carries out Boolean calculation and adds operation, obtain wide area city overall model.
3. the multiple dimensioned wind environment method for numerical simulation in city according to claim 1, wherein the step S200 is also wrapped
It includes:
Step S210: establishing height is 5-6H rotary table, wherein H is urban architecture maximum height；
Step S220: it carries out Boolean calculation to subtract operation including: to subtract City Terrain model and city with the rotary table that step S210 is established
Build group model in city；
Step S230: rectangular area delimited, city scope is confined；
Step S240: the result that step S220 and step S230 are obtained carries out Boolean calculation, by round outside but spuare inside form
The wide area city overall model of segmentation；
Step S250: the wide area city overall model after step S240 segmentation is simulated.
4. the multiple dimensioned wind environment method for numerical simulation in city according to claim 1, wherein the step S300 includes:
Step S310: the calculating boundary of local block model is defined；
Step S320: the interpolation in the overall model of city is calculated, and interpolation is mounted in the model boundary of local block；
Step S330: using the good boundary condition of interpolation, simulated conditions are set, local area block model carries out numerical simulation.
5. the multiple dimensioned wind environment method for numerical simulation in city according to claim 2, wherein the step S200 is also wrapped
It includes:
Step S260: grid dividing is carried out to the wide area city overall model that step S250 is obtained.
6. the multiple dimensioned wind environment method for numerical simulation in city according to claim 3, wherein the step S300 is also wrapped
It includes:
Step S340: grid dividing is carried out to the local block model that step S330 is obtained.
7. the multiple dimensioned wind environment method for numerical simulation in city according to claim 2, wherein generated in the step S110
Circular City Terrain model centered on target cities, further includes:
Step S111: the smallest circumference that target cities can be included is established；
Step S112: the area of the obtained smallest circumference of step S111 is reciprocal multiplied by maximum obstruction rate, it obtains and step S111
In the concentric circumferential area of smallest circumference, as City Terrain model Calculation Plane range；Wherein, the maximum obstruction rate is small
In 5%.
8. the multiple dimensioned wind environment method for numerical simulation in city according to claim 3, wherein the step S250 is also wrapped
It includes:
Step S260: the wide area city overall model that step S250 is obtained is arranged using common flow rate boundary.
9. the multiple dimensioned wind environment method for numerical simulation in city according to claim 4, wherein inserting in the step S320
Value is one or more in velocity field, pressure field and K E.
10. the multiple dimensioned wind environment method for numerical simulation in city according to claim 3 or 9, wherein the step S250 or
In the step S330, using K- ε two-equation model, the field flow of computational domain is initialized, and carry out numerical simulation calculation；The K-
ε two-equation model includes:
Continuity equation
The equation of momentum
K equation
E equation
Wherein,uiFor the fast component under cartesian coordinate system；xiFor under cartesian coordinate system
Coordinate components；T is time term；Re is reynolds number Re=UD/ ν；U a and D are the characteristic length of speed and simulation；P is pressure term；
giComponent for acceleration of gravity in the direction i；ρ is density；ν molecular viscosity coefficient is (μ=ρ ν)；νt(or μt=ρ νt) it is turbulent flow
The turbulent flow viscosity of model；K is turbulent fluctuation kinetic energy；ε is turbulent fluctuation kinetic energy consumption rate；C1ε=1.44, C2ε=1.92, Cμ
=0.09, σk=1.0, Cμ=0.09 and σε=1.3 be model constants.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910179920.7A CN109948214B (en) | 2019-03-11 | 2019-03-11 | Urban multi-scale wind environment numerical simulation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910179920.7A CN109948214B (en) | 2019-03-11 | 2019-03-11 | Urban multi-scale wind environment numerical simulation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109948214A true CN109948214A (en) | 2019-06-28 |
| CN109948214B CN109948214B (en) | 2023-01-31 |
Family
ID=67008685
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910179920.7A Active CN109948214B (en) | 2019-03-11 | 2019-03-11 | Urban multi-scale wind environment numerical simulation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109948214B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110531107A (en) * | 2019-09-09 | 2019-12-03 | 中国科学技术大学 | The monitoring support component and monitoring test equipment of large scale brattice wind field |
| CN111767597A (en) * | 2020-06-18 | 2020-10-13 | 软通动力信息技术有限公司 | City model verification method, device, equipment and storage medium |
| WO2022047959A1 (en) * | 2020-09-04 | 2022-03-10 | 东南大学 | Artificially intelligent urban design morphological layout method for enhancing wind environment |
| CN115329691A (en) * | 2022-08-01 | 2022-11-11 | 上海师范大学 | CFD and GIS-based super-large city wind environment simulation method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030117397A1 (en) * | 2001-12-21 | 2003-06-26 | Hubrecht Alain Yves Nestor | Systems and methods for generating virtual reality (VR) file(s) for complex virtual environments |
| US20160203245A1 (en) * | 2015-01-14 | 2016-07-14 | State Grid Corporation Of China | Method for simulating wind field of extreme arid region based on wrf |
| CN109033664A (en) * | 2018-08-07 | 2018-12-18 | 北京建筑大学 | Based on the considerations of the architectural wind environment appraisal procedure of CFD building body draining effect |
| CN109063361A (en) * | 2018-08-17 | 2018-12-21 | 温州建正节能科技有限公司 | A kind of calculation method that simulation solid particulate matter moves in urban environment |
-
2019
- 2019-03-11 CN CN201910179920.7A patent/CN109948214B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030117397A1 (en) * | 2001-12-21 | 2003-06-26 | Hubrecht Alain Yves Nestor | Systems and methods for generating virtual reality (VR) file(s) for complex virtual environments |
| US20160203245A1 (en) * | 2015-01-14 | 2016-07-14 | State Grid Corporation Of China | Method for simulating wind field of extreme arid region based on wrf |
| CN109033664A (en) * | 2018-08-07 | 2018-12-18 | 北京建筑大学 | Based on the considerations of the architectural wind environment appraisal procedure of CFD building body draining effect |
| CN109063361A (en) * | 2018-08-17 | 2018-12-21 | 温州建正节能科技有限公司 | A kind of calculation method that simulation solid particulate matter moves in urban environment |
Non-Patent Citations (6)
| Title |
|---|
| "《土木建筑与环境工程》2013年第35卷第1~6期总第170~176期总目次", 《土木建筑与环境工程》 * |
| 刘丽等: "基于CFD技术的河谷型城市风环境模拟――以兰州市城关区为例", 《中国沙漠》 * |
| 吴珍珍等: "基于CFD模拟技术的深圳市城市风环境分析", 《工程质量》 * |
| 孙磊等: "基于RANS模型的球型结构风荷载数值模拟研究", 《建筑科学》 * |
| 崔桂香等: "城市大气环境的大涡模拟研究进展", 《力学进展》 * |
| 杨俊宴: "全数字化城市设计的理论范式探索", 《国际城市规划》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110531107A (en) * | 2019-09-09 | 2019-12-03 | 中国科学技术大学 | The monitoring support component and monitoring test equipment of large scale brattice wind field |
| CN111767597A (en) * | 2020-06-18 | 2020-10-13 | 软通动力信息技术有限公司 | City model verification method, device, equipment and storage medium |
| CN111767597B (en) * | 2020-06-18 | 2023-10-24 | 北京软通智慧科技有限公司 | City model verification method, device, equipment and storage medium |
| WO2022047959A1 (en) * | 2020-09-04 | 2022-03-10 | 东南大学 | Artificially intelligent urban design morphological layout method for enhancing wind environment |
| CN115329691A (en) * | 2022-08-01 | 2022-11-11 | 上海师范大学 | CFD and GIS-based super-large city wind environment simulation method |
| CN115329691B (en) * | 2022-08-01 | 2024-03-12 | 上海师范大学 | CFD and GIS-based ultra-large city wind environment simulation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109948214B (en) | 2023-01-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109948214A (en) | The multiple dimensioned wind environment method for numerical simulation in city | |
| CN112100720B (en) | Artificial intelligence city design form layout method for improving wind environment | |
| CN107133427A (en) | A kind of construction method of the flood risk analysis model based on 2DGIS platforms | |
| Samsonov et al. | Object-oriented approach to urban canyon analysis and its applications in meteorological modeling | |
| Yue et al. | An adaptive method of high accuracy surface modeling and its application to simulating elevation surfaces | |
| Li | Global transport on a spherical multiple-cell grid | |
| CN115329691B (en) | CFD and GIS-based ultra-large city wind environment simulation method | |
| CN111767597A (en) | City model verification method, device, equipment and storage medium | |
| CN115994496B (en) | Urban park high resolution atmospheric CO 2 Numerical simulation method of concentration three-dimensional field | |
| Hussein et al. | Influences of wind flow over heritage sites: A case study of the wind environment over the Giza Plateau in Egypt | |
| CN114357571A (en) | Inversion method and system for atmospheric boundary layer wind field characteristics in constructed building environment | |
| CN108052704A (en) | Mesoscale photochemical pollution simulation and forecast algorithm with Grid Nesting function | |
| CN114580310A (en) | Method for realizing scale reduction processing of WRF (hand-wrenching simulation) wind field based on PALM (PALM fiber laser) | |
| García-Sánchez et al. | The impact of level of detail in 3D city models for CFD-based wind flow simulations | |
| Li et al. | Exploring urban space quantitative indicators associated with outdoor ventilation potential | |
| CN105893590A (en) | Automatic processing method for real-situation cases of DTA (Digital Terrain Analysis) modelling knowledge | |
| Camelli et al. | Generating seamless surfaces for transport and dispersion modeling in GIS | |
| CN110765694B (en) | Urban surface water flow numerical simulation method based on simplified shallow water equation set | |
| Sayehbani et al. | Numerical modeling of wave and current patterns of beris port in east of chabahar-iran | |
| Houda et al. | A CFD model for simulating urban flow in complex morphological street network | |
| Jędrzejewski et al. | The problem of airflow around building clusters in different configurations | |
| CN104199887B (en) | Analysis of water environment method and system based on GIS-Geographic Information System Yu Delft3D FLOW computing engines | |
| CN114861263A (en) | Method and system for finely simulating mountainous terrain wind speed field | |
| KR102537488B1 (en) | Wind speed analyzing system in urban area based on inflow wind speed and morphological parameters, and wind speed analyzing method using the same | |
| Kwon et al. | Improvement of building-construction algorithm for using GIS data and analysis of flow and dispersion around buildings |
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 |