Simulation device, simulation method, and simulation program

Abstract

To efficiently simulate the phenomenon of melting a solid heated in a melting furnace. [Solution] The system includes an initial value setting unit that sets initial values ​​for particle method data showing air and metal particles in the melting furnace, and Euler method data representing the space of the melting furnace including air and metal divided into a mesh, based on structural information of the melting furnace and the metal placed in the melting furnace; and an update unit that updates the particle method data of metal and air based on the particle method, updates the Euler method data based on the Euler method, and updates the particle method data of air based on the Euler method data.

Description

Technical Field

[0001] The present invention relates to a simulation device and the like.

Background Art

[0002] In metal processing such as casting, in the process of melting the metal in the melting furnace, heating by air may be performed.

[0003] FIG. 41 is a diagram showing an example of melting metal in a melting furnace. As shown in FIG. 41, a hot air outlet 5a is provided in the melting furnace 5. In step S1, the metal 6 before melting is arranged, and the release of hot air from the outlet 5a is started. In step S2 after the passage of time from step S1, the metal 6 is melting.

[0004] It is required to make the melting time of the metal as short as possible, and as a simulation technique for predicting the melting time of the metal, a particle method, an Euler method, etc. are used.

[0005] The particle method divides substances such as fluids, solids, and gases into a plurality of particles for modeling. On the other hand, the Euler method divides the simulation space into a plurality of meshes and predicts the change in physical quantities within each mesh.

Prior Art Documents

Patent Documents

[0006]

Patent Document 1

Patent Document 2

Summary of the Invention

Problems to be Solved by the Invention

[0007] However, predicting the melting time of a metal using the particle method results in long computation times. Since melting time is known to depend on multiple parameters such as the metal's shape and wind direction, it is difficult to investigate the conditions for parameters that can shorten the melting time while changing these multiple parameters using the particle method.

[0008] While the Euler method can shorten calculation time when predicting metal melting time compared to the particle method, its prediction accuracy is inferior.

[0009] In one aspect, the present invention aims to provide a simulation device, a simulation method, and a simulation program that can efficiently simulate the phenomenon of melting a solid heated in a melting furnace. [Means for solving the problem]

[0010] In the first proposal, the simulation device has an initial value setting unit that sets initial values ​​for particle method data showing air and metal particles in the melting furnace, and Euler method data representing the space of the melting furnace including air and metal divided into a mesh, based on structural information of the melting furnace and the metal placed in the melting furnace; and an update unit that updates the particle method data of metal and air based on the particle method, updates the Euler method data based on the Euler method, and updates the particle method data of air based on the Euler method data. [Effects of the Invention]

[0011] This technology allows for efficient simulation of the phenomenon of melting a solid heated in a melting furnace. [Brief explanation of the drawing]

[0012] [Figure 1] Figure 1 is a diagram illustrating the particle method. [Figure 2] Figure 2 shows an example of the data structure of particle method data. [Figure 3]FIG. 3 is a flowchart showing a processing procedure for one-hour step calculation of central coordinates by the particle method. [Figure 4] FIG. 4 is a flowchart showing the specific processing procedure shown in step S14. [Figure 5] FIG. 5 is a diagram showing an example of a list of nozzle outlet particle numbers. [Figure 6] FIG. 6 is a flowchart showing a processing procedure for one-hour step calculation of temperature T by the particle method. [Figure 7] FIG. 7 is a diagram for explaining the Euler method. [Figure 8] FIG. 8 is a diagram showing an example of the data structure of Euler method data. [Figure 9] FIG. 9 is a diagram showing an example of a list of nozzle outlet mesh numbers. [Figure 10] FIG. 10 is a flowchart showing the calculation procedure by the Euler method. [Figure 11] FIG. 11 is a functional block diagram showing the configuration of the simulation device according to this embodiment. [Figure 12] FIG. 12 is a diagram for explaining the stl format. [Figure 13] FIG. 13 is a flowchart showing a processing procedure for initializing particle method data. [Figure 14] FIG. 14 is a diagram (1) for explaining the processing of the initial setting unit. [Figure 15] FIG. 15 is a diagram (2) for explaining the processing of the initial setting unit. [Figure 16] FIG. 16 is a diagram (3) for explaining the processing of the initial setting unit. [Figure 17] FIG. 17 is a diagram (4) for explaining the processing of the initial setting unit. [Figure 18] FIG. 18 is a flowchart showing the initial setting process of temperature. [Figure 19] FIG. 19 is a diagram (5) for explaining the processing of the initial setting unit. [Figure 20]Figure 20 is a flowchart showing the processing procedure for initializing Euler method data. [Figure 21] Figure 21 is a diagram (6) for explaining the processing of the initialization section. [Figure 22] Figure 22 is a diagram (7) for explaining the processing of the initialization section. [Figure 23] Figure 23 is a flowchart showing the initial setting process of the flow velocity. [Figure 24] Figure 24 is a diagram (8) for explaining the processing of the initialization section. [Figure 25] Figure 25 is a flowchart showing the processing procedure of the first update process. [Figure 26] Figure 26 is a flowchart showing the processing procedure for selecting particle number j. [Figure 27] Figure 27 is a diagram for supplementarily explaining the first update process. [Figure 28] Figure 28 is a diagram for supplementarily explaining the second update process. [Figure 29] Figure 29 is a flowchart showing the processing procedure of the third update process. [Figure 30] Figure 30 is a diagram (1) for supplementarily explaining the third update process. [Figure 31] Figure 31 is a flowchart showing the process of updating the particle method data based on the Euler method data. [Figure 32] Figure 32 is a diagram showing an example of the metal shape data. [Figure 33] Figure 33 is a diagram showing an example of the data structure of the metal shape data. [Figure 34] Figure 34 is a diagram for supplementarily explaining the fourth update process. [Figure 35] Figure 35 is a flowchart showing the processing procedure of the fifth update process. [Figure 36] Figure 36 is a diagram for supplementarily explaining the fifth update process. [Figure 37] Figure 37 is a diagram for explaining the change amount of the metal shape. [Figure 38]Figure 38 is a flowchart showing the processing procedure for the sixth update process. [Figure 39] Figure 39 is a flowchart showing the processing procedure of the simulation device according to this embodiment. [Figure 40] Figure 40 shows an example of a computer hardware configuration that achieves similar functionality to the simulation device described in the embodiment. [Figure 41] Figure 41 shows an example of melting metal in a melting furnace. [Modes for carrying out the invention]

[0013] The following describes in detail, with reference to the drawings, embodiments of the simulation apparatus, simulation method, and simulation program disclosed in this application. However, this embodiment does not limit the present invention. [Examples]

[0014] Before describing this embodiment, we will explain the particle method and the Euler method in turn.

