JP2020191033A - Simulation result visualization system - Google Patents

Abstract

To provide a simulation result visualization system that allows a designer to easily determine the validity of a result from a three-dimensional fluid analysis simulation result.SOLUTION: A simulation result visualization system includes: a piping system design information holding unit that holds information on start coordinates and end coordinates of a pipe axis, pipe axial length, and a diameter of each piping component in a piping system consisting of a plurality of piping components; a fluid analysis simulation unit that calculates a physical quantity related to a fluid phenomenon at each time based on a physical formula that simulates the fluid phenomenon in the piping system; and a visualization unit that summarizes information on the physical quantity, draws it, and outputs drawn contents. The visualization unit includes: a physical statistic calculation unit that calculates physical statistics from the physical quantity; a pipe axis coordinate calculation unit that calculates pipe axis coordinates in the piping system; and a graph drawing unit that draws a physical statistic visualization graph with the horizontal axis taking the pipe axis coordinates and the vertical axis taking the physical statistics.SELECTED DRAWING: Figure 1

Description

The present invention relates to a simulation result visualization system related to piping design.

When designing the installation position of the flow meter in a piping system such as a nuclear power plant, if the flow meter cannot be installed at a position that satisfies the standards related to safety standards, the actual flow test or 3D as an alternative to the actual flow test It is necessary to perform a fluid analysis simulation and use the results to evaluate the effect of the flow field in the piping system on the accuracy of the flowmeter.

As a technique for evaluating the measurement accuracy of the flow meter in such a piping system, for example, Patent Document 1 describes the spool structure data of the spool of the ultrasonic flow meter and the piping flow path of the flow path through which the fluid to be measured flows. The actual machine condition storage means for storing the data, the fluid condition data of the fluid, the temperature condition data of the fluid, the spool structure data, the water supply pipe flow path data, the fluid condition data, and the temperature condition. A coupled analysis means of turbulence analysis / structural analysis for calculating the flow velocity distribution in the spool and the displacement of the spool under the conditions of the respective data at the time of flow measurement using the data, and the flow velocity distribution. An analysis result storage means for storing the displacement data of the spool, an ultrasonic propagation analysis means for performing ultrasonic propagation analysis of each measurement line of the ultrasonic flowmeter using the analysis result, and each measurement. From the ultrasonic propagation analysis result of the line, the deviation is added to at least one kind of data of the flow rate measurement accuracy analysis means for calculating the flow rate measurement accuracy, the spool structure data, the piping flow path data, the fluid condition data, and the temperature condition data. It is provided with a sensitivity analysis means for performing the ultrasonic propagation analysis and the flow rate analysis and performing a sensitivity analysis for each parameter, and a design / manufacturing support data output means for outputting the design / manufacturing support data based on the sensitivity analysis result. A design and manufacturing support system for a sonic flow meter is disclosed.

Further, Patent Document 2 is a method for verifying an ultrasonic water supply flow meter of a nuclear power plant, and sets conditions for setting conditions for a simulated test in which water at a flow rate equivalent to that of an actual machine is flowed using a flow path that simulates the conditions of an actual machine. The measurement accuracy is verified by the process and the mock test verification step for verifying the measurement accuracy of the ultrasonic flowmeter under the conditions of the mock test set in the condition setting step, and the flow path and the operation of the actual machine are simulated. The flow velocity distribution calculation step for calculating the flow velocity distribution by flow analysis under the conditions, the sensor position calculation step for calculating the sensor position under the actual operating conditions, and the ultrasonic propagation are calculated in consideration of the flow velocity distribution and the sensor position. The measurement accuracy is verified by the ultrasonic propagation calculation step and the analysis verification step of verifying the flow rate measurement accuracy from the ultrasonic propagation calculation result calculated by the ultrasonic propagation calculation step, and further by the mock test verification step. A verification method of an ultrasonic water flow meter including a double verification step of comparing the verification result of a mock test with the analysis verification result of the analysis verification step and determining the validity of the verification result is disclosed.

Japanese Unexamined Patent Publication No. 2012-013416 JP 2011-112533

By the way, in the above-mentioned conventional technique, while it is possible to confirm the measurement accuracy of the flow meter, another method is required to investigate the error factor of the flow meter from the entire piping system. One of the methods is, for example, to perform a fluid analysis simulation and use the result to determine an error factor. However, since the fluid analysis simulation is a three-dimensional analysis and the amount of information is very large, there is a problem that it is difficult for the designer to judge the validity of the result.

The present invention has been made in view of the above, and an object of the present invention is to provide a simulation result visualization system in which a person in charge of design can easily judge the validity of a result from a three-dimensional fluid analysis simulation result.

The present application includes a plurality of means for solving the above problems. For example, in a piping system composed of a plurality of piping components, the start coordinates, end coordinates, and axial length of each piping component of each piping component. In addition, a piping system design information holding section that holds caliber information, a fluid analysis simulation section that calculates the physical quantity related to the fluid phenomenon at each time based on a physical formula that simulates the fluid phenomenon in the piping system, and information on the physical quantity. In a simulation result visualization system having a visualization unit that summarizes and outputs the drawing contents, the visualization unit includes a physical statistic calculation unit that calculates a physical statistic from the physical quantity, and a pipe axis coordinate in the piping system. It is assumed that it has a tube axis coordinate calculation unit for calculating the above, and a graph drawing section for drawing a physical quantity visualization graph in which the tube axis coordinates are taken on the horizontal axis and the physical statistics are taken on the vertical axis.

According to the present invention, the person in charge of design can easily judge the validity of the result from the result of the three-dimensional fluid analysis simulation.

It is a figure which shows the hardware configuration of the simulation result visualization system in the 1st Embodiment. It is a functional block diagram which shows the processing function of the simulation result visualization system. It is a functional block diagram which shows the processing function of a piping system design part. It is a functional block diagram which shows the processing function of the 3D fluid analysis simulation part. It is a functional block diagram which shows the processing function of a piping route information management part. It is a functional block diagram which shows the processing function of the tube axis coordinate calculation part. It is a functional block diagram which shows the processing function of a physical statistics calculation part. It is a functional block diagram which shows the processing function of a graph output part. It is a functional block diagram which shows the processing function of the 3D data output part. It is a figure explaining the processing flow of the simulation visualization system. It is a flowchart which shows the processing content of the piping system design processing. It is a flowchart which shows the processing content of the piping system design processing. It is a flowchart which shows the processing content of the fluid analysis simulation execution processing. It is a flowchart which shows the processing content of the part information extraction / additional setting processing. It is a flowchart which shows the processing content of the section definition processing. It is a flowchart which shows the processing content of the tube axis coordinate calculation process among the physical statistics / tube axis coordinate calculation process. It is a flowchart which shows the processing content of the physical statistic calculation processing among the physical statistic / tube axis coordinate calculation processing. It is a flowchart which shows the processing content of the graph output processing among the graph 3D data output processing. Graph ・ It is a flowchart which shows the processing content of 3D data output processing among 3D data output processing. It is a figure which shows the data structure of various data of the recording apparatus in a simulation result visualization apparatus. It is a figure which shows the data structure of a component mesh information, a component list, and a piping composition list. It is a figure which shows the data structure of a part information. It is a figure which shows the data structure of the piping composition information. It is a figure which shows the data structure of a pipe route information and a pipe route list. It is a figure which shows the data structure of a simulation condition and a model list. It is a figure which shows the data structure of the simulation result. It is a figure which shows the data structure of the cross-sectional information in each part on the path. It is a figure which shows the data structure of the cross-section unit simulation result. It is a figure which shows the data structure of the cross-section physical statistic information. It is a figure which shows the screen display example on a display. It is a figure which shows the hardware configuration of the simulation result visualization system in the modification.

Hereinafter, embodiments of the present invention will be described with reference to the drawings.

<First Embodiment>
The first embodiment of the present invention will be described with reference to FIGS. 1 to 30.

FIG. 1 is a diagram showing a hardware configuration of a simulation result visualization system according to the present embodiment.

In FIG. 1, the simulation result visualization system includes a simulation result visualization device D100, a display output unit D105, and an operation input unit D106.

The simulation result visualization device D100 outputs a graph of the simulation result, 3D data of the piping system, and the like to the display output unit D105 according to the user operation input from the operation input unit D106.

The simulation result visualization device D100 is composed of a CPU (Central Processor Unit) D101, a RAM (Random Access Memory) D102, a recording device D103, and an I / F (Interface) D104. The CPUD 101 executes arithmetic processing according to a program for visualizing the simulation result. The RAM_D102 temporarily stores data in a program for visualizing the simulation result. The recording device D103 records a program, data, or the like for visualizing the simulation result, and corresponds to a hard disk drive or the like.

The recording device D103 stores a piping system design program, a 3D fluid analysis simulation program, and a visualization program as programs for visualizing the simulation results. In addition, various data D1 used when executing these programs are also stored.

In the piping system program, parts of the piping system are created according to the operation input from the user, the piping configuration is constructed by connecting the parts, and the piping route is constructed by specifying the route on the piping configuration, and the piping route information. It is stored in the recording device as. The 3D fluid analysis simulation program specifies piping route information, simulation conditions, and a simulation model according to operation input from the user, executes simulation calculation, and saves the simulation result in the recording device D103.

