A cartesian mesh generation method, apparatus, device and medium

Abstract

The application discloses a Cartesian grid generation method and device, equipment and medium, relates to the technical field of scientific calculation and numerical simulation, and comprises the following steps: determining the type of a grid generation instruction; if the type is to generate a new grid, determining geometric file information and construction parameter information, generating a Cartesian grid, determining configuration information, determining fluid marker points and grid points, performing intersection point quantity calculation on the fluid marker points, obtaining intersection point quantity, judging whether the grid points match the geometric file information, updating the construction parameter information if the grid points match the geometric file information, marking the Cartesian grid based on the updated construction parameter information and the geometric file information, and obtaining a new grid; and if the type is to adjust a historical grid, screening a target historical grid, determining grid adjustment parameters, and adjusting the target historical grid. The application can improve the Cartesian grid generation efficiency, realize fast grid adjustment, and provide support for a wind field numerical model at any time during ship operation.

Description

Technical Field

[0001] This invention relates to the fields of scientific computing and numerical simulation technology, and in particular to a method, apparatus, device and medium for generating Cartesian grids. Background Technology

[0002] Currently, wind field environmental forecasting, as a component of digital battlefield construction, has important guiding significance for ship combat and carrier-based aircraft take-off and landing. Affected by the dynamic changes in the equipment layout on the ship deck and the position and number of aircraft, the grid for dynamic wind field simulation needs to be changed frequently. When simulating wind fields on ships, the first step is to generate the grid required for calculation based on the geometric model. However, ships have complex shapes and a large number of equipment on the deck (such as carrier-based aircraft, island superstructure, measuring instruments and equipment, etc.). Therefore, the manual cost of conventional grid generation methods is high. Furthermore, the equipment layout on the deck may change at any time during ship operations (such as the number and position of carrier-based aircraft, etc.), making the cost of regenerating the grid very high and unable to quickly calculate the wind field results.

[0003] As can be seen from the above, how to improve the efficiency of Cartesian mesh generation, achieve rapid mesh adjustment without manipulating geometric files, thereby avoiding repeated mesh generation, reducing the manual and time costs of mesh generation, and providing a guarantee for the continuous generation of wind field numerical models during ship operations is a problem to be solved in this field. Summary of the Invention

[0004] In view of this, the purpose of this invention is to provide a Cartesian mesh generation method, apparatus, device, and medium that can improve the efficiency of Cartesian mesh generation, achieve rapid mesh adjustment without manipulating geometric files, thereby avoiding repeated mesh generation, reducing the manual and time costs of mesh generation, and providing a guarantee for the on-the-spot numerical modeling of wind fields during ship operations. The specific solution is as follows:

[0005] In a first aspect, this application discloses a method for generating Cartesian grids, including:

[0006] Obtain the mesh generation instruction and determine the type of the mesh generation instruction; wherein, the type includes generating a new mesh and updating the mesh;

[0007] If the type is to generate a new mesh, then the geometry file information and construction parameter information are determined based on the mesh generation instruction. A Cartesian mesh is generated based on the construction parameter information, the configuration information of the Cartesian mesh is determined, fluid marker points and mesh points are determined, the number of intersections of the fluid marker points is calculated to obtain the number of intersections, and the number of intersections is used to determine whether the mesh points match the geometry file information. If the mesh points match the geometry file information, the construction parameter information is updated to obtain updated construction parameter information. The Cartesian mesh is marked based on the updated construction parameter information and the geometry file information to obtain a new mesh.

[0008] If the type of the mesh generation instruction is to adjust the historical mesh, then the target historical mesh that matches the mesh generation instruction is selected from the historical meshes, the mesh adjustment parameters are determined based on the mesh generation instruction, and the target historical mesh is adjusted based on the mesh adjustment parameters to obtain the adjusted target historical mesh.

[0009] Optionally, determining the geometric file information and construction parameter information based on the mesh generation instructions, and generating a Cartesian mesh based on the construction parameter information, includes:

[0010] Based on the mesh generation instructions, geometric file information and construction parameter information are determined, and the number of geometric files is determined; wherein, the construction parameter information includes the coordinate range and precision information of the Cartesian mesh;

[0011] A Cartesian grid is generated within the coordinate range in a uniformly discrete manner, and the total number of grid points is recorded.

[0012] Optionally, determining the configuration information of the Cartesian grid, and determining the fluid marker points and grid points, includes:

[0013] Attribute labeling and equipment type labeling are configured for all grid points in the Cartesian grid to obtain the configuration information;

[0014] The watershed area is determined, and fluid markers are randomly selected within the watershed area, along with grid points.

[0015] Optionally, the step of calculating the number of intersections of the fluid marker points to obtain the number of intersections includes:

[0016] The geometric file model is obtained, and the device environment of the Turing architecture graphics processor is initialized to obtain the initialized Turing architecture graphics processor. The environment of the geometric file model is traversed using the initialized Turing architecture graph. Then, the geometric file model is divided into bounding boxes. The number of intersections between the line connecting each mesh and fluid marker point and the geometric triangle face in the geometric file model is calculated using RTCore hardware to obtain the number of intersections.

[0017] Optionally, determining whether the grid point matches the geometry file information based on the number of intersections, and updating the construction parameter information if the grid point matches the geometry file information, includes:

[0018] Determine the type of the number of intersection points; wherein the type includes odd and even numbers;

[0019] If the number of intersection points is odd, the grid point is determined to match the geometric file information, and the construction parameter information is then updated. If the number of intersection points is even, the grid point is determined not to match the geometric file information, and the process jumps directly to the step of marking the Cartesian grid based on the updated construction parameter information and the geometric file information.

[0020] Optionally, determining the mesh adjustment parameters based on the mesh generation command and adjusting the target historical mesh based on the mesh adjustment parameters includes:

