CN106372320A - Method for performing large eddy simulation on highway tunnel turbulence by using sub-filtering scale model - Google Patents

Method for performing large eddy simulation on highway tunnel turbulence by using sub-filtering scale model Download PDF

Info

Publication number
CN106372320A
CN106372320A CN201610789060.5A CN201610789060A CN106372320A CN 106372320 A CN106372320 A CN 106372320A CN 201610789060 A CN201610789060 A CN 201610789060A CN 106372320 A CN106372320 A CN 106372320A
Authority
CN
China
Prior art keywords
model
sub
eddy
simulation
flow
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
CN201610789060.5A
Other languages
Chinese (zh)
Other versions
CN106372320B (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.)
Taizhou Vocational and Technical College
Original Assignee
Taizhou Vocational and Technical College
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 Taizhou Vocational and Technical College filed Critical Taizhou Vocational and Technical College
Priority to CN201610789060.5A priority Critical patent/CN106372320B/en
Publication of CN106372320A publication Critical patent/CN106372320A/en
Application granted granted Critical
Publication of CN106372320B publication Critical patent/CN106372320B/en
Expired - Fee Related 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
    • 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

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 provides a method for performing large eddy simulation on highway tunnel turbulence by using a sub-filtering scale model. By building a correction function D, an eddy-viscosity coefficient v<Tau> near the railway tunnel wall surface is corrected; after the method for performing large eddy simulation on highway tunnel turbulence by using the sub-filtering scale model is used, the dissipation quantity of the model near the wall surface is greatly improved, so that the dissipation of the small-dimension eddy near the wall surface is balanced; therefore the result conformity of the improved average flow speed profile and the DNS (domain name system) is better.

Description