The visualization program creates an image including a graph of pipe axis coordinates and physical statistics and a 3D data shape of the piping system from the simulation result according to the operation input from the user, and outputs the image to the display output unit D105. The visualization program consists of a physical statistics calculation module, a graph output module, a pipe axis coordinate calculation module, a three-dimensional data output module, and a piping route information management module.

The physical statistic calculation module extracts simulation results on a cross section designated by the piping route information management unit and calculates physical statistics such as turning strength. The graph output module creates and generates a graph image from physical statistics and tube axis coordinates. The pipe axis coordinate calculation module calculates the pipe axis coordinates in the cross section specified by the piping route information management unit. The three-dimensional data output module generates a three-dimensional shape image of a pipe from a simulation result according to a pipe route information management unit. The piping route information management module defines the cross section of the piping from the piping route information.

The I / F_D104 receives data related to the user's operation input from the operation input unit, and transmits an image including a graph of pipe axis coordinates and physical statistics and a 3D data shape of the piping system to the display output unit D105.

The display output unit D105 is for outputting an image including a graph of tube axis coordinates and physical statistics and a 3D data shape of the piping system to the user, and refers to a display or a printer. The operation input unit D106 converts the operation input from the user into an electric signal and transmits it to the simulation result visualization device, and corresponds to a keyboard, a touch panel, a mouse, or the like.

FIG. 2 is a functional block diagram showing a processing function of the simulation result visualization system.

In FIG. 2, the simulation result visualization system is composed of a piping system design unit B1, a 3D fluid analysis simulation unit B2, and a visualization unit B3.

The piping system design unit B1 constructs a piping route according to the piping system program and the operation input from the user, and stores the piping route information in the recording device. The 3D fluid analysis simulation unit B2 executes the simulation calculation according to the 3D fluid analysis simulation program and the operation input from the user, and saves the simulation result in the recording device. The visualization unit B3 creates an image including a graph of pipe axis coordinates and physical statistics and a 3D data shape of the piping system from the simulation results according to the visualization program, and outputs the image to the display output unit. The visualization unit B3 is composed of a piping route information management unit B31, a physical statistic calculation unit B32, a pipe axis coordinate calculation unit B33, a graph output unit B34, and a three-dimensional data output unit B35.

The piping route information management unit B31 defines the cross section of the piping from the piping route information according to the piping route information management module. The physical statistic calculation unit B32 calculates physical statistics such as turning strength according to the physical statistic calculation module. The pipe axis coordinate calculation unit B33 calculates the pipe axis coordinates in the cross section designated by the pipe route information management unit according to the pipe axis coordinate calculation module. The graph output unit B34 creates and generates a graph image G1 from physical statistics and tube axis coordinates according to the graph output module. The three-dimensional data output unit B35 generates a three-dimensional shape image B35 of the pipe from the simulation result according to the pipe route information management unit according to the three-dimensional data output module.

FIG. 10 is a diagram illustrating a processing flow of the simulation visualization system.

In FIG. 10, the piping system design unit B1 executes the piping system design process F1, creates piping route information, and transmits it to the 3D fluid analysis simulation unit B2 and the visualization unit B3. The 3D fluid analysis simulation unit B2 performs the fluid analysis simulation execution process F2 using the piping route information, generates a simulation result, and transmits it to the visualization unit B3. On the other hand, when the visualization unit B3 receives the piping route information from the piping system design unit B1, the visualization unit B3 acquires the type of the component from the piping route information and specifies the display color of the component by performing the component information extraction / additional setting process F3. Next, the visualization unit B3 defines the cross section of the pipe from the pipe route information by performing the route cross section definition process F4. When the visualization unit B3 receives the simulation result from the 3D fluid analysis simulation unit B2, the visualization unit B3 calculates the pipe axis coordinates from the cross-sectional information on the pipe route information by executing the physical statistic / pipe axis coordinate calculation process F5, and cross-sections. Extract the above simulation results and calculate the physical statistics. Finally, the visualization unit B3 generates a graph image from the physical statistics and the pipe axis coordinates by performing graph / 3D data output processing, and generates a pipe three-dimensional shape image from the path cross-sectional information and the simulation result. The three-dimensional shape of the pipe is displayed on the display or output to the printer.

FIG. 3 is a functional block diagram showing a processing function of the piping system design unit.

In FIG. 3, the piping system design unit B31 is composed of a component creation unit B11, a component storage unit B12, a component selection / connection unit B13, a piping configuration storage unit B14, a piping route information generation unit B15, and a piping route information storage unit B16. ing. The component creation unit B11 creates component shapes such as straight pipes, curved pipes, branch parts, and valves according to operation input from the user. The component storage unit B12 stores the component shape created by the component creation unit B11 in the recording device as component mesh information. The component selection / connection unit B13 selects a plurality of component mesh information stored in the component storage unit B12 according to the operation input from the user, and connects the component mesh information to each other to generate piping configuration information. Is. The pipe configuration storage unit B14 stores the pipe configuration information generated by the component selection / connection unit in the storage device. The piping route information generation unit B15 selects one from a plurality of piping configuration information stored in the piping configuration storage unit B14, specifies the direction and route through which the fluid flows in the piping configuration by input from the user. Piping route information is generated by adding to the piping configuration information. The piping route information storage unit B16 stores the piping route information generated by the piping route information generation unit B15 in the recording device.

11 and 12 are flowcharts showing the processing contents of the piping system design process.

In FIGS. 11 and 12, in the piping system design process F1, the component creation unit B11 first creates component mesh information by executing the component creation process, and transmits the component mesh information to the component storage unit B12. The component storage unit B12 stores the component mesh information in the recording device by executing the component information storage process, and registers the component mesh ID and the component name given to the component mesh information in the component list.

The component selection / connection unit B13 receives the component list from the component storage unit B12 by transmitting the component list request signal. By executing the component selection process, the component selection / connection unit B13 selects a component from the component list according to an operation input from the user, and acquires a component mesh ID of the component. The component selection / connection unit B13 receives the component mesh information from the component storage unit B12 by transmitting the component mesh information request signal together with the component mesh ID to the component storage unit B12. The component selection / connection unit B13 executes the component connection process after acquiring the mesh information of the two components, thereby rotating or moving each component according to the operation input from the user, or connecting the two components. Or The component selection / connection unit B13 generates piping configuration information from the result of connecting the components and transmits the piping configuration information to the piping configuration storage unit B14. The pipe configuration storage unit B14 saves the received pipe configuration information in the recording device, and registers the pipe configuration ID and the pipe configuration name given to the pipe configuration information in the pipe configuration list.

The pipe route information generation unit B15 receives the pipe configuration list from the pipe configuration storage unit B14 by transmitting the pipe configuration list request to the pipe configuration storage unit B14. By executing the piping configuration selection process, the piping route information generation unit B15 selects a piping configuration from the piping configuration list and acquires a piping configuration ID. The pipe route information generation unit B15 receives the pipe configuration information from the pipe configuration storage unit B14 by transmitting the pipe configuration information request including the pipe configuration ID to the pipe configuration storage unit B14. By executing the piping route designation process, the piping route information generation unit B15 specifies the direction and route in which the fluid flows in the piping configuration according to the operation input from the user, and adds the piping route information to the piping configuration information. Generate. The pipe route information storage unit B16 saves the pipe route information received from the pipe route information generation unit B15 in the recording device by executing the pipe route information storage process, and stores the pipe route ID in the pipe route information in the pipe route list. to register.

FIG. 4 is a functional block diagram showing the processing function of the 3D fluid analysis simulation unit.

In FIG. 4, the 3D fluid analysis simulation unit B2 is composed of a simulation condition setting unit B21, a simulation model selection unit B22, a piping route information reading unit B23, a simulation calculation unit B24, and a simulation result storage unit B25.

The simulation condition setting unit B21 sets calculation conditions for performing simulation calculation, such as initial values of flow velocity and pressure, calculation time interval, and integral calculation method. In this embodiment, all the setting items are put together in one simulation condition file and selected from a plurality of simulation condition files. A method of setting each setting item independently may be used. The simulation model selection unit B22 prepares a plurality of physical models of the fluid necessary for performing the simulation calculation, and selects from them. The piping route information reading unit 1_B23 reads the piping route information data created by the piping system design unit B1 from the recording device.

The simulation calculation unit B24 sets the piping route information according to the simulation conditions set by the simulation condition setting unit B21, the simulation model selected by the simulation model selection unit B22, and the piping route information read by the piping route information reading unit 1_B23. The simulation calculation of the fluid analysis in the piping system shown in is executed and the simulation result is generated. The simulation result storage unit B25 stores the simulation result generated by the simulation calculation unit B24 in the recording device.

FIG. 13 is a flowchart showing the processing contents of the fluid analysis simulation execution process.