[0021] Obtain a grid copy corresponding to the target historical grid, and open the grid copy in a read-write manner;

[0022] Based on the mesh generation command, mesh adjustment parameters are determined, and the mesh copy is adjusted according to the mesh adjustment parameters.

[0023] Optionally, determining the mesh adjustment parameters based on the mesh generation instruction and adjusting the mesh copy according to the mesh adjustment parameters includes:

[0024] The grid adjustment parameters are determined based on the grid generation instructions; wherein, the grid adjustment parameters include the number of devices to be adjusted, the Cartesian grid point number, and the attribute marker;

[0025] The grid copy is adjusted according to the grid adjustment parameters to obtain the adjusted grid copy, and the grid copy is saved to obtain the adjusted target historical grid.

[0026] Secondly, this application discloses a Cartesian grid generation apparatus, comprising:

[0027] The instruction acquisition module is used to acquire mesh generation instructions and determine the type of the mesh generation instructions; wherein, the type includes generating a new mesh and updating a mesh;

[0028] The new mesh generation module is used to, if the type is to generate a new mesh, determine the geometry file information and construction parameter information based on the mesh generation instruction, generate a Cartesian mesh based on the construction parameter information, determine the configuration information of the Cartesian mesh, determine fluid marker points and mesh points, calculate the number of intersections of the fluid marker points to obtain the number of intersections, determine whether the mesh points match the geometry file information based on the number of intersections, update the construction parameter information to obtain updated construction parameter information, and mark the Cartesian mesh based on the updated construction parameter information and the geometry file information to obtain a new mesh;

[0029] The historical grid adjustment module is used to, if the type of the grid generation instruction is to adjust the historical grid, select a target historical grid that matches the grid generation instruction from the historical grid, determine the grid adjustment parameters based on the grid generation instruction, and adjust the target historical grid based on the grid adjustment parameters to obtain the adjusted target historical grid.

[0030] Thirdly, this application discloses an electronic device, including:

[0031] Memory, used to store computer programs;

[0032] A processor is used to execute the computer program to implement the aforementioned Cartesian grid generation method.

[0033] Fourthly, this application discloses a computer storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the steps of the aforementioned disclosed Cartesian grid generation method.

[0034] As can be seen, this application provides a Cartesian mesh generation method, including obtaining a mesh generation instruction and determining the type of the mesh generation instruction; wherein, the type includes generating a new mesh and updating a mesh; if the type is generating a new mesh, then based on the mesh generation instruction, geometric file information and construction parameter information are determined, a Cartesian mesh is generated based on the construction parameter information, the configuration information of the Cartesian mesh is determined, fluid marker points and mesh points are determined, the number of intersections of the fluid marker points is calculated to obtain the number of intersections, and based on the number of intersections, it is determined whether the mesh points match the geometric file information; if the mesh points match the geometric file information, the construction parameter information is updated to obtain updated construction parameter information, and the Cartesian mesh is marked based on the updated construction parameter information and the geometric file information to obtain a new mesh; if the type of the mesh generation instruction is adjusting a historical mesh, then a target historical mesh that matches the mesh generation instruction is selected from the historical meshes, mesh adjustment parameters are determined based on the mesh generation instruction, and the target historical mesh is adjusted based on the mesh adjustment parameters to obtain the adjusted target historical mesh. This application utilizes the lattice Boltzmann method for wind field calculations, primarily employing Cartesian grids to more accurately capture flow field structures. Addressing the rapid response requirements of digital maritime warfare, it enables the rapid generation of grids needed for wind field calculations with minimal human intervention. The software calculation process leverages Turing architecture GPU hardware, allowing for quick adjustment and creation of new geometries without modifying geometry files or recalculating grid generation. This reduces the manual and time costs of grid generation, making it suitable for wind field calculations on naval vessels. It minimizes the impact of newly added geometric equipment and minor changes in equipment shape on the wind field, ensuring the ability to perform wind field numerical models at any time during ship operations. Attached Figure Description

[0035] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0036] Figure 1 This is a flowchart of a Cartesian mesh generation method disclosed in this application;

[0037] Figure 2 Here is an example diagram of the two-dimensional Cartesian mesh disclosed in this application;

[0038] Figure 3 This is a schematic diagram of the automatic labeling algorithm disclosed in this application;

[0039] Figure 4 This is a flowchart of the hardware-accelerated alternative to software computing disclosed in this application;

[0040] Figure 5 This is a diagram of the OptiX program relationships disclosed in this application;

[0041] Figure 6 This is an example diagram of the STL format geometry disclosed in this application;

[0042] Figure 7 This is a flowchart illustrating a Cartesian grid generation method disclosed in this application.

[0043] Figure 8 This is a schematic diagram of a Cartesian grid generation device disclosed in this application;

[0044] Figure 9 This application provides a structural diagram of an electronic device. Detailed Implementation

[0045] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0046] Currently, wind field environmental forecasting, as an integral part of digital battlefield construction, plays a crucial guiding role in naval operations and carrier-based aircraft takeoffs and landings. However, due to the dynamic changes in the equipment layout on the ship's deck and the position and number of aircraft, the mesh for dynamic wind field simulation needs frequent changes. When simulating wind fields on ships, the first step is to generate the necessary mesh based on the geometric model. However, ships have complex shapes and numerous pieces of equipment on their decks (such as carrier-based aircraft, the island superstructure, and measuring instruments), making conventional mesh generation methods labor-intensive and costly. Furthermore, the equipment layout on the deck can change at any time during ship operations (e.g., the number and position of carrier-based aircraft), making mesh regeneration extremely expensive and hindering the rapid calculation of wind field results. Therefore, improving the efficiency of Cartesian mesh generation, enabling rapid mesh adjustments without manipulating geometric files, avoiding repeated mesh generation, reducing the labor and time costs of mesh generation, and providing a guarantee for real-time wind field numerical modeling during ship operations are problems that need to be solved in this field.

