CN116542178A - Constant flow numerical simulation method, device, equipment and medium - Google Patents

Constant flow numerical simulation method, device, equipment and medium Download PDF

Info

Publication number
CN116542178A
CN116542178A CN202310591185.7A CN202310591185A CN116542178A CN 116542178 A CN116542178 A CN 116542178A CN 202310591185 A CN202310591185 A CN 202310591185A CN 116542178 A CN116542178 A CN 116542178A
Authority
CN
China
Prior art keywords
flow field
updated
iteration
iteration step
round
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202310591185.7A
Other languages
Chinese (zh)
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.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
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 Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202310591185.7A priority Critical patent/CN116542178A/en
Publication of CN116542178A publication Critical patent/CN116542178A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • 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
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application discloses a steady flow numerical simulation method, a device, equipment and a medium, and relates to the field of computational fluid mechanics research, wherein the method comprises the following steps: sequentially updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain updated computing units and updated current flow fields in the iteration step of the round; judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result. Through the scheme, the flow field meeting the convergence condition can be rapidly obtained in the steady flow numerical simulation process.

Description

Constant flow numerical simulation method, device, equipment and medium
Technical Field
The invention relates to the field of computational fluid dynamics research, in particular to a steady flow numerical simulation method, a steady flow numerical simulation device, steady flow numerical simulation equipment and a steady flow numerical simulation medium.
Background
Explicit format is a commonly used time format in computational fluid dynamics iterative type methods. When the explicit format advances the time, the used variables are all known values of the current moment, so that the variable value of the next time step can be obtained by directly advancing. The explicit format has the advantages of simple principle, easy implementation and small memory. However, the advance time step of the explicit format is usually small, limited by the numerical stability, and for steady flow problems, a large number of iterations are required to reach steady state, and thus the computational efficiency is low.
The iterative convergence process of the steady flow field is the development process of the initial field to the final solution under the constraint of the boundary condition. The standard explicit format uses a synchronous pushing method when time pushing, and all flow field units directly use the flow field value before the current iteration update to perform unified calculation, and the flow field value of the updated unit in the iteration process is not adopted. In general, the flow field values of the updated elements contain more fully developed information, whereas the standard explicit format does not take this important information into account during the iteration process, thus not favoring the rapid convergence of the steady flow.
In summary, how to quickly obtain a flow field meeting the convergence condition in the steady flow numerical simulation process is a problem to be solved in the field.
Disclosure of Invention
Accordingly, the present invention is directed to a method, apparatus, device and medium for simulating steady flow values, which can rapidly obtain a flow field satisfying a convergence condition during steady flow value simulation. The specific scheme is as follows:
in a first aspect, the present application discloses a steady flow numerical simulation method, comprising:
sequentially updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain updated computing units and updated current flow fields in the iteration step of the round;
judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result.
Optionally, before updating each un-updated calculation unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, the method further includes:
setting a surface boundary condition of an aircraft to be simulated and a boundary condition of the periphery of a calculation domain to obtain a target boundary condition; the boundary conditions of the periphery of the calculation domain comprise far field conditions, inflow conditions, outflow conditions and symmetrical plane conditions.
Optionally, the determining whether the preset iteration stop condition is met currently or not, if not, starting a next iteration step, including:
judging whether the current total iteration times are larger than a preset iteration threshold, and if not, starting the next iteration step.
Optionally, the determining whether the preset iteration stop condition is met currently or not, if not, starting a next iteration step, including:
judging whether the updated current flow field in the iteration step of the round meets a preset convergence condition, and if not, starting the next iteration step.
Optionally, the step of determining whether the updated current flow field in the current iteration step meets a preset convergence condition, and if not, starting a next iteration step includes:
acquiring the residual error of the current flow field updated in the iteration step of the round, and judging whether the residual error meets a preset convergence condition;
if not, starting the next iteration step.
Optionally, the updating each un-updated calculation unit in the iteration step of the present round sequentially based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain an updated calculation unit in the iteration step of the present round and the updated current flow field, including:
judging whether updated computing units and non-updated computing units exist in the iteration steps of the round, if so, determining the current computing unit to be updated in the iteration steps of the round from the non-updated computing units of the iteration steps of the round;
updating the current computing unit to be updated in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field values of all computing units in the peripheral computing template range of the current computing unit to be updated in the iteration step of the round, and then re-jumping to the step of judging whether updated computing units and non-updated computing units exist in the iteration step of the round or not until all computing units in the iteration step of the round are updated.
Optionally, before updating the current calculation unit to be updated in the iteration step of this round based on the target boundary condition, the target flow field unit updating function and the flow field value of the updated calculation unit, the method further includes:
and determining an initial flow field value corresponding to each calculation unit in the initial flow field, and determining a current calculation unit to be updated of an initial iteration step from each calculation unit.
In a second aspect, the present application discloses a steady flow numerical simulation device comprising:
the flow field iteration module is used for sequentially updating each non-updated calculation unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field so as to obtain an updated calculation unit in the iteration step of the round and the updated current flow field;
the simulation result obtaining module is used for judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step in turn based on the target boundary condition, the target flow field unit updating function and the current flow field until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result.
In a third aspect, the present application discloses an electronic device comprising:
a memory for storing a computer program;
and a processor for executing the computer program to implement the steps of the steady flow numerical simulation method disclosed above.
In a fourth aspect, the present application discloses a computer-readable storage medium for storing a computer program; wherein the computer program when executed by a processor implements the steps of the steady flow numerical simulation method disclosed previously.
The beneficial effects of the application are that: sequentially updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain updated computing units and updated current flow fields in the iteration step of the round; judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result. Therefore, when the flow field is updated, each non-updated computing unit in the iteration step of the round needs to be updated based on the current flow field, specifically, when the non-updated computing units are sequentially updated in the iteration of the round, if other computing units located in the computing template range of each unit to be updated have completed the iteration updating of the round, the latest flow field value of the corresponding unit in the computing template range needs to be adopted for computing, so that updated flow field information can be fully utilized, the target flow field can be acquired more quickly, and the effect of accelerating convergence is achieved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only embodiments of the present invention, and that other drawings may be obtained according to the provided drawings without inventive effort to a person skilled in the art.
FIG. 1 is a schematic flow diagram of a steady flow numerical simulation method disclosed in the present application;
FIG. 2 is a schematic flow diagram of a specific steady flow numerical simulation method disclosed in the present application;
FIG. 3 is a flow chart of another exemplary steady flow numerical simulation method disclosed in the present application;
FIG. 4 is a schematic diagram of a particular steady flow numerical simulation disclosed herein;
FIG. 5 is a schematic diagram of a specific residual convergence curve disclosed herein;
FIG. 6 is a schematic diagram of a steady flow numerical simulation device disclosed in the present application;
fig. 7 is a block diagram of an electronic device disclosed in the present application.
Detailed Description
The following description of the technical solutions in the embodiments of the present application will be made clearly and completely with reference to the drawings in the embodiments of the present application, and it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Explicit format is a commonly used time format in computational fluid dynamics iterative type methods. When the explicit format advances the time, the used variables are all known values of the current moment, so that the variable value of the next time step can be obtained by directly advancing. The explicit format has the advantages of simple principle, easy implementation and small memory. However, the advance time step of the explicit format is usually small, limited by the numerical stability, and for steady flow problems, a large number of iterations are required to reach steady state, and thus the computational efficiency is low.
The iterative convergence process of the steady flow field is the development process of the initial field to the final solution under the constraint of the boundary condition. The standard explicit format uses a synchronous pushing method when time pushing, and all flow field units directly use the flow field value before the current iteration update to perform unified calculation, and the flow field value of the updated unit in the iteration process is not adopted. In general, the flow field values of the updated elements contain more fully developed information, whereas the standard explicit format does not take this important information into account during the iteration process, thus not favoring the rapid convergence of the steady flow.
Therefore, the application correspondingly provides a steady flow numerical simulation scheme which can rapidly acquire the flow field meeting the convergence condition in the steady flow numerical simulation process.
Referring to fig. 1, an embodiment of the present application discloses a steady flow numerical simulation method, including:
step S11: and updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence to obtain an updated computing unit and an updated current flow field in the iteration step of the round.
In this embodiment, before the updating of each non-updated computing unit in the iteration step of this round based on the target boundary condition, the target flow field unit updating function and the current flow field sequentially, the method further includes: setting a surface boundary condition of an aircraft to be simulated and a boundary condition of the periphery of a calculation domain to obtain a target boundary condition; the boundary conditions of the periphery of the calculation domain comprise far field conditions, inflow conditions, outflow conditions and symmetrical plane conditions. It should be noted that, before updating the computing unit, firstly, for the steady flow problem to be solved, the physical space in the computing domain is discretized by means of commercial software and the like to generate a computing grid file containing space discrete units, and grid information is obtained by means of file reading; secondly, boundary conditions required for flow simulation are set according to the characteristics of the solved problem, and may include an object plane boundary condition, a far field boundary condition, an inflow boundary condition, an outflow boundary condition, a symmetry plane boundary condition and the like.
The current flow field is the latest flow field, that is, if the current iteration is the first iteration, the latest flow field is the initialized flow field, for example, the 1 st calculation unit of the first iteration is updated based on the target boundary condition and the target flow field unit updating function on the assumption that the total number of calculation units is 5; when updating the 1 st computing unit, when updating the 2 nd computing unit, the first iteration is ended based on the updated 1 st computing unit until all 5 computing units are updated.
Step S12: judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result.
It can be appreciated that in the process of obtaining the target flow field, a plurality of iterations are required, and an iteration stop condition can be preset to stop starting the next iteration step. If the preset iteration stop condition is not met, a new round of updating is needed to be carried out on each computing unit again. Because the latest flow field is referred to each time the computing unit is updated, that is, the updated flow field value generally contains more fully developed information and is closer to the final solution, the target flow field, that is, the corresponding flow field simulation result, can be obtained more quickly.
The beneficial effects of the application are that: sequentially updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain updated computing units and updated current flow fields in the iteration step of the round; judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result. Therefore, when the flow field is updated, each non-updated computing unit in the iteration step of the round needs to be updated based on the current flow field, specifically, when the non-updated computing units are sequentially updated in the iteration of the round, if other computing units located in the computing template range of each unit to be updated have completed the iteration updating of the round, the latest flow field value of the corresponding unit in the computing template range needs to be adopted for computing, so that updated flow field information can be fully utilized, the target flow field can be acquired more quickly, and the effect of accelerating convergence is achieved.
Referring to fig. 2, an embodiment of the present application discloses a specific steady flow numerical simulation method, which includes:
step S21: and updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence to obtain an updated computing unit and an updated current flow field in the iteration step of the round.
In this embodiment, the updating each non-updated calculating unit in the iteration step of the present round based on the target boundary condition, the target flow field unit updating function and the current flow field sequentially to obtain an updated calculating unit and an updated current flow field in the iteration step of the present round includes: judging whether updated computing units and non-updated computing units exist in the iteration steps of the round, if so, determining the current computing unit to be updated in the iteration steps of the round from the non-updated computing units of the iteration steps of the round; updating the current computing unit to be updated in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field values of all computing units in the peripheral computing template range of the current computing unit to be updated in the iteration step of the round, and then re-jumping to the step of judging whether updated computing units and non-updated computing units exist in the iteration step of the round or not until all computing units in the iteration step of the round are updated.
It may be understood that, in a preset sequence, each un-updated computing unit in the iteration step of the present round is sequentially taken as a current computing unit to be updated in the iteration step of the present round, the current computing unit to be updated and surrounding units thereof are selected as computing templates, for all computing units within the range of the computing templates, if the updating of the iteration step of the present round has been completed, the current computing unit to be updated in the iteration step of the present round is updated by using the flow field value updated in the iteration step of the present round, and if the updating of the iteration step of the present round has not been completed, the current computing unit to be updated in the iteration step of the present round is updated by using the flow field value updated in the previous iteration step of the present round, where the target flow field unit updating function may be expressed as:
in the method, in the process of the invention,representing the flow field value of the I-th calculation unit at the nth iteration step, namely the flow field value of the calculation unit currently to be solved,/for>Representing flow field values of computation units for which the I-th computation unit periphery has completed the n-th iteration step update,flow field value of calculation unit indicating that the I-th calculation unit has not completed the n-th iteration step update,/->Flow field value f representing calculation unit whose periphery is not completed with nth iteration step update Explicit (. Cndot.) is an iterative formula for flow field elements in a standard explicit format.
For example, there are 5 computing units in total, the iteration step of this round is the 3 rd iteration, and in the 3 rd iteration, the 1 st computing unit and the 2 nd computing unit have been updated, i.e. n is 3, i is 3, and the iteratively updated template includes 2 units on the left and right sides in addition to the 3 rd unit, so that the updating function of the target flow field unit is as follows:
it can be understood that the argument of the update function of the target flow field unit is the flow field value of the unit to be updated and its peripheral units at the current iteration step, and the argument is the flow field value of the updated unit, that is,the update has been completed in this round of iteration, thus +.>Is a known quantity, +.>And->The update has been completed in the previous iteration, so +.>Andthe flow field value is also known, so that independent variables in the updating function of the target flow field unit are known, and the current flow field value to be updated in the iteration step of the round can be directly solved.
In this embodiment, before updating the current calculation unit to be updated in the present iteration step based on the target boundary condition, the target flow field unit updating function, and the flow field value of the updated calculation unit, the method further includes: and determining an initial flow field value corresponding to each calculation unit in the initial flow field, and determining a current calculation unit to be updated of an initial iteration step from each calculation unit. For example, 5 calculation units are numbered in advance to obtain a unit update sequence, for example, the unit update sequence can be defined according to the sequence from small unit number to large unit number, the sequence from large unit number to small unit number, and the unit update sequence can be defined.
Step S22: judging whether the current total iteration times are larger than a preset iteration threshold.
Step S23: if the current flow field is not greater than the preset iteration stop condition, starting a next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result.
For example, the preset iteration threshold is set to 100 times, which means that the maximum allowed iteration number is 100, and if 99 iterations are completed currently, 1 iteration can be performed; if 100 iterations have been completed, the iteration is stopped and the current flow field is output as the target flow field.
Therefore, when the current computing unit of the iteration step of the current round is updated, the latest flow field value of the computing unit in the template range of the computing unit to be updated currently is needed to be utilized instead of the flow field value of the previous round, that is, the current computing unit to be updated can be updated by fully utilizing the developed information, and the updated flow field value is closer to the final solution, so that the fast convergence of steady flow is facilitated.
Referring to fig. 3, an embodiment of the present application discloses a specific steady flow numerical simulation method, which includes:
step S31: and updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence to obtain an updated computing unit and an updated current flow field in the iteration step of the round.
In this embodiment, the first computing unit is updatedWhen this cell is a boundary cell, the boundary condition can be used to obtain a virtual cell +.>And->And updating according to the flow field updating method; update last calculation Unit->In this case, the cell is also a boundary cell, and the boundary condition can be used to obtain a virtual cell +.>And->And updated according to the flow field updating method.
Step S32: and judging whether the current flow field updated in the iteration step of the round meets a preset convergence condition.
In this embodiment, the determining whether the updated current flow field in the iteration step of the present round meets a preset convergence condition, and if not, starting a next iteration step includes: and acquiring the residual error of the current flow field updated in the iteration step of the round, and judging whether the residual error meets a preset convergence condition. The calculation formula for obtaining the residual error of the current flow field is as follows:
wherein N represents the total number of flow field units, q i Representing the calculated value, q, of the calculation unit i exact Representing the exact value of the computational unit i, e.g. q exact ＝1。
Step S33: if the current flow field is not satisfied, starting a next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step in turn based on the target boundary condition, the target flow field unit updating function and the current flow field until the preset iteration stop condition is currently satisfied, and determining the current flow field which is currently output as the target flow field to obtain a corresponding flow field simulation result.
It can be understood that whether to stop the new iteration according to whether the updated flow field in the iteration step of the present round meets the preset convergence condition, i.e. whether to stop the current flow field updated in the iteration step of the present roundWhether the residual meets a preset convergence condition, e.g. the preset convergence condition is that the residual is less than 10 -6 The iteration may be stopped.
Therefore, in the iteration process, after each calculation unit completes the iteration, the numerical value of the calculation unit is replaced by the flow field value after the iteration from the flow field value of the previous iteration step, so that the updated flow field value can be used when other calculation units update the iteration. Compared with the standard explicit format, the method of the application utilizes the information of the updated unit, and can play a role in accelerating the convergence of the flow field.
The present application will be described with reference to a specific steady flow numerical simulation schematic diagram shown in fig. 4. If the one-dimensional linear convection problem is calculated by adopting a time first-order forward difference and a space first-order windward format, the formula of a control equation is as follows:
where q represents a scalar flow field value, t represents time, x represents a spatial coordinate, and a represents a propagation velocity.
1) Reading in a calculation grid, setting target boundary conditions, and uniformly distributing 51 points in a closed zone of [ -10, 10] as the calculation grid; the left side of the one-dimensional calculation domain is set as a velocity inflow boundary, the unit value of the inlet boundary is defined as 1, the right side of the calculation domain is set as an outlet boundary, and the unit value of the outlet boundary is always equal to the unit value of the adjacent inner point.
2) Initializing a one-dimensional flow field, wherein the flow field value of a left side interval [ -10,0] is 1, the flow field value of a right side interval (0, 10) is 2, and an initialization formula is as follows:
where q represents the spatial distribution function of the initial scalar field, q L Represents the flow field value of the left side, q R Represents the flow field value on the right side, x representsAnd (5) space coordinates.
3) The flow field unit for performing the iterative step of this round is updated by updating all the calculation units according to a preset updating sequence, for example, dispersing the control equation by adopting a space-time first-order formatCan be expressed in the following form:
wherein CFL is an abbreviation of the number Courant-Friedrichs-Lewy and has a value in the range of (0, 1]Taking 0.9 in this embodiment, the above formula is performed on all cells, e.g., in calculating the 10 th cellWhen the number 9 unit has completed updating, then->Substituting the known quantity into the above formula yields an updated value for the 10 th cell, as shown in the following formula:
updating all the computing units according to the above flow 3), judging whether all the computing units finish updating, if so, indicating that the iteration step of the round is finished.
4) After the iteration step of the round is completed, whether a preset iteration stopping condition is met currently is required to be judged, for example, whether the updated current flow field in the iteration step of the round meets a preset convergence condition is judged, namely, the residual error of the updated current flow field in the iteration step of the round is required to be counted, if the residual error meets the preset convergence condition, the updated current flow field in the iteration step of the round is judged to be a target flow field, the target flow field is saved, and the iteration is terminated; if the residual error does not meet the preset convergence condition, the process jumps back to the process 3), and the next iteration step is started until the residual error meets the preset convergence condition, so as to obtain the target flow field.
Fig. 5 is a schematic diagram of a specific residual convergence curve, where a standard explicit curve corresponds to a residual convergence curve based on a steady flow numerical simulation of a standard explicit format, and NewExplicit corresponds to a residual convergence curve based on a steady flow numerical simulation of a new explicit format of the present application, as can be seen from fig. 5, taking a calculation condition of cfl=0.9 as an example, the number of convergence steps of the standard explicit format is 55, that is, the number of convergence steps of the explicit calculation method of the present application needs to iterate 55 times to satisfy the preset convergence condition, that is, the number of convergence steps of the explicit calculation method of the present application needs to iterate only 30 times to satisfy the preset convergence condition, so that it can be seen that, when the preset convergence conditions are the same, the convergence speed of the explicit calculation method of the present application is approximately doubled compared with the convergence speed of the standard explicit format.
Referring to fig. 6, an embodiment of the present application discloses a steady flow numerical simulation device, including:
the flow field iteration module 11 is used for sequentially updating each non-updated calculation unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain an updated calculation unit in the iteration step of the round and the updated current flow field;
the simulation result obtaining module 12 is configured to, if the preset iteration stop condition is currently met, start a next iteration step to re-jump to the step of sequentially updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field until the preset iteration stop condition is currently met, and determine the current flow field currently output as the target flow field to obtain a corresponding flow field simulation result.
The beneficial effects of the application are that: sequentially updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain updated computing units and updated current flow fields in the iteration step of the round; judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result. Therefore, when the flow field is updated, each un-updated calculation unit in the iteration step of the round needs to be updated based on the current flow field, and when the calculation template where the unit to be updated is located contains the updated unit of the round, the updated flow field value of the unit is adopted to calculate, so that updated flow field information can be fully utilized, the target flow field can be acquired more quickly, and the effect of accelerating convergence is achieved.
Further, the embodiment of the application also provides electronic equipment. Fig. 7 is a block diagram of an electronic device 20, according to an exemplary embodiment, and the contents of the diagram should not be construed as limiting the scope of use of the present application in any way.
Fig. 7 is a schematic structural diagram of an electronic device according to an embodiment of the present application. Specifically, the method comprises the following steps: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input output interface 25, and a communication bus 26. The memory 22 is used for storing a computer program, and the computer program is loaded and executed by the processor 21 to implement relevant steps in the steady flow numerical simulation method executed by the electronic device as disclosed in any of the foregoing embodiments.
In this embodiment, the power supply 23 is configured to provide an operating voltage for each hardware device on the electronic device; the communication interface 24 can create a data transmission channel between the electronic device and an external device, and the communication protocol to be followed is any communication protocol applicable to the technical solution of the present application, which is not specifically limited herein; the input/output interface 25 is used for acquiring external input data or outputting external output data, and the specific interface type thereof may be selected according to the specific application requirement, which is not limited herein.
Processor 21 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 21 may be implemented in at least one hardware form of DSP (Digital Signal Processing ), FPGA (Field-Programmable Gate Array, field programmable gate array), PLA (Programmable Logic Array ). The processor 21 may also comprise a main processor, which is a processor for processing data in an awake state, also called CPU (Central Processing Unit ); a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 21 may integrate a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content required to be displayed by the display screen. In some embodiments, the processor 21 may also include an AI (Artificial Intelligence ) processor for processing computing operations related to machine learning.
The memory 22 may be a carrier for storing resources, such as a read-only memory, a random access memory, a magnetic disk, or an optical disk, and the resources stored thereon include an operating system 221, a computer program 222, and data 223, and the storage may be temporary storage or permanent storage.
The operating system 221 is used for managing and controlling various hardware devices on the electronic device and the computer program 222, so as to implement the operation and processing of the processor 21 on the mass data 223 in the memory 22, which may be Windows, unix, linux. The computer program 222 may further include a computer program that can be used to perform other specific tasks in addition to the computer program that can be used to perform the steady flow numerical simulation method performed by the electronic device as disclosed in any of the previous embodiments. The data 223 may include, in addition to data received by the electronic device and transmitted by the external device, data collected by the input/output interface 25 itself, and so on.
Further, the embodiment of the application also discloses a computer readable storage medium, wherein the storage medium stores a computer program, and when the computer program is loaded and executed by a processor, the method steps executed in the steady flow numerical simulation process disclosed in any embodiment are realized.
Finally, it is further noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
The foregoing has described in detail a method, apparatus, device and medium for steady flow numerical simulation provided by the present invention, and specific examples are applied herein to illustrate the principles and embodiments of the present invention, and the above examples are only used to help understand the method and core idea of the present invention; meanwhile, as those skilled in the art will have variations in the specific embodiments and application scope in accordance with the ideas of the present invention, the present description should not be construed as limiting the present invention in view of the above.