In FIG. 13, in the fluid analysis simulation execution process F2, the pipe route information reading unit 1_B23 executes the pipe route information reading process to read the pipe route information data from the recording device and transmit it to the simulation execution unit B24. The simulation condition setting unit B21 sets the simulation condition and transmits it to the simulation execution unit B24 by executing the condition setting process independently of the pipe route information reading process. The model selection unit B22 selects a simulation model and transmits it to the simulation execution unit B24 by executing the model selection process independently of the piping route information reading process and the condition setting process. The simulation calculation unit B24 receives the simulation model, the simulation conditions, and the piping route information, and after all of them are prepared, performs the simulation calculation execution process to generate the simulation result and transmit it to the simulation result storage unit B25. The simulation result storage unit B25 receives the simulation result and stores it in the recording device by performing ten simulation result storage processes.

FIG. 5 is a functional block diagram showing a processing function of the piping route information management unit.

In FIG. 5, the piping route information management unit B31 includes the piping route information reading unit 2_B311, the route component information extraction / additional setting unit B312, the route component information storage unit B313, the route cross section definition unit B314, and the internal cross section of each component on the route. It is composed of an information storage unit B and 315. The piping route information reading unit 2_B311 reads the piping route information data created by the piping system design unit B1 from the recording device.

The in-route component information extraction / additional setting unit B312 extracts component information including the component type from the piping route information and sets the component color according to the operation input from the user. The in-route component information extraction / additional setting unit B312 is composed of a component information extraction unit B3121, a component type extraction unit B3122, and a component color setting unit B3123. The component information extraction unit B3121 extracts information about the component from the piping route information. The component type extraction unit B3122 extracts information on the type of component, such as a straight pipe or a curved pipe, from the component information. The part color setting unit B3123 manually or automatically sets the color of the part in the information about the part. The in-route component information storage unit B313 constructs in-route component information based on information about the component, the type of the component, and the color of the component, and stores the information in the route in the recording device. The route cross-section definition unit B314 defines the position of the pipe cross section for which the physical statistics are to be calculated in the pipe specified by the route component information generated by the route component information extraction / additional setting unit B312, and inside the route. It is added to the product information to form the cross-sectional information inside each part on the route. The cross-section information storage unit B315 in each part on the route stores the part information on the route configured by the route cross-section definition unit B314 in the recording device.

FIG. 14 is a flowchart showing the processing contents of the component information extraction / additional setting process.

In FIG. 14, in the component information extraction / additional setting process F3, the piping route information reading unit 2_B311 first executes the piping route information reading process to transmit the piping route information from the recording device to the component information extraction unit B3121. The component information extraction unit B3121 receives the piping route information and executes the component information extraction process to extract the component unit information in the route on the piping from the piping route information, and sets the component type extraction unit B3122 and the component color. Information on individual components in the route is transmitted to unit B3123. The component type extraction unit B3122 extracts the component type from the component unit information in the route, and the component color setting unit B3123 sets the display color of the component for the component unit information in the route. Regarding the display color of parts, the colors of adjacent parts must be different. For example, there is a method in which the hues of adjacent parts are always separated by 90 degrees or more. Specifically, when the piping route is composed of parts other than branches within 256 points, the color of the part number i is (R (i), G (i), B (i)), and R (i). ) Is represented by the following (Equation 1) to (Equation 4), G (i) is represented by the following (Equation 5) to (Equation 8), and B (i) is represented by the following (Equation 9) to (Equation 12). Can be done.

Further, the floor function floor (x, y) constituting each of the above equations is defined by the following (Equation 13), the ceiling function ceil (x, y) is represented by the following (Equation 14), and the remainder function mod (x, y) is defined. In the following (Equation 15), the function fz (x) is the following (Equation 16), the step size m of the color of the part is the following (Equation 17), and the maximum value of the number of colors of the part is the following (Equation 18). Each can be represented.

After that, the in-route part information storage unit B313 receives the part information / part type and the part color from the part type extraction unit B3122 and the part color setting unit B3123, respectively, and transmits the in-path part information by the part information / part type and the part color. Build and save to recording device.

FIG. 15 is a flowchart showing the processing contents of the cross-section definition processing.

In FIG. 15, in the cross-section definition process F3, the in-route component information extraction / additional setting unit B312 transmits the in-route component information to the route cross-section definition unit B315. The path cross-section definition unit B315 extracts information on the length of each part in the pipe axis direction. Next, the route cross-section definition unit B315 sets whether to divide the cross-section of the route equally or arbitrarily for the component by the operation input from the user. If it is divided into equal parts, the number of divisions is set by the operation input from the user, and the cross-sectional distance is calculated by dividing the length of the part in the pipe axis direction by (the number of divisions + 1). On the other hand, in any case, the cross-sectional distance is specified within the range of the length of the component in the pipe axis direction by the operation input from the user. This is repeated for all parts. Then, the route cross-section definition unit B315 embeds the cross-section distance in the part information in the route and constitutes the cross-section information in each part on the route. The cross-section information storage unit B316 of each part on the route receives the cross-section information of each part on the route from the route cross-section definition unit B315 and stores it in the recording device.

FIG. 6 is a functional block diagram showing a processing function of the tube axis coordinate calculation unit.

In FIG. 6, the pipe axis coordinate calculation unit B33 is composed of a section information reading unit 1_B331 in each part on the route, each section tube axis coordinate calculation unit B332, and a section information storage unit B333 in each part on the route with pipe axis coordinates. ing. The section 1_B331 for reading the cross-section information inside each component on the route reads the cross-section information inside each component on the route stored in the recording device. Each cross-section tube axis coordinate calculation unit B332 calculates the tube axis coordinates, the center coordinates, and the normal vector in each cross section by using the cross-section information in each part on the path. The section B333 for storing the cross-section information inside each part on the route with the pipe axis coordinates adds the pipe axis coordinates, the cross-section center coordinates, and the normal vector to the cross-section information inside each part on the route, and the cross-section information inside each part on the route with the pipe axis coordinates. It is summarized as, and stored in the recording device.

FIG. 16 is a flowchart showing the processing contents of the tube axis coordinate calculation process among the physical statistics / tube axis coordinate calculation processes.

In FIG. 16, in the pipe axis coordinate calculation process of the physical statistics / tube axis coordinate calculation process F5, the section information reading unit 1_B331 in each component on the route reads the cross section information in each component on the route stored in the recording device. It is transmitted to each cross-section tube axis coordinate calculation unit B332. Each cross-section tube axis coordinate calculation unit B332 first resets the tube axis coordinates to 0. Next, the start point coordinates and the end point coordinates are extracted for each component, and the tube axis coordinates LS (j) of the start point are calculated using the following (Equation 19).

Here, LP (j) is the length of the part number j in the pipe axis direction.

After that, by integrating the lengths in the pipe axis direction for each cross section, the pipe axis coordinates Lz (i, j) of the cross section number i of the part number j are calculated by the following (Equation 20).

Here, L (i, j) indicates the length in the pipe axial direction from the cross-section number i of the part number j to the cross-section number i + 1.

Next, the cross-section center coordinates (xO (i, j), yO (i, j), zO (i, j)) and the normal vector ez (j) = (ezx (j), in the cross-section number i of the part number j), Ezy (j) and ezz (j)) are calculated and obtained respectively. For example, when the part is a straight pipe and the cross section is created by equally dividing the part by the number of divisions imax, the center coordinates of the cross section are calculated by the following (Equation 21) and the normal vector is calculated by the following (Equation 22). be able to.

Here, (xs (j), ys (j), zs (j)) is the center coordinate of the component start point, and (xe (j), ye (j), ze (j)) is the center coordinate of the component end point. , Each shown. Each cross-section tube axis coordinate calculation unit B332 executes the above processing for all cross sections in all parts, and stores the tube axis coordinates, the cross-section center coordinates, and the normal vector in each part on the path with the tube axis coordinates. Send to B333. The cross-section information storage unit B333 of each part on the route with the pipe axis coordinates receives the cross-section information of each part on the route, the pipe axis coordinates, the cross-section center coordinates, and the normal vector, and the pipe axis coordinates are added to the cross-section information of each part on the route. And the cross-section center coordinates and the normal vector are added and saved in the recording device as the cross-section information in each part on the path with the tube axis coordinates.

FIG. 7 is a functional block diagram showing a processing function of the physical statistics calculation unit.

In FIG. 7, the physical statistics calculation unit B32 includes a simulation result reading unit 1_B321, a section information reading unit 2_B322 in each part on the route, a simulation result section extraction unit B323, a cell unit physical statistics calculation unit B324, and a cell unit physical statistics. It is composed of a calculation unit B325, a cross-sectional unit physical statistic integration calculation unit B326, a physical statistic time average / standard deviation calculation unit B326, and a cross-sectional unit physical statistic storage unit B327.

The simulation result reading unit 1_B321 reads the simulation result calculated by the simulation calculation execution unit B2 from the storage device. The route information reading unit 2_B322 reads the cross-section information of each component on the route from the storage device, which is generated by the piping route information management unit B31. The simulation result cross-section extraction unit B323 extracts the simulation result at the cross-section position defined by the cross-section information in each part on the path from the simulation results.