[0047] See Figure 1 As shown in the figure, an embodiment of the present invention discloses a method for generating Cartesian grids, which may specifically include:

[0048] Step S11: Obtain the mesh generation instruction and determine the type of the mesh generation instruction; wherein, the type includes generating a new mesh and updating the mesh.

[0049] Step S12: If the type is to generate a new mesh, then the geometry file information and construction parameter information are determined based on the mesh generation instruction. A Cartesian mesh is generated based on the construction parameter information, the configuration information of the Cartesian mesh is determined, fluid marker points and mesh points are determined, the number of intersections of the fluid marker points is calculated to obtain the number of intersections, and the mesh points are judged to match the geometry file information based on the number of intersections. If the mesh points match the geometry file information, the construction parameter information is updated to obtain updated construction parameter information. The Cartesian mesh is marked based on the updated construction parameter information and the geometry file information to obtain a new mesh.

[0050] In this embodiment, if the type is to generate a new mesh, then the geometric file information and construction parameter information are determined based on the mesh generation instruction, and the number of geometric files is determined. The construction parameter information includes the coordinate range and precision information of the Cartesian mesh. A Cartesian mesh is generated within the coordinate range in a uniformly discrete manner, and the total number of mesh points is recorded. Then, attribute labeling and equipment type labeling configurations are performed on all mesh points in the Cartesian mesh to obtain the configuration information. The watershed range is determined, and fluid marker points of the fluid domain are randomly determined within the watershed range. Mesh points are also determined. Then, the geometric file model is obtained, and the device environment of the Turing architecture graphics processor is initialized to obtain the initial... The initialized Turing architecture graphics processor performs environment traversal on the geometric file model using the initialized Turing architecture graph. Then, it divides the geometric file model into bounding boxes and uses RTCore hardware to calculate the number of intersections between the line connecting each mesh and fluid marker point and the geometric triangle face in the geometric file model. Based on this number of intersections, it determines whether the mesh point matches the geometric file information. If the mesh point matches the geometric file information, the construction parameter information is updated to obtain updated construction parameter information. Based on the updated construction parameter information and the geometric file information, the Cartesian mesh is marked to obtain a new mesh.

[0051] This application primarily uses Cartesian grids for wind field calculations using the Lattice Boltzmann Method (LBM). Its advantages lie mainly in the generation of spatial non-body grids with low dependence on the object surface. This generation method makes the grid generation process simple, fast, and highly automated. Secondly, the Cartesian grid's ease of subdivision combined with the cross-tree data structure facilitates adaptive behavior, thus enabling more accurate capture of the flow field structure. The calculations typically require labeling fluid, solid, and boundary regions. A common method is to mark the grid using coordinate ranges or manually, such as... Figure 2 As shown in the figure, 0 represents fluid points and 1 represents solid points, indicating the flow domain and solid region. In the figure, (1) is the geometric image of the two-dimensional Cartesian mesh example, (2) is the generated Cartesian mesh image, and (3) is the attribute label image. This method is easy to implement in two-dimensional or regular boundaries, but for irregular geometry in three-dimensional real-world complex applications, when the number of geometry is very large and large-scale parallel computing is required, the overhead of this manual labeling method is very large. In addition, the equipment layout scheme on the ship deck, the position and number of aircraft change frequently, so it is necessary to re-determine and regenerate, which brings a great obstacle to the rapid numerical modeling during operation. When performing external flow field and wind field calculations in any scenario, this application needs to use a sufficiently large flow domain to surround the obstacles to set the boundary of the external flow field. In order to ensure the correct calculation of flow details, the boundary distance from the obstacle usually needs to be sufficient to ensure the full development of the fluid. This application adopts an automatic labeling algorithm, and the specific algorithm principle is as follows. Figure 3 As shown, for wind field calculations involving obstacles of arbitrary shapes, a marker point (the sun marker in the figure) within the fluid domain can be randomly selected at the upper boundary of the adjacent computational domain grid. For grid points within the computational domain, the type of grid point can be determined by calculating the number of intersections between the line connecting the marker point and the grid point and the geometric surface, based on this fluid domain marker point. For example... Figure 3 As shown at point A, for a grid point in the fluid domain, the number of intersections between the line connecting it to the marker point and the geometric surface must be even. For a solid point, represented by point B, the number of intersections must be odd. For a solid point, if its adjacent point in any direction is a fluid point, it is marked as a boundary point. This application transforms the marking of the fluid domain into a problem of the intersection of light rays and triangles. The sun marker point is the starting point of the light ray, and the line connecting each grid point and the marker serves as the direction and endpoint of a light ray.

[0052] To quickly solve the problem of the intersection of light rays and triangles, the intersection determination in ship calculations can be mainly divided into two parts: (1) traversal determination of bounding boxes: For ships, we can build complete STL models of the ship, each carrier-based aircraft, and each piece of equipment separately. First, we need to calculate the interaction between a light ray and these large entities. This step is to divide the space where the ship and equipment are located into multiple bounding boxes and first determine whether the light ray intersects with these large bounding boxes. This step is the intersection determination of light rays and cuboids, which is relatively easy to handle. We only need to compare the coordinates of the endpoints of the light rays with the coordinate range of the cuboids to determine whether the points on the light rays may fall into the cuboids. If so, we can further divide the bounding box of the large entity into N smaller bounding boxes and calculate the intersection again. This process is repeated until the bounding box of the triangle that intersects with the light ray is found, and then the final intersection determination of the light ray and triangle is performed on this triangle. This method can quickly filter out non-intersecting triangles. (2) For a ray in space, let the starting point be O and its direction be D. According to the parameter formula of the ray, any point P on it (that is, the intersection point to be found) is as follows:

[0053] P = O + tD

[0054] For a spatial triangle with vertices V1, V2, and V3, for any point inside the triangle, the following parametric equation holds:

[0055] P = (1-uv)V1 + uV2 + vV3

[0056] u and v are the weights of V2 and V3, 1-uv is the weight of V1, and u≥0, v≥0, u+v≤1 are satisfied.

[0057] Combining the formulas for rays and triangles, we get: (1-uv)V1+uV2+vV3=O+tD

[0058] Clearly, u, v, and t are all unknowns. By rearranging and simplifying, we obtain the following system of linear equations:

[0059]

[0060] Solving this system of equations allows us to determine the intersection of the spatial ray and the triangle, and to obtain the location of the intersection point.

[0061] This application employs a Turing architecture GPU (Graphics Processing Unit) ray tracing structure to achieve automated mesh marking; such as Figure 4As shown, compared to previous generations of GPUs, the latest Turing architecture GPUs, in order to better perform 3D graphics rendering, have added dedicated RTCores (Ray-Tracing Cores) to the Streaming Multiprocessor (SM) for ray tracing. The functions of traversing and intersecting bounding boxes and determining intersections with triangles are integrated as fixed hardware units within the RTCore. Therefore, the RTCore allows software implementations to be transformed into hardware implementations, resulting in greater speed and efficiency. GPU hardware supporting this functionality includes, but is not limited to: NVIDIA GeForce RTX 20 series, NVIDIA GeForce RTX 30 series, NVIDIA Quadro RTX, and AMD Radeon RX 6000 series. Taking NVIDIA's OptiX ray tracing engine as an example, this ray tracing engine is free for commercial use and provides GPU-accelerated ray tracing through the inherent RTCores of RTX GPUs. Its basic architecture is as follows... Figure 5 As shown, Ray generation, Continuation callable, Miss, Closest-hit, Any-hit, and Direct callable require user-defined implementations. These parts allow for additional programmable functionality within the standard set of OptiX programs, thus also enabling utilization in non-graphics rendering domains; while Scratch traversal and Intersection are implemented by OptiX. Furthermore, the entire logic layer is also handled by OptiX. The ray generation program generally consists of three parts: creating rays, initiating scene traversal, and storing ray tracing results. The intersection procedure implements ray geometry intersection testing. When a ray does not intersect any geometric object in the scene, the Miss procedure is executed, indicating no intersection with the geometry. When the tracing ray finds a new intersection, Any-hit is called, indicating that the ray intersects with geometry in the scene. This application can achieve rapid generation, automated discrimination, easy storage, and efficient adjustment of 3D Cartesian meshes based on stereolithography format geometry files. STL (STereoLithography) is a file format (i.e., geometric file) created by 3D Systems software company, originally used for stereolithography computer-aided design software. An example of an STL format geometry is shown below. Figure 6 As shown, STL is composed of triangular faces and can only be used to represent closed faces or volumes.

[0062] This application utilizes the efficient storage and rapid adjustment of mesh information to create new meshes. When using the mesh, different equipment (such as carrier-based aircraft and mobile measuring instruments) on the ship are modeled separately. Each individual mobile instrument uses a separate STL format file. Since ships typically do not add new equipment during operations at sea, all equipment at the time of departure can be individually labeled. Assuming the ship's main body is labeled 0, carrier-based aircraft 1, and a mobile measuring instrument 2, the solid point labeled 1 can be subdivided into solid points inside the ship's hull (1,0), solid points inside the carrier-based aircraft (1,1), and points inside the mobile measuring instrument (1,2), etc., where the first digit 1 indicates a solid area, and the second digit indicates the equipment type. After the mesh is generated and labeled, the 3D coordinates, attribute labels, and equipment type labels of each point are recorded and formatted into a file for storage. When reusing the mesh, if mesh adjustments are needed, a copy of the original mesh information can be made, and a new mesh can be quickly created by retrieving the equipment type label. For example, when a carrier-based aircraft takes off and leaves the ship, the calculation only requires finding the grid point marked (1,1) in the copy, changing its corresponding coordinates to (0, -1), and then re-finding and changing it back to (1,1) when it returns to its docking position. When the position of movable equipment changes, only the attribute and equipment type markings of the grid points at the initial and final positions need to be modified. For example, when a movable measuring equipment moves 10 meters in the X direction, only the corresponding grid point marked (1,2) needs to be found, its attribute marking set to 0, and the corresponding grid point at 10 meters in the X direction needs to be marked (1,2). Changes in direction and rotation are handled similarly. This method avoids adjusting the geometry file and generates the required mesh without recalculation.

[0063] Step S13: If the type of the mesh generation instruction is to adjust the historical mesh, then select the target historical mesh that matches the mesh generation instruction from the historical mesh, determine the mesh adjustment parameters based on the mesh generation instruction, and adjust the target historical mesh based on the mesh adjustment parameters to obtain the adjusted target historical mesh.