Claims (10)

1. A steady flow numerical simulation method, characterized by comprising:
sequentially updating each non-updated computing unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain updated computing units and updated current flow fields in the iteration step of the round;
judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step of the current iteration step based on the target boundary condition, the target flow field unit updating function and the current flow field in sequence, until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result.
2. The steady flow numerical simulation method according to claim 1, wherein before updating each un-updated calculation unit in the iteration step of the present round in sequence based on the target boundary condition, the target flow field unit updating function and the current flow field, the method further comprises:
setting a surface boundary condition of an aircraft to be simulated and a boundary condition of the periphery of a calculation domain to obtain a target boundary condition; the boundary conditions of the periphery of the calculation domain comprise far field conditions, inflow conditions, outflow conditions and symmetrical plane conditions.
3. The steady flow numerical simulation method according to claim 1, wherein the determining whether the preset iteration stop condition is satisfied at present, if not, starting the next iteration step, includes:
judging whether the current total iteration times are larger than a preset iteration threshold, and if not, starting the next iteration step.
4. A steady flow numerical simulation method according to any one of claims 1 to 3, characterized in that the determining whether the preset iteration stop condition is currently met or not, and if not, starting a next iteration step, comprises:
judging whether the updated current flow field in the iteration step of the round meets a preset convergence condition, and if not, starting the next iteration step.
5. The method for modeling stationary flow values according to claim 4, wherein said determining whether the current flow field updated in the current iteration step satisfies a preset convergence condition, if not, starting a next iteration step, comprises:
acquiring the residual error of the current flow field updated in the iteration step of the round, and judging whether the residual error meets a preset convergence condition;
if not, starting the next iteration step.
6. The steady flow numerical simulation method according to claim 1, wherein the sequentially updating each non-updated calculation unit in the iteration step of the present round based on the target boundary condition, the target flow field unit updating function and the current flow field to obtain an updated calculation unit in the iteration step of the present round and the updated current flow field comprises:
judging whether updated computing units and non-updated computing units exist in the iteration steps of the round, if so, determining the current computing unit to be updated in the iteration steps of the round from the non-updated computing units of the iteration steps of the round;
updating the current computing unit to be updated in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field values of all computing units in the peripheral computing template range of the current computing unit to be updated in the iteration step of the round, and then re-jumping to the step of judging whether updated computing units and non-updated computing units exist in the iteration step of the round or not until all computing units in the iteration step of the round are updated.
7. The steady flow numerical simulation method according to claim 6, wherein before updating the current calculation unit to be updated in the present iteration step based on the target boundary condition, the target flow field unit update function, and the flow field value of the updated calculation unit, further comprises:
and determining an initial flow field value corresponding to each calculation unit in the initial flow field, and determining a current calculation unit to be updated of an initial iteration step from each calculation unit.
8. A steady flow numerical simulation device, comprising:
the flow field iteration module is used for sequentially updating each non-updated calculation unit in the iteration step of the round based on the target boundary condition, the target flow field unit updating function and the current flow field so as to obtain an updated calculation unit in the iteration step of the round and the updated current flow field;
the simulation result obtaining module is used for judging whether the preset iteration stop condition is met currently, if not, starting the next iteration step to jump to the step of updating each un-updated calculation unit in the iteration step in turn based on the target boundary condition, the target flow field unit updating function and the current flow field until the preset iteration stop condition is met currently, and determining the current flow field which is output currently as the target flow field to obtain a corresponding flow field simulation result.
9. An electronic device, comprising:
a memory for storing a computer program;
a processor for executing the computer program to implement the steps of the steady flow numerical simulation method as claimed in any one of claims 1 to 7.
10. A computer-readable storage medium storing a computer program; wherein the computer program when executed by a processor implements the steps of the steady flow numerical simulation method of any of claims 1 to 7.
CN202310591185.7A 2023-05-24 2023-05-24 Constant flow numerical simulation method, device, equipment and medium Pending CN116542178A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310591185.7A CN116542178A (en) 2023-05-24 2023-05-24 Constant flow numerical simulation method, device, equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310591185.7A CN116542178A (en) 2023-05-24 2023-05-24 Constant flow numerical simulation method, device, equipment and medium