A kind of large eddy simulation is carried out to vcehicular tunnel turbulent flow using sub- filter scale model Method
Technical field
The present invention relates to vcehicular tunnel field, especially, specifically one kind adopts sub- filter scale model to highway The method that turbulent flow big whirlpool in tunnel is simulated.
Background technology
Turbulent flow, is a kind of flow phenomenon of irregular.Its scrambling not only shows that the physical quantitys such as speed, pressure exist When aerial random distribution, more show the not repeated of it.Simple turbulent flow can obtain it with the method for mathematical analyses Main flow behavior, but the unpredictable complicated turbulent flow of Mathematical Method.For complicated turbulent, real except carrying out physics Beyond testing, numerical simulation is the Main Means it studied and is predicted.
For in theory, with direct Numerical (dns), it is possible to obtain the full detail of field of turbulent flow.But in fact, The complicated turbulent that direct Numerical is applied to high reynolds number (re) calculates, and needs huge computing resource.In theory Estimate, the amount of calculation realizing direct Numerical is proportional to re3.
Reynolds average method (rans) is to be asked with numerical method under the boundary condition of given mean motion and initial condition Solution Reynolds equation.Equal physical quantity when being that Reynolds average method is given, its major advantage is that amount of calculation is little.But its accuracy is relatively Difference, especially the corner in windward side easily occur the excessive phenomenon of Turbulent Kinetic cite { jincheng2011, cuncfd }.
Large eddy simulation method (les) is between this two classes method above-mentioned.Large eddy simulation method will be whole by filtering operation Individual Turbulent Flow Field partition of the scale is to filter distinguishable yardstick and the not distinguishable yardstick of sub- filtering.For the fluid filtering distinguishable yardstick Motion, is calculated with direct Numerical (dns), and the fluid motion for sub- filter scale, then pass through the sub- filter scale of construction Model to be simulated, the closed mode in similar Reynolds average method.Little yardstick motion due to high reynolds number turbulent flow has , there is the possibility that construction does not rely on concrete flowing, pervasive sub- filter scale model in theory in universality, thus big whirlpool mould Intend being considered to have the ability improving turbulent flow precision of prediction, it will substitute commonly used thunder in current engineering in the near future Promise averaging method, becomes the main method of turbulent flow engineering calculation.
However, due to the understanding enough to turbulent flow little yardstick poverty of movement, existing all kinds of Asias filter scale model is also not Sub- filter scale (sfs) motion can be described exactly to filtering the impact that distinguishable yardstick (rfs) moves, thus drawing in the calculation Enter model error.On the other hand, because mostly the turbulent flow handled by large eddy simulation is the unsteady nonlinear many of complexity Nanoscale systems, directly apply some successful numerical methods in traditional permanent, laminar flow and turbulent flow Reynolds average calculate can carry Carry out larger numerical error.The non-linear dynamics process that numerical error and model error pass through complexity influences each other so that mesh Front large eddy simulation calculates and shows sizable uncertainty.
Content of the invention
The purpose of the present invention be calculate for current large eddy simulation vcehicular tunnel turbulent flow is shown sizable not true Qualitatively problem, proposes a kind of method carrying out large eddy simulation to vcehicular tunnel turbulent flow using sub- filter scale model.Based on Upper thought, wishes to propose new eddy viscosity computational methods, so that sfs model can cutting by force near wall herein Appropriate kinetic energy dissipation rate is provided in the calculating cutting turbulent flow.
The technical scheme is that
A kind of method using the big whirlpool simulation to vcehicular tunnel turbulent flow of sub- filter scale model, it includes following step Rapid:
S1, set up correction function d, aforementioned correction function d is used for the eddy viscosity υ to vcehicular tunnel near wallτCarry out Revise;
D=d1×d2
d1=1-exp [- (y+/a+)3]
d 2 = c 1 e x p &lsqb; ( y + / y 0 + ) 4 &rsqb; + 1
Wherein, a+=25, c1=200,y+Represent is the dimensionless vertical height of wall;
System is glued in s2, the whirlpool of the subfilter model being adopted according to following formula calculating vcehicular tunnel turbulent flow large eddy simulation Number:
υτ=(cd δ)2|s|
Wherein, c is constant, δ=(δ x δ y δ z)1/3, wherein δ x, δ y and δ z are in the x, y and z directions respectively Mesh width;S is local train rate, and s is defined asWherein,
The c of the present invention is 0.1.
Beneficial effects of the present invention:
The present invention using sub- filter scale model vcehicular tunnel turbulent flow is carried out large eddy simulation method application after, model Dissipation amount near wall is greatly improved, and so that the dissipation near wall small size whirlpool is balanced, average after then improving Flow velocity profile is met preferably with the result of dns.
Brief description
Fig. 1 is velocity flow profile and distribution of shear stress schematic diagram in the vcehicular tunnel of the present invention.
Fig. 2 is using one of mean flow generalized section during standard smagorinsky modeling vcehicular tunnel turbulent flow.
Fig. 3 is two using mean flow generalized section during standard smagorinsky modeling vcehicular tunnel turbulent flow.
Fig. 4 is the distribution schematic diagram of the correction function d that the present invention is given.
Fig. 5 is the function curve near wall for the correction function d that is given of the present invention.
Fig. 6 is the correction function d that is given of the present invention in the function curve remotely from wall.
Fig. 7 is the distribution curve of eddy stress.
1st, vcehicular tunnel top surface;2nd, Tunnel Pavement;3rd, the elevation of zero point in vcehicular tunnel flow field;4th, vcehicular tunnel flow field Middle part flex point height;5th, the distribution of time average velocity in vcehicular tunnel is illustrated;6th, the turbulent shear stress in vcehicular tunnel is illustrated.
Specific embodiment
The present invention is further illustrated with reference to the accompanying drawings and examples.
Existing Asia filter scale (sfs) model can be roughly classified into feature simulation according to its building method (functional modeling) and constructivity simulate (struetral modeling) two classes.
The model of feature simulation, with eddy viscosity models as representative, filters distinguishable speed by constructionFunction come straight Connect the impact of simulation sfs stress.This class model mainly considers the actual effect on the motion impact of rfs yardstick for the sfs yardstick motion, no It should be understood that the details of sfs yardstick motion.Its Typical Representative has: smagorinsky model, composes eddy viscosity models, structure function Model (strueture-function model).Kolmogorov eddy viscosity models of v.c.wong and d.k.lilly etc..Mesh The front dynamic smagorins material model obtaining immense success obviously falls within this class.
The model of constructivity simulation constructs part or all of sfs speed by certain approach from the distinguishable velocity field of filtering Field u', then byAnd the definition of sfs stress tensor directly obtains the model of sfs stress tensor.This class model is main There is chaotic maps model, sub- filter scale estimates model and Similarity Model etc..
Air flow in vcehicular tunnel is the complicated turbulent that vehicle travels and causes in the confined space.And due to motion , as the presence of flow field disturbance factor, tunnel flow field turbulence level is strong for vehicle.Under uniformly continuous wagon flow, in tunnel when equal Velocity flow profile assumes the s type distribution of double-log along tunnel height.In a small thickness on tunnel inner wall surface (about 20mm), turbulent shear stress is a steady state value.In other regions in tunnel flow field, turbulent shear stress is high with vcehicular tunnel Degree assumes parabola change.At the middle part (h=(h+hv)/2) in tunnel, shear stress reaches its minima.Tunnel flow field velocity The flex point of distribution also here height.The maximum of turbulent shear stress concentrates on top and the near-bottom region in tunnel.Highway Velocity flow profile in tunnel and distribution of shear stress are as shown in Figure 1.1st, vcehicular tunnel top surface;2nd, Tunnel Pavement;3rd, highway The elevation of zero point in tunnel flow field, can be taken as 0.05m;4th, flex point height in the middle part of vcehicular tunnel flow field, can be taken as 2.525m.
When using standard smagorinsky modeling vcehicular tunnel turbulent flow, due near wall Strong shear stress Exist so that calculated mean flow section is higher, result is not good enough.Horiutik large eddy simulation calculating in, using dynamic Eddy viscosity models (dynamic smagorinsky model) and dynamic mixed model (dynamic mixed model), calculate The mean flow section obtaining is substantially higher (Fig. 2).The calculating of gullbrand j also indicates that, dynamic eddy viscosity models and dynamically mixing Matched moulds type can cause the mean flow section of tunnel wall surface flow field substantially higher (Fig. 3).Also carried out with dynamic eddy viscosity models herein Calculate, the result obtaining also has a substantially higher problem (Fig. 2, Fig. 3).
It is substantially energy effect that sub- filter scale motion can divide the effect of pigtail yardstick motion to filtering, thus (sub- filtering Scale Model) as long as being capable of the energy transmission in this Liang Ge yardstick area of correct balance it is possible to fully describe sub- filter scale motion The impact of (pigtail yardstick can be divided to move filtering).The model of feature simulation is simulated by eddy viscosity from the angle of the energy balance Sfs dissipates, and accurate mean kinetic energy generally can be provided in the case of free turbulence to dissipate.However, due to eddy viscosity Dependency between model and sfs stress is very low, and this class model is difficult to provide accurate sfs stress, thus in sfs stress lifting The situation such as occasion to be acted on, such as turbulent closure scheme, wall flowing, result is just not ideal enough.This research is paid close attention to tunnel top and is opened Flow field characteristic at mouthful, so obtain relatively accurate flow field parameter at wall with regard to inevitable requirement.
Real sfs stress while filtering distinguishable yardstick extracts kinetic energy, also anti-to filtering distinguishable yardstick To transmission kinetic energy, i.e. so-called inverse level string.It is true that the research to turbulent flow direct Numerical result shows, although averagely dynamic Energy level string always positive (transmitting to little yardstick from large scale), the instantaneous kinetic energy transmission of each local space position then both may Be forward direction be also likely to be reverse, and the positive level string of energy and inverse level string in magnitude quite, both ratio mean kinetic energy levels String is much greater.In order to provide correct Turbulent Mean Kinetic Energy Dissipation Rate, the model of the so pure dissipativeness of eddy viscosity models, its magnitude Must match with Turbulent Mean Kinetic Energy Dissipation Rate, and Turbulent Mean Kinetic Energy Dissipation Rate is much smaller than the instantaneous energy transmission of partial points, this just determines Determine eddy viscosity models to be difficult to match with sfs stress in partial points.Dynamic smagorinsky model utilize local similarity Dynamic decision model coefficient, can obtain inverse level string effect by producing negative sfs viscosity, however such against can transmit be with very Real level string effect is inconsistent.These shortcomings of eddy viscosity models make it be difficult to become good " universality " model, It is true that the universal constant different pattern of flow not being all suitable for is a major defect of eddy viscosity models.
The model of constructivity simulation to obtain sfs stress by constructing sfs yardstick and moving, and some of them achieve huge Success.For example, Similarity Model, inverse convolution model and velocity gradient model etc., all show the height and true stress between Degree dependency, and inverse level string effect can be automatically generated.However, these models and some other explicitly consider energy The sfs model of amount balance, all shows and can not produce the problem that enough mean kinetic energies dissipate, and needs and eddy viscosity models cooperation Use, obtain so-called mixed model.
In the energy level string mechanism of turbulent flow, the momentum-transport of turbulent flow and Turbulent Mean Kinetic Energy Dissipation Rate are all by Large Scale Motion institute Determine, Large Scale Motion determines and extracts kinetic energy from mean flow and be transmitted to the speed of little yardstick motion, and the motion of little yardstick is then led to Cross and produce the small-scale structure accordingly enriching degree, with the kinetic energy dissipation providing the kinetic energy input with Large Scale Motion to match Rate.The fact that it is meant that Large Scale Motion is insensitive to little yardstick motion details, as long as thus large eddy simulation calculate sfs Model provides appropriate kinetic energy dissipation rate, and the result of calculation of Large Scale Motion just will have ballpark statistical property, and The details of sfs model dissipation mechanism then affects less on result.Based on above thought, wish that proposing a new whirlpool glues system herein Number calculating method, so that sfs model can provide appropriate kinetic energy dissipation in the calculating of the Strong shear turbulent flow of near wall Rate.
The most basic problem of large eddy simulation seeks to provide suitable subscale grid model.
The present invention proposes new sub- filter scale model, is allowed to be applied to the flow pattern of vcehicular tunnel.New sub- filter Ripple Scale Model is mainly reflected in the algorithms of different of the correction function d of near wall eddy viscosity.Correction factor in the present invention Comprise two parts, concrete form is as follows:
D=d1×d2
d1=1-exp [- (y+/a+)3]
d 2 = c 1 e x p &lsqb; ( y + / y 0 + ) 4 &rsqb; + 1
Wherein, a+=25, c1=200,
Fig. 4 is the distribution of correction function d given herein.At wall, correction function d is zero.Mistake in velocity profile Cross correction function d in layer to increase rapidly, whenWhen, new correction function d reaches its peak value, is gradually lowered afterwards again, and becomes To in 1.As seen from the figure, eddy viscosity is strengthened near the transition zone of velocity profile.Fig. 6 exists for new correction function d The function curve of near wall, Fig. 7 is new correction function d in the function curve remotely from wall.
Compliance test result and Test of accuracy
Fig. 7 gives the contrast of the mean flow velocity profile before and after improvement and dns result.Abscissa is with wall parameter no The normal direction coordinate of dimension, empty circles are the result of calculation before improving, and solid circles are the result of calculation after improving, wherein Rer=387, solid line directly simulates the result of calculating, rer=395 for moser r d.It can be seen that before improvement, Because the dissipation near wall small size whirlpool is larger, and the dissipation amount of model is also not enough to be allowed to balance, then mean flow speed Section is higher.After improvement, dissipation amount near wall for the model is greatly improved, and so that the dissipation near wall small size whirlpool is obtained Balance, after then improving, mean flow velocity profile and the result of dns meet preferably.
Part that the present invention does not relate to is all same as the prior art or can be realized using prior art.

