CN111783277B - Fluid-solid interface decoupling algorithm, device and equipment - Google Patents
Fluid-solid interface decoupling algorithm, device and equipment Download PDFInfo
- Publication number
- CN111783277B CN111783277B CN202010498716.4A CN202010498716A CN111783277B CN 111783277 B CN111783277 B CN 111783277B CN 202010498716 A CN202010498716 A CN 202010498716A CN 111783277 B CN111783277 B CN 111783277B
- Authority
- CN
- China
- Prior art keywords
- fluid
- solid
- time step
- state parameter
- field state
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a fluid-solid interface decoupling algorithm, a device and equipment, which are applied to the technical field of fluid-solid coupling mechanics, wherein the method comprises the steps of initializing a flow field and a solid structure field based on a flow field state parameter and a solid structure field state parameter of the last time step; calculating flow field state parameters and solid structure field state parameters in the current time step of the initialized flow field and solid structure field according to a preset algorithm; and decoupling the fluid-solid interface according to the flow field state parameter and the solid structure field state parameter in the current time step based on a third-order Dragon-Gregory tower iterative algorithm to obtain a state value of the fluid-solid interface of the current time step.
Description
Technical Field
The invention relates to the technical field of fluid-solid coupling mechanics, in particular to a fluid-solid interface decoupling algorithm, a fluid-solid interface decoupling device and fluid-solid interface decoupling equipment.
Background
Fluid-solid interactions are a physical phenomenon in almost all areas of fluid mechanics, marine engineering and aerospace. The impact of the fluid on the solids and the resulting response are of great significance for the design optimization of the engineering equipment. Thus, to date, there are a number of numerical algorithms applied to simulate fluid-solid interactions, including mainly any of the Lagrangian-Eulerian (ALE) algorithms and the two-way loose coupling algorithms.
When the Euler Lagrangian algorithm is applied to fluid-solid coupling problem simulation, the grid is required to be cut and reconstructed at a fluid-solid interface, so that the problem of low efficiency exists in calculation.
The bidirectional loose coupling algorithm considers the response of the solid to the fluid, but the response calculation of the fluid and the solid is not in one time step, so that the algorithm has the defect of obvious calculation accuracy in the case that the transient deformation is large due to the impact of the solid.
Disclosure of Invention
In view of the foregoing, the present invention provides a fluid-solid interface decoupling algorithm, apparatus and device, which overcomes at least some of the problems associated with the related art.
In order to solve the technical problems, the invention adopts the following technical scheme:
in a first aspect, a fluid-solid interface decoupling algorithm includes:
initializing a flow field and a solid structure field based on the flow field state parameter and the solid structure field state parameter of the last time step;
calculating flow field state parameters and solid structure field state parameters in the current time step of the initialized flow field and solid structure field according to a preset algorithm;
and decoupling the fluid-solid interface according to the flow field state parameter and the solid structure field state parameter in the current time step based on a third-order Dragon-Gregory tower iterative algorithm to obtain a state value of the fluid-solid interface of the current time step.
Optionally, the initializing the flow field and the solid structure field based on the flow field state parameter and the solid structure field state parameter of the previous time step includes:
acquiring a solid structure field state parameter of the last time step and a field state parameter of the last time step;
and inputting the solid structure field state parameter and the flow field state parameter into a grid which is generated at a fluid-solid interface in advance.
Optionally, the calculating the flow field state parameter and the solid structure field state parameter in the current time step of the initialized flow field and the initialized solid structure field according to a preset algorithm includes:
calculating flow field state parameters of the initialized flow field in the current time step according to a preset fluid control equation;
and calculating the solid structure field state parameters of the initialized solid structure field in the current time step according to a preset solid control equation.
Optionally, the preset solid control equation includes:
wherein x is 1 And x 2 Respectively represent M 1 And M 2 Variation of (u) 1 And u 2 Respectively represent M 1 And M 2 Is a velocity response of (2); m is m 1 And m 2 Represents M 1 And M 2 Mass, p I And p atm Respectively representing the pressure and the atmospheric pressure of the fluid-solid interface, l 0 The initial length of the spring in the equilibrium state is represented by k, the spring rate of the spring, and g, the gravitational acceleration.
Optionally, the preset fluid control equation is a control equation of the fluid under the euler coordinate system, and specifically includes:
wherein u= [ ρ, ρu, E] T ,F(U)＝[ρu,ρu 2 +p,(E+p)u] T Indicating the density of the fluid, u indicating the velocity of the fluid, p indicating the pressure of the fluid, e=e+0.5ρu 2 Representing total energy, e representing internal energy, t representing time, and x representing space.
Optionally, the third-order longge tower iterative algorithm includes:
U (1) ＝U n +ΔtL(U n )
U＝[x 1 ,x 2 ,p I ,u 1 ,u 2 ] T
wherein Δt represents a time step, ρ IL Representing the density at the left side of the fluid-solid interface, c IL Representing the sound velocity to the left of the fluid-solid interface and n represents the last time step.
Optionally, the fluid-solid structural field comprises: two mass solids and a spring connecting the two mass solids;
the solid structure field state parameters comprise deformation and speed response of two mass solids and deformation of the spring.
Optionally, the flow field state parameters include: the density of the fluid, the velocity of the fluid, the pressure of the fluid, the total energy of the fluid.
In a second aspect, a fluid-solid interface decoupling apparatus includes:
the initialization module is used for initializing the flow field and the solid structure field based on the flow field state parameter and the solid structure field state parameter of the last time step;
the parameter calculation module is used for calculating flow field state parameters and solid structure field state parameters in the current time step of the initialized flow field and solid structure field according to a preset algorithm;
and the decoupling module is used for decoupling the fluid-solid interface according to the flow field state parameter and the solid structure field state parameter in the current time step based on a third-order Dragon-Gregory tower iterative algorithm to obtain the state value of the fluid-solid interface of the current time step.
In a third aspect, a fluid-solid interface decoupling apparatus, comprising:
a processor, and a memory coupled to the processor;
the memory is used for storing a computer program;
the processor is configured to invoke and execute the computer program in the memory to perform the fluid-solid interface decoupling algorithm of the first aspect.
In a fourth aspect, a storage medium stores a computer program which, when executed by a processor, implements a fluid-solid interface decoupling algorithm according to any one of the first aspects of the present invention.
The invention adopts the technical scheme, and can realize the following technical effects: firstly initializing a flow field and a solid structure field based on flow field state parameters and solid structure field state parameters of the previous time step, and then calculating flow field state parameters and solid structure field state parameters in the current time step of the initialized flow field and solid structure field according to a preset algorithm; and decoupling the fluid-solid interface according to the flow field state parameter and the solid structure field state parameter in the current time step based on a third-order Dragon-Gregory tower iterative algorithm to obtain a state value of the fluid-solid interface of the current time step. Therefore, after the flow field and the solid field parameters of the previous time step are initialized, the parameters of the current time step are calculated, and the grid of the fluid-solid interface is not required to be cut and reconstructed, so that the calculation efficiency is improved; in addition, in the same time step, the state parameter values of the current time step flow field and the solid structure field are calculated firstly, and then the state value of the fluid-solid interface is calculated, so that the calculation error is reduced, and the calculation accuracy is improved.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of a fluid-solid interface decoupling algorithm according to an embodiment of the present invention;
FIG. 2 is a flow chart of a fluid-solid interface decoupling algorithm according to an embodiment of the present invention;
FIG. 3 is a flow chart of a fluid-solid interface decoupling algorithm according to another embodiment of the present invention;
FIG. 4 shows the experimental results and results of the present inventionA comparison verification result diagram of the reference solution;
FIG. 5 is a schematic structural view of a fluid-solid interface decoupling apparatus according to an embodiment of the present invention;
fig. 6 is a schematic structural diagram of a fluid-solid interface decoupling apparatus according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail below. It will be apparent that the described embodiments are only some, but not all, embodiments of the invention. All other embodiments, based on the examples herein, which are within the scope of the invention as defined by the claims, will be within the scope of the invention as defined by the claims.
For a better understanding of the solution provided in this application, it is necessary to understand the following:
in the fluid-solid coupling problem, fluid-solid coupling interface processing and data transmission determine success or failure of fluid-solid coupling problem calculation. The interaction of the fluid and the solid is mainly reflected at the moving interface, and when the interface moves and changes, the fluid-solid coupling phenomenon is generated. The fluid-solid coupling interface is a boundary between two different media, unlike boundary problems in solid mechanics or fluid mechanics.
Examples
In this application, one-dimensional algorithms for analyzing complete coupling of compressible fluids in elastic solid structures are implemented, and are described in summary based on the following cases.
Referring to fig. 1, a plane wave impacts an elastic solid consisting of two masses connected by springs in water to a ship. This case constitutes the most simplified problem of one-dimensional compressible fluid interacting with elastic solids. Wherein M1 and M2 represent two elastic solids, K is the elastic coefficient of the spring, u 1 And u 2 Respectively represent M 1 And M 2 Is provided.
Fig. 2 is a flow chart of a fluid-solid interface decoupling algorithm according to an embodiment of the present invention. As shown in fig. 2, the present embodiment provides a fluid-solid interface decoupling algorithm, including:
step 201, initializing a flow field and a solid structure field based on a flow field state parameter and a solid structure field state parameter of the last time step;
in some embodiments, the flow field state parameter and the solid structure field state parameter of the previous time step may be obtained in the calculation of the previous time step, and if the previous time step is the initial time step, the state parameter of the previous time step becomes the initial parameter. Initializing a flow field and a solid structure field, wherein the initialization comprises the steps of obtaining solid structure field state parameters of the last time step and flow field state parameters of the last time step; and inputting the solid structure field state parameter and the flow field state parameter into a grid which is generated at a fluid-solid interface in advance.
Step 202, calculating flow field state parameters and solid structure field state parameters in the current time step of the initialized flow field and solid structure field according to a preset algorithm;
in some embodiments, the preset algorithm includes a preset fluid control equation and a preset solid control equation. Specifically, step 202 includes: calculating flow field state parameters of the initialized flow field in the current time step according to a preset fluid control equation; and calculating the solid structure field state parameters of the initialized solid structure field in the current time step according to a preset solid control equation.
And 203, decoupling the fluid-solid interface based on a third-order Dragon-Gregory tower iterative algorithm according to the flow field state parameter and the solid structure field state parameter in the current time step to obtain a state value of the fluid-solid interface of the current time step.
In the method, after the flow field and the solid field parameters of the previous time step are initialized, the parameters of the current time step are calculated, and the grid of the fluid-solid interface is not required to be cut and reconstructed, so that the calculation efficiency is improved; in addition, in the same time step, the state parameter values of the current time step flow field and the solid structure field are calculated firstly, and then the state value of the fluid-solid interface is calculated, so that the calculation error is reduced, and the calculation accuracy is improved. And by completely decoupling the fluid-solid interface, iterative computation is performed on fluid computation and solid response in one time step, so that complete fluid-solid interaction coupling computation is achieved. The method can greatly improve the calculation accuracy, especially in the case that the solid is subjected to strong impact of fluid to generate instantaneous large deformation. When the solid is subjected to strong impact, large deformation can be generated in a time step, if the deformation exceeds a grid length, the calculation is not converged, and meanwhile, under the condition of large deformation, the flow field and the change of the structural field can be accurately calculated only through iterative calculation of fluid and solid responses. In the related art, basically all fluid-solid coupling algorithms cannot effectively solve such problems, especially the physical problems of severe response of shock waves generated when the structure is subjected to shock waves and cavitation collapse. The application provides an algorithm structure and a theoretical basis under the one-dimensional problem, and is a core algorithm for solving the actual engineering problem.
Fig. 3 is a schematic flow chart of a fluid-solid interface decoupling algorithm according to another embodiment of the present invention. As shown in fig. 3, the present embodiment provides a fluid-solid interface decoupling algorithm, including:
step 301, acquiring a solid structure field state parameter of a previous time step and a field state parameter of the previous time step;
based on the above related embodiments, the fluid-solid structured field of the present application comprises: two mass solids and a spring connecting the two mass solids. The solid structure field state parameters comprise deformation and speed response of two mass solids and deformation of the spring; the flow field state parameters include: the density of the fluid, the velocity of the fluid, the pressure of the fluid, the total energy of the fluid.
Step 302, inputting the solid structure field state parameter and the flow field state parameter into a grid which is generated at a fluid-solid interface in advance;
in some embodiments, grid generation is to divide a continuous space area into a limited number of subareas, and determine node distribution in each subarea, where the quality of the grid directly affects the accuracy and convergence of fluid-solid coupling calculation. Generally, grid generation can be classified into a structured grid method and an unstructured grid method. In the structured grid method, the connection relation between each node and the adjacent point is fixed, the connection between the nodes and the adjacent point is not required to be confirmed by setting data, and the number of grids is easy to control. In the unstructured grid method, the numbers of the cells and the nodes have no rule. In this embodiment, the grid may be generated by a generation method in the prior art, which is not limited herein.
Step 303, calculating flow field state parameters of the initialized flow field in the current time step according to a preset fluid control equation;
in some embodiments, the preset fluid control equation is a control equation of the fluid in the euler coordinate system, and specifically includes:
wherein u= [ ρ, ρu, E] T ,F(U)＝[ρu,ρu 2 +p,(E+p)u] T Indicating the density of the fluid, u indicating the velocity of the fluid, p indicating the pressure of the fluid, e=e+0.5ρu 2 Representing total energy, e representing internal energy, t representing time, and x representing space.
And step 304, calculating the solid structure field state parameters of the initialized solid structure field in the current time step according to a preset solid control equation.
In some embodiments, the predetermined solids control equation comprises:
wherein x is 1 And x 2 Respectively represent M 1 And M 2 Variation of (u) 1 And u 2 Respectively represent M 1 And M 2 Is a velocity response of (2); m is m 1 And m 2 Represents M 1 And M 2 Mass, p I And p atm The pressure and the atmospheric pressure of the fluid-solid interface are indicated, respectively. l (L) 0 The initial length of the spring in the equilibrium state is represented by k, the spring rate of the spring, and g, the gravitational acceleration.
Step 305, decoupling the fluid-solid interface based on the third-order Dragon-Gregory tower iterative algorithm according to the flow field state parameter and the solid structure field state parameter in the current time step, and obtaining the state value of the fluid-solid interface in the current time step.
Specifically, the third-order longge tower iterative algorithm includes:
U (1) ＝U n +ΔtL(U n )
U＝[x 1 ,x 2 ,p I ,u 1 ,u 2 ] T
wherein Δt represents a time step, ρ IL Representing the density at the left side of the fluid-solid interface, c IL Representing the sound velocity to the left of the fluid-solid interface and n represents the last time step.
Further, embodiments of the present invention may be described by the following case, as shown in fig. 1, in which the inside of the tube is compressible water.
The calculation region is a one-dimensional problem, where m 2 /m 1 ＝5,p max ＝161.5barL＝0.6345m,m 1 ＝8.678kg,m 2 ＝43.392kg,k＝42825kg/s 2 ,c f ＝1500m/s,p 0 ＝1.0bar,ρ 0 ＝1025kg/m 3 。
In this implementation, 1200 uniform grids are distributed over the calculation region (-10 m,5.08 m), with the fluid-solid interface at x=0.0. The calculation time is 150ms.
FIG. 4 shows the experimental results and results of the present inventionThe comparative verification of the baseline solution is schematically shown in fig. 4, where the fluid undergoes cavitation and collapses twice, which results in a multiple change of direction of velocity, positive-negative-zero. The results show that the velocity response and the historical data of cavitation zone size calculated by the present application are +.>The solution is basically identical. The first cavitation collapse (closed) spatiotemporal information captured herein is (x, t) = (1.25 m,49 ms), and the second cavitation collapse (closed) spatiotemporal information is (x, t) = (0.12 m,121 ms). Meanwhile, the result shows that the calculation speed of the technology exceeds +.>The model is more than 10 times.
The technology can calculate the parameter values of the flow field and the solid structure field under the interaction of fluid and solid through the flow, and the technology carries out simulation analysis on the problem of multiple physical fields through a complete two-way coupling method. Is an effective technique for researching strong fluid impact and solid strong response. Besides the fluid-solid interaction complete coupling algorithm of the patent, a scheme comprising ALE, unidirectional coupling and bidirectional sparse coupling can realize fluid-solid coupling simulation analysis to a certain extent. However, in the engineering problem that the solid is subjected to strong impact to generate instantaneous large deformation, the technical scheme of the patent is the scheme with the most reliability, the most stability and the highest precision.
Fig. 5 is a schematic structural diagram of a fluid-solid interface decoupling apparatus according to an embodiment of the present invention. As shown in fig. 5, the present embodiment provides a fluid-solid interface decoupling apparatus, including:
an initialization module 501, configured to initialize a flow field and a solid structure field based on a flow field state parameter and a solid structure field state parameter of a previous time step;
the parameter calculation module 502 is configured to calculate a flow field state parameter and a solid structure field state parameter in a current time step of the initialized flow field and solid structure field according to a preset algorithm;
and the decoupling module 503 is configured to decouple the fluid-solid interface according to the flow field state parameter and the solid structure field state parameter in the current time step based on a third-order Dragon-Gregory tower iterative algorithm, so as to obtain a state value of the fluid-solid interface in the current time step.
The specific implementation of this embodiment may refer to the related descriptions in the fluid-solid interface decoupling algorithm and the method embodiment described in the foregoing embodiments, which are not repeated herein.
Fig. 6 is a schematic structural diagram of a fluid-solid interface decoupling apparatus according to an embodiment of the present invention. Referring to fig. 6, an embodiment of the present application provides a fluid-solid interface decoupling apparatus, including:
a processor 601 and a memory 602 coupled to the processor;
the memory 602 is used to store a computer program;
the processor 601 is operative to invoke and execute a computer program in the memory 602 to perform the fluid-solid interface decoupling algorithm as in the embodiments described above.
The specific implementation of this embodiment may refer to the related descriptions in the fluid-solid interface decoupling algorithm and the method embodiment described in the foregoing embodiments, which are not repeated herein.
The embodiment of the invention provides a storage medium, which stores a computer program, and when the computer program is executed by a processor, the steps in a fluid-solid interface decoupling algorithm are realized.
For a specific implementation of this embodiment, reference may be made to the description related to the above embodiments of the fluid-solid interface decoupling algorithm, which is not repeated herein.
It is to be understood that the same or similar parts in the above embodiments may be referred to each other, and that in some embodiments, the same or similar parts in other embodiments may be referred to.
It should be noted that in the description of the present invention, the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. Furthermore, in the description of the present invention, unless otherwise indicated, the meaning of "plurality" means at least two.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and further implementations are included within the scope of the preferred embodiment of the present invention in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present invention.
It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution device. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
Those of ordinary skill in the art will appreciate that all or a portion of the steps carried out in the method of the above-described embodiments may be implemented by a program to instruct related hardware, where the program may be stored in a computer readable storage medium, and where the program, when executed, includes one or a combination of the steps of the method embodiments.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing module, or each unit may exist alone physically, or two or more units may be integrated in one module. The integrated modules may be implemented in hardware or in software functional modules. The integrated modules may also be stored in a computer readable storage medium if implemented in the form of software functional modules and sold or used as a stand-alone product.
The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, or the like.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.
Claims (6)
1. A method of fluid-solid interface decoupling comprising:
initializing a flow field and a solid structure field based on the flow field state parameter and the solid structure field state parameter of the last time step;
calculating the flow field state parameters of the initialized flow field in the current time step according to a preset fluid control equation;
calculating solid structure field state parameters of the initialized solid structure field in the current time step according to a preset solid control equation;
decoupling the fluid-solid interface according to the flow field state parameter and the solid structure field state parameter in the current time step based on a third-order Dragon-Gregory tower iterative algorithm to obtain a state value of the fluid-solid interface of the current time step;
wherein the preset solid control equation includes:
；
；
wherein,and->Respectively indicate->And->Is (are) modified>And->Respectively indicate->And->Is a velocity response of (2); />Andrepresentation->And->Quality of (1)>And->Representing two elastic solids, < >>And->Respectively the pressure and the atmospheric pressure of the fluid-solid interface, < ->Indicating the initial length of the spring in the state of equilibrium, < >>Indicating the spring rate of the spring, +.>Representing gravitational acceleration;
the preset fluid control equation includes:
；
wherein,,/>，/>indicating the density of the fluid, u indicating the velocity of the fluid, p indicating the pressure of the fluid, +.>Representing total energy, e representing internal energy, t representing time, x representing space;
the third-order Dragon-Gregory tower iterative algorithm comprises the following steps:
wherein,represents a time step,/->Represents the density on the left side of the fluid-solid interface, +.>Representing the sound velocity to the left of the fluid-solid interface and n represents the last time step.
2. The fluid-solid interface decoupling method of claim 1, wherein initializing the flow field and the solid structure field based on the flow field state parameter and the solid structure field state parameter of the last time step comprises:
acquiring a solid structure field state parameter of the last time step and a field state parameter of the last time step;
and inputting the solid structure field state parameter and the flow field state parameter into a grid which is generated at a fluid-solid interface in advance.
3. The fluid solid interface decoupling method of claim 1 or 2, wherein the fluid solid structural field comprises: two mass solids and a spring connecting the two mass solids;
the solid structure field state parameters comprise deformation and speed response of two mass solids and deformation of the spring.
4. The fluid solid interface decoupling method of claim 1 or 2, wherein the flow field state parameters comprise: the density of the fluid, the velocity of the fluid, the pressure of the fluid, the total energy of the fluid.
5. A fluid-solid interface decoupling apparatus, comprising:
the initialization module is used for initializing the flow field and the solid structure field based on the flow field state parameter and the solid structure field state parameter of the last time step;
the parameter calculation module is used for calculating the flow field state parameter of the initialized flow field at the current time step according to a preset fluid control equation;
calculating solid structure field state parameters of the initialized solid structure field in the current time step according to a preset solid control equation;
the decoupling module is used for decoupling the fluid-solid interface according to the flow field state parameter and the solid structure field state parameter in the current time step based on a third-order Dragon-Gregory tower iterative algorithm to obtain a state value of the fluid-solid interface of the current time step;
wherein the preset solid control equation includes:
；
；
wherein,and->Respectively indicate->And->Is (are) modified>And->Respectively indicate->And->Is a velocity response of (2); />Andrepresentation->And->Quality of (1)>And->Representing two elastic solids, < >>And->Respectively the pressure and the atmospheric pressure of the fluid-solid interface, < ->Indicating the initial length of the spring in the state of equilibrium, < >>Indicating the spring rate of the spring, +.>Representing gravitational acceleration;
the preset fluid control equation includes:
；
wherein,,/>，/>indicating the density of the fluid, u indicating the velocity of the fluid, p indicating the pressure of the fluid, +.>Representing total energy, e representing internal energy, t representing time, x representing space;
the third-order Dragon-Gregory tower iterative algorithm comprises the following steps:
wherein,represents a time step,/->Represents the density on the left side of the fluid-solid interface, +.>Representing the sound velocity to the left of the fluid-solid interface and n represents the last time step.
6. A fluid-solid interface decoupling apparatus, comprising:
a processor, and a memory coupled to the processor;
the memory is used for storing a computer program;
the processor is configured to invoke and execute the computer program in the memory to perform the fluid-solid interface decoupling method of any of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010498716.4A CN111783277B (en) | 2020-06-04 | 2020-06-04 | Fluid-solid interface decoupling algorithm, device and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010498716.4A CN111783277B (en) | 2020-06-04 | 2020-06-04 | Fluid-solid interface decoupling algorithm, device and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111783277A CN111783277A (en) | 2020-10-16 |
| CN111783277B true CN111783277B (en) | 2024-04-12 |
Family
ID=72753645
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010498716.4A Active CN111783277B (en) | 2020-06-04 | 2020-06-04 | Fluid-solid interface decoupling algorithm, device and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111783277B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102913464A (en) * | 2012-09-05 | 2013-02-06 | 江苏大学 | Method for predicting transient fluid-solid coupling characteristic of centrifugal pump rotor |
| NO20140079A1 (en) * | 2014-01-24 | 2015-07-27 | Ledaflow Technologies Da | Method for Determining Flow Parameters for a Multiphase Flow in a Pipeline Section |
| CN109902376A (en) * | 2019-02-25 | 2019-06-18 | 北京理工大学 | A high-precision numerical simulation method of fluid-structure interaction based on continuum mechanics |
| CN110765651A (en) * | 2019-11-12 | 2020-02-07 | 中国电子科技集团公司第二十九研究所 | Modeling and intelligent design method of microstrip direct coupling filter |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7430500B2 (en) * | 2002-12-27 | 2008-09-30 | Riken | Method and device for numerical analysis of flow field of incompressible viscous fluid, directly using V-CAD data |
| US8214185B2 (en) * | 2005-05-13 | 2012-07-03 | Seiko Epson Corporation | Stability performance of the coupled algorithms for viscoelastic ink jet simulations |
| CN102203782B (en) * | 2010-09-09 | 2013-03-13 | 天津空中代码工程应用软件开发有限公司 | Numerical method for solving eulerian equation of lagrange type |
-
2020
- 2020-06-04 CN CN202010498716.4A patent/CN111783277B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102913464A (en) * | 2012-09-05 | 2013-02-06 | 江苏大学 | Method for predicting transient fluid-solid coupling characteristic of centrifugal pump rotor |
| NO20140079A1 (en) * | 2014-01-24 | 2015-07-27 | Ledaflow Technologies Da | Method for Determining Flow Parameters for a Multiphase Flow in a Pipeline Section |
| CN109902376A (en) * | 2019-02-25 | 2019-06-18 | 北京理工大学 | A high-precision numerical simulation method of fluid-structure interaction based on continuum mechanics |
| CN110765651A (en) * | 2019-11-12 | 2020-02-07 | 中国电子科技集团公司第二十九研究所 | Modeling and intelligent design method of microstrip direct coupling filter |
Non-Patent Citations (2)
| Title |
|---|
| Application of a coupled Eulerian–Lagrangian method to simulate interactions between deformable composite structures and compressible multiphase flow;Wenfeng Xie et al;《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING》;1497-1519 * |
| Simulations of Compressible Two-Medium Flow by Runge-Kutta Discontinuous GalerkinMethodswith the Ghost Fluid Method;Jianxian Qiu;《COMMUNICATIONS IN COMPUTATIONAL PHYSICS》;479-504 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111783277A (en) | 2020-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Song | The scaled boundary finite element method in structural dynamics | |
| CN114757070B (en) | A new WENO format construction method for numerical simulation based on trigonometric function framework | |
| Witteveen et al. | An adaptive stochastic finite elements approach based on Newton–Cotes quadrature in simplex elements | |
| CN109740182A (en) | A kind of mesh free physical deformation emulation mode based on regeneration nuclear particle | |
| CN112001109A (en) | Realization of Structural Impact Dynamics Simulation Method by Regenerative Kernel Particle Algorithm | |
| Nichita et al. | A level set method coupled with a volume of fluid method for modeling of gas-liquid interface in bubbly flow | |
| Funes et al. | An efficient dynamic formulation for solving rigid and flexible multibody systems based on semirecursive method and implicit integration | |
| CN111783277B (en) | Fluid-solid interface decoupling algorithm, device and equipment | |
| Camara | A method for coupled arch dam‐foundation‐reservoir seismic behaviour analysis | |
| Emamzadeh et al. | Dynamic adaptive finite element analysis of acoustic wave propagation due to underwater explosion for fluid-structure interaction problems | |
| CN111801672A (en) | System and method for designing and manufacturing stable lattice structures | |
| Colella | Volume-of-fluid methods for partial differential equations | |
| Liu et al. | Numerical simulation of fluid-structure interaction using modified ghost fluid method and Naviers equations | |
| Fanelli et al. | Fluid-structure interaction problem of a deformable lamina solved with an original OpenFOAM code | |
| Yagawa | Free Mesh Method: fundamental conception, algorithms and accuracy study | |
| CN116595748B (en) | A high-precision numerical simulation method for strong shock compressible multiphase fluids | |
| Park et al. | Output-adaptive tetrahedral cut-cell validation for sonic boom prediction | |
| Bigarella et al. | A study of convective flux computation schemes for aerodynamic flows | |
| CN119106597B (en) | Ultrasonic suspension simulation method and device, electronic equipment and storage medium | |
| DOLEJŠÍ et al. | Anisotropic mesh adaptation for transonic and supersonic flow simulation | |
| Duan et al. | Recent progresses on a meshless Euler solver for compressible flows | |
| Cerroni et al. | Fluid structure interaction solver coupled with volume of fluid method for two-phase flow simulations | |
| Zeng et al. | A systematic procedure for achieving higher-order spatial accuracy in ghost fluid and other embedded boundary methods for fluid-structure interaction problems | |
| Kochi et al. | Shock capturing with discontinuous Galerkin Method using Overset grids for two-dimensional Euler equations | |
| Kaszynski et al. | Harmonic Convergence Estimation Through Strain Energy Superconvergence |
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 |