CN114722745A - Turbulent flow wall surface distance calculation method and device, computer equipment and storage medium - Google Patents
Turbulent flow wall surface distance calculation method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN114722745A CN114722745A CN202210650302.8A CN202210650302A CN114722745A CN 114722745 A CN114722745 A CN 114722745A CN 202210650302 A CN202210650302 A CN 202210650302A CN 114722745 A CN114722745 A CN 114722745A
- Authority
- CN
- China
- Prior art keywords
- coordinate
- wall surface
- parallel
- cartesian
- plane
- 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/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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Image Generation (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
The application relates to a turbulent wall surface distance calculation method, a device, computer equipment and a storage medium, wherein the method comprises the following steps: establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane; determining the number of wall surface units parallel to a Cartesian coordinate surface in a geometric object plane; determining a rotation coordinate axis according to the number of the parallel wall surface units corresponding to each coordinate surface; in a Cartesian three-dimensional rectangular coordinate system, rotating a geometric object plane by a preset angle around a rotation coordinate axis, and storing data of a wall surface unit of the rotated geometric object plane into an ADT data structure; in a Cartesian three-dimensional rectangular coordinate system, rotating a space point by a preset angle around a rotating coordinate axis, and establishing a search box for the temporary space point obtained by rotation; searching wall units falling into the range of the search box in the ADT data structure, and recording the wall units as target units; and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the turbulent wall surface distance. By adopting the method, the calculation efficiency of the turbulent wall surface distance can be improved.
Description
Technical Field
The present application relates to the technical field of computational fluid mechanics, and in particular, to a method and an apparatus for calculating a distance between turbulent walls, a computer device, and a storage medium.
Background
Computational Fluid Dynamics (CFD) is a cross discipline for performing numerical simulation analysis on Fluid mechanics problems by using computers and numerical algorithms, and the purpose of research is to obtain good application in numerous industrial fields represented by aerospace and solve numerous key aerodynamic problems in aerospace aircraft development. The real fluid motion has two flow states with different properties of laminar flow and turbulent flow, and the flow in nature is mostly turbulent flow. The use of CFD to simulate laminar flow has become mature, but faces the problem of turbulence in actual flow, and due to the physical high complexity, it can only be handled in a simplified manner. At present, three main methods for turbulence numerical simulation exist: direct Numerical Simulation (DNS), Large Eddy Simulation (LES), and Reynolds-average N-S equation method (RANS). In these turbulence simulation methods, the wall distance, which is a parameter critical to the simulation accuracy, exists, i.e., the minimum distance from the center of a grid cell to the surface of the aircraft is calculated.
A common geometric method for calculating the wall distance in engineering applications is the ADT method (Alternating Digital Tree, a special binary Tree data structure). Because the ADT method does not need to traverse all wall units when calculating the wall distance, the efficiency is extremely high, and the ADT method is widely applied. However, when the geometric object plane is parallel to the cartesian coordinate plane, employing ADT search will face a very embarrassing problem: since the geometric object plane is parallel to the cartesian coordinate plane, the boundaries of the geometric box of spatial points are parallel to the geometric object plane. When ADT search is carried out, if a given search radius r1 is too small, a wall surface unit cannot be searched; on the other hand, if the search radius r2 is enlarged, too many wall cells are searched. In both cases, the ADT method is degenerated to a direct search method, and the efficiency is greatly reduced.
Disclosure of Invention
In view of the above, there is a need to provide a turbulent wall distance method, apparatus, computer device and storage medium capable of improving the efficiency of calculating the wall distance in the case where the geometric object plane is parallel to the cartesian coordinate plane.
A method of turbulence wall distance calculation, the method comprising:
establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane;
determining the number of wall surface units of a coordinate surface parallel to the Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
determining a rotation coordinate axis according to the number of parallel wall surface units corresponding to each coordinate surface;
in the Cartesian three-dimensional rectangular coordinate system, the geometric object surface is rotated by a preset angle around the rotation coordinate axis, and the data of the wall surface unit of the rotated geometric object surface is stored in an ADT data structure;
in the Cartesian three-dimensional rectangular coordinate system, rotating the space points by the preset angle around the rotation coordinate axis, and establishing a search box for the temporary space points obtained by rotation;
searching the ADT data structure for the wall surface unit falling in the range of the search box, and recording the wall surface unit as a target unit;
and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
In one embodiment, the determining the number of wall elements in the geometric object plane parallel to the coordinate plane of the cartesian three-dimensional rectangular coordinate system includes: traversing each wall surface unit on the geometric object surface; judging whether the normal vector of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system or not; and if so, adding one to the number of the parallel wall surface units of the coordinate surface perpendicular to the coordinate axis to obtain the number of the parallel wall surface units corresponding to each coordinate surface.
In one embodiment, the determining the rotation coordinate axis according to the number of the parallel wall units corresponding to each coordinate plane includes: judging whether the number of parallel wall surface units corresponding to at least one coordinate surface is larger than a preset value in the number of parallel wall surface units corresponding to the three coordinate surfaces; if so, identifying a coordinate plane with the number of parallel wall surface units larger than a preset value and the maximum number as a target coordinate plane; and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
In one embodiment, before rotating the geometric object plane by a preset angle around the rotation coordinate axis in the cartesian three-dimensional rectangular coordinate system, the method further includes: determining a quantity interval to which the quantity of the parallel wall surface units corresponding to the target coordinate surface belongs; and determining the preset angle according to the preset rotation angle associated with the number interval.
In one embodiment, the creating a search box for the rotated temporary spatial point includes: determining a search diameter; on a plane parallel to the target coordinate plane, establishing a square by taking the temporary space point obtained by rotation as the center and the search diameter as the side length; taking the square as a search box of the temporary space point; two adjacent edges of the search box are respectively parallel to two coordinate axes in the target coordinate plane.
A turbulent wall distance calculation apparatus, the apparatus comprising:
the parallel checking module is used for establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane; determining the number of wall surface units of a coordinate surface parallel to the Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
the rotating module is used for determining a rotating coordinate axis according to the number of the parallel wall surface units corresponding to each coordinate surface; in the Cartesian three-dimensional rectangular coordinate system, the geometric object surface is rotated by a preset angle around the rotation coordinate axis, and the data of the wall surface unit of the rotated geometric object surface is stored in an ADT data structure; in the Cartesian three-dimensional rectangular coordinate system, rotating the space points by the preset angle around the rotation coordinate axis, and establishing a search box for the temporary space points obtained by rotation;
a turbulent wall distance calculation module for searching the ADT data structure for the wall units falling within the range of the search box and recording as target units; and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
In one embodiment, the parallelism checking module is further configured to traverse each wall element on the geometry plane; judging whether the normal vector of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system or not; and if so, adding one to the number of the parallel wall surface units of the coordinate surface perpendicular to the coordinate axis to obtain the number of the parallel wall surface units corresponding to each coordinate surface.
In an embodiment, the rotation module is further configured to determine whether there is at least one of the number of parallel wall units corresponding to the three coordinate surfaces, and if so, determine a coordinate surface with the number of parallel wall units greater than a preset value and the maximum number as the target coordinate surface; and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane;
determining the number of wall surface units of a coordinate surface parallel to the Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
determining a rotation coordinate axis according to the number of parallel wall surface units corresponding to each coordinate surface;
in the Cartesian three-dimensional rectangular coordinate system, the geometric object surface is rotated by a preset angle around the rotation coordinate axis, and the data of the wall surface unit of the rotated geometric object surface is stored in an ADT data structure;
in the Cartesian three-dimensional rectangular coordinate system, rotating the space points by the preset angle around the rotation coordinate axis, and establishing a search box for the temporary space points obtained by rotation;
searching the ADT data structure for the wall surface unit falling in the range of the search box, and recording the wall surface unit as a target unit;
and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
A computer-readable storage medium, on which a computer program is stored which, when executed by a processor, carries out the steps of:
establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane;
determining the number of wall surface units of a coordinate surface parallel to the Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
determining a rotation coordinate axis according to the number of parallel wall surface units corresponding to each coordinate surface;
in the Cartesian three-dimensional rectangular coordinate system, the geometric object surface is rotated by a preset angle around the rotation coordinate axis, and the data of the wall surface unit of the rotated geometric object surface is stored in an ADT data structure;
in the Cartesian three-dimensional rectangular coordinate system, rotating the space points by the preset angle around the rotation coordinate axis, and establishing a search box for the temporary space points obtained by rotation;
searching the ADT data structure for the wall surface unit falling in the range of the search box, and recording the wall surface unit as a target unit;
and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
The method, the device, the computer equipment and the storage medium for calculating the turbulent wall surface distance are suitable for efficient turbulent wall surface distance calculation parallel to a Cartesian coordinate plane, when the wall surface distance is calculated, if a geometric object surface is parallel to the Cartesian coordinate plane, the wall surface unit and the space point are rotated, so that the rotated object surface is not parallel to the Cartesian coordinate plane any more, then the wall surface distance is calculated, the problem that the wall surface unit searched by an original ADT method is not available or too many, the problem that the original ADT method does not work is solved, and the wall surface distance calculation efficiency is greatly improved.
Drawings
FIG. 1 is a diagram of an embodiment of an application of the method for calculating the distance to a turbulent wall.
FIG. 2 is a flow chart illustrating a method for calculating a turbulent wall distance in one embodiment.
FIG. 3 is a schematic diagram of a geometry plane parallel to a Cartesian coordinate plane in one embodiment.
FIG. 4 is a diagram illustrating an angle of rotation of a geometric object plane about a rotational axis in one embodiment.
FIG. 5 is a diagram illustrating an angle of rotation of a spatial point about a rotational axis in one embodiment.
FIG. 6 is a schematic representation of the relative positions of the rotated geometry plane and the temporary spatial points in one embodiment.
Fig. 7 is a block diagram showing a structure of a turbulent wall distance calculating device in one embodiment.
FIG. 8A is a diagram illustrating an internal structure of a computer device in one embodiment.
Fig. 8B is an internal structural diagram of a computer device in another 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.
The turbulent wall surface distance calculation method provided by the application can be applied to the application environment shown in fig. 1. Wherein the terminal 102 communicates with the server 104 via a network. The turbulent wall surface distance calculation method provided by the application is executed by a computer device, specifically, the turbulent wall surface distance calculation method can be executed by the terminal 102, the server 104, or the terminal 102 and the server 104 cooperatively. The terminal 102 may be, but not limited to, various personal computers, notebook computers, smart phones, tablet computers, and portable wearable devices, and the server 104 may be implemented by an independent server or a server cluster formed by a plurality of servers.
In one embodiment, as shown in fig. 2, a turbulent wall distance calculation method is provided, which is described by taking the method as an example applied to the terminal in fig. 1, and includes the following steps:
step 202, a cartesian three-dimensional rectangular coordinate system of the geometric object plane is established.
Computational mesh refers to the discretization of a continuous flow field into a finite number of discrete elements (triangles, tetrahedrons, etc.), these elements and the geometric topological relationships between them, called computational mesh. The geometry plane refers to the surface of a solid contained within the computational domain of the CFD simulation. Wall cells refer to discrete cells on a geometric object plane in a computational grid. The cartesian coordinate system is a generic term referring to a rectangular coordinate system and a diagonal coordinate system. Two coordinate axes intersecting the origin constitute a planar affine coordinate system. If the measurement units on the two coordinate axes are equal, the affine coordinate system is called as a cartesian coordinate system. A cartesian coordinate system with two mutually perpendicular coordinate axes is called a cartesian rectangular coordinate system, otherwise called a cartesian oblique coordinate system. The Cartesian coordinate systems adopted by the application are all Cartesian three-dimensional rectangular coordinate systems.
And step 204, determining the number of wall surface units of a coordinate plane parallel to the Cartesian three-dimensional rectangular coordinate system in the geometric object plane.
The coordinate axes of the Cartesian three-dimensional rectangular coordinate system comprise an X axis, a Y axis and a Z axis which are perpendicular to each other in pairs. The coordinate plane of the cartesian three-dimensional rectangular coordinate system is a plane formed by two coordinate axes in the cartesian three-dimensional rectangular coordinate system, and specifically comprises three XOY planes, YOZ planes and XOZ planes.
In one embodiment, determining the number of wall elements in the geometry plane parallel to a coordinate plane of a cartesian three-dimensional rectangular coordinate system comprises: traversing each wall surface unit on the geometric object surface; judging whether the normal vector of the surface of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system or not; and if so, adding one to the number of the parallel wall units of the coordinate plane perpendicular to the coordinate axis to obtain the number of the parallel wall units corresponding to each coordinate plane.
Specifically, the computer device traverses each wall surface unit on the geometric object surface, and judges whether the surface normal vector of the wall surface unit is parallel to the cartesian coordinate axis, that is, whether the wall surface unit is parallel to the cartesian coordinate plane. For example, if the normal vector of the wall element is parallel to the X-axis, the wall element is parallel to the YOZ plane of the cartesian coordinate plane. As shown in fig. 3, if the wall element is parallel to the cartesian coordinate plane, the computer device saves the total number of wall elements parallel to the YOZ plane of the cartesian coordinate plane, recorded as nFaceX. Similarly, the computer device records the total number of wall elements nFaceY, nFaceZ parallel to the XOZ, XOY plane of the cartesian coordinate plane. In fig. 3, each dot represents a wall element.
When nFaceX, nFaceY, nFaceZ are less than or equal to the preset values, the computer device assumes that the geometric object plane is not parallel to the Cartesian coordinate plane. If there is a condition of greater than a preset value for nFaceX, nFaceY, nFaceZ, the computer device assumes that the geometric object plane is parallel to the Cartesian coordinate plane.
And step 206, determining a rotation coordinate axis according to the number of the parallel wall surface units corresponding to each coordinate surface.
In one embodiment, determining the rotation coordinate axis according to the number of the parallel wall units corresponding to each coordinate plane includes: judging whether the number of parallel wall surface units corresponding to at least one coordinate surface is larger than a preset value in the number of parallel wall surface units corresponding to the three coordinate surfaces; if so, identifying a coordinate plane with the number of parallel wall surface units larger than a preset value and the maximum number as a target coordinate plane; and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
Specifically, if nFaceX > nFaceY and nFaceX > nFaceZ, the computer device determines that the geometric object plane is parallel to the YOZ plane of the Cartesian coordinate plane, i.e., the target coordinate plane is the YOZ plane and the rotation coordinate axis is the X axis; if nFaceY > nFaceX and nFaceY > nFaceZ, the computer device determines that the geometric object plane is parallel to the XOZ plane of the Cartesian coordinate plane, namely the target coordinate plane is the XOZ plane and the rotation coordinate axis is the Y axis; if nFaceZ > nFaceX and nFaceZ > nFaceY, the computer device assumes that the geometry plane is parallel to the XOY plane of the Cartesian coordinate plane, i.e., the target coordinate plane is the XOY plane and the rotation coordinate axis is the Z axis.
And step 208, in a Cartesian three-dimensional rectangular coordinate system, rotating the geometric object plane by a preset angle around a rotation coordinate axis, and storing the data of the wall surface unit of the rotated geometric object plane into an ADT data structure.
Specifically, as shown in fig. 4, the coordinate values of the wall surface units are rotated around the rotation coordinate axis by an angle θ, and the coordinate values of the wall surface units on the rotated geometric object plane are stored in the ADT.
And step 210, rotating the space points by a preset angle around a rotation coordinate axis in a Cartesian three-dimensional rectangular coordinate system, and establishing a search box for the temporary space points obtained by rotation.
The wall surface distance refers to the closest distance from any point in space to the geometric object surface.
Specifically, when calculating the wall surface distance of any point in the space, the computer device rotates the coordinates of the space point to generate a temporary space point for calculating the wall surface distance. As shown in fig. 5, the rotation operation of the spatial point is kept in agreement with the rotation operation of the wall surface unit.
Further, determining a proper search radius, and constructing a search box of the temporary space point. The search box may be a circular area centered on the temporary space point and having a preset search radius as a radius.
In one embodiment, establishing a search box for the rotated temporary spatial points comprises: determining a search diameter; on a plane parallel to a target coordinate plane, establishing a square by taking a temporary space point obtained by rotation as a center and a search diameter as a side length; taking the square as a search box of the temporary space point; two adjacent edges of the search box are respectively parallel to two coordinate axes in the target coordinate plane.
In one embodiment, before rotating the geometric object plane by a preset angle around the rotation coordinate axis in the cartesian three-dimensional rectangular coordinate system, the method further includes: determining a quantity interval to which the quantity of the parallel wall surface units corresponding to the target coordinate surface belongs; and determining the preset angle according to the preset rotation angle associated with the number interval. In other words, the preset angle may be a preset fixed angle (for example, 45 ° clockwise), or may be an angle dynamically determined based on a mapping relationship between preset intervals of the number of parallel wall units and the rotation angle, or in other manners.
At step 212, the ADT data structure is searched for wall cells that fall within the search box, denoted as target cells.
Specifically, the computer device searches the ADT data structure for the wall cells contained by the search box based on the constructed search box. Fig. 6 shows the relative positions of the geometric object plane and the temporary spatial point after the geometric object plane and the spatial point are synchronously rotated for a certain angle. In fig. 6, each dot represents a wall cell, the dotted square represents a search box for a temporary space point, the light dots represent wall cells that fall into the search box, and the dark dots represent wall cells that do not fall into the search box.
And 214, calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
Specifically, the computer device calculates the straight-line distances of the temporary spatial points to the center points of the wall cells. And taking the minimum value of all the linear distances to obtain the wall surface distance from the temporary space point to the rotated geometric object surface. In the embodiment, the wall surface units are rotated, and only the rotated wall surface units are stored in the ADT, so that the original calculation grid cannot be changed; and rotating the space points, and only performing wall surface distance calculation by using the rotated temporary space points without changing the original calculation grid. And performing the same rotation operation on the wall surface unit and the space point, so that the relative position of the temporary space point to the wall surface unit in the ADT is consistent with the relative position of the real space point to the real wall surface unit, and the correctness of the wall surface distance calculation result is ensured. Therefore, the calculated wall distance is also the wall distance from the real space point to the real geometric object plane.
According to the high-efficiency turbulence wall surface distance calculation method suitable for the parallel Cartesian coordinate surface, when the wall surface distance is calculated, if the geometric object surface is parallel to the Cartesian coordinate surface, the wall surface unit and the space point are rotated, so that the rotated object surface is not parallel to the Cartesian coordinate surface, then the wall surface distance is calculated, the problems that the wall surface unit searched by an original ADT method is not available or too many, the original ADT method does not work are solved, and the wall surface distance calculation efficiency is greatly improved.
It should be understood that, although the steps in the flowchart of fig. 2 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. 2 may include multiple steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed in turn or alternately with other steps or at least a portion of the other steps or stages.
In one embodiment, as shown in fig. 7, there is provided a turbulent wall distance calculation apparatus 700, including: a parallelism checking module 702, a rotation module 704, and a turbulence wall distance calculation module 706, wherein:
a parallel checking module 702, configured to establish a cartesian three-dimensional rectangular coordinate system of a geometric object plane; determining the number of wall surface units of a coordinate surface parallel to a Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
a rotation module 704, configured to determine a rotation coordinate axis according to the number of parallel wall units corresponding to each coordinate plane; in a Cartesian three-dimensional rectangular coordinate system, rotating a geometric object plane by a preset angle around a rotation coordinate axis, and storing data of a wall surface unit of the rotated geometric object plane into an ADT data structure; in a Cartesian three-dimensional rectangular coordinate system, rotating a space point by a preset angle around a rotating coordinate axis, and establishing a search box for the temporary space point obtained by rotation;
a turbulent wall distance calculation module 706 for searching the ADT data structure for wall cells falling within the search box range, denoted as target cells; and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
In one embodiment, the parallelism checking module is further configured to traverse each wall element on the geometry plane; judging whether the normal vector of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system; and if so, adding one to the number of the parallel wall surface units of the coordinate surface vertical to the coordinate axis to obtain the number of the parallel wall surface units corresponding to each coordinate surface.
In one embodiment, the rotation module is further configured to determine whether there is at least one of the number of parallel wall units corresponding to the coordinate surfaces that is greater than a preset value, and if so, determine a coordinate surface with the largest number of parallel wall units that is greater than the preset value as the target coordinate surface; and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
In one embodiment, the rotation module is further configured to determine a number interval to which the number of parallel wall surface units corresponding to the target coordinate surface belongs; and determining the preset angle according to the preset rotation angle associated with the number interval.
In one embodiment, the turbulent wall distance calculation module is further configured to determine a search diameter; on a plane parallel to a target coordinate plane, establishing a square by taking a temporary space point obtained by rotation as a center and a search diameter as a side length; taking the square as a search box of the temporary space point; two adjacent edges of the search box are respectively parallel to two coordinate axes in the target coordinate plane.
According to the high-efficiency turbulence wall surface distance calculation method suitable for the parallel Cartesian coordinate surface, when the wall surface distance is calculated, if the geometric object surface is parallel to the Cartesian coordinate surface, the wall surface unit and the space point are rotated, so that the rotated object surface is not parallel to the Cartesian coordinate surface, then the wall surface distance is calculated, the problems that the wall surface unit searched by an original ADT method is not available or too many, the original ADT method does not work are solved, and the wall surface distance calculation efficiency is greatly improved.
For specific definition of the turbulent wall distance calculation device, reference may be made to the above definition of the turbulent wall distance calculation method, which is not described herein again. The modules in the turbulent wall distance calculation device 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 server, and its internal structure diagram may be as shown in fig. 8A. The computer device includes a processor, a memory, and a network interface 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, a computer program, and a database. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The database of the computer device is used to store ADT structure data. 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 turbulent wall distance calculation method.
In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as shown in fig. 8B. The computer device includes a processor, a memory, a communication 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 operating system and the computer program to run on the non-volatile storage medium. The communication interface of the computer device is used for carrying out wired or wireless communication with an external terminal, and the wireless communication can be realized through WIFI, an operator network, NFC (near field communication) or other technologies. The computer program is executed by a processor to implement a turbulent wall distance calculation method. 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 the shell of the computer equipment, an external keyboard, a touch pad or a mouse and the like.
Those skilled in the art will appreciate that the configurations shown in fig. 8A and 8B are merely block diagrams of some configurations relevant to the present disclosure, and do not constitute a limitation on the computing devices to which the present disclosure may be applied, and that a particular computing device may include more or less components than those shown, or combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory and a processor, the memory having a computer program stored therein, the processor implementing the following steps when executing the computer program:
establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane;
determining the number of wall surface units of a coordinate surface parallel to a Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
determining a rotation coordinate axis according to the number of parallel wall surface units corresponding to each coordinate surface;
in a Cartesian three-dimensional rectangular coordinate system, rotating a geometric object plane by a preset angle around a rotation coordinate axis, and storing data of a wall surface unit of the rotated geometric object plane into an ADT data structure;
in a Cartesian three-dimensional rectangular coordinate system, rotating a space point by a preset angle around a rotating coordinate axis, and establishing a search box for the temporary space point obtained by rotation;
searching wall surface units falling in the range of the search box in the ADT data structure, and recording the wall surface units as target units;
and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
In one embodiment, the processor, when executing the computer program, further performs the steps of: determining the number of wall elements in the geometric object plane parallel to a coordinate plane of a cartesian three-dimensional rectangular coordinate system, comprising: traversing each wall surface unit on the geometric object surface; judging whether the normal vector of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system; and if so, adding one to the number of the parallel wall surface units of the coordinate surface vertical to the coordinate axis to obtain the number of the parallel wall surface units corresponding to each coordinate surface.
In one embodiment, the processor, when executing the computer program, further performs the steps of: determining a rotation coordinate axis according to the number of the parallel wall units corresponding to each coordinate surface, wherein the method comprises the following steps: judging whether the number of parallel wall surface units corresponding to at least one coordinate surface is larger than a preset value in the number of parallel wall surface units corresponding to the three coordinate surfaces; if so, identifying a coordinate plane with the number of parallel wall surface units larger than a preset value and the maximum number as a target coordinate plane; and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
In one embodiment, the processor, when executing the computer program, further performs the steps of: in the three-dimensional rectangular coordinate system of cartesian, before rotating the geometry face around the rotation coordinate axis and predetermineeing the angle, still include: determining a quantity interval to which the quantity of the parallel wall surface units corresponding to the target coordinate surface belongs; and determining the preset angle of the preset rotation angle associated with the number interval.
In one embodiment, the processor when executing the computer program further performs the steps of: establishing a search box for the rotated temporary spatial point, comprising: determining a search diameter; on a plane parallel to a target coordinate plane, establishing a square by taking a temporary space point obtained by rotation as a center and a search diameter as a side length; taking the square as a search box of the temporary space point; two adjacent edges of the search box are respectively parallel to two coordinate axes in the target coordinate plane.
According to the high-efficiency turbulence wall surface distance calculation method suitable for the parallel Cartesian coordinate surface, when the wall surface distance is calculated, if the geometric object surface is parallel to the Cartesian coordinate surface, the wall surface unit and the space point are rotated, so that the rotated object surface is not parallel to the Cartesian coordinate surface, then the wall surface distance is calculated, the problems that the wall surface unit searched by an original ADT method is not available or too many, the original ADT method does not work are solved, and the wall surface distance calculation efficiency is greatly improved.
In one embodiment, a computer-readable storage medium is provided, having a computer program stored thereon, which when executed by a processor, performs the steps of:
establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane;
determining the number of wall surface units of a coordinate surface parallel to a Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
determining a rotation coordinate axis according to the number of parallel wall surface units corresponding to each coordinate surface;
in a Cartesian three-dimensional rectangular coordinate system, rotating a geometric object plane by a preset angle around a rotation coordinate axis, and storing data of a wall surface unit of the rotated geometric object plane into an ADT data structure;
in a Cartesian three-dimensional rectangular coordinate system, rotating a space point by a preset angle around a rotating coordinate axis, and establishing a search box for the temporary space point obtained by rotation;
searching wall surface units falling in the range of the search box in the ADT data structure, and recording the wall surface units as target units;
and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
In one embodiment, the computer program when executed by the processor further performs the steps of: determining the number of wall elements in the geometric object plane parallel to a coordinate plane of a cartesian three-dimensional rectangular coordinate system, comprising: traversing each wall surface unit on the geometric object surface; judging whether the normal vector of the surface of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system or not; and if so, adding one to the number of the parallel wall surface units of the coordinate surface vertical to the coordinate axis to obtain the number of the parallel wall surface units corresponding to each coordinate surface.
In one embodiment, the computer program when executed by the processor further performs the steps of: determining a rotation coordinate axis according to the number of the parallel wall units corresponding to each coordinate surface, wherein the method comprises the following steps: judging whether the number of parallel wall surface units corresponding to at least one coordinate surface is larger than a preset value in the number of parallel wall surface units corresponding to the three coordinate surfaces; if so, identifying a coordinate plane with the number of parallel wall surface units larger than a preset value and the maximum number as a target coordinate plane; and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
In one embodiment, the computer program when executed by the processor further performs the steps of: in the three-dimensional rectangular coordinate system of cartesian, before rotating the geometry face around the rotation coordinate axis and predetermineeing the angle, still include: determining a quantity interval to which the quantity of the parallel wall surface units corresponding to the target coordinate surface belongs; and determining the preset angle according to the preset rotation angle associated with the number interval.
In one embodiment, the computer program when executed by the processor further performs the steps of: establishing a search box for the rotated temporary spatial points, comprising: determining a search diameter; on a plane parallel to a target coordinate plane, taking a temporary space point obtained by rotation as a center and a search diameter as a side length, establishing a square, and taking the square as a search box of the temporary space point; two adjacent edges of the search box are respectively parallel to two coordinate axes in the target coordinate plane.
According to the high-efficiency turbulence wall surface distance calculating method suitable for the parallel Cartesian coordinate surface, when the wall surface distance is calculated, if the geometric object surface is parallel to the Cartesian coordinate surface, the wall surface unit and the space point are rotated, so that the rotated object surface is not parallel to the Cartesian coordinate surface, then the wall surface distance is calculated, the problems that the wall surface unit searched by the original ADT method is not available or too many, the original ADT method does not work are solved, and the wall surface distance calculating efficiency is greatly improved.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware related to instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can 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 can include at least one of non-volatile and volatile memory. Non-volatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical storage, or the like. Volatile Memory can include Random Access Memory (RAM) or external cache Memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.
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 examples 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 method for calculating a distance to a turbulent wall surface, the method comprising:
establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane;
determining the number of wall surface units of a coordinate surface parallel to the Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
determining a rotation coordinate axis according to the number of the parallel wall surface units corresponding to each coordinate surface;
in the Cartesian three-dimensional rectangular coordinate system, the geometric object surface is rotated by a preset angle around the rotation coordinate axis, and the data of the wall surface unit of the rotated geometric object surface is stored in an ADT data structure;
in the Cartesian three-dimensional rectangular coordinate system, rotating the space points by the preset angle around the rotation coordinate axis, and establishing a search box for the temporary space points obtained by rotation;
searching the ADT data structure for the wall surface unit falling in the range of the search box, and recording the wall surface unit as a target unit;
and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the turbulent wall surface distance.
2. The method of claim 1, wherein determining the number of wall elements in the geometric object plane that are parallel to the coordinate plane of the cartesian three-dimensional rectangular coordinate system comprises:
traversing each wall surface unit on the geometric object surface;
judging whether the normal vector of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system or not;
and if so, adding one to the number of the parallel wall surface units of the coordinate surface perpendicular to the coordinate axis to obtain the number of the parallel wall surface units corresponding to each coordinate surface.
3. The method of claim 2, wherein determining the rotation axis according to the number of parallel wall elements corresponding to each coordinate plane comprises:
judging whether the number of parallel wall surface units corresponding to at least one coordinate surface is larger than a preset value in the number of parallel wall surface units corresponding to the three coordinate surfaces;
if so, identifying a coordinate plane with the number of parallel wall surface units larger than a preset value and the maximum number as a target coordinate plane;
and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
4. The method of claim 3, wherein prior to rotating the geometric object plane by a preset angle around the rotation coordinate axis in the Cartesian three-dimensional rectangular coordinate system, further comprises:
determining a quantity interval to which the quantity of the parallel wall surface units corresponding to the target coordinate surface belongs;
and determining the preset angle according to the preset rotation angle associated with the number interval.
5. The method of claim 3, wherein the building a search box for the rotated temporary spatial points comprises:
determining a search diameter;
on a plane parallel to the target coordinate plane, establishing a square by taking the temporary space point obtained by rotation as the center and the search diameter as the side length;
taking the square as a search box of the temporary space point; two adjacent edges of the search box are respectively parallel to two coordinate axes in the target coordinate plane.
6. A turbulent wall distance calculation apparatus, comprising:
the parallel checking module is used for establishing a Cartesian three-dimensional rectangular coordinate system of a geometric object plane; determining the number of wall surface units of a coordinate surface parallel to the Cartesian three-dimensional rectangular coordinate system in the geometric object surface;
the rotating module is used for determining a rotating coordinate axis according to the number of the parallel wall surface units corresponding to each coordinate surface; in the Cartesian three-dimensional rectangular coordinate system, the geometric object plane is rotated by a preset angle around the rotation coordinate axis, and data of a wall surface unit of the rotated geometric object plane are stored in an ADT data structure; in the Cartesian three-dimensional rectangular coordinate system, rotating the space points by the preset angle around the rotation coordinate axis, and establishing a search box for the temporary space points obtained by rotation;
a turbulent wall distance calculation module for searching the ADT data structure for the wall units falling within the range of the search box and recording as target units; and calculating the minimum distance between the rotated space point and the central point of each target unit to obtain the distance of the turbulent wall surface.
7. The apparatus of claim 6, wherein the parallelism inspection module is further configured to traverse each wall element on the geometry plane; judging whether the normal vector of the wall surface unit is parallel to the coordinate axis of the Cartesian three-dimensional rectangular coordinate system or not; and if so, adding one to the number of the parallel wall surface units of the coordinate surface perpendicular to the coordinate axis to obtain the number of the parallel wall surface units corresponding to each coordinate surface.
8. The apparatus of claim 7, wherein the rotation module is further configured to determine whether there is at least one of the number of parallel wall units corresponding to the three coordinate planes, and if so, determine a coordinate plane having the largest number of parallel wall units and the number of parallel wall units larger than a predetermined value as the target coordinate plane; and determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor when executing the computer program performs the steps of the method according to any of claims 1 to 5.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210650302.8A CN114722745B (en) | 2022-06-10 | 2022-06-10 | Turbulent flow wall surface distance calculation method and device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210650302.8A CN114722745B (en) | 2022-06-10 | 2022-06-10 | Turbulent flow wall surface distance calculation method and device, computer equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114722745A true CN114722745A (en) | 2022-07-08 |
| CN114722745B CN114722745B (en) | 2022-08-26 |
Family
ID=82232598
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210650302.8A Active CN114722745B (en) | 2022-06-10 | 2022-06-10 | Turbulent flow wall surface distance calculation method and device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114722745B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115357849A (en) * | 2022-10-24 | 2022-11-18 | 中国空气动力研究与发展中心计算空气动力研究所 | Method and device for calculating wall surface distance under Cartesian grid |
| CN117217130A (en) * | 2023-11-07 | 2023-12-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Layer-by-layer advancing wall distance determining method, device and medium |
| CN117390775A (en) * | 2023-12-13 | 2024-01-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Viscous Cartesian grid generation method, viscous Cartesian grid generation equipment, viscous Cartesian grid medium and viscous Cartesian grid generation system |
| CN117521563A (en) * | 2024-01-08 | 2024-02-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Pneumatic data processing method based on impeller mechanical turbulence wall distance calculation |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030002049A1 (en) * | 2001-06-20 | 2003-01-02 | Michael Kuchel | Scanning interferometer for aspheric surfaces and wavefronts |
| CN101882177A (en) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | Aeroelastic stability fluid-structure interaction prediction method of turbo-machine changed interblade phase angles |
| CN102682144A (en) * | 2011-11-30 | 2012-09-19 | 天津空中代码工程应用软件开发有限公司 | Flight icing numerical value simulation method of helicopter rotor wing |
| CN103737426A (en) * | 2013-12-24 | 2014-04-23 | 西安交通大学 | Numerical control machine tool rotating shaft geometric error three-wire measurement method |
| CN103900167A (en) * | 2014-03-21 | 2014-07-02 | 美的集团股份有限公司 | Air outlet protection hood of fan, fan and air conditioner outdoor machine |
| RU157012U1 (en) * | 2015-06-23 | 2015-11-20 | федеральное государственное бюджетное образовательное учреждение высшего образования "Национальный исследовательский университет "МЭИ" (ФГБОУ ВО "НИУ "МЭИ") | LONGITUDINAL PROBE FOR MEASURING LOCAL CHARACTERISTICS OF TURBULENT FLOWS IN PIPES |
| CN107391950A (en) * | 2017-08-03 | 2017-11-24 | 中国空气动力研究与发展中心计算空气动力研究所 | A kind of turbulent flow wall is apart from quick calculation method |
| CN109948002A (en) * | 2019-03-14 | 2019-06-28 | 空气动力学国家重点实验室 | The nearest wall surface of unstrctured grid based on balance KD tree is apart from method for solving |
| CN110476558A (en) * | 2019-08-26 | 2019-11-22 | 南京理工大学 | Transplanting hole tray seedling machine based on weathering steel framework |
| CN111501658A (en) * | 2020-04-20 | 2020-08-07 | 水利部交通运输部国家能源局南京水利科学研究院 | Hydrodynamic lifting device based on vertical shaft plug flow and performance test method thereof |
| CN112100937A (en) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | Efficient overlapped grid parallel assembly method based on wall surface distance |
-
2022
- 2022-06-10 CN CN202210650302.8A patent/CN114722745B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030002049A1 (en) * | 2001-06-20 | 2003-01-02 | Michael Kuchel | Scanning interferometer for aspheric surfaces and wavefronts |
| CN101882177A (en) * | 2010-06-18 | 2010-11-10 | 北京航空航天大学 | Aeroelastic stability fluid-structure interaction prediction method of turbo-machine changed interblade phase angles |
| CN102682144A (en) * | 2011-11-30 | 2012-09-19 | 天津空中代码工程应用软件开发有限公司 | Flight icing numerical value simulation method of helicopter rotor wing |
| CN103737426A (en) * | 2013-12-24 | 2014-04-23 | 西安交通大学 | Numerical control machine tool rotating shaft geometric error three-wire measurement method |
| CN103900167A (en) * | 2014-03-21 | 2014-07-02 | 美的集团股份有限公司 | Air outlet protection hood of fan, fan and air conditioner outdoor machine |
| RU157012U1 (en) * | 2015-06-23 | 2015-11-20 | федеральное государственное бюджетное образовательное учреждение высшего образования "Национальный исследовательский университет "МЭИ" (ФГБОУ ВО "НИУ "МЭИ") | LONGITUDINAL PROBE FOR MEASURING LOCAL CHARACTERISTICS OF TURBULENT FLOWS IN PIPES |
| CN107391950A (en) * | 2017-08-03 | 2017-11-24 | 中国空气动力研究与发展中心计算空气动力研究所 | A kind of turbulent flow wall is apart from quick calculation method |
| CN109948002A (en) * | 2019-03-14 | 2019-06-28 | 空气动力学国家重点实验室 | The nearest wall surface of unstrctured grid based on balance KD tree is apart from method for solving |
| CN110476558A (en) * | 2019-08-26 | 2019-11-22 | 南京理工大学 | Transplanting hole tray seedling machine based on weathering steel framework |
| CN111501658A (en) * | 2020-04-20 | 2020-08-07 | 水利部交通运输部国家能源局南京水利科学研究院 | Hydrodynamic lifting device based on vertical shaft plug flow and performance test method thereof |
| CN112100937A (en) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | Efficient overlapped grid parallel assembly method based on wall surface distance |
Non-Patent Citations (3)
| Title |
|---|
| FRANCESCO CAPIZZANO: "Turbulent Wall Model for Immersed Boundary Methods", 《AIAA JOURNAL》 * |
| 沈志伟: "基于混合笛卡尔网格方法的非定常流动问题研究", 《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》 * |
| 赵钟: "湍流模拟壁面距离MPI/OpenMP混合并行计算方法", 《空气动力学学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115357849A (en) * | 2022-10-24 | 2022-11-18 | 中国空气动力研究与发展中心计算空气动力研究所 | Method and device for calculating wall surface distance under Cartesian grid |
| CN117217130A (en) * | 2023-11-07 | 2023-12-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Layer-by-layer advancing wall distance determining method, device and medium |
| CN117217130B (en) * | 2023-11-07 | 2024-03-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Layer-by-layer advancing wall distance determining method, device and medium |
| CN117390775A (en) * | 2023-12-13 | 2024-01-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Viscous Cartesian grid generation method, viscous Cartesian grid generation equipment, viscous Cartesian grid medium and viscous Cartesian grid generation system |
| CN117390775B (en) * | 2023-12-13 | 2024-02-27 | 中国空气动力研究与发展中心计算空气动力研究所 | Viscous Cartesian grid generation method, viscous Cartesian grid generation equipment, viscous Cartesian grid medium and viscous Cartesian grid generation system |
| CN117521563A (en) * | 2024-01-08 | 2024-02-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Pneumatic data processing method based on impeller mechanical turbulence wall distance calculation |
| CN117521563B (en) * | 2024-01-08 | 2024-03-15 | 中国空气动力研究与发展中心计算空气动力研究所 | Pneumatic data processing method based on impeller mechanical turbulence wall distance calculation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114722745B (en) | 2022-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114722745B (en) | Turbulent flow wall surface distance calculation method and device, computer equipment and storage medium | |
| Neiva et al. | A scalable parallel finite element framework for growing geometries. Application to metal additive manufacturing | |
| Sieger et al. | On shape deformation techniques for simulation-based design optimization | |
| Brix et al. | Parallelisation of multiscale-based grid adaptation using space-filling curves | |
| Jafari et al. | Developing an algorithm for reconstruction blocky systems in discontinuous media: three‐dimensional analysis | |
| CN113987666B (en) | BIM (building information modeling) model examination method, device, equipment and storage medium | |
| Lee et al. | Topology-oriented incremental algorithm for the robust construction of the Voronoi diagrams of disks | |
| CN110008387B (en) | Flow field visualization implementation method and device and electronic equipment | |
| Fang et al. | A contact detection algorithm for multi‐sphere particles by means of two‐level‐grid‐searching in DEM simulations | |
| Xing et al. | Efficient collision detection based on AABB trees and sort algorithm | |
| Ji et al. | A robust and efficient hybrid cut-cell/ghost-cell method with adaptive mesh refinement for moving boundaries on irregular domains | |
| CN112528428A (en) | Method and device for displaying physical parameters of engineering structure and computer equipment | |
| Vidović et al. | Convex combinations for diffusion schemes | |
| Nie et al. | Development of an object-oriented finite element program with adaptive mesh refinement for multi-physics applications | |
| Hu et al. | Voronoi cells of non-general position spheres using the GPU | |
| TWI406189B (en) | Method for constructing triangular grids of point clouds | |
| CN115859749A (en) | Constraint building method and device of three-dimensional model, electronic equipment and storage medium | |
| Tafuri et al. | Automatic collision free path planning in hybrid triangle and point models: a case study | |
| Novikov et al. | Layer-by-layer partitioning of finite element meshes for multicore architectures | |
| CN112989193A (en) | Data processing method and device, electronic equipment and computer storage medium | |
| Cui et al. | Point-in-polyhedra test with direct handling of degeneracies | |
| de Magalhães et al. | Exact fast parallel intersection of large 3-D triangular meshes | |
| Soukal et al. | Walking algorithms for point location in TIN models | |
| Zapata et al. | A GPU parallel finite volume method for a 3D Poisson equation on arbitrary geometries | |
| Gao et al. | Flip to regular triangulation and convex hull |
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 |