Claims (2)

1. a kind of using the sub- filter scale model method that big whirlpool is simulated to vcehicular tunnel turbulent flow, it is characterized in that it include with Lower step:
S1, set up correction function d, aforementioned correction function d is used for the eddy viscosity υ to vcehicular tunnel near wallτIt is modified;
D=d1×d2
d1=1-exp [- (y+/a+)3]
d 2 = c 1 e x p &lsqb; ( y + / y 0 + ) 4 &rsqb; + 1
Wherein, a+=25, c1=200,y+Represent is the dimensionless vertical height of wall;
S2, the eddy viscosity of the subfilter model being adopted according to following formula calculating vcehicular tunnel turbulent flow large eddy simulation:
υτ=(cd δ)2|s|
Wherein, c is constant, δ=(δ x δ y δ z)1/3, wherein δ x, δ y and δ z are grid in the x, y and z directions respectively Width;S is local train rate, and s is defined asWherein,
2. the method using the big whirlpool simulation to vcehicular tunnel turbulent flow of sub- filter scale model according to claim 1, It is characterized in that described c is 0.1.
CN201610789060.5A 2016-08-31 2016-08-31 Method for performing large vortex simulation on highway tunnel turbulence by adopting sub-filtering scale model Expired - Fee Related CN106372320B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610789060.5A CN106372320B (en) 2016-08-31 2016-08-31 Method for performing large vortex simulation on highway tunnel turbulence by adopting sub-filtering scale model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610789060.5A CN106372320B (en) 2016-08-31 2016-08-31 Method for performing large vortex simulation on highway tunnel turbulence by adopting sub-filtering scale model