[0015] First, let's explain an example of the particle method. Figure 1 is a diagram illustrating the particle method. In the particle method, metals and air are represented as aggregates of tiny particles. Furthermore, in the particle method, the change in the shape of a metal due to melting is represented by the movement of the particles that make up the metal. For example, the phenomenon of the shape of metal 7a changing to the shape of metal 7b is treated in the particle method as the particles in particle method data 8a, which represent an aggregate of particles, moving and changing to particle method data 8b.

[0016] Although not shown in Figure 1, the particle method data 8a and 8b represent not only the metals 7a and 7b, but also the air contained in the molten furnace as aggregates of particles.

[0017] Figure 2 shows an example of the data structure of particle method data. As shown in Figure 2, particle method data includes a particle number, center coordinates x, y, z, temperature T, metal flag M, and flow velocity v. The particle number is a number that identifies the particle. The center coordinates x, y, z indicate the center position of the particle. Temperature T indicates the temperature of the particle. The metal flag M is a flag that indicates whether the particle is metallic or not. If the particle is metallic, the metal flag M is "true", and if the particle is air, the metal flag M is "false". Flow velocity v indicates the flow velocity of the particle. A value is set for flow velocity v when the particle is air. If the particle is metallic, a value such as NULL is set to indicate emptiness.

[0018] In the particle method, the particle's center coordinates and temperature are calculated at predetermined time steps, and the particle method data is updated. Below, the calculation of the center coordinates for one hour step and the calculation of the temperature for one hour step using the particle method are explained in order. For convenience of explanation, the device that performs the particle method will be referred to as "device." The time width dt for the one-hour step is predetermined.

[0019] First, we will explain the calculation of the central coordinates for one hour step using the particle method. The apparatus performs the calculation using the particle method data shown in Figure 2. Figure 3 is a flowchart showing the processing procedure for the calculation of the central coordinates for one hour step using the particle method. In Figure 3, particles for which the metal flag M is true are referred to as "metal particles," and particles for which the metal flag M is false are referred to as "air particles." The melting point Tm explained in Figure 3 is set in advance.

[0020] As shown in Figure 3, the device sets i to the particle number of the first particle in the particle method data (step S10). If the metal flag M for particle number i is not true (step S11, No), the device proceeds to step S16. On the other hand, if the metal flag M for particle number i is true (step S11, Yes), the device proceeds to step S12.

[0021] If the temperature T of particle number i is not equal to or greater than the melting point Tm (step S12, No), the apparatus proceeds to step S16. On the other hand, if the temperature T of particle number i is equal to or greater than the melting point Tm (step S12, Yes), the apparatus proceeds to step S13.

[0022] The device calculates the change in the central coordinates (x, y, z) of the metal particle with particle number i based on the particle method (step S13). If an air particle is present at the destination of the metal particle, the device moves the air particle to the source of the metal particle (step S14).

[0023] The device updates the center coordinates of the metal particle with particle number i in the particle method data (step S15). If particle number i reaches the last particle number in the particle method data (step S16, Yes), the device terminates the calculation of the center coordinates for one hour step using the particle method.

[0024] If particle number i has not reached the last particle number in the particle method data (step S16, No), the device sets particle number i to the next particle number (step S17) and proceeds to step S11.

[0025] Next, the processing procedure shown in step S14 of Figure 3 will be explained in more detail. Figure 4 is a flowchart showing the specific processing procedure shown in step S14. The "i" shown in Figure 4 corresponds to the particle number i explained in Figure 3. The particle diameter Lp explained in Figure 4 is set in advance.

[0026] As shown in Figure 4, the device sets the particle number of the first particle in the particle method data to j (step S20). If the metal flag M for particle number j is not false (step S21, No), the device proceeds to step S25. On the other hand, if the metal flag M for particle number j is false (step S21, Yes), the device proceeds to step S22.

[0027] The device calculates the interparticle distance between the center coordinates of particle number i and the center coordinates of particle number j (step S22). If the interparticle distance is not greater than or equal to the particle diameter Lp (step S23), the device updates the center coordinates of particle number i in the particle method data with the center coordinates of particle number j (step S24).

[0028] We will now proceed to the explanation of step S25. The device terminates processing if particle number j reaches the last particle number in the particle method data (step S25, Yes). On the other hand, if particle number j has not reached the last particle number in the particle method data (step S25, No), the device sets particle number j to the next particle number (step S26) and proceeds to step S21.

[0029] Next, we will explain the calculation of temperature T for a one-hour step using the particle method. When performing such calculations, the apparatus uses the outlet particle number list shown in Figure 5. Figure 5 is a diagram showing an example of the outlet particle number list. As shown in Figure 5, the outlet particle number list 15 is set with particle numbers that identify each particle present in the outlet 20b when the entire melting furnace 20 containing the metal 20a is represented as an aggregate of particles. In the example shown in Figure 5, for convenience, the particle numbers set in the outlet particle number list 15 are p_out1, p_out2, p_out3, etc., but the same particle is assigned the same particle number as the particle number in Figure 2.

[0030] Figure 6 is a flowchart showing the processing procedure for calculating the temperature T for a 1-hour step using the particle method. As shown in Figure 6, the device sets i to the particle number of the first particle in the particle method data (step S30). If the metal flag M for particle number i is true (step S31, Yes), the device calculates the temperature T of particle number i using the center coordinates and temperature T of the particle method data based on the particle method (step S32), and proceeds to step S35.

[0031] On the other hand, if the metal flag M for particle number i is not true (step S31, No), the device determines whether particle number i matches the particle number in the outlet particle number list 15 (step S33). If particle number i matches the particle number in the outlet particle number list 15 (step S33, Yes), the device proceeds to step S35.

[0032] On the other hand, if particle number i does not match the particle number in the outlet particle number list 15 (step S33, No), the device calculates the temperature T of particle number i based on the particle method, using the center coordinates, temperature T, and flow velocity v of the particle method data (step S34).

[0033] The device updates the temperature T of particle number i in the particle method data (step S35). If particle number i reaches the last particle number in the particle method data (step S36, Yes), the device terminates the process. On the other hand, if particle number i has not reached the last particle number in the particle method data (step S36, No), the device sets particle number i to the next particle number (step S37) and proceeds to step S31.

[0034] This concludes the explanation of the particle method.

[0035] Next, an example of the Euler method will be described. Figure 7 is a diagram illustrating the Euler method. As shown in Figure 7, the Euler method divides the space to be simulated (the entire melting furnace 20 including the metal 20a and the outlet 20b) into a mesh, and calculates the flow velocity and pressure of each mesh. In this invention, both the metal and the air are meshed.