[0064] In this embodiment, a mesh generation instruction is obtained, and the type of the mesh generation instruction is determined. The type includes generating a new mesh and updating a mesh. If the type is generating a new mesh, geometric file information and construction parameter information are determined based on the mesh generation instruction. A Cartesian mesh is generated based on the construction parameter information, and the configuration information of the Cartesian mesh is determined. Fluid marker points and mesh points are determined, and the number of intersections of the fluid marker points is calculated. Based on the number of intersections, it is determined whether the mesh points match the geometric file information. If the mesh points match the geometric file information, the construction parameter information is updated to obtain updated construction parameter information. The Cartesian mesh is marked based on the updated construction parameter information and the geometric file information to obtain a new mesh. If the type of the mesh generation instruction is adjusting a historical mesh, a target historical mesh matching the mesh generation instruction is selected from the historical meshes. Mesh adjustment parameters are determined based on the mesh generation instruction, and the target historical mesh is adjusted based on the mesh adjustment parameters to obtain the adjusted target historical mesh. This application utilizes the lattice Boltzmann method for wind field calculations, primarily employing Cartesian grids to more accurately capture flow field structures. Addressing the rapid response requirements of digital maritime warfare, it enables the rapid generation of grids needed for wind field calculations with minimal human intervention. The software calculation process leverages Turing architecture GPU hardware, allowing for quick adjustment and creation of new geometries without modifying geometry files or recalculating grid generation. This reduces the manual and time costs of grid generation, making it suitable for wind field calculations on naval vessels. It minimizes the impact of newly added geometric equipment and minor changes in equipment shape on the wind field, ensuring the ability to perform wind field numerical models at any time during ship operations.

[0065] See Figure 7 As shown in the figure, an embodiment of the present invention discloses a method for generating Cartesian grids, which may specifically include:

[0066] Step S21: Obtain the mesh generation instruction and determine the type of the mesh generation instruction; wherein the type includes generating a new mesh and updating the mesh.

[0067] Step S22: If the type is to generate a new mesh, then the geometry file information and construction parameter information are determined based on the mesh generation instruction. A Cartesian mesh is generated based on the construction parameter information. The configuration information of the Cartesian mesh is determined. Fluid marker points and mesh points are determined. The number of intersections of the fluid marker points is calculated to obtain the number of intersections. Based on the number of intersections, it is determined whether the mesh points match the geometry file information. If the mesh points match the geometry file information, the construction parameter information is updated to obtain updated construction parameter information. The Cartesian mesh is marked based on the updated construction parameter information and the geometry file information to obtain a new mesh.

[0068] In this embodiment, the specific process for determining whether the grid point matches the geometric file information is as follows: determine the type of the number of intersection points; wherein, the type includes odd and even numbers; if the type of the number of intersection points is odd, then it is determined that the grid point matches the geometric file information, and then the construction parameter information is updated; if the type of the number of intersection points is even, then it is determined that the grid point does not match the geometric file information, and then directly jump to the step of marking the Cartesian mesh based on the updated construction parameter information and the geometric file information.

[0069] Step S23: If the type of the mesh generation instruction is to adjust the historical mesh, then select the target historical mesh that matches the mesh generation instruction from the historical mesh, and determine the mesh adjustment parameters based on the mesh generation instruction.

[0070] Step S24: Obtain a grid copy corresponding to the target historical grid, open the grid copy in read-write mode, determine the grid adjustment parameters based on the grid generation instruction, and adjust the grid copy according to the grid adjustment parameters to obtain the adjusted target historical grid.

[0071] In this embodiment, the specific process of adjusting the grid copy is as follows: determining the grid adjustment parameters based on the grid generation instruction; wherein, the grid adjustment parameters include the number of devices to be adjusted, the Cartesian grid point number, and the attribute marker; adjusting the grid copy according to the grid adjustment parameters to obtain the adjusted grid copy, and saving the grid copy to obtain the adjusted target historical grid.

[0072] The specific working steps of this application are as follows: (1) According to the mesh generation instruction, determine whether a new mesh needs to be generated. If so, proceed to (2) the process of generating a new mesh. If not, proceed to (29) the process of adjusting the generation history mesh. (2) Read in the storage location of the geometry file and the coordinate range and precision of the Cartesian mesh to be generated. Count the number of geometry files as I. (3) Generate Cartesian meshes of corresponding precision uniformly within the coordinate range. Record the total number of mesh points N. Set the attribute markers and equipment type markers of all mesh points to (0, -1) by default. Where 0 represents the fluid region and -1 represents no equipment. (4) Randomly select a reference marker of the fluid region near the inlet of the watershed. (5) Initialize i = 0, where i represents the number of the geometry file; (6) Load the i-th 3D STL format geometry file; (7) Initialize the device environment of the Turing architecture GPU; (8) Initialize n = 0, where n represents the number of the Cartesian mesh point; (9) Generate ray n, the starting point of the ray is the reference marker point selected in (4), and the ending point of the ray is the position of the n-th mesh point; (10) The GPU performs environment traversal on the geometry within the range; (11) The RTCore hardware automatically performs hierarchical bounding box traversal judgment; (12) The RTCore hardware automatically performs spatial ray and triangle intersection judgment; (13) Count the number of hits, and judge the parity of the number of hits. (14) If not inside a solid, proceed to (15), otherwise proceed to (16); (15) If determined to be a fluid point, do not change the mesh label; (16) If determined to be a solid point, set the attribute label and equipment type label to (1, n); (17) Let n = n + 1; (18) Determine if n is less than the total number of mesh points N. If yes, proceed to (9), otherwise proceed to step (19); (19) Unload the i-th geometry file and let i = i + 1; (20) Determine if i is less than the total number of geometry files I. If yes, proceed to (6), otherwise proceed to (21); (21) Create and open the mesh file and write the formatting information; (22) Reinitialize n = 0, n still represents the Cartesian mesh. Grid point number; (23) Determine if the nth grid point is a solid point and its neighboring grid points are fluid points. If so, execute (24), otherwise execute (25); (24) Determine if it is a boundary point, modify the attribute marker of the nth grid point to 2, and leave the equipment type marker unchanged; (25) Write the three-dimensional coordinates, attribute marker, and equipment type marker of the nth grid point into the grid file; (26) Let n = n + 1; (27) Determine if n is less than the total number of grid points N. If so, execute (22), otherwise execute (28); (28) Close the grid file and terminate the process; Step 29: Read the grid file address, copy a grid file copy, open the original grid file in read-only mode, and open the grid file copy in read-write mode;Determine whether mesh adjustment is needed based on the input. If so, execute (31), operate and modify the mesh file copy; otherwise, execute (46). (31) Determine whether devices need to be reduced / added based on the input. If so, read in the number of devices to be reduced / added, S, and the corresponding device number list, and execute (32); otherwise, execute (41). (32) Initialize s = 0, where s represents the s-th device to be reduced / added. (33) Adjust the attribute flag of the s-th device using a NOT gate operation. If it is 1 or 2, set it to 0; if it is 0, set it to 1. (34) Let s = s + 1. (35) Determine whether s is less than the number of devices to be reduced / added, S. If so, execute (33); otherwise, execute (36). (36) Initialize n = 0, where n represents the Cartesian mesh point number. (37) Determine whether the n-th mesh point is a solid point and its neighboring mesh is a fluid point. If so, execute ( 38), otherwise execute step; (38) if it is a boundary point, modify the attribute mark of the nth grid point to 2, and keep the equipment type mark unchanged; (39) let n = n + 1; (40) determine if n is less than the total number of grid points N, if yes, execute (37), otherwise execute (41); (41) determine whether the equipment needs to be translated / rotated according to the input, if yes, read the number of equipment to be translated / rotated S and the corresponding equipment number list, execute step 42, otherwise execute step Y; (42) initialize s = 0, s represents the sth equipment to be translated / rotated; (43) modify the attribute and equipment type mark of the initial and final positions of the sth grid point; (44) let s = s + 1; (45) determine if s is less than the number of equipment to be translated / rotated S, if yes, execute step (43), otherwise execute (46); (46) save the file copy, close the original grid file, and the process terminates.

