CN115526091A - Separated coupling numerical simulation method and device for multi-physical-field application - Google Patents
Separated coupling numerical simulation method and device for multi-physical-field application Download PDFInfo
- Publication number
- CN115526091A CN115526091A CN202211469835.2A CN202211469835A CN115526091A CN 115526091 A CN115526091 A CN 115526091A CN 202211469835 A CN202211469835 A CN 202211469835A CN 115526091 A CN115526091 A CN 115526091A
- Authority
- CN
- China
- Prior art keywords
- coupling
- physical
- region
- solid
- library
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application relates to a separated coupling numerical simulation method and a separated coupling numerical simulation device for multi-physical-field application. The method comprises the following steps: in numerical simulation of multi-physical-field application of flow-solid coupling, a fluid area is solved through a lattice Boltzmann method, a solid area is solved through a finite volume method, a consistent data mapping scheme provided by an open source coupling library is adopted to perform data interpolation and mapping on a discrete space at the coupling area, and parameters and physical field information of the coupling area are transmitted in real time through the open source coupling library based on time step; and carrying out iterative solution on the discrete equation of the subarea according to a coupling scheme package provided by the open source coupling library to obtain coupling response information applied by the multiple physical fields. The invention combines a macroscopic finite volume method and a mesoscopic lattice Boltzmann method, realizes the real-time interaction of multi-region parameters and physical information in the iterative solving process, greatly facilitates the analysis of coupled multi-physical field application, and can provide efficient and reliable reference for engineering practice.
Description
Technical Field
The application relates to the field of multi-physical field numerical simulation, in particular to a separated coupling numerical simulation method and device for multi-physical field application.
Background
In essence, the real world is a world with multiple physical fields coupled, and the coupled physical fields are affected with each other, which is complicated and complicated. The discovery and research of multi-physics applications suggests that people are reconsidering the world from a perspective that is deeper and closer to the natural essence. With the development of science and technology, in the process of researching the application of multiple physical fields, the traditional research method based on observation and experiment and constructed on the basis of the thought of single physical field analysis cannot deal with the interaction between complex coupled multiple physical fields. Sometimes, it is not a difficult matter to obtain experimental results, but it is a real challenge how to give convincing theoretical explanations. In recent years, the numerical simulation technology has the advantages of low cost, capability of reproducing special experimental conditions and the like, and has been developed into the visual field of people. At present, the combination of experimental research and numerical simulation technology, whether scientific research or product development, is in the trend.
The numerical simulation technology applied to the multi-physical field is not only to simply superpose the simulation of a plurality of single physical fields, but in the field of numerical simulation, the coupling of the multi-physical field can be understood from three angles of multi-field, multi-region and multi-scale. Wherein, the multi-field refers to the excitation and response of a plurality of physical fields existing in the system at the same time; multi-region refers to a system in which multiple continuum with different characteristics interact directly through boundaries; multi-scale refers to the continuous span of behavior from micro-to macro-scale in a system. Common multi-physical field couplings are fluid-solid couplings, electrical-thermal couplings, thermal-structural couplings, thermal-electrical-structural couplings, acoustic-structural couplings, fluid-thermal couplings, and the like.
At present, numerical simulations for multi-physics applications are mainly classified into a direct coupling method and a split coupling method according to a solving method. The separation type coupling method comprises the steps of firstly analyzing each physical field related to multi-physical-field application, then transferring a certain physical field in a coupling area to other physical fields through interpolation and mapping for calculation, and obtaining the response of a coupling system through iterative solution.
At present, respective numerical algorithms have been developed in various scientific and engineering fields, but in a fluid-solid coupling scenario, different algorithms cannot realize real-time coupling of multiple physical fields, which greatly limits the application range of multi-physical field numerical simulation. Although the macro scale numerical algorithms can ensure the conservation characteristic of the system and easily process irregular boundaries, the algorithms need to continuously calculate global residuals and correct intermediate results, so the algorithms are not efficient and have poor expandability. The micro-scale numerical algorithm can fully describe the system details, but the simulation scale is limited to the atomic and molecular level, the simulation time step is small, and therefore the simulation process needs great computing performance. The mesoscopic algorithm with the scale between the two is simple, good in locality and easy to parallelize, but due to the introduction of statistical errors, the precision of the simulation is affected.
Disclosure of Invention
In view of the foregoing, there is a need to provide a separated coupling numerical simulation method, apparatus, computer device and storage medium for multi-physical-field application, which can implement real-time coupling of multiple physical fields in a flow-solid coupling scenario.
A split coupled numerical simulation method for multi-physics applications, the method comprising:
establishing a corresponding example shape according to a sub-region applied by the multi-physics field of the fluid-solid coupling, and setting a boundary condition; the sub-region comprises a fluid region and a solid region, and a coupling region is formed by surface contact between the fluid region and the solid region;
constructing a discrete equation of the fluid region by a lattice Boltzmann method, constructing a discrete equation of the solid region by a finite volume method, and determining a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring in the sub-region;
when the multi-physical-field application is subjected to coupling simulation according to the first solver and the second solver, real-time transmission is carried out on parameters and physical field information of the coupling area through a preset open source coupling library based on a time step;
and iteratively solving the discrete equation of the fluid region and the discrete equation of the solid region according to a coupling scheme package provided by the open source coupling library until convergence, so as to obtain coupling response information applied to the multi-physics field.
In one embodiment, the method further comprises: aligning the different discrete spaces involved in the coupling region.
In one embodiment, the method further comprises: determining physical parameters involved in a coupling simulation process;
and converting the representation forms of the physical parameters in the lattice boltzmann method and the finite volume method according to equivalent condition information of the physical parameters in the lattice boltzmann method and the finite volume method so as to ensure parameter equivalence of the physical parameters in two types of solving algorithms.
In one embodiment, the method further comprises: the preset open source coupling library is an open source coupling library precCE.
In one embodiment, the method further comprises: and determining the time steps of the lattice Boltzmann method and the finite volume method according to the requirements of the time interpolation scheme provided by the open source coupling library, and establishing the equivalent correlation of the time steps of the two algorithms.
In one embodiment, the method further comprises: selecting a corresponding numerical coupling mode from a coupling scheme package provided by the open source coupling library according to the characteristics of the multi-physical-field application and the requirement on numerical solving precision;
and carrying out iterative solution on the discrete equation of the fluid region and the discrete equation of the solid region according to the selected numerical coupling mode until convergence, and obtaining coupling response information applied by the multi-physical field.
In one embodiment, the method further comprises: the multi-physical field application of the flow-solid coupling is a conjugate heat transfer multi-physical field.
A split coupled numerical simulation apparatus for multi-physics applications, the apparatus comprising:
the initialization module is used for establishing a corresponding example shape according to a sub-area applied by the flow-solid coupled multi-physical field and setting a boundary condition; the sub-region comprises a fluid region and a solid region, and a coupling region is formed by surface contact between the fluid region and the solid region;
a solver determination module, configured to construct a discrete equation of the fluid region by a lattice boltzmann method, construct a discrete equation of the solid region by a finite volume method, determine a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring within the partitioned region;
the real-time interaction module is used for transmitting the parameters and the physical field information of the coupling area in real time based on a time step through a preset open source coupling library when the coupling simulation is carried out on the multi-physical field application according to the first solver and the second solver;
and the iterative solution module is used for iteratively solving the discrete equation of the fluid region and the discrete equation of the solid region according to the coupling scheme package provided by the open source coupling library until convergence, so as to obtain the coupling response information applied by the multi-physical field.
A computer device comprising a memory storing a computer program and a processor implementing the following steps when the computer program is executed:
establishing corresponding example shapes according to the subareas applied by the flow-solid coupled multi-physical field, and setting boundary conditions; the sub-region comprises a fluid region and a solid region, and a coupling region is formed by surface contact between the fluid region and the solid region;
constructing a discrete equation of the fluid region by a lattice Boltzmann method, constructing a discrete equation of the solid region by a finite volume method, and determining a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring in the sub-region;
when the multi-physical-field application is subjected to coupling simulation according to the first solver and the second solver, real-time transmission is carried out on parameters and physical field information of the coupling area through a preset open source coupling library based on a time step;
and carrying out iterative solution on the discrete equation of the fluid region and the discrete equation of the solid region according to a coupling scheme package provided by the open source coupling library until convergence, so as to obtain coupling response information applied by the multi-physical field.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
establishing corresponding example shapes by regions according to the multi-physical field application of flow-solid coupling, and setting boundary conditions; the sub-region comprises a fluid region and a solid region, and a coupling region is formed by surface contact between the fluid region and the solid region;
constructing a discrete equation of the fluid region by a lattice boltzmann method, constructing a discrete equation of the solid region by a finite volume method, and determining a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring in the sub-region;
when the multi-physical-field application is subjected to coupling simulation according to the first solver and the second solver, parameters and physical field information of the coupling area are transmitted in real time through a preset open-source coupling library based on a time step;
and iteratively solving the discrete equation of the fluid region and the discrete equation of the solid region according to a coupling scheme package provided by the open source coupling library until convergence, so as to obtain coupling response information applied to the multi-physics field.
In the separated coupling numerical simulation method, the separated coupling numerical simulation device, the computer equipment and the storage medium for the multi-physical-field application, in the numerical simulation of the flow-solid coupling multi-physical-field application, a fluid area is solved by a lattice Boltzmann method, a solid area is solved by a finite volume method, for the coupling area, a consistent data mapping scheme provided by an open source coupling library is adopted to interpolate and map data in a discrete space at the coupling area, and in the coupling numerical simulation process, parameters and physical field information of the coupling area are transmitted in real time based on time step length through a preset open source coupling library; and carrying out iterative solution on the discrete equation of the subarea according to a coupling scheme package provided by the open source coupling library until convergence, and obtaining the coupling response information applied by the multi-physical field of the flow-solid coupling. According to the characteristics of multiple physical fields and multiple regions, a macroscopic finite volume method and a mesoscopic lattice Boltzmann method are combined, real-time interaction of parameters and physical information of the multiple regions is realized in an iterative solving process, and coupling of a non-gridding method and a structured grid method in a numerical simulation process of the multiple physical fields is realized, so that the analysis of coupling multiple physical field application is greatly facilitated, and efficient and reliable reference can be provided for engineering practice; in addition, in the simulation process of the fluid, compared with the traditional finite volume method, the lattice boltzmann method has higher efficiency and better expandability, and is more suitable for simulation solution on a high-performance computer.
Drawings
FIG. 1 is a diagram of an application scenario of a split coupling numerical simulation approach for multi-physics applications in an embodiment;
FIG. 2 is a schematic illustration of the geometry and boundary conditions for a conjugate heat transfer coupling multiphysics application in an example embodiment;
FIG. 3 is a schematic diagram illustrating spatial dispersion of a conjugate heat transfer coupling multi-physics field application in one embodiment;
FIG. 4 is a diagram illustrating a consistency mapping scheme provided by precCE in an exemplary embodiment;
FIG. 5 is a schematic diagram illustrating a dimensional transformation of a design in an embodiment;
FIG. 6 is a schematic diagram of a time stepping scheme in an embodiment;
FIG. 7 is a schematic diagram of a block Cheng Ouge embodiment;
FIG. 8 is a schematic diagram of an implicit coupling scheme in an embodiment;
FIG. 9 is a block diagram of a split coupled numerical simulation apparatus for multi-physics-oriented applications according to an embodiment;
FIG. 10 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
In one embodiment, as shown in fig. 1, a separated coupling numerical simulation method for multi-physics applications is provided, comprising the following steps:
102, establishing corresponding example shapes according to the multi-physics field application subarea of the fluid-solid coupling, and setting boundary conditions.
The sub-region comprises a fluid region and a solid region, and the fluid region and the solid region are in surface contact to form a coupling region, so that the coupling region is formed by surface coupling.
The coupling of multiple physical fields can be understood from three angles of multiple fields, multiple areas and multiple scales, and at present, the corresponding numerical coupling mode is determined based on the characteristics of different physical fields in different areas and the requirement of numerical simulation precision.
The application of the multi-physical field of the fluid-solid coupling in the embodiment is the application of the conjugate heat transfer coupling multi-physical field. In industry and life, the conjugated heat transfer process is widely existed. The process of heat transfer between a solid and a contacting fluid is studied by conjugate heat transfer. The solids are heat transferred by a conduction process; while fluid heat transfer is dominated by convection. And establishing corresponding example shapes according to the selected conjugate heat transfer coupling multi-physical field application subarea, and applying corresponding boundary conditions on the boundaries of the examples.
And 104, constructing a discrete equation of the fluid region by a lattice Boltzmann method, constructing a discrete equation of the solid region by a finite volume method, and determining a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring in the region.
The numerical simulation process is based on spatial dispersion, and continuous simulation areas need to be spatially dispersed according to different numerical algorithm requirements before simulation begins. The fluid area is solved by the lattice Boltzmann method and the solid area is solved by the finite volume method by combining the characteristics of the macroscopic finite volume method and the mesoscopic lattice Boltzmann method.
Generally, a finite volume method or a finite difference method is usually selected to develop a numerical solver to simulate the motion condition of a flow field, and the deformation condition of a solid is simulated by the finite element method. For the selected application of the conjugated heat transfer multi-physical field, based on the basic characteristics of the physical field, the solid heat transfer part is calculated by a Laplace solver based on a finite volume method, the fluid part is solved by selecting a lattice Boltzmann method, and the coupling of the non-lattice Boltzmann method and the finite volume method based on the structured lattice in the numerical simulation process of the multi-physical field is realized by combining the non-lattice Boltzmann method and the finite volume method based on the structured lattice. In addition, in the simulation process of the fluid, compared with the traditional finite volume method, the lattice boltzmann method has higher efficiency and better expandability, and is more suitable for simulation solution on a high-performance computer.
And 106, when the coupling simulation is carried out on the multi-physical-field application according to the first solver and the second solver, real-time transmission is carried out on the parameters and the physical field information of the coupling area through a preset open source coupling library based on the time step length.
Aiming at multi-region coupling multi-physical field application, the invention combines a finite volume method and a lattice Boltzmann method, respectively carries out dispersion based on different spatial positions, in the coupling solving process, a solver of a fluid region and a solid region needs to transmit the spatial dispersion condition of a coupling region to an open source coupling library according to requirements, and the transmission of physical field data at the coupling region is completed through a consistency mapping scheme provided by the open source coupling library.
Specifically, the open source coupled library of the present embodiment is preCICE, and the physical field data at the coupled region, which is transferred by the interpolation and mapping scheme provided by preCICE, is temperature and heat flux.
And 108, iteratively solving the discrete equation of the fluid region and the discrete equation of the solid region according to the coupling scheme provided by the open source coupling library until convergence, so as to obtain coupling response information applied by the multiple physical fields.
When a finite volume method and a lattice boltzmann method jointly simulate a coupled multi-physical field application, time step lengths of two simulation algorithms are determined according to requirements of a time interpolation scheme provided by a precCE based on time step transmission parameters and physical field information, and equivalent correlation of the time step lengths of the two simulation algorithms is established.
The coupling scheme of the numerical equation determines the coupling mode of each separation solver, the precICE provides many different coupling schemes for the separation solver, including explicit and implicit, serial and parallel and the like, and four possible combination modes include serial explicit, serial implicit, parallel explicit and parallel implicit and the like. And determining a corresponding numerical coupling scheme according to the selected application characteristics of the conjugated heat transfer coupling multi-physical field and the requirement of the application on numerical solving precision, wherein the coupling scheme comprises data information corresponding to the coupling scheme provided for the open source coupling library.
In the separated coupling numerical simulation method for multi-physical-field application, in the numerical simulation of the multi-physical-field application of flow-solid coupling, a fluid area is solved by a lattice boltzmann method, a solid area is solved by a finite volume method, for the coupling area, a consistent data mapping scheme provided by an open source coupling library is adopted to interpolate and map data in a discrete space at the coupling area, and in the coupling numerical simulation process, parameters and physical field information of the coupling area are transmitted in real time on the basis of time step length through a preset open source coupling library; and carrying out iterative solution on the discrete equation of the subarea according to a coupling scheme package provided by the open source coupling library until convergence, and obtaining the coupling response information applied by the multi-physical field of the flow-solid coupling. According to the characteristics of multiple physical fields and multiple regions, a macroscopic finite volume method and a mesoscopic lattice Boltzmann method are combined, real-time interaction of parameters and physical information of the multiple regions is realized in an iterative solving process, and coupling of a non-gridding method and a structured grid method in a numerical simulation process of the multiple physical fields is realized, so that the coupling analysis of the application of the multiple physical fields is greatly facilitated, and efficient and reliable reference can be provided for engineering practice; in addition, in the simulation process of the fluid, compared with the traditional finite volume method, the lattice boltzmann method has higher efficiency and better expandability, and is more suitable for simulation solution on a high-performance computer.
In one embodiment, the method further comprises: the different discrete spaces involved in the coupling region are aligned.
Since the data mapping scheme provided by the open source coupling library precce has corresponding requirements on the discrete space at the coupling region, different discrete spaces involved in the coupling region need to be aligned first, and then consistency of data mapping and interpolation of the coupling region is ensured.
In one embodiment, the method further comprises: determining physical parameters involved in a coupling simulation process; and converting the representation forms of the physical parameters in the lattice boltzmann method and the finite volume method according to equivalent condition information of the physical parameters in the lattice boltzmann method and the finite volume method so as to ensure parameter equivalence of the physical parameters in two solving algorithms.
The finite volume method and the lattice boltzmann method respectively belong to two different numerical solving mechanisms of macroscopic continuous flow and mesoscopic dynamics, the meanings and the expression forms of physical parameters in the two solving systems are different, the solving processes based on the lattice boltzmann and the finite volume method in the coupling simulation process are respectively calculated aiming at respective dimension systems, but the two solving algorithms need to exchange physical field information at a coupling area. Therefore, the representation forms of the basic physical parameters in different numerical algorithms need to be converted according to the equivalent conditions of the parameters in the two methods, so as to ensure the parameter equivalence in the finite volume method and the lattice boltzmann solution method.
It should be understood that, although the steps in the flowchart of fig. 1 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least a portion of the steps in fig. 1 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternately with other steps or at least a portion of the sub-steps or stages of other steps.
In one embodiment, the coupling response of a coupled thermal coupling multiphysics application is numerically simulated by the method proposed by the present invention.
FIG. 2 is a graphical representation of the geometry and boundary conditions for selected applications of the embodiment of the conjugate heat transfer coupling multiphysics. The application is as follows: at a temperature of
At a certain flow rate
Flows through a length of
A bottom of the plate is heated
Heating (
) According to the research, the thickness of the flat plate can influence the temperature conduction process, and the thickness of the flat plate is selected to be the thickness
. The fluid part establishes an example shape and applies boundary conditions through a Palabos platform based on a lattice Boltzmann method, and the solid part establishes an example shape as shown in the figure and applies boundary conditions through an OpenFOAM platform based on a finite volume method.
FIG. 3 is a schematic diagram illustrating the spatial dispersion of selected applications of the coupled multi-physics field. The fluid part adopts a lattice boltzmann method to carry out discrete solution, the solid part adopts a finite volume method to carry out discrete solution, and a square box in the figure marks an applied coupling interface. In the coupling solution process, the solver of the fluid and solid regions needs to transfer the spatial discrete condition of the coupling region to preCICE as required, and the transfer of the physical field (temperature and heat flux) data at the coupling region is completed through a consistency mapping scheme provided by the preCICE as shown in fig. 4.
Fig. 5 is a schematic diagram of dimension conversion designed by the embodiment. The simulation system may be classified into a physical system, a dimensionless system, and a Lattice Boltzmann (LBM) system according to a type difference of the physical parameter dimension. A solver based on a finite volume method carries out numerical simulation on a physical process in a physical or dimensionless system, and the solver based on a lattice Boltzmann method carries out simulation on the physical process of the lattice Boltzmann system. In the coupling simulation process, two solvers need to exchange physical field information. The calculation process of the lattice boltzmann needs to use a temperature field calculated by a finite volume method at a coupling interface to solve the problem of the finite volume methodThe solution requires the use of a heat flux calculated by lattice boltzmann at the coupling interface. The unit conversion process from the finite volume method to the lattice Boltzmann method can be divided into two stages, namely a dimensionless process from a physical system to a dimensionless system and a scaling process from the dimensionless system to the lattice Boltzmann system. Characteristic length from physical system to dimensionless system
Characteristic time
Specific density of
Characteristic temperature of
On the basis of discrete physical quantities of length from dimensionless system to lattice Boltzmann system
Time discrete physical quantity
And temperature discrete physical quantity
Based on the fact that the physical parameters can be unified by the Reynolds numbers in different systems
And prandtl number
To perform the engraving.
Fig. 6 is a schematic diagram of the time stepping scheme of the present embodiment. preCICE provides two time stepping schemes: the fixed time window and the first participant specify a time step. The present invention selects a fixed time window as shown in fig. 6A step-by-step scheme. In a fixed time window, each participating solver solves at a time step no greater than the time window, and the coupled data is communicated only when the time window is completed each time. After each solver finishes a time step, the used time step is sent to the preconCE, and the preconCE returns the time from the end of the current time window (the beginning of the next time window)
Each solver calculates the next time step, if greater than
E.g. dark grey part of fig. 6, then use
(of dark grey portion in FIG. 6)
) As the next time step; if less than
The time step calculated using the solver (FIG. 6 light grey part), as in FIG. 6 light grey part
) The solution is performed as the next time step.
FIG. 7 is a schematic diagram of the equation coupling scheme of this example. preCICE provides some different coupling schemes for numerical equations. According to the selected conjugate heat transfer multi-physical field application characteristics, a serial implicit coupling scheme is determined. In the implicit coupling scheme, the output of one solver is used as the input of the other solver, and an iterative algorithm is called to solve at a fixed point until the numerical solution converges. The Acc part in fig. 7 refers to the fixed-point iterative algorithm. In a serial implicit coupling scheme:
both equations use
The value of the time, the first equation using
After the completion of the sub-iteration
The second equation uses the first
After the iteration is completed
。
The converged serial scheme is as follows:
the general form is as follows:
wherein the content of the first and second substances,
，
(Picard iteration).
As shown in FIG. 8, in the implicit coupling scheme, a new time window is entered, and the solutions for coupling are solvedThe device will first save the state of the current physical field. If it is used
And
without convergence, precce is recalculated using the new boundary values, but the values of the internal fields have been affected, and thus the previously saved physical field states are reloaded and iteration continues with the new boundary values until convergence. The precICE provides a plurality of methods for solving Picard iteration, and in the embodiment, an IQN-ILS scheme is selected for iterative computation.
In one embodiment, as shown in fig. 9, there is provided a separated coupling numerical simulation apparatus for multi-physics applications, comprising: an initialization module 902, a solver determination module 904, a real-time interaction module 906, and an iterative solution module 908, wherein:
an initialization module 902, configured to establish a corresponding example shape according to a sub-region applied by a flow-solid coupled multi-physics field, and set a boundary condition; the sub-regions comprise a fluid region and a solid region, and the fluid region and the solid region are in surface contact to form a coupling region;
a solver determining module 904, configured to construct a discrete equation of the fluid region by a lattice boltzmann method, construct a discrete equation of the solid region by a finite volume method, and determine a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring in the partitioned region;
the real-time interaction module 906 is configured to transmit parameters of the coupling area and physical field information in real time based on a time step through a preset open source coupling library when the coupling simulation is performed on the multi-physical field application according to the first solver and the second solver;
and the iterative solution module 908 is configured to iteratively solve the discrete equations of the fluid region and the discrete equations of the solid region according to the coupling scheme provided by the open source coupling library until convergence, so as to obtain coupling response information applied to the multi-physical field.
The real-time interaction module 906 is also used to align the different discrete spaces involved within the coupling region.
The real-time interaction module 906 is further configured to determine physical parameters involved in the coupling simulation process; and converting the representation forms of the physical parameters in the lattice boltzmann method and the finite volume method according to equivalent condition information of the physical parameters in the lattice boltzmann method and the finite volume method so as to ensure that the parameters of the physical parameters are equivalent in two types of solving algorithms.
The real-time interaction module 906 is further configured to determine time step lengths of the lattice boltzmann method and the finite volume method according to a requirement of a time interpolation scheme provided by the open source coupling library, and establish an equivalent correlation of the time step lengths of the two algorithms.
The real-time interaction module 906 is further configured to select a corresponding numerical coupling mode from a coupling scheme package provided by the open source coupling library according to the characteristics of the multi-physics field application and the requirement on numerical solution accuracy; and carrying out iterative solution on the discrete equation of the fluid region and the discrete equation of the solid region according to the selected numerical coupling mode until convergence, and obtaining coupling response information applied by the multiple physical fields.
For specific limitations of the split coupling numerical simulation apparatus for multi-physical-field application, refer to the above limitations of the split coupling numerical simulation method for multi-physical-field application, and are not described herein again. The modules in the separated coupling numerical simulation device for multi-physical-field application can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 10. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a split coupling numerical simulation method for multi-physics applications. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, a key, a track ball or a touch pad arranged on a shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
It will be appreciated by those skilled in the art that the configuration shown in fig. 10 is a block diagram of only a portion of the configuration associated with the present application, and is not intended to limit the computing device to which the present application may be applied, and that a particular computing device may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory storing a computer program and a processor implementing the steps of the method embodiments described above when the processor executes the computer program.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which, when executed by a processor, performs the steps of the above-described method embodiment.
It will be understood by those skilled in the art that all or part of the processes of the method of the embodiments described above may be implemented by hardware instructions of a computer program, which may be stored in a non-volatile computer-readable storage medium, and when executed, may include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), rambus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A separated coupling numerical simulation method for multi-physics application is characterized by comprising the following steps:
establishing corresponding example shapes according to the subareas applied by the flow-solid coupled multi-physical field, and setting boundary conditions; the sub-region comprises a fluid region and a solid region, and a coupling region is formed by surface contact between the fluid region and the solid region;
constructing a discrete equation of the fluid region by a lattice Boltzmann method, constructing a discrete equation of the solid region by a finite volume method, and determining a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring in the sub-region;
when the multi-physical-field application is subjected to coupling simulation according to the first solver and the second solver, real-time transmission is carried out on parameters and physical field information of the coupling area through a preset open source coupling library based on a time step;
and carrying out iterative solution on the discrete equation of the fluid region and the discrete equation of the solid region according to a coupling scheme package provided by the open source coupling library until convergence, so as to obtain coupling response information applied by the multi-physical field.
2. The method of claim 1, further comprising, before real-time transferring the parameters and physical field information of the coupling region based on a time step through a preset open source coupling library, the method further comprising:
aligning the different discrete spaces involved in the coupling region.
3. The method of claim 1, further comprising, before real-time transferring the parameters and physical field information of the coupling region based on a time step through a preset open source coupling library, the method further comprising:
determining physical parameters involved in a coupling simulation process;
and converting the representation forms of the physical parameters in the lattice boltzmann method and the finite volume method according to equivalent condition information of the physical parameters in the lattice boltzmann method and the finite volume method so as to ensure parameter equivalence of the physical parameters in two kinds of solving algorithms.
4. The method of claim 1, wherein the preset open-source coupled library is an open-source coupled library preCICE.
5. The method of claim 4, before real-time transferring the parameters and the physical field information of the coupling region based on the time step through a preset open source coupling library, further comprising:
and according to the requirement of a time interpolation scheme provided by the open source coupling library, determining the time step of the lattice boltzmann method and the finite volume method, and establishing the equivalent correlation of the time step of the two algorithms.
6. The method of claim 4, wherein iteratively solving the discrete equations for the fluid region and the discrete equations for the solid region according to a coupling scheme package provided by the open-source coupling library until convergence to obtain coupling response information for the multi-physics application comprises:
selecting a corresponding numerical coupling mode from a coupling scheme package provided by the open source coupling library according to the characteristics of the multi-physical-field application and the requirement on numerical solving precision;
and carrying out iterative solution on the discrete equation of the fluid region and the discrete equation of the solid region according to the selected numerical coupling mode until convergence, and obtaining coupling response information applied by the multi-physical field.
7. The method according to any one of claims 1 to 6, wherein the application of the multi-physical field of flow-solid coupling is a conjugate heat transfer multi-physical field.
8. A split coupled numerical simulation apparatus for multi-physics applications, the apparatus comprising:
the initialization module is used for establishing corresponding example shapes in a partitioned mode according to the multi-physical field application of flow-solid coupling and setting boundary conditions; the sub-region comprises a fluid region and a solid region, and a coupling region is formed by surface contact between the fluid region and the solid region;
a solver determination module, configured to construct a discrete equation of the fluid region by a lattice boltzmann method, construct a discrete equation of the solid region by a finite volume method, determine a first solver of the discrete equation of the fluid region and a second solver of the discrete equation of the solid region according to a physical process occurring within the partitioned region;
the real-time interaction module is used for transmitting the parameters and the physical field information of the coupling area in real time through a preset open source coupling library based on a time step when the multi-physical-field application is subjected to coupling simulation according to the first solver and the second solver;
and the iterative solution module is used for iteratively solving the discrete equation of the fluid region and the discrete equation of the solid region according to a coupling scheme package provided by the open source coupling library until convergence, so as to obtain the coupling response information applied to the multi-physics field.
9. The apparatus of claim 8, wherein the real-time interaction module is further configured to:
aligning the different discrete spaces involved in the coupling region.
10. The apparatus of claim 8, wherein the real-time interaction module is further configured to:
determining physical parameters involved in a coupling simulation process;
and converting the representation forms of the physical parameters in the lattice boltzmann method and the finite volume method according to equivalent condition information of the physical parameters in the lattice boltzmann method and the finite volume method so as to ensure parameter equivalence of the physical parameters in two types of solving algorithms.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211469835.2A CN115526091B (en) | 2022-11-22 | 2022-11-22 | Separated coupling numerical simulation method and device for multi-physics field application |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211469835.2A CN115526091B (en) | 2022-11-22 | 2022-11-22 | Separated coupling numerical simulation method and device for multi-physics field application |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115526091A true CN115526091A (en) | 2022-12-27 |
| CN115526091B CN115526091B (en) | 2023-03-17 |
Family
ID=84705026
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211469835.2A Active CN115526091B (en) | 2022-11-22 | 2022-11-22 | Separated coupling numerical simulation method and device for multi-physics field application |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115526091B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004092992A (en) * | 2002-08-30 | 2004-03-25 | Mitsubishi Heavy Ind Ltd | Device and method for calculating heat flux of furnace heat transfer tube |
| CN102592016A (en) * | 2011-12-31 | 2012-07-18 | 英特工程仿真技术（大连）有限公司 | Engineering multiphysics coupling analysis method |
| CN104866695A (en) * | 2015-06-24 | 2015-08-26 | 武汉大学 | GPU-accelerated fluid-structure coupling simulation method through immersion boundary and lattice Boltzmann methods |
| CN111695309A (en) * | 2020-06-02 | 2020-09-22 | 上海科技大学 | High-performance large-scale fluid-solid coupling fluid simulation method based on statistical dynamics |
| CN113946905A (en) * | 2021-09-13 | 2022-01-18 | 华南理工大学 | CFD-FEM-SPH four-way coupled carrier-liquid ship hydro-elastic response simulation method |
| FR3114669A1 (en) * | 2020-09-25 | 2022-04-01 | Safran | Method for simulating the flow of a fluid in contact with a mobile solid |
| CN115238611A (en) * | 2022-08-09 | 2022-10-25 | 南京航空航天大学 | Multiphase flow simulation efficiency optimization method based on multiphase lattice Boltzmann flux method |
-
2022
- 2022-11-22 CN CN202211469835.2A patent/CN115526091B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004092992A (en) * | 2002-08-30 | 2004-03-25 | Mitsubishi Heavy Ind Ltd | Device and method for calculating heat flux of furnace heat transfer tube |
| CN102592016A (en) * | 2011-12-31 | 2012-07-18 | 英特工程仿真技术（大连）有限公司 | Engineering multiphysics coupling analysis method |
| CN104866695A (en) * | 2015-06-24 | 2015-08-26 | 武汉大学 | GPU-accelerated fluid-structure coupling simulation method through immersion boundary and lattice Boltzmann methods |
| CN111695309A (en) * | 2020-06-02 | 2020-09-22 | 上海科技大学 | High-performance large-scale fluid-solid coupling fluid simulation method based on statistical dynamics |
| FR3114669A1 (en) * | 2020-09-25 | 2022-04-01 | Safran | Method for simulating the flow of a fluid in contact with a mobile solid |
| CN113946905A (en) * | 2021-09-13 | 2022-01-18 | 华南理工大学 | CFD-FEM-SPH four-way coupled carrier-liquid ship hydro-elastic response simulation method |
| CN115238611A (en) * | 2022-08-09 | 2022-10-25 | 南京航空航天大学 | Multiphase flow simulation efficiency optimization method based on multiphase lattice Boltzmann flux method |
Non-Patent Citations (8)
| Title |
|---|
| LIU Y , GUO X , LI C , ET AL.: "The Communication-Overlapped Hybrid Decomposition Parallel Algorithm for Multi-Scale Fluid Simulations", 《 THE 48TH INTERNATIONAL CONFERENCE ON PARALLEL PROCESSING》 * |
| SEN ZHANG, XIAOWEI GUO, CHAO LI, YI LIU, RAN ZHAO, CANQUN YANG: "Numerical Study of Fluid-Structure Interaction Dynamics under High-explosive Detonation on Massively Parallel Computers", 《IEEE 22ND INTERNATIONAL CONFERENCE ON HIGH PERFORMANCE COMPUTING AND COMMUNICATIONS》 * |
| SEN ZHANG, XIAO-WEI GUO, CHAO LI, YI LIU, SIJIANG FAN, RAN ZHAO, CANQUN YANG: "A large scale parallel fluid‐structure interaction computing platform for simulating structural responses to a detonation shock" * |
| 曾非同等: "基于多尺度多物理场的油浸式变压器流动-传热数值研究", 《电工技术学报》 * |
| 栾辉宝等: "有限体积法与LBM分区耦合模拟方腔自然对流", 《西安交通大学学报》 * |
| 栾辉宝等: "有限容积法与格子Boltzmann方法耦合模拟传热流动问题", 《科学通报》 * |
| 段宝岩等: "电子装备机电热多场耦合问题初探", 《电子机械工程》 * |
| 皮兴才等: "有限体积格子Boltzmann方法用于近空间连续流区绕流模拟", 《载人航天》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115526091B (en) | 2023-03-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108573089B (en) | System and method for automatically creating fluid model by using 3D CAD point cloud | |
| Bove et al. | Modeling solid oxide fuel cell operation: Approaches, techniques and results | |
| Geyer et al. | Automated metamodel generation for Design Space Exploration and decision-making–A novel method supporting performance-oriented building design and retrofitting | |
| US20220027539A1 (en) | Method and system for manufacturing a heat exchanger for supercritical pressure fluid | |
| Zhu et al. | Model reductions for multiscale stochastic optimization of cooling water system equipped with closed wet cooling towers | |
| Bobič et al. | Dynamic behaviour of a plate heat exchanger: Influence of temperature disturbances and flow configurations | |
| CN114077796A (en) | High-adaptability multi-phase particle dispersion type fuel element temperature field calculation method | |
| CN105760586A (en) | Fuzzy temperature response subordinating degree function solving method based on point collocation theory | |
| CN115659908A (en) | Three-unit unbalanced porous medium method of printed circuit board heat exchanger | |
| CN116595884A (en) | Method for realizing ultra-real-time calculation of miniature nuclear reactor system | |
| CN115526091B (en) | Separated coupling numerical simulation method and device for multi-physics field application | |
| CN112230679B (en) | Group coupling system cooperative control method and device based on time delay | |
| Cao et al. | Reducio: Model reduction for data center predictive digital twins via physics-guided machine learning | |
| Sun et al. | An intelligent plate fin-and-tube heat exchanger design system through integration of CFD, NSGA-II, ANN and TOPSIS | |
| CN105677995A (en) | Method for numerical solution of fuzzy steady state heat conduction problem based on full grid point collocation theory | |
| Tao et al. | Lattice Boltzmann simulation of complex thermal flows via a simplified immersed boundary method | |
| CN116882253A (en) | Modelica-based thermal member modeling method, device and medium | |
| CN112347638B (en) | Dual-unit method based three-dimensional integrated micro-system electric-thermal coupling analysis method | |
| Sharma et al. | Surrogate modeling approach towards coupling computational fluid dynamics and energy simulations for analysis and design of energy efficient attics | |
| Chen et al. | Design methodology of large-scale thermoelectric generation: A hierarchical modeling approach | |
| Reuter et al. | Computational models for predicting cooling tower fill performance in cross-counterflow configuration | |
| Hillary et al. | Guidelines for developing efficient thermal conduction and storage models within building energy simulations | |
| Pejman et al. | Topology optimization of microvascular composites for active-cooling applications using a geometrical reduced-order model | |
| CN113962130A (en) | Calculation method for three-dimensional contact heat transfer | |
| Bautista et al. | Transient heat conduction in a solid slab using multiple-scale analysis |
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 |