The cell-based physical statistic calculation unit B324 calculates the cell-based physical statistic from the simulation result on the cross section. Here, the cell is a polar coordinate system having the center coordinate of the cross section as the origin on the cross-sectional plane of the tube axis coordinate z and the elapsed time t in the simulation, (r, θ), (r + dr, θ), (r, θ + dθ). ), (R + dr, θ + dθ) refers to the region composed of four points. Note that (r, θ) refers to radial coordinates and rotational coordinates in the polar coordinate system. The cell-based physical statistic calculation unit B324 is composed of an angular momentum calculation unit, a momentum calculation unit, and other calculation units. The angular momentum calculation unit calculates the angular momentum of the fluid in cell units using the following (Equation 23). The momentum calculation unit calculates the momentum of the fluid in cell units using the following (Equation 24).

Here, ρ indicates the density of the fluid, Ur (r, θ, z, t) indicates the flow velocity in the r-axis direction, and Uθ (r, θ, z, t) indicates the flow velocity in the θ-axis direction. The other calculation unit calculates the parameters necessary for calculating other physical statistics.

The section unit physical statistic integration calculation unit B325 integrates the cell unit physical statistics and calculates the section unit physical statistics. The cross-section unit physical statistics integral calculation unit B325 is composed of each time turning strength calculation unit and other calculation units. Each time turning strength calculation unit calculates the turning strength m (z, t) at the tube axis coordinate z and the simulation elapsed time t by integrating the angular momentum and the momentum of the fluid in each cell by the following (Equation 25). ..

Here, Rmax indicates the radius of the pipe. The other calculation unit also calculates other physical statistics by integrating and calculating various parameters for each cell.

The physical statistic time average / standard deviation calculation unit B226 calculates the time average and standard deviation of the physical statistic obtained in time units. The physical statistic time average / standard deviation calculation unit B226 is composed of a turning intensity time average calculation unit, a turning strength standard deviation calculation unit, and other calculation units. The turning intensity time average calculation unit calculates the average value mmean (z) of the turning strength at the tube axis coordinate z from the turning strength at each time by using the following (Equation 26). Further, the turning strength standard deviation calculation unit calculates the standard deviation dm (z) of the turning strength from the turning strength at each time and its average value by using the following (Equation 27).

Here, Ts indicates the time average calculation start time in the simulation, and T indicates the elapsed time in the simulation. The other calculation unit calculates the time average and standard deviation of the physical statistics obtained in time units as well as the turning intensity.

The cross-section unit physical statistic storage unit B327 stores the calculated cross-section unit physical statistic in the recording device. The cross-section unit physical statistic storage unit B327 includes a turning strength storage unit and other storage units. The turning strength storage unit stores the time average and standard deviation of the turning strength in each cross section in the recording device. The other storage unit stores the time average and standard deviation of other physical statistics in each cross section in the recording device.

FIG. 17 is a flowchart showing the processing contents of the physical statistic calculation process among the physical statistic / tube axis coordinate calculation processes.

In FIG. 17, in the physical statistic calculation process of the physical statistic / tube axis coordinate calculation process F5, the simulation result reading unit 1_B321 reads the simulation result from the recording device and transmits it to the simulation result section extraction unit B323. On the other hand, at the same time, the cross-section information reading unit B322 of each part on the route reads the cross-section information of each part on the route from the recording device and transmits the simulation result cross-section extraction unit B323. The simulation result cross-section extraction unit B323 reads the simulation result and the cross-section information of each part on the route, and cross-sections on the pipe of the simulation result according to the pipe axis coordinates, the normal vector, and the cross-section center coordinates indicated by the cross-section information of each part on the route. Is created, the simulation result on the cross section is extracted, and is transmitted to the cell unit physical statistics calculation unit B324 as the cross section unit simulation result.

The cell unit physical statistic calculation unit B324 calculates the cell unit physical statistic from the simulation result on the cross section and transmits it to the cross section unit physical statistic integration calculation unit B325. The section unit physical statistic integration calculation unit B325 integrates the cell unit physical statistics, calculates the section unit physical statistics at each time, and transmits the physical statistics to the time average / standard deviation calculation unit B226. The physical statistic time average / standard deviation calculation unit B226 calculates the time average and standard deviation at each tube axis coordinate for the physical statistic obtained in time units, and transmits the time average and standard deviation to the section unit physical statistic storage unit B327. The cross-section unit physical statistic storage unit B327 stores the calculated cross-section unit physical statistic in the recording device.

FIG. 8 is a functional block diagram showing a processing function of the graph output unit.

In FIG. 8, the graph output unit B34 includes the cross-section information reading unit 1_B341 in each part on the path with pipe axis coordinates, the cross-section unit physical statistics reading unit B342, and the physical statistics / pipe axis coordinates (horizontal axis vertical axis) on each part cross section. ) It is composed of an extraction unit B343, a graph drawing unit B344, a graph setting unit B345, and a graph display output unit B346. The section 1_B341 for reading the cross-section information inside each component on the path with the tube axis coordinates reads the cross-section information inside each component on the path with the tube axis coordinates stored in the recording device. The cross-section unit physical statistic reading unit B342 reads the cross-section unit physical statistic stored in the recording device. Physical statistics on each part cross section / pipe axis coordinates (horizontal axis vertical axis) Extraction unit B343 is the pipe axis coordinates corresponding to each cross section from the cross section information in each part and the cross section unit physical statistics on the route with each pipe axis coordinates. And the time average and standard deviation of the part number and the physical statistic are extracted and linked by the section identifier. In addition, the part number, part type name, and part color are extracted from the cross-sectional information inside each part on the path with the coordinates of each pipe axis.

The graph drawing unit B344 draws a graph on the RAM with the horizontal axis as the tube axis coordinates and the vertical axis as the parameters related to the physical statistics. The graph drawing unit B344 is composed of a physical statistic time average drawing unit and a physical statistic time average ± standard deviation drawing unit. The physical statistic time average drawing unit draws a graph in which the horizontal axis is the tube axis coordinates and the vertical axis is the time average of the physical statistic. The physical statistics time average ± standard deviation drawing unit draws a graph in which the horizontal axis is the tube axis coordinates and the vertical axis is the value obtained by subtracting the standard deviation from the physical statistics time average and adding the value.

The graph setting unit B345 sets a color or adds a part type name to the graph according to the part number, the part type name, and the part color extracted from the cross-sectional information inside each part on the path with the coordinates of each pipe axis. The graph setting unit B345 includes a component type display unit and a component color setting unit. The component type display unit writes the name of the component in the area indicated by each component in the graph image data drawn on the RAM. The component color setting unit paints the graph portion of the region portion indicated by each component in the designated color in the image data of the graph drawn on the RAM. The graph display output unit B346 displays the image data of the graph drawn on the RAM on a display or prints it on paper with a printer.

FIG. 18 is a flowchart showing the processing contents of the graph output processing among the graph / 3D data output processing.

In FIG. 18, in the graph output process of the graph / 3D data output process F6, the cross-section information reading unit 1_B341 in each part on the path with the tube axis coordinates has the cross section inside each part on the path with the tube axis coordinates stored in the recording device. The information is read and transmitted to the physical statistics / tube axis coordinates (horizontal axis vertical axis) extraction unit B343 on each part cross section. At the same time, the cross-section unit physical statistic reading unit B342 reads the cross-section unit physical statistic stored in the recording device to the cross-section physical statistic / tube axis coordinate (horizontal axis vertical axis) extraction unit B343 of each part. Send. Physical statistics on each part cross section / pipe axis coordinates (horizontal axis vertical axis) Extractor B343 receives cross-section information inside each part and cross-section unit physical statistics on each route with pipe axis coordinates, and the pipe corresponding to each cross section. The axis coordinates, part numbers, time averages and standard deviations of physical statistics are extracted and linked by cross-section identifiers. In addition, the part number, part type name, and part color are extracted from the cross-sectional information inside each part on the path with the coordinates of each pipe axis.

Then, the physical statistics, the tube axis coordinates, the part number, the part type name, and the part color are transmitted to the graph drawing unit B344. The graph drawing unit B344 receives the physical statistics, the pipe axis coordinates, the part number, the part type name, and the part color, the horizontal axis is the pipe axis coordinates, and the vertical axis is the time average of the physical statistics and the standard deviation is subtracted from the time average. A graph obtained by adding the standard deviation to the time average and the time average is drawn on the RAM, and the graph image, the part number, the part type name, and the part color are transmitted to the graph setting unit. The graph setting unit B345 sets a color in the received graph image according to the received part number, part type name, and part color, adds the part type name, and transmits the decorated graph image to the graph display output unit B346. The graph display output unit B346 receives the decorated graph image, displays the image data of the graph drawn on the RAM on the display, or prints it on paper with a printer.

FIG. 9 is a functional block diagram showing a processing function of the 3D data output unit.

