CN115630592B - Pressure initialization method and device for flow field - Google Patents
Pressure initialization method and device for flow field Download PDFInfo
- Publication number
- CN115630592B CN115630592B CN202211659815.1A CN202211659815A CN115630592B CN 115630592 B CN115630592 B CN 115630592B CN 202211659815 A CN202211659815 A CN 202211659815A CN 115630592 B CN115630592 B CN 115630592B
- 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.)
- Active
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 device 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; determining whether the target flow field contains at least one boundary of a boundary pressure type that is a specified pressure type; if the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type, acquiring a reference point corresponding to the specified pressure type 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 designated pressure types, acquiring a reference point selected by a user and the reference pressure of the reference point; determining an initial pressure boundary condition of a 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 initialization pressure of the target flow field. The technical problem that only the initial pressure of the 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 apparatus for initializing pressure of a flow field.
Background
In common CFD (Computational Fluid Dynamics ) software, such as Fluent and OpenFOAM, the pressure initiation field of the target flow field is typically set to a uniform field, thereby solving the Navier-Stokes equation (Navier-Stokes). In practice, however, the initial field of pressure in the target flow field may not be a uniform field, which may result in slow convergence of the calculated wiener-stokes equation or may result in non-physical calculation results if the initial field of pressure in the target flow field is set to be a uniform field.
Disclosure of Invention
In view of this, the present application aims at providing at least a method and a device for initializing the pressure of a flow field, which solve the technical problem of initializing the pressure of a target flow field as a uniform field without considering the boundary pressure type of the target flow field by calculating a pressure initialization formula according to the pressure form of the target flow and the boundary pressure types of all boundaries, and achieve the technical effect of improving the accuracy of initializing the pressure and thus the accuracy of solving the nano-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, the method including: 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 a specified pressure or not; determining whether the target flow field contains at least one boundary of a boundary pressure type that is a specified pressure type; if the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type, acquiring a reference point corresponding to the specified pressure type 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 designated pressure types, acquiring a reference point selected by a user and the reference pressure of the reference point; determining an initial pressure boundary condition of a 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 initialization 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 types, acquiring the reference point selected by the user and the reference pressure of the reference point includes: if the boundary pressure types of all the boundaries do not contain the designated pressure type, 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 the reference point selected by the user in the target surface and the reference pressure of the reference point.
Optionally, the flow field morphology comprises: the first form and the second form, wherein the first form is used for describing that the boundary surface of the target flow field is a wall surface, and the second form is used for describing that the target flow field contains a plurality of fluids with different densities; according to the flow field morphology of the target flow field, determining a target surface with a zero height vector in the target flow field comprises the following steps: if the flow field shape is the first shape, setting a 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 between any two adjacent fluids in the target flow field is set as a target surface.
Optionally, the pressure form comprises a raw form, a first transformed form, and a second transformed form; boundary pressure types include: a specified pressure type and a non-specified pressure type, the non-specified pressure type being used to describe that the force of the boundary surface of the target flow field is not a 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 the pressure form being 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:
If the boundary pressure types of all the boundaries of the target flow field are non-designated 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 include the unspecified pressure type and the specified pressure type, determining the initial pressure boundary condition of the target flow field includes:
in the above-mentioned formula(s),for fluid density->Acceleration of gravity, ++>Is the reference point (x 0 ,y 0 ,z 0 ) Reference pressure of>Is the fluid pressure at any point (x, y, z) in the original form in the target flow field, +.>For the distance from any point (x, y, z) in the target flow field to the reference point, +.>Refers to the fluid pressure of the boundary surface, +.>Refers to the normal of the boundary surface.
Optionally, determining an initial pressure boundary condition of the target flow field according to the pressure form, the boundary pressure type of all boundaries, the reference point and the reference pressure, including: in the case of a first variant of the pressure 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-designated 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 include the specified pressure type and the unspecified pressure type, determining the initial pressure boundary condition of the target flow field includes:
in the above-mentioned formula(s),for fluid density->Fluid pressure field for the target flow field, +.>Acceleration of gravity, ++>For the reference point (x 0 ,y 0 ,z 0 ) Reference pressure value of>For the fluid pressure value at any point in the target flow field, is->For the distance from any point (x, y, z) in the target flow field to the reference point, n refers to normal,/o>Refers to the height of the boundary surface, h refers to the vertical height from any point (x, y, z) within 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 type of all boundaries, the reference point and the reference pressure, including: in the case of a second variant of the pressure 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-designated 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 include the specified pressure type and the unspecified pressure type, determining the initial pressure boundary condition of the target flow field includes:
In the above formula, ρ is the fluid density,is of preset density (I)>Is the density of boundary surface->The fluid pressure field quantity for the target flow field, n being the normal, < >>Acceleration of gravity, ++>For the reference point (x 0 ,y 0 ,z 0 ) Reference pressure value of>For 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.
Optionally, determining the initialization pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form includes: calculating a first target pressure of the target flow field according to an initial pressure boundary condition and a pressure initial formula corresponding to a 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 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 respectively belong to the corresponding difference value ranges, taking the second target pressure as the initializing pressure of the target flow field; if the first difference value and the second difference value do not belong to the corresponding difference value ranges, taking the second target pressure as the new first target pressure, jumping to calculate the second target density of the target flow field according to the first target pressure and the pressure initial formula, and continuing to execute.
In a second aspect, embodiments of the present application further provide a pressure initializing device for a flow field, where the device 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 the designated pressure or not; a first determining module for determining whether the target flow field contains at least one boundary of a boundary pressure type that 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 containing at least one boundary is the specified pressure type; the third acquisition module is used for acquiring the reference point selected by the 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 designated pressure types; the second determining module is used for 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 the third determining module is used for determining the initialization pressure of the target flow field according to the initial pressure boundary condition and a pressure initial formula corresponding to the pressure form.
In a third aspect, embodiments of the present application further provide an electronic device, including: a processor, a memory and a bus, the memory storing machine readable instructions executable by the processor, the processor and the memory in communication via the bus when the electronic device is in operation, the machine readable instructions when executed by the processor perform the steps of the method for pressure initialization of a flow field in any one of the possible embodiments of the first aspect or the first aspect described above.
In a fourth aspect, the embodiments of the present application further provide a computer readable storage medium having a computer program stored thereon, which when executed by a processor performs the step of pressure initialization of the flow field in the first aspect or any of the possible embodiments 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 a specified pressure or not; determining whether the target flow field contains at least one boundary of a boundary pressure type that is a specified pressure type; if the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type, acquiring a reference point corresponding to the specified pressure type 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 designated pressure types, acquiring a reference point selected by a user and the reference pressure of the reference point; determining an initial pressure boundary condition of a 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 initialization 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, so that the technical problem that in the prior art, the target flow field is used as a uniform field only for pressure initialization without considering the boundary pressure types of the target flow field is solved, and the technical effect of improving the accuracy of pressure initialization and thus the accuracy of solving the Navier-Stokes equation is achieved.
In order to make the above 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 needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered limiting the scope, and that other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 shows a flow chart of a method for initializing the pressure of a flow field according to an embodiment of the present application.
Fig. 2 shows a schematic diagram of all boundaries of a target flow field provided by an embodiment of the present application, including at least one boundary of a specified pressure type.
Fig. 3 shows a schematic view of a target flow field provided in an embodiment of the present application in a first configuration.
Fig. 4 shows a schematic view of a target flow field provided in an embodiment of the present application in a second configuration.
Fig. 5 shows a functional block diagram of a method for initializing the pressure of a flow field according to an embodiment of the present application.
Fig. 6 shows a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more clear, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application, and it should be understood that the accompanying drawings in the present application are only for the purpose of illustration and description, and are not intended to limit the protection scope of the present application. In addition, it should be understood that the schematic drawings are not drawn to scale. A flowchart, as used in this application, illustrates operations implemented according to some embodiments of the present application. It should be appreciated that the operations of the flow diagrams may be implemented out of order and that steps without logical context may be performed in reverse order or concurrently. Moreover, one or more other operations may be added to the flow diagrams and one or more operations may be removed from the flow diagrams as directed by those skilled in the art.
In addition, the described embodiments are only some, but not all, of the embodiments of the present application. The components of the embodiments of the present application, which are generally described and illustrated in the figures herein, may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, as provided in the accompanying drawings, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present application without making any inventive effort, are intended to be within the scope of the present application.
In the prior art, the initial pressure fields of the target flow fields are all set to be uniform fields, which are inconsistent with the actual situation, and the convergence of the Navier-Stokes equation is slow after the uniform fields are initialized, so that non-physical calculation results can be caused.
Based on the above, the embodiment of the application provides a method and a device for initializing the pressure of a flow field, which solve the technical problem that in the prior art, the pressure initialization is carried out by taking the target flow field as a uniform field only without considering the boundary pressure type of the target flow field by calculating a pressure initialization formula according to the pressure form of the target flow and the boundary pressure types of all boundaries, thereby achieving the technical effect of improving the accuracy of the pressure initialization and further improving the accuracy of solving a Navier-Stokes equation. The method comprises the following steps:
referring to fig. 1, fig. 1 is a flowchart of a method for initializing a flow field according to an embodiment of the present application. As shown in fig. 1, the method for initializing the pressure of the flow field provided in the embodiment of the present application includes the following steps:
s101: the pressure form of the target flow field and the boundary pressure types of all boundaries are acquired.
The boundary pressure type is used to describe whether the force applied to the boundary surface of the target flow field is a specified pressure. Boundary pressure types include: 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 the stress of the boundary surface of the target flow field as the non-specified pressure. All of the boundaries of the target flow field may be boundaries of the specified pressure type, boundaries of the non-specified pressure type, boundaries of both the specified pressure type and boundaries of the non-specified pressure type.
The boundary pressure type of the boundary being a specified pressure type means that the boundary is a specified pressure boundary, i.e. the force of the boundary is a specified pressure value. The non-designated pressure type comprises a wall 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 boundary of the wall surface 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 neither of the three pressure forms is influenced by the calculation result.
Acquiring the pressure form of the target flow field and the boundary pressure types of all boundaries comprises: acquiring a pressure form of a target process selected by a user; acquiring a pretreatment file of computer aided engineering (computer aided engineering, CAE) simulation corresponding to a target flow field, wherein the pretreatment file contains information such as grids, total number of boundaries, boundary pressure type and the like of the target flow field; the boundaries of the specified pressure type with the target number and the boundaries of the unspecified pressure type with the target number subtracted from the total number of the boundaries are obtained from the pretreatment file.
S102: it is determined whether the boundary pressure type of the target flow field contains at least one boundary is a specified pressure type.
That is, it is determined whether all boundaries of the target flow field include boundaries of which the boundary pressure type is the specified pressure type.
S103: and acquiring the reference point corresponding to the specified pressure type and the reference pressure of the reference point.
And if the boundary pressure type of the target flow field containing 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 a specified pressure type, the reference point and the reference pressure of the reference point need not be specified by the user. At this time, one 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.
Referring to fig. 2, fig. 2 is a schematic diagram of a target flow field provided in an embodiment of the present application, where at least one boundary is a specified pressure type. As shown in fig. 2, the target flow 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 One point is selected as a reference point in the boundary, and the height of the reference point is 0, namely h 0 Reference pressure p of the reference point =0 0 =p ref . Since the upper surface of the hatched portion is set as a boundary designating the pressure type, the vertical height of any point of the target flow field from the reference point is-h.
S104: the reference point selected by the user and the reference pressure of the reference point are obtained.
And if the boundary pressure types of all the boundaries of the target flow field are not the designated pressure types, acquiring the reference point selected by the 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 designated pressure types, acquiring the reference point selected by the user and the reference pressure of the reference point, wherein the method comprises the following steps: if the boundary pressure types of all the boundaries of the target flow field are not the designated 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 the reference point selected by the user in the target surface and the reference pressure of the reference point.
The flow field comprises the following forms: the first form and the second form, wherein the first form is used for describing that the boundary surface of the target flow field is a wall surface, and the second form is used for describing that the target flow field contains multiple fluids with different densities.
According to the flow field morphology of the target flow field, determining a target surface with a zero height vector in the target flow field comprises the following steps: if the flow field shape is the first shape, setting a 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 between any two adjacent fluids in the target flow field is set as a target surface.
Referring to fig. 3, fig. 3 is a schematic view of a target flow field in a first configuration according to an embodiment of the present application. As shown in FIG. 3, the boundary surfaces of the target flow field are all wall surfaces, no boundary is of a specified pressure type, at this time, any one of the lowest point of the target flow or the lowest surface with the height of 0 is taken as a reference point, and the height of the reference point is set to 0, namely h 0 =0, and sets the reference pressure p of the reference point 0 =p ref I.e., the reference pressure is set by the user, where the height value h of any point in the target flow field is the vertical height from the lowest point or from the lowest plane.
Referring to fig. 4, fig. 4 is a schematic view of a target flow field in a second configuration according to an embodiment of the present application. As shown in fig. 4, the target flow field contains two fluids, and the densities of the two fluids are different, 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 taken as a reference point, and the height of the point is set to be 0, namely h 0 Reference pressure p of reference point is set =0 0 =p ref . Thus, at this time, the target flow field is at the upper part of the interfaceThe vertical height from any point of the interface or the reference point is h, and the vertical height from any point of the object flow field, which is positioned at the lower part of the interface, to the interface or the reference point is-h.
After step S103 or step S104 is performed, step S105 is performed: 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 initialization 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 an initial pressure boundary condition and a pressure initial formula corresponding to a pressure form, determining an initial pressure of a target flow field comprises the following steps: calculating a first target pressure of the target flow field according to an initial pressure boundary condition and a pressure initial formula corresponding to a 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 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 respectively belong to the corresponding difference value ranges, taking the second target pressure as the initializing pressure of the target flow field; if the first difference value and the second difference value do not belong to the corresponding difference value ranges, taking the second target pressure as the new first target pressure, jumping to calculate the second target density of the target flow field according to the first target pressure and the pressure initial formula, and continuing to execute.
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 momentum conservation equation for the Navie-Stokes is as follows:
in the formula (1),for the fluid density of the target area, t is time, < >>For the fluid velocity of the target area, +.>For the fluid pressure of the target area, +.>For the first fluid viscosity coefficient of the target area, +.>For the fluid temperature of the target area, +.>Is the second fluid viscosity coefficient of the target area and +.> ,/>Gravitational acceleration.
Equation (1) can also be written as:
when the boundary pressure type is a non-specified pressure type, the momentum conservation equation can be reduced to:
if the boundary pressure types of all the boundaries of the target flow field are non-designated pressure types, determining the initial pressure boundary condition of the target flow field comprises:
when the boundary pressure type is the specified pressure type, the velocity gradient of the boundary is still 0, so that equation (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:
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, there is a boundary of both the non-specified pressure type and the specified pressure type in all the boundaries of the target flow field), the determining the initial pressure boundary condition of the target flow field includes:
in the formulas (5) to (7),for fluid density->Acceleration of gravity, ++>Is the reference point (x 0 ,y 0 ,z 0 ) Reference pressure of>Is the fluid pressure at any point (x, y, z) in the original form in the target flow field, +.>For the distance from any point (x, y, z) in the target flow field to the reference point, +.>Refers to the fluid pressure of the boundary surface, +.>Refers to the normal of 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,as 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 of the target flow field.
Equation (8) can be derivedWriting it in the form of a Laplace equation yields a first pressure initial equation. The first pressure initial formula is:
According to the initial pressure boundary condition, namely formula (5) or formula (6) or formula (7), the first pressure can be solved by combining the first pressure initial formulaIs a distribution of (a); 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 a second target density of the target flow field according to the second pressure and the first pressure initial formula; determining whether a difference between the second pressure and the first pressure belongs to a preset pressure difference range and determining whether a difference between the second target density and the first density belongs to a preset density difference range; if the difference value between the second pressure and the first pressure belongs to a preset pressure difference value range and the difference value between the second target density and the first density belongs to a preset density difference value range, the second pressure is used as an initialization pressure; if the difference between the second pressure and the first pressure does not belong to the preset pressure difference range, or the difference between the second target density and the first density does not belong to the preset density difference range, or the difference between the second pressure and the first pressure does not belong to the preset pressure difference range and the second target density and the first density And if the difference value does not belong to the preset density difference value range, taking the second pressure as a new first pressure, and calculating the first density of the target flow field again according to the first pressure and the first pressure initial formula.
In the case where the pressure form is the first variation, the fluid pressure of the target flow field may be defined as:
the momentum conservation equation for the Navie-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 liquid crystal display device,the boundary condition of (2) is no longer normal zero gradient and the force of the boundary is the same as the force of the interior. Thus, 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 includes:
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 include the specified pressure type and the unspecified pressure type, determining the initial pressure boundary condition of the target flow field includes:
In the formulas (10) to (16),for fluid density->For the fluid pressure value of the target flow field, +.>Acceleration of gravity, ++>For the reference point (x 0 ,y 0 ,z 0 ) Reference pressure value of>For the fluid pressure value at any point in the target flow field, is->For the distance from any point (x, y, z) in the target flow field to the reference point, n refers to normal,/o>Refers to the height of the boundary surface, h refers to the vertical height from any point (x, y, z) within the target flow field to the reference point.
In the case of the first transformation of the pressure form, the target flow field is an incompressible density field, which can be considered asIs a uniform field; for a compressible density field of the target flow field, < > is>The derivation can be achieved>Writing the expression in the form of a Laplace equation yields a second pressure initiation formula.
The second pressure initial formula is:
according to the initial pressure boundary condition, namely the formula (14) or the formula (15) or the formula (16), the third pressure can be solved by combining the second pressure initial formulaIs a distribution of (a); 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 belongs to a preset pressure difference range and determining whether a difference between the fourth density and the third density belongs to 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, the fourth pressure is used as an initialization pressure; if the difference value between the fourth pressure and the third pressure does not belong to the preset pressure difference range, or the difference value between the fourth density and the third density does not belong to the preset density difference range, or the difference value between the fourth pressure and the third pressure does not belong to the preset pressure difference range and the difference value between the fourth density and the third density does not belong to the preset density difference range, the fourth pressure is taken as a new third pressure, and the third density of the target flow field is calculated again according to the third pressure and the second pressure initial formula.
In the case where the pressure form is the second variation, the fluid pressure of the target flow field may be defined as:
the momentum conservation equation for the Navie-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 processing unit is used for processing the data,the boundary condition of (2) is no longer a normal zero gradient, such that if the boundary pressure type of all boundaries of the target flow field are of a non-specified pressure type, 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 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 include the specified pressure type and the unspecified pressure type, determining the initial pressure boundary condition of the target flow field includes:
in the formulas (18) to (23), ρ is the fluid density,is of preset density (I)>Is the density of boundary surface->The fluid pressure field quantity for the target flow field, n being the normal, < >>Acceleration of gravity, ++>For the reference point (x 0 ,y 0 ,z 0 ) Reference pressure value of >For 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. />
Can be obtained according to the formula (18)Can obtainThe alignment derivative may be obtained: />The expression is written in the form of a Laplace equation to yield a third pressure initiation formula.
The third pressure initial formula is:
according to the boundary condition of the initial pressure, namely the formula (22) or the formula (23) or the formula (24), the fifth pressure can be solved by combining the third pressure initial formulaIs a distribution of (a); 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 the fifth density and the 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 belongs to a preset pressure difference range and determining whether a difference between the sixth density and the fifth density belongs to a preset density difference range; if the difference value between the sixth pressure and the fifth pressure belongs to the preset pressure difference value range and the difference value between the sixth density and the fifth density belongs to the preset density difference value range, the sixth pressure is used as the initializing pressure; if the difference between the sixth pressure and the fifth pressure does not belong to the preset pressure difference range, or the difference between the sixth density and the fifth density does not belong to the preset density difference range, or the difference between the sixth pressure and the fifth pressure does not belong to the preset pressure difference range and the difference between the sixth density and the fifth density does not belong to the preset density difference range, the sixth pressure is taken as a new fifth pressure, and the fifth density of the target flow field is calculated again according to the fifth pressure and the third pressure initial formula.
That is, in the first variation, the fluid pressure of the target flow field defines the medium densityIn the second variant, the fluid pressure of the target flow field defines the medium density +.>Is a fixed value.
Based on the same application conception, the embodiment of the present application further provides a flow field pressure initializing device corresponding to the flow field pressure initializing method provided in the foregoing embodiment, and since the principle of solving the problem by the device in the embodiment of the present application is similar to that of the flow field pressure initializing method in the foregoing embodiment of the present application, the implementation of the device may refer to the implementation of the method, and the repetition is omitted.
As shown in fig. 5, fig. 5 is a functional block diagram of a pressure initializing device for a flow field according to an embodiment of the present application. The pressure initializing device 10 of the flow field includes: the first acquisition module 101, the first determination module 102, the second acquisition module 103, the third acquisition module 104, the second determination module 105, and the third determination module 106.
A first obtaining module 101, 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 a stress of a boundary surface of the target flow field is a specified pressure; a first determining module 102, configured to determine whether the target flow field includes at least one boundary having a boundary pressure type that is a specified pressure type; a second obtaining module 103, configured to obtain, if the boundary pressure type of the target flow field including at least one boundary is a specified pressure type, a reference point corresponding to the specified pressure type and a reference pressure of the reference point; a third obtaining module 104, configured to obtain the reference point selected by the 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; 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 an initial pressure boundary condition and a pressure initial formula corresponding to the pressure form.
Based on the same application concept, referring to fig. 6, which is a schematic structural diagram of an electronic device provided in an embodiment of the present application, the electronic device 20 includes: processor 201, memory 202, and bus 203, memory 202 storing machine-readable instructions executable by processor 201, which when executed by processor 201 perform the steps of a method for initializing pressure of a flow field as in any of the embodiments described above, when electronic device 20 is in operation, processor 201 and memory 202 communicate via bus 203.
Specifically, machine readable instructions, when executed by 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 a specified pressure or not; determining whether the target flow field contains at least one boundary of a boundary pressure type that is a specified pressure type; if the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type, acquiring a reference point corresponding to the specified pressure type 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 designated pressure types, acquiring a reference point selected by a user and the reference pressure of the reference point; determining an initial pressure boundary condition of a 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 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, the embodiment of the present application further provides 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 pressure initialization method of the flow field provided by the above embodiment are executed.
Specifically, the storage medium can be a general storage medium, such as a mobile disk, a hard disk, and the like, and when a computer program on the storage medium is run, the pressure initialization method of the flow field can be executed, and by calculating a 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 by taking the target flow field as a uniform field without considering the boundary pressure type of the target flow field is solved, and the technical effect of improving the accuracy of the pressure initialization is achieved, thereby improving the accuracy of solving the Navier-Stokes equation.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described system and apparatus may refer to corresponding procedures in the foregoing method embodiments, which are not described herein again. In the several embodiments provided in this application, it should be understood that the disclosed systems, devices, and methods may be implemented in other manners. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be other manners of division in actual implementation, and for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some communication interface, device or unit indirect coupling or communication connection, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in each embodiment of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in 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 essence or a part contributing to the prior art or a part of the technical solutions, or in the form of a software product, which is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
The foregoing is merely a specific embodiment of the present application, but the protection scope of the present application is not limited thereto, and any person skilled in the art can easily think about changes or substitutions within the technical scope of the present application, and the changes or substitutions are covered in the protection 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 a specified pressure or not; the boundary pressure type includes: a specified pressure type for describing the stress of the boundary surface of the target flow field as a specified pressure and a non-specified pressure type for describing the stress of the boundary surface of the target flow field as a non-specified pressure;
determining whether the target flow field contains at least one boundary of a boundary pressure type that is a specified pressure type;
if the boundary pressure type of the target flow field containing at least one boundary is a specified pressure type, acquiring a reference point corresponding to the specified pressure type and the reference pressure of the reference point;
If the boundary pressure types of all the boundaries of the target flow field are non-designated pressure types, acquiring a reference point selected by a user and the 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 initialization 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 of claim 1, wherein if the boundary pressure types of all boundaries of the target flow field are non-specified pressure types, obtaining the user-selected reference point and the reference pressure of the reference point comprises:
if the boundary pressure types of all the boundaries of the target flow field are non-designated 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 the reference pressure of the reference point.
3. The method of claim 2, wherein the flow field morphology comprises: a first state and a second state, wherein the first state is used for describing that boundary surfaces of the target flow field are wall surfaces, and the second state is used for describing that the target flow field contains a plurality of fluids with different densities;
The determining the target surface with the height vector of zero in the target flow field according to the flow field shape of the target flow field comprises the following steps:
if the flow field shape is the first shape, 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 between 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;
said determining an initial pressure boundary condition of said target flow field from said pressure pattern, said boundary pressure type of all boundaries, said reference point and said reference pressure, comprising:
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 condition of the target flow field includes:
if the boundary pressure types of all the boundaries of the target flow field are non-designated pressure types, determining the initial pressure boundary condition of the target flow field includes:
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, determining the initial pressure boundary condition of the target flow field includes:
In the above-mentioned formula(s),for fluid density->Acceleration of gravity, ++>Is the reference point (x 0 ,y 0 ,z 0 ) Reference pressure of>Is the fluid pressure at any point (x, y, z) in the original form in the target flow field, +.>For any point (x, y)Z) distance to the reference point, +.>Refers to the fluid pressure of the boundary surface, +.>Refers to the normal of the boundary surface.
5. The method of claim 4, wherein said determining initial pressure boundary conditions for the target flow field as a function of the pressure form, boundary pressure type for all boundaries, the reference point, and the reference pressure comprises:
in the case where the pressure form is the first conversion 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 condition of the target flow field includes:
if the boundary pressure types of all the boundaries of the target flow field are non-designated pressure types, determining the initial pressure boundary condition of the target flow field includes:
if the boundary pressure types of all the boundaries of the target flow field include a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field includes:
In the above-mentioned formula(s),for fluid density->Fluid pressure field for the target flow field, +.>The acceleration of the gravity is that,for the reference point (x 0 ,y 0 ,z 0 ) Reference pressure value of>For the fluid pressure value at any point in the target flow field, is->For the distance from any point (x, y, z) in the target flow field to the reference point, n refers to normal,/o>Refers to the height of the boundary surface, h refers to the vertical height from any point (x, y, z) within the target flow field to the reference point.
6. The method of claim 4, wherein said determining initial pressure boundary conditions for the target flow field as a function of the pressure form, boundary pressure type for all boundaries, the reference point, and the reference pressure comprises:
in the case where the pressure form is the second variation:
if the boundary pressure types of all the boundaries of the target flow field are the designated pressure types, determining the initial pressure boundary condition of the target flow field includes:
if the boundary pressure types of all the boundaries of the target flow field are non-designated pressure types, determining the initial pressure boundary condition of the target flow field includes:
if the boundary pressure types of all the boundaries of the target flow field include a specified pressure type and a non-specified pressure type, determining the initial pressure boundary condition of the target flow field includes:
In the above formula, ρ is the fluid density,is of preset density (I)>Is the density of boundary surface->The fluid pressure field quantity for the target flow field, n being the normal, < >>Acceleration of gravity, ++>For the reference point (x 0 ,y 0 ,z 0 ) Is used for the reference pressure value of (a),for 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 the initialization pressure of the target flow field according to the initial pressure boundary condition and the 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, 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 the corresponding difference value ranges, the second target pressure is used as the initializing 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 ranges, taking the second target pressure as a new first target pressure, jumping to calculate 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 initialisation apparatus for a flow field, the apparatus comprising:
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 the designated pressure or not; the boundary pressure type includes: a specified pressure type for describing the stress of the boundary surface of the target flow field as a specified pressure and a non-specified pressure type for describing the stress of the boundary surface of the target flow field as a non-specified pressure;
A first determining module for determining whether the target flow field contains at least one boundary of a boundary pressure type that is a specified pressure type;
the second acquisition module is used for acquiring a reference point corresponding to the appointed pressure type and the reference pressure of the reference point if the boundary pressure type of the target flow field containing at least one boundary is the appointed pressure type;
the third acquisition module is used for 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 non-designated 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 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 a pressure initial formula corresponding to the pressure form.
9. An electronic device, comprising: a processor, a memory and a bus, said memory storing machine readable instructions executable by said processor, said processor and said memory in communication via said bus when the electronic device is running, said machine readable instructions when executed by said processor performing the steps of the pressure initialization method of the flowfield of any one of claims 1 to 7.
10. A computer readable storage medium, characterized in that it has stored thereon a computer program which, when executed by a processor, performs the steps of the pressure initialization method of a flow field according to any 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 CN115630592A (en) | 2023-01-20 |
CN115630592B true 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 (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114282448A (en) * | 2021-11-16 | 2022-04-05 | 北京百度网讯科技有限公司 | Flow field information acquisition method, model training method and device and electronic equipment |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8437990B2 (en) * | 2011-03-11 | 2013-05-07 | Aerion Corporation | Generating a simulated fluid flow over an aircraft surface using anisotropic diffusion |
CN111368487B (en) * | 2020-03-17 | 2023-07-18 | 广西师范大学 | 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 |
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 (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114282448A (en) * | 2021-11-16 | 2022-04-05 | 北京百度网讯科技有限公司 | Flow field information acquisition method, model training method and device and electronic equipment |
Also Published As
Publication number | Publication date |
---|---|
CN115630592A (en) | 2023-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zahedi et al. | Spurious currents in finite element based level set methods for two‐phase flow | |
Dhondt | Calculix crunchix user’s manual version 2.12 | |
JP6434512B2 (en) | Temperature coupling algorithm for hybrid hot-lattice Boltzmann method | |
Vacondio et al. | Smoothed Particle Hydrodynamics: Approximate zero‐consistent 2‐D boundary conditions and still shallow‐water tests | |
Elgeti et al. | Deforming fluid domains within the finite element method: five mesh-based tracking methods in comparison | |
JP2019125102A (en) | Fluid analysis device, fluid analysis method, and fluid analysis program | |
Choi | A least‐square weighted residual method for level set formulation | |
Mesri et al. | Advanced parallel computing in material forming with CIMLib | |
Viré et al. | An immersed-shell method for modelling fluid–structure interactions | |
Chou et al. | Extrapolated local radial basis function collocation method for shallow water problems | |
CN109783935B (en) | Implementation method for improving splash fluid stability based on ISPH | |
AU2017227323B2 (en) | Particle simulation device, particle simulation method, and particle simulation program | |
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 | |
Engberg et al. | Numerical simulation of rising droplets in liquid–liquid systems: A comparison of continuous and sharp interfacial force models | |
Capodaglio et al. | A particle tracking algorithm for parallel finite element applications | |
CN115630592B (en) | Pressure initialization method and device for flow field | |
EP3179390A1 (en) | Method and apparatus for modeling movement of air bubble based on fluid particles | |
Hilger et al. | A novel approach to fluid-structure interaction simulations involving large translation and contact | |
Lee | Mathematical model of contractile ring-driven cytokinesis in a three-dimensional domain | |
Jassim et al. | Numerical analysis of Newtonian flows based on artificial compressibility AC method | |
US9690545B2 (en) | Floating-point calculation apparatus, program, and calculation apparatus | |
KR20220012662A (en) | Fluid analysis simulation method and fluid simulation apparatus | |
Joubert et al. | Boundary condition enforcement for renormalised weakly compressible meshless Lagrangian methods | |
CN115906709A (en) | Fluid pressure field prediction method and device | |
KR102181981B1 (en) | Apparatus, method and computer program for sph based fluid analysis simulation |
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 |