[0073] In this embodiment, a mesh generation instruction is obtained, and the type of the mesh generation instruction is determined. The type includes generating a new mesh and updating a mesh. If the type is generating a new mesh, geometric file information and construction parameter information are determined based on the mesh generation instruction. A Cartesian mesh is generated based on the construction parameter information, and the configuration information of the Cartesian mesh is determined. Fluid marker points and mesh points are determined, and the number of intersections of the fluid marker points is calculated. Based on the number of intersections, it is determined whether the mesh points match the geometric file information. If the mesh points match the geometric file information, the construction parameter information is updated to obtain updated construction parameter information. The Cartesian mesh is marked based on the updated construction parameter information and the geometric file information to obtain a new mesh. If the type of the mesh generation instruction is adjusting a historical mesh, a target historical mesh matching the mesh generation instruction is selected from the historical meshes. Mesh adjustment parameters are determined based on the mesh generation instruction, and the target historical mesh is adjusted based on the mesh adjustment parameters to obtain the adjusted target historical mesh. This application utilizes the lattice Boltzmann method for wind field calculations, primarily employing Cartesian grids to more accurately capture flow field structures. Addressing the rapid response requirements of digital maritime warfare, it enables the rapid generation of grids needed for wind field calculations with minimal human intervention. The software calculation process leverages Turing architecture GPU hardware, allowing for quick adjustment and creation of new geometries without modifying geometry files or recalculating grid generation. This reduces the manual and time costs of grid generation, making it suitable for wind field calculations on naval vessels. It minimizes the impact of newly added geometric equipment and minor changes in equipment shape on the wind field, ensuring the ability to perform wind field numerical models at any time during ship operations.

[0074] See Figure 8 As shown, this embodiment of the invention discloses a Cartesian grid generation device, which may specifically include:

[0075] The instruction acquisition module 11 is used to acquire mesh generation instructions and determine the type of the mesh generation instructions; wherein, the type includes generating a new mesh and updating a mesh;

[0076] The new mesh generation module 12 is used to, if the type is to generate a new mesh, determine the geometry file information and construction parameter information based on the mesh generation instruction, generate a Cartesian mesh based on the construction parameter information, determine the configuration information of the Cartesian mesh, determine fluid marker points and mesh points, calculate the number of intersections of the fluid marker points to obtain the number of intersections, determine whether the mesh points match the geometry file information based on the number of intersections, update the construction parameter information to obtain updated construction parameter information, and mark the Cartesian mesh based on the updated construction parameter information and the geometry file information to obtain a new mesh;

[0077] The historical grid adjustment module 13 is used to, if the type of the grid generation instruction is to adjust the historical grid, select a target historical grid that matches the grid generation instruction from the historical grid, determine the grid adjustment parameters based on the grid generation instruction, and adjust the target historical grid based on the grid adjustment parameters to obtain the adjusted target historical grid.

[0078] In this embodiment, a mesh generation instruction is obtained, and the type of the mesh generation instruction is determined. The type includes generating a new mesh and updating a mesh. If the type is generating a new mesh, geometric file information and construction parameter information are determined based on the mesh generation instruction. A Cartesian mesh is generated based on the construction parameter information, and the configuration information of the Cartesian mesh is determined. Fluid marker points and mesh points are determined, and the number of intersections of the fluid marker points is calculated. Based on the number of intersections, it is determined whether the mesh points match the geometric file information. If the mesh points match the geometric file information, the construction parameter information is updated to obtain updated construction parameter information. The Cartesian mesh is marked based on the updated construction parameter information and the geometric file information to obtain a new mesh. If the type of the mesh generation instruction is adjusting a historical mesh, a target historical mesh matching the mesh generation instruction is selected from the historical meshes. Mesh adjustment parameters are determined based on the mesh generation instruction, and the target historical mesh is adjusted based on the mesh adjustment parameters to obtain the adjusted target historical mesh. This application utilizes the lattice Boltzmann method for wind field calculations, primarily employing Cartesian grids to more accurately capture flow field structures. Addressing the rapid response requirements of digital maritime warfare, it enables the rapid generation of grids needed for wind field calculations with minimal human intervention. The software calculation process leverages Turing architecture GPU hardware, allowing for quick adjustment and creation of new geometries without modifying geometry files or recalculating grid generation. This reduces the manual and time costs of grid generation, making it suitable for wind field calculations on naval vessels. It minimizes the impact of newly added geometric equipment and minor changes in equipment shape on the wind field, ensuring the ability to perform wind field numerical models at any time during ship operations.

