CN115630592A - Pressure initialization method and device for flow field - Google Patents
Pressure initialization method and device for flow field Download PDFInfo
- Publication number
- CN115630592A CN115630592A CN202211659815.1A CN202211659815A CN115630592A CN 115630592 A CN115630592 A CN 115630592A CN 202211659815 A CN202211659815 A CN 202211659815A CN 115630592 A CN115630592 A CN 115630592A
- Authority
- CN
- China
- Prior art keywords
- pressure
- flow field
- boundary
- target flow
- target
- 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
-
- 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 provides a pressure initialization method and a pressure initialization device for a flow field, wherein the method comprises the following steps: acquiring a pressure form of a target flow field and boundary pressure types of all boundaries; determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type; if the boundary pressure type of the target flow field including at least one boundary is the designated pressure type, acquiring a reference point corresponding to the designated pressure type and reference pressure of the reference point; if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, acquiring a reference point selected by a user and reference pressure of the reference point; determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure; and determining a corresponding pressure initial formula according to the initial pressure boundary condition and the pressure form, and calculating the initial pressure of the target flow field. The technical problem that only the initial pressure of a target flow field is uniformly initialized in the prior art is solved.
Description
Technical Field
The present disclosure relates to the field of fluid mechanics, and in particular, to a method and an apparatus for initializing pressure of a flow field.
Background
In common CFD (Computational Fluid Dynamics) software, such as Fluent and OpenFOAM, the initial field of pressure in the target flow field is usually set as a uniform field, and then the Navier-Stokes equation is solved. However, in practical cases, the pressure initial field in the target flow field may not be a uniform field, and if the pressure initial field of the target flow field is set to be a uniform field, the convergence of the calculated navier-stokes equation may be slow, or a non-physical calculation result may be caused.
Disclosure of Invention
In view of this, an object of the present application is to provide at least a method and an apparatus for initializing pressure of a flow field, which calculate a pressure initialization formula according to a pressure form of a target process and boundary pressure types of all boundaries, so as to solve the technical problem in the prior art that pressure initialization is performed only by using a target flow field as a uniform field without considering the boundary pressure type of the target flow field, and achieve the technical effect of improving accuracy of pressure initialization and thus improving accuracy of solving a navier-stokes equation.
The application mainly comprises the following aspects:
in a first aspect, an embodiment of the present application provides a method for initializing a pressure of a flow field, where the method includes: acquiring a pressure form of a target flow field and boundary pressure types of all boundaries, wherein the boundary pressure types are used for describing whether the stress of a boundary surface of the target flow field is designated pressure or not; determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type; if the boundary pressure type of the target flow field including at least one boundary is the designated pressure type, acquiring a reference point corresponding to the designated pressure type and reference pressure of the reference point; if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, acquiring a reference point selected by a user and reference pressure of the reference point; determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure; and determining the initial pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form.
Optionally, if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure type, obtaining a reference point selected by a user and a reference pressure of the reference point, including: if the boundary pressure types of all the boundaries do not contain the specified pressure type, determining a target surface with a zero height vector in the target flow field according to the flow field form of the target flow field; and acquiring a reference point selected by a user in the target surface and a reference pressure of the reference point.
Optionally, the flow field morphology comprises: the device comprises a first form and a second form, wherein the first form is used for describing that the boundary surfaces of a target flow field are all wall surfaces, and the second form is used for describing that the target flow field contains a plurality of fluids with different densities; determining a target surface with a zero height vector in a target flow field according to the flow field shape of the target flow field, comprising: if the flow field form is the first form, setting the surface with the height of 0 of the target flow field as a target surface; if the flow field shape is the second shape, the interface of any two adjacent fluids in the target flow field is set as the target surface.
Optionally, the pressure form comprises an original form, a first transformed form, and a second transformed form; the boundary pressure types include: the method comprises the steps of specifying a pressure type and a non-specified pressure type, wherein the non-specified pressure type is used for describing that the stress of a boundary surface of a target flow field is not specified pressure; determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure, wherein the initial pressure boundary condition comprises the following steps: in the case of a pressure form in its original form: if the boundary pressure types of all the boundaries of the target flow field are the specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
if the boundary pressure types of all the boundaries of the target flow field comprise a non-specified pressure type and a specified pressure type, determining the initial pressure boundary condition of the target flow field comprises the following steps:
in the above-mentioned formula,in order to be the density of the fluid,in order to be the acceleration of the gravity,as a reference point (x) in its original form 0 ,y 0 ,z 0 ) Is measured with respect to the reference pressure of (c),is the fluid pressure at any point (x, y, z) within the target flow field in its original form,is the distance from any point (x, y, z) in the target flow field to the reference point,it is referred to the fluid pressure at the boundary surface,which refers to the normal to the boundary surface.
Optionally, determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all the boundaries, the reference point and the reference pressure, including: in the case where the pressure form is the first transform form:
if the boundary pressure types of all the boundaries of the target flow field are the specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
if the boundary pressure types of all boundaries of the target flow field comprise a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field comprises the following steps:
in the above-mentioned formula,in order to be the density of the fluid,is the fluid pressure field magnitude of the target flow field,in order to be the acceleration of the gravity,is a reference point (x) 0 ,y 0 ,z 0 ) Is measured with respect to the reference pressure value of (c),the fluid pressure value at any point in the target flow field,for any point in the target stream field (x, y, z)The distance to the reference point, n referring to the normal,refers to the height of the boundary surface, and h refers to the vertical height from any point (x, y, z) in the target flow field to the reference point.
Optionally, determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all the boundaries, the reference point and the reference pressure, including: in the case where the pressure form is the second transformation form:
if the boundary pressure types of all the boundaries of the target flow field are the specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
if the boundary pressure types of all the boundaries of the target flow field comprise a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field comprises the following steps:
in the above formula, ρ is the fluid density,in order to have a predetermined density, the density of the coating,is the density of the boundary surface(s),the fluid pressure field magnitude, n refers to the normal,in order to be the acceleration of the gravity,is a reference point (x) 0 ,y 0 ,z 0 ) Is measured with respect to the reference pressure value of (c),the fluid pressure value of any point (x, y, z) in the target flow field,is the distance length from any point (x, y, z) in the target stream field to the reference point.
Optionally, determining the initial pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form, including: calculating a first target pressure of a target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form; calculating a first target density of the target flow field according to a first target pressure and a pressure initial formula; solving a second target pressure according to the first target density and the pressure initial formula; calculating a second target density of the target flow field according to a second target pressure and a pressure initial formula; calculating a first difference between the second target pressure and the first target pressure, and a second difference between the second target density and the first target density; determining whether the first difference value and the second difference value respectively belong to corresponding difference value ranges; if the first difference value and the second difference value belong to the corresponding difference value ranges respectively, the second target pressure is used as the initialization pressure of the target flow field; and if the first difference value and the second difference value do not belong to the corresponding difference value range respectively, taking the second target pressure as a new first target pressure, skipping to the step of calculating a second target density of the target flow field according to the first target pressure and the pressure initial formula, and continuing to execute the step.
In a second aspect, an embodiment of the present application further provides a pressure initialization apparatus for a flow field, where the apparatus includes: the first acquisition module is used for acquiring the pressure form of the target flow field and the boundary pressure types of all boundaries, wherein the boundary pressure types are used for describing whether the stress of the boundary surface of the target flow field is designated pressure or not; the first determining module is used for determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type; the second acquisition module is used for acquiring a reference point corresponding to the specified pressure type and reference pressure of the reference point if the boundary pressure type of the target flow field including at least one boundary is the specified pressure type; the third acquisition module is used for acquiring a reference point selected by a user and reference pressure of the reference point if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types; the second determination module is used for determining the initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure; and the third determining module is used for determining the initialization pressure of the target flow field according to the initial pressure boundary condition and the pressure initial formula corresponding to the pressure form.
In a third aspect, an embodiment of the present application further provides an electronic device, including: a processor, a memory and a bus, the memory storing processor-executable machine-readable instructions, the processor and the memory communicating via the bus when the electronic device is running, the machine-readable instructions being executed by the processor to perform the steps of the method for initializing a pressure of a flow field in the first aspect or any one of the possible embodiments of the first aspect.
In a fourth aspect, this application provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to perform the step of initializing the pressure of the flow field in the first aspect or any one of the possible implementations of the first aspect.
The embodiment of the application provides a method and a device for initializing the pressure of a flow field, wherein the method comprises the following steps: acquiring a pressure form of a target flow field and boundary pressure types of all boundaries, wherein the boundary pressure types are used for describing whether the stress of a boundary surface of the target flow field is designated pressure or not; determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type; if the boundary pressure type of the target flow field including at least one boundary is the designated pressure type, acquiring a reference point corresponding to the designated pressure type and reference pressure of the reference point; if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, acquiring a reference point selected by a user and reference pressure of the reference point; determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure; and determining the initial pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form. According to the method and the device, the pressure initial formula is calculated according to the pressure form of the target flow and the boundary pressure types of all boundaries, the technical problem that in the prior art, the boundary pressure type of the target flow field is not considered, and only the target flow field is used as a uniform field for pressure initialization is solved, so that the technical effect of improving the accuracy of pressure initialization and further improving the accuracy of solving the Navier-Stokes equation is achieved.
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.
Fig. 1 is a flowchart illustrating a pressure initialization method for a flow field according to an embodiment of the present disclosure.
Fig. 2 shows a schematic diagram of a target flow field provided by an embodiment of the present application, which includes at least one boundary of all boundaries as a specified pressure type.
Fig. 3 is a schematic diagram illustrating a target flow field in a first configuration according to an embodiment of the present application.
Fig. 4 is a schematic diagram illustrating a target flow field in a second configuration according to an embodiment of the present application.
Fig. 5 is a functional block diagram illustrating a pressure initialization method of a flow field provided in an embodiment of the present application.
Fig. 6 shows a schematic structural diagram of an electronic device provided in an embodiment of the present application.
Detailed Description
In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it should be understood that the drawings in the present application are for illustrative and descriptive purposes only and are not intended to limit the scope of the present application. Additionally, it should be understood that the schematic drawings are not necessarily drawn to scale. The flowcharts used in this application illustrate operations implemented according to some embodiments of the present application. It should be understood that the operations of the flow diagrams may be performed out of order, and that steps without logical context may be reversed in order or performed concurrently. One skilled in the art, under the guidance of this application, may add one or more other operations to, or remove one or more operations from, the flowchart.
In addition, the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, presented in the accompanying drawings, is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present application without making any creative effort, shall fall within the protection scope of the present application.
In the prior art, the pressure initial fields of the target flow field are all set as uniform fields, which are inconsistent with the actual situation, and the convergence of the calculated Navier-Stokes equation after the uniform fields are initialized is slow, which may also result in non-physical calculation results.
Based on this, the embodiment of the application provides a method and a device for initializing pressure of a flow field, and by calculating a pressure initial formula according to a pressure form of a target process and boundary pressure types of all boundaries, the technical problem that pressure initialization is performed only by using a target flow field as a uniform field without considering the boundary pressure type of the target flow field in the prior art is solved, and the technical effect of improving the accuracy of pressure initialization and thus improving the accuracy of solving a navier-stokes equation is achieved. The method comprises the following specific steps:
referring to fig. 1, fig. 1 is a flowchart illustrating a pressure initialization method for a flow field according to an embodiment of the present disclosure. As shown in fig. 1, a pressure initialization method for a flow field provided in an embodiment of the present application includes the following steps:
s101: and acquiring the pressure form of the target flow field and the boundary pressure types of all boundaries.
The boundary pressure type is used for describing whether the stress of the boundary surface of the target flow field is a specified pressure. The boundary pressure types include: and specifying a pressure type and a non-specified pressure type, wherein the specified pressure type is used for describing the stress of the boundary surface of the target flow field as the specified pressure, and the non-specified pressure type is used for describing that the stress of the boundary surface of the target flow field is not the specified pressure. All boundaries of the target flow field may be boundaries containing designated pressure types, boundaries containing non-designated pressure types, or boundaries containing both designated pressure types and non-designated pressure types.
The boundary pressure type of the boundary is the specified pressure type, which means that the boundary is the specified pressure boundary, that is, the stress of the boundary is the specified pressure value. The unspecified pressure type comprises a wall surface boundary and an open boundary, wherein the open boundary is a far field boundary at the boundary of the target flow field, and the boundary is in an undisturbed state, namely the velocity gradient is 0; the wall surface boundary is the boundary of the target flow field, and the boundary is in an undisturbed state, namely the velocity gradient is 0.
The pressure form comprises an original form, a first transformation form and a second transformation form, wherein the first transformation form and the second transformation form are derived from the original form, and a user can randomly select one of the three pressure forms, and the selection of the pressure form does not influence the calculation result.
Acquiring the pressure form of the target flow field and the boundary pressure types of all boundaries comprises the following steps: acquiring a pressure form of a target process selected by a user; acquiring a preprocessing file of Computer Aided Engineering (CAE) simulation corresponding to a target flow field, wherein the preprocessing file comprises information such as grids of the target flow field, the total number of boundaries, the type of boundary pressure and the like; and acquiring boundaries with a target number of specified pressure types and boundaries with a total number of boundaries minus a target number of non-specified pressure types from all boundaries of the target flow field from the pre-processing file.
S102: determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type.
That is, it is determined whether a boundary pressure type including all boundaries among all boundaries of the target flow field is a boundary of the specified pressure type.
S103: and acquiring a reference point corresponding to the specified pressure type and a reference pressure of the reference point.
And if the boundary pressure type of the target flow field including at least one boundary is the specified pressure type, acquiring a reference point corresponding to the specified pressure type and the reference pressure of the reference point.
That is, if the boundary pressure type of the target flow field including at least one boundary is the specified pressure type, the reference point and the reference pressure of the reference point do not need to be specified by the user. At this time, a boundary is randomly selected from at least one boundary of the specified pressure type, any point on the boundary is taken as a reference point, and a specified pressure value corresponding to the boundary is taken as a reference pressure of the reference point.
For example, referring to fig. 2, fig. 2 is a schematic diagram illustrating that at least one of all boundaries of a target flow field provided by an embodiment of the present application is a specified pressure type. As shown in FIG. 2, the target streamThe field is a shaded portion, the upper surface of the shaded portion is a boundary of a specified pressure type, and the specified pressure value of the boundary is p ref Optionally, a point in the boundary is used as a reference point, and the height of the reference point is 0, i.e. h 0 =0, reference pressure p of the reference point 0 =p ref . Since the upper surface of the shaded portion is set as the boundary of the specified pressure type, the vertical height of any point of the target flow field from the reference point is-h.
S104: a reference point selected by a user and a reference pressure of the reference point are acquired.
And if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, acquiring a reference point selected by a user and the reference pressure of the reference point.
If the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, acquiring a reference point selected by a user and reference pressure of the reference point, wherein the reference pressure comprises the following steps: if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, determining a target surface with a zero height vector in the target flow field according to the flow field shape of the target flow field; and acquiring a reference point selected by a user in the target surface and a reference pressure of the reference point.
The flow field morphology includes: the device comprises a first form and a second form, wherein the first form is used for describing that boundary surfaces of a target flow field are wall surfaces, and the second form is used for describing that the target flow field contains a plurality of fluids with different densities.
Determining a target surface with a zero height vector in a target flow field according to the flow field shape of the target flow field, comprising: if the flow field form is the first form, setting the surface of the target flow field with the height of 0 as a target surface; if the flow field shape is the second shape, the interface of any two adjacent fluids in the target flow field is set as the target surface.
For example, referring to fig. 3, fig. 3 is a schematic diagram of a target flow field in a first form according to an embodiment of the present application. As shown in fig. 3, all the boundary surfaces of the target flow field are wall surfaces, and none of the boundary surfaces is a designated pressure type, and at this time, any point of the lowest point of the target flow or the lowest surface with the height of 0 is taken as a reference point, andsetting the height of the reference point to 0, h 0 =0, and sets the reference pressure p of the reference point 0 =p ref That is, the reference pressure is set by the user, and the height value h of any point in the target flow field at this time is the vertical height from the lowest point or from the lowest surface.
For example, referring to fig. 4, fig. 4 is a schematic diagram of a target flow field in a second form according to an embodiment of the present application. As shown in fig. 4, the target flow field contains two fluids, and the two fluids have different densities, so that an interface is formed between the two fluids in the initial state of the target flow field, any point selected by a user in the interface is used as a reference point, and the height of the point is set to be 0, that is, h 0 =0, setting a reference pressure p of the reference point 0 =p ref . Therefore, at this time, the vertical height from any point on the upper part of the interface to the interface or the reference point in the target flow field is h, and the vertical height from any point on the lower part of the interface to the interface or the reference point in the target flow field is-h.
After step S103 or step S104 is executed, step S105 is executed: and determining the initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure.
S105: and determining the initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure.
S106: and determining the initial pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form.
Determining the initial pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form, wherein the pressure initial formula comprises the following steps: calculating a first target pressure of a target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form; calculating a first target density of the target flow field according to a first target pressure and a pressure initial formula; solving a second target pressure according to the first target density and the pressure initial formula; calculating a second target density of the target flow field according to a second target pressure and a pressure initial formula; calculating a first difference between the second target pressure and the first target pressure, and a second difference between the second target density and the first target density; determining whether the first difference value and the second difference value respectively belong to corresponding difference value ranges; if the first difference value and the second difference value belong to the corresponding difference value ranges respectively, the second target pressure is used as the initialization pressure of the target flow field; and if the first difference value and the second difference value do not belong to the corresponding difference value range respectively, taking the second target pressure as a new first target pressure, skipping to the step of calculating a second target density of the target flow field according to the first target pressure and the pressure initial formula, and continuing to execute the step.
Determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure, wherein the initial pressure boundary condition comprises the following steps:
in the case of the pressure form being the original form, the conservation of momentum equation for the Navier-Stokes is as follows:
in the formula (1), the first and second groups,is the fluid density of the target area, t is time,is the fluid velocity of the target area and,is the fluid pressure in the target area and,is the first fluid viscosity coefficient of the target area,is the temperature of the fluid in the target area,is the second fluid viscosity coefficient of the target area and ,is the acceleration of gravity.
Equation (1) can also be written as:
when the boundary pressure type is a non-specified pressure type, the momentum conservation equation can be simplified as follows:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
when the boundary pressure type is the designated pressure type, the velocity gradient of the boundary is still 0, and thus the formula (4) can still be obtained.
That is, in the case where the pressure form is the original form: if the boundary pressure types of all the boundaries of the target flow field are the designated pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
if the boundary pressure types of all the boundaries of the target flow field include a non-specified pressure type and a specified pressure type (that is, all the boundaries of the target flow field have boundaries of both the non-specified pressure type and the specified pressure type), linearly combining the non-specified pressure type and the specified pressure type, and determining the initial pressure boundary condition of the target flow field includes:
in the formulas (5) to (7),in order to be the density of the fluid,in order to be the acceleration of the gravity,as a reference point (x) in its original form 0 ,y 0 ,z 0 ) The reference pressure of (a) is set,is the fluid pressure at any point (x, y, z) within the target flow field in its original form,is the distance from any point (x, y, z) in the target flow field to the reference point,of fingersIs the pressure of the fluid at the boundary surface,which refers to the normal to the boundary surface.
In the case where the pressure form is the original form, the fluid pressure of the target flow field may be defined as:
in equation (8), p refers to the fluid pressure of the target flow field,is the reference pressure of the reference point, h is the vertical distance from any point in the target flow field to the reference point,is the pressure field magnitude of the target flow field.
Equation (8) can be obtainedWriting it in the form of laplace's equation yields a first pressure initialization formula. The first pressure initial equation is:
the first pressure (can be solved in accordance with the initial pressure boundary condition, i.e., equation (5) or equation (6) or equation (7), in combination with the first pressure initial equationDistribution of (d); calculating a first density of the target flow field according to the first pressure and a first pressure initial formula; solving a second pressure according to the first density and the first pressure initial formula; calculating the target flow according to the second pressure and the first pressure initial formulaA second target density of fields; determining whether a difference between the second pressure and the first pressure falls within a preset pressure difference range and determining whether a difference between the second target density and the first density falls within a preset density difference range; if the difference value between the second pressure and the first pressure belongs to the preset pressure difference value range and the difference value between the second target density and the first density belongs to the preset density difference value range, taking the second pressure as the initialization pressure; and if the difference value between the second pressure and the first pressure does not belong to the preset pressure difference value range, or the difference value between the second target density and the first density does not belong to the preset density difference value range, or the difference value between the second pressure and the first pressure does not belong to the preset pressure difference value range, and the difference value between the second target density and the first density does not belong to the preset density difference value range, calculating the first density of the target flow field by taking the second pressure as the new first pressure according to the first pressure and the first pressure initial formula again.
In the case where the pressure form is the first alternate form, the fluid pressure of the target flow field may be defined as:
the conservation of momentum equation for Navier-Stokes can be written as:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, the following formula (11) can be obtained:
at this time, the process of the present invention,the boundary condition of (2) is no longer normal zero gradient, and the force on the boundary is the same as the force on the inside. Therefore, if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field are the specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
If the boundary pressure types of all the boundaries of the target flow field comprise a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field comprises the following steps:
in the formulas (10) to (16),in order to be the density of the fluid,is the fluid pressure value of the target flow field,in order to be the acceleration of the gravity,is a reference point (x) 0 ,y 0 ,z 0 ) Is measured with respect to the reference pressure value of (c),the fluid pressure value at any point in the target flow field,is the distance from any point (x, y, z) in the target flow field to a reference point, n refers to the normal,refers to the height of the boundary surface, and h refers to the vertical height from any point (x, y, z) in the target flow field to the reference point.
In the case of a pressure profile of the first transformation profile, it can be assumed that the target flow field is an incompressible density fieldA uniform field is formed; for a compressible density field for the target flow field,the derivation can be obtainedWriting the expression in the form of laplace's equation yields a second pressure initialization formula.
The second pressure initial formula is:
the third pressure (can be solved by combining the second pressure initialization formula with the initial pressure boundary condition, i.e., formula (14) or formula (15) or formula (16) (14))Distribution of (c); calculating a third density of the target flow field according to a third pressure and a second pressure initial formula; solving a fourth pressure according to the third density and the second pressure initial formula; calculating a fourth density of the target flow field according to a fourth pressure and a second pressure initial formula; determining whether a difference between the fourth pressure and the third pressure falls within a preset pressure difference range and determining whether a difference between the fourth density and the third density falls within a preset density difference range; if the difference value between the fourth pressure and the third pressure belongs to the preset pressure difference value range and the difference value between the fourth density and the third density belongs to the preset density difference value range, taking the fourth pressure as the initialization pressure; and if the difference value between the fourth pressure and the third pressure does not belong to the preset pressure difference value range, or the difference value between the fourth density and the third density does not belong to the preset density difference value range, or the difference value between the fourth pressure and the third pressure does not belong to the preset pressure difference value range, and the difference value between the fourth density and the third density does not belong to the preset density difference value range, taking the fourth pressure as a new third pressure, and calculating the third density of the target flow field again according to the third pressure and the second pressure initial formula.
In the case where the pressure form is the second alternate form, the fluid pressure of the target flow field may be defined as:
the conservation of momentum equation for navier-stokes can be written as:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, the following formula (19) can be obtained:
further, the first and second liquid crystal display panels,the boundary condition of (2) is no longer a normal zero gradient, so that if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary condition of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field are the specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises the following steps:
If the boundary pressure types of all the boundaries of the target flow field comprise a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field comprises the following steps:
in the equations (18) to (23), ρ is the fluid density,in order to set the density to a predetermined value,is the density of the boundary surface(s),the fluid pressure field magnitude, n refers to the normal,is the acceleration of the gravity, and the acceleration is the acceleration of the gravity,is a reference point (x) 0 ,y 0 ,z 0 ) Is measured with respect to the reference pressure value of (c),is the fluid pressure value at any point (x, y, z) in the target flow field,is the distance length from any point (x, y, z) in the target stream field to the reference point.
According to the formula (18), theCan obtainThe alignment derivative may result in:writing the expression in the form of laplace's equation may result in a third pressure initialization formula.
The third pressure initial formula is:
depending on the initial pressure boundary condition, i.e., equation (22) or equation (23) or equation (24), in combination with the third pressure initial equation, a fifth pressure(s) (24) may be solvedDistribution of (d); calculating a fifth density of the target flow field according to a fifth pressure and a third pressure initial formula; solving a sixth pressure according to a fifth density and a third pressure initial formula; calculating a sixth density of the target flow field according to a sixth pressure and a third pressure initial formula; determining whether a difference between the sixth pressure and the fifth pressure falls within a preset pressure difference range and determining whether a difference between the sixth density and the fifth density falls within a preset density difference range; if the difference value between the sixth pressure and the fifth pressure belongs to a preset pressure difference value range and the difference value between the sixth density and the fifth density belongs to a preset density difference value range, taking the sixth pressure as an initialization pressure; and if the difference value between the sixth pressure and the fifth pressure does not belong to a preset pressure difference value range, or the difference value between the sixth density and the fifth density does not belong to a preset density difference value range, or the difference value between the sixth pressure and the fifth pressure does not belong to a preset pressure difference value range, and the difference value between the sixth density and the fifth density does not belong to a preset density difference value range, taking the sixth pressure as a new fifth pressure, and calculating the fifth density of the target flow field again according to the fifth pressure and a third pressure initial formula.
That is, in the first variation, the fluid pressure of the target flow field defines the medium densityVariable, in a second alternative form, the fluid pressure of the target flow field defines the medium densityIs a fixed value.
Based on the same application concept, the embodiment of the present application further provides a pressure initialization apparatus for a flow field corresponding to the pressure initialization method for a flow field provided in the foregoing embodiment, and since the principle of the apparatus in the embodiment of the present application for solving the problem is similar to the pressure initialization method for a flow field in the foregoing embodiment of the present application, the implementation of the apparatus may refer to the implementation of the method, and repeated details are not described.
As shown in fig. 5, fig. 5 is a functional block diagram of a pressure initialization apparatus for a flow field according to an embodiment of the present disclosure. The pressure initialization apparatus 10 of the flow field includes: a first obtaining module 101, a first determining module 102, a second obtaining module 103, a third obtaining module 104, a second determining module 105 and a third determining module 106.
The first obtaining module 101 is configured to obtain a pressure form of the target flow field and boundary pressure types of all boundaries, where the boundary pressure types are used to describe whether stress on a boundary surface of the target flow field is specified pressure; a first determining module 102, configured to determine whether a boundary pressure type of the target flow field including at least one boundary is a specified pressure type; a second obtaining module 103, configured to obtain a reference point corresponding to a specified pressure type and a reference pressure of the reference point if a boundary pressure type of the target flow field including at least one boundary is the specified pressure type; a third obtaining module 104, configured to obtain a reference point selected by a user and a reference pressure of the reference point if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure type; a second determining module 105, configured to determine an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point, and the reference pressure; the third determining module 106 is configured to determine an initialization pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form.
Based on the same application concept, referring to fig. 6, a schematic structural diagram of an electronic device provided in an embodiment of the present application is shown, where the electronic device 20 includes: a processor 201, a memory 202 and a bus 203, wherein the memory 202 stores machine-readable instructions executable by the processor 201, and when the electronic device 20 is operated, the processor 201 communicates with the memory 202 via the bus 203, and the machine-readable instructions are executed by the processor 201 to perform the steps of the pressure initialization method for the flow field according to any one of the embodiments.
In particular, the machine readable instructions, when executed by the processor 201, may perform the following: acquiring a pressure form of a target flow field and boundary pressure types of all boundaries, wherein the boundary pressure types are used for describing whether the stress of a boundary surface of the target flow field is designated pressure or not; determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type; if the boundary pressure type of the target flow field including at least one boundary is the specified pressure type, acquiring a reference point corresponding to the specified pressure type and reference pressure of the reference point; if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, acquiring a reference point selected by a user and reference pressure of the reference point; determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure; and determining the initial pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form.
Based on the same application concept, embodiments of the present application further provide a computer-readable storage medium, where a computer program is stored, and when the computer program is executed by a processor, the steps of the method for initializing a pressure of a flow field provided by the above embodiments are performed.
Specifically, the storage medium can be a general storage medium, such as a mobile disk, a hard disk, or the like, when a computer program on the storage medium is run, the pressure initialization method of the flow field can be executed, and by calculating the pressure initialization formula according to the pressure form of the target flow and the boundary pressure types of all boundaries, the technical problem that in the prior art, the pressure initialization is performed only by using the target flow field as a uniform field without considering the boundary pressure type of the target flow field is solved, so that the technical effect of improving the accuracy of the pressure initialization and thus improving the accuracy of solving the navier-stokes equation is achieved.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again. In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions when actually implemented, and for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or units through some communication interfaces, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one position, or may be distributed on multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solutions of the present application may be embodied in the form of a software product, which is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The above description is only for the specific embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily think of the changes or substitutions within the technical scope of the present application, and shall cover the scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.
Claims (10)
1. A method of pressure initialization of a flow field, the method comprising:
acquiring a pressure form of a target flow field and boundary pressure types of all boundaries, wherein the boundary pressure types are used for describing whether the stress of a boundary surface of the target flow field is designated pressure or not;
determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type;
if the boundary pressure type of at least one boundary in the target flow field is a specified pressure type, acquiring a reference point corresponding to the specified pressure type and reference pressure of the reference point;
if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure types, acquiring a reference point selected by a user and reference pressure of the reference point;
determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all the boundaries, the reference point and the reference pressure;
and determining the initial pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form.
2. The method according to claim 1, wherein if the boundary pressure types of all the boundaries of the target flow field are not the specified pressure type, acquiring a reference point selected by a user and a reference pressure of the reference point, comprising:
if the boundary pressure types of all boundaries of the target flow field are not the specified pressure types, determining a target surface with a zero height vector in the target flow field according to the flow field shape of the target flow field;
and acquiring a reference point selected by a user in the target surface and a reference pressure of the reference point.
3. The method of claim 2, wherein the flow field morphology comprises: the device comprises a first form and a second form, wherein the first form is used for describing that boundary surfaces of the target flow field are wall surfaces, and the second form is used for describing that the target flow field contains a plurality of fluids with different densities;
determining a target surface with a zero height vector in the target flow field according to the flow field shape of the target flow field, including:
if the flow field form is a first form, setting a surface with the height of 0 of the target flow field as a target surface;
and if the flow field shape is the second shape, setting the interface of any two adjacent fluids in the target flow field as a target surface.
4. The method of claim 1, wherein the pressure form comprises a raw form, a first transformed form, and a second transformed form; the boundary pressure types include: the method comprises the steps of specifying a pressure type and a non-specified pressure type, wherein the non-specified pressure type is used for describing that the stress of a boundary surface of a target flow field is not specified pressure;
determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all the boundaries, the reference point and the reference pressure, including:
in the case where the pressure form is the original form:
if the boundary pressure types of all the boundaries of the target flow field are the specified pressure types, determining the initial pressure boundary conditions of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary condition of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field comprise a non-specified pressure type and a specified pressure type, determining the initial pressure boundary condition of the target flow field comprises:
in the above-mentioned formula,in order to be the density of the fluid,in order to be the acceleration of the gravity,as a reference point (x) in its original form 0 ,y 0 ,z 0 ) The reference pressure of (a) is set,is the fluid pressure at any point (x, y, z) within the target flow field in its original form,is the distance from any point (x, y, z) in the target flow field to the reference point,it is referred to the fluid pressure at the boundary surface,which refers to the normal to the boundary surface.
5. The method according to claim 4, wherein the determining an initial pressure boundary condition of the target flow field as a function of the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure comprises:
in the case where the pressure form is the first transform form:
if the boundary pressure types of all the boundaries of the target flow field are the designated pressure types, determining the initial pressure boundary conditions of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary condition of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field comprise a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field comprises:
in the above-mentioned formula,in order to be the density of the fluid,is the fluid pressure field magnitude of the target flow field,in order to be the acceleration of the gravity,is a reference point (x) 0 ,y 0 ,z 0 ) Is measured with respect to the reference pressure value of (c),is the fluid pressure value at any point within the target flow field,is the distance from any point (x, y, z) in the target flow field to a reference point, n refers to the normal,refers to the height of the boundary surface, and h refers to the vertical height from any point (x, y, z) in the target flow field to the reference point.
6. The method according to claim 4, wherein the determining an initial pressure boundary condition of the target flow field as a function of the pressure form, the boundary pressure types of all boundaries, the reference point and the reference pressure comprises:
in the case where the pressure form is the second transformation form:
if the boundary pressure types of all the boundaries of the target flow field are the designated pressure types, determining the initial pressure boundary conditions of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field are non-specified pressure types, determining the initial pressure boundary condition of the target flow field comprises:
if the boundary pressure types of all the boundaries of the target flow field comprise a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field comprises:
in the above formula, ρ is the fluid density,in order to have a predetermined density, the density of the coating,is the density of the boundary surface(s),the fluid pressure field magnitude, n refers to the normal,in order to be the acceleration of the gravity,is a reference point (x) 0 ,y 0 ,z 0 ) Is measured with respect to the reference pressure value of (c),is the fluid pressure value at any point (x, y, z) in the target flow field,is the distance length from any point (x, y, z) in the target flow field to the reference point.
7. The method of claim 1, wherein determining an initialization pressure of the target flow field according to the initial pressure boundary condition and a pressure initialization formula corresponding to the pressure form comprises:
calculating a first target pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form;
calculating a first target density of the target flow field according to the first target pressure and the pressure initial formula;
solving a second target pressure according to the first target density and the pressure initial formula;
calculating a second target density of the target flow field according to the second target pressure and the pressure initial formula;
calculating a first difference between the second target pressure and the first target pressure, and a second difference between the second target density and the first target density;
determining whether the first difference value and the second difference value respectively belong to corresponding difference value ranges;
if the first difference value and the second difference value respectively belong to corresponding difference value ranges, the second target pressure is used as the initialization pressure of the target flow field;
and if the first difference value and the second difference value do not belong to the corresponding difference value range respectively, taking the second target pressure as a new first target pressure, skipping to calculating a second target density of the target flow field according to the first target pressure and the pressure initial formula, and continuing to execute.
8. A pressure initialization apparatus for a flow field, the apparatus comprising:
the first acquisition module is used for acquiring a pressure form of a target flow field and boundary pressure types of all boundaries, wherein the boundary pressure types are used for describing whether the stress of a boundary surface of the target flow field is designated pressure or not;
the first determining module is used for determining whether the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type;
the second acquisition module is used for acquiring a reference point corresponding to the specified pressure type and the reference pressure of the reference point if the boundary pressure type of the target flow field including at least one boundary is the specified pressure type;
the third acquisition module is used for acquiring a reference point selected by a user and reference pressure of the reference point if the boundary pressure types of all boundaries of the target flow field are not the specified pressure types;
a second determining module, configured to determine an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure types of all the boundaries, the reference point, and the reference pressure;
and the third determining module is used for determining the initialization pressure of the target flow field according to the initial pressure boundary condition and the pressure initial formula corresponding to the pressure form.
9. An electronic device, comprising: a processor, a memory and a bus, the memory storing machine-readable instructions executable by the processor, the processor and the memory communicating over the bus when the electronic device is operating, the machine-readable instructions when executed by the processor performing the steps of the pressure initialization method of the flow field according to any one of claims 1 to 7.
10. A computer-readable storage medium, characterized in that a computer program is stored thereon, which computer program, when being executed by a processor, performs the steps of the pressure initialization method of a flow field according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211659815.1A CN115630592B (en) | 2022-12-23 | 2022-12-23 | Pressure initialization method and device for flow field |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211659815.1A CN115630592B (en) | 2022-12-23 | 2022-12-23 | Pressure initialization method and device for flow field |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115630592A true CN115630592A (en) | 2023-01-20 |
CN115630592B CN115630592B (en) | 2023-04-28 |
Family
ID=84910774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211659815.1A Active CN115630592B (en) | 2022-12-23 | 2022-12-23 | Pressure initialization method and device for flow field |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115630592B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120232860A1 (en) * | 2011-03-11 | 2012-09-13 | Desktop Aeronautics, Inc. | Generating a simulated fluid flow over an aircraft surface using anisotropic diffusion |
CN111368487A (en) * | 2020-03-17 | 2020-07-03 | 广西师范大学 | Flow field processing method for simulating periodic movement of particles based on lattice Boltzmann model |
CN111931306A (en) * | 2020-07-31 | 2020-11-13 | 上海交通大学四川研究院 | Method and system for regulating and controlling interaction process of main vortex pair based on influence of auxiliary vortex pair |
CN114282448A (en) * | 2021-11-16 | 2022-04-05 | 北京百度网讯科技有限公司 | Flow field information acquisition method, model training method and device and electronic equipment |
CN115422772A (en) * | 2022-09-27 | 2022-12-02 | 哈尔滨工程大学 | Flow-pressure coupling transient calculation method based on Brayton cycle volume equipment model |
-
2022
- 2022-12-23 CN CN202211659815.1A patent/CN115630592B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120232860A1 (en) * | 2011-03-11 | 2012-09-13 | Desktop Aeronautics, Inc. | Generating a simulated fluid flow over an aircraft surface using anisotropic diffusion |
CN111368487A (en) * | 2020-03-17 | 2020-07-03 | 广西师范大学 | Flow field processing method for simulating periodic movement of particles based on lattice Boltzmann model |
CN111931306A (en) * | 2020-07-31 | 2020-11-13 | 上海交通大学四川研究院 | Method and system for regulating and controlling interaction process of main vortex pair based on influence of auxiliary vortex pair |
CN114282448A (en) * | 2021-11-16 | 2022-04-05 | 北京百度网讯科技有限公司 | Flow field information acquisition method, model training method and device and electronic equipment |
CN115422772A (en) * | 2022-09-27 | 2022-12-02 | 哈尔滨工程大学 | Flow-pressure coupling transient calculation method based on Brayton cycle volume equipment model |
Also Published As
Publication number | Publication date |
---|---|
CN115630592B (en) | 2023-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zahedi et al. | Spurious currents in finite element based level set methods for two‐phase flow | |
JP6825038B2 (en) | Mass exchange model for relative penetration simulation | |
Zheng et al. | Incompressible SPH method based on Rankine source solution for violent water wave simulation | |
Shin et al. | Accurate representation of surface tension using the level contour reconstruction method | |
JP7102741B2 (en) | Fluid analyzer, fluid analysis method, and fluid analysis program | |
US10489937B2 (en) | Paintbrush and liquid simulation | |
JP6671064B2 (en) | Particle simulation apparatus, particle simulation method, and particle simulation program | |
CN109783935B (en) | Implementation method for improving splash fluid stability based on ISPH | |
Mesri et al. | Advanced parallel computing in material forming with CIMLib | |
Mukherjee et al. | Simulating liquid droplets: A quantitative assessment of lattice Boltzmann and Volume of Fluid methods | |
CN113361161B (en) | Shale oil-water two-phase flow simulation method and device based on digital core and storage medium | |
Boroomand et al. | On the use of Laplace's equation for pressure and a mesh-free method for 3D simulation of nonlinear sloshing in tanks | |
Fierz et al. | Maintaining large time steps in explicit finite element simulations using shape matching | |
Gatica et al. | An L p spaces-based formulation yielding a new fully mixed finite element method for the coupled Darcy and heat equations | |
Sheu et al. | Development of a dispersively accurate conservative level set scheme for capturing interface in two-phase flows | |
CN110188407B (en) | Method and device for determining liquid flow parameters in porous medium | |
Darbani et al. | Meshless method for shallow water equations with free surface flow | |
Thömmes et al. | A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method | |
Mao et al. | A conservative and consistent Lagrangian gradient smoothing method for simulating free surface flows in hydrodynamics | |
Chiocchetti et al. | An exactly curl-free staggered semi-implicit finite volume scheme for a first order hyperbolic model of viscous two-phase flows with surface tension | |
Zografos et al. | A viscoelastic two-phase solver using a phase-field approach | |
CN114065593A (en) | Method and device for underground water coupling simulation in multi-scale fracture matrix system | |
CN115630592A (en) | Pressure initialization method and device for flow field | |
JP5221603B2 (en) | Ground deformation analysis device, ground deformation analysis method, program | |
KR102436658B1 (en) | Fluid analysis simulation method and fluid simulation apparatus |
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 |