Publications (1)

Publication Number Publication Date
CN116542178A true CN116542178A (en) 2023-08-04

Family

ID=87446851

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310591185.7A Pending CN116542178A (en) 2023-05-24 2023-05-24 Constant flow numerical simulation method, device, equipment and medium

Country Status (1)

Country Link
CN (1) CN116542178A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116738892A (en) * 2023-08-08 2023-09-12 中国空气动力研究与发展中心计算空气动力研究所 Flow field numerical simulation data protection method, system, equipment and medium

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116738892A (en) * 2023-08-08 2023-09-12 中国空气动力研究与发展中心计算空气动力研究所 Flow field numerical simulation data protection method, system, equipment and medium
CN116738892B (en) * 2023-08-08 2023-10-20 中国空气动力研究与发展中心计算空气动力研究所 Flow field numerical simulation data protection method, system, equipment and medium

Similar Documents

Publication Publication Date Title
CN114862656B (en) Multi-GPU-based acquisition method for training cost of distributed deep learning model
JP7350923B2 (en) Deep learning framework operator registration method, equipment, device and storage medium
CN116542178A (en) Constant flow numerical simulation method, device, equipment and medium
CN114968559B (en) LSF-based multi-host multi-GPU distributed arrangement deep learning model method
CN112527281B (en) Operator upgrading method and device based on artificial intelligence, electronic equipment and medium
CN115033391B (en) Data flow method and device for neural network calculation
JP7435951B2 (en) Floating point number generation method, apparatus, electronic device, storage medium and computer program for integrated circuit chip verification
WO2018152319A1 (en) Spreadsheet recalculation algorithm for directed acyclic graph processing
CN112749809A (en) Method and device for constructing quantum simulation system
CN116933676A (en) Flow numerical simulation method, device, equipment and medium
CN117077294B (en) Parameterization-based structural grid generation method, apparatus, device and medium
US20200272901A1 (en) Optimization apparatus, optimization method, and non-transitory computer readable medium
CN116167445B (en) Quantum measurement mode processing method and device and electronic equipment
WO2023216500A1 (en) Computing power resource deployment method and apparatus for intelligent computing center, and device and storage medium
CN105653355A (en) Method and system for calculating Hadoop configuration parameters
EP3968202A1 (en) Customizable reinforcement of learning column placement in structural design
US10467120B2 (en) Software optimization for multicore systems
CN109492086A (en) A kind of answer output method, device, electronic equipment and storage medium
US11941545B2 (en) Systems and methods for generating a boundary of a footprint of uncertainty for an interval type-2 membership function based on a transformation of another boundary
CN114492251A (en) Low-speed flow field divergence processing method, device, equipment and medium in supercomputing environment
CN110378037B (en) CFD simulation data storage method and device based on Ceph and server
CN113570067A (en) Synchronization method, device and program product of distributed system
CN117827619B (en) Time-consuming prediction simulation method, device, equipment, medium and system for heterogeneous calculation force
KR102457154B1 (en) Method and system for generating intermediate representation for program for execution on accelerator
EP4216097A1 (en) Method, system, equipment and medium for modifying the layering layer information of finite element model unit

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