CN109614638A - A kind of city wind environment CFD analogy method of indirect modeling - Google Patents

A kind of city wind environment CFD analogy method of indirect modeling Download PDF

Info

Publication number
CN109614638A
CN109614638A CN201811273522.3A CN201811273522A CN109614638A CN 109614638 A CN109614638 A CN 109614638A CN 201811273522 A CN201811273522 A CN 201811273522A CN 109614638 A CN109614638 A CN 109614638A
Authority
CN
China
Prior art keywords
city
wind
studied
building
wind environment
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
Application number
CN201811273522.3A
Other languages
Chinese (zh)
Other versions
CN109614638B (en
Inventor
刘震卿
张冲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201811273522.3A priority Critical patent/CN109614638B/en
Publication of CN109614638A publication Critical patent/CN109614638A/en
Application granted granted Critical
Publication of CN109614638B publication Critical patent/CN109614638B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention belongs to urban wind environmental areas, and disclose a kind of city wind environment CFD analogy method of indirect modeling.This method includes the following steps: that (a) chooses city to be studied and determine its size, obtains the building data in the city to be studied;The computational domain for constructing the city to be studied carries out grid dividing to computational domain;(b) grid of building position adds sources of resistance in computational domain, carries out computational fluid dynamics emulation to the computational domain after addition sources of resistance, the wind field distribution situation in the city is obtained with this, to realize the wind environment simulation in the city.Through the invention, workload is reduced, the simulation to different regions city wind environment is particularly suitable for, reduces repetitive operation, improves working efficiency.

Description

A kind of city wind environment CFD analogy method of indirect modeling
Technical field
The invention belongs to urban wind environmental areas, more particularly, to a kind of city wind environment CFD mould of indirect modeling Quasi- method.
Background technique
In recent years, urban architecture density continues to increase, construction style is different, space layout is changeable, has to city wind environment Larger impact.And the area that city is most intensive as population activity, wind environment is to human comfort, structure security implication weight Greatly, Gucheng city wind environment CFD analogue simulation very it is necessary to.
Wu Zhenzhen etc. carries out grid dividing by taking Shenzhen as an example, to urban architecture Direct Modeling and to model, is based on Fluent fluid calculation software carries out CFD simulation to city wind environment using k- ε model;Similar research is also: Liu Li Jun etc., Money justice etc. carries out wind environment research simulation by establishing actual cities model.When being simulated to city wind environment, if to city Build Direct Modeling and grid division, it is necessary first to establish urban architecture geometrical model, then carry out network topology, workload is huge Greatly, when and for different regions Urban Simulation, model and grid division need to be re-established, repetitive operation is big, inefficiency.
Summary of the invention
Aiming at the above defects or improvement requirements of the prior art, the present invention provides a kind of city vanes of indirect modeling Sources of resistance is added in border CFD analogy method, the position by first establishing the computational domain in city, and in computational domain where building, With this simcity building to the inhibition of wind field, wherein carry out net after avoiding the geometrical model for establishing City Building LF Topology reduces workload, the simulation to different regions city wind environment is particularly suitable for, without constantly re-establishing model Grid division reduces repetitive operation, improves working efficiency.
To achieve the above object, it is proposed, according to the invention, provide a kind of city CFD Simulation of Wind Environment side of indirect modeling Method, which is characterized in that this method includes the following steps:
(a) it chooses city to be studied and determines the size in the city to be studied, obtain the building in the city to be studied Data;The computational domain that the city to be studied is constructed according to the size in the city to be studied carries out grid to the computational domain and draws Point;
(b) grid of the building position described in the computational domain adds sources of resistance, is built with this simcity Object is built to the inhibition of wind field, the computational domain after addition sources of resistance is subjected to computational fluid dynamics emulation, is obtained with this The wind field distribution situation in the city is obtained, to realize the wind environment simulation in the city.
It is further preferred that the sources of resistance preferably uses following expression formula in step (b),
Wherein,It is the resistance source item in the direction i, i represents the direction x, y or z CdIt is resistance coefficient,It is wind speed size It is the speed in the direction i.
It is further preferred that the building data include building plane angle point data and building in step (a) Altitude information.
It is further preferred that the calculating domain model is cuboid in step (a).
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, then the present invention carries out grid dividing to computational domain, with tradition by establishing the model about city computational domain Urban wind environmental simulation need to establish urban architecture geometrical model then subdivided grid compare, the present invention does not need to establish Build geometrical model, it is only necessary to cuboid computational domain grid division, efficiency be greatly improved, avoid and directly build to urban architecture The lengthy and jumbled process of mould grid division, reduce workload, especially to multiple and different Urban Simulations when, due to building in each city Difference need to repeatedly establish building model and carry out grid dividing, and workload is huge, by using method provided by the invention, nothing Grid dividing need to be carried out to building, a kind of pair of more convenient and fast simulation means of city wind environment are provided;
2, the present invention is by there are the grid cells at situation of building to add resistance source item automatically to simulate in computational domain Urban architecture is to the inhibition of wind field, and convenient and efficient, more efficient to carry out analogue simulation to different cities wind environment, simulation is tied Structure is more in line with reality.
Detailed description of the invention
Fig. 1 is the flow chart of indirect modeling Simulation city wind environment constructed by preferred embodiment according to the invention;
Fig. 2 is middle urban architecture schematic diagram constructed by preferred embodiment according to the invention;
Fig. 3 (a) is the axonometric drawing schematic diagram of urban architecture Direct Modeling model meshes in the prior art;
Fig. 3 (b) is the sectional drawing schematic diagram of urban architecture Direct Modeling model meshes in the prior art;
Fig. 4 (a) is model meshes axonometric drawing schematic diagram constructed by preferred embodiment according to the invention;
Fig. 4 (b) is model meshes axonometric drawing schematic diagram constructed by preferred embodiment according to the invention;
Fig. 5 is the wind field streamline and speed cloud atlas of Direct Modeling calculated result in the prior art;
Fig. 6 is the wind field streamline and speed cloud atlas of calculated result constructed by preferred embodiment according to the invention.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention 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 Not constituting a conflict with each other can be combined with each other.
The city wind environment CFD simulated technological process of a kind of indirect modeling of the present invention, as shown in Figure 1, specific implementation Mode is as follows:
1, research urban area is chosen, by taking certain imaginary zonule city as an example, urban architecture distribution is as shown in Fig. 2, determine Urban area size obtains urban architecture data;
(1) research urban area size is chosen;
(2) it obtains selected urban area and builds data, which includes: architectural plane angle point data and building height Data.
2, calculating domain sizes of the setting comprising research urban area are established and calculate domain model, and carry out grid to computational domain It divides.
(1) according to target cities area size, building height, determine calculate domain sizes (the long 10km of computational domain in embodiment, Wide 5km, high 1.5km);
(2) according to domain sizes are calculated, computational domain geometrical model is established, carries out grid dividing to domain model is calculated.
3, it is based on fluid calculation software, by grid cell in program automatic searching computational domain, and to by building geometry The grid cell addition resistance source item of encirclement builds the inhibition to wind field with simcity, and calculates Urban Wind Field.
(1) resistance is added in computational domain according to City Terrain data dot file automatically at the grid cell of right urban architecture Source item, with simcity building to the inhibition of wind field, resistance source item is calculate by the following formula:
Wherein,It is the resistance source item in the direction i, i represents the direction x, y or z CdIt is resistance coefficient, simulant building is to wind field Obstruction take CdFor infinity,It is wind speed size,It is the speed in the direction i.
(2) Simulation of Wind Field simulation is carried out to target cities based on fluid calculation software, uses standard k- ε turbulent flow in embodiment Model selects stable state to solve.
In order to embody indirect modeling proposed by the invention city wind environment CFD analogy method validity, will be with Urban architecture Direct Modeling calculated result is compared, Fig. 3 (a) is the axis of urban architecture Direct Modeling model meshes in the prior art Mapping schematic diagram;Fig. 3 (b) is the sectional drawing schematic diagram of urban architecture Direct Modeling model meshes in the prior art;Such as Fig. 3 (a) (b) shown in, city model grid dividing is as shown in the figure, and Fig. 4 (a) is mould constructed by preferred embodiment according to the invention Type grid axonometric drawing schematic diagram, Fig. 4 (b) are the signals of model meshes axonometric drawing constructed by preferred embodiment according to the invention Figure;As shown in Fig. 4 (a) and (b), white portion is addition resistance source item in figure.Two kinds of models are imported into fluid calculation respectively Software (embodiment use Fluent software) carries out Simulation of Wind Field simulation, and when simulation uses the uniform wind speed that becomes a mandarin, size 10m/s.
Fig. 5 is the wind field streamline and speed cloud atlas of Direct Modeling calculated result in the prior art, and Fig. 6 is according to the invention The wind field streamline and speed cloud atlas of calculated result constructed by preferred embodiment, the result in Fig. 5 and Fig. 6 are coincide good, equal table When revealing wind by urban architecture due to inhibition and there is a phenomenon where stream, it was demonstrated that it is proposed by the present invention by adding automatically Resistance source item is added to consider validity of the urban architecture to the inhibition of wind field, this indirect modeling method simcity vane Border is reliable.
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, any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should all include Within protection scope of the present invention.

Claims (4)

1. a kind of city wind environment CFD analogy method of indirect modeling, which is characterized in that this method includes the following steps:
(a) it chooses city to be studied and determines the size in the city to be studied, obtain the building data in the city to be studied; The computational domain in the city to be studied is constructed according to the size in the city to be studied, and grid then is carried out to the computational domain and is drawn Point;
(b) sources of resistance is added to the grid of the building position in the computational domain, with this simcity building The computational domain after addition sources of resistance is carried out computational fluid dynamics emulation, is obtained with this by object to the inhibition of wind field The wind field streamline and speed cloud distribution situation in the city, to realize the wind environment simulation in the city.
2. a kind of city wind environment CFD analogy method of indirect modeling as described in claim 1, which is characterized in that in step Suddenly in (b), the sources of resistance preferably uses following expression formula,
Wherein,It is the resistance source item in the direction i, i represents the direction x, y or z, CdIt is resistance coefficient,It is wind speed size,It is The speed in the direction i.
3. a kind of city wind environment CFD analogy method of indirect modeling as claimed in claim 1 or 2, which is characterized in that In step (a), the building data include building plane angle point data and depth of building data.
4. a kind of city wind environment CFD analogy method of indirect modeling as described in any one of claims 1-3, feature exist In in step (a), the calculating domain model is cuboid.
CN201811273522.3A 2018-10-30 2018-10-30 Indirect modeling urban wind environment CFD simulation method Active CN109614638B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811273522.3A CN109614638B (en) 2018-10-30 2018-10-30 Indirect modeling urban wind environment CFD simulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811273522.3A CN109614638B (en) 2018-10-30 2018-10-30 Indirect modeling urban wind environment CFD simulation method

Publications (2)

Publication Number Publication Date
CN109614638A true CN109614638A (en) 2019-04-12
CN109614638B CN109614638B (en) 2023-04-07

Family

ID=66002575

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811273522.3A Active CN109614638B (en) 2018-10-30 2018-10-30 Indirect modeling urban wind environment CFD simulation method

Country Status (1)

Country Link
CN (1) CN109614638B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110532666A (en) * 2019-08-26 2019-12-03 中国石油大学(华东) Urban micro-climate evolution analysis method based on wind environment variation
CN111767597A (en) * 2020-06-18 2020-10-13 软通动力信息技术有限公司 City model verification method, device, equipment and storage medium
CN113343421A (en) * 2021-05-04 2021-09-03 东南大学 Urban building energy consumption simulation method and system coupled with urban vegetation form
CN114021225A (en) * 2021-10-20 2022-02-08 中邮通建设咨询有限公司 Method, device, equipment and medium for simulating wind environment of base station by Web end
CN115329691A (en) * 2022-08-01 2022-11-11 上海师范大学 CFD and GIS-based super-large city wind environment simulation method
CN116151141A (en) * 2022-12-12 2023-05-23 中广核风电有限公司 Urban wind environment CFD simulation area selection method and device

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN108563867A (en) * 2018-04-13 2018-09-21 华中科技大学 A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN108563867A (en) * 2018-04-13 2018-09-21 华中科技大学 A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110532666A (en) * 2019-08-26 2019-12-03 中国石油大学(华东) Urban micro-climate evolution analysis method based on wind environment variation
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
CN113343421A (en) * 2021-05-04 2021-09-03 东南大学 Urban building energy consumption simulation method and system coupled with urban vegetation form
CN113343421B (en) * 2021-05-04 2022-10-25 东南大学 Urban building energy consumption simulation method and system coupled with urban vegetation form
CN114021225A (en) * 2021-10-20 2022-02-08 中邮通建设咨询有限公司 Method, device, equipment and medium for simulating wind environment of base station by Web end
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
CN116151141A (en) * 2022-12-12 2023-05-23 中广核风电有限公司 Urban wind environment CFD simulation area selection method and device
CN116151141B (en) * 2022-12-12 2024-01-30 中广核风电有限公司 Urban wind environment CFD simulation area selection method and device

Also Published As

Publication number Publication date
CN109614638B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
CN109614638A (en) A kind of city wind environment CFD analogy method of indirect modeling
CN107862738B (en) One kind carrying out doors structure three-dimensional rebuilding method based on mobile laser measurement point cloud
CN105513133B (en) A kind of city wind environment numerical map makes and display methods
CN108563867A (en) A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM
CN111768502A (en) Non-structural grid two-dimensional flood simulation system based on GPU acceleration technology
CN104318622B (en) Triangular mesh modeling method of indoor scene inhomogeneous three dimension point cloud data
CN102509339B (en) Method for simplifying vertex clustering of three-dimensional models with texture constraint
CN105160707B (en) Threedimensional model quick visualization method based on viewpoint index
CN109033664A (en) Based on the considerations of the architectural wind environment appraisal procedure of CFD building body draining effect
CN108664705B (en) OpenFOAM-based method for simulating surface roughness of complex terrain
CN111767597A (en) City model verification method, device, equipment and storage medium
CN102903139B (en) Accelerated rendering method for contours
CN109882164A (en) A kind of large scale acidification analogy method of fractured carbonate rock oil reservoir
CN104360896A (en) Parallel fluid simulation acceleration method based on GPU (Graphics Processing Unit) cluster
CN103679803A (en) Method for dynamic modeling for sea ice in navigation simulator scene
CN116050689A (en) Intelligent searching method, system, terminal and medium for wide-area space railway line
CN114417467A (en) Method, device, medium and equipment for establishing ventilation corridor of urban dense area
CN105528423A (en) Self-adaptive co-location pattern obtaining method and apparatus considering spatial instance distance weight
CN109635317A (en) A kind of high altitude localities CFD artificially generated terrain topological method
CN110162903A (en) A kind of urban architecture windward side density calculating method and system parallel based on grid
CN106909730A (en) Three-dimensional model building emulation mode and system based on homotopy mapping algorithm
CN113012286B (en) Method for constructing road DEM based on dense point cloud data
Zhang et al. Conformal adaptive hexahedral-dominant mesh generation for CFD simulation in architectural design applications
Menzel Evolvable free-form deformation control volumes for evolutionary design optimization
KR20220095964A (en) Apparatus and method for generating learning data based on map

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