CN116070071A - Wall function selection method, device, terminal equipment and storage medium - Google Patents
Wall function selection method, device, terminal equipment and storage medium Download PDFInfo
- Publication number
- CN116070071A CN116070071A CN202310347874.3A CN202310347874A CN116070071A CN 116070071 A CN116070071 A CN 116070071A CN 202310347874 A CN202310347874 A CN 202310347874A CN 116070071 A CN116070071 A CN 116070071A
- Authority
- CN
- China
- Prior art keywords
- flow type
- wall function
- object plane
- flow
- compressible
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- 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
Abstract
The application discloses a method, a device, a terminal device and a storage medium for selecting a wall function, which are based on the followingy + Standard wall function and basis of (2)y * Is based on the fact that the compressible correction wall function of the system is simulated and the object plane normal partial derivative of the tangential velocity of the object plane at grid points in the boundary layer near the object plane is calculatedy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined; according to the flow type of the turbulent flow, the target wall function corresponding to the flow type is selected, and through the embodiment of the invention, the calculation result can be more accurate, the embodiment of the invention can automatically judge the flow type when the compressible flow is simulated by using the wall function, and the more adaptive wall function type is selected, so that the more accurate simulation result is obtained.
Description
Technical Field
The application belongs to the technical field of hydrodynamics, and particularly relates to a method and a device for selecting a wall function, terminal equipment and a storage medium.
Background
The RANS (reynolds average equation) method is designed to reduce the dependence of the numerical solution (in particular frictional resistance, heat flow, etc.) on the location of the first layer mesh nodes adjacent to the wall by providing a boundary condition at the solid wall when simulating turbulent flow. The turbulent boundary layer simulation adopts a wall function method, so that the scale of a first layer of grid of the wall surface can be greatly widened, the numerical rigidity problem can be reduced, the convergence effect in the calculation process can be accelerated, and the stability of iterative calculation can be improved.
At present, four common and widely applied wall functions are available, one is based ony + The standard wall function of (2), based on the teaching of Laundery * Third, nichols introduces a compressible logarithmic form of White and Christoph based on a uniform effective velocity wall form to obtain a uniform effective, compressible and thermally conductive-based wall-approaching layery + Four is based on Nichols's continuous wall lawy * Is a continuous compressible wall law. In the actual numerical simulation process, since different compressible flow calculation examples use different wall functions to simulate effects, how to select proper wall functions can be calculated and accurately is a problem which needs to be solved at present.
Disclosure of Invention
The invention aims to provide a method, a device, a terminal device and a storage medium for selecting a wall function, so as to solve the defects in the prior art.
In a first aspect, an embodiment of the present invention provides a method for selecting a wall function, where the method includes:
respectively based ony + Standard wall function and basis of (2)y * Performing simulation and calculation on the compressible correction wall function to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane;
according to the basisy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined;
and selecting a target wall function corresponding to the flow type according to the flow type of the turbulent flow.
Optionally, the method is based ony + Andy * determining the flow type of the turbulent flow, comprising:
if based ony + Andy * if the simulation results of the turbulence are inconsistent, the flow type of the turbulence is determined to be a highly compressible flow type.
Optionally, the selecting a target wall function corresponding to the flow type according to the flow type of the turbulent flow includes:
if the flow type of the turbulent flow is a flow type with strong compressibility, calculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane;
if the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with strong compressibility is a separation flow type with strong compressibility;
and selecting a target wall function corresponding to the flow type according to whether the flow type with strong compressibility is a separation flow type.
Optionally, the selecting a target wall function corresponding to the flow type according to whether the flow type with strong compressibility is a separation flow type includes:
if the compressibility is strongThe flow type is a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type;
if the highly compressible flow type is not a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type.
Optionally, the selecting a target wall function corresponding to the flow type according to the flow type of the turbulent flow includes:
if the flow type of the turbulent flow is a flow type with weak compressibility, calculating an object plane normal partial derivative of the tangential velocity of an object plane on grid points in a boundary layer near the object plane;
if the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with low compressibility is a separation flow type with low compressibility;
and selecting a target wall function corresponding to the flow type according to whether the flow type with low compressibility is a separation flow type.
Optionally, the selecting a target wall function corresponding to the flow type according to whether the flow type with low compressibility is a separation flow type includes:
if the non-compressible flow type is a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type;
if the non-compressible flow type is not a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type.
In a second aspect, an embodiment of the present invention provides a device for selecting a wall function, where the device includes:
acquisition modules for respectively using the basey + Standard wall function and basis of (2)y * Performing simulation and calculation on the compressible correction wall function to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane;
a determining module for according to the basey + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined;
and the selection module is used for determining a target wall function corresponding to the flow type according to the flow type of the turbulent flow.
Optionally, the determining module is configured to:
if based ony + Andy * if the simulation results of the turbulence are inconsistent, the flow type of the turbulence is determined to be a highly compressible flow type.
If the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with strong compressibility is a separation flow type with strong compressibility;
optionally, the selecting module is configured to:
and selecting a target wall function corresponding to the flow type according to whether the flow type with low compressibility is a separation flow type.
If the highly compressible flow type is a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type;
if the highly compressible flow type is not a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type.
Optionally, the determining module is configured to:
if based ony + Andy * if the simulation results of the turbulence are consistent, the flow type of the turbulence is determined to be the flow type with weak compressibility.
If the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with low compressibility is a separation flow type with low compressibility;
optionally, the selecting module is configured to:
and selecting a target wall function corresponding to the flow type according to whether the flow type with low compressibility is a separation flow type.
If the non-compressible flow type is a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type;
if the non-compressible flow type is not a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type.
In a third aspect, an embodiment of the present invention provides a terminal device, including: at least one processor and memory;
the memory stores a computer program; the at least one processor executes the computer program stored in the memory to implement the method of selecting a wall function provided in the first aspect.
In a fourth aspect, an embodiment of the present invention provides a computer readable storage medium having stored therein a computer program which, when executed, implements the method for selecting a wall function provided in the first aspect.
The embodiment of the invention has the following advantages:
the method, the device, the terminal equipment and the storage medium for selecting the wall function provided by the embodiment of the invention are based on the following steps of respectively usingy + Standard wall function and basis of (2)y * Performing simulation and calculation on the compressible correction wall function to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane; according to the basis ofy + Andy * consistency of simulation results of (c) and boundary layer inner net near object planeDetermining the flow type of turbulent flow by the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point; according to the flow type of turbulent flow, a target wall function corresponding to the flow type is selected, and the flow is based on the flow typey + Andy * the simulation result consistency of the flow type and the object plane normal partial derivative of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane are judged to determine the flow type, and then the target wall function corresponding to the flow type is selected, so that the calculation result is more accurate.
Drawings
In order to more clearly illustrate the embodiments or prior art solutions of the present application, the drawings that are required for the description of the embodiments or prior art will be briefly described below, it being apparent that the drawings in the following description are only some of the embodiments described in the present application, and that other drawings may be obtained according to these drawings without inventive faculty for a person skilled in the art.
FIG. 1 is a flow chart of a method for selecting a wall function according to an embodiment of the present application;
FIG. 2 is a flow chart of a method of selecting a wall function according to an embodiment of the present application;
FIG. 3 is a graph showing the results of modeling a compressible flat turbulent boundary layer using different wall functions in one embodiment of the present application;
FIG. 4 is a graph showing the results of simulating a supersonic groove separation flow using different wall functions in one embodiment of the present application;
FIG. 5 is a block diagram of an embodiment of a wall function selection device of the present invention;
fig. 6 is a schematic structural diagram of a terminal device of the present invention.
Detailed Description
For the purposes, technical solutions and advantages of the present application, the technical solutions of the present application will be clearly and completely described below with reference to specific embodiments and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
An embodiment of the invention provides a method for selecting a wall function, which is used for determining a proper wall function. The selecting device of the embodiment, in which the execution body is a wall function, is provided on a terminal device, for example, the terminal device includes at least a computer terminal and the like.
Referring to fig. 1, there is shown a flow chart of steps of an embodiment of a method of selecting a wall function according to the present invention, which may include the steps of:
s101, respectively using the basey + Standard wall function and basis of (2)y * Performing simulation and calculation on the compressible correction wall function to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane;
specifically, in the embodiment of the present invention, terminal devices respectively use basey + Standard wall function and basis of (2)y * And (3) carrying out simulation and calculation on the compressible correction wall function to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane.
S102, according to the basisy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined;
specifically, the terminal device pair is based ony + Andy * and (3) judging the consistency of simulation results and the normal partial derivative of the object plane of tangential velocity of the object plane on grid points in the boundary layer near the object plane, and determining the flow type according to the judgment result.
As an alternative embodiment, if based ony + Andy * if the simulation results of (a) are inconsistent, determining the flow type of the turbulent flow as being possibleA highly compressible flow type.
As another alternative, if based ony + Andy * if the simulation results are consistent, the turbulent flow type is determined to be a flow type with weak compressibility.
S103, selecting a target wall function corresponding to the flow type according to the flow type of the turbulent flow.
Specifically, according to different flow types, the terminal equipment selects a target wall function corresponding to the flow type as an optional implementation mode, and if the flow type of the turbulent flow is a flow type with strong compressibility, the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane is calculated;
if the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with strong compressibility is a separation flow type with strong compressibility;
according to whether the flow type with strong compressibility is a separation flow type, selecting a target wall function corresponding to the flow type.
Specifically, if the highly compressible flow type is a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type;
if the highly compressible flow type is not a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type.
As another alternative, if the turbulent flow type is not compressible, calculating the object plane normal partial derivative of the object plane tangential velocity at grid points in the boundary layer near the object plane;
if the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with low compressibility is a separation flow type with low compressibility;
according to whether the flow type with low compressibility is a separation flow type, selecting a target wall function corresponding to the flow type.
Specifically, if the non-compressible flow type is a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type;
if the non-compressible flow type is not a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type.
The wall function selection method provided by the embodiment of the invention is respectively based ony + Standard wall function and basis of (2)y * Performing simulation and calculation to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane, and determining the flow type of the turbulence; according to the flow type of turbulent flow, a target wall function corresponding to the flow type is selected, and the flow is based on the flow typey + Andy * the simulation result consistency of the flow type and the object plane normal partial derivative of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane are judged to determine the flow type, and then the target wall function corresponding to the flow type is selected, so that the calculation result is more accurate.
In a further embodiment of the present invention, a method for selecting a wall function provided in the above embodiment is further described.
FIG. 2 is a flow chart of a method of selecting a wall function according to an embodiment of the present application; as shown in fig. 2, the embodiment of the invention provides a method for selecting a wall function in compressible flow, which can automatically judge the flow type when the compressible flow is simulated by using the wall function, and a more suitable wall function type is selected to obtain a more accurate simulation result.
In order to achieve the technical purpose, the invention adopts the following technical scheme.
Four wall functions are implemented in the simulation program code to form independent wall function modules.
A wall function selecting module is realized in the simulation program code, and a plurality of monitoring points can be arranged in the flow field in the selecting module for judging the flow type and selecting the corresponding wall function.
By utilizing the method, the wall function capable of realizing compressible flow is automatically selected, and the method sequentially comprises the following steps of:
(1) Preparing a simulation example, namely a program with a wall function and a selected module, a grid and a parameter file;
(2) Starting calculation and automatically running a wall function selecting module;
(3) The wall function selecting modules firstly respectively use the basey + Standard wall function and basis of (2)y * Is simulated and calculated to obtain all grid points in boundary layer near object plane
Values.
Wherein u is the flow direction speed,yin the normal direction of the wall, ∂ is the sign of the partial derivative,
the partial derivative of the normal direction of the flow velocity at the wall surface;
(4) If based ony + Andy * if the simulation results of the two are inconsistent, judging that the flow with strong compressibility is the flow with weak compressibility, otherwise, judging that the flow with strong compressibility is the flow with weak compressibility;
(5) If it is determined that the flow with high compressibility is detected, then if there are grid points
Judging that the flow with strong compressibility is the separation flow with strong compressibility, otherwise, judging that the flow with strong compressibility is not the separation flow with strong compressibility;
(6) If the compressibility is strongThe flow is a highly compressible split flow, then automatic enablement is based ony * Is a continuous compressible wall law;
(7) If the above-described highly compressible flow is not a highly compressible separation flow, then the automatic enablement is based ony * Is used for compressing and correcting the wall function;
(8) If it is determined that the flow with low compressibility is not strong, then if there are grid points
Judging that the flow with the incompressibility is the separation flow with the incompressibility, otherwise, judging that the flow with the incompressibility is not the separation flow with the incompressibility;
(9) If the non-compressible flow is a non-compressible separation flow, the automatic activation is based ony + Is a continuous compressible wall law;
(10) If the non-compressible flow is not a non-compressible separation flow, the automatic activation is based ony + Is a standard wall function of (2);
(11) And carrying out subsequent calculation processing according to the selected target wall surface function.
Compared with the prior art, the embodiment of the invention provides a method for selecting the wall function in the compressible flow, and a more accurate simulation result can be obtained when the wall function is used for simulating the compressible flow.
In the embodiment of the invention, specifically:
(1) Preparing a simulation example, namely a program with a wall function and a selected module, a grid and a parameter file;
(2) Starting calculation and automatically running a wall function selecting module;
(3) Respectively based ony + Standard wall function and basis of (2)y * Is based on the simulation of the compressible correction wall functiony + Andy * the simulation results of (2) are consistent, and the flow with weak compressibility is judged;
(4) Calculating grid points in boundary layers near all object planes
A value;
(5) With grid points
Further, it is determined that the separation flow is not strong in compressibility;
(6) Automatic enablement based ony + Is a continuous compressible wall law.
(7) And continuing the subsequent calculation processing.
FIG. 3 is a graph showing the results of modeling a compressible flat turbulent boundary layer using different wall functions in one embodiment of the present application; indicating use is based ony * The compressible modified wall function (labeled wf 1) better simulates the attached turbulent boundary layer flow without separation.
FIG. 4 is a graph showing the results of simulating a supersonic groove separation flow using different wall functions in one embodiment of the present application; indicating that the use of a continuously compressible wall law, labeled wf6, better simulates the split flow.
Letter and corresponding name in the drawings:C f as a coefficient of friction, the coefficient of friction,xas the coordinate axis, wf0 is the unopened wall function.
It should be noted that, for simplicity of description, the method embodiments are shown as a series of acts, but it should be understood by those skilled in the art that the embodiments are not limited by the order of acts, as some steps may occur in other orders or concurrently in accordance with the embodiments. Further, those skilled in the art will appreciate that the embodiments described in the specification are presently preferred embodiments, and that the acts are not necessarily required by the embodiments of the invention.
The wall function selection method provided by the embodiment of the invention is respectively based ony + Standard wall function and basis of (2)y * Is simulated by the compressible correction wall functionCalculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane; according to the basis ofy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined; according to the flow type of turbulent flow, a target wall function corresponding to the flow type is selected, and the flow is based on the flow typey + Andy * the simulation result consistency of the flow type and the object plane normal partial derivative of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane are judged to determine the flow type, and then the target wall function corresponding to the flow type is selected, so that the calculation result is more accurate.
Another embodiment of the present invention provides a device for selecting a wall function, configured to execute the method for selecting a wall function provided in the foregoing embodiment.
Referring to fig. 5, there is shown a block diagram of an embodiment of a wall function selection device of the present invention, which may include the following modules: an acquisition module 501, a determination module 502, and a selection module 503, wherein:
the acquisition module 501 is configured to use the base respectivelyy + Standard wall function and basis of (2)y * Performing simulation and calculation on the compressible correction wall function to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane;
the determining module 502 is used for determining that the data is based ony + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined;
the selection module 503 is configured to select a target wall function corresponding to a flow type according to the flow type of the turbulent flow.
The wall function selecting device provided by the embodiment of the invention is respectively realized byUsage based ony + Standard wall function and basis of (2)y * Performing simulation on the compressible modified wall function and calculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane; according to the basis ofy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined; according to the flow type of turbulent flow, a target wall function corresponding to the flow type is selected, and the flow is based on the flow typey + Andy * the simulation result consistency of the flow type and the object plane normal partial derivative of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane are judged to determine the flow type, and then the target wall function corresponding to the flow type is selected, so that the calculation result is more accurate.
In a further embodiment of the present invention, the wall function selecting device provided in the above embodiment is further described in addition.
Optionally, the determining module is configured to:
according to the basis ofy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined;
and the selection module is used for determining a target wall function corresponding to the flow type according to the flow type of the turbulent flow.
Optionally, the determining module is configured to:
if based ony + Andy * if the simulation results of the turbulence are inconsistent, the flow type of the turbulence is determined to be a highly compressible flow type.
If the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with strong compressibility is a separation flow type with strong compressibility;
optionally, the selecting module is configured to:
and selecting a target wall function corresponding to the flow type according to whether the flow type with strong compressibility is a separation flow type.
If the highly compressible flow type is a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type;
if the highly compressible flow type is not a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type.
Optionally, the determining module is configured to:
if based ony + Andy * if the simulation results of the turbulence are consistent, the flow type of the turbulence is determined to be the flow type with weak compressibility.
If the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with low compressibility is a separation flow type with low compressibility;
optionally, the selecting module is configured to:
and selecting a target wall function corresponding to the flow type according to whether the flow type with low compressibility is a separation flow type.
If the non-compressible flow type is a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type;
if the non-compressible flow type is not a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type.
For the device embodiments, since they are substantially similar to the method embodiments, the description is relatively simple, and reference is made to the description of the method embodiments for relevant points.
The wall function selecting device provided by the embodiment of the invention is based on the following modesy + Standard wall function and basis of (2)y * Performing simulation on the compressible modified wall function and calculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane; according to the basis ofy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined; according to the flow type of turbulent flow, a target wall function corresponding to the flow type is selected, and the flow is based on the flow typey + Andy * the simulation result consistency of the flow type and the object plane normal partial derivative of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane are judged to determine the flow type, and then the target wall function corresponding to the flow type is selected, so that the calculation result is more accurate.
A further embodiment of the present invention provides a terminal device, configured to execute the method for selecting a wall function provided in the foregoing embodiment.
Fig. 6 is a schematic structural view of a terminal device of the present invention, as shown in fig. 6, the terminal device includes: at least one processor 601 and memory 602;
the memory stores a computer program; the at least one processor executes the computer program stored in the memory to implement the method for selecting a wall function provided in the above embodiment.
The terminal device provided in this embodiment is based on the following by using respectivelyy + Standard wall function and basis of (2)y * Performing simulation on the compressible modified wall function and calculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane; according to the basis ofy + Andy * is a simulated knot of (2)Determining the flow type of turbulent flow according to the consistency of the result and the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane; according to the flow type of turbulent flow, a target wall function corresponding to the flow type is selected, and the flow is based on the flow typey + Andy * the simulation result consistency of the flow type and the object plane normal partial derivative of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane are judged to determine the flow type, and then the target wall function corresponding to the flow type is selected, so that the proper wall function can be selected, the calculation result is more accurate, the flow type can be automatically judged when the compressible flow is simulated by using the wall function, and the more adaptive wall function type is selected to obtain the more accurate simulation result.
A further embodiment of the present application provides a computer readable storage medium having a computer program stored therein, which when executed implements the method for selecting a wall function provided in any of the above embodiments.
According to the computer-readable storage medium of the present embodiment, by using the base stations respectivelyy + Standard wall function and basis of (2)y * Performing simulation on the compressible modified wall function and calculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane; according to the basis ofy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined; according to the flow type of turbulent flow, a target wall function corresponding to the flow type is selected, and the flow is based on the flow typey + Andy * the consistency of the simulation result of the flow type and the normal partial derivative of the object plane on the grid point in the boundary layer near the object plane are judged to determine the flow type, and then the target wall function corresponding to the flow type is selected, so that the calculation result is more accurate, the embodiment of the invention can automatically judge the flow type when the wall function is used for simulating compressible flow, and a more adaptive wall function is selectedThe number type, obtain more accurate simulation result.
It should be noted that the foregoing detailed description is exemplary and is intended to provide further explanation of the invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments in accordance with the present application. As used herein, the singular is intended to include the plural unless the context clearly indicates otherwise. Furthermore, it will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, steps, operations, devices, components, and/or groups thereof.
It should be noted that the terms "first," "second," and the like in the description and claims of the present application and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the terms so used are interchangeable under appropriate circumstances such that the embodiments of the application described herein are capable of operation in sequences other than those illustrated or otherwise described herein.
Furthermore, the terms "comprise" and "have," as well as any variations thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those elements but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
Spatially relative terms, such as "above … …," "above … …," "upper surface at … …," "above," and the like, may be used herein for ease of description to describe one device or feature's spatial location relative to another device or feature as illustrated in the figures. It will be understood that the spatially relative terms are intended to encompass different orientations in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "above" or "over" other devices or structures would then be oriented "below" or "beneath" the other devices or structures. Thus, the exemplary term "above … …" may include both orientations of "above … …" and "below … …". The device may also be positioned in other different ways, such as rotated 90 degrees or at other orientations, and the spatially relative descriptors used herein interpreted accordingly.
In the above detailed description, reference is made to the accompanying drawings, which form a part hereof. In the drawings, like numerals typically identify like components unless context indicates otherwise. The illustrated embodiments described in the detailed description, drawings, and claims are not meant to be limiting. Other embodiments may be utilized, and other changes may be made, without departing from the spirit or scope of the subject matter presented herein.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A method of selecting a wall function, the method comprising:
respectively based ony + Standard wall function and basis of (2)y * Performing simulation on the compressible modified wall function and calculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane;
according to the basisy + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined;
and selecting a target wall function corresponding to the flow type according to the flow type of the turbulent flow.
2. The method according to claim 1, wherein said step of based ony + Andy * determining the flow type of the turbulent flow, comprising:
if based ony + Andy * if the simulation results of the turbulence are inconsistent, the flow type of the turbulence is determined to be a highly compressible flow type.
3. The method of claim 2, wherein selecting a target wall function corresponding to the flow type based on the flow type of the turbulent flow comprises:
if the flow type of the turbulent flow is a flow type with strong compressibility, calculating the normal partial derivative of the object plane of the tangential velocity of the object plane on grid points in the boundary layer near the object plane;
if the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with strong compressibility is a separation flow type with strong compressibility;
and selecting a target wall function corresponding to the flow type according to whether the flow type with strong compressibility is a separation flow type.
4. A method according to claim 3, wherein said selecting a target wall function corresponding to said flow type based on whether said highly compressible flow type is a split flow type comprises:
if the highly compressible flow type is a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type;
if the highly compressible flow type is not a highly compressible separation flow type, the selection is based ony * Is a target wall function corresponding to the flow type.
5. The method according to claim 1, wherein said step of based ony + Andy * determining the flow type of the turbulent flow, comprising:
if based ony + Andy * if the simulation results of the turbulence are consistent, the flow type of the turbulence is determined to be the flow type with weak compressibility.
6. The method of claim 5, wherein selecting a target wall function corresponding to the flow type based on the flow type of the turbulent flow comprises:
if the flow type of the turbulent flow is a flow type with weak compressibility, calculating an object plane normal partial derivative of the tangential velocity of an object plane on grid points in a boundary layer near the object plane;
if the partial derivative of at least one grid point is smaller than or equal to a preset value, determining that the flow type with low compressibility is a separation flow type with low compressibility;
and selecting a target wall function corresponding to the flow type according to whether the flow type with low compressibility is a separation flow type.
7. The method of claim 6, wherein selecting the target wall function corresponding to the flow type based on whether the flow type that is not highly compressible is a split flow type comprises:
if the non-compressible flow type is a non-compressible separation flow type, the selection is based ony + Is a target wall function corresponding to the flow type;
if the non-compressible flow type is not a separate flow type that is not compressible,the selection is based ony + Is a target wall function corresponding to the flow type.
8. A wall function selection device, the device comprising:
acquisition modules for respectively using the basey + Standard wall function and basis of (2)y * Performing simulation and calculation on the compressible correction wall function to obtain the normal partial derivative of the object plane of the tangential velocity of the object plane on the grid point in the boundary layer near the object plane;
a determining module for determining the basis of the basey + Andy * the consistency of simulation results of the flow meter and the tangential velocity of the object plane on grid points in the boundary layer near the object plane are determined;
and the selection module is used for selecting a target wall surface function corresponding to the flow type according to the flow type of the turbulent flow.
9. A terminal device, comprising: at least one processor and memory;
the memory stores a computer program; the at least one processor executes the computer program stored by the memory to implement the method of selecting a wall function of any one of claims 1-7.
10. A computer readable storage medium, characterized in that the computer readable storage medium has stored therein a computer program which, when executed, implements the method of selecting a wall function according to any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310347874.3A CN116070071B (en) | 2023-04-04 | 2023-04-04 | Wall function selection method, device, terminal equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310347874.3A CN116070071B (en) | 2023-04-04 | 2023-04-04 | Wall function selection method, device, terminal equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116070071A true CN116070071A (en) | 2023-05-05 |
| CN116070071B CN116070071B (en) | 2023-06-27 |
Family
ID=86173513
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310347874.3A Active CN116070071B (en) | 2023-04-04 | 2023-04-04 | Wall function selection method, device, terminal equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116070071B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103218507A (en) * | 2012-12-13 | 2013-07-24 | 中国电力科学研究院 | Two-dimensional numerical simulation method for icing process of power transmission line |
| CN106682398A (en) * | 2016-12-12 | 2017-05-17 | 华南理工大学 | RANS/LES method based on tubulence energy and turbulence length scale truncation |
| CN106777823A (en) * | 2017-01-25 | 2017-05-31 | 天津大学 | Construction schemes of underground structure group progress simulation optimization method based on ventilation numerical simulation |
| CN110276090A (en) * | 2019-04-09 | 2019-09-24 | 广州市香港科大霍英东研究院 | A kind of eddy simulation of turbulent flows method based on coherent structure |
| CN110489709A (en) * | 2019-08-01 | 2019-11-22 | 中国空气动力研究与发展中心计算空气动力研究所 | The method for numerical simulation of parsing Wall-function based on compressible flows |
| CN112052512A (en) * | 2020-07-23 | 2020-12-08 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for judging layering of turbulent boundary layer |
| CN113158338A (en) * | 2021-04-13 | 2021-07-23 | 中国空气动力研究与发展中心计算空气动力研究所 | Rapid turbulence wall function aerodynamic force prediction method based on coarse grid |
| CN115862779A (en) * | 2022-11-25 | 2023-03-28 | 兰州大学 | Liquid metal conjugate heat transfer calculation method based on two-equation thermal turbulence model |
-
2023
- 2023-04-04 CN CN202310347874.3A patent/CN116070071B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103218507A (en) * | 2012-12-13 | 2013-07-24 | 中国电力科学研究院 | Two-dimensional numerical simulation method for icing process of power transmission line |
| CN106682398A (en) * | 2016-12-12 | 2017-05-17 | 华南理工大学 | RANS/LES method based on tubulence energy and turbulence length scale truncation |
| CN106777823A (en) * | 2017-01-25 | 2017-05-31 | 天津大学 | Construction schemes of underground structure group progress simulation optimization method based on ventilation numerical simulation |
| CN110276090A (en) * | 2019-04-09 | 2019-09-24 | 广州市香港科大霍英东研究院 | A kind of eddy simulation of turbulent flows method based on coherent structure |
| CN110489709A (en) * | 2019-08-01 | 2019-11-22 | 中国空气动力研究与发展中心计算空气动力研究所 | The method for numerical simulation of parsing Wall-function based on compressible flows |
| CN112052512A (en) * | 2020-07-23 | 2020-12-08 | 中国空气动力研究与发展中心计算空气动力研究所 | Method for judging layering of turbulent boundary layer |
| CN113158338A (en) * | 2021-04-13 | 2021-07-23 | 中国空气动力研究与发展中心计算空气动力研究所 | Rapid turbulence wall function aerodynamic force prediction method based on coarse grid |
| CN115862779A (en) * | 2022-11-25 | 2023-03-28 | 兰州大学 | Liquid metal conjugate heat transfer calculation method based on two-equation thermal turbulence model |
Non-Patent Citations (7)
| Title |
|---|
| ABDUSSAMET SUBASI 等: "Performance evaluation of RANS-based turbulence models in simulating a honeycomb heat sink", 《HEAT AND MASS TRANSFER》, pages 2435 - 2443 * |
| L.PATRUNO: "Accuracy of numerically evaluated flutter derivatives of bridge deck sections using RANS: Effects on the flutter onset velocity", 《ENGINEERING STRUCTURES》, pages 49 - 65 * |
| SERGEI V. UTYUZHNIKOV: "Generalized Wall-functions for High-Reynolds-number Turbulence Models", 《ARXIV》, pages 1 - 21 * |
| 伊伊爸: "Fluent湍流模型中壁面函数选择建议", Retrieved from the Internet <URL:360doc.com/content/21/0918/13/15913066_996114587.shtml> * |
| 毛枚良 等: "可压缩湍流边界层壁面函数方法综述", 《空气动力学学报》, pages 1 - 11 * |
| 王新光 等: "壁面函数在超声速湍流模拟中的应用", 《航空学报》, pages 126153 - 1 * |
| 赵会超 等: "CFD 用于流量测量辅助研究中的仿真精度", 《清华大学学报》, pages 1052 - 1056 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116070071B (en) | 2023-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108647668A (en) | The construction method of multiple dimensioned lightweight Face datection model and the method for detecting human face based on the model | |
| CN104217433B (en) | A kind of method and device for analyzing image | |
| CN107767375A (en) | One kind point cloud method for evaluating quality and device | |
| CN111428918B (en) | Soil pollution range prediction method and system for heavy metal attenuation pollution source | |
| CN113792478B (en) | Heat source layout temperature field reliability analysis method based on integrated depth proxy model | |
| US20100145668A1 (en) | Method for dynamic repartitioning in adaptive mesh processing | |
| US20150073730A1 (en) | Mechanical strain gauge simulation | |
| CN106556877A (en) | A kind of earth magnetism Tonghua method and device | |
| CN116070071B (en) | Wall function selection method, device, terminal equipment and storage medium | |
| CN107016218A (en) | A kind of method and apparatus for determining that finite element point load is distributed in airplane wingtip winglet aerofoil | |
| CN107544098B (en) | Surface roughness generation method and device, storage medium and processor | |
| CN108280265A (en) | Distortion of the mesh interpolation method, system and medium based on novel radial basis function | |
| CN105260249B (en) | Space computational domain calculates strength characteristic extracting method | |
| CN116227045B (en) | Local stress strain field construction method and system for structural test piece | |
| CN115272594A (en) | Iso-surface generation method based on geotools | |
| CN109684649A (en) | A kind of wind speed revision method and system based on landform | |
| CN116702571B (en) | Numerical simulation method and device based on multiple smoothness measurement factors | |
| CN110334459B (en) | Rapid and refined modeling system and method for power transmission tower line system | |
| CN108564308B (en) | Method and device for evaluating total radiation change characteristics of photovoltaic power station | |
| CN113496097A (en) | Aircraft fuel tank fuel oil shaking simulation analysis method based on SPH | |
| CN112016956A (en) | BP neural network-based ore grade estimation method and device | |
| CN111160530A (en) | Compression processing method and device of model and computer equipment | |
| CN110174865A (en) | A kind of method and device controlling target information acquisition rate | |
| CN108510778A (en) | A kind of implementation method of the vehicle polymerization display based on recorder management platform | |
| CN109032325A (en) | A kind of acquisition methods, system, device and the readable storage medium storing program for executing of controller power consumption |
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 |