[0079] In some specific embodiments, the instruction acquisition module 11 may specifically include:

[0080] The information determination module is used to determine geometric file information and construction parameter information based on the mesh generation instructions, and to determine the number of geometric files; wherein, the construction parameter information includes the coordinate range and precision information of the Cartesian mesh;

[0081] The Cartesian grid generation module is used to generate a Cartesian grid in a uniformly discrete manner within the coordinate range and record the total number of grid points.

[0082] In some specific embodiments, the new mesh generation module 12 may specifically include:

[0083] The configuration information determination module is used to perform attribute marking configuration and equipment type marking configuration on all grid points in the Cartesian grid to obtain the configuration information;

[0084] The grid point determination module is used to determine the watershed range, randomly determine the fluid marker points of the fluid domain within the watershed range, and determine the grid points.

[0085] In some specific embodiments, the new mesh generation module 12 may specifically include:

[0086] An initialization module is used to acquire a geometric file model, initialize the device environment of the Turing architecture graphics processor to obtain the initialized Turing architecture graphics processor, perform environment traversal on the geometric file model using the initialized Turing architecture graph, then divide the geometric file model into bounding boxes, and use RTCore hardware to calculate the number of intersection points between the line connecting each mesh and fluid marker point and the geometric triangle face in the geometric file model to obtain the number of intersection points.

[0087] In some specific embodiments, the new mesh generation module 12 may specifically include:

[0088] A type determination module is used to determine the type of the number of intersection points; wherein the type includes odd and even numbers;

[0089] The construction parameter information update module is used to determine if the type of the number of intersection points is odd, and then update the construction parameter information; if the type of the number of intersection points is even, and then determine if the grid points do not match the geometric file information, then directly jump to the step of marking the Cartesian mesh based on the updated construction parameter information and the geometric file information.

[0090] In some specific embodiments, the historical grid adjustment module 13 may specifically include:

[0091] The grid copy acquisition module is used to acquire a grid copy corresponding to the target historical grid and open the grid copy in a read-write manner.

[0092] The grid copy adjustment module is used to determine grid adjustment parameters based on the grid generation command, and adjust the grid copy according to the grid adjustment parameters.

[0093] In some specific embodiments, the historical grid adjustment module 13 may specifically include:

[0094] The grid adjustment parameter determination module is used to determine the grid adjustment parameters based on the grid generation command; wherein, the grid adjustment parameters include the number of devices to be adjusted, the Cartesian grid point number, and the attribute marker;

[0095] The grid copy adjustment module is used to adjust the grid copy according to the grid adjustment parameters to obtain the adjusted grid copy, and save the grid copy to obtain the adjusted target historical grid.

[0096] Figure 9 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the Cartesian grid generation method performed by the electronic device as disclosed in any of the foregoing embodiments.

[0097] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0098] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored thereon include operating system 221, computer program 222 and data 223, etc., and the storage method can be temporary storage or permanent storage.

[0099] The operating system 221 manages and controls the various hardware devices on the electronic device 20 and the computer program 222 to enable the processor 21 to perform calculations and processing on the data 223 in the memory 22. The operating system 221 can be Windows, Unix, Linux, etc. The computer program 222, in addition to including a computer program capable of performing the Cartesian grid generation method executed by the electronic device 20 as disclosed in any of the foregoing embodiments, may further include computer programs capable of performing other specific tasks. The data 223 may include data received by the Cartesian grid generation device from external devices, as well as data collected by its own input / output interface 25.

[0100] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0101] Furthermore, embodiments of this application also disclose a computer-readable storage medium storing a computer program, which, when loaded and executed by a processor, implements the Cartesian grid generation method steps disclosed in any of the foregoing embodiments.

[0102] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0103] The present invention provides a detailed description of a Cartesian grid generation method, apparatus, device, and storage medium. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for generating Cartesian grids, characterized in that, include: Obtain the mesh generation instruction and determine the type of the mesh generation instruction; wherein, the type includes generating a new mesh and updating the mesh; If the type is to generate a new mesh, then the geometry file information and construction parameter information are determined based on the mesh generation instruction. A Cartesian mesh is generated based on the construction parameter information, the configuration information of the Cartesian mesh is determined, fluid marker points and mesh points are determined, the number of intersections of the fluid marker points is calculated to obtain the number of intersections, and based on the number of intersections, it is determined whether the mesh points match the geometry file information. If the mesh points match the geometry file information, the construction parameter information is updated to obtain updated construction parameter information. Based on the updated construction parameter information and the geometry file information, the Cartesian mesh is marked to obtain a new mesh; wherein, the Cartesian mesh is the mesh used when performing wind field calculations using the lattice Boltzmann method; If the type of the mesh generation instruction is to adjust the historical mesh, then the target historical mesh that matches the mesh generation instruction is selected from the historical mesh, the mesh adjustment parameters are determined based on the mesh generation instruction, and the target historical mesh is adjusted based on the mesh adjustment parameters to obtain the adjusted target historical mesh. The configuration information of the Cartesian grid is determined, including the determination of fluid marker points and grid points. This includes: configuring attribute markings and equipment type markings for all grid points in the Cartesian grid to obtain configuration information; determining the watershed range, randomly determining fluid marker points within the watershed range, and determining grid points; where grid points are points in the Cartesian grid; and fluid marker points are randomly selected marker points within the watershed range that are located within the fluid watershed. The number of intersections of fluid markers is calculated, including: acquiring the geometric file model; initializing the device environment of the Turing architecture graphics processor to obtain the initialized Turing architecture graphics processor; using the initialized Turing architecture graphics processor to traverse the environment of the geometric file model; then dividing the geometric file model into bounding boxes; and using RTCore hardware to calculate the number of intersections between the line connecting each mesh and fluid marker and the geometric triangle face in the geometric file model. An intersection is defined as the point formed by the intersection of the line connecting the mesh point and the fluid marker with the geometric triangle face in the geometric file model. Based on the mesh generation command, mesh adjustment parameters are determined, and the target historical mesh is adjusted according to the mesh adjustment parameters, including: obtaining a mesh copy corresponding to the target historical mesh and opening the mesh copy in read-write mode; based on the mesh generation command, mesh adjustment parameters are determined, and the mesh copy is adjusted according to the mesh adjustment parameters; The grid adjustment parameters are determined based on the grid generation command, and the grid copy is adjusted according to the grid adjustment parameters, including: determining the grid adjustment parameters based on the grid generation command; wherein, the grid adjustment parameters include the number of devices to be adjusted, the Cartesian grid point number, and the attribute marker; adjusting the grid copy according to the grid adjustment parameters to obtain the adjusted grid copy, and saving the grid copy to obtain the adjusted target historical grid.