Publications (2)

Publication Number Publication Date
CN106372320A true CN106372320A (en) 2017-02-01
CN106372320B CN106372320B (en) 2020-03-20

Family

ID=57899838

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610789060.5A Expired - Fee Related CN106372320B (en) 2016-08-31 2016-08-31 Method for performing large vortex simulation on highway tunnel turbulence by adopting sub-filtering scale model

Country Status (1)

Country Link
CN (1) CN106372320B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108133079A (en) * 2017-12-01 2018-06-08 上海理工大学 Field of turbulent flow coherent structure extracting method in IC engine cylinder
CN109858148A (en) * 2019-01-30 2019-06-07 南京航空航天大学 A kind of turbulent flow calculation method based on part filtering
CN112131800A (en) * 2020-07-20 2020-12-25 中国科学院力学研究所 Novel large vortex simulation method and device based on energy flow similarity
CN113111610A (en) * 2021-05-10 2021-07-13 中国空气动力研究与发展中心计算空气动力研究所 Novel sub-lattice scale model establishing method
CN115859760A (en) * 2022-04-29 2023-03-28 中广核风电有限公司 Method and device for simulating flow field turbulence

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102867094A (en) * 2012-09-19 2013-01-09 西安交通大学 Establishment method for free surface flow model in moving particle semi-implicit algorithm

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102867094A (en) * 2012-09-19 2013-01-09 西安交通大学 Establishment method for free surface flow model in moving particle semi-implicit algorithm

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JANG-OH MO ET AL.: ""Effects of wind speed changes on wake instability of a wind turbine in a virtual wind tunel using large eddy simulation"", 《JOURNAL OF WIND ENGINEERING AND INDUSTRIAL AERODYNAMICS》 *
刘克同 等: ""基于MRT-LBM大涡模拟的桥梁气动参数数值仿真"", 《四川大学学报(工程科学版)》 *
金斯科 等: ""单向连续车流下公路隧道紊流分区结构与流场特性"", 《南京工业大学学报(自然科学版)》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108133079A (en) * 2017-12-01 2018-06-08 上海理工大学 Field of turbulent flow coherent structure extracting method in IC engine cylinder
CN109858148A (en) * 2019-01-30 2019-06-07 南京航空航天大学 A kind of turbulent flow calculation method based on part filtering
CN109858148B (en) * 2019-01-30 2023-06-09 南京航空航天大学 Turbulence calculation method based on partial filtering
CN112131800A (en) * 2020-07-20 2020-12-25 中国科学院力学研究所 Novel large vortex simulation method and device based on energy flow similarity
CN112131800B (en) * 2020-07-20 2024-04-12 中国科学院力学研究所 Novel large vortex simulation method and device based on energy flow similarity
CN113111610A (en) * 2021-05-10 2021-07-13 中国空气动力研究与发展中心计算空气动力研究所 Novel sub-lattice scale model establishing method
CN113111610B (en) * 2021-05-10 2022-10-14 中国空气动力研究与发展中心计算空气动力研究所 Sub-lattice scale model establishing method
CN115859760A (en) * 2022-04-29 2023-03-28 中广核风电有限公司 Method and device for simulating flow field turbulence