In FIG. 9, the 3D data output unit B35 includes a cross-section information reading unit 2_B351 in each component, a simulation result reading unit 2_B352, each component mesh extraction unit B353, a 3D data drawing unit B354, and a drawing setting unit B355 on the path with tube axis coordinates. It is composed of a 3D data display / output unit B356. The section 2_B351 for reading the cross-section information inside each component on the path with the coordinates of the tube axis reads the cross-section information inside each component on the path with the coordinates of the tube axis stored in the recording device. The simulation result reading unit 2_B352 reads the simulation result stored in the recording device. Each component mesh extraction unit B353 extracts the number, component color, and mesh data of each component from the cross-sectional information inside each component and the simulation result on the path with the coordinates of each pipe axis.

The 3D data drawing unit B354 draws an image of the three-dimensional shape of the pipe to be simulated and the mesh of the parts in the pipe on the RAM. The 3D data drawing unit B354 is composed of a 3D data outer frame drawing unit and each component mesh drawing unit. The 3D data outer frame drawing unit draws a three-dimensional image of the mesh data that becomes the outer frame of the pipe from the mesh data in the simulation result on the RAM. Each component mesh drawing unit draws a three-dimensional image of the mesh data of each component on the RAM.

The drawing setting unit B355 sets the color of the graph or sets the drawing angle according to the part color associated with the part number, and redraws the graph. The drawing setting unit B355 is composed of a component color setting unit and a drawing angle setting unit. The component color setting unit paints the mesh of each component in a designated color in the three-dimensional image of the pipe and the component drawn on the RAM. The drawing angle setting unit changes the angle of the viewpoint and redraws the three-dimensional image of the pipes and parts drawn on the RAM. The 3D data display / output unit B356 displays the three-dimensional image data of the pipes and parts drawn on the RAM on a display or prints them on paper with a printer.

FIG. 19 is a flowchart showing the processing contents of the 3D data output processing among the graph / 3D data output processing.

In FIG. 19, in the 3D data output process of the graph / 3D data output process F6, the cross-section information reading unit 2_B351 in each component on the path with tube axis coordinates is in each component on the path with tube axis coordinates stored in the recording device. The cross-section information is read and transmitted to each component mesh extraction unit B353. At the same time, the simulation result reading unit 2_B352 reads the simulation result stored in the recording device and transmits it to each component mesh extraction unit B353. Each component mesh extraction unit B353 receives the cross-sectional information inside each component and the simulation result on the path with the coordinates of each pipe axis, and extracts the mesh data of each component in the pipe. In addition, the part number, part type name, and part color are extracted from the cross-sectional information inside each part on the path with the coordinates of each pipe axis. Then, the simulation result, the part number, the color, and the mesh are transmitted to the 3D data drawing unit B354.

The 3D data drawing unit B354 receives the simulation result, the part number, the color, and the mesh, draws the 3D data which is the three-dimensional shape of the pipe and the part on the RAM, and draws the 3D data image, the part number, the part type name, and the part. Send the color to the graph setting section. The drawing setting unit B355 sets a color for the received 3D data image according to the received part number, part type name, and part color, sets the drawing angle, redraws, and displays the decorated 3D data in the 3D data display output unit B356. Send the image. The 3D data display output unit B356 receives the decorated 3D data image and displays the image data of the 3D data drawn on the RAM on a display or prints it on paper with a printer.

FIG. 20 is a diagram showing a data structure of various data of the recording device in the simulation result visualization device.

In FIG. 20, various data D1 includes data D11 for a piping system design program, a plurality of piping route information D12, a number of piping route information, a piping route list D13, data D14 for a 3D fluid analysis simulation program, and a plurality of simulation calculation time samples. It is composed of the simulation result D15 of the above, the data D16 for the visualization program, and other variables.

The piping system design program data D11 is a collection of constants and variables used by the piping system design unit B1 to execute processing according to the piping system design program. The data for the piping system design program is composed of a plurality of component mesh information D111, the number of component mesh information and the component list D112, a plurality of piping configuration information D113, the number of piping configuration information, the piping configuration list D114, and other temporary variables and constants. ing.

The part mesh information D111 is a collection of information such as points and cells constituting the parts in each part in the piping system. The part list D112 is a list in which a plurality of designed part mesh information D111s are registered, and is necessary for the user to display and select which part and which part are to be connected when connecting the parts. .. The pipe configuration information D113 is information describing a plurality of component mesh information D111 constituting the pipe and a connection status between the component meshes. The piping configuration list D114 is a list in which a plurality of designed piping configuration information D113s are registered, and is necessary for the user to display and select which piping configuration information D113 is to be set when setting the piping route. is there.

The piping route information D12 is information indicating which path and in which direction the fluid flows in the piping indicated by the piping configuration information D113 and the piping configuration information D113. The piping route list D13 is a list in which a plurality of designed piping route information D12s are registered, and is necessary for the user to display a list and select which piping route information D12 to use when executing the simulation. ..

The data D14 for the 3D fluid analysis simulation program is a set of constants and variables used by the 3D fluid analysis simulation unit B2 to execute the process according to the 3D fluid analysis simulation program. The data D14 for the 3D fluid analysis simulation program is composed of a simulation condition / model list D141, a simulation result time sample number, and other temporary variables. The simulation condition / model list D141 is a list containing a plurality of simulation conditions and simulation calculation formula models used when executing simulation calculation, and the user can select which simulation condition and simulation model to use by displaying the list. It is necessary to make it.

The simulation calculation time sample number is data indicating how much the step size of the elapsed time in the simulation should be when the simulation is executed. The simulation result D15 shows the value of the physical quantity at each position as a result of the simulation calculation using the simulation condition, the model, and the number of calculation time samples in the piping route information D12.

The visualization program data D16 is a collection of constants and variables used by the visualization unit B3 to execute processing according to the visualization program. The data D16 for the visualization program is composed of cross-sectional information D161 in each part on the path, cross-sectional unit simulation result D162, cross-sectional physical statistic information D163, graph image, piping 3D data image, and other temporary variables.

The cross-section information D161 in each part on the route is a collection of information on each part on the piping route, pipe axis coordinates, cross-section center coordinates, normal vector, and the like of each cross-section in the part. The cross-section unit simulation result D162 is an extraction of the simulation results for each cross-section on the piping path. The cross-section physical statistic information D163 is a physical statistic for each cross section on the piping path. The graph image is a graph image in which the horizontal axis is the tube axis coordinates and the vertical axis is the physical statistics. The pipe 3D data image is an image depicting the three-dimensional shape of the pipe and the parts inside the pipe.

FIG. 21 is a diagram showing a component mesh information, a component list, and a data configuration of a piping configuration list.

In FIG. 21, the component mesh information D111 is composed of a component mesh ID, component information D1111, a plurality of point information, a number of points in the mesh, a plurality of cell information, and a number of cells in the mesh. The part mesh ID is an identifier of the part mesh. The component information D1111 is information showing an outline of the component such as the type and name of the component, the start point and the end point. The point information shows information about the points that make up the mesh. The point information is composed of a point ID, a point X coordinate, a Y coordinate, a Z coordinate, and a local point ID.

The point ID indicates the ID of a plurality of points in the part mesh information. The local point ID indicates the ID of the point in the pipe route to be calculated by the simulation, and if the route in the pipe is not set, a value indicating not set is entered. The number of points in the mesh indicates the number of point information in the part mesh information. The cell information indicates which of the plurality of point information is used to form the cell that is the smallest unit of the three-dimensional shape. The cell information is composed of a cell ID, a plurality of point IDs, a number of cell constituent points, and a local cell ID.

The cell ID indicates an identifier of a plurality of cell information in the component mesh information. The plurality of point IDs describe the point IDs indicating specific point information. The number of cell constituent points indicates the number of point information constituting the cell. The local cell ID indicates the ID of the cell in the pipe route that is the target of the simulation calculation, and if the route in the pipe is not set, a value indicating that the route is not set is entered. The number of cells in the mesh indicates the number of cell information in the part mesh information.

The parts list D112 is composed of a plurality of parts list elements and the number of parts list elements. The parts list element shows an outline of parts that can be selected for creating a piping configuration, and is composed of a part mesh ID, a part type, and a part name. The part type indicates the major classification of parts such as straight pipes, curved pipes, and branches, and the part name indicates the minor classification of parts such as the model number and specifications of the parts. The pipe configuration list D114 is composed of a plurality of pipe list components and the number of pipe list components. The piping configuration list element shows an outline of a piping configuration that can be selected as a target for setting a piping route, and is composed of a piping configuration ID that is an identifier of the piping configuration and a name of the piping configuration.

FIG. 22 is a diagram showing a data structure of component information.

In FIG. 22, the component information D1111 includes a component type, a component name, a plurality of start point information, a plurality of start point information, a plurality of end point information, a plurality of end point information, a plurality of way point information, a number of way points, a representative point information, and a plurality of path element information. It consists of the total number of route elements, component route information, and the number of component routes.

The start point information shows the center coordinates of the cross section at the uppermost stream of the part. However, the definitions of upstream and downstream are tentatively added for convenience, and the upstream and downstream directions may be reversed as soon as the piping route is set. The start point information is composed of a start / end point ID, an X coordinate of the start point, a Y coordinate, and a Z coordinate. The start / end point ID indicates identifiers of a plurality of start points, end points, and waypoints. The reason why there are a plurality of starting points is to correspond to a branch component having two upstream points. The end point information shows the center coordinates of the cross section at the most downstream of the component, and is composed of the start / end point ID, the X coordinate, the Y coordinate, and the Z coordinate of the end point.