[0036] The Euler method uses the Euler method data shown in Figure 8. Figure 8 is a diagram illustrating an example of the data structure of the Euler method data. As shown in Figure 8, the Euler method data includes a mesh number, center coordinates x, y, z, flow velocity v, pressure p, temperature T, fluid resistance r, and an air flag a. The mesh number is a number that identifies the mesh. The center coordinates x, y, z indicate the center position of the mesh. Flow velocity v, pressure p, temperature T, and fluid resistance r are the flow velocity, pressure, temperature, and fluid resistance of the mesh, respectively.

[0037] The air flag 'a' indicates whether the mesh with the corresponding mesh number should be treated as air or metal. The air flag is 'true' when the mesh is treated as air, and 'false' when it is treated as metal. For example, a mesh is treated as metal if its center coordinates are within the coordinates of one of the particles in the particle method data that make up the metal (a sphere with a particle radius at the particle's center coordinates).

[0038] The flow velocity v is 0 when the mesh is metal. The fluid resistance r is 0 when the mesh corresponds to air, and when it is metal, it is a sufficiently large value (for example, around 10 to the power of 9) that prevents air from moving.

[0039] As will be described later, in this invention, with respect to the temperature T of Euler method data, if the center coordinates of the mesh are contained within any particle (a sphere with a particle radius at the particle's center coordinates) of the particle collection constituting the metal in the particle method data, the temperature of that particle is set as the mesh temperature T.

[0040] In the Euler method data, the temperature T and flow velocity v of the mesh corresponding to the outlet are set to predetermined fixed values. For example, the mesh corresponding to the outlet is the mesh corresponding to outlet 20b shown in Figure 7. General Euler methods have a function (boundary condition) to set such fixed values, and this function is used.

[0041] In the calculation using the Euler method, the outlet mesh number list shown in Figure 9 is used. Figure 9 is a diagram showing an example of the outlet mesh number list. As shown in Figure 9, the outlet mesh number list 16 is set to a mesh number that identifies each mesh present in the outlet 20b when the entire melting furnace 20, including the metal 20a, is represented as a collection of meshes. In the example shown in Figure 9, for convenience, the mesh numbers set in the outlet mesh number list 16 are mb1, mb2, mb3, etc., but the same mesh is assigned the same mesh number as the mesh number in Figure 8.

[0042] Next, we will explain the procedure for calculation using the Euler method. For the sake of explanation, the device that performs the Euler method will be referred to as the "device". Figure 10 is a flowchart showing the procedure for calculation using the Euler method. Based on the outlet mesh number list 16, the device sets a fixed value for the flow velocity v corresponding to the outlet mesh number among the mesh numbers of the Euler method data (step S40). The following processing steps will not be explained, but the fixed value set for the flow velocity v corresponding to the outlet mesh number is not updated.

[0043] The device identifies the shape of the metal based on the air flag in the Euler method data (step S41). The device calculates the steady-state flow velocity v and pressure p for the entire mesh space in the shape of the metal using the Euler method (step S41). Based on the calculation results, the device updates the flow velocity and pressure for each mesh number in the Euler method data (step S42).

[0044] The device sets i to the first mesh number in the Euler method data (step S43). The device updates the flow velocity v and pressure p for mesh number i in the Euler method data (step S44). If mesh number i reaches the last mesh number in the Euler method data (step S45, Yes), the device terminates the process. On the other hand, if mesh number i has not reached the last mesh number in the Euler method data (step S45, No), the device sets mesh number i to the next mesh number (step S46) and proceeds to step S44.

[0045] This concludes the explanation of Euler's method.

[0046] Next, the simulation apparatus according to this embodiment will be described. The simulation apparatus according to this embodiment performs initial settings for particle method data and Euler method data based on structural information regarding the melting furnace and the metal inside the melting furnace. The simulation apparatus updates the information of the particles corresponding to the metal among the multiple particles in the particle method data based on the particle method, and updates the information of the particles corresponding to air among the multiple particles in the particle method data based on the mesh information of the updated Euler method data based on the Euler method. In this way, by using the particle method and the Euler method appropriately, the phenomenon of melting a solid heated in a melting furnace can be efficiently simulated.

[0047] An example of the configuration of the simulation device according to this embodiment will be described. Figure 11 is a functional block diagram showing the configuration of the simulation device according to this embodiment. As shown in Figure 11, the simulation device 100 has a communication unit 110, an input unit 120, a display unit 130, a storage unit 140, and a control unit 150.

[0048] The communication unit 110 communicates with external devices via a network. The communication unit 110 is implemented using a NIC (Network Interface Card) or the like.

[0049] The input unit 120 is an input device for inputting various types of information to the control unit 150 of the simulation device 100. The input unit 120 can be a keyboard, mouse, touch panel, etc.

[0050] The display unit 130 is a display device that displays information output from the control unit 150.

[0051] The storage unit 140 includes particle method data 141, Euler method data 142, outlet particle number list 143, outlet mesh number list 144, metal STL file 145, melting furnace STL file 146, outlet STL file 147, and metal shape table 148. The storage unit 140 is a storage device such as an HDD (Hard Disk Drive) or memory.

[0052] Particle method data 141 is data used in the particle method. The explanation of particle method data 141 is the same as the explanation of particle method data in Figure 2.

[0053] Euler method data 142 is data used in the Euler method. The explanation of Euler method data 142 is the same as the explanation of Euler method data in Figure 3.

[0054] The outlet particle number list 143 shows the particle numbers from the particle method data 141 that correspond to the particles at the outlet. The explanation for the outlet particle number list 143 is the same as the explanation for the outlet particle number list 15 in Figure 5.

[0055] The outlet mesh number list 144 shows the mesh numbers corresponding to the outlet meshes among the mesh numbers in the Euler method data. The explanation for outlet mesh number list 144 is the same as the explanation for outlet mesh number list 16 in Figure 9.

[0056] Metal STL file 145 is data that defines the surface of a three-dimensional shape of a metal. Melting furnace STL file 146 is data that defines the surface of a three-dimensional shape of a melting furnace. For example, metal STL file 145 and melting furnace STL file 146 are STL format files generated by 3D CAD. Figure 12 is a diagram illustrating the STL format. For example, an STL format file is data in which the surface of a three-dimensional shape 17 is divided into multiple triangles, and the vertex coordinates of each triangle are recorded in the file.

[0057] The outlet STL file 147 is data that defines the three-dimensional surface shape of the outlet of the melting furnace. Like the metal STL file 145 and the melting furnace STL file 146, the outlet STL file 147 is also an STL file generated by 3D CAD.

[0058] For example, metal STL file 145 and melting furnace STL file 146 correspond to "structural information".

[0059] The metal shape table 148 holds the metal shape data calculated from the particle method data 141. For example, the metal shape data is generated by scanning the center coordinates x, y, and z of the particle method data 141 where the metal flag M is "true", and extracting the outermost particle (particle number). The metal shape data will be explained later.

[0060] We will now move on to the explanation of the control unit 150 in Figure 11. The control unit 150 has an initial setting unit 151, an update unit 152, and an output unit 153. The control unit 150 is a CPU (Central Processing Unit), etc.

[0061] The initial setup unit 151 performs the initial setup of particle method data 141 and Euler method data 142. The process by which the initial setup unit 151 performs the initial setup of particle method data 141 and Euler method data 142 will be described in order below.

[0062] The process by which the initial setup unit 151 initializes the particle method data 141 will be described. First, the processing procedure of the initial setup unit 151 will be explained, and then supplementary information will be provided for each process. Figure 13 is a flowchart of the processing procedure for initializing the particle method data. As shown in Figure 13, the initial setup unit 151 reads the metal STL file 145 and the melting furnace STL file 146 (step S50).

[0063] The initial setup unit 151 determines the center coordinates of the particles using the metal stl file 145 and records them in the particle method data 141 in association with the particle number (step S51). The initial setup unit 151 sets the metal flag M to true and the flow velocity v to NULL for the particle number recorded in the particle method data (step S52).

[0064] The initial setup unit 151 newly determines the particle number and center coordinates of the particles using the metal STL file 145 and the melting furnace STL file 146, and adds the newly determined particle number and center coordinates to the particle method data 141 (step S53).

[0065] The initial setup unit 151 sets the metal flag M to false and the flow velocity v to v0 for the particle numbers newly recorded in the particle method data 141 (step S54). The initial setup unit 151 sets the temperature T for all particle numbers in the particle method data 141 to the initial value (step S55).

[0066] Next, we will provide a supplementary explanation of the process shown in Figure 13. We will provide a supplementary explanation of the process in step S51 of Figure 13. The initial setup unit 151 uses the metal stl file 145 to determine the central coordinates (x, y, z) so that the entire region inside the metal is filled with particles of diameter Lp. The initial setup unit 151 uses the functions set in general particle method software to determine the central coordinates (x, y, z) of each particle. The initial setup unit 151 records the central coordinates (x, y, z) of each determined particle in the particle method data 141 and assigns particle numbers in the order in which they are recorded.

[0067] Figure 14 is a diagram (1) illustrating the processing of the initial setup unit. In Figure 14, the metal surface 18 is identified based on the metal stl file 145. For example, the initial setup unit 151 sets n particles inside the metal surface 18 and determines the center coordinates (x, y, z) of each particle. The initial setup unit 151 assigns particle numbers R1, R2, ..., Rn to the n particles and records them in the particle method data 141 in correspondence with the center coordinates (x, y, z).

[0068] Next, we will provide a supplementary explanation of the process in step S52 of Figure 13. Figure 15 is Figure (2) for explaining the process of the initial setup unit. As shown in Figure 15, the initial setup unit 151 sets the metal flag M to true for particle numbers R1, R2, ..., Rn recorded in the particle method data 141, and sets the flow velocity v to NULL.

[0069] Next, a supplementary explanation will be given regarding step S53 in Figure 13. Figure 16 is Figure (3) for explaining the processing of the initial setup unit. In Figure 16, the metal surface 18 is identified based on the metal stl file 145. The melting furnace surface 19 is identified based on the melting furnace stl file 146. The initial setup unit 151 targets the region 19a from inside the melting furnace surface 19, excluding the inside of the metal surface 18, and identifies the central coordinates (x, y, z) so that the entire region 19a is filled with particles of diameter Lp. The initial setup unit 151 identifies the central coordinates (x, y, z) of each particle using the functions set in general particle method software. The initial setup unit 151 stores the central coordinates (x, y, z) of each particle in the identified region 19a in the particle method data 141 and assigns a new particle number in the order in which they are recorded. For example, the initial setup unit 151 assigns particle numbers Rn+1, Rn+2, ..., Rn+m to m particles, and records them in the particle method data 141, corresponding to the center coordinates (x, y, z).

[0070] Next, we will provide supplementary explanations for step S54 in Figure 13. Figure 17 is Figure (4) for explaining the processing of the initial setup unit. The initial setup unit 151 sets the metal flag M to false and the flow velocity v to v0 for the particle numbers Rn+1, Rn+2, ..., Rn+m newly recorded in the particle method data 141.

[0071] Next, we will provide supplementary explanations regarding step S55 in Figure 13. For example, as an explanation of step S55 in Figure 13, the initial setting unit 151 more specifically performs the temperature initial setting process shown in Figure 18. Figure 18 is a flowchart of the temperature initial setting process.

[0072] As shown in Figure 18, the initial setup unit 151 reads the outlet stl file 147 (step S60). The initial setup unit 151 sets i to the particle number of the first particle in the particle method data (step S61). If the center coordinates of particle number i are not located inside the outlet (step S62, No), the initial setup unit 151 sets the temperature T of particle number i in the particle method data 141 to the value of T0 (step S63) and proceeds to step S66.

[0073] On the other hand, if the center coordinates of particle number i are located inside the outlet (step S62, Yes), the initial setup unit 151 records particle number i in the outlet particle number list 143 (step S64). The initial setup unit 151 sets Tout to the temperature T of particle number i in the particle method data 141 (step S65).

[0074] The initial setup unit 151 terminates processing if particle number i reaches the last particle number in the particle method data 141 (step S66, Yes). On the other hand, if particle number i has not reached the last particle number in the particle method data (step S66, No), the initial setup unit 151 sets particle number i to the next particle number (step S67) and proceeds to step S62.

[0075] For example, the initial setup unit 151 updates the particle method data 141 as shown in Figure 19 by performing the process shown in Figure 18. Figure 19 is a diagram (5) illustrating the process of the initial setup unit. For example, if the center coordinates of particle numbers R1, R2, Rn, Rn+1, Rn+m are not located inside the outlet, the initial setup unit 151 sets the temperature T of particle numbers R1, R2, Rn, Rn+1, Rn+m in the particle method data 141 to T0.

[0076] The initial setup unit 151 sets Tout to the temperature T corresponding to particle number Rn+2 in the particle method data 141 if the center coordinates of particle number Rn+2 are located inside the outlet. The initial setup unit 151 also records particle number Rn+2 in the outlet particle number list 143.

[0077] The above describes the process by which the initial setup unit 151 initializes the particle method data 141.

[0078] Next, we will explain the process by which the initial setup unit 151 initializes the Euler method data 142. First, we will explain the processing procedure of the initial setup unit 151, and then provide supplementary information for each process. Figure 20 is a flowchart of the processing procedure for initializing the Euler method data. As shown in Figure 20, the initial setup unit 151 reads the melting furnace stl file 146 (step S70).

[0079] The initial setup unit 151 uses the melting furnace stl file 146 to divide the entire melting furnace into multiple meshes (rectangular parallelepipeds), determines the center coordinates of each mesh, and records the center coordinates in Euler method data 142 in correspondence with the mesh number (step S71).

[0080] The initial setup unit 151 sets the pressure p corresponding to all mesh numbers in the Euler method data 142 to p0 (step S72). The initial setup unit 151 determines and records the flow velocity v corresponding to all mesh numbers in the Euler method data 142 (step S73).

[0081] Next, we will provide a supplementary explanation of the process shown in Figure 20. We will provide a supplementary explanation of the process in step S71 of Figure 20. The initial setup unit 151 divides the entire melting furnace into multiple meshes (rectangular parallelepipeds) of sizes dx, dy, and dz using the melting furnace stl file 146. General mesh generation software can be used for meshing the entire melting furnace. The initial setup unit 151 records the center coordinates (x, y, z) of each mesh in the Euler method data 142 and assigns mesh numbers in the order in which they are recorded.

[0082] Figure 21 is a diagram (6) illustrating the processing of the initial setup unit. In Figure 21, the melting furnace surface 19 is identified based on the melting furnace stl file 146. For example, the initial setup unit 151 divides the melting furnace surface 19 into multiple meshes (rectangular parallelepipeds) of sizes dx, dy, and dz. The initial setup unit 151 records the center coordinates (x, y, z) of each mesh in Euler method data 142 and assigns mesh numbers m1, m2, m3, ..., mn in the order in which they are recorded.

[0083] Next, we will provide a supplementary explanation of the process in step S72 of Figure 20. Figure 22 is a diagram (7) illustrating the process of the initial setup unit. In Figure 22, the initial setup unit 151 sets pO for the pressures p of mesh numbers m1, m2, m3, ..., mn recorded in the Euler method data 142.

[0084] Next, we will provide a supplementary explanation regarding the process in step S73 of Figure 20. For example, the initial setup unit 151, as the process in step S73 of Figure 20, more specifically, performs the initial flow velocity setting process shown in Figure 23. Figure 23 is a flowchart of the initial flow velocity setting process.

[0085] As shown in Figure 23, the initial setup unit 151 reads the outlet stl file 147 (step S80). The initial setup unit 151 sets i to the first mesh number of the Euler method data 142 (step S81). If the center coordinates of mesh number i are not inside the outlet (step S82, No), the initial setup unit 151 sets the value of v0 to the velocity v of mesh number i in the Euler method data 142 (step S83) and proceeds to step S86.

[0086] On the other hand, if the center coordinates of mesh number i are located inside the outlet (step S82, Yes), the initial setup unit 151 records mesh number i in the outlet mesh number list 144 (step S84). The initial setup unit 151 sets vout to the velocity v of mesh number i in the Euler method data 142 (step S85).

[0087] The initial setup unit 151 terminates processing if the mesh number i reaches the last mesh number in the Euler method data 142 (step S86, Yes). On the other hand, if the initial setup unit 151 has not reached the last mesh number in the Euler method data (step S86, No), it sets the mesh number i to the next mesh number (step S87) and proceeds to step S82.

[0088] For example, the initial setup unit 151 updates the Euler method data 142 as shown in Figure 24 by performing the process shown in Figure 23. Figure 24 is a diagram (8) illustrating the process of the initial setup unit. For example, if the center coordinates of mesh numbers m1, m3, mn are not located inside the outlet, the initial setup unit 151 sets the velocity v of mesh numbers m1, m3, mn in the Euler method data 142 to v0.

[0089] The initial setup unit 151 sets vout to the velocity v corresponding to mesh number m2 in the Euler method data 142 if the center coordinates of mesh number m2 are located inside the outlet. The initial setup unit 151 also records mesh number m2 in the outlet mesh number list 144.

[0090] The above describes the process by which the initial setup unit 151 performs the initial setup of the particle method data 141 and the Euler method data 142.

[0091] Let's return to the explanation of Figure 11. The update unit 152 executes the first update process, the second update process, the third update process, the fourth update process, the fifth update process, and the sixth update process. The following explanation will describe the first to sixth update processes in order.

[0092] First, the first update process performed by the update unit 152 will be described. The first update process is a process that updates the Euler method data 142 based on the particle method data 141. For example, with respect to the temperature T of the Euler method data 142, if the center coordinates of the mesh are contained within any particle (a sphere with a particle radius at the particle's center coordinates) of the particle method data 141 that constitutes the metal, the update unit 152 sets the temperature of that particle as the mesh temperature T.

[0093] Figure 25 is a flowchart showing the processing procedure for the first update process. As shown in Figure 25, the update unit 152 sets i to the starting mesh number of the Euler method data 142 (step S90).

[0094] The update unit 152 searches the particle method data 141 for a particle with a center coordinate nearest to the center coordinate of mesh number i, and sets the particle number of the found particle to j (step S91). If the metal flag for particle number j is not true (step S92, No), the update unit 152 sets the temperature T of particle number j to the temperature T of mesh number i, sets the fluid resistance r to 0, sets the air flag a to true (step S93), and proceeds to step S95.

[0095] On the other hand, if the metal flag for particle number j is true (step S92, Yes), the update unit 152 sets the flow velocity v for mesh number i to 0, sets the temperature T to the value of temperature T for particle number j, sets the fluid resistance r to R0, and sets the air flag a to false (step S94).

[0096] The update unit 152 terminates processing if the mesh number i reaches the last mesh number in the Euler method data 142 (step S95, Yes). On the other hand, if the mesh number i has not reached the last mesh number in the Euler method data 142 (step S95, No), the update unit 152 sets the mesh number i to the next mesh number (step S96) and proceeds to step S91.

[0097] Here, the process of selecting the particle number j shown in step S91 of FIG. 25 will be described more specifically. FIG. 26 is a flowchart showing the procedure for selecting the particle number j. As shown in FIG. 26, the update unit 152 sets a sufficiently large real value (for example, 10 to the 9th power) for lmin (step S100). The update unit 152 sets 0 for jtmp (step S101). The update unit 152 sets k to the first particle number of the particle method data 141 (step S102).

[0098] The update unit 152 calculates the distance ltmp between the center coordinates of the mesh number i and the center coordinates of the particle number k (step S103). If the update unit 152 does not satisfy the condition of ltmp < lmin (step S104, No), it proceeds to step S107.

[0099] On the other hand, if the update unit 152 satisfies the condition of ltmp < lmin (step S104, Yes), it sets the value of lmin to the value of ltmp (step S105). The update unit 152 sets k to the value of jtmp (step S106).

[0100] If the update unit 152 reaches the last particle number of the particle method data 141 (step S107, Yes), it sets the value of jtmp to j (step S108). On the other hand, if the update unit 152 has not reached the last particle number of the particle method data 141 (step S107, No), it sets the particle number k to the next particle number (step S109) and proceeds to step S103.

[0101] Next, we will provide supplementary information about the first update process described in Figure 25. Figure 27 is a diagram that provides supplementary explanation of the first update process. In the example shown in Figure 27, we will explain the case where the determination in step 92 of Figure 25 is "Yes" in the processing of mesh number m3 of the Euler method data 142. In this case, the update unit 152 sets the flow velocity v corresponding to mesh number m3 to "0 (zero)", the fluid resistance r to "R0", and the air flag a to "false". The update unit 152 also sets the temperature T corresponding to mesh number m3 to the value of the temperature T of the particle (particle number) with the center coordinates nearest to the center coordinates of mesh number m3 (for example, T3).

[0102] Next, the second update process performed by the update unit 152 will be explained. The second update process updates the flow velocity v and pressure p of the Euler method data 142 by executing the calculation procedure of the Euler method described in Figure 10.

[0103] Figure 28 is a diagram that provides supplementary explanation of the second update process. The update unit 152 updates the values ​​of flow velocity v and pressure p by performing calculations using the Euler method on the Euler method data 142. Note that since the mesh with mesh number m2 is inside the outlet, the flow velocity v for mesh number m2 is vout. Since the fluid resistance r of the mesh with mesh number m2 is R0, the flow velocity remains 0 (zero).

[0104] Next, the third update process performed by the update unit 152 will be described. The third update process is a process that updates the particle method data 141 based on the Euler method data 142. Figure 29 is a flowchart showing the processing procedure of the third update process. As shown in Figure 29, the update unit 152 sets i to the particle number of the first particle in the particle method data 141 (step S200).

[0105] If the metal flag M for particle number i is not false (step S201, No), the update unit 152 proceeds to step S204. On the other hand, if the metal flag M for particle number i is false (step S201, Yes), the update unit 152 proceeds to step S202.

[0106] The update unit 152 sets the mesh number j to the one whose center coordinates are closest to the center coordinates of particle number i among the records in the Euler method data 142 where the air flag a is true (step S202). The update unit 152 records the value of the flow velocity v of mesh number j as the flow velocity v of particle number i (step S203).

[0107] If particle number i reaches the last particle number in the particle method data (step S204, Yes), the update unit 152 generates metal shape data based on the particle method data 141 (step S205).

[0108] On the other hand, if particle number i has not reached the last particle number in particle method data 141 (step S204, No), the update unit 152 sets particle number i to the next particle number (step S206) and proceeds to step S201.

[0109] Figure 30 is a diagram that provides supplementary explanation of the third update process. In the example shown in Figure 30, the update unit 152 sets the flow velocity v value of the corresponding mesh number for the flow velocity v corresponding to particle numbers Rn+1, Rn+2, and Rn+3 in the particle method data 141 records where the metal flag M is false. For example, vn+1 is set for the flow velocity v corresponding to particle number Rn+1. vn+2 is set for the flow velocity v corresponding to particle number Rn+2. vn+3 is set for the flow velocity v corresponding to particle number Rn+3.

[0110] Here, the process described in step S202 of FIG. 29 will be described more specifically with reference to FIG. 31. FIG. 31 is a flowchart showing the process of updating particle method data based on Euler method data. As shown in FIG. 31, the update unit 152 sets a sufficiently large real value (for example, 10 to the 9th power) for lmin (step S300). The update unit 152 sets 0 for jtmp (step S301). The update unit 152 sets the mesh number of the head of the Euler method data 142 for k (step S302).

[0111] If the air flag a of the mesh number k is not true (step S303, No), the update unit 152 proceeds to step S308. On the other hand, if the air flag a of the mesh number k is true (step S303, Yes), the update unit 152 proceeds to step S304.

[0112] The update unit 152 calculates the distance ltmp between the center coordinates of the particle number i and the center coordinates of the mesh number k (step S304). If the condition ltmp < lmin is not satisfied (step S305, No), the update unit 152 proceeds to step S308.

[0113] On the other hand, if the condition ltmp < lmin is satisfied (step S305, Yes), the update unit 152 sets the value of lmin to the value of ltmp (step S306). The update unit 152 sets k to the value of jtmp (step S307).

[0114] If the mesh number k has reached the last mesh number of the Euler method data 142 (step S308, Yes), the update unit 152 sets the value of jtmp to j (step S309). On the other hand, if the mesh number k has not reached the last mesh number of the Euler method data 142 (step S308, No), the update unit 152 sets the mesh number k to the next particle number (step S310) and proceeds to step S303.

[0115] Next, we will explain in more detail the process for generating the metal shape data, as described in step S205 of Figure 29.

[0116] The update unit 152 extracts metal shape data based on the particle method data 141. Figure 32 shows an example of metal shape data. For example, the update unit 152 scans the center coordinates x, y, z of the particle method data 141 where the metal flag M is "true" and extracts the outermost particle (particle number) to generate metal shape data 9b.

[0117] Figure 33 shows an example of the data structure of metal shape data. As shown in Figure 33, the metal shape data includes a particle number and central coordinates x, y, and z. The particle number in Figure 33 is a number that identifies a particle, and is a number that identifies a particle located on the metal surface. In particle method data 141 and the metal shape data in Figure 33, the same particle is assigned the same particle number. The central coordinates x, y, and z indicate the central position of the particle.

[0118] The update unit 152 generates metal shape data and registers the generated metal shape data in the metal shape table 148. The metal shape data generated by the update unit 152 through the third update process is referred to as the "first metal shape data".

[0119] Next, the fourth update process executed by the update unit 152 will be explained. The fourth update process updates the "temperature T" of the particle method data 141 by executing the calculation procedure of the particle method described in Figure 3. For example, in the fourth update process, the change in temperature T of the particle method data 141 for a one-hour step is calculated using the particle method.

[0120] Figure 34 is a diagram to supplement the explanation of the fourth update process. When the update unit 152 calculates the temperature for a one-hour step using the particle method for the particle method data 141, the temperature T portion of the particle method data 141 is updated accordingly. In the example in Figure 34, the particle with particle number Rn+2 is located inside the outlet, so the temperature remains fixed at Tout.

[0121] Next, the fifth update process executed by the update unit 152 will be described. The fifth update process calculates the "center coordinates x, y, z" for one hour step of the particle method data 141 by executing the calculation procedure of the particle method described in Figure 4, and further generates metal shape data. Figure 35 is a flowchart showing the processing procedure of the fifth update process.

[0122] As shown in Figure 35, the update unit 152 calculates the center coordinates of the particle method data 141 for one hour step using the particle method (step S400) and updates the particle method data 141 (step S401). Based on the particle method data 141, the update unit 152 generates metal shape data (step S402).

[0123] The update unit 152 refers to the metal shape data generated by the fifth update process as "second metal shape data". The update unit 152 records the second metal shape data in the metal shape table 148.

[0124] Figure 36 is a diagram that provides supplementary explanation of the fifth update process. In the example shown in Figure 36, the update unit 152 updates the central coordinates x, y, and z corresponding to particle numbers R1 to Rn+m of the particle method data 141. The update unit 152 scans the central coordinates x, y, and z corresponding to the updated particle numbers R1 to Rn+m for which the metal flag M is "true", and generates the second metal shape data by extracting the outermost particle (particle number).

[0125] Next, the sixth update process executed by the update unit 152 will be described. In the sixth update process, the update unit 152 calculates the amount of deformation of the metal shape based on the first metal shape data and the second metal shape data recorded in the metal shape table 148.

[0126] Figure 37 is a diagram illustrating the change in metal shape. A particle corresponding to a certain "particle number RN (where N is a natural number)" is sometimes simply referred to as "particle RN".

[0127] In Figure 37, metal shape data 10a is described as the first metal data, and metal shape data 10b is described as the second metal shape data. The update unit 152 selects a particle (here, particle R101) from the metal shape data 10b. The update unit 152 also selects a particle (here, particle R101) from the metal shape data 10a. The update unit 152 calculates the distance between the center coordinates of particle R101 in the metal shape data 10a and the center coordinates of particle R101 in the metal shape data 10b as the change amount.

[0128] The update unit 152 calculates the total amount of change in all particles while changing the pairs of particles in the metal shape data 10a, and identifies the maximum amount of change among the calculated amounts of change as the change amount dl.

[0129] Figure 38 is a flowchart showing the processing steps of the sixth update process. As shown in Figure 38, the update unit 152 sets Imax to 0.0 (step S500). The update unit 152 sets Itmp to 0.0 (step S501). Imax and Itmp are real number variables. The update unit 152 sets i to 0 (step S502). i is a variable.

[0130] The update unit 152 calculates the distance between the center coordinates of the i-th particle in the first metal shape data and the center coordinates of the i-th particle in the second metal shape data, and sets it to the real variable Itmp (step S503).

[0131] If the update unit 152 does not satisfy the condition Itmp > Imax (step S504, No), it proceeds to step S506. On the other hand, if the update unit 152 satisfies the condition Itmp > Imax (step S504, Yes), it proceeds to step S505.

[0132] The update unit 152 updates the value of Imax based on the value of Itmp (Imax = Itmp) (step S505). The update unit 152 increments i by 1 (step S506). If the condition i > N is not met (step S507, No), the update unit 152 proceeds to step S503. N is the maximum number of records of the metal shape data.

[0133] On the other hand, if the condition i > N is met (step S507, Yes), the update unit 152 sets the change amount dl to the value of Imax (dl = Imax) (step S508).

[0134] The update unit 152 calculates the change amount dl by executing the process shown in Figure 38. If the change amount dl is greater than or equal to a threshold, the update unit 152 sequentially executes the following: "the process of updating Euler method data 142 based on particle method data 141", "the process of performing calculations on Euler method data 142 using the Euler method", and "the process of updating particle method data 141 based on Euler method data". After that, it updates the first metal shape data with the metal shape data generated in the "process of updating particle method data 141 based on Euler method data", and returns to the fourth update process described above.

[0135] The process of "updating Euler method data 142 based on particle method data 141" is the same as the first update process described above. The process of "performing Euler method calculations on Euler method data 142" is the same as the second update process described above. The process of "updating particle method data 141 based on Euler method data" is the same as the third update process described above.

[0136] Furthermore, if the change amount dl is less than or equal to the threshold, the update unit 152 returns to the fourth update process described above.

[0137] The update unit 152 repeatedly executes the above process in a time width of 1 hour step dt until the time reaches a predetermined value.

[0138] The output unit 153 outputs the particle method data 141 and Euler method data 142, which have been updated by the update unit 152, to the display unit 130 for display as simulation results.

[0139] The above describes an example of the configuration of the simulation device 100 according to this embodiment.

[0140] Next, an example of the overall processing procedure of the simulation device 100 according to this embodiment will be described. Figure 39 is a flowchart showing the processing procedure of the simulation device according to this embodiment. Note that the values ​​for the step time dt and the end time tend shown in Figure 39 are assumed to be predetermined. As shown in Figure 39, the simulation device 100 sets the time t to an initial value (for example, 0) (step S600).

[0141] The initial setup unit 151 of the simulation device 100 initializes the particle method data 141 (step S601). The specific process in step S601 corresponds to the processing procedure described in Figure 13.

[0142] The initial setup unit 151 initializes the Euler method data (step S602). The specific processing procedure for step S602 corresponds to the processing procedure described in Figure 20.

[0143] The update unit 152 of the simulation device 100 updates the temperature T, fluid resistance r, and air flag a of the Euler method data 142 based on the particle method data 141 (step S603). The specific processing procedure of step S603 corresponds to the processing procedure of the first update process described in Figure 25.

[0144] The update unit 152 performs calculations using the Euler method and updates the flow velocity v and pressure p of the Euler method data (step S604). The specific processing procedure for step 604 corresponds to the processing procedure for the second update process described in Figure 28.

[0145] The update unit 152 updates the particle method data 141 based on the Euler method data 142 (step S605). The specific processing procedure for step S605 corresponds to the processing procedure for the third update process described in Figure 29.

[0146] The update unit 152 generates first metal shape data based on the particle method data 141 (step S606).

[0147] The update unit 152 calculates the change in temperature T of the particle method data 141 for a one-hour step using the particle method and updates the particle method data 141 (step S607). The specific processing procedure for step S607 corresponds to the fourth update process. The fourth update process is a process that updates the "temperature T" of the particle method data 141 by executing the calculation procedure using the particle method described in Figure 6.

[0148] The update unit 152 calculates the change in the center coordinates of the particle method data 141 for a one-hour step using the particle method and updates the particle method data (step S608). The specific processing procedure for step S608 corresponds to the fifth update process. The fifth update process is a process that updates the "center coordinates x, y, z" of the particle method data 141 by executing the calculation procedure using the particle method described in Figure 3.

[0149] The update unit 152 generates second metal shape data based on the particle method data (step S609). If second metal shape data already exists in step S609, it is updated with the newly generated second metal shape data. The update unit 152 calculates the amount of change based on the first metal shape data and the second metal shape data (step S610).

[0150] If the amount of change is greater than or equal to the threshold (step S611, Yes), the update unit 152 proceeds to step S612. On the other hand, if the amount of change is not greater than or equal to the threshold (step S611, No), the update unit 152 proceeds to step S616.

[0151] The update unit 152 updates the temperature T, fluid resistance r, and air flag a of the Euler method data 142 based on the particle method data 141 (step S612). The specific processing procedure of step S612 corresponds to the processing procedure of the first update process described in Figure 25.

[0152] The update unit 152 performs calculations using the Euler method to update the flow velocity v and pressure p data using the Euler method (step S613). The specific processing procedure for step 613 corresponds to the processing procedure for the second update process described in Figure 28. Since the update of the flow velocity v and pressure p using the Euler method is performed only when the amount of change in the metal shape is greater than or equal to a threshold, it is possible to significantly reduce the simulation time without reducing the prediction accuracy.

[0153] The update unit 152 updates the particle method data 141 based on the Euler method data 142 (step S614). The specific processing procedure for step S614 corresponds to the processing procedure for the third update process described in Figure 29.

[0154] The update unit 152 generates first metal shape data based on the particle method data 141, updates the first metal shape data (step S615), and proceeds to step S616.

[0155] The update unit 152 outputs the simulation result (step S617) if the condition t > tend is met (step S616, Yes). On the other hand, if the update unit 152 does not meet the condition t > tend (step S616, No), it updates the time t by t = t + dt (step S618) and proceeds to step S607.

[0156] Next, the effects of the simulation device 100 according to this embodiment will be explained. The simulation device 100 performs initial settings for particle method data 141 and Euler method data 142 based on structural information regarding the melting furnace and the metal inside the melting furnace. The simulation device 100 updates the particle method data 141 of the metal based on the particle method, updates the Euler method data 142 based on the Euler method, and updates the particle method data 141 of the air based on the Euler method data 142. In this way, by using the particle method and the Euler method appropriately, the phenomenon of melting a solid heated in a melting furnace can be efficiently simulated.

[0157] For example, focusing on computation time, while the conventional technology requires a computation time of "4005.85h," the simulation device 100 according to this embodiment requires a computation time of "9.85h," representing an expected speed improvement of more than 407 times. However, this assumes a real time of "6000s" for the particle method, a time step of "0.0015s," and a computation time of 3.6s for one step in the Euler method.

[0158] Furthermore, the simulation device 100 repeatedly performs a process to calculate the temperature and position of particles in the particle method data 141 based on the particle method at predetermined time intervals. When the amount of change in the metal based on the particle method data 141 exceeds a threshold, it updates the flow velocity of the mesh in the Euler method data based on the Euler method, and updates the flow velocity of the particle method data with the updated mesh flow velocity. This allows the particle method data and the Euler method data 142 to be updated at appropriate timings.

[0159] Furthermore, the simulation device 100 updates the temperature of the mesh corresponding to the metal among the multiple meshes included in the Euler method data using the temperature of the particles corresponding to the metal in the particle method data. This allows the Euler method data 142 to be updated appropriately.

[0160] Next, an example of a computer hardware configuration that achieves the same functions as the simulation device 100 shown in the above embodiment will be described. Figure 40 is a diagram showing an example of a computer hardware configuration that achieves the same functions as the simulation device according to the embodiment.

[0161] As shown in Figure 40, the computer 200 includes a CPU 201 that performs various calculations, an input device 202 that receives data input from the user, and a display 203. The computer 200 also includes a communication device 204 that communicates with external devices via a wired or wireless network, and an interface device 205 that connects to other devices. The computer 200 also includes a RAM 206 for temporarily storing various information and a hard disk drive 207. Each of the devices 201 to 207 is connected to a bus 208.

[0162] The hard disk drive 207 contains an initial setup program 207a, an update program 207b, and an output program 207c. The CPU 201 reads each of the programs 207a to 207c and loads them into RAM 206.

[0163] The initial setup program 207a functions as the initial setup process 206a. The update program 207b functions as the update process 206b. The output program 207c functions as the output process 206c.

[0164] The processing of the initial setup process 206a corresponds to the processing of the initial setup unit 151. The processing of the update process 206b corresponds to the processing of the update unit 152. The processing of the output process 206c corresponds to the processing of the output unit 153.

[0165] Furthermore, programs 207a to 207c do not necessarily have to be stored on the hard disk drive 207 from the beginning. For example, each program could be stored on a "portable physical medium" such as a flexible disk (FD), CD-ROM, DVD, magneto-optical disk, or IC card inserted into the computer 200. Then, the computer 200 could read and execute each program 207a to 207c. [Explanation of Symbols]

[0166] 100 Simulation devices 120 Input section 130 Display section 140 Storage section 141 Particle method data 142 Euler method data 143 List of particle numbers for air outlets 144 Outlet Mesh Number List 145 Metal STL File 146 Melting furnace STL file 147 Air outlet stl file 148 Metal Shaped Table 150 Control Unit 151 Initial setting section 152 Update section 153 Output section

Claims

1. An initial value setting unit sets initial values ​​for particle method data showing particles of air and metal in the melting furnace, and Euler method data representing the space of the melting furnace including the air and metal, based on structural information of the melting furnace and the metal placed in the melting furnace. An update unit updates the particle method data for the metal and the air based on the particle method, updates the Euler method data based on the Euler method, and updates the particle method data for the air based on the Euler method data. A simulation device characterized by having the following features.

2. The particle information in the aforementioned particle method data includes temperature, position, and flow velocity. The simulation apparatus according to claim 1, characterized in that the update unit repeatedly performs a process to calculate the temperature and position of particles in the particle method data at predetermined time intervals, updates the flow velocity of the mesh in the Euler method data when the amount of change in the particle method data of the metal exceeds a threshold, and further performs a process to update the flow velocity of the particle method data using the updated flow velocity of the mesh.

3. The simulation apparatus according to claim 2, wherein the update unit further performs a process to calculate the amount of change based on the position of the particles in the particle method data calculated at a certain time and the position of the particles in the particle method data calculated at the time following the certain time.

4. The mesh information in the aforementioned Euler method data includes temperature and flow velocity. The simulation apparatus according to claim 2, characterized in that the update unit further performs a process of updating the temperatures of a plurality of meshes included in the Euler method data with the temperatures of the particle method data.

5. Based on structural information of the melting furnace and the metal placed in the melting furnace, initial values ​​are set for particle method data showing the air and metal particles in the melting furnace, and for Euler method data representing the space of the melting furnace, including the air and the metal, divided into a mesh. Based on the particle method, the particle method data for the metal and the air were updated. The Euler method data is updated based on the Euler method. A simulation method characterized in that a computer performs a process to update the particle method data of air based on the Euler method data.

6. Based on structural information of the melting furnace and the metal placed in the melting furnace, initial values ​​are set for particle method data showing the air and metal particles in the melting furnace, and for Euler method data representing the space of the melting furnace, including the air and the metal, divided into a mesh. Based on the particle method, the particle method data for the metal and the air were updated. The Euler method data is updated based on the Euler method. A simulation program characterized by causing a computer to perform a process to update the particle method data of air based on the Euler method data.

Images