2. The Cartesian mesh generation method according to claim 1, characterized in that, The process of determining geometric file information and construction parameter information based on the mesh generation instructions, and generating a Cartesian mesh based on the construction parameter information, includes: Based on the mesh generation instructions, geometric file information and construction parameter information are determined, and the number of geometric files is determined; wherein, the construction parameter information includes the coordinate range and precision information of the Cartesian mesh; A Cartesian grid is generated within the coordinate range in a uniformly discrete manner, and the total number of grid points is recorded.

3. The Cartesian grid generation method according to claim 1, characterized in that, The step of determining whether the grid point matches the geometry file information based on the number of intersections, and updating the construction parameter information if the grid point matches the geometry file information, includes: Determine the type of the number of intersection points; wherein the type includes odd and even numbers; If the number of intersection points is odd, the grid point is determined to match the geometric file information, and the construction parameter information is then updated. If the number of intersection points is even, the grid point is determined not to match the geometric file information, and the process jumps directly to the step of marking the Cartesian grid based on the updated construction parameter information and the geometric file information.

4. A Cartesian grid generation device, characterized in that, include: The instruction acquisition module is used to acquire mesh generation instructions and determine the type of the mesh generation instructions; wherein, the type includes generating a new mesh and updating a mesh; The new mesh generation module is used to, if the type is to generate a new mesh, determine the geometric file information and construction parameter information based on the mesh generation instruction, generate a Cartesian mesh based on the construction parameter information, determine the configuration information of the Cartesian mesh, determine fluid marker points and mesh points, calculate the number of intersections of the fluid marker points to obtain the number of intersections, determine whether the mesh points match the geometric file information based on the number of intersections, update the construction parameter information to obtain updated construction parameter information, and mark the Cartesian mesh based on the updated construction parameter information and the geometric file information to obtain a new mesh; wherein, the Cartesian mesh is the mesh used when performing wind field calculations using the lattice Boltzmann method; The historical grid adjustment module is used to, if the type of the grid generation instruction is to adjust the historical grid, select a target historical grid that matches the grid generation instruction from the historical grid, determine the grid adjustment parameters based on the grid generation instruction, and adjust the target historical grid based on the grid adjustment parameters to obtain the adjusted target historical grid. The configuration information of the Cartesian grid is determined, including the determination of fluid marker points and grid points. This includes: configuring attribute markings and equipment type markings for all grid points in the Cartesian grid to obtain configuration information; determining the watershed range, randomly determining fluid marker points within the watershed range, and determining grid points; where grid points are points in the Cartesian grid; and fluid marker points are randomly selected marker points within the watershed range that are located within the fluid watershed. The number of intersections of fluid markers is calculated, including: acquiring the geometric file model; initializing the device environment of the Turing architecture graphics processor to obtain the initialized Turing architecture graphics processor; using the initialized Turing architecture graphics processor to traverse the environment of the geometric file model; then dividing the geometric file model into bounding boxes; and using RTCore hardware to calculate the number of intersections between the line connecting each mesh and fluid marker and the geometric triangle face in the geometric file model. An intersection is defined as the point formed by the intersection of the line connecting the mesh point and the fluid marker with the geometric triangle face in the geometric file model. Based on the mesh generation command, mesh adjustment parameters are determined, and the target historical mesh is adjusted according to the mesh adjustment parameters, including: obtaining a mesh copy corresponding to the target historical mesh and opening the mesh copy in read-write mode; based on the mesh generation command, mesh adjustment parameters are determined, and the mesh copy is adjusted according to the mesh adjustment parameters; The grid adjustment parameters are determined based on the grid generation command, and the grid copy is adjusted according to the grid adjustment parameters, including: determining the grid adjustment parameters based on the grid generation command; wherein, the grid adjustment parameters include the number of devices to be adjusted, the Cartesian grid point number, and the attribute marker; adjusting the grid copy according to the grid adjustment parameters to obtain the adjusted grid copy, and saving the grid copy to obtain the adjusted target historical grid.

5. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the Cartesian grid generation method as described in any one of claims 1 to 3.

6. A computer-readable storage medium, characterized in that, Used to store a computer program; wherein, when the computer program is executed by a processor, it implements the Cartesian grid generation method as described in any one of claims 1 to 3.

Images