CN116629162B - Unsteady flow field data processing method and device, electronic equipment and storage medium - Google Patents
Unsteady flow field data processing method and device, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN116629162B CN116629162B CN202310886830.8A CN202310886830A CN116629162B CN 116629162 B CN116629162 B CN 116629162B CN 202310886830 A CN202310886830 A CN 202310886830A CN 116629162 B CN116629162 B CN 116629162B
- Authority
- CN
- China
- Prior art keywords
- flow field
- grid
- data
- unsteady flow
- field data
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 48
- 238000012800 visualization Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 31
- 230000003993 interaction Effects 0.000 claims abstract description 12
- 230000000007 visual effect Effects 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000009877 rendering Methods 0.000 description 14
- 238000005516 engineering process Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000002452 interceptive effect Effects 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- RKTYLMNFRDHKIL-UHFFFAOYSA-N copper;5,10,15,20-tetraphenylporphyrin-22,24-diide Chemical group [Cu+2].C1=CC(C(=C2C=CC([N-]2)=C(C=2C=CC=CC=2)C=2C=CC(N=2)=C(C=2C=CC=CC=2)C2=CC=C3[N-]2)C=2C=CC=CC=2)=NC1=C3C1=CC=CC=C1 RKTYLMNFRDHKIL-UHFFFAOYSA-N 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
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/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Processing Or Creating Images (AREA)
Abstract
The application discloses an unsteady flow field data processing method, an unsteady flow field data processing device, electronic equipment and a storage medium, which are applied to the technical field of scientific simulation. The method comprises the steps of loading target unsteady flow field data of required time steps into a memory according to user requirements; a surface mesh storage structure is constructed that includes a plurality of cell containers for storing each dimension of mesh information and a sequence of structures for accessing the corresponding cell containers. Based on the characteristics of grid cells with different dimensions, respectively splitting the flow field grids of the target unsteady flow field data according to the dimension types of the grid cells, and storing the obtained surface grid data of each dimension grid cell into a surface grid storage structure; finally, performing animation visualization interaction processing on all the surface flow field data of the non-fixed-length time steps read from the surface grid storage structure, so that the playing efficiency of the non-fixed-length flow field animation visualization can be effectively improved.
Description
Technical Field
The present application relates to the field of scientific simulation technology, and in particular, to a method and apparatus for processing unsteady flow field data, an electronic device, and a readable storage medium.
Background
The unsteady flow field data is taken as important data generated by scientific simulation, and is a research object for computational fluid dynamics and scientific visualization. Since most scientific phenomena have time dependence, scientific computation requires sampling of a spatial flow field over a period of time to generate spatial volume data that varies with time, i.e., variant data, also referred to as unsteady flow field data. Unsteady flow field data is a collection of points describing one or more scalar values in three-dimensional European space, which contains many variables and features, and simultaneously records time and space information, and compared with unsteady flow field data, the physical quantity distribution and grid position of the unsteady flow field data are changed along with time. In actual storage, a series of data files are typically named chronologically, with a collection of stationary flow fields being described in different time steps. In other words, the unsteady flow field data is expressed by a plurality of transient unsteady flow field data which are continuous in time, and the unsteady flow field data of each time step are processed and rendered and played in an animation mode, so that the characteristics and rules of the flow field changing along with time can be vividly and particularly reflected.
In the related art, rendering and drawing are generally carried out on a steady flow field of a single time step, and then the whole data set is continuously played in an animation mode according to a time sequence to analyze unsteady flow field data. The method has good universality and can be combined with the visual interaction technology to assist analysis, but a large amount of memory space is required to be consumed for large-scale and high-time-precision unsteady data, and the requirement on hardware resources is high; when the grid is complex in a single time step, the data processing and rendering are long, so that the playing delay is large, and the animation is not smooth.
In view of this, improving the playing efficiency of the animation visualization of the unsteady flow field is a technical problem that needs to be solved by those skilled in the art.
Disclosure of Invention
The application provides a method and a device for processing unsteady flow field data, electronic equipment and a readable storage medium, which effectively improve the playing efficiency of animation visualization of an unsteady flow field.
In order to solve the technical problems, the application provides the following technical scheme:
in one aspect, the present application provides a method for processing unsteady flow field data, including:
pre-constructing a surface grid storage structure; the surface grid storage structure comprises a plurality of unit containers for storing each dimension grid information and a structure sequence for accessing the corresponding unit containers;
Loading target unsteady flow field data of the required time step into a memory according to the requirement of a user;
based on the characteristics of grid cells with different dimensions, respectively splitting the flow field grids of the target unsteady flow field data according to the dimension types of the grid cells, and storing the obtained surface grid data of each dimension grid cell into the surface grid storage structure;
and carrying out animation visual interaction processing on the surface flow field data of all the non-fixed-length time steps read in the surface grid storage structure.
Optionally, the loading the target unsteady flow field data of the required time step into the memory according to the user requirement includes:
all unsteady flow field data to be processed are read into an auxiliary memory in advance;
determining animation reading data parameters according to the overview information of the unsteady flow field data to be processed;
calculating the number of bytes occupied by the target unsteady flow field data in the file corresponding to the single time step required by the user according to the vertex, the number of units and the variable type corresponding to the animation read data parameter;
and directly reading the target unsteady flow field data from the memory based on the byte number.
Optionally, the determining the animation read data parameter according to the overview information of the unsteady flow field data to be processed includes:
Reading and visualizing the full data of the first time step of the unsteady flow field data to be processed to obtain overview information of the unsteady flow field data to be processed;
creating a character string container for storing variable names of the animation read data parameters based on the overview information;
according to the user demand instruction, a corresponding value is given to the pre-created zone bit; the zone bit is used for identifying whether to extract the three-dimensional surface of the flow field.
Optionally, the calculating the number of bytes occupied by the target unsteady flow field data in the file corresponding to the single time step of the user demand through the vertex and the number of units of the flow field and the variable type corresponding to the animation read data parameter includes:
reading vertexes and units forming a flow field according to a preset sequence to construct a geometric topological structure of the flow field;
reading variable data corresponding to the flow field based on the animation reading data parameters, and storing the variable data to the vertex or the unit to which the variable data belongs;
and calculating the area of the file inside occupied by each variable data according to the number of top points, the number of units and the variable type corresponding to each variable data.
Optionally, the reading the target unsteady flow field data directly from the memory based on the byte number includes:
Determining the number of vertexes, the number of units and the number of bytes occupied by one variable of the flow field according to the overview information;
calculating the unit data offset according to the number of the vertexes or the number of units and the number of bytes occupied by one variable;
when the file pointer reaches the unselected variable region, the file pointer is repositioned based on the unit data offset to jump directly to the selected variable region.
Optionally, the splitting the flow field grid of the target unsteady flow field data according to the dimension type of the grid unit, and storing the obtained surface grid data of each dimension grid unit into the surface grid storage structure, including:
and for the zero-dimensional grid cells and the one-dimensional grid cells in the target unsteady flow field data, dividing each zero-dimensional grid cell and each one-dimensional grid cell into a plurality of independent vertex cells, inserting each independent vertex cell into a cell container with corresponding dimension in the surface grid storage structure, and copying corresponding vertex variables and cell variables.
Optionally, the splitting the flow field grid of the target unsteady flow field data according to the dimension type of the grid unit, and storing the obtained surface grid data of each dimension grid unit into the surface grid storage structure, including:
For the two-dimensional grid cells in the target unsteady flow field data, if the current grid cells are linear cells, obtaining a vertex list of the current grid cells in the target unsteady flow field data, and inserting the vertex list into a cell container with corresponding dimensionality in the surface grid storage structure according to a preset sequence;
if the current grid unit is a nonlinear unit, replacing the current grid unit with a plurality of target linear units, and inserting a vertex list of the intersection point of each target linear unit into a unit container of corresponding dimension in the surface grid storage structure according to a preset sequence; and each target linear unit is spliced with the current grid unit.
Optionally, the splitting the flow field grid of the target unsteady flow field data according to the dimension type of the grid unit, and storing the obtained surface grid data of each dimension grid unit into the surface grid storage structure, including:
for linear three-dimensional grid cells in the target unsteady flow field data, inserting each face of the current grid cell into a cell list of a corresponding shape created in advance based on the type of the grid shape;
and after all the linear three-dimensional grid cells in the target unsteady flow field data are processed, inserting the surface data which do not belong to the inner face of the flow field in each cell list into a two-dimensional cell container, and copying corresponding strain values.
Optionally, the building the surface mesh storage structure includes:
generating three cell containers based on the grid cell dimensions, the total number of grid cells, and the total number of vertices contained in each grid cell;
and generating a corresponding structure sequence for each cell container according to the dimension corresponding to each grid cell in the current cell container and the cell data offset corresponding to each grid cell in the array.
Another aspect of the present application provides an unsteady flow field data processing apparatus, including:
the flow field network construction module is used for constructing a surface grid storage structure in advance; the surface grid storage structure comprises a plurality of unit containers for storing each dimension grid information and a structure sequence for accessing the corresponding unit containers;
the data loading module is used for loading the target unsteady flow field data of the required time step into the memory according to the user requirement;
the grid analysis and extraction module is used for respectively splitting the flow field grids of the target unsteady flow field data according to the dimension types of the grid cells based on the characteristics of the grid cells with different dimensions, and storing the obtained surface grid data of the grid cells with different dimensions into the surface grid storage structure;
And the visualization module is used for carrying out animation visualization interaction processing on the surface flow field data of all the non-fixed-length time steps read in the surface grid storage structure.
The application also provides an electronic device comprising a processor for implementing the steps of the unsteady flow field data processing method according to any one of the preceding claims when executing a computer program stored in a memory.
The application finally provides a readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the unsteady flow field data processing method of any one of the preceding claims.
The technical scheme provided by the application has the advantages that the flow field grids are simplified by constructing the surface grid storage structure, the surface grid storage structure only needs to store the surface grid data of the unsteady flow field data, and the redundant internal grid units of the unsteady flow field data and the vertex data thereof are deleted, so that after single time step data are loaded into the memory, the surface grid is analyzed and extracted as visual input, the data for visual processing do not contain redundant geometric and attribute data, the grid quantity is effectively reduced, the processing and rendering time of the single time step of the unsteady data is reduced, the accuracy and fluency of the visualization are ensured, the problem of the inter-frame drawing delay caused by the complex grid structure of the single time step is solved, and the playing efficiency and the interactive frame rate of the unsteady flow field animation visualization are effectively improved; memory resources are also effectively saved, and memory optimization management is realized.
In addition, the application also provides a corresponding implementation device, electronic equipment and a readable storage medium for the unsteady flow field data processing method, so that the method has more practicability, and the device, the electronic equipment and the readable storage medium have corresponding advantages.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application as claimed.
Drawings
For a clearer description of the present application or of the technical solutions related thereto, the following brief description will be given of the drawings used in the description of the embodiments or of the related art, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained from these drawings without the inventive effort of a person skilled in the art.
FIG. 1 is a schematic flow chart of an unsteady flow field data processing method provided by the application;
FIG. 2 is a schematic diagram of a surface mesh storage structure according to the present application;
FIG. 3 is a schematic flow chart of another unsteady flow field data processing method provided by the application;
FIG. 4 is a schematic representation of a nonlinear two-dimensional grid cell in an illustrative example provided by the present application;
FIG. 5 is a block diagram of an embodiment of an unsteady flow field data processing apparatus provided by the present application;
fig. 6 is a block diagram of an embodiment of an electronic device according to the present application.
Detailed Description
In order to better understand the aspects of the present application, the present application will be described in further detail with reference to the accompanying drawings and detailed description. It will be apparent that the described embodiments are only some, but not all, embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
The terms first, second, third, fourth and the like in the description and in the claims and in the above drawings are used for distinguishing between different objects and not necessarily for describing a particular sequential or chronological order. Furthermore, the terms "comprise" and "have," as well as any variations of the two, are intended to cover a non-exclusive inclusion. For example, a process, method, system, article, or apparatus that comprises a list of steps or elements is not limited to only those steps or elements but may include other steps or elements not expressly listed. Various non-limiting embodiments of the present application are described in detail below.
Referring first to fig. 1, fig. 1 is a schematic flow chart of a unsteady flow field data processing method provided by the present application, where the present application may include the following matters:
s101: the surface mesh storage structure is pre-built.
In this embodiment, the surface grid storage structure is a simplified flow field grid, and the flow field grid is a space structure for describing a flow field, and the grid cells can be divided into a zero-dimensional grid cell, a one-dimensional grid cell, a two-dimensional grid cell and a three-dimensional grid cell according to dimensions. Although the surface rendering visualization method only needs to use flow field surface data, complete grid information is still needed to acquire the three-dimensional surface of the flow field. In order to improve the visualization efficiency, the application only stores the unsteady time step data required in the animation visualization process by utilizing the surface grid storage structure, deletes the redundant internal grid units and the geometric data attribute data such as the vertexes thereof in the original data, reduces the grid quantity, and ensures that the animation visualization is applicable to larger unsteady flow field data. The surface mesh storage structure of the present embodiment includes a cell container (may be referred to as CellArray) and a junction uniquely corresponding theretoThe cell container is used for storing the grid information of each dimension, namely, the cell container is used for storing a plurality of grid cells in the original data after being split. The grid cells are arranged in the memory by arrays In the form of a storage, wherein,cellthe units are represented by the terms of a unit,cell 1 for the identification number of the first unit,cell n is the firstnThe identification number of the individual cell(s),pthe vertex is represented as such,nas a total number of units,m n represent the firstnIndividual unitscell n The number of the vertices included in the interior,m、nare all positive integers. />The identification number of the first vertex of the first unit,/for the first unit>Identification number for the second vertex of the first unit,/for the second vertex of the first unit>The first unit is the firstmIdentification number of each vertex->Is the firstnIdentification number of the first vertex of the individual element, for example>Is the firstnThe identification number of the second vertex of the cell,is the firstnFirst of unitsmThe identification number of the vertex. Since the number of top points of constituent cells is uncertain, a certain cell cannot be randomly accessed. In order to increase the access speed, a structure sequence of the auxiliary access unit container, that is to say a structure sequence for accessing the corresponding unit container, is simultaneously formed at the read-in unit.
S102: and loading the target unsteady flow field data of the required time step into the memory according to the user requirement.
In order to process large-scale unsteady flow field data in a mode of completely reading in the data, when the unsteady flow field data needing animation visualization processing, namely unsteady flow field data to be processed, is read, loading is carried out according to the requirement, wherein the target unsteady flow field data is the unsteady flow field data to be loaded to a memory at present, and the target unsteady flow field data is only a small part of the original unsteady flow field data to be processed.
S103: based on the characteristics of grid cells with different dimensions, respectively splitting the flow field grids of the target unsteady flow field data according to the dimension types of the grid cells, and storing the obtained surface grid data of the grid cells with different dimensions into a surface grid storage structure.
After the required time step data is loaded into the memory in the previous step, the surface mesh is analyzed and extracted as an animated visual input. It can be understood that grid cells with different dimensions have different characteristics, when analyzing the target unsteady flow field data and extracting the surface data, the splitting modes of the grid cells with different dimensions are different, the splitting purpose is to extract the surface grid data, namely, the finally stored data only keep the grid data necessary for animation drawing, at the moment, all flow field data necessary for animation visualization are loaded in a memory, redundant geometric and attribute data are removed, the accuracy and fluency of the visualization are ensured, the memory space necessary for storing the unsteady data is greatly reduced, memory resources are saved, and the optimized storage is performed by adopting the mode of extracting the surface grid on the memory management. When all grid cells of the unsteady flow field data to be processed are processed according to the steps, the original flow field data, namely the unsteady flow field data to be processed and auxiliary storage data are released, and all unsteady time steps are read and processed according to user parameters from user demands to obtain a set of surface flow field data. Furthermore, the extraction process of the grid on the surface of the flow field is relatively independent, so that the task can be divided to accelerate in a multithreading mode for improving the visual efficiency of the animation, in other words, the zero-dimensional grid unit, the one-dimensional grid unit, the two-dimensional grid unit and the three-dimensional grid unit can be processed in parallel.
S104: and carrying out animation visual interaction processing on all the surface flow field data of the non-fixed-length time steps read in the surface grid storage structure.
After the data to be processed is intensively read and processed in the steps, animation visualization is performed through the steps, so that the play delay between frames can be remarkably reduced, and the real-time interaction experience of users is improved. By way of example, the surface data corresponding to each time step can be mapped into the graphic primitive and rendered for display in sequence, so that the unsteady flow field animation is formed. The graphic objects calculated from the data of each time step, such as the isosurface, the characteristic region and the like, are stored as a plurality of data sets, participate in playing of the animation, and can adjust the properties of the visibility, the transparency and the like of the drawn image through the graphic objects. The user can adjust the playing speed by modifying the timer interval through interface interaction, and can also switch the timer state so as to pause or continue playing. After the pause, the current time step data can be analyzed by using a visualization technology of a steady flow field. In addition, under the support of geometric and topological data, the display window can apply camera interaction in a three-dimensional space, such as operations of field zooming, rotation, translation and the like when playing animation, so that flow field information is explored more comprehensively.
In the technical scheme provided by the application, the flow field grids are simplified by constructing the surface grid storage structure, the surface grid storage structure only needs to store the surface grid data of the unsteady flow field data, and the redundant internal grid units of the unsteady flow field data and the vertex data thereof are deleted, so that after single time step data are loaded into the memory, the surface grids are analyzed and extracted as visual input, the data for visual processing do not contain redundant geometric and attribute data, the grid quantity is effectively reduced, the processing and rendering time consumption of the single time step of the unsteady data is reduced, the accuracy and fluency of the visualization are ensured, the problem of the inter-frame drawing delay caused by the complex grid structure of the single time step is solved, and the playing efficiency and the interactive frame rate of the unsteady flow field animation visualization are effectively improved; memory resources are also effectively saved, and memory optimization management is realized.
The above embodiment does not limit how to construct the surface mesh storage structure, and this embodiment also provides an optional construction manner, which may include the following:
generating three cell containers based on the grid cell dimensions, the total number of grid cells, and the total number of vertices contained in each grid cell;
And generating a corresponding structure sequence for each cell container according to the dimension corresponding to each grid cell in the current cell container and the cell data offset corresponding to each grid cell in the array.
In this embodiment, the structure of the cell container and the corresponding structure sequence is shown in fig. 2, where the cell container includes a zero-dimensional cell container storing zero-dimensional grid information, a one-dimensional cell container storing one-dimensional grid information, and a two-dimensional cell container storing two-dimensional grid information (including two-dimensional grid cells and two-dimensional grid cells converted by three-dimensional grid cells), and the total number of grid cells refers to the total number of grid cells included in the cell container. The CellTypes order table is used to record the grid cell type, i.e., the grid cell dimension, and the offset of the cell within the array, i.e., the cell data offset.
In the embodiment, three groups of grid unit data structures can be used for storing zero-dimensional, one-dimensional and two-dimensional grid information after grid simplification processing, so that the grid quantity is effectively reduced, and the memory optimization management is realized.
It can be understood that the time dimension needs to be additionally considered when the unsteady flow field data is visualized, early research is constrained by hardware performance in the visual analysis process of the unsteady flow field data, and under the limit of limited computing resources, the flow field data cannot be loaded at one time, but a mode of reading in the data and rendering is used, so that the real-time interaction and rendering frame rate of the visualization are limited, and the full visual analysis of the large-scale unsteady flow field data is difficult. With the development of flow field visualization technology and the updating of computer hardware, the mode of firstly reading all data into processing and then rendering and playing becomes a feasible choice, and the most time-consuming step of reading data is finished in advance, so that the method can effectively analyze the dynamic and evolution of the phenomenon hidden in the data for the data with moderate scale, and the unsteady animation visualization can obviously reduce the delay between animation frames and improve the interactive response speed. However, for large-scale unsteady data with high time and spatial resolution, it is still difficult to directly load all the data into the memory for visualization, and the playing efficiency of the animation visualization of the unsteady flow field is still low. Based on this, the present application also provides the following embodiments, which are used to solve the problem that in the unsteady flow field animation visualization, data cannot be completely loaded into the memory due to large data size, as shown in fig. 3, and may include the following contents:
All unsteady flow field data to be processed are read into an auxiliary memory in advance; determining animation reading data parameters according to overview information of unsteady flow field data to be processed; calculating the number of bytes occupied by the target unsteady flow field data corresponding to the single time step of the user demand in the file according to the vertex, the number of units and the animation reading data parameter corresponding variable types of the flow field; and directly reading the target unsteady flow field data from the memory based on the byte number.
It can be understood that the variable values and grid positions of the unsteady flow field data can change with time, but the variable types and the grid number are always the same, so that only flow field variables required by research can be read and surface flow field grids necessary for visualization can be extracted when data are intensively loaded according to the focus of a user. In order to enable the user to initially explore unsteady flow field data, research data of interest is found. In order to determine the position of the target unsteady flow field data in the storage, it is necessary to determine animation read data parameters, which can be obtained by primarily analyzing the data set information, and by interactively setting the animation read data parameters, i.e. variable names of required variables, such as speed, density, temperature, pressure, etc., as animation visualizations for providing data input.
As an optional implementation manner, this embodiment also provides a setting manner of the parameters of the animation read data, which may include: reading and visualizing the full data of the first time step of the unsteady flow field data to be processed to obtain overview information of the unsteady flow field data to be processed; creating a character string container for storing variable names of animation read data parameters based on the overview information; according to the user demand instruction, a corresponding value is given to the pre-created zone bit; the flag bit is used for identifying whether to extract the three-dimensional surface of the flow field.
In this embodiment, the first time step full amount of the unsteady flow field data input by the user can be read and visualized to provide overview information of the unsteady flow field. Under the support of complete flow field data, visualization technology including volume drawing, isosurface and cloud image maturation can be used for assisting user analysis, and other time steps can be switched to and the same operation can be performed. For example, a string container vector < string > may be created that holds the variable name of the desired load for skipping some of the redundant flow field variables during subsequent I/O (input/output) processing. And setting a zone bit, and selecting whether to extract the three-dimensional surface of the flow field by a user, thereby simplifying the flow field grid.
As an alternative implementation manner, the embodiment also provides an implementation manner of skipping redundant flow field variables, which may include determining the size of the data to be read in the file, and how to skip from the data to be read, that is, the unselected variable data, to the data to be read, that is, the selected variable data, which may include the following contents:
reading vertexes and units forming the flow field according to a preset sequence to construct a geometric topological structure of the flow field; reading variable data corresponding to the flow field based on the animation reading data parameters, and storing the variable data to the vertex or the unit to which the variable data belongs; and calculating the area of the file inside occupied by each variable data according to the number of top points, the number of units and the variable type corresponding to each variable data.
Determining the number of vertexes, the number of units and the number of bytes occupied by one variable of the flow field according to the overview information; calculating the unit data offset according to the number of vertexes or the number of units and the number of bytes occupied by one variable; when the file pointer reaches the unselected variable region, the file pointer is repositioned based on the unit data offset to jump directly to the selected variable region.
In this embodiment, the hierarchical structure of the storage system can know that the process of loading data from the auxiliary memory into the memory is very time-consuming, so a visualization scheme is used that all data is read in first, and then a single time step is directly obtained from the memory and rendered. In order to accelerate the data reading speed and reduce the memory occupation, the embodiment uses the animation obtained in the last step to read the data parameters such as variable names, and calculates the number of bytes occupied by the variable data in the file through the vertexes, the number of units and the variable types of the flow field when the data file is read: when loading flow field data, firstly, the vertexes and units forming the flow field are required to be read in sequentially, so that the geometric and topological structures of the flow field are constructed, each vertex corresponds to a unique vertex number, and each grid unit corresponds to a unique unit number. And then reading variable data of the flow field, storing the variable data on the vertex or the unit to which the variable belongs, wherein the header of the file variable data area contains variable description information such as variable data type, variable name and whether the variable is distributed on the vertex or the unit. Each variable occupies a continuous area in the data file, and the size of the area can be calculated through the number of top points, the number of units and the variable data type. When the file pointer reads an unwanted variable, a unit data offset can be added to the pointer offsetThus skipping this portion of content, the unit data offset may be calculated in the following manner:offset=n×size，nthe number of vertices or cells, because both vertices and cells may contain flow field attribute data,sizethe number of bytes occupied for a variable. Derived from the first time-step complete data read innAndsizefor subsequent data reading, the file pointer reaching the non-selected variable area and then being usedoffsetAnd repositioning, so that the optimization of the data reading algorithm is completed, and the performance is improved at the data I/O level.
In this embodiment, the variable offset is used to skip redundant variable data during data reading, thereby realizing selective reading. The flow field data on-demand loading technology is adopted, so that the memory overhead of the loaded data can be reduced while the visual result is kept correct, and the one-time total loading of a plurality of time steps of data of an unsteady flow field is realized. Of course, those skilled in the art may delete the redundant flow field variables after reading the complete data to achieve the purpose of on-demand loading without adopting the above-mentioned method of on-demand loading variable data.
The above embodiment does not limit how to analyze and extract the surface mesh data of the non-stationary flow field data, and based on this, the present application also provides an alternative implementation, which may include the following:
For a zero-dimensional grid unit and a one-dimensional grid unit in the target unsteady flow field data, dividing each zero-dimensional grid unit and each one-dimensional grid unit into a plurality of independent vertex units, inserting each independent vertex unit into a unit container with corresponding dimension in a surface grid storage structure, and copying corresponding vertex variables and unit variables.
For two-dimensional grid cells in the target unsteady flow field data, if the current grid cells are linear cells, obtaining a vertex list of the current grid cells in the target unsteady flow field data, and inserting the vertex list into cell containers with corresponding dimensions in a surface grid storage structure according to a preset sequence; if the current grid unit is a nonlinear unit, replacing the current grid unit with a plurality of target linear units, and inserting a vertex list of the intersection point of each target linear unit into a unit container of corresponding dimension in the surface grid storage structure according to a preset sequence; each target linear cell concatenates the current grid cell.
For linear three-dimensional grid cells in the target unsteady flow field data, inserting each face of the current grid cell into a cell list of a corresponding shape created in advance based on the type of the grid shape; after all the linear three-dimensional grid cells in the target unsteady flow field data are processed, the surface data which do not belong to the inner face of the flow field in each cell list are inserted into a two-dimensional cell container, and corresponding strain values are copied.
In this embodiment, in order to increase the data processing speed, the extraction of the original grid surface is performed by traversing the auxiliary access sequence table, and the grid cells are processed respectively according to different types, so that unnecessary indirect access operations are reduced. For zero-dimensional grid cells of the vertex and multi-vertex type, the multi-vertex cell can be split into a plurality of independent vertices, and then all the vertex cells are inserted into a zero-dimensional cell container of an output grid, namely a surface grid storage structure, and the corresponding vertices and cell variables are copied. Similarly, the line segment and polyline one-dimensional elements are split and inserted.
All two-dimensional grid cells in the original data participate in forming a flow field surface, but the two-dimensional grid cells have the difference of linearity and nonlinearity, and the nonlinearity grid cells need to be additionally divided into a plurality of linear grid cells so as to facilitate graphic mapping, so that different processing modes are also needed according to the types of the grid cells:
if the linear two-dimensional grid unit is: for example, for linear grid cells such as triangles and polygons, a vertex list of the grid cells in the original data can be obtained and inserted into a two-dimensional cell container of an output two-dimensional cell, namely, a surface grid storage structure according to a preset sequence. The preset sequence is the same as the storage sequence of the points inside the grid unit, for example, the points inside the grid unit need to be stored in a certain sequence to be normally used in subsequent calculation. These involved units are typically stored in a counter-clockwise order, and correspondingly, the vertex list of the acquisition unit in the raw data is inserted into the output two-dimensional unit in a counter-clockwise order.
If the non-linear two-dimensional grid unit is: the subdivision of the nonlinear two-dimensional cells is to accurately replace the original nonlinear grid cells with a plurality of linear two-dimensional grid cells, in other words, to subdivide the nonlinear grid cells into a plurality of linear two-dimensional grid cells, which can be spliced into nonlinear grid cells. For example, as shown in fig. 4, the quadratic tetragonal unit in the figure is a nonlinear two-dimensional grid unit, and may be divided into 4 quadrilaterals by introducing one central line on each of four sides, and the four quadrilaterals are linear grid units, and central line intersection points may be added to the vertex list. Similarly, a secondary triangle may be represented by four triangles.
For three-dimensional grid cells, only the two-dimensional faces that are flow field surfaces need to be preserved, i.e., faces that are contained by only one three-dimensional cell. Since the two-dimensional surfaces of the three-dimensional units are not stored separately, the normal extraction method searches for the vertices constituting the two-dimensional surfaces, each vertex obtains a list of all units including the vertex, and each list is intersected, and if the result is the current unit only, the two-dimensional surfaces can be confirmed as the surfaces to be reserved. The judging mode can cause larger data retrieval cost, and the embodiment can be used for classifying and processing according to the types of the three-dimensional units, so that the two-dimensional surface process of extracting the three-dimensional units is obviously accelerated.
For linear three-dimensional grid cells, a plurality of cell lists of shapes such as three cell lists of triangles, quadrilaterals and polygons can be created first, and each face of the linear three-dimensional cell is inserted into the corresponding list by type. The insertion rule may be: and inserting vertex numbers forming the faces according to a preset sequence, searching whether a current face exists in the currently inserted list, setting the grid cell number to-1 if the current face exists, skipping the face, and setting the grid cell number to the number where the grid cell in the original data exists if the current face does not exist. The preset sequence is the same as the storage sequence of the points inside the grid unit, for example, the points inside the grid unit need to be stored in a certain sequence to be normally used in subsequent calculation. The units concerned are usually stored in a counterclockwise order, and the vertex numbers constituting the faces are inserted in a counterclockwise order, accordingly. After all grid cells of the original data are processed, traversing the shape lists, wherein the surface with the cell number of-1 is the inner surface of the flow field, and the flow field is skipped directly. And the number of the non-1 is the effective surface, the effective surface is inserted into a two-dimensional cell container of the surface data, and variable values are copied from the original data according to the grid cell number.
The nonlinear three-dimensional grid unit can not accelerate operation by adopting the processing mode of the linear three-dimensional unit, and still needs to obtain a nonlinear surface unit by intersection according to a conventional method, then subdivide the nonlinear surface unit into a plurality of linear two-dimensional grid units, and insert the linear two-dimensional grid units into a two-dimensional unit container of output data, namely a surface grid storage structure. In the embodiment, when the grid on the surface of the flow field is extracted, an auxiliary storage structure is created for the three-dimensional grid unit, so that the processing speed is improved, and the visual playing efficiency of the unsteady flow field animation is further improved.
As can be seen from the above, the related art adopts loading a single time step data file from an external memory to an internal memory, and storing a geometric structure by using a vertex list according to the vertex and grid unit numbers forming the flow field in the storage structure in the internal memory, so as to store the spatial distribution of the flow field. And then reading the variables recorded on the vertexes and the grid units into the memory in the form of a plurality of arrays, so that the variables correspond to the serial numbers of the vertexes or the units where the variables are positioned, and obtaining the physical quantity distribution of the flow field. Thus, the animation drawing of each time step in the related art needs to completely load the steady flow field data of the step, and when the data size is large, the time cost of the data I/O can cause high play delay. Compared with the related art, the method adopts the flow field data on-demand loading technology in the data reading stage, and can reduce the memory overhead of loading data while keeping the visual result correct. The related art selects a data processing and rendering mode according to flow field information of interest in the visual study. Because of the high computational cost of volume rendering, surface rendering is generally used for animation playback smoothness, intermediate geometric primitives are constructed in a three-dimensional data field, and then the rendering of pictures is realized by conventional computer graphics technology. Because the data input required by the visualization is only a part of the original data, the prior art can cause the waste of memory resources when loading the complete original data, thereby reducing the unsteady data scale capable of carrying out animation visualization. The application adopts the surface grid storage structure to simplify the flow field grid, can effectively improve the processing and rendering efficiency of unsteady flow field data and improve the playing efficiency of unsteady flow field animation visualization.
It should be noted that, in the present application, there is no strict sequence of execution among the steps, so long as the sequence accords with the logic sequence, the steps may be executed simultaneously, or may be executed according to a certain preset sequence, and fig. 1 and fig. 2 are only a schematic manner, and do not represent that only such execution sequence is possible.
The application also provides a corresponding device for the unsteady flow field data processing method, so that the method has more practicability. Wherein the device may be described separately from the functional module and the hardware. In the following description, the unsteady flow field data processing apparatus provided by the present application is used to implement the unsteady flow field data processing method provided by the present application, in this embodiment, the unsteady flow field data processing apparatus may include or be divided into one or more program modules, where the one or more program modules are stored in a storage medium and executed by one or more processors, to complete the unsteady flow field data processing method disclosed in the first embodiment. Program modules in the present application are directed to a series of computer program instruction segments capable of performing particular functions, and are more suitable than programs themselves for describing the execution of an unsteady flow field data processing apparatus in a storage medium. The following description will specifically describe the functions of each program module in this embodiment, and the unsteady flow field data processing apparatus described below and the unsteady flow field data processing method described above may be referred to correspondingly.
Based on the angle of the functional module, referring to fig. 5, fig. 5 is a block diagram of an unsteady flow field data processing apparatus provided by the present application under an embodiment, where the apparatus may include:
a flow field network construction module 501, configured to pre-construct a surface grid storage structure; the surface grid storage structure comprises a plurality of unit containers for storing each dimension grid information and a structure sequence for accessing the corresponding unit containers;
the data loading module 502 is configured to load the target unsteady flow field data of the required time step into the memory according to the user requirement;
the grid analysis extraction module 503 is configured to split flow field grids of the target unsteady flow field data according to the dimension types of the grid cells based on characteristics of the grid cells with different dimensions, and store the obtained surface grid data of the grid cells with different dimensions into a surface grid storage structure;
and the visualization module 504 is used for performing animation visualization interaction processing on the surface flow field data of all the non-fixed-length time steps read in the surface grid storage structure.
Optionally, in some implementations of this embodiment, the data loading module 502 may be further configured to: all unsteady flow field data to be processed are read into an auxiliary memory in advance; determining animation reading data parameters according to overview information of unsteady flow field data to be processed; calculating the number of bytes occupied by the target unsteady flow field data corresponding to the single time step of the user demand in the file according to the vertex, the number of units and the animation reading data parameter corresponding variable types of the flow field; and directly reading the target unsteady flow field data from the memory based on the byte number.
As an alternative implementation of the foregoing embodiment, the foregoing data loading module 502 may be further configured to: reading and visualizing the full data of the first time step of the unsteady flow field data to be processed to obtain overview information of the unsteady flow field data to be processed; creating a character string container for storing variable names of animation read data parameters based on the overview information; according to the user demand instruction, a corresponding value is given to the pre-created zone bit; the flag bit is used for identifying whether to extract the three-dimensional surface of the flow field.
As another alternative to the above embodiment, the data loading module 502 may be further configured to: reading vertexes and units forming the flow field according to a preset sequence to construct a geometric topological structure of the flow field; reading variable data corresponding to the flow field based on the animation reading data parameters, and storing the variable data to the vertex or the unit to which the variable data belongs; and calculating the area of the file inside occupied by each variable data according to the number of top points, the number of units and the variable type corresponding to each variable data.
As yet another alternative implementation of the foregoing embodiment, the foregoing data loading module 502 may be further configured to: determining the number of vertexes, the number of units and the number of bytes occupied by one variable of the flow field according to the overview information; calculating the unit data offset according to the number of vertexes or the number of units and the number of bytes occupied by one variable; when the file pointer reaches the unselected variable region, the file pointer is repositioned based on the unit data offset to jump directly to the selected variable region.
Alternatively, in other implementations of this embodiment, the grid analysis extraction module 503 may include a first splitting unit, a second splitting unit, and a third splitting unit.
Wherein the first splitting unit is operable to: for a zero-dimensional grid unit and a one-dimensional grid unit in the target unsteady flow field data, dividing each zero-dimensional grid unit and each one-dimensional grid unit into a plurality of independent vertex units, inserting each independent vertex unit into a unit container with corresponding dimension in a surface grid storage structure, and copying corresponding vertex variables and unit variables.
Wherein the second splitting unit is operable to: for two-dimensional grid cells in the target unsteady flow field data, if the current grid cells are linear cells, obtaining a vertex list of the current grid cells in the target unsteady flow field data, and inserting the vertex list into cell containers with corresponding dimensions in a surface grid storage structure according to a preset sequence; if the current grid unit is a nonlinear unit, replacing the current grid unit with a plurality of target linear units, and inserting a vertex list of the intersection point of each target linear unit into a unit container of corresponding dimension in the surface grid storage structure according to a preset sequence; each target linear cell concatenates the current grid cell.
Wherein the third splitting unit is operable to: for linear three-dimensional grid cells in the target unsteady flow field data, inserting each face of the current grid cell into a cell list of a corresponding shape created in advance based on the type of the grid shape; after all the linear three-dimensional grid cells in the target unsteady flow field data are processed, the surface data which do not belong to the inner face of the flow field in each cell list are inserted into a two-dimensional cell container, and corresponding strain values are copied.
Optionally, in still other implementations of this embodiment, the flow field network building module 501 may be further configured to: generating three cell containers based on the grid cell dimensions, the total number of grid cells, and the total number of vertices contained in each grid cell; and generating a corresponding structure sequence for each cell container according to the dimension corresponding to each grid cell in the current cell container and the cell data offset corresponding to each grid cell in the array.
The functions of each functional module of the unsteady flow field data processing device of the present application may be specifically implemented according to the method in the above method embodiment, and the specific implementation process may refer to the related description of the above method embodiment, which is not repeated herein.
From the above, the present embodiment can effectively improve the playing efficiency of the animation visualization of the unsteady flow field.
The unsteady flow field data processing device mentioned above is described from the perspective of functional module, and further, the application also provides an electronic device, which is described from the perspective of hardware. Fig. 6 is a schematic structural diagram of an electronic device according to an embodiment of the present application. As shown in fig. 6, the electronic device comprises a memory 60 for storing a computer program; a processor 61 for implementing the steps of the unsteady flow field data processing method as mentioned in any of the embodiments above when executing a computer program.
Processor 61 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and processor 61 may also be a controller, microcontroller, microprocessor, or other data processing chip, among others. The processor 61 may be implemented in at least one hardware form of DSP (Digital Signal Processing ), FPGA (Field-Programmable Gate Array, field programmable gate array), PLA (Programmable Logic Array ). The processor 61 may also include a main processor, which is a processor for processing data in an awake state, also called a CPU (Central Processing Unit ), and a coprocessor; a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 61 may be integrated with a GPU (Graphics Processing Unit, image processor) for taking care of rendering and drawing of the content that the display screen is required to display. In some embodiments, the processor 61 may also include an AI (Artificial Intelligence ) processor for processing computing operations related to machine learning.
Memory 60 may include one or more computer-readable storage media, which may be non-transitory. Memory 60 may also include high-speed random access memory as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. The memory 60 may in some embodiments be an internal storage unit of the electronic device, such as a hard disk of a server. The memory 60 may in other embodiments also be an external storage device of the electronic device, such as a plug-in hard disk provided on a server, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash Card (Flash Card) or the like. Further, the memory 60 may also include both internal storage units and external storage devices of the electronic device. The memory 60 may be used to store not only application software installed in an electronic device, but also various types of data, such as: code or the like that executes a program during the unsteady flow field data processing method may also be used to temporarily store data that has been output or is to be output. In this embodiment, the memory 60 is at least used for storing a computer program 601, where the computer program can implement the relevant steps of the unsteady flow field data processing method disclosed in any of the foregoing embodiments after being loaded and executed by the processor 61. In addition, the resources stored in the memory 60 may further include an operating system 602, data 603, and the like, where the storage manner may be transient storage or permanent storage. The operating system 602 may include Windows, unix, linux, among other things. The data 603 may include, but is not limited to, data corresponding to the unsteady flow field data processing result, and the like.
In some embodiments, the electronic device may further include a display 62, an input/output interface 63, a communication interface 64, or referred to as a network interface, a power supply 65, and a communication bus 66. Among other things, the display 62, input output interface 63 such as a Keyboard (Keyboard) pertain to a user interface, which may optionally also include standard wired interfaces, wireless interfaces, etc. Alternatively, in some embodiments, the display may be an LED display, a liquid crystal display, a touch-sensitive liquid crystal display, an OLED (Organic Light-Emitting Diode) touch, or the like. The display may also be referred to as a display screen or display unit, as appropriate, for displaying information processed in the electronic device and for displaying a visual user interface. Communication interface 64 may optionally include a wired interface and/or a wireless interface, such as a WI-FI interface, a bluetooth interface, etc., typically used to establish a communication connection between an electronic device and other electronic devices. The communication bus 66 may be a peripheral component interconnect standard (peripheral component interconnect, PCI) bus, or an extended industry standard architecture (extended industry standard architecture, EISA) bus, among others. The bus may be classified as an address bus, a data bus, a control bus, etc. For ease of illustration, only one thick line is shown in fig. 6, but not only one bus or one type of bus.
Those skilled in the art will appreciate that the configuration shown in fig. 6 is not limiting of the electronic device and may include more or fewer components than shown, for example, may also include a sensor 67 that performs various functions.
The functions of each functional module of the electronic device according to the present application may be specifically implemented according to the method in the above method embodiment, and the specific implementation process may refer to the relevant description of the above method embodiment, which is not repeated herein.
From the above, the present embodiment can effectively improve the playing efficiency of the animation visualization of the unsteady flow field.
It will be appreciated that the unsteady flow field data processing method of the above embodiments may be stored in a computer readable storage medium if implemented in the form of software functional units and sold or used as a stand alone product. Based on such understanding, the technical solution of the present application may be embodied essentially or in part or all of the technical solution contributing to the related art, or may be embodied in the form of a software product stored in a storage medium, which performs all or part of the steps of the methods of the various embodiments of the present application. And the aforementioned storage medium includes: a U disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrically erasable programmable ROM, registers, a hard disk, a multimedia card, a card-type Memory (e.g., SD or DX Memory, etc.), a magnetic Memory, a removable disk, a CD-ROM, a magnetic disk, or an optical disk, etc., that can store program code.
Based on this, the present application also provides a readable storage medium storing a computer program which, when executed by a processor, performs the steps of the unsteady flow field data processing method according to any one of the embodiments above.
In this specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, so that the same or similar parts between the embodiments are referred to each other. For the hardware including the device and the electronic equipment disclosed in the embodiments, the description is relatively simple because the hardware includes the device and the electronic equipment corresponding to the method disclosed in the embodiments, and relevant places refer to the description of the method.
Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative elements and steps are described above generally in terms of functionality in order to clearly illustrate the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
The application provides a method and a device for processing unsteady flow field data, electronic equipment and a readable storage medium. The principles and embodiments of the present application have been described herein with reference to specific examples, the description of which is intended only to facilitate an understanding of the method of the present application and its core ideas. It should be noted that it will be apparent to those skilled in the art that various modifications and adaptations of the application can be made without departing from the principles of the application and these modifications and adaptations are intended to be within the scope of the application as defined in the following claims.
Claims (12)
1. A method for unsteady flow field data processing, comprising:
pre-constructing a surface grid storage structure; the surface grid storage structure comprises a plurality of unit containers for storing each dimension grid information and a structure sequence for accessing the corresponding unit containers;
loading target unsteady flow field data of the required time step into a memory according to the requirement of a user;
based on the characteristics of grid cells with different dimensions, respectively splitting the flow field grids of the target unsteady flow field data according to the dimension types of the grid cells, and storing the obtained surface grid data of each dimension grid cell into the surface grid storage structure;
And carrying out animation visual interaction processing on the surface flow field data of all the non-fixed-length time steps read in the surface grid storage structure.
2. The unsteady flow field data processing method according to claim 1, wherein loading the target unsteady flow field data of the required time step into the memory according to the user's requirement comprises:
all unsteady flow field data to be processed are read into an auxiliary memory in advance;
determining animation reading data parameters according to the overview information of the unsteady flow field data to be processed;
calculating the number of bytes occupied by the target unsteady flow field data in the file corresponding to the single time step required by the user according to the vertex, the number of units and the variable type corresponding to the animation read data parameter;
and directly reading the target unsteady flow field data from the memory based on the byte number.
3. The unsteady flow field data processing method according to claim 2, wherein determining animation read data parameters according to overview information of the unsteady flow field data to be processed includes:
reading and visualizing the full data of the first time step of the unsteady flow field data to be processed to obtain overview information of the unsteady flow field data to be processed;
Creating a character string container for storing variable names of the animation read data parameters based on the overview information;
according to the user demand instruction, a corresponding value is given to the pre-created zone bit; the zone bit is used for identifying whether to extract the three-dimensional surface of the flow field.
4. The unsteady flow field data processing method according to claim 2, wherein the calculating the number of bytes occupied by the target unsteady flow field data corresponding to the single time step of the user demand in the file through the vertex, the number of units and the animation read data parameter corresponding variable type of the flow field includes:
reading vertexes and units forming a flow field according to a preset sequence to construct a geometric topological structure of the flow field;
reading variable data corresponding to the flow field based on the animation reading data parameters, and storing the variable data to the vertex or the unit to which the variable data belongs;
and calculating the area of the file inside occupied by each variable data according to the number of top points, the number of units and the variable type corresponding to each variable data.
5. The unsteady flow field data processing method according to claim 2, wherein the reading the target unsteady flow field data directly from the memory based on the byte count includes:
Determining the number of vertexes, the number of units and the number of bytes occupied by one variable of the flow field according to the overview information;
calculating the unit data offset according to the number of the vertexes or the number of the units and the number of bytes occupied by one variable;
when the file pointer reaches the unselected variable region, the file pointer is repositioned based on the unit data offset to jump directly to the selected variable region.
6. The unsteady flow field data processing method according to claim 1, wherein the splitting the flow field grid of the target unsteady flow field data according to the dimension type of the grid unit, and storing the obtained surface grid data of each dimension grid unit in the surface grid storage structure includes:
and for the zero-dimensional grid cells and the one-dimensional grid cells in the target unsteady flow field data, dividing each zero-dimensional grid cell and each one-dimensional grid cell into a plurality of independent vertex cells, inserting each independent vertex cell into a cell container with corresponding dimension in the surface grid storage structure, and copying corresponding vertex variables and cell variables.
7. The unsteady flow field data processing method according to claim 1, wherein the splitting the flow field grid of the target unsteady flow field data according to the dimension type of the grid unit, and storing the obtained surface grid data of each dimension grid unit in the surface grid storage structure includes:
For the two-dimensional grid cells in the target unsteady flow field data, if the current grid cells are linear cells, obtaining a vertex list of the current grid cells in the target unsteady flow field data, and inserting the vertex list into a cell container with corresponding dimensionality in the surface grid storage structure according to a preset sequence;
if the current grid unit is a nonlinear unit, replacing the current grid unit with a plurality of target linear units, and inserting a vertex list of the intersection point of each target linear unit into a unit container of corresponding dimension in the surface grid storage structure according to a preset sequence; and each target linear unit is spliced with the current grid unit.
8. The unsteady flow field data processing method according to claim 1, wherein the splitting the flow field grid of the target unsteady flow field data according to the dimension type of the grid unit, and storing the obtained surface grid data of each dimension grid unit in the surface grid storage structure includes:
for linear three-dimensional grid cells in the target unsteady flow field data, inserting each face of the current grid cell into a cell list of a corresponding shape created in advance based on the type of the grid shape;
And after all the linear three-dimensional grid cells in the target unsteady flow field data are processed, inserting the surface data which do not belong to the inner face of the flow field in each cell list into a two-dimensional cell container, and copying corresponding strain values.
9. The unsteady flow field data processing method according to any one of claims 1 to 8, wherein the constructing a surface mesh storage structure includes:
generating three cell containers based on the grid cell dimensions, the total number of grid cells, and the total number of vertices contained in each grid cell;
and generating a corresponding structure sequence for each cell container according to the dimension corresponding to each grid cell in the current cell container and the cell data offset corresponding to each grid cell in the array.
10. An unsteady flow field data processing apparatus, comprising:
the flow field network construction module is used for constructing a surface grid storage structure in advance; the surface grid storage structure comprises a plurality of unit containers for storing each dimension grid information and a structure sequence for accessing the corresponding unit containers;
the data loading module is used for loading the target unsteady flow field data of the required time step into the memory according to the user requirement;
The grid analysis and extraction module is used for respectively splitting the flow field grids of the target unsteady flow field data according to the dimension types of the grid cells based on the characteristics of the grid cells with different dimensions, and storing the obtained surface grid data of the grid cells with different dimensions into the surface grid storage structure;
and the visualization module is used for carrying out animation visualization interaction processing on the surface flow field data of all the non-fixed-length time steps read in the surface grid storage structure.
11. An electronic device comprising a processor and a memory, the processor being configured to implement the steps of the unsteady flow field data processing method of any one of claims 1 to 9 when executing a computer program stored in the memory.
12. A readable storage medium, wherein a computer program is stored on the readable storage medium, which computer program, when being executed by a processor, implements the steps of the unsteady flow field data processing method of any one of claims 1 to 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310886830.8A CN116629162B (en) | 2023-07-19 | 2023-07-19 | Unsteady flow field data processing method and device, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310886830.8A CN116629162B (en) | 2023-07-19 | 2023-07-19 | Unsteady flow field data processing method and device, electronic equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116629162A CN116629162A (en) | 2023-08-22 |
| CN116629162B true CN116629162B (en) | 2023-10-20 |
Family
ID=87636870
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310886830.8A Active CN116629162B (en) | 2023-07-19 | 2023-07-19 | Unsteady flow field data processing method and device, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116629162B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226634A (en) * | 2013-04-19 | 2013-07-31 | 华南理工大学 | Computing method for unsteady flow field of rotary jet pump based on three-dimensional dynamic mesh |
| WO2017084106A1 (en) * | 2015-11-20 | 2017-05-26 | 田川 | System and method for numerical simulation of aircraft flow field |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10901118B2 (en) * | 2016-10-28 | 2021-01-26 | Exxonmobil Upstream Research Company | Method and system for enhancing meshes for a subsurface model |
-
2023
- 2023-07-19 CN CN202310886830.8A patent/CN116629162B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226634A (en) * | 2013-04-19 | 2013-07-31 | 华南理工大学 | Computing method for unsteady flow field of rotary jet pump based on three-dimensional dynamic mesh |
| WO2017084106A1 (en) * | 2015-11-20 | 2017-05-26 | 田川 | System and method for numerical simulation of aircraft flow field |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116629162A (en) | 2023-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Telea | Data visualization: principles and practice | |
| CN105261066B (en) | A kind of three-dimensional geographic information system real-time rendering multithreading distribution and control method | |
| US11436800B2 (en) | Interactive system and method providing real-time virtual reality visualization of simulation data | |
| JP5218109B2 (en) | Visualization data processing device, visualization data processing device control method, and visualization data processing device control program | |
| US20080246766A1 (en) | Numerical analysis mesh generating method and apparatus | |
| CN105517681A (en) | Chart conversion system using metadata and method therefor | |
| CN115798654A (en) | Model material processing method, device, equipment and medium | |
| CN114882188A (en) | Implementation method of small map based on digital twin | |
| US8339401B2 (en) | Curved surface importing device and method | |
| CN116629162B (en) | Unsteady flow field data processing method and device, electronic equipment and storage medium | |
| Stolte et al. | Parallel spatial enumeration of implicit surfaces using interval arithmetic for octree generation and its direct visualization | |
| She et al. | An efficient method for rendering linear symbols on 3D terrain using a shader language | |
| CN111402369A (en) | Interactive advertisement processing method and device, terminal equipment and storage medium | |
| Li et al. | Investigations on visualization and interaction of ship structure multidisciplinary finite element analysis data for virtual environment | |
| EP2715564B1 (en) | Mesh sequence file format | |
| Saquib | Mathematica Data Visualization | |
| CN114327395A (en) | Webpage design method, device, computer equipment, storage medium and product | |
| Ahrens et al. | A modular extensible visualization system architecture for culled prioritized data streaming | |
| CN116958385B (en) | Material texture dynamic updating method supporting mass monomer models, storage medium and equipment | |
| Griffith et al. | Feature Tracking in VR for Cumulus Cloud Life-Cycle Studies. | |
| CN114708140B (en) | Real-time interactive sectioning visualization method for numerical calculation model based on unit | |
| CN113627456B (en) | Method and device for generating reference data | |
| US11954802B2 (en) | Method and system for generating polygon meshes approximating surfaces using iteration for mesh vertex positions | |
| Li et al. | Mass models dynamic loading method research and implementation based on OSG | |
| US20230394767A1 (en) | Method and system for generating polygon meshes approximating surfaces using root-finding and iteration for mesh vertex positions |
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 |