The waypoint information indicates the waypoint of the flow of parts, and is composed of a start / end point ID, an X coordinate of the end point, a Y coordinate, and a Z coordinate. The representative point information indicates a representative point of the component, and is composed of the X coordinate, the Y coordinate, and the Z coordinate of the representative point. In the case of a straight pipe, the representative point is often defined as the midpoint between the start point and the end point, and is basically a waypoint, but there is a representative point in the outer frame of the part like a curved pipe. Sometimes.

The path element information is information that divides the path of the fluid passing through the component for each element. The route element information is composed of a route element ID, a plurality of start / end IDs, a number of start / end points, a length in the pipe axis direction, an outer diameter, and an inner diameter. The route element ID indicates an identifier in a plurality of route element information. The length in the pipe axis direction indicates the length in the pipe axis direction of the path element. Inner diameter and outer diameter indicate the outer diameter and inner diameter of parts on the path. The route element is composed of a plurality of start points, end points, and waypoints. In the case of a straight pipe, there is only one combination of a start point and an end point, and in the case of a curved pipe, there is only one combination of a start point, an end point and a plurality of waypoints. In the case of branching, there are a plurality of route elements such as start point 1 and way point 1, start point 2 and way point 1, way point 1 and way point 2, way point 2 and end point 1, and way point 2 and end point 2.

The component element information is information indicating the entire path of the fluid passing from the specific start point to the specific end point in the component. The component route information is composed of the component route ID, the pipe axial length, the route element ID, and the number of route elements in the route. The component route ID indicates an identifier of a plurality of component routes. The length in the pipe axis direction indicates the length of the path from the specific start point to the specific end point of the component. The route element ID indicates the ID of the route element that constitutes the entire route. For example, in the case of the above-mentioned branch, the component route information includes the route elements such as the start point 1 and the way point 1, the way point 1 and the way point 2, and the way point 2 and the end point 1 as the part route information, so that the start point 1 or the way point 1 Alternatively, a component path such as a waypoint 2 or an end point 1 is formed. Since there are two starting points and two ending points, 2 × 2 = 4 component path information is formed.

FIG. 23 is a diagram showing a data configuration of piping configuration information.

In FIG. 23, the pipe configuration information D113 is composed of a pipe configuration ID, a pipe configuration name, a plurality of component mesh local information, a number of component mesh local information, a plurality of component connection information, and a number of component connections. The pipe configuration ID indicates a plurality of identifiers of pipe configuration information. The pipe configuration name indicates the character string of the name given to the pipe configuration by the user. The component mesh local information is obtained by adding information that is positioned in the piping configuration to the component mesh information D111.

The part mesh local information is composed of the part mesh information D111, the part mesh local ID, the part color, the representative point offset, and the normal vector. The component mesh local ID is an identifier of the component mesh local information in the piping configuration information. The part color indicates the display color of the part. The representative point offset indicates the position where the representative point of the component mesh is displayed in the piping configuration, and is composed of the X coordinate, the Y coordinate, and the Z coordinate. The normal vector indicates the angle of the component mesh, and is composed of an X-direction component, a Y-direction component, and a Z-direction component.

The component connection information shows the connection relationship between the component meshes. The component connection information is composed of a component connection information ID, a connection source component mesh local ID, a connection destination component mesh local ID, a connection source component mesh component path ID, a connection destination component mesh component path ID, and connection surface center coordinates. The component connection information ID is an identifier of a plurality of component connection information in the piping configuration information. The connection source component mesh local ID and the connection destination component mesh local ID indicate the local IDs of the connection source and connection destination component meshes. The connection source component mesh component path ID and the connection destination component mesh component path ID indicate the component path IDs of the connection source and the connection destination, and by specifying the component path ID, which component path is to be connected can be specified. Can be done. The connection surface center coordinates indicate the center coordinates on the connection surface between the parts, and are composed of the X coordinate, the Y coordinate, and the Z coordinate.

FIG. 24 is a diagram showing a data structure of piping route information and a piping route list.

In FIG. 24, the pipe route information D12 is composed of a pipe route ID, a pipe route name, a pipe configuration information, a plurality of pipe route single route information, and a number of pipe route single routes. The piping route ID indicates an identifier of a plurality of piping route information. The pipe route name indicates the character string of the name given to the pipe route by the user. The pipe configuration information D113 indicates a pipe configuration for which a route is specified. The piping route single route information specifies the piping route, and is composed of the piping route single route ID, a plurality of component connection direction information, and the number of component connection direction information. The pipe route single road ID is an identifier of a plurality of pipe route single road information. The component connection direction information specifies a flow direction for specific component connection information, and is composed of a component connection information ID and a flow direction. The flow direction is numerically identified as to whether the flow is from the temporary upstream to the temporary downstream or vice versa.

The piping route list D13 is composed of a plurality of piping route list elements and the number of piping route list elements. The piping route list element is a combination of a piping route ID and a piping route name for displaying and selecting a list.

FIG. 25 is a diagram showing a data structure of a simulation condition / model list.

In FIG. 25, the simulation condition / model list D141 is composed of a plurality of simulation conditions, a number of simulation conditions, a plurality of simulation models, a number of simulation models, a plurality of simulation condition constant information, and a plurality of simulation model constant information.

A simulation condition is a collection of a plurality of simulation condition constants, and is composed of a simulation condition ID, a simulation condition name, and a plurality of simulation condition constants. The simulation condition ID is an identifier in a plurality of simulation conditions. The simulation condition name is a character string in which the name of the simulation condition is described. The simulation condition constant is one of a group of physical constants representing a simulation condition, and is composed of an ID and a constant. The ID is an identifier of a plurality of simulation condition constants, and the value indicates the value of the physical constant indicated by the identifier. For example, when the ID indicates the initial value of the pipe pressure, a specific numerical value indicating the initial value of the pipe pressure is entered in the value.

A simulation model is a collection of a plurality of simulation model constants, and is composed of a simulation model ID, a simulation model name, and a plurality of simulation model constants. The simulation model ID is an identifier in a plurality of simulation models. The simulation model name is a character string in which the name of the simulation model is described. The simulation model constant is one of the coefficient groups of the mathematical formula representing the simulation model, and is composed of an ID and a constant. The ID is an identifier of a plurality of simulation model constants, and the value indicates the value of the coefficient of the mathematical formula indicated by the identifier. For example, when the ID indicates the influence coefficient of the pressure at the current time with respect to the pressure at the next time, a specific numerical value indicating the influence coefficient is entered in the value.

The simulation condition constant information indicates which physical constant is contained in which ID in the simulation condition constant, and is composed of an ID, a name, and a constant name. The ID is an identifier of a plurality of simulation condition constants, and the name is the name of the simulation condition constant. The constant name is the constant name of the simulation condition constant in the simulation calculation program. For example, if the ID indicates the initial value of the pipe pressure, the character string of "pipe pressure initial value" is entered in the name, and the constant name of the pipe pressure initial value used in the simulation calculation program is entered in the variable name. ..

The simulation model constant information indicates which ID contains the coefficient of which mathematical expression in the simulation model constant, and is composed of an ID, a name, and a constant name. The ID is an identifier of a plurality of simulation model constants, and the name is the name of the simulation model constant. The constant name is the constant name of the simulation model constant in the simulation calculation program. For example, if the ID indicates the coefficient of influence of the current time pressure on the pressure at the next time, the name will contain "the coefficient of influence of the current time pressure on the next time" and the variable name will be used in the simulation calculation program. Enter the constant name of the influence coefficient.

FIG. 26 is a diagram showing a data structure of the simulation result.

In FIG. 26, the simulation result D15 is composed of a result ID, a model ID, a condition ID, a piping route information D12, a plurality of point-based physical quantities, a total number of mesh points, a plurality of cell-based physical quantities, a total number of mesh cells, and time. There is. The result ID is an identifier of a plurality of simulation results D15. The model ID is the simulation model ID, and the condition ID is the simulation condition ID. The piping route information D12 refers to the piping route information D12 that is the target of the simulation calculation. Time is the elapsed time in the simulation results.

The point-based physical quantity indicates what value the physical quantity indicates at the coordinates of the points constituting the mesh on the pipe. The point unit physical quantity is composed of a local point ID, a plurality of physical quantities, and a number of physical quantity types. The local point ID is an ID for associating which point information among the piping route information D12 internal piping configuration information D113 internal component mesh information D111 internal point information corresponds to the point unit physical quantity. The physical quantity indicates the physical quantity at the coordinates indicated by the point information associated with the local point ID, and is composed of a number of dimensions and a plurality of physical quantity elements. For example, if it is a physical quantity indicating the scalar amount of the pipe pressure, the number of dimensions is 1, and the value of the scalar amount of the pipe pressure is entered in the physical quantity element. On the other hand, if the physical quantity indicates a three-dimensional vector value of the flow velocity, the number of dimensions is 3, and the X-axis direction component, the Y-axis direction component, and the Z-axis direction component of the flow velocity are contained in each of the three physical quantity elements.

