CN107133418A - Fluids material stratospheric transport analogy method based on alternately TVD difference algorithms - Google Patents
Fluids material stratospheric transport analogy method based on alternately TVD difference algorithms Download PDFInfo
- Publication number
- CN107133418A CN107133418A CN201710383303.XA CN201710383303A CN107133418A CN 107133418 A CN107133418 A CN 107133418A CN 201710383303 A CN201710383303 A CN 201710383303A CN 107133418 A CN107133418 A CN 107133418A
- Authority
- CN
- China
- Prior art keywords
- transport
- fluids
- fluids material
- tvd
- stratospheric
- 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
- 239000000463 material Substances 0.000 title claims abstract description 85
- 239000012530 fluid Substances 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 36
- 230000008569 process Effects 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 8
- 238000009792 diffusion process Methods 0.000 claims description 18
- 210000003205 muscle Anatomy 0.000 claims description 4
- 230000032258 transport Effects 0.000 description 21
- 230000000694 effects Effects 0.000 description 8
- 238000007711 solidification Methods 0.000 description 7
- 230000008023 solidification Effects 0.000 description 7
- 238000004088 simulation Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 239000006185 dispersion Substances 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000003467 diminishing effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005194 fractionation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses the advection transport equation that a kind of Fluids material stratospheric transport analogy method based on alternately TVD difference algorithms includes building Fluids material, and then numerical discretization is carried out to the advection transport equation, obtain the discrete equation of Fluids material conveying;The discrete equation that Fluids material is conveyed disassembles into the step number equal with its dimension;Obtain the flow field data of Fluids material and input and iterative numerical calculating is carried out in the discrete equation conveyed to Fluids material, predict the Fluids material concentration in each default iteration time section;During iterative numerical is calculated, different types of total variance abatement restricted function is used alternatingly according to the step number equal with its dimension in each dimension respectively in different default iteration time sections to calculate Fluids material concentration in next period, the change of Fluids material stratospheric transport process medium fluid material concentration is finally simulated and predict.
Description
Technical field
Calculated the invention belongs to fluid dynamic numerical simulation and calculating field, more particularly to a kind of alternately TVD difference that is based on
The Fluids material stratospheric transport analogy method of method.
Background technology
Numerical model is research and the important work for forecasting the motion such as river, ocean and its interior transport of substances in nature
Tool, with the development of computer technology and social economy, numerical model is in scientific research, operational forecast and Hu Ku, ocean etc.
The effect become more and more important is played in the industry such as environmental management and planning.
Advection is that the main process that material carries out space conveying is carried in Fluids, and advection process is represented by equation
(1).It is generally using the method for difference that the progress of (1) formula is discrete in numerical model, obtain such as the advection difference side of formula (2) form
Journey, then carries out iterative numerical to formula (2), obtains material concentration c spatial distribution and change over time.
Wherein, c transports material concentration for certain, and u, v, w are respectively the flow velocity in x, y, z direction.
Wherein,For the material concentration at half storey grid, Δ t is discrete
Time step, Δ x, Δ y, Δ z are respectively the grid spacing in x, y, z direction, and label subscript n represents time step number, subscript i,
J, k represent the grid label in x, y, z direction respectively.
In difference equation (2), the material concentration at half storey grid how is calculatedIt is the main difference of different difference schemes
Where not.For example, upstreame scheme uses upstream grid concentration, middle poor form uses being averaged for upstream and downstream grid
Concentration.Although these algorithms calculate easy, the problem of there is low precision or strong numerical solidification, cause analog result distortion compared with
It is strong and the stability of model can be influenceed.
Total variance abatement (Total Variation Diminishing, abbreviation TVD) difference scheme is the high-precision nothing of a class
The advection difference algorithm of frequency dispersion, it is widely used in ocean and Atmospheric Numerical Model and correlative study.TVD algorithms are mainly
The high-precision item of back-diffusion is added on the basis of original low precision upstreame scheme, the precision of advection form is improved, draws simultaneously
Enter restricted function to control the intensity of high-precision item added, with this, keep in whole integral process forward, numeric format it is total
It is deteriorated and keeps not increasing, so as to effectively suppresses the numerical solidification that high accurate scheme generally has.
The algorithm of TVD advection forms is following (by taking x directions as an example):
As shown in formula (3), calculatingWhen, plus high-precision correction term on the basis of upstreame scheme, and before correction term
Add restricted function ψ (ri+1/2), thus obtain:
Wherein,
Add after restricted function and ensure the total variance of result of calculationHolding does not increase, i.e.,:
TV(cn+1)≤TV(cn) (4)
Restricted function need to be caused to meet following conditions:
Wherein, r represents concentration gradient ratio, i.e.,Or
Although TVD restricted functions have a variety of, a kind of restricted function of single use can not be controlled in TVD forms well
The intensity of the high-precision item of addition, thereby resulted in above TVD numeric formats still have certain numerical d ispersion (e.g., Minmod,
Van Leer, MUSCL, HSIMT) and the pseudo- effect of back-diffusion (e.g., Superbee).In addition, TVD forms can ensure one-dimensional before
In the case of total variance do not increase, it is and then barely satisfactory in the case of multidimensional, however it remains numerical solidification.Therefore, existing TVD is utilized
When advection difference scheme to the transport of substances in Fluids simulate for a long time, unsatisfactory simulation knot is still suffered from
Really.
The content of the invention
In order to solve the deficiencies in the prior art, the invention provides a kind of Fluids based on alternately TVD difference algorithms
Material stratospheric transport analogy method, it provides a kind of effectively decrease numerical d ispersion and the pseudo- effect of back-diffusion and solves multidimensional frequency dispersion and puts down
Difference method-the TVDal (TVD with Alternating Limiters) of stream format, it is right in Fluids to be applied to
The simulation and forecast of transport of substances.
A kind of Fluids material stratospheric transport analogy method based on alternately TVD difference algorithms of the present invention, including:
The advection transport equation of Fluids material is built, and then numerical discretization is carried out to the advection transport equation, is obtained
The discrete equation conveyed to Fluids material;The dimension of the discrete equation of the Fluids material conveying is at least two dimension;
The discrete equation that Fluids material is conveyed disassembles into the step number equal with its dimension;
Obtain the flow field data of Fluids material and input in the discrete equation conveyed to Fluids material and enter line number
Value iterative calculation, predicts the Fluids material concentration in each default iteration time section;
During iterative numerical is calculated, in different default iteration time sections respectively in each dimension according to
The equal step number of its dimension is used alternatingly different types of total variance abatement restricted function to calculate the earth in next period
Flowing material concentration, finally simulates and predicts the change of Fluids material stratospheric transport process medium fluid material concentration.
Further, the building method based on the abatement form that is deteriorated, the advection transport equation to earth flowing material is carried out
Numerical discretization.
Further, total variance abatement restricted function includes the restricted function and numerical value back-diffusion property of numerical d ispersion property
Restricted function.
Further, Fluids material concentration is the function on time and space.
Further, the process for asking for the advection transport equation of Fluids material is:
Fluids material concentration is respectively to the derived function for the directioin parameter for characterizing space, then the side with corresponding sign space
After being multiplied to the flow velocity of parameter, the Fluids material concentration that finally adds up obtains Fluids material to the derived function of time
Advection transport equation.
Wherein, the restricted function of the numerical d ispersion property is Minmod functions.
The restricted function of the numerical d ispersion property is Van Leer functions.
The restricted function of the numerical d ispersion property is MUSCL functions.
The restricted function of the numerical d ispersion property is HSIMT functions.
The restricted function of the numerical value back-diffusion property includes Superbee functions.
Wherein, actually calculate in, frequently with several classical TVD restricted functions mainly have:
Superbee:
ψ(r)Superbee=max [0, min (2r, 1), min (r, 2)] (6)
Van Leer(or Harmonic):
ψ(r)VanLeer=(r+ | r |)/(r+1) (7)
MUSCL:
Minmod:
ψ(r)Minmod=max [0, min (r, 1)] (9)
HSIMT:
ψ(r)HSIMT=max [0, min (2r, 2, β)] (10)
Wherein,
Compared with prior art, the beneficial effects of the invention are as follows:
(1) this method does not change original formula increase extra computation amount, but is used alternatingly not in different time step
Connatural algorithm, largely reduces the numerical d ispersion and the pseudo- effect of back-diffusion of traditional TVD advections difference algorithm so that
The error of analog result is smaller.
(2) alternately restricted function is split with discrete equation in this method and be combined, it is flat without frequency dispersion based on obtained multidimensional
Flow algorithm, effectively eliminates the TVD advection forms numerical solidification problem that transport of substances is simulated in the case of Multi-dimension calculation.
Brief description of the drawings
The Figure of description for constituting the part of the application is used for providing further understanding of the present application, and the application's shows
Meaning property embodiment and its illustrate be used for explain the application, do not constitute the improper restriction to the application.
Fig. 1 is the schematic flow sheet of the present invention;
Fig. 2 (a) is initial time material concentration distribution map;
Fig. 2 (b) is the Superbee analog results based on traditional TVD algorithms;
Fig. 2 (c) is the HSIMT analog results based on traditional TVD algorithms;
Fig. 2 (d) is the Superbee analog results that traditional TVD algorithms combination discrete equation splits algorithm;
Fig. 2 (e) is the HSIMT analog results that traditional TVD algorithms combination discrete equation splits algorithm;
Fig. 2 (f) is the analog result of the inventive method.
Embodiment
It is noted that described further below is all exemplary, it is intended to provide further instruction to the application.Unless another
Indicate, all technologies used herein and scientific terminology are with usual with the application person of an ordinary skill in the technical field
The identical meanings of understanding.
It should be noted that term used herein above is merely to describe embodiment, and be not intended to restricted root
According to the illustrative embodiments of the application.As used herein, unless the context clearly indicates otherwise, otherwise singulative
It is also intended to include plural form, additionally, it should be understood that, when in this manual using term "comprising" and/or " bag
Include " when, it indicates existing characteristics, step, operation, device, component and/or combinations thereof.
The dimension of discrete equation of Fluids material conveying is at least two dimension, below with Fluids material convey from
Exemplified by dissipating the dimension of equation and being two dimension:
As shown in figure 1, the Fluids material stratospheric transport analogy method based on alternately TVD difference algorithms of the present invention,
Including:
Step 1:Build Fluids material advection transport equation, and then to the advection transport equation carry out numerical value from
Dissipate, obtain the discrete equation of Fluids material conveying.
The advection transport equation for building Fluids material is formula (1):
Numerical discretization is carried out to the advection transport equation, the discrete equation for obtaining the conveying of Fluids material is formula
(2):
Wherein,For the material concentration at half storey grid, Δ t is discrete time step,
Δ x, Δ y are respectively the grid spacing in x, y direction, and label subscript n represents time step number, and subscript i, j represents x, y direction respectively
Grid label.
Step 2:The discrete equation that Fluids material is conveyed disassembles into the step number equal with its dimension.
Step 3:Obtain the flow field data of Fluids material and input in the discrete equation conveyed to Fluids material
Iterative numerical calculating is carried out, the Fluids material concentration in each default iteration time section is predicted;
During iterative numerical is calculated, in different default iteration time sections respectively in each dimension according to
The equal step number of its dimension is used alternatingly different types of total variance abatement restricted function to calculate the earth in next period
Flowing material concentration, finally simulates and predicts the change of Fluids material stratospheric transport process medium fluid material concentration.
Specifically, flow velocity u, v information is inputted, the parity of iterative steps is iterated and judged to discrete equation, if
Odd number is calculated then to be calculated using below step a), if even number is then calculated using step b).
Step a):
Discrete equation formula (2) is disassembled as formula (3-1) and the calculating of the step of formula (3-2) two:
Wherein,Subscript com and dif represent respectively using numerical value back-diffusion property (for example, Superbee) and numerical value expansion
The TVD restricted functions of property (for example, Minmod) are dissipated to calculate material concentration at half storey grid
Discrete equation is disassembled into (3-1) and (3-2) two step, advection process in x and y both directions is calculated successively:First
When step calculates advection process on x directions, first using a kind of TVD restricted functions of numerical value back-diffusion property, calculating obtains middle change
Measure c*;Intermediate variable c* is substituted into y directions difference equation by second step, is calculated advection process on y directions, is used a kind of numerical d ispersion
The TVD restricted functions of property calculate the concentration c for obtaining new time stepn+1。
Step b):
Equally it is divided to (4-1) and (4-2) two step to calculate, but the first step first calculates y directions advection process, it is anti-using a kind of numerical value
The TVD restricted functions of scattering nature, calculating obtains intermediate variable c*;Intermediate variable c* is substituted into x directions difference equation by second step,
X equation advection processes are calculated, using a kind of TVD restricted functions of numerical d ispersion property, the material that calculating obtains new time step is dense
Degree.
Wherein,Subscript com and dif represent respectively using numerical value back-diffusion property (for example, Superbee) and numerical value expansion
The TVD restricted functions of property (for example, Minmod) are dissipated to calculate material concentration at half storey grid
Finally judge to integrate the setting value whether duration reaches pattern, no then return to step 2 continues, terminated if reaching
Calculate.
Present disclosure can be further illustrated by following analog result.
1st, simulated experiment:Using the inventive method and other method, the material advection process to giving a certain concentration distribution
Simulated.
2nd, analog result
Exist in a horizontal two-dimension marine site, in it with 12 hours for reciprocal flowing of the cycle along 45 ° of directions, initial time
Material concentration distribution is given to be distributed as shown in Fig. 2 (a), and knot is simulated after 100 cycles obtained using different advection difference algorithms
Fruit is respectively as shown in Fig. 2 (b)-Fig. 2 (f).Wherein, Fig. 2 (b) and 2 (c) are based on traditional TVD algorithms Superbee and HSIMT
Analog result, Fig. 2 (d) and Superbee and the HSIMT simulation that Fig. 2 (e) is that tradition TVD algorithm combinations discrete equation splits algorithm
As a result;Fig. 2 (f) is the analog result of the inventive method.
From Fig. 2 (a)-Fig. 2 (f) as can be seen that the analog result based on traditional TVD algorithms has obvious numerical value frequency
Dissipate, splitting numerical solidification after algorithm using discrete equation is improved, but still suffers from the pseudo- effect of numerical value back-diffusion and diffusion, is made
Then maintain the distribution of true solution substantially with the inventive method, reduce the pseudo- effect of numerical value back-diffusion and diffusion of traditional TVD algorithms,
Result of calculation achieves gratifying analog result without numerical solidification simultaneously.
To sum up, this method proposed by the present invention effectively weakens numerical d ispersion and the pseudo- effect of back-diffusion and solves multidimensional frequency dispersion
Difference algorithm-the TVDal of advection form material advection analogy method.The method of the present invention is conversely pseudo- by the way that property is used alternatingly
Traditional TVD restricted functions of effect, effectively reduce numerical d ispersion and back-diffusion, reduce the error of advection analog result,
Meanwhile, Multidimensional Discrete equation is subjected to fractionation calculating, it is ensured that multidimensional advective simulation is that no numerical solidification is produced, and the above is improved logical
The analog result for crossing the experiment of material advection can be seen that.
Although above-mentioned the embodiment of the present invention is described with reference to accompanying drawing, not to present invention protection model
The limitation enclosed, one of ordinary skill in the art should be understood that on the basis of technical scheme those skilled in the art are not
Need to pay various modifications or deform still within protection scope of the present invention that creative work can make.
Claims (10)
1. a kind of Fluids material stratospheric transport analogy method based on alternately TVD difference algorithms, it is characterised in that including:
The advection transport equation of Fluids material is built, and then numerical discretization is carried out to the advection transport equation, obtains ground
The discrete equation of ball flowing material conveying;The dimension of the discrete equation of the Fluids material conveying is at least two dimension;
The discrete equation that Fluids material is conveyed disassembles into the step number equal with its dimension;
Obtain the flow field data of Fluids material and input progress numerical value in the discrete equation conveyed to Fluids material and change
In generation, calculates, and predicts the Fluids material concentration in each default iteration time section;
During iterative numerical is calculated, tieed up respectively in each dimension according to it in different default iteration time sections
The equal step number of number is used alternatingly different types of total variance abatement restricted function to calculate Fluids in next period
Material concentration, finally simulates and predicts the change of Fluids material stratospheric transport process medium fluid material concentration.
2. the Fluids material stratospheric transport analogy method as claimed in claim 1 based on alternately TVD difference algorithms, it is special
Levy and be, based on the building method for the abatement form that is deteriorated, the advection transport equation to earth flowing material carries out numerical discretization.
3. the Fluids material stratospheric transport analogy method as claimed in claim 1 based on alternately TVD difference algorithms, it is special
Levy and be, total variance abatement restricted function includes the restricted function of numerical d ispersion property and the limitation letter of numerical value back-diffusion property
Number.
4. the Fluids material stratospheric transport analogy method as claimed in claim 1 based on alternately TVD difference algorithms, it is special
Levy and be, Fluids material concentration is the function on time and space.
5. the Fluids material stratospheric transport analogy method as claimed in claim 4 based on alternately TVD difference algorithms, it is special
Levy and be, the process for asking for the advection transport equation of Fluids material is:
Fluids material concentration is joined respectively to the derived function for the directioin parameter for characterizing space, then with the corresponding direction for characterizing space
After several flow velocitys is multiplied, the Fluids material concentration that finally adds up obtains the advection of Fluids material to the derived function of time
Transport equation.
6. the Fluids material stratospheric transport analogy method as claimed in claim 3 based on alternately TVD difference algorithms, it is special
Levy and be, the restricted function of the numerical d ispersion property is for Minmod functions or from his scattering nature restricted function.
7. the Fluids material stratospheric transport analogy method as claimed in claim 3 based on alternately TVD difference algorithms, it is special
Levy and be, the restricted function of the numerical d ispersion property is for Van Leer functions or from his scattering nature restricted function.
8. the Fluids material stratospheric transport analogy method as claimed in claim 3 based on alternately TVD difference algorithms, it is special
Levy and be, the restricted function of the numerical d ispersion property is for MUSCL functions or from his scattering nature restricted function.
9. the Fluids material stratospheric transport analogy method as claimed in claim 3 based on alternately TVD difference algorithms, it is special
Levy and be, the restricted function of the numerical d ispersion property is for HSIMT functions or from his scattering nature restricted function.
10. the Fluids material stratospheric transport analogy method as claimed in claim 3 based on alternately TVD difference algorithms, its
It is characterised by, the restricted function of the numerical value back-diffusion property includes Superbee functions or limits letter from his back-diffusion property
Number.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710383303.XA CN107133418B (en) | 2017-05-26 | 2017-05-26 | Earth fluid material advection transportation simulation method based on alternative TVD differential algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710383303.XA CN107133418B (en) | 2017-05-26 | 2017-05-26 | Earth fluid material advection transportation simulation method based on alternative TVD differential algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107133418A true CN107133418A (en) | 2017-09-05 |
CN107133418B CN107133418B (en) | 2020-04-21 |
Family
ID=59733351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710383303.XA Expired - Fee Related CN107133418B (en) | 2017-05-26 | 2017-05-26 | Earth fluid material advection transportation simulation method based on alternative TVD differential algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107133418B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114840974A (en) * | 2022-03-25 | 2022-08-02 | 中国气象局地球***数值预报中心 | Advection mode system suitable for complex terrain and operation method thereof |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101484922A (en) * | 2006-04-05 | 2009-07-15 | 财团法人Seoul大学校产学协力财团 | Method of simulating detailed movements of fluids using derivative particles |
CN101615219A (en) * | 2009-08-04 | 2009-12-30 | 北京师范大学 | A kind of high-precision difference method of simulating transport diffusion problem |
CN102830121A (en) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | Real-time detection method of soft abrasive flow abrasive group |
CN104239625A (en) * | 2014-09-05 | 2014-12-24 | 中国矿业大学 | Corrective fluid motion equation linear iteration-based steady state solution method |
CN104636566A (en) * | 2015-03-13 | 2015-05-20 | 中国矿业大学 | Dynamic mesh numerical solution method based on modified fluid motion equation |
CN104732093A (en) * | 2015-03-30 | 2015-06-24 | 西安交通大学 | FCT-FDM forward simulation method based on dispersion viscosity wave equation |
US9430870B2 (en) * | 2003-02-19 | 2016-08-30 | California Institute Of Technology | Level set surface editing operators |
CN106294282A (en) * | 2016-08-01 | 2017-01-04 | 中国石油天然气股份有限公司 | Black oil reservoir simulation method and device |
-
2017
- 2017-05-26 CN CN201710383303.XA patent/CN107133418B/en not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9430870B2 (en) * | 2003-02-19 | 2016-08-30 | California Institute Of Technology | Level set surface editing operators |
CN101484922A (en) * | 2006-04-05 | 2009-07-15 | 财团法人Seoul大学校产学协力财团 | Method of simulating detailed movements of fluids using derivative particles |
CN101615219A (en) * | 2009-08-04 | 2009-12-30 | 北京师范大学 | A kind of high-precision difference method of simulating transport diffusion problem |
CN102830121A (en) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | Real-time detection method of soft abrasive flow abrasive group |
CN104239625A (en) * | 2014-09-05 | 2014-12-24 | 中国矿业大学 | Corrective fluid motion equation linear iteration-based steady state solution method |
CN104636566A (en) * | 2015-03-13 | 2015-05-20 | 中国矿业大学 | Dynamic mesh numerical solution method based on modified fluid motion equation |
CN104732093A (en) * | 2015-03-30 | 2015-06-24 | 西安交通大学 | FCT-FDM forward simulation method based on dispersion viscosity wave equation |
CN106294282A (en) * | 2016-08-01 | 2017-01-04 | 中国石油天然气股份有限公司 | Black oil reservoir simulation method and device |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114840974A (en) * | 2022-03-25 | 2022-08-02 | 中国气象局地球***数值预报中心 | Advection mode system suitable for complex terrain and operation method thereof |
CN114840974B (en) * | 2022-03-25 | 2022-11-25 | 中国气象局地球***数值预报中心 | Advection mode system suitable for complex terrain and operation method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN107133418B (en) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104268943B (en) | Fluid simulation method based on Eulerian-Lagrangian coupling method | |
Wang | A Quadtree-based adaptive Cartesian/Quad grid flow solver for Navier-Stokes equations | |
Morgenthal et al. | An immersed interface method for the Vortex-In-Cell algorithm | |
Hartmann et al. | Discontinuous Galerkin methods for computational aerodynamics—3D adaptive flow simulation with the DLR PADGE code | |
CN113850008B (en) | Self-adaptive grid disturbance domain updating acceleration method for aircraft aerodynamic characteristic prediction | |
CN111859529B (en) | Multi-grid disturbance domain updating acceleration method for aircraft streaming numerical simulation | |
CN113609598B (en) | RANS/LES disturbance domain updating method for aircraft aerodynamic characteristic simulation | |
CN107657075B (en) | The Region Decomposition FInite Element of Darcy velocity at simulated groundwater dielectric interface | |
Bermejo et al. | A conservative quasi-monotone semi-Lagrangian scheme | |
Giles et al. | Progress in adjoint error correction for integral functionals | |
Suchde et al. | Point cloud movement for fully Lagrangian meshfree methods | |
Thömmes et al. | A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method | |
Miller et al. | Surface heterogeneity effects on regional-scale fluxes in the stable boundary layer: Aerodynamic roughness length transitions | |
Moukalled et al. | Performance comparison of the NWF and DC methods for implementing high-resolution schemes in a fully coupled incompressible flow solver | |
CN107133418A (en) | Fluids material stratospheric transport analogy method based on alternately TVD difference algorithms | |
Okusanya et al. | Algebraic multigrid for stabilized finite element discretizations of the Navier–Stokes equations | |
Liu et al. | An unconditionally stable spline difference scheme of O (k2+ h4) for solving the second-order 1D linear hyperbolic equation | |
Parent | Positivity-preserving flux difference splitting schemes | |
Sonar | Classical finite volume methods | |
Jiang et al. | A new turbulence model for separated flows | |
Shanbi et al. | The simulation of natural gas gathering pipeline network | |
CN108072899B (en) | Self-adaptive implementation method of discontinuous Galerkin finite element seismic numerical simulation algorithm | |
Hong et al. | A sharp interface immersed edge-based smoothed finite element method with extended fictitious domain scheme | |
Shi et al. | An adaptive semi-lagrangian level-set method for convection-diffusion equations on evolving interfaces | |
Gillow et al. | Computation of currents at microelectrodes using hp-DGFEM |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200421 |