CN114722745B - 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
- CN114722745B CN114722745B CN202210650302.8A CN202210650302A CN114722745B CN 114722745 B CN114722745 B CN 114722745B CN 202210650302 A CN202210650302 A CN 202210650302A CN 114722745 B CN114722745 B CN 114722745B
- Authority
- CN
- China
- Prior art keywords
- wall surface
- coordinate
- cartesian
- parallel
- 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.)
- Active
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, which 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 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. By adopting the method, the calculation efficiency of the turbulent wall surface distance can be improved.
Description
Technical Field
The application relates to the technical field of computational fluid mechanics, in particular to a method and a device for calculating a distance between turbulent wall surfaces, computer equipment 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 are used for turbulent flow numerical simulation: 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, it is necessary to provide a turbulent wall distance method, apparatus, computer device and storage medium capable of improving efficiency of calculating a wall distance in a case where a geometric object plane is parallel to a 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 rotating coordinate axis, and establishing a search box for the temporary space points obtained through 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 the 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 of 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 rotating coordinate axis, and establishing a search box for the temporary space points obtained through rotation;
a turbulent wall distance calculation module for searching the ADT data structure for the wall units falling in 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 a surface normal vector of a wall surface unit is parallel to a 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 largest number as a 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 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 rotating coordinate axis, and establishing a search box for the temporary space points obtained through 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 representation of a geometry plane parallel to a Cartesian coordinate plane in one embodiment.
FIG. 4 is a schematic illustration of an angle of rotation of a geometric surface 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 an internal structural diagram 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 coordinate axes perpendicular to each other is called a cartesian rectangular coordinate system, otherwise it is 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 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; 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.
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 the nFaceX, nFaceY, nFaceZ conditions are greater than the preset values, the computer device determines 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, determining 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 208, rotating the geometric object plane by a preset angle around the rotation coordinate axis in a Cartesian three-dimensional rectangular coordinate system, 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 axis by an angle θ, and the coordinate values of the wall surface units on the rotated geometric object surface are stored in the ADT.
And step 210, in a Cartesian three-dimensional rectangular coordinate system, rotating the space points by a preset angle around a rotation coordinate axis, and establishing a search box for the temporary space points obtained through rotation.
The wall distance refers to the closest distance from any point in space to the geometric object plane.
Specifically, when calculating the wall surface distance at any point in 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 coincident 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 with a temporary spatial point as a center and 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 of temporary spatial points, the light dots represent wall cells falling into the search box, and the dark dots represent wall cells not falling 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 completely same rotation operation is carried out on the wall surface unit and the space point, so that the relative position from the temporary space point to the wall surface unit in the ADT is consistent with the relative position from 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 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 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 check 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 a geometric object plane;
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, a geometric object plane is rotated by a preset angle around a rotation coordinate axis, and data of a wall surface unit of the rotated geometric object plane are stored in 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 turbulent wall surface distance.
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 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 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 of the preset rotation angle associated with the number interval.
In one embodiment, the turbulence wall distance calculation module is further configured to determine a search diameter; on a plane parallel to a target coordinate surface, 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 a 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.
The specific definition of the turbulence wall distance calculation device can be referred to the definition of the turbulence wall distance calculation method in the foregoing, and will not be described in detail here. The modules in the turbulence wall distance calculating 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 operation of an operating system and computer programs in the non-volatile storage medium. The communication interface of the computer device is used for communicating with an external terminal in a wired or wireless manner, and the wireless manner 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 portions of 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 fewer components than shown, or may 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 turbulent wall surface distance.
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 units of the coordinate plane perpendicular to the coordinate axis to obtain the number of the parallel wall units corresponding to each coordinate plane.
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 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 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 point, 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 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 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 Memory, 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 application shall be subject to the appended claims.
Claims (8)
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;
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;
determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis;
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.
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 1, wherein before rotating the geometric object plane by a preset angle around the rotation coordinate axis in the cartesian three-dimensional rectangular coordinate system, further comprising:
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.
4. The method of claim 1, 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.
5. 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 judging whether the number of the parallel wall surface units corresponding to at least one coordinate surface is larger than a preset value or not in the number of the parallel wall surface units corresponding to the three coordinate surfaces, and determining one coordinate surface with the number of the parallel wall surface units larger than the preset value and the maximum number as a target coordinate surface if the number of the parallel wall surface units corresponding to at least one coordinate surface is larger than the preset value; determining a coordinate axis perpendicular to the target coordinate plane in the Cartesian three-dimensional rectangular coordinate system as a rotation coordinate axis; 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 rotating coordinate axis, and establishing a search box for the temporary space points obtained through 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.
6. The apparatus of claim 5, wherein the parallelization inspection module is further configured to traverse each wall element on the geometry 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.
7. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 4 when executing the computer program.
8. 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 4.
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 CN114722745A (en) | 2022-07-08 |
| CN114722745B true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115357849B (en) * | 2022-10-24 | 2023-03-28 | 中国空气动力研究与发展中心计算空气动力研究所 | Method and device for calculating wall surface distance under Cartesian grid |
| CN117217130B (en) * | 2023-11-07 | 2024-03-12 | 中国空气动力研究与发展中心计算空气动力研究所 | Layer-by-layer advancing wall distance determining method, device and medium |
| 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 |
| CN117521563B (en) * | 2024-01-08 | 2024-03-15 | 中国空气动力研究与发展中心计算空气动力研究所 | Pneumatic data processing method based on impeller mechanical turbulence wall distance calculation |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6781700B2 (en) * | 2001-06-20 | 2004-08-24 | Kuechel Michael | Scanning interferometer for aspheric surfaces and wavefronts |
| CN101882177B (en) * | 2010-06-18 | 2012-09-26 | 北京航空航天大学 | Aeroelastic stability fluid-structure interaction prediction method of turbo-machine changed interblade phase angles |
| CN102682144B (en) * | 2011-11-30 | 2014-02-19 | 天津空中代码工程应用软件开发有限公司 | Flight icing numerical value simulation method of helicopter rotor wing |
| CN103737426B (en) * | 2013-12-24 | 2015-10-21 | 西安交通大学 | A kind of Digit Control Machine Tool rotating shaft geometric error three line mensuration |
| CN103900167B (en) * | 2014-03-21 | 2017-06-06 | 美的集团股份有限公司 | The air-out protective cover of blower fan, blower fan and air-conditioner outdoor unit |
| 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 |
| CN111501658B (en) * | 2020-04-20 | 2021-10-08 | 水利部交通运输部国家能源局南京水利科学研究院 | Hydrodynamic lifting device based on vertical shaft plug flow and performance test method thereof |
| CN112100937B (en) * | 2020-09-11 | 2021-05-28 | 南京航空航天大学 | Efficient overlapped grid parallel assembly method based on wall surface distance |
-
2022
- 2022-06-10 CN CN202210650302.8A patent/CN114722745B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114722745A (en) | 2022-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114722745B (en) | Turbulent flow wall surface distance calculation method and device, computer equipment and storage medium | |
| Wang et al. | Computational algorithms for tracking dynamic fluid–structure interfaces in embedded boundary methods | |
| Bercea et al. | A structure-exploiting numbering algorithm for finite elements on extruded meshes, and its performance evaluation in Firedrake | |
| Sieger et al. | On shape deformation techniques for simulation-based design optimization | |
| Capizzano | Automatic generation of locally refined Cartesian meshes: Data management and algorithms | |
| Brix et al. | Parallelisation of multiscale-based grid adaptation using space-filling curves | |
| CN113221371A (en) | Method and device for determining critical sliding surface of side slope and terminal equipment | |
| Ji et al. | A robust and efficient hybrid cut-cell/ghost-cell method with adaptive mesh refinement for moving boundaries on irregular domains | |
| US11417073B2 (en) | System and method for generating hierarchical level-of-detail measurements for runtime calculation and visualization | |
| CN113987666A (en) | BIM (building information modeling) model examination method, device, equipment and storage medium | |
| Vidović et al. | Convex combinations for diffusion schemes | |
| CN112528428A (en) | Method and device for displaying physical parameters of engineering structure and computer equipment | |
| Hu et al. | Voronoi cells of non-general position spheres using the GPU | |
| Tafuri et al. | Automatic collision free path planning in hybrid triangle and point models: a case study | |
| Xing et al. | Efficient collision detection based on AABB trees and sort algorithm | |
| Novikov et al. | Layer-by-layer partitioning of finite element meshes for multicore architectures | |
| Gao et al. | Flip to regular triangulation and convex hull | |
| CN112989193A (en) | Data processing method and device, electronic equipment and computer storage medium | |
| Soukal et al. | Walking algorithms for point location in TIN models | |
| Nie et al. | Development of an object-oriented finite element program with adaptive mesh refinement for multi-physics applications | |
| Zapata et al. | A GPU parallel finite volume method for a 3D Poisson equation on arbitrary geometries | |
| JP3683948B2 (en) | Finite element model processing system and method | |
| Ubach et al. | On the interpolation of normal vectors for triangle meshes | |
| Fan et al. | Schwarzschild Harmonic Path Planning and Tetrahedral Mesh Smoothing | |
| CN116028679B (en) | Data processing method, device and equipment |
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 |