The cell unit physical quantity indicates what value the physical quantity shows in the center coordinates of the cells constituting the mesh on the pipe. The cell unit physical quantity is composed of a local cell ID, a plurality of physical quantities, and a number of physical quantity types. The local cell ID is an ID for associating which cell information among the cell information in the piping route information D12 internal piping configuration information D113 internal component mesh information D111 corresponds to the point unit physical quantity. The physical quantity indicates the physical quantity in the representative coordinates of the cell indicated by the cell information associated with the local cell ID, and is composed of a number of dimensions and a plurality of physical quantity elements.

FIG. 27 is a diagram showing a data structure of cross-sectional information inside each component on the route.

In FIG. 27, the cross-sectional information D161 inside the parts on the route is composed of the piping route ID, the information on a plurality of individual parts in the route, and the number of parts in the route. The piping route ID is an identifier associated with the piping route ID described in the piping route information D12.

The individual component information in the route is information about the individual component in the route, which is necessary for drawing a graph of tube axis coordinates and physical statistics on the route. The single part information in the route includes the part mesh local ID, the front part mesh local ID, the rear part mesh local ID, the part type, the part color, the length in the pipe axis direction, the starting point pipe axis coordinates, the number of passing routes, and the number of passing routes. It is composed of route information, the number of cross sections, and cross section information for the number of cross sections. The component mesh local ID is a component mesh local ID associated with the piping route information D12 internal piping configuration information D113 internal component mesh information D111. The front part mesh local ID is the local ID of the part mesh directly connected to the upstream side of the corresponding part mesh. The rear part mesh local ID is the local ID of the part mesh directly connected to the downstream side of the corresponding part mesh.

The part type is information about the type of the part mesh, and indicates a straight pipe, a curved pipe, a branch, and the like. The part color is the color specified for the part and is used to identify the part in graphs and 3D data. The pipe axial length is the value of the length of the component in the pipe axial direction. The starting point tube axis coordinates are the values of the tube axis coordinates of the most upstream cross section of the component in the entire path. The passage route information is information indicating which pipe route single route the component passes through, and is composed of a pipe route single route ID and a flow direction number. The pipe route single road ID is an identifier indicating which of the pipe route single road information in the pipe route information D12 is the pipe route single road information. The flow direction number indicates which flow direction the single pipe route is in the pipe route.

The cross-section information is information indicating the attributes of each cross-section in the part. The section information is composed of a section ID, a section distance, a tube axis coordinate, a section center coordinate, and a normal vector. The cross-section ID is an identifier of a plurality of cross-section information. The cross-sectional distance indicates the distance from the previous cross section in the pipe axis direction. The pipe axis coordinates indicate the coordinates in the pipe axis direction in the entire piping route, and the actual value is entered in this item in the cross-section information inside each part on the route with the pipe axis coordinates, and the cross-section information inside each part on the route. Is a value indicating an exception in this item, for example, a physically impossible value of -1E + 10. The cross-section center coordinates indicate the center coordinates of the circle when the cross section is a circle, and are composed of the X coordinate, the Y coordinate, and the Z coordinate. The normal vector indicates in which direction the cross section is oriented in the coordinate system of the entire piping path, and is composed of an X direction component, a Y direction component, and a Z direction component. The in-path component information shown in FIG. 5 is the in-route component unit information without the cross-section information, and 0 (zero) is entered in the value of the number of cross-sections.

FIG. 28 is a diagram showing a data structure of a cross-section unit simulation result.

In FIG. 28, the cross-sectional unit simulation result D162 includes a result ID, a model ID, a condition ID, a time, a plurality of point information, a plurality of point unit physical quantities, a total number of all mesh points, a plurality of cell information, and a plurality of cell unit physical quantities + physical quantity. It consists of statistics and the total number of mesh cells.

The result ID is an identifier of a plurality of simulation results D15. The model ID is the simulation model ID, and the condition ID is the simulation condition ID. Time is the elapsed time in the simulation results.

The point information and cell information are regions in which the three-dimensional shape of the piping path composed of the piping route information D12 internal piping configuration information D11 internal component mesh is sliced into round slices in the cross section defined in each component internal cross section information D161 on the path. It is the information of points and cells in the mesh of the above, and the data structure is the same as that of the part mesh information D111. The point-based physical quantity is a mesh of a region in which the three-dimensional shape of the piping path composed of the piping route information D12 internal piping configuration information D11 internal component mesh is sliced in a cross section defined by each component internal cross-section information D161 on the path. The physical quantity at the points constituting the above is shown, and the data structure is the same as that of the simulation result D15.

The cell unit physical quantity + physical statistic is obtained by cutting the three-dimensional shape of the piping route composed of the piping route information D12 internal piping configuration information D11 internal component mesh into round slices with the cross section defined in each component internal cross-section information D161 on the route. It shows the physical quantities and physical statistics in the cells that make up the mesh of the area. Cell unit physical quantity + physical statistic is local cell ID, part mesh local ID, section ID, multiple physical quantities and physical quantity types, multiple cell unit physical statistics, multiple physical statistic integration values, and physical statistic types. It is composed of and. The local cell ID is an ID for associating which cell information among the cell information in the piping route information D12 internal piping configuration information D113 internal component mesh information D111 corresponds to the point unit physical quantity.

The component mesh local ID is a component mesh local ID associated with the piping route information D12 internal piping configuration information D113 internal component mesh information D111. The cross-section ID is an ID for associating which cross-section information of the cross-section information D161 inner cross-section information of each part corresponds to on the route. The physical quantity indicates the physical quantity at the representative coordinates of the cell indicated by the cell information associated with the local point ID, and is composed of a number of dimensions and a plurality of physical quantity elements. The cell-based physical statistic is the physical statistic in each cell before integration, for example, the angular momentum in the cell or the momentum in the cell. The physical statistic integral value indicates a physical statistic such as a turning intensity obtained as a result of integrating the angular momentum of each cell and the momentum of each cell in all cells.

FIG. 29 is a diagram showing a data structure of cross-sectional physical statistics information.

In FIG. 29, the cross-section physical statistic information D163 includes a result ID, a model ID, a condition ID, a pipe axis coordinate maximum value, a pipe axis coordinate minimum value, a cross section number, a total time, a pipe route single route number, and a cross section unit physical statistic. It consists of information.

The result ID is an identifier of a plurality of simulation results D15. The model ID is the simulation model ID, and the condition ID is the simulation condition ID. The total time is the elapsed time obtained by subtracting the start time from the end time of the simulation result. The pipe axis coordinate maximum and minimum values are the maximum and minimum values of the pipe section coordinates in the piping route. The number of single piping routes indicates the number of single piping routes in the piping route.

The cross-section unit physical statistic information is a collection of information necessary for drawing a graph of the physical statistic with respect to the tube axis coordinates in the cross-section unit. The cross-section unit physical statistic information includes the cross-section ID, the part mesh local ID, the physical statistic total information for the number of multiple physical statistic types, the part type, the part color, the tube axis coordinates, the plurality of passage route information, and the number of passage routes. It is composed of. The cross-section ID is an ID for associating which cross-section information of the cross-section information D161 inner cross-section information of each part corresponds to on the route. The component mesh local ID is a component mesh local ID associated with the piping route information D12 internal piping configuration information D113 internal component mesh information D111.

The physical statistic summary information is a summary of the time average and time standard deviation of each physical statistic over time. The part type is information about the type of the part mesh, and indicates a straight pipe, a curved pipe, a branch, and the like. The part color is the color specified for the part and is used to identify the part in graphs and 3D data. The pipe axial length is the value of the length of the component in the pipe axial direction. The tube axis coordinates are the values of the tube axis coordinates of the cross section in the entire path. The passage route information is information indicating which pipe route single route the component passes through, and is composed of a pipe route single route ID and a flow direction number. The pipe route single road ID is an identifier indicating which of the pipe route single road information in the pipe route information D12 is the pipe route single road information. The flow direction number indicates which flow direction the single pipe route is in the pipe route.

FIG. 30 is a diagram showing an example of screen display on a display.

In FIG. 30, the screen display on the display is composed of a pipe axis coordinate-physical statistic graph output screen G1 and a piping system 3D data output screen G2. The tube axis coordinate-physical statistic graph output screen G1 displays a graph (physical statistic visualization graph) of the turning intensity with respect to the tube axis coordinates in each of the path 1 and the path 2, with the path 1 in the upper row and the path in the lower row. 2 is shown. The horizontal axis shows the pipe axis coordinates, and the vertical axis shows the turning strength. The time average of turning strength, the average value of turning strength plus the standard deviation, and the average value of turning strength minus the standard deviation. Is drawing. Each part is color-coded, and the part name is described in the part area.

The piping system 3D data output screen G2 is drawn by projecting the three-dimensional shape of the piping system in two dimensions from a specific viewpoint, and each component is color-coded and drawn. In addition, the character string of the part name is drawn near each part, and the arrow of the guideline and the route name are also drawn for the route.

In the present embodiment configured as described above, since the physical statistics visualization graph with the tube axis coordinates on the horizontal axis and the physical statistics on the vertical axis is drawn, a three-dimensional fluid analysis is performed. The designer can easily judge the validity of the result from the simulation result.

<Modified example of the first embodiment>
Hereinafter, this modification will be described with reference to FIG. 31.

This modification is a case where the hardware configuration is different from that of the first embodiment.

FIG. 31 is a diagram showing the hardware configuration of the simulation result visualization system in this modified example. In the figure, the same members as those in the first embodiment are designated by the same reference numerals, and the description thereof will be omitted.

In FIG. 31, the simulation result visualization system D100A includes a piping system design unit D200, a piping system design operation input unit D106A, a piping system design display unit D105A, a 3D fluid analysis simulation unit D300, and a simulation operation input unit D106B. It is composed of a simulation display unit D105B, a visualization unit D400, a visualization operation input unit D106C, a visualization display unit D105C, and a network D500.

As described above, the piping system design unit D200, the 3D fluid analysis simulation unit D300, and the visualization unit D400 have operation input units D106A to D106C and display units D105A to D105C, respectively. Each operation input unit D106A to D106C converts an operation input from a user into an electric signal and transmits it to each core unit, and corresponds to a keyboard, a touch panel, a mouse, or the like. Further, the display units D105A to D105C are for outputting the image output from each core unit to the user, and refer to a display or a printer.

The piping system design unit D200, 3D fluid analysis simulation unit D300 and visualization unit D400 are CPU (Central Processor Unit) D201 to D401, RAM (Random Access Memory) D202 to D402, recording devices D203 to D403, and I / F, respectively. (Interface) It is composed of D204 to D404. The CPUs D201 to D401 execute arithmetic processing according to a program for visualizing the simulation result. The RAM_D202 to D402 temporarily store data in a program for visualizing the simulation result. The recording devices D203 to D403 record programs, data, and the like for processing each part, and correspond to a hard disk drive or the like.

The recording device D203 of the piping system design unit D200 has a piping system design program, the recording device D303 of the 3D fluid analysis simulation unit D300 has a 3D fluid analysis simulation program, and the recording device D403 of the visualization unit D400 has a visualization program. It is stored.

The functions of the piping system program, the 3D fluid analysis simulation program, and the visualization program are the same as those in the first embodiment. However, data is exchanged between each part through a network.

Other configurations are the same as in the first embodiment.

In the present modification configured as described above, the same effect as that of the first embodiment can be obtained.

<Additional notes>
The present invention is not limited to the above-described embodiment, and includes various modifications and combinations within a range that does not deviate from the gist thereof. Further, the present invention is not limited to the one including all the configurations described in the above-described embodiment, and includes the one in which a part of the configurations is deleted. Further, each of the above configurations, functions and the like may be realized by designing a part or all of them by, for example, an integrated circuit. Further, each of the above configurations, functions, and the like may be realized by software by the processor interpreting and executing a program that realizes each function.

B1 ... Piping system design unit, B11 ... Parts creation unit, B12 ... Parts storage unit, B13 ... Parts selection / connection unit, B14 ... Piping configuration storage unit, B15 ... Piping route information generation unit, B16 ... Piping route information storage unit, B2 ... 3D fluid analysis simulation unit, B21 ... simulation condition setting unit, B22 ... simulation model selection unit, B23 ... piping route information reading unit 1, B24 ... simulation calculation unit, B25 ... simulation result storage unit, B3 ... visualization unit, B31 ... Piping route information management unit, B311 ... Piping route information reading unit 2, B312 ... In-route parts information extraction / additional setting unit, B3121 ... Parts information extraction unit, B3122 ... Parts type extraction unit, B3123 ... Parts color setting unit, B313 ... In-route part information storage unit, B314 ... Path cross-section definition unit, B315 ... In-route cross-section information storage unit, B32 ... Physical statistics calculation unit, B321 ... Simulation result reading unit 1, B322 ... In each part on the route Sectional information reading section, B323 ... Simulation result section section extraction section, B324 ... Cell unit physical statistic calculation section, B325 ... Sectional unit physical statistic integration calculation section, B326 ... Physical statistics time average / standard deviation calculation section, B327 ... Section Unit physical statistics storage unit, B33 ... Pipe axis coordinate calculation unit, B331 ... Section information reading unit in each part on the route, B332 ... Each section pipe axis coordinate calculation unit, B333 ... Section information in each part on the route with pipe axis coordinates Storage unit, B34 ... Graph output unit, B341 ... Sectional information reading unit for each part on the route with tube axis coordinates 1, B342 ... Sectional unit physical statistic reading unit, B343 ... Physical statistics on each part section / tube axis coordinates ( Horizontal axis vertical axis) Extraction unit, B344 ... Graph drawing unit, B345 ... Graph setting unit, B346 ... Graph display output unit, B35 ... Three-dimensional data output unit, B351 ... Section information reading unit in each part on the route with tube axis coordinates 2, B352 ... Simulation result reading unit 2, B353 ... Each component mesh extraction unit, B354 ... 3D data drawing unit, B355 ... Drawing setting unit, B356 ... 3D data display output unit, D1 ... Various data, D11 ... Piping system design program Data, D111 ... Parts mesh information, D1111 ... Parts information, D112 ... Parts list, D113 ... Piping configuration information, D114 ... Piping configuration list, D12 ... Piping route information, D13 ... Piping route list, D14 ... 3D fluid analysis simulation program Data, D141 ... Simulation condition / model list, D15 ... Simulation result, D16 ... Data for visualization program , D161 ... Cross section information of each part on the path, D162 ... Cross section unit simulation result, D163 ... Cross section physical statistic information, F1 ... Piping system design process, F2 ... Fluid analysis simulation execution process, F3 ... Part information extraction / additional setting process , F4 ... Path cross section definition processing, F5 ... Physical statistics / tube axis coordinate calculation processing, F6 ... Graph / 3D data output processing, G1 ... Pipe axis coordinates-physical statistics graph output screen, G2 ... Piping system 3D data output screen

Claims (10)

The piping system design information holding unit that holds information on the start coordinates, end coordinates, length in the pipe axis direction, and diameter of each piping component in a piping system consisting of multiple piping components, and the fluid phenomenon in the piping system. In a simulation result visualization system having a fluid analysis simulation unit that calculates a physical quantity related to a fluid phenomenon at each time based on a simulated physical formula, and a visualization unit that summarizes and draws the information of the physical quantity and outputs the drawn contents.
The visualization unit
A physical statistic calculation unit that calculates a physical statistic from the physical quantity,
A pipe axis coordinate calculation unit that calculates pipe axis coordinates in the piping system,
A simulation result visualization system characterized by having a graph drawing unit for drawing a physical statistic visualization graph with the tube axis coordinates on the horizontal axis and physical statistics on the vertical axis. A piping system design information holding device that holds information on the start coordinates, end coordinates, length in the pipe axis direction, and diameter of each piping component in a piping system consisting of multiple piping components, and fluid phenomena in the piping system. A fluid analysis simulation device that calculates the physical quantity related to the fluid phenomenon at each time based on the simulated physical formula, a visualization device that summarizes and draws the information of the physical quantity and outputs the drawn contents, and data transmission between the devices. In a simulation result visualization system with a network device to perform
The visualization device is
A physical statistic calculation unit that calculates a physical statistic from the physical quantity,
A pipe axis coordinate calculation unit that calculates pipe axis coordinates in the piping system,
A simulation result visualization system characterized by having a graph drawing unit for drawing a physical statistic visualization graph with the tube axis coordinates on the horizontal axis and physical statistics on the vertical axis. In the simulation result visualization system according to claim 1 or 2.
The simulation result visualization system is characterized in that a graph having pipe axis coordinates on the horizontal axis and physical statistics on the vertical axis is displayed for each path of the piping system in the physical statistics visualization graph. In the simulation result visualization system according to claim 1 or 2.
The simulation result visualization system is characterized in that the three-dimensional shape of the piping system is drawn together with the graph having the pipe axis coordinates on the horizontal axis and the physical statistics on the vertical axis in the physical statistics visualization graph. In the simulation result visualization system according to claim 4,
A simulation result visualization system characterized in that the three-dimensional shape of the piping system is color-coded for each part. In the simulation result visualization system according to claim 1 or 2.
A simulation result visualization system characterized in that a component name is drawn in a display area of each component in the three-dimensional shape of the piping system. In the simulation result visualization system according to claim 1 or 2.
A simulation result visualization system characterized in that the physical statistic is a physical statistic calculated from a swirl number. In the simulation result visualization system according to claim 1 or 2.
A simulation result visualization system characterized in that the physical statistic is a time average of the physical statistic. In the simulation result visualization system according to claim 1 or 2.
A simulation result visualization system characterized in that the physical statistic is a time average, a maximum value, and a minimum value of the physical statistic. In the simulation result visualization system according to claim 1 or 2.
The simulation result visualization system is characterized in that the physical statistic is a time average related to the physical statistic, a time average plus a time standard deviation, and a time average minus the time standard deviation.

Images