Also Published As

Publication number Publication date
CN106372320B (en) 2020-03-20

Similar Documents

Publication Publication Date Title
CN106372320A (en) Method for performing large eddy simulation on highway tunnel turbulence by using sub-filtering scale model
CN103538068A (en) Fuzzy sliding mode trajectory tracking control method for SCARA robot
CN103473386B (en) A kind of method determining downburst wind profile of horizontal movement
CN109918821B (en) Windward conservation type river overflow outflow numerical simulation method
CN105047057B (en) A kind of high through street Macro-traffic Flow emulation mode of multi-user&#39;s multilane
CN106599456B (en) It is a kind of to distinguish the domatic and Geomorphologic Instantaneous Unit Hydrograph construction method of raceway groove conflux networks difference
CN103433924A (en) High-accuracy position control method for serial robot
CN107600176A (en) A kind of intelligent vehicle active steering control method theoretical based on extension control
CN105787959A (en) Method for multi-agent network object tracking based on improved adaptive particle filtering
CN106870620A (en) A kind of nonparametric dynamics calculation method of magnetic rheological liquid damper damping force
CN110110371A (en) Three-dimensional slope safety coefficient based on lower bound limit analysis theorem iteratively solves method
CN106325075A (en) H-infinity control method for time-delay linear parameter variation discrete system
CN104298797A (en) Method for determining remaining oil of fracture and hole type oil reservoir under encirclement of high flow guide channel
Braun et al. Finite element simulation of the wind action over bridge sectional models: Application to the Guamá River Bridge (Pará State, Brazil)
Tullis Size-scale effects of labyrinth weir hydraulics
Ho et al. Application of computational fluid dynamics to evaluate hydraulic performance of spillways in Australia
CN106295018A (en) A kind of transformator phantom construction method and system
CN103777521B (en) A kind of low speed control method of vehicle based on fuzzy control
Derickson et al. Development of a powerful hybrid tool for evaluating wind power in complex terrain: atmospheric numerical models and wind tunnels
Patel et al. Aerodynamic analysis of the utility truck with the morphing boom equipment
CN103150450B (en) Method for establishing dynamics real-time model of rail transit vehicle
CN116108965A (en) Method and device for calculating urban flood prediction model
KR101910274B1 (en) Method for simultaneous consideration of overland and underground space in flood modelling
CN111625987A (en) Turbulence model generation method based on deep learning
Wingstedt et al. Numerical simulations of particle dispersion in an urban area

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 210009 Zhongshan North Road, Jiangsu, No. 200,

Applicant after: Jin Sike

Applicant after: Taizhou Vocational & Technical College

Address before: 318001 Xueyuan Road, Jiaojiang, Zhejiang, No. 788,

Applicant before: Jin Sike

Applicant before: Taizhou Vocational & Technical College

COR Change of bibliographic data
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200320