CN113688554A - Flow field calculation method combining discontinuous finite element and immersion boundary method - Google Patents
Flow field calculation method combining discontinuous finite element and immersion boundary method Download PDFInfo
- Publication number
- CN113688554A CN113688554A CN202111071637.6A CN202111071637A CN113688554A CN 113688554 A CN113688554 A CN 113688554A CN 202111071637 A CN202111071637 A CN 202111071637A CN 113688554 A CN113688554 A CN 113688554A
- Authority
- CN
- China
- Prior art keywords
- flow field
- unit
- eta
- wall surface
- wall
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Complex Calculations (AREA)
Abstract
The invention discloses a flow field calculation method combining an interrupted finite element and an immersed boundary method, which comprises the steps of S1, generating right-angle grid units for the whole flow field; s2, constructing an anisotropic shape function in the right-angle grid unit; s3, initializing a flow field; and S4, performing time-marching circulation on the flow field by using an immersion boundary method based on the constructed anisotropic shape function to obtain flow field data and realize flow field analysis. The method can greatly improve the flow numerical simulation efficiency of various complex geometric areas in the practical application of engineering, and has better application prospect.
Description
Technical Field
The invention belongs to the technical field of fluid mechanics, and particularly relates to a flow field calculation method combining an interrupted finite element method and an immersion boundary method.
Background
In recent years, immersion Boundary methods (Immersed Boundary methods) in the field of Computational Fluid Dynamics (CFD) have attracted much attention in the academic and application fields, and in the immersion Boundary methods, rectangular grids are used for the whole flow field calculation, and a finite difference Method discrete control Method is matched. In the vicinity of the solid wall surface boundary, aiming at the intersection condition of the solid surface and the rectangular grid, the immersion boundary method applies boundary conditions by adopting an extrapolation method or an equivalent external force method. The step of generating the grid is omitted, and the problem of coordinate conversion from a physical space to a calculation space does not need to be processed in the calculation, so that the calculation efficiency can be greatly improved, and the method becomes a hotspot of research in the field of CFD at present.
As shown in fig. 1, the direct mesh of the submerged boundary method cannot resolve the extremely thin boundary layer of the surface of an object in high reynolds number flow, since it is not possible to avoid normal encryption alone in a skin mesh. If the thin boundary layer is directly distinguished by encrypting the isotropic right-angle grid, which is equivalent to performing unnecessary encryption on the flow direction and the spreading direction while encrypting the normal grid, the grid number is increased sharply, the calculation amount is huge, and the numerical simulation can not be completed under reasonable calculation resources (calculation power, memory and time consumption). Therefore, the immersion boundary method is only used for simulating low Reynolds number flow initially, because the boundary layer of the flow is thick and can be resolved by a thicker right-angle grid, as shown in FIG. 2, and the high Reynolds number attached boundary layer is simulated later by combining the wall function technology, because the wall function technology abandons the direct resolution of the internal structure of the boundary layer by the grid and adopts an approximate model to estimate the influence of the boundary layer indirectly, thereby effectively reducing the requirement on the size of the grid. However, this approach actually circumvents the difficulties, and the immersion boundary method still fails in the convection field analysis.
Disclosure of Invention
Aiming at the defects in the prior art, the flow field calculation method combining the discontinuous finite element and the immersion boundary method solves the problem that the flow of high Reynolds number wall surfaces is difficult to flow respectively in the process of performing numerical simulation on the boundary layer structure by using the conventional boundary entering method.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: a flow field calculation method combining an interrupted finite element method and an immersion boundary method comprises the following steps:
s1, generating right-angle grid units for the whole flow field;
s2, constructing shape functions of all units in the flow field based on the generated right-angle grid units;
s3, initializing a flow field;
and S4, performing time advancing circulation on each unit of the flow field by using an immersion boundary method based on the shape function of each unit of the structure to obtain flow field calculation data.
Further, in step S1, in the process of generating the rectangular grid cells, the encryption process is performed near the wall surface.
Further, the step S2 is specifically:
s21, in the rectangular grid cells, marking the rectangular cells near the wall surface as near-wall cells, and defining a xi-eta coordinate system of the near-wall cells according to the direction of the local wall surface;
wherein xi is a wall surface tangential direction, and eta is a wall surface normal direction;
s22, respectively selecting polynomials of corresponding times in xi and eta directions in the near wall unit, constructing corresponding anisotropic shape functions, and forming an approximate space corresponding to a near wall unit solution, namely an approximate space of a discontinuous finite element;
and simultaneously, adopting an isotropic shape function of an x-y space in a non-near-wall unit to complete the shape function construction of all units in the flow field.
Further, in step S21, the ξ - η coordinate system of the near-wall unit is determined specifically as follows:
in the rectangular grid unit, for the rectangular unit of the coarse rectangular grid in the x-y space, a xi-eta coordinate system is set at the centroid of the unit, and the xi-eta is obtained by arbitrarily rotating the x-y coordinate system of the coordinate system.
Further, in step S22, when the gradient of the tangential velocity field in the vicinity of the wall surface in the wall surface normal direction is larger than the gradient of the wall surface tangential direction in the viscous flow, the degree of the η -direction corresponding polynomial is larger than the degree of the ξ -direction corresponding polynomial.
Further, the step S22 is specifically:
establishing a local coordinate system xi-eta by taking the unit centroid as an origin (0,0), and setting the distribution of the flow field variable u (xi, eta) in the unit to be approximately obtained by the linear combination of corresponding shape functions, namely:
in the formula, B(m,n)(xi, eta) is a shape function, the highest degree of xi direction is n times, the highest degree of eta direction is m times, the polynomial of xi is p times, the polynomial of eta is q times,coefficients that are shape functions;
the degree of freedom and the shape function of the (p +1) × (q +1) term pair are obtained based on u (xi, η) and are used as an approximate space of the near wall unit solution.
Further, the step S4 is specifically:
s41, applying physical boundary conditions;
s42, based on the applied physical boundary conditions, using the immersion boundary method to correct the freedom degree of the wall surface cut in the right-angle grid unit;
s43, calculating residual errors of discontinuous finite elements of the units according to the shape functions of the constructed units;
s44, performing time advancing by adopting a set time dispersion mode, and updating the degree of freedom of solution variables in the units according to the residual errors corresponding to each unit;
wherein each cell comprises a near-wall cell and a non-near-wall cell;
s45, calculating current flow field data according to the degree of freedom of solution variables in each current unit and outputting the data;
s46, repeating the steps S41-S45 until the set time step length is reached;
and S47, completing time advancing circulation, storing output flow field data and completing flow field calculation.
Further, the flow field data in the step S45 includes an aerodynamic coefficient, a fluid density, a fluid pressure, and a fluid temperature.
The invention has the beneficial effects that:
the method can greatly improve the flow numerical simulation efficiency of various complex geometric areas in practical engineering application, has good application prospect and is embodied in the following points:
(1) the degree of freedom configuration of the anisotropic shape function discontinuous finite element adopted in the invention does not depend on the shape of the grid, and has better flexibility;
(2) in the invention, the high Reynolds number boundary layer of the irregular wall surface can be really distinguished only by adopting the right-angle grids;
(3) the invention adopts the right-angle grid, which can greatly save the manpower, time and computing resources required by grid generation and can easily realize the step of grid automatic generation in CFD.
Drawings
Fig. 1 is a schematic diagram of a cylindrical skin grid structure grid (left) and a cylindrical cartesian grid (right) in the background art provided by the present invention.
Fig. 2 is a schematic diagram illustrating the boundary layer resolving power comparison between the conformal grid (left) and the rectangular grid (right) of the NACA0012 airfoil structure in the background art provided by the present invention.
FIG. 3 is a flow chart of a flow field calculation method combining discontinuous finite element and dipping boundary methods according to the present invention.
FIG. 4 is a schematic view of an anisotropic flow field under a coarse grid of right angles provided by the present invention.
FIG. 5 is a schematic diagram of a ξ - η coordinate system of a rectangular unit provided by the present invention.
Fig. 6 is a schematic diagram of a unit cut by a wall surface according to the present invention.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
As shown in fig. 3, a flow field calculation method combining discontinuous finite element and dipping boundary method includes the following steps:
s1, generating right-angle grid units for the whole flow field;
s2, constructing shape functions of all units in the flow field based on the generated right-angle grid units;
s3, initializing a flow field;
and S4, performing time advancing circulation on each unit of the flow field by using an immersion boundary method based on the shape function of each unit of the structure to obtain flow field calculation data.
In step S1 of the present embodiment, in the process of generating the rectangular grid cells, an encryption process is performed near the wall surface to improve the approximation to the irregular wall surface.
In step S2 of this embodiment, an anisotropic shape function discontinuous finite element is constructed, and decoupling of the degree of freedom configuration and the mesh is achieved, that is, the degree of freedom configuration no longer depends on the shape, size, and other factors of the mesh.
The degrees of freedom in the conventional immersion boundary method are arranged at the vertices of the rectangular meshes or the cell centers, and are tightly bound to the mesh shapes, so that since the meshes cannot decompose the wall geometry, the degrees of freedom defined by the coarse rectangular meshes cannot match the wall geometry to the boundary layer anisotropy resolution form. The invention combines the anisotropic shape function discontinuous finite element in the immersion boundary method, and can distinguish the irregular wall surface shape and the thin boundary layer in the coarse right-angle grid.
The discontinuous finite element of the anisotropic function in this embodiment is obtained by sampling any anisotropic shape function and improving on the basis of the taylor unfolding discontinuous finite element method, and based on this, the step S2 specifically includes:
s21, in the rectangular grid cells, marking the rectangular cells near the wall surface as near-wall cells, and defining a xi-eta coordinate system of the near-wall cells according to the direction of the local wall surface;
wherein xi is a wall surface tangential direction, and eta is a wall surface normal direction;
s22, respectively selecting polynomials of corresponding times in xi and eta directions in the near wall unit, constructing corresponding anisotropic shape functions, and forming an approximate space corresponding to a near wall unit solution, namely an approximate space of a discontinuous finite element;
and simultaneously, adopting an isotropic shape function of an x-y space in a non-near-wall unit to complete the shape function construction of all units in the flow field.
In step S21, the method for determining the ξ - η coordinate system of the near-wall unit specifically includes:
in the rectangular grid unit, for the rectangular unit of the coarse rectangular grid in the x-y space, a xi-eta coordinate system is set at the centroid of the unit, and the xi-eta is obtained by arbitrarily rotating the x-y coordinate system of the coordinate system.
As the shape function of the near-wall unit is based on the xi-eta coordinate system, the shape function and the corresponding freedom degree can be matched with the designated direction at will, and the decoupling of the freedom degree configuration and the rectangular grid is realized.
The step S22 is specifically:
establishing a local coordinate system xi-eta by taking the unit centroid as an origin (0,0), and setting the distribution of the flow field variable u (xi, eta) in the unit to be approximately obtained by the linear combination of corresponding shape functions, namely:
in the formula, B(m,n)(xi, eta) is a shape function, the highest degree of xi direction is n times, the highest degree of eta direction is m times, the polynomial of xi is p times, the polynomial of eta is q times,coefficients that are shape functions;
obtaining the degree of freedom and the anisotropic shape function of (p +1) × (q +1) term pairing based on u (xi, eta), taking the degree of freedom and the anisotropic shape function as an approximate space of a near-wall unit solution, and performing spatial dispersion on a flow field control equation; since p and q are not necessarily the same, the approximation capability of the approximation space of these shape functions to form the near-wall cell solution is different in both ξ and η directions, i.e., anisotropic for the ξ - η space, in order to better resolve the drastic changes in the normal direction within the wall bounding layer.
Based on the obtained anisotropic shape function discontinuous finite element, in a coarse right-angle grid of the immersion boundary method, the anisotropy characteristics of the irregular wall surface shape and the wall surface boundary layer in the right-angle grid can be distinguished through flexible configuration freedom.
As shown in fig. 4, in step S22 described above, when the gradient of the tangential velocity field in the vicinity of the wall surface in the wall surface normal direction is larger than the gradient of the wall surface tangential direction in the viscous flow, the degree of the η -direction corresponding polynomial expression is larger than the degree of the ξ -direction corresponding polynomial expression. When the variable has larger change in the tangential direction, the tangential polynomial degree can also be properly increased, in short, the polynomial degrees in different directions are self-adaptively adjusted according to the local flow field change to form an anisotropic shape function.
Specifically, the flow field near the wall surface has a larger gradient in the normal direction of the wall surface and a less variable flow field in the tangential direction. In order to analyze the flow field, enough degrees of freedom need to be arranged in the wall surface normal direction, for a right-angle unit close to the wall surface, the wall surface tangential direction and the wall surface normal direction can be selected as a shape function space, for a right-angle unit close to the wall surface, and the wall surface tangential direction and the wall surface normal direction can be selected as a xi-eta coordinate system of the shape function space. Meanwhile, in order to distinguish the drastic change of the flow field near the wall surface in the normal direction (eta direction), a higher-order polynomial can be adopted for the shape function in the eta direction, and a lower-order polynomial can be adopted for the shape function along the flow direction (xi direction) (as shown in fig. 5). In this way, the identification of the wall surface geometric characteristics and the analysis of the anisotropic boundary layer under the coarse right-angle grid can be realized.
Step S4 of this embodiment specifically includes:
s41, applying physical boundary conditions;
s42, based on the applied physical boundary conditions, using the immersion boundary method to correct the freedom degree of the wall surface cut in the right-angle grid unit;
as shown in fig. 6, for a cell cut by a wall surface (in the figure, a dark color is a flow field area, and a light color is a solid area), the value of the degree of freedom is corrected, so that the value of the flow area is almost unchanged, and at the same time, a function is ensured to meet boundary conditions, thereby avoiding boundary application in the sense of an immersion boundary method is completed, and specific embodiments refer to boundary measures of a conventional immersion method, such as cut-cell method, ghost-cell method, grid-less method, and the like;
s43, calculating residual errors of discontinuous finite elements of the units according to the shape functions of the constructed units;
s44, performing time advancing by adopting a set time dispersion mode, and updating the degree of freedom of solution variables in the units according to the residual errors corresponding to each unit;
wherein each cell comprises a near-wall cell and a non-near-wall cell;
s45, calculating current flow field data according to the degree of freedom of solution variables in each current unit and outputting the data;
s46, repeating the steps S41-S45 until the set time step length is reached;
and S47, completing time advancing circulation, storing output flow field data and completing flow field calculation.
The flow field data in the above step S45 includes an aerodynamic coefficient, a fluid density, a fluid pressure, and a fluid temperature.
In the flow field calculation process, the key point of distinguishing the high Reynolds number boundary layer is that enough degrees of freedom are arranged along the normal direction of the wall surface, so that the step length of the normal close-fit radius layer grid is very small, and the step lengths in other directions can be larger, so that the thin boundary layer with high anisotropy can be economically and effectively distinguished.
The invention provides a novel method for realizing flow field analysis by combining an anisotropic shape function and an immersion boundary method, which has the capability of realizing more freedom degree in the wall surface normal direction in a thicker right-angle grid so as to ensure that a high Reynolds number thin boundary layer is distinguished in a thick right-angle grid, and simultaneously reserves the advantages of simple generation and easy realization of automation of the right-angle grid in the immersion boundary method.
Claims (8)
1. A flow field calculation method combining an interrupted finite element method and an immersed boundary method is characterized by comprising the following steps:
s1, generating right-angle grid units for the whole flow field;
s2, constructing shape functions of all units in the flow field based on the generated right-angle grid units;
s3, initializing a flow field;
and S4, performing time advancing circulation on each unit of the flow field by using an immersion boundary method based on the shape function of each unit of the structure to obtain flow field calculation data.
2. The method for calculating a flow field according to claim 1, wherein in step S1, the method includes performing encryption processing near the wall surface during the process of generating the rectangular grid cells.
3. The method for calculating a flow field by combining discontinuous finite element and dipping boundary method according to claim 1, wherein the step S2 is specifically as follows:
s21, in the rectangular grid cells, marking the rectangular cells near the wall surface as near-wall cells, and defining a xi-eta coordinate system of the near-wall cells according to the direction of the local wall surface;
wherein xi is a wall surface tangential direction, and eta is a wall surface normal direction;
s22, respectively selecting polynomials of corresponding times in xi and eta directions in the near wall unit, constructing corresponding anisotropic shape functions, and forming an approximate space corresponding to a near wall unit solution, namely an approximate space of a discontinuous finite element;
and simultaneously, adopting an isotropic shape function of an x-y space in a non-near-wall unit to complete the shape function construction of all units in the flow field.
4. A flow field calculation method incorporating discontinuous finite element and dipping boundary method as claimed in claim 3, wherein in step S21, the xi- η coordinate system of the near wall element is determined by:
in the rectangular grid unit, for the rectangular unit of the coarse rectangular grid in the x-y space, a xi-eta coordinate system is set at the centroid of the unit, and the xi-eta is obtained by arbitrarily rotating the x-y coordinate system of the coordinate system.
5. The flow field calculation method combining the discontinuous finite element method and the dipping boundary method according to claim 3, wherein in the step S22, when the gradient of the tangential velocity field near the wall surface in the normal direction of the wall surface is larger than that in the tangential direction of the wall surface under the viscous flow condition, the degree of the η -direction corresponding polynomial is larger than that in the ξ -direction corresponding polynomial.
6. The method for calculating a flow field by combining an interrupted finite element method and an immersed boundary method according to claim 4, wherein the step S22 is specifically as follows:
establishing a local coordinate system xi-eta by taking the unit centroid as an origin (0,0), and setting the distribution of the flow field variable u (xi, eta) in the unit to be approximately obtained by the linear combination of corresponding shape functions, namely:
in the formula, B(m,n)(xi, eta) is a shape function, the highest degree of xi direction is n times, the highest degree of eta direction is m times, the polynomial of xi is p times, the polynomial of eta is q times,coefficients that are shape functions;
the degree of freedom and the shape function of the (p +1) × (q +1) term pair are obtained based on u (xi, η) and are used as an approximate space of the near wall unit solution.
7. The method for calculating a flow field by combining discontinuous finite element and dipping boundary method according to claim 5, wherein the step S4 is specifically as follows:
s41, applying physical boundary conditions;
s42, based on the applied physical boundary conditions, using the immersion boundary method to correct the freedom degree of the wall surface cut in the right-angle grid unit;
s43, calculating residual errors of discontinuous finite elements of the units according to the shape functions of the constructed units;
s44, performing time advancing by adopting a set time dispersion mode, and updating the degree of freedom of solution variables in the units according to the residual errors corresponding to each unit;
wherein each cell comprises a near-wall cell and a non-near-wall cell;
s45, calculating current flow field data according to the degree of freedom of solution variables in each current unit and outputting the data;
s46, repeating the steps S41-S45 until the set time step length is reached;
and S47, completing time advancing circulation, storing output flow field data and completing flow field calculation.
8. The method of claim 7, wherein the flow field data in step S45 includes an aerodynamic coefficient, a fluid density, a fluid pressure, and a fluid temperature.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111071637.6A CN113688554B (en) | 2021-09-14 | 2021-09-14 | Flow field calculation method combining discontinuous finite element and immersion boundary method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111071637.6A CN113688554B (en) | 2021-09-14 | 2021-09-14 | Flow field calculation method combining discontinuous finite element and immersion boundary method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113688554A true CN113688554A (en) | 2021-11-23 |
CN113688554B CN113688554B (en) | 2023-03-21 |
Family
ID=78586261
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111071637.6A Active CN113688554B (en) | 2021-09-14 | 2021-09-14 | Flow field calculation method combining discontinuous finite element and immersion boundary method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113688554B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115114689A (en) * | 2022-07-08 | 2022-09-27 | 中国航空发动机研究院 | Geometric figure boundary marking method and device, electronic equipment and readable storage medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103853921A (en) * | 2014-02-24 | 2014-06-11 | 昆明理工大学 | Method for predicting flow-induced vibration characteristic of large-deformation super-elastic structure |
CN103970989A (en) * | 2014-04-15 | 2014-08-06 | 昆明理工大学 | Immersing boundary flow field calculation method based on fluid/solid interface consistency |
CN111159853A (en) * | 2019-12-10 | 2020-05-15 | 北京航空航天大学 | Boundary layer high-precision processing method for high Reynolds number viscous flow problem |
CN113051634A (en) * | 2019-12-27 | 2021-06-29 | 河海大学 | Fluid-solid coupling numerical calculation method for bank slope landslide surge |
-
2021
- 2021-09-14 CN CN202111071637.6A patent/CN113688554B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103853921A (en) * | 2014-02-24 | 2014-06-11 | 昆明理工大学 | Method for predicting flow-induced vibration characteristic of large-deformation super-elastic structure |
CN103970989A (en) * | 2014-04-15 | 2014-08-06 | 昆明理工大学 | Immersing boundary flow field calculation method based on fluid/solid interface consistency |
CN111159853A (en) * | 2019-12-10 | 2020-05-15 | 北京航空航天大学 | Boundary layer high-precision processing method for high Reynolds number viscous flow problem |
CN113051634A (en) * | 2019-12-27 | 2021-06-29 | 河海大学 | Fluid-solid coupling numerical calculation method for bank slope landslide surge |
Non-Patent Citations (2)
Title |
---|
WEI-XI HUANG 等: "Three-dimensional simulation of elastic capsules in shear flow by the penalty immersed boundary method" * |
屈崑 等: "CFD数学模型的线性化方法及其应用" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115114689A (en) * | 2022-07-08 | 2022-09-27 | 中国航空发动机研究院 | Geometric figure boundary marking method and device, electronic equipment and readable storage medium |
CN115114689B (en) * | 2022-07-08 | 2024-03-12 | 中国航空发动机研究院 | Geometric border marking method, device, electronic equipment and readable storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN113688554B (en) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | A hybrid boundary node method | |
Nakahashi et al. | Building-cube method for large-scale, high resolution flow computations | |
Hansen et al. | Simulation of wake effects between two wind farms | |
Li et al. | Dynamic wind farm wake modeling based on a Bilateral Convolutional Neural Network and high-fidelity LES data | |
CN101887478A (en) | Sequence radial basis function agent model-based high-efficiency global optimization method | |
Gang et al. | Mesh deformation on 3D complex configurations using multistep radial basis functions interpolation | |
CN113688554B (en) | Flow field calculation method combining discontinuous finite element and immersion boundary method | |
CN114638048A (en) | Three-dimensional spray pipe flow field rapid prediction and sensitivity parameter analysis method and device | |
CN113076678B (en) | Frequency domain two-degree body weight abnormity rapid numerical simulation method and device | |
CN105224726B (en) | The method that structured grid Dynamic mesh is used for unstrctured grid flow field calculation device | |
CN117854643A (en) | MEMS membrane simulation method and system based on graphic neural network | |
Zhang et al. | Adaptation for hybrid unstructured grid with hanging node method | |
CN111339688B (en) | Method for solving rocket simulation model time domain equation based on big data parallel algorithm | |
CN110287632B (en) | Method for simulating numerical value of cathode contact resistance of solid oxide fuel cell | |
CN116191404A (en) | Distributed power generation prediction method and device, electronic equipment and storage medium | |
CN106650293B (en) | Underground water model evaluation method based on AM nested sampling algorithm | |
CN108233382B (en) | Method for extracting Jacobian matrix of rectangular coordinate tidal current equation | |
Ivan et al. | Hyperbolic conservation laws on three-dimensional cubed-sphere grids: a parallel solution-adaptive simulation framework | |
Malik et al. | Complex electrode microstructure simulations using a smoothed boundary method with adaptive mesh refinement | |
CN104037756B (en) | A kind of Power System Stability Assessment method containing Complex Power device model | |
Zhao et al. | Diffeomorphism Neural Operator for various domains and parameters of partial differential equations | |
CN107957974B (en) | Method for extracting polar coordinate tidal current equation Jacobi matrix | |
CN117077298B (en) | Aircraft robust optimization design method based on gradient enhancement random Co-Kriging model | |
El Marghichi et al. | Modelling photovoltaic modules with enhanced accuracy using particle swarm clustered optimization | |
Kim et al. | High-density mesh flow computations by building-cube method |
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 |