KR20130042055A - Temperature estimation method and device for fluid system - Google Patents

Abstract

The temperature estimation device acquires the flow field of the fluid system. Subsequently, the temperature estimating apparatus is a region including each of the temperature measurement sites A to D, the endothermic region E, and the inflow / exit sites F and G in the fluid system, and the regions E31 to E34, E, which do not overlap each other, Set F and G. Subsequently, the temperature estimating apparatus occupies a proportion of the fluid generated in the flow diffusion phenomenon due to the flow field in the total fluid of the temperature estimation point of the fluid that has reached the temperature estimation point without passing through other regions. Obtained as a downstream force for every G. And a temperature estimating apparatus estimates the temperature of a temperature estimation point using the downstream side force for every site | part in a temperature estimation point based on the known temperature of each site | part A-G.

Description

TEMPERATURE ESTIMATION METHOD AND DEVICE FOR FLUID SYSTEM

The present invention provides a method for estimating a temperature of a fluid system, a method for estimating a temperature distribution of a fluid system, a method for monitoring a temperature distribution of a fluid system, a temperature estimating apparatus, a method for controlling a molten zinc temperature in a hot dip galvanizing port, a hot dip galvanized steel sheet, and a tundish It relates to a molten steel temperature control method in tundish.

In order to know directly the temperature distribution of the fluid system flowing through the walls separated by the wall of the fluid equipment or the building used in the industrial process, the temperature measuring device is arranged in a number and arrangement sufficient to characterize the temperature distribution in the fluid system. Need to install However, the fluid installation, the inside of a building, etc. often have a complicated shape, and there exist some places where a temperature measuring device cannot be installed. In addition, when the fluid system is high temperature or highly corrosive, the measurement of the temperature using the temperature measuring device may be limited, and the temperature measuring device may not be installed in a number and arrangement sufficient to know the temperature distribution. many. In order to make up for the lack of the temperature measurement site in the case where the temperature measuring device cannot be provided in this way, it is necessary to estimate the temperature distribution of the entire fluid system from the measured temperature measurement value and interpolate.

In the case where the object for estimating the temperature distribution is a homogeneous solid, the temperature correlation between the temperature measurement point and the temperature estimation point in the geometrically close position is large, so that it is relatively easy to use a well-known interpolation method such as spline interpolation. It is possible to estimate the temperature distribution and interpolate. For example, a method called an inverse distance weighting method is known which estimates various values including the temperature of the estimated point and the like based on the distance between the measured part and the estimated point. (For example, refer nonpatent literature 1). The method of nonpatent literature 1 calculates the distance between the position of a measured site | part, and an estimated point, weights as the weight becomes smaller, so that the value of the measured site with a calculated distance is large, and as a weighted average The value at the estimation point is estimated. Specifically, the method described in Non-Patent Literature 1 estimates the temperature of a temperature estimation point, for example, by using the following equation (1) with a weight of the power of the inverse of the distance (l / l _i ). do. Here, Te ^l is the estimated temperature in the temperature estimation point, l _i is the distance between the position of the temperature measurement site i and the temperature estimation point, and T _i is the temperature in the temperature measurement site i Actual value. U is an interpolation parameter that takes a positive value.

[1]

On the other hand, when the object for estimating the temperature distribution is a fluid system, heat transfer occurs due to convection, so that even if the temperature measurement site and the temperature estimation point are geometrically close to each other, the temperature correlation is not large. . As a method of estimating the temperature distribution of such a fluid system, for example, Patent Document 1 discloses a temperature distribution in a storage tank based on the moving speed of hot water and the history of measured temperatures in a hot-water cylinder. A temperature distribution estimation system for estimating is disclosed. Patent Document 2 also superimposes tidal vector data obtained by a tidal current meter together with water temperature data obtained by a track and a water temperature meter which are position data obtained by a navigation apparatus. A water temperature distribution display device is disclosed. In addition, Patent Literature 3 establishes boundary conditions on the basis of actual values of sensors installed indoors, and expresses environmental conditions such as temperature, humidity, and carbon dioxide concentration at designated locations, or formulas related to thermal conductivity or Navier-Stokes. Disclosed is an air conditioning sensor system estimated using the following equation.

On the other hand, as a technique relating to the flow field of a fluid system, for example, Non-Patent Document 2, Non-Patent Document 3, and Patent Document 4, the power range of the outlet and the suction port used as one of the indicators of the ventilation efficiency The concept is described. When the power range of the outlet described in this nonpatent literature 2, the nonpatent literature 3, and the patent document 4 paid attention to any specific point of the room provided with the some outlet, it is the airflow from the outlet to consider. Indicates how far is reaching the point. In addition, the range of the force of an intake port represents the distribution state in each indoor point of the air discharged | emitted through the intake port made into examination object. Moreover, the nonpatent literature 2, the nonpatent literature 3, and the patent document 4 also describe the calculation method of the range of forces by a numerical analysis.

Japanese Laid-Open Patent Publication No. 2006-214622 Japanese Laid-Open Patent Publication No. 61-151428 Japanese Laid-Open Patent Publication No. 2008-75973 Japanese Laid-Open Patent Publication No. 2004-101058

 Shepard, D .: A two-dimensional interpolating function for irregularly spaced data. Proc. ACM. nat. Conf., 517e24, 1968.  Shuka Murakami: Environmental design engineering of architecture, city by CFD, Tokyo society publishing society  S. Kato, S. Murakami, H. Kobayashi: Newscales for evaluating ventilation efficiency as affected by supply and exhaust openings based on spatial distribution of contaminant, Proceedings of the 12th ISCC, 341-348, 1994

However, in the reverse distance weighting method described in Non-Patent Literature 1, since the temperature is estimated, for example, by a weighted average by a weight based only on the distance l _i between the temperature measurement site and the temperature estimation point, the temperature obtained. The estimation result does not reflect the influence of the flow of the fluid. For this reason, in a fluid system in which the contribution of the flow to heat transport such as an actual fluid process is very large, the same temperature distribution is estimated even though the temperature distribution is significantly different in the case where the flow velocity is large and small. Therefore, it is difficult to apply to the fluid system in which heat transport by the flow of the fluid is dominant.

Moreover, the system of patent document 1 applies only a 1-dimensional fluid flow, and it is difficult to apply to the fluid system which has a 3-dimensional fluid flow field. In addition, although Patent Document 2 describes that there is a correlation between the distribution of water temperature, the direction of the algae, and its speed, it is not disclosed how to specifically obtain the correlation to estimate the water temperature distribution, and thus requires a fluid system. The accuracy of estimation of the non-measurement point to be obtained cannot be obtained. In addition, since the water temperature data and the tidal current vector data are assumed only for the case of two-dimensional sea surface, application to a fluid system having a three-dimensional fluid flow field is difficult.

In addition, as in the system described in Patent Literature 3, in the method using a formula related to heat conduction or a Navy Stokes formula, the sensor is located at both the most upstream positions of the flow of air, such as a window into which the air enters the room and an outlet of the air conditioning equipment. It is essential to install. For this reason, when a sensor cannot be installed in all of the most upstream positions, it becomes impossible to estimate a temperature distribution in the whole fluid system. Therefore, it is difficult to apply to a fluid system in which the sensor cannot be installed at the inflow position of the fluid.

In addition, the methods described in Non-Patent Literature 2 and Non-Patent Literature 3 focus only on visualizing the presence position of the fluid inflow and outflow from the intake and intake ports of air, and assume application to temperature estimation. I'm not doing it. Similarly, the method described in Patent Document 4 focuses on calculating an air age spatial distribution, and does not assume application to temperature estimation.

The present invention has been made in view of the above, and a temperature estimation method of a fluid system capable of realizing a high precision temperature estimation in consideration of heat transportation due to flow of a fluid without restricting the arrangement of the temperature measuring device, An object of the present invention is to provide a temperature distribution estimation method of a fluid system, a temperature distribution monitoring method of a fluid system, and a temperature estimation device.

In addition, another object of the present invention is a molten zinc produced using a molten zinc temperature control method in a hot dip galvanizing port capable of manufacturing a hot dip galvanized steel sheet without surface defects and a molten zinc temperature control method in the hot dip galvanizing port. It is to provide a coated steel sheet. Moreover, another object of this invention is to provide the molten steel temperature control method in a tundish which can suppress the refractory damage of a tundish.

In order to solve the above-mentioned problems and to achieve the object, the temperature estimation method of the fluid system according to the present invention includes a temperature at an arbitrary temperature estimation point of a fluid system having two or more temperature known regions. A method of estimating a temperature of a fluid system, the method of estimating a temperature, wherein the temperature is passed through a temperature known area by using position information of the temperature known area and information about a flow field of a fluid system indicating a flow of fluid in the entire fluid system. The ratio of the fluid generated in the temperature base region to the temperature estimated point of the fluid which reaches the temperature estimation point without passing through the other temperature base region is the ratio of the force of the temperature base region at the temperature estimation point. The temperature using the force acquiring step acquired as a function, and information on the forces in the temperature known areas and the forces at the temperature estimation points. It characterized in that it comprises a temperature estimation step of estimating the temperature at the apex.

In order to solve the above problems and to achieve the object, the temperature distribution estimation method of a fluid system according to the present invention is a temperature distribution estimation method of a fluid system having a temperature distribution, using the above-described invention throughout the entire fluid system. It is characterized by estimating the temperature of the temperature estimation point set in the above, and estimating the temperature estimated for each of the temperature estimation points as the temperature distribution of the fluid system.

MEANS TO SOLVE THE PROBLEM In order to achieve the objective and to achieve the objective, the temperature distribution monitoring method of the fluid system which concerns on this invention is a temperature distribution monitoring method of the fluid system which has a temperature distribution, Comprising: The said fluid system estimated using said invention. Based on the temperature distribution of, the temperature distribution in any cross section of the fluid system is visualized and displayed on the screen.

In order to solve the above problems and achieve the object, the temperature estimating apparatus according to the present invention is a temperature estimating apparatus for estimating the temperature at an arbitrary temperature estimation point of a fluid system having two or more temperature known regions. Among the fluids having passed through the temperature known area or generated in the temperature known area by using the position information of the temperature known area and information about the flow field of the fluid system indicating the flow of the fluid throughout the fluid system, Force acquisition means which acquires the ratio which the fluid which reached the said temperature estimation point in the whole fluid of the temperature estimation point as a force of the temperature known area in a temperature estimation point, without passing through a temperature base area, and each temperature base The temperature which estimates the temperature in the said temperature estimation point using the information of the temperature of an area | region and the force in the said temperature estimation point. An estimation means is provided.

In order to solve the above problems and to achieve the object, the molten zinc temperature control method in the hot dip galvanizing port according to the present invention, the molten zinc in the hot dip galvanizing port estimated by the temperature estimation method of the fluid system according to the present invention A temperature extraction step of extracting the temperature of the molten zinc in a predetermined region in the hot dip galvanizing port from the zinc temperature data, a determination step of determining whether the extracted temperature is within a predetermined threshold range, and the determination In the step, if it is determined that the extracted temperature is out of the threshold range, a control step of manipulating the output of the heating means of the hot dip galvanizing port so that the extracted temperature is within the threshold range is characterized in that it comprises a.

In order to solve the above-mentioned problems and to achieve the object, the hot-dip galvanized steel sheet according to the present invention is produced using a hot-dip zinc temperature control method in the hot-dip galvanizing port according to the present invention.

In order to solve the above-mentioned problems and achieve the object, the molten steel temperature control method in the tundish according to the present invention includes the molten steel temperature data in the tundish estimated by the temperature estimation method of the fluid system according to the present invention. A temperature extraction step of extracting the temperature of the molten steel in a predetermined region in the tundish; a determination step of determining whether the extracted temperature is within a predetermined threshold range; and in the determination step, the extracted temperature is a threshold And when it is determined to be out of the value range, a control step of operating the output of the heating means of the tundish so that the extracted temperature is within a threshold range.

According to the present invention, it is possible to realize a high-precision temperature estimation in consideration of heat transportation due to the flow of fluid without restricting the arrangement of the temperature measuring device. According to this invention, a hot-dip galvanized steel plate without surface defects can be manufactured. According to the present invention, refractory damage of the tundish can be suppressed.

1 is a block diagram illustrating the concept of the present invention.
2 is a block diagram showing an apparatus configuration for implementing the present invention.
3A is a diagram illustrating an example of a fluid system.
3B is a diagram showing a downstream force distribution in the temperature known region R1.
3C is a diagram showing the downstream force distribution in the temperature known region R2.
3D is a diagram showing a downstream force distribution in the temperature known region R3.
3E is a diagram showing the upstream force distribution of the temperature known region R1.
3F is a diagram showing the upstream force distribution in the temperature known region R2.
3G is a diagram showing the upstream force distribution of the temperature known region R3.
4 is a flowchart showing a processing procedure of a weight calculation process in the case where the temperature field of the fluid system does not change with time.
FIG. 5 is a flowchart showing a processing procedure of a weight calculation process following FIG. 4.
FIG. 6 is a flowchart showing a processing procedure of weight calculation processing following FIG. 5.
7 is a flowchart showing a processing procedure of a temperature estimation process.
8 is a flowchart showing a processing procedure of a delivery time calculating process.
9 is a flowchart showing a processing procedure of a weight calculation process in the case where the temperature field of the fluid system can change over time.
FIG. 10 is a flowchart showing a processing procedure of a weight calculation process following FIG. 9.
FIG. 11 is a flowchart showing a processing procedure of a weight calculation process subsequent to FIG. 10.
12 is a flowchart showing a processing procedure of the temperature estimation process in the case where the temperature field of the fluid system can change with time.
FIG. 13 is a schematic view showing an interior of a room to be applied in Embodiment 1 from above. FIG.
14 is a block diagram showing the functional configuration of the temperature estimating apparatus according to the first embodiment.
15 is a schematic diagram showing an actual flow field of the room of FIG. 13.
It is a figure which shows the downstream force of the temperature measurement site | part A. FIG.
It is a figure which shows the downstream force of the temperature measurement site | part B.
It is a figure which shows the downstream force of the temperature measurement site | part C.
It is a figure which shows the downstream force of the temperature measurement site | part D. FIG.
It is a figure which shows the force downstream of the endothermic part E. FIG.
It is a figure which shows the downstream force of the inflow-exit part F. FIG.
It is a figure which shows the downstream force of the inflow-exit part G. FIG.
It is a figure which shows the upstream force of the temperature measurement site | part A. FIG.
It is a figure which shows the upstream force of the temperature measurement site | part B.
It is a figure which shows the upstream force of the temperature measurement site | part C.
It is a figure which shows the upstream force of the temperature measurement site | part D.
It is a figure which shows the upstream force of the endothermic area E. FIG.
It is a figure which shows the upstream force of the inflow-exit part F. FIG.
It is a figure which shows the upstream force of the inflow-exit part G. FIG.
It is a figure which shows the weight of the temperature measurement site | part A. FIG.
It is a figure which shows the weight of the temperature measurement site | part B.
It is a figure which shows the weight of the temperature measured part C.
It is a figure which shows the weight of the temperature measurement site | part D. FIG.
18E is a diagram showing the weight of the endothermic portion E. FIG.
It is a figure which shows the weight of the inflow_out | pouring part F. FIG.
It is a figure which shows the weight of the inflow_out | portion part G.
19A is another diagram illustrating the weight of the temperature measured part A. FIG.
19B is another diagram illustrating the weight of the temperature measured part B. FIG.
19C is another diagram illustrating a weight of the temperature measured part C. FIG.
19D is another diagram illustrating a weight of the temperature measured part D. FIG.
19E is another diagram showing the weight of the endothermic portion E;
19F is another diagram illustrating the weight of inflow and outflow region F. FIG.
19G is another diagram illustrating the weight of the inflow and outflow region G. FIG.
20 is a diagram illustrating the estimation result of Experimental Example 1 according to the first embodiment.
FIG. 21 is a diagram showing an estimation result of Experimental Example 2 in Embodiment 1. FIG.
22 is a diagram illustrating the estimation result of the comparative example in the first embodiment.
FIG. 23 is a diagram showing the actual temperature distribution in the room of FIG.
It is a schematic diagram which showed the inside of the water tank made into the application object in Embodiment 2 from the side.
FIG. 25 is a schematic view showing the inside of the water tank of FIG. 24 from above. FIG.
Fig. 26 is a block diagram showing the functional configuration of the temperature estimating apparatus according to the second embodiment.
FIG. 27 is a diagram showing an estimation result of Experimental Example 1 in Embodiment 2. FIG.
FIG. 28 is a diagram showing an estimation result of Experimental Example 2 in Embodiment 2. FIG.
It is a figure which shows the estimation result of the comparative example in Embodiment 2. FIG.
It is a figure which shows the additional installation position of the thermometer in a water tank.
FIG. 31A is a diagram showing a time course of temperature measured at position P41 in the water tank. FIG.
FIG. 31B is a diagram showing a time course of temperature measured at position P42 in the water tank. FIG.
FIG. 31C is a diagram showing the time course of the temperature measured at the position P43 in the water tank. FIG.
FIG. 31D is a diagram showing the time course of the temperature measured at the position P44 in the water tank.
FIG. 31E is a diagram showing a time course of temperature measured at position P45 in the water tank. FIG.
FIG. 31F shows the time course of temperature measured at position P46 in the water bath. FIG.
It is a figure which shows the temperature distribution in the horizontal cross section which passes through the center of the water tank of FIG. 25 (after 1 minute from the inflow water temperature change).
FIG. 32B is a diagram showing the temperature distribution in the horizontal cross section passing through the center of the water tank of FIG. 25 (after 2 minutes from the influent temperature change).
FIG. 32C is a diagram showing the temperature distribution in the horizontal cross section passing through the center of the water tank of FIG. 25 (after 3 minutes from the influent water temperature change).
FIG. 32D is a diagram showing the temperature distribution in the horizontal cross section passing through the center of the water tank of FIG. 25 (after 4 minutes from the influent water temperature change).
FIG. 32E is a diagram showing the temperature distribution in the horizontal cross section passing through the center of the water tank of FIG. 25 (after 5 minutes from the influent water temperature change).
FIG. 32F is a diagram showing the temperature distribution in the horizontal cross section passing through the center of the water tank of FIG. 25 (after 6 minutes from the influent water temperature change).
33 is a schematic view showing the inside of a hot-dip galvanizing port to be applied in Embodiment 4 from the side.
34 is a block diagram showing the functional configuration of the temperature estimating apparatus according to the fourth embodiment.
35 is a diagram illustrating the estimation result in the fourth embodiment.
36 is a perspective view schematically showing a configuration of a tundish to be applied in Embodiment 5;
FIG. 37 is a diagram illustrating a mounting position of a thermocouple installed in the tundish of the fifth embodiment; FIG.
38 is a block diagram showing the functional configuration of the temperature estimation device according to the fifth embodiment.
FIG. 39 is a diagram illustrating the estimation result in the fifth embodiment. FIG.

(Mode for carrying out the invention)

EMBODIMENT OF THE INVENTION Hereinafter, with reference to drawings, the form for implementing the temperature estimation method of the fluid system of this invention, the temperature distribution estimation method of a fluid system, the temperature distribution monitoring method of a fluid system, and a temperature estimation apparatus is demonstrated. In addition, this invention is not limited by this embodiment. In addition, in description of drawing, the same code | symbol is attached | subjected to the same part.

[Concept of the present invention]

1 is a functional block diagram for explaining the concept of the present invention. As shown in FIG. 1, this invention estimates the temperature in the arbitrary temperature estimation point in the fluid with two or more temperature known areas. In detail, this invention uses each coordinate in a temperature estimation point using the coordinate which is the positional information of two or more temperature known areas, and the information about the flow field of the fluid system which shows the flow of the fluid in the whole fluid system. Acquire information on the forces in the known area and estimate the temperature at any temperature estimated point in the fluid using the temperature measured value (base temperature) of each temperature known area and the information on the forces in the temperature estimation point. do. The forces in each temperature base region are fluids which have flowed from the temperature base region to the advective diffusion phenomenon due to the flow field or the inverted flow field among the total fluids at the temperature estimation point, and in the temperature base region It means the proportion of fluid (contribution rate) that has reached the temperature estimate without passing through other temperature known areas.

[Configuration of Temperature Estimation Device]

2 is a block diagram showing an example of an apparatus configuration for implementing the present invention. The temperature estimating apparatus 1 shown in FIG. 2 is connected with the 1 or more temperature measuring apparatus 2 provided in the predetermined temperature measurement site | part in the fluid system of the estimation object which estimates temperature. The temperature estimating apparatus 1 includes various IC memories such as a ROM and a RAM such as a CPU and a flash memory, a hard disk, a storage device such as various storage media, a communication device, an output device such as a display device or a printing device, an input device, or the like. It can be realized by a well-known hardware configuration provided with, for example, a general purpose computer such as a workstation or a computer can be used.

This temperature estimating apparatus 1 makes a fluid system containing a temperature known area | region into at least two places, and estimates the temperature in the predetermined temperature estimation point in a fluid system based on the temperature (base temperature) of a temperature known area | region. Estimate The position and number of temperature estimate points can be set suitably. As a typical temperature known area | region, there exists a temperature measuring area | region where the temperature measurement apparatus was arrange | positioned and the temperature was measured directly. At that time, the temperature measured value becomes a known temperature. 2 is a device configuration diagram when a temperature measuring device is used. Although the temperature measurement area | region was made into the temperature known area | region in the said description, as long as temperature is known, any area may be sufficient and it is not limited to a temperature measurement area | region.

When the temperature field is regarded as almost steady, the instantaneous value of the temperature near the time point at which temperature estimation is performed may be referred to as known temperature. If the temperature field can change in time, time series temperature data in which the time observed as the known temperature is stored in order is required, so that the temperature and time in the temperature known area are stored in time series in correspondence with each other. The temperature of the temperature known region at any time may be extracted, interpolated or extrapolated as appropriate, and output. When the temperature can be obtained by some means, such as when the temperature is fixed, for example, stored in advance in the storage device, this temperature is also treated equally to the temperature in the temperature known area, and the temperature and time are dealt with. Even if it is preserve | saved in time series, what is necessary is just to be able to extract, interpolate, or extrapolate and output the value of temperature in arbitrary time suitably.

The temperature known region can be considered to be classified into a temperature measuring region, a heat generating and absorbing region, and an inflow and outflow region by the thermofluid characteristics of the region. The temperature measurement region refers to a region including a temperature measurement site or a temperature measurement site and its vicinity. The endothermic region refers to an endothermic region or an area including the vicinity of the endothermic region. The inflow and outflow area refers to an inflow and outflow area, or an area including the inflow and outflow area and its vicinity. The temperature measurement part is a point in a fluid system in which the temperature is known by means such as actually measuring the temperature, and in which the fluid flows in and out of the fluid, or heat generation and endotherm does not occur in the area including the area or its vicinity. Or say area. The temperature measurement site is not necessarily limited to the site where the temperature is directly measured. For example, the site | part where temperature is indirectly known, such as the site | part which can convert temperature from another parameter by a model formula, etc., the site which controls temperature with a control apparatus, etc. are also included. The position and number of the temperature measurement site | part in a fluid system, for example, the installation position and number of the temperature measuring device 2 can be set suitably. An endothermic end part means the point, surface, or area | region in the fluid system in which exotherm or endotherm generate | occur | produces. An inflow-out part means the point, surface, or area | region where the inflow of the fluid into a system or the outflow of the fluid to the outside of a system generate | occur | produces. Regarding the heat-repelling portion and the inflow / exit portion, the temperature may also include an unknown portion, whereby it is possible to improve the reliability of the temperature estimation.

[Estimate principle of temperature]

First, the principle of estimation of the temperature at the temperature estimation point will be described. In addition, in the following description, the fluid system of estimation object is K temperature measurement site | part i (i = 1-K), and the heat absorption site | part i (i = K + 1-K + L) which is L heat generation site | part or an endothermic site. And an inflow / outflow site i (i = K + L + 1 to K + L + M), which is an M inflow site or an outflow site, and the temperature at the N temperature estimation point j (j = 1 to N) set in the fluid system. It is assumed to be.

The temperature estimating apparatus 1 of this embodiment uses the flow field of a fluid system, and the force of the fluid in the temperature estimation point j, specifically, the temperature measurement part i of the fluid in the temperature estimation point j, The forces of the endothermic end portion i and the inflow and outflow portion i are acquired. And the temperature estimating apparatus 1 estimates the temperature of the temperature estimation point j which considered the advection spread of the heat in the flow field by using the force of this fluid as an index value. In the following, the temperature estimation device (1), a force of the above-described, obtained downstream-side forces R _1ij and upstream forces two kinds of powers of _{R 2ij (R 1ij, R 2ij} ) , and used as an indicator value of the temperature presumption .

With reference to FIGS. 3A-3G, the concept of the downstream force and the upstream force is demonstrated. 3A is an example of a fluid system, and includes a fluid 100, a container 101, and a separator plate 102. It is assumed that there is a flow field F in the fluid system shown in FIG. 3A, and the flow circulates in the direction of the arrow in the figure. For simplicity, the fluid system of FIG. 3A does not generate or lose fluid due to inflow from the outside, outflow to the outside, chemical reactions, or the like. In the fluid system shown in FIG. 3A, there are three temperature known regions R1, R2, and R3, and are shown in circles. At this time, the definition of the force on the downstream side of the temperature known region R1 at any temperature estimation point in the fluid system is as follows.

That is, the downstream force of the temperature base region R1 at the temperature estimation point is the fluid which flows from the temperature base region R1 by the flow diffusion phenomenon by the flow field F among all the fluids at the temperature estimation point. It is defined as the ratio of the fluid which reaches from the known area R1 to another temperature known area, in this example, to the temperature estimation point without passing through the temperature known area R2 and the temperature known area R3. By this definition, the downstream force in the temperature known region R1 can be calculated for all the places in the fluid system. Similarly, the downstream forces in the temperature known region R2 and the downstream forces in the temperature known region R3 can also be calculated. As a result, the downstream force distributions I11, I12, and I13 for the temperature known regions R1, R2, and R3 are shown in Figs. 3B, 3C, and 3D, respectively. As is apparent from FIGS. 3B, 3C, and 3D, the downstream force distributions I11, I12, and I13 for the temperature known regions R1, R2, and R3 follow the flow from each region of the temperature known regions R1, R2, R3. It is a distribution extending downstream. Moreover, if there exists another temperature-known area | region on the way, it will become the distribution of the form which avoided the area | region. Since the regions indicated by the downstream side force distributions of the temperature base zones R1, R2, and R3 are areas in which a large amount of fluid flows from the temperature base zones R1, R2, and R3, respectively, the known temperatures of the temperature base zones R1, R2, and R3 are respectively. Has a strong temperature correlation with

On the other hand, the definition of the upstream force in the temperature known region R1 at the temperature estimation point is as follows. First, with respect to the flow field F, the magnitude of the flow velocity vector is the same, and a flow field in which only the directions are inverted (in the present specification, called an inversion flow field) is obtained. Then, among the entire fluids at the temperature estimation point, the fluid flows from the temperature known region R1 in accordance with the advection diffusion phenomenon due to the inversion flow field, and is different from the temperature known region R1 in the temperature known region, in this example, the temperature known region R2 and the temperature. The ratio of the fluid which reached the estimated point without passing through the known area R3 is defined as the upstream force of the temperature known area R1 at the temperature estimated point. By the above definition, the upstream force of the temperature known region R1 can be calculated for all places in the fluid system. Similarly, the upstream side force of the temperature base region R2 and the upstream side force of the temperature base region R3 can be calculated. As a result, the upstream force distributions I21, I22, and I23 for the temperature known regions R1, R2, and R3 are shown in Figs. 3E, 3F, and 3G, respectively. As is apparent from FIGS. 3E, 3F, and 3G, the upstream force distribution for the temperature known regions R1, R2, and R3 is from the respective regions of the temperature known regions R1, R2, and R3 to the upstream side opposite to the flow. It becomes an extended distribution. Moreover, if there exists another temperature-known area | region on the way, it will become the distribution of the form which avoided the area | region. Since most of the fluid in the region indicated by the upstream side force distributions of the temperature base zones R1, R2, and R3 flows into the temperature base zones R1, R2, and R3, respectively, the upstream side force distribution of the temperature base zones R1, R2, and R3. There is a strong temperature correlation between the temperature of the region indicated by and the known temperature of the temperature known regions R1, R2, and R3.

In order to acquire the downstream side force _R1ij and the upstream side force _R2ij , in this embodiment, in each part i of the K temperature measurement site | part i, the L end part endothermic site | part, and the M inflow-out part i, On the other hand, each part i is included and finite area | region i (i = 1-K + L + M) which does not overlap with each other is set. Specifically, the temperature measurement region i is set as the region i corresponding to the temperature measured portion i, the heat absorbing region i is set as the region i corresponding to the heat absorbing portion i, and the region i corresponding to the inflow and outflow portion i. Set the inflow and outflow area i. The shape of the temperature measurement area | region i, the heat absorption area | region i, and the inflow-out area | region i which are set may be any shape as long as it contains the temperature measurement site | part i, the heat absorption area | region i, or the inflow-outflow part i. That is, each area i may be a point, a line, or a surface, for example, and may be an area having a three-dimensional finite volume.

In order to be applicable to a wide range of objects, the temperature measuring region i, the heat absorbing region i, and the temperature measuring region i, the heat absorbing region i, and the inflow-out region i corresponding to the temperature measuring region i, the endothermic region i, and the inflow-out region i using the following method It is good to set. For example, the temperature measurement site | part i is an installation position of the temperature measuring device 2 provided in a fluid system. For this reason, when the temperature measurement area | region i sets the spherical region of radius r centering on the installation position of the temperature measuring device 2 which is this temperature measurement site | part i, and sets it as the temperature measurement area | region i, good. If the radius r of the sphere area to be set is a large value, the estimated temperature distribution of the fluid system becomes a steep temperature distribution, and a small value results in a smoothed temperature distribution. The value of the specific radius r differs from an optimum value according to the flow characteristic of a fluid system. For example, when 1 m of water tank is applied to one side and 6 temperature measuring apparatuses 2 are installed in this tank in 6 places, and 6 temperature measuring site | parts i are installed, the temperature measuring area i is, for example, It is preferable to set it as the sphere area | region of radius r = 0.05m. The radius r of this sphere region may be any value as long as the region i of the temperature measurement region i, the heat absorbing region i, and the inflow-out region i do not overlap, but is set for each temperature measurement region i, respectively. It is preferable that all the radii of the temperature measurement region i are the same.

When an area in which the heat generation or endotherm occurs due to heating by a heating device, endotherm by a heat absorbing device, a chemical reaction, or the like is included in the fluid system, the area becomes an endothermic portion i. In this case, this area is referred to as the heat absorption region i. For example, heat generation and endotherm cause a step of the fluid system, specifically, a wall surface such as a facility for defining a flow of the fluid system to be estimated, or a bath surface of the fluid system to be estimated ( When generated on a bath surface, this wall surface or bath surface is the heat absorption region i. For example, when a solid immersed in a fluid system generates heat or endothermic, when a substance which generates heat or endotherm due to a chemical reaction or the like is immersed in the fluid system, the surface of the solid is absorbed by the heat-absorbing region i. Shall be. When heat generation or endotherm occurs in a part of the fluid system, for example, when an induction heating apparatus is installed in the fluid system, the area in the fluid system to which heating energy is applied is referred to as the heat absorption region i.

If there is an inflow of fluid into the system or an outflow of fluid to the outside of the system, the inflow region or the outflow region becomes the inflow / outflow region i. In this case, this area is referred to as the inflow / outflow area i. For example, when fluid flows in from the boundary surface which divides the flow of the fluid system of estimation object, or fluid flows out from this interface surface, the said boundary surface is made into the inflow-out area | region i.

However, when temperature measurement area | region i, heat absorption area | region i, and area | region i and i 'made into inflow-out area | region i and i' (i = 1-K + L + M, i '= 1-K + L + M, i ≠ i') overlap, it is a downstream side. The forces R _1ij and the upstream forces R _2ij cannot be obtained. For this reason, the temperature measurement area i, the heat absorption area i, and so that each area i necessarily belongs to any one of the temperature measurement area i, the heat absorption region i, and the inflow-out area i, that is, each region i does not overlap. In this case, it is necessary to determine the shape and size of the inflow and outflow area i.

A fluid that has passed through one region of interest (notice region) i or that has been created within the region of interest i and has reached a temperature estimate j without passing through another region i ' It is defined as the fluid component of the region of interest i in FIG. Then, the ratio of the fluid component of the region of interest i to the total fluid in the temperature estimation point j is defined as the force of the corresponding region i in the temperature estimation point j, and the flow field of the fluid system (hereinafter, "actually Force obtained by using the downstream force R _1ij of the corresponding part i in the temperature estimation point j, and the force acquired by using the inverse flow field of the fluid system at the temperature estimation point j. It is defined as the upstream force R _2ij of the corresponding site i.

That is, downstream force _R1ij is acquired using the flow field (henceforth "real flow field") of a fluid system. The actual flow field is calculated using, for example, numerical simulations, actual equipment, experimental apparatus that simulates the actual equipment, and the like. For example, a flow velocity vector of the entire fluidic system to be estimated, specifically, a flow velocity vector representing a direction and a flow rate of the fluid in each region where the entire fluidic system is divided in the same size, is obtained to be an actual flow field. .

Using this actual flow field, the proportion of the fluid component at the temperature estimation point j is calculated for all the regions i, and is obtained as the downstream force R _1ij . Specifically, for each of the regions of the temperature measurement region i, the heat absorbing region i, and the inflow-out region i, the fluid has passed through the corresponding region i or generated in the region i, and the other temperature measurement region. The ratio of the fluid (fluid component) which reached to the temperature estimation point j without passing through i ', the heat absorption zone i' and the inflow and exit area i 'with respect to the total fluid of the temperature estimation point j, _Let it be the force R _1ij .

On the other hand, the upstream force R _2ij is obtained using the inversion flow field of the fluid system. This inversion flow field is obtained by inverting all the velocity vectors obtained as the actual flow fields. By using this inversion flow field, the ratio of the fluid component in the temperature estimation point j is calculated in all the regions i, and is obtained as the upstream force R _2ij . Specifically, for each region i of the temperature measuring region i and the inflow / outflow region i, it is a fluid which has passed through the corresponding region i or is generated in this region i, and the other temperature measuring region i ', the heat absorption region The ratio of the fluid (fluid component) to the total fluid of the temperature estimation point j of the fluid (fluid component) which has reached the temperature estimation point j without passing through i 'or the inflow / exit area i' is calculated to be the upstream force R _2ij . In addition, about the heat absorption area | region i, the ratio of a fluid component is made into "0", and it is set to the upstream force _R2ij (it is set to the upstream force _R2ij = 0). This is because there is a temperature correlation between the endothermic end portion i and the position downstream of the endothermic end portion i, but there is no temperature correlation between the endothermic end portion i and the position upstream of the endothermic end portion i.

Then, the downstream-side forces obtained as described above R _1ij and upstream forces R _2ij of the pair forces (R _1ij, R _2ij) on the basis of, and even a weight W _ij, the details of each portion i of the temperature estimation point j The weight W _ij for weighting the known temperature of each part i is calculated. For example, the weight W _ij is calculated as W _ij = W (R _1ij , R _2ij ) using a weighting function W (R ₁ , R ₂ ) which becomes a monotone non-decreasing function. .

In addition, in Non-Patent Document 2, the force ranges of the outlet and the suction port are defined. This is similar to the downstream and upstream forces described in the present invention, but with a different concept. In other words, the force ranges of the outlet and the suction port of the non-patent document 2 are a method that can be applied only to the inflow site and the outflow site of the fluid system, that is, the site that can be set as the boundary condition. Therefore, the force range of the blowout port and the suction port cannot be defined for the measured portion and the heat absorbing portion. On the other hand, the downstream force and the upstream force newly devised in the present invention can be defined not only for the boundary but also for the measured portion and the endothermic portion existing inside the fluid system. In the present invention, it is very important to consider the measured portion or the endothermic portion, and it is essential to use the concepts of the downstream force and the upstream force newly devised in the present invention.

If the flow sheets regarded as almost normal for, and then, downstream of the forces obtained as described above R _1ij and upstream force pair, the power of R _2ij based on the (R _1ij, R _2ij), a temperature estimation part for the point j the weight W _ij, i Details of each calculates a weight W _ij for performing weighting for the known temperature of the parts i. For example, the weight W _ij is calculated as W _ij = W (R _1ij , R _2ij ) using the weighting functions W (R _1ij , R _2ij ) _serving as the monotonic non-reducing function.

Then, by the weighted average using the known temperature T _i of each site i of the temperature measurement site i, the endothermic heat site i, and the inflow and outflow site i, and the weight W _ij for each site i with respect to the calculated temperature estimation point j , The estimated temperature at the temperature estimation point j is calculated. Base temperature of each part of i T _i, and a temperature estimation estimates the temperature of the point j with the relationship Te _j is indicated by the food (2). Temperature base temperature of the measured area i T _i is a temperature actually measured value that is measured by the temperature measuring device (2). The known temperature T _i of the heat-absorbing end region i and the inflow-out region i uses the value when the temperature of the corresponding region i is known. In the case where the temperature of the endothermic end portion i or the inflow and outflow portion i is unknown, the temperature estimate point j is estimated according to the following equation (2) after replacing the value of the weight W _ij for the corresponding portion i with "0". The temperature Te _j is calculated.

 [Number 2]

On the other hand, the present invention can be used even when the flow field is considered to be almost normal, in which case the temperature field may change in time. In this case, the transfer time mentioned later together with the said force is acquired, and temperature is estimated using time series temperature data. The time series temperature data is a sequence of the measured actual temperature i, the endothermic region i, and the known actual temperature observed in each of the inlet or outlet regions i, or the known temperature determined by other means, and the observed time in order. Can be output by interpolating and extrapolating the temperature T _i (t) at an arbitrary time t from the recording of the measured temperature and time. If the measured temperature is not observed at time t due to a measuring instrument failure or the like, the temperature may be treated as unknown, but when the measured temperature is observed at another time, interpolation and extrapolation of the data in the vicinity One temperature may be output as T _i (t).

The transfer time is the time required for the fluid to travel by advection diffusion between each site i and the temperature estimate point j. Specifically, the delivery time consists of the downstream delivery time and the upstream delivery time, and the time required for the fluid to move from each site i to the temperature estimate point j is the downstream delivery time τ _{1 ij} , the fluid is the temperature The time required to move from the estimated point j to each site i is the upstream propagation time τ _{2 ij} . Then, using time-series temperature data and time t ₀ for estimating temperature and propagation time τ _{1 ij} , τ _{2 ij} , the temperature at a point in time in the past or future from the time for which temperature is to be estimated by the transfer time minute Calculate and let it be known temperature. Specifically, based on the time t ₀ to estimate the temperature, the temperature of the site i observed in the past by the downstream transmission time τ _{1 ij} minutes at a past point in time is the downstream known location of the site i at the temperature estimation point j. Let it be temperature. In other words, the actual output temperature _{T i (t 0 -τ 1 ij} ) at the time t ₀ -τ _{1 ij} from the time-series temperature data, may be to a downstream base temperature. Similarly, the temperature of the site i observed at a future time point by the upstream delivery time τ _{2 ij} minutes on the basis of the time t ₀ for which temperature is to be estimated is taken as the upstream known temperature of the site i at the temperature estimation point j. . In other words, the actual output temperature _{T i (t 0 + τ 2} ij) at the time t ₀ + τ _{2 ij} from the time-series temperature data, may be to the upstream base temperature.

Downstream forces R _1ij and upstream forces R _2ij pair forces (R _1ij, R _2ij) on the basis of the downstream-side weight W downstream base temperature of _1ij and upstream weights W _2ij, specifically, each portion i T _The weights W _1ij and W _2ij respectively corresponding to _i (t ₀ -τ _{1 ij} ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ) are calculated. The downstream side for the weights W _1ij weights W _1ij using a minor non-decreasing function of the downstream side forces R _1ij weighting functions _{W 1 (R 1ij, R 2ij} ) that for any upstream forces R _2ij W _1ij = W ₁ (R _1ij, R _2ij), the upstream-side weights for W _2ij weighting function where the upstream forces minor non-decreasing function of R ₂ with respect to any of the downstream side forces R _1ij W ₂ and (R _1ij, R _2ij) using a it may be calculated by a weight W _2ij in _{W 2ij = W 2 (R 1ij} , R 2ij).

Downstream weights W _{1ij for} each site i with respect to the temperature measurement site i, the endothermic area i, and the downstream known temperature T _i (t ₀ -τ _{1 ij} ) of each site i of the inflow and outflow site i and the temperature estimation point j and calculates the upstream-side base temperature _{t i (t 0 + τ 2} ij) and a temperature estimation estimates the temperature of the point j in the, time t ₀ by a weighted average using a weight W upstream _2ij. The estimated temperature Te _j (t ₀ ) of the temperature estimation point j at time t ₀ is represented by the following equation (3).

[Number 3]

The downstream known temperature T _i (t ₀ -τ _{1 ij} ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ) of the heat absorption end region i and the inflow and outflow region i are known when the temperature of the corresponding region i is known. Use that value in. When the to the heat absorbing portion temperature i or inflow output area i of the unknown is, after replacing the downstream weight W _1ij the upstream value of the weight W _2ij for the area i to "0", the formula (3) The estimated temperature Te _j (t ₀ ) of the temperature estimation point j is calculated accordingly.

Next, the process sequence which the temperature estimation apparatus 1 performs is demonstrated with reference to FIGS. When the temperature field of the target fluid system is regarded as almost normal, the processing sequence of FIGS. 4 to 7 is used, and when the temperature field of the target fluid system can change over time, the processing sequence of FIGS. 8 to 12 is used. . The temperature estimating apparatus 1 performs a process according to the processing sequence shown in FIGS. 4-7 or 8-12, and the temperature estimation method of a fluid system, the temperature distribution estimation method of a fluid system, and the temperature distribution of a fluid system Implement the monitoring method. The processing described herein can be realized by storing a program for realizing this processing in a storage device of the temperature estimating apparatus 1, for example, by reading and executing the program.

First, the processing sequence in the case where the temperature field of the target fluid system is regarded as almost normal will be described with reference to FIGS. 4 to 7. First, the procedure of the process (weight value calculation process) which the temperature estimation apparatus 1 performs in order to calculate said weight W _ij is demonstrated. 4-6 is a flowchart which shows the processing sequence of a weight calculation process. Here, for the method for obtaining the forces (R _1ij, R _2ij) by the temperature distribution analysis using fluid numerical simulation example illustrates a processing procedure of the case of calculating the weight W _ij.

In the weight calculation process shown in FIG. 4 and FIG. 5, first, as shown in FIG. 4, as a flow field acquisition process, typical flow field calculation conditions of the fluid system which are estimation objects are set using a numerical fluid simulation (step S1). ), The normal flow field is calculated based on the set flow field calculation conditions, and is set as the actual flow field (step S3). Here, calculation of the flow field of a fluid is performed using a well-known technique. Specifically, any fluid analysis solver capable of obtaining the flow and temperature fields of the fluid may be used, including a commercially available product. For example, an ANSYS FLUENT (registered trademark) or the like may be used. Calculate the flow field of. Moreover, although the method of calculating a two-dimensional flow field and the method of calculating a three-dimensional flow field conventionally are known, according to the characteristic of the fluid system to be estimated, a two-dimensional flow field or a three-dimensional flow field is known. It is good also as a method of selecting and using suitably.

Subsequently, as the region setting step, the temperature measurement region i, the endothermic region i and the temperature measurement region i, the endothermic region i and the inflow-out region i corresponding to the respective portions i of the temperature measurement region i, the endothermic region i in the fluid system i i (i = 1 to K + L + M) is set (step S5). Moreover, as an estimation point setting process, the temperature estimation point j (j = 1-N) is set in a fluid system (step S7). Next, the site | part i which acquires the downstream side force _R1ij is designated from the temperature measurement site | part i, the heat absorption part i, and each site | part i of the inflow-out | portion part i (step S9). The processing here can be realized by increasing the value of i sequentially in the range of 1 to K + L + M for each iteration of steps S9 to S23.

Next, the temperature estimation point j is specified (step S11). The processing here can be realized by increasing the value of j in order in the range of 1 to N every time the steps S11 to S21 are repeated. Subsequently, as a downstream force acquiring step, first, a boundary condition necessary for a numerical fluid simulation is given, but since the force is calculated by the temperature distribution analysis, the value of the temperature in each region i is used as the boundary condition. To give. Specifically, the boundary condition which fixes the temperature of the area | region i corresponding to the designated site | part i to "1" is provided, and the boundary condition which fixes the temperature of another area | region i '(i ≠ i') to "0" is provided. (Step S13). Thereafter, a numerical fluid simulation of the advection diffusion phenomenon is carried out using the actual flow field, and the temperature distribution analysis is performed under the given boundary condition. Specifically, the normal temperature distribution is calculated (step S15), and the temperature value at the temperature estimation point j is obtained according to the obtained normal temperature distribution (step S17). This temperature value corresponds to the ratio of the fluid component of the area | region i in the temperature estimation point j in an actual flow field. And the acquired temperature value is made into the value of the downstream force _R1ij of the designated part i in temperature estimation point j (step S19).

Subsequently, it is determined whether or not the downstream force R _1ij is acquired for all the temperature estimated values j. When _there exists the temperature estimated value j whose downstream side force _R1ij is not acquired (step S21: No), it returns to step S11 and repeats the above process. Subsequently, it is determined whether or not the downstream side force R _1ij is obtained for the temperature actual measurement site i, the heat absorbing heat site i, and all the sites i of the inflow and outflow site i. If there is a portion i in which the downstream force R _1ij is not acquired (step S23: No), the process returns to step S9 and the above-described processing is repeated. If downstream force R _1ij is acquired about all the parts i (step S23: Yes), then, as shown in FIG. It calculates as a chapter (step S25). And the site | part i which acquires an upstream side force _R2ij is specified from the temperature measurement site | part i, the heat absorption part i, and each site | part i of the inflow-out | portion part i (step S27). In the same manner as in step S9, the value of i may be increased in order in the range of 1 to K + L + M for each iteration of steps S27 to S45.

Next, the temperature estimation point j which acquires the upstream force _R2ij is specified (step S29). Similarly to step S11, every time the steps S29 to S43 are repeated, the value of j may be increased in the order of 1 to N in order. Subsequently, the treatment branches for the case where the designated portion i is the temperature measurement region i or the inflow-out region i and the case where the endothermic heat region i is used. That is, when the designated part i is the temperature measured part i or the inflow-out part i (step S31: Yes), as an upstream force acquisition process, first, the area i corresponding to this designated temperature measured part i or the inflow-out part i The boundary condition which fixes the temperature of to is set to "1", and the boundary condition which fixes the temperature of another area | region i '(i ≠ i') to "0" is given (step S33). Then, the numerical fluid simulation is performed using the inverted flow field, and the temperature distribution analysis is performed under the given boundary condition. Specifically, the normal temperature distribution is calculated (step S35), and the temperature value at the temperature estimation point j is obtained according to the obtained normal temperature distribution (step S37). This temperature value corresponds to the ratio of the fluid component of the area | region i in the temperature estimation point j in an inversion flow field. And the acquired temperature value is made into the value of the upstream-side force _R2ij of the designated site | part i in the temperature estimation point j (step S39), and it transfers to step S43 after that.

On the other hand, when the designated part i is not the temperature actual measured part i or the inflow / exited part i but the heat absorbing end part i (step S31: No), the value of the upstream force R _2ij of the designated part i is set to "0" ( Step S41), the flow then advances to step S43. In step S43, it is determined whether or not the upstream force R _2ij is acquired for all the temperature estimation points j. If there is a temperature estimation point j in which the upstream side force R _2ij is not acquired (step S43: No), the process returns to step S29 and the above-described processing is repeated. If the upstream forces R _2ij are acquired for all the temperature estimation points j (step S43: Yes), it is determined whether or not the upstream forces R _2ij are acquired for all the sites i. If there is a part i whose upstream force R _2ij is not acquired (step S45: No), the process returns to step S27 and the above-described processing is repeated.

The acquisition method of the downstream force R _1ij and upstream forces R _2ij, the whole with respect to the specified temperature estimation point j, but that acquires downstream forces R _1ij and upstream forces R _2ij, the present method is readily fluid zone It is possible to obtain the downstream force distribution and the upstream force distribution of the entire fluid region. Specifically, the temperature estimate points j are placed at positions of all calculation grids j 'of the arrangement, for example, numerical fluid simulations, that cover the fluid region sufficiently thoroughly, and the downstream forces R for all the temperature estimate points j. _{By obtaining 1ij} and the upstream side force R _2ij , the downstream side force distribution and the upstream side force distribution of the entire fluid region can be obtained.

If upstream side force _R2ij was acquired about all the site | parts i (step S45: Yes), then, as shown in FIG. 6, the site | part i which calculates the weight W _ij first from each site | part i as a weight calculation process first. (Step S47). In the same manner as in step S9, the value of i may be increased in order in the range of 1 to K + L + M for each iteration of steps S47 to S59. Next, the temperature estimation point j which calculates weight W _ij is specified (step S49). In the same manner as in step S11, the value of j may be increased in order in the range of 1 to N for each iteration of steps S49 to S57.

Next, it is determined whether the temperature of the designated part i is known or not. In general, when the designated part i is a temperature measured part i, the temperature is known. However, although the temperature observation may be temporarily unavailable due to a malfunction of the measuring instrument, in this case, the temperature of the measured portion i may be unknown. On the other hand, temperature may be unknown about the heat absorption part i or inflow-out part i. For this reason, when the temperature of the designated part i is known (step S51: Yes), the temperature is determined using the weighting function W (R _1ij , R _2ij ) based on the forces R _1ij and R _2ij of the designated part i. The weight W _ij of the designated part i in the estimation point j is calculated (step S53).

The weighting function W (R _1ij, R _2ij) is, and is a monotone non-decreasing function of R _1ij for any R _2ij As shown in the food 4, the forging ratio decrease in R _2ij for any R _1ij If the function is to be a function and the downstream force R _1ij and the upstream force R _2ij are both "0", that is, a function that becomes "0" when R _1ij = R _2ij = 0, This can also be applied.

 [4]

The weighting functions W (R _1ij , R _2ij ) _{vary depending on the} spatial scale, the flow rate scale, the interval between the temperature measurement sites, and the like, but are relatively simple and can be widely used even when a certain fluid system is to be estimated. a weighting function, may be mentioned the food (5) downstream of the forces R _1ij and upstream power weighting function for calculating an average value for each portion i of _{R 2ij W (R 1ij, R} 2ij) as shown in. For example, the weighting functions W (R _1ij , R _2ij ) shown in the following _equation (5) are not able to grasp the presence or absence of the endothermic end portion i in the fluid system to be estimated, or Although the endothermic part i is included, it is suitable for the case where the exact position cannot be understood.

 [Number 5]

The temperature of all the endothermic end portions i and the inlet and outlet regions i is known, or the temperature measurement region i exists in the vicinity of all the endothermic end portions i and the inlet and outlet portion i (all endothermic regions i and the inlet and outlet regions i Temperature measurement site i is present within a predetermined predetermined distance range of and < RTI ID = 0.0 > and < / RTI > _{Since the precision} of the downstream force R _1ij is higher than that of _2ij , the weight function W (R _1ij , R _2ij ) shown in the following _equation (6) using only the downstream force R _1ij as the weighting function W (R _1ij , R _2ij) ) Can be used.

 [Number 6]

The temperature of all the endothermic region i and the inlet / outlet region i is known, or the temperature measurement region i exists in the vicinity of all the endothermic region i and the inlet / outlet region i, and all of the endothermic region i and the inlet / outlet region i the temperature in the case where the downstream side of the bore flow from the measured region i, downstream forces R due _1ij than the upstream forces R _2ij side is a high precision, the weight function W (R _1ij, R _2ij) as the upstream forces R _The weighting functions W (R _1ij , R _2ij ) shown in the following expression (7) using only _2ij may be used.

 [Numeral 7]

If all of the positions in the fluid system of the endothermic end region i or the inlet / outlet region i which contribute significantly to the temperature distribution are specified, and some or all of the temperature of the endothermic end region i and the inlet / outlet region i are unknown, The weighting functions W (R _1ij , R _2ij ) shown in 8) may be used. S _1j is the sum of the temperature measurement site i, the endothermic heat generating site i, and the downstream force R _1ij with respect to the inflow and outflow site i whose temperature is known, and S _2j is the temperature measurement site i and the heat absorbing site where the temperature is known The sum of i and the upstream forces R _2ij for the inflow and outflow region i.

 [Numeral 8]

The location in the fluid system of the endothermic end region i or the inlet / outlet region i, which contributes greatly to the temperature distribution, is all specified, some or all of the temperature of the endothermic end region i and the inlet / outlet region i are unknown, and the temperature estimate j In the case where both the downstream force R _1ij and the upstream force R _2ij become small values, the weighting functions W (R _1ij , R _2ij ) shown in the following expression (9) may be used. S _1j is the sum of the temperature measurement site i, the endothermic heat generating site i, and the downstream force R _1ij with respect to the inflow and outflow site i, whose temperature is known, and S _avej is the temperature measurement site i and the heat absorbing site where the temperature is known. It is the sum total of the average value (1/2) x (R _1ij + R _2ij ) between the downstream and upstream forces with respect to i and the inflow-exit site | part i.

 [Number 9]

The weighting function W (R _1ij , R _2ij ) may apply the same thing uniformly to all the temperature estimation points j set in step S7, and for each temperature estimation point j, the appropriate weighting functions W (R _1ij , R _2ij ) may be selectively used.

Returning to FIG. 6, if the weight W _ij is calculated as described above, the flow proceeds to step S57. In addition, when the temperature of the designated part i is unknown (step S51: No), the weight W _ij of the designated part i in the temperature estimation point j is set to "0" (step S55), and then the process proceeds to step S59. do. In step S57, it is determined whether the weight W _ij is computed with respect to all the temperature estimation points j. If there is a temperature estimation point j whose weight W _ij is not calculated (step S57: No), the process returns to step S49 and the above-described processing is repeated. And if the weight W _ij is computed with respect to all the temperature estimate points j (step S57: Yes), it transfers to step S59.

In step S59, it is determined whether or not the weight W _ij has been calculated for all the portions i. If there is a part i whose weight W _ij is not calculated (step S59: No), the process returns to step S47 and the above-described processing is repeated. And if the weight W _ij is computed about all the site | parts i (step S59: Yes), the weight W _ij for each site | part i in the calculated temperature estimation point j is preserve | saved in a memory | storage device (step S61), and weight calculation process Ends. The estimated temperature Te _j of any of the temperature estimation points j in the fluid system can be estimated by calculating the weight W _ij for each site i in this temperature estimation point j in the above-described order, but the temperature distribution of the entire fluid system is estimated. In order to visualize, it is necessary to set the temperature estimation point j over the whole fluid system and calculate the weight W _ij for all the set temperature estimation points j. In this case, it is preferable to calculate the weight W _ij for each site i for all the temperature estimation points j in advance, and save it as a database (weight database) in the storage device.

Next, the procedure of the process (temperature estimation process) for estimating the temperature of arbitrary temperature estimation point j using the weight W _ij for each site | part i with respect to the temperature estimation point j computed as mentioned above is demonstrated. 7 is a flowchart showing a processing procedure of a temperature estimation process.

In a temperature estimation process, as shown in FIG. 7, first, the temperature measurement part i, the heat absorption part i, and the known temperature T _i of the inflow-and-exit part _i are acquired (step S71). For the actual measurement temperature region i, actually measured value of the temperature input from the temperature measured area i temperature measuring device (2) installed on the base is obtained as a temperature T _i. About the heat-absorbing end portion i and the inflow / exit portion i, a temperature measuring device is installed in the corresponding portion i, and the temperature is measured. If it can be acquired by any means, such as when it is stored in advance in the device, it is acquired.

Next, the temperature estimation point j which estimates temperature is specified (step S73). The processing here can be realized by increasing the value of j in order in the range of 1 to N every time the steps S73 to S81 are repeated. Subsequently, the weight W _ij for each site i with respect to the designated temperature estimation point j is read and obtained from the storage device (step S75). For example, the weight W _ij for the designated temperature estimation point j is obtained from the above weight database. And, according to the above equation (2), subjected to a weighted average process using the weight W _ij obtained at the base temperature T _i and the step S75 of the parts i obtained in step S71, the estimated temperature Te _j of the temperature estimation point j It calculates (step S77). Thereafter, the estimated temperature Te _j of the calculated temperature estimation point _j is stored in the storage device (step S79).

Thereafter, it is determined whether or not the estimated temperature Te _j is calculated for all the temperature estimation points j. If there is a temperature estimation point j whose estimated temperature Te _j is not calculated (step S81: No), the process returns to step S73 and the above-described processing is repeated. On the other hand, if the calculated estimated temperature Te _j for all temperature estimation point j (step S81: Yes), the end temperature estimation process.

Next, the processing sequence in the case where the temperature field of the target fluid system can fluctuate over time will be described with reference to FIGS. 8 to 12. In the case where the temperature field can fluctuate, in addition to the forces, acquisition of time-series temperature data, calculation of propagation time, calculation of downstream weight and upstream weight, calculation of downstream known temperature and upstream known temperature are necessary. .

First, the calculation of the propagation time, that is, the downstream propagation time τ _{1 ij} and the upstream propagation time τ _{2 ij} , will be described. 8 is a flowchart showing a processing procedure of a delivery time calculating process. The downstream delivery time τ _{1 ij} means the time required for the fluid to move from the temperature measurement site i, the endothermic site i, and the inlet and outlet site i to the temperature estimate point j by the advection diffusion, and the upstream delivery The time τ _{2 ij} means the time required for the fluid to move from the temperature estimation point j to the temperature measurement region i, the endothermic region i, and the inflow and outflow region i. τ _{1 ij} is the time required for the fluid to move in the direction of the temperature estimate point j on the downstream side of the flow from the temperature measurement site i, the endothermic site i, and the inflow and outflow site i, so that the downstream delivery time Similarly, since τ _{2 ij} is the time required for the fluid to move from the direction of the temperature estimation point j on the upstream side of the flow as seen from the temperature measurement site i, the endothermic site i, and the inflow and outflow site i This is called upstream propagation time. Hereinafter, the pair (τ _1ij , τ _{2 ij} ) of the downstream transfer time and the upstream transfer time is called a transfer time.

Hereinafter, as an example of the transfer time calculation method, a calculation method of the transfer time τ _{1 ij} , τ _{2 ij} using a numerical fluid simulation of temperature will be described.

First, a typical boundary condition of a fluid system is set using a numerical fluid simulation (step S101), and then a flow field is calculated based on the set boundary condition (step S103). Since this flow field is the same as the actual flow field calculated | required in the order (step S3 of FIG. 4) when a temperature field is considered almost normal, you may use the said actual flow field as it is.

Next, after setting the temperature measurement part i of the fluid system, the heat absorption part i, and the inflow-outflow part i (i = 1-K + L + M) (step S105), the temperature estimation point j (j = 1-N) is set. (Step S107). From the set temperature measurement site i, endothermic heat site i, and inflow-outflow site i (i = 1 to K + L + M) and the temperature estimation point j (j = 1 to N), the transfer time (τ _{1 ij} , τ _{2 ij} ) is determined. The site | part i to calculate and the temperature estimation point j are specified, respectively (step S109 and step S111). The processing here increases the value of i sequentially in the range of 1 to K + L + M at the time of repetition of step S109 to step S131, and the value of j at the time of repetition of step S109 to step S129 in the range of 1 to N. This can be achieved by increasing in order from.

Then, setting the heating condition of the oil system, the initial temperature T ₀ to the total (unit K) and a given hereinafter (step S113), part i heating value S (unit W) at a position along the (step S115). Under these conditions, an abnormality calculation of the temperature distribution is performed (step S117), and the temperature rise behavior at the temperature estimation point j is calculated. When the temperature at the temperature estimation point j reaches the threshold temperature T _C (unit K), the time τ _{1 ij} taken until the temperature reaches T _C from T ₀ is calculated (step S119). τ _{1 ij} becomes the downstream propagation time. Since the initial temperature T ₀ is a value that does not affect the passing of time, may grant any value. Regarding the calorific value S and the threshold temperature T _C , the optimum value differs depending on the target fluid system. For example, in the case of a hot dip galvanizing port, a hot metal retaining furnace and a tundish, S = 2,000KW and T _C = T ₀ + 1K may be sufficient.

In the same manner, after the initial temperature T _{0 is applied} to the whole fluid system (step S121), the calorific value S is given to the position of the temperature estimation point j (step S123), and abnormal calculation of the temperature distribution is performed (step S125). the position of the temperature of the i and calculates the time τ _{2 ij} jammed until it reaches the T _C from T ₀ (step S127). τ _{2 ij} is the upstream propagation time. The propagation time τ _{1 ij} , τ _{2 ij} is the time and temperature required for the fluid to move from the temperature measurement site i, the endothermic region i, and the inflow and outflow site i to the temperature estimate point j by advection diffusion. Any indicator may be used as long as it is an indicator corresponding to the time required for the fluid to move from the estimated point j to the temperature measurement site i, the endothermic heat site i, and the inflow and outflow site i, and the definition method is not particularly limited.

In step S129, it is determined whether the propagation time has been calculated for all the temperature estimation points j. If transmission time is _(ij τ _1, τ _{2 ij)} already calculated temperature estimated point j (the step S129: No), and repeats the above processing and returns to step S111. If the transfer time has been calculated for all the temperature estimation points j (step S129: Yes), the process proceeds to step S131. In step S131, it is determined whether or not the transfer time is calculated for all sites i. If there is a portion i in which the propagation times τ _{1 ij} and τ _{2 ij} are not calculated (step S131: No), the process returns to step S109 and the above-described processing is repeated. Then, if the transfer times τ _{1 ij} and τ _{2 ij} have been calculated for all sites i (step S131: Yes), the downstream transfer time τ _{1 ij} and the upstream transfer time τ _{2 ij} for each calculated site i are calculated. It saves to a memory device (step S133), and complete | finishes a delivery time calculation process. In order to estimate and visualize the temperature distribution of the entire fluid system, the downstream transfer time τ _{1 ij} and the upstream transfer time τ _{2 ij} for each site i are calculated in advance for all the temperature estimation points j, and the data is stored in the storage device. It is preferable to save as a base.

Next, the flow field acquisition process, the area | region setting process, the downstream side force acquisition process, and the upstream side force acquisition process shown in FIG. 9 and FIG. 10 are performed. 9 and 10 are flowcharts showing the processing procedure of the calculation process of the downstream weight W _1ij and the upstream weight W _2ij . These should just be the same procedure as the said procedure (step S1-step S45 of FIG. 4 and FIG. 5) described about the case where the temperature field of the target fluid system is regarded as almost normal.

Next, the procedure of the process (weighting calculation process) which the temperature estimation apparatus 1 performs in order to calculate a weight is demonstrated. When the temperature field can fluctuate over time, the downstream weight W _1ij and / or the upstream weight W _2ij are calculated as the weight. If the downstream side force R _1ij and the upstream side force R _2ij were acquired about all the site | parts i (step S201-S245, FIG. 9 and FIG. 10), then, as shown in FIG. 11, it is downstream from each site i. The site | part i which calculates side weight _W1ij and upstream weight _W2ij is specified (step S247). Herein, the value of i may be increased in order in the range of 1 to K + L + M at each iteration.

Next, the temperature estimation point j which calculates downstream weight _W1ij and upstream weight _W2ij is specified (step S249). For each iteration, the value of j should be increased in order from 1 to N. Next, it is determined whether or not the temperature of the designated part i is known (step S251). In general, when the designated part i is a temperature measured part i, the temperature is known. However, when temperature observation is temporarily unavailable due to a malfunction of the measuring instrument or the like, the temperature of the measured portion i may be unknown. On the other hand, temperature may be unknown about the heat absorption part i or inflow-out part i. For this reason, when the temperature of the designated part i is known (step S251: Yes), the downstream weighting function W ₁ (R _1ij , R _2ij ) and the basis of the forces R _1ij and R _2ij of the designated part i and and it calculates the upstream weighting function _{W 2 (R 1ij, R 2ij} ) downstream and upstream weights W _{1ij 2ij} weight W of the designated area i of the temperature estimation point j using (step S253 and step S255).

Downstream of the weighting function _{W 1 (R 1ij, R 2ij} ) and upstream of the weighting function _{W 2 (R 1ij, R 2ij} ) is, and as shown in the food 10, the R _1ij for any R _2ij forging ratio It becomes a decreasing function and _becomes a monotonic non-decreasing function of R _2ij with respect to any R _1ij , and when both the downstream force R _1ij and the upstream force R _2ij are "0", that is, R _1ij = Any function can be applied as long as it is a function that becomes "0" when R _2ij = 0.

 [Number 10]

The downstream weighting functions W ₁ (R _1ij , R _2ij ) and the upstream weighting functions W ₂ (R _1ij , R _2ij ) change the optimal functional type depending on the spatial scale, the flow rate scale, and the interval between the temperature measurement sites. It is relatively simple, and the downstream weighting functions W ₁ (R _1ij , R _2ij ) and the upstream weighting functions W ₂ (R _1ij , R _2ij ) that can be widely used even when any fluid system is to be estimated, 11a) and downstream weighting functions W ₁ (R _1ij , R _2ij ) with 0.5 times the downstream force as the downstream weight and 0.5 times the upstream force as the upstream weight, as shown in (11b) and And an upstream weighting function W ₂ (R _1ij , R _2ij ). For example, this following formula (11) is not able to grasp the presence or absence of the endothermic end part i in the fluid system of estimation object, or although the fluid system contains the endothermic end part i, the exact position is carried out. It is suitable for cases where it is impossible to grasp.

 [Number 11]

The temperature of all the endothermic end portions i and the inlet and outlet regions i is known, or the temperature measurement region i exists in the vicinity of all the endothermic end portions i and the inlet and outlet portion i (all endothermic regions i and the inlet and outlet regions i In the case where the temperature measurement site i is present within a predetermined predetermined distance range of and each of the endothermic end portions i and the inflow and outflow area i is upstream of the flow as viewed from the temperature measurement site i, the following equation (12a) ) And (12b), the downstream weighting function W ₁ (R _1ij , R _2ij ) may be set to the downstream force R _1ij , and the upstream weighting function W ₂ (R _1ij , R _2ij ) may be zero.

 [Number 12]

The temperature of all the endothermic areas i and the inlet and outlet parts i is known, or the temperature measurement part i exists in the vicinity of all the endothermic areas i and the inlet and outlet part i, and all the endothermic areas i and the inlet and outlet parts i Is at the downstream side of the flow as seen from the temperature measurement site i, the downstream weighting functions W ₁ (R _1ij , R _2ij ) are set to 0 and upstream weighting as shown in the following equations (13a) and (13b). The function W ₂ (R _1ij , R _2ij ) may be set to the upstream force R _2ij .

 [Num. 13]

If all of the positions in the fluid system of the endothermic end portion i or the inlet / outlet portion i which greatly contribute to the temperature distribution are specified, and some or all of the temperature of the endothermic end portion i and the inlet / outlet portion i are unknown, As shown in 14a) and (14b), the downstream weighting function W ₁ (R _1ij , R _2ij ) is the downstream force R _1ij , and the upstream weighting function W ₂ (R _1ij , R _2ij ) is represented by the formula (14b): ) Can be used. S _1j is the sum of the temperature measurement site i, the endothermic heat generating site i, and the downstream force R _1ij with respect to the inflow and outflow site i, whose temperature is known, and S _avej is the temperature measurement site i and the heat absorbing site where the temperature is known. It is the sum total of the average value (1/2) x (R _1ij + R _2ij ) of the downstream and upstream forces with respect to i and the inflow-exit site | part i.

 [Number 14]

The downstream weighting functions W ₁ (R _1ij , R _2ij ) and the upstream weighting functions W ₂ (R _1ij , R _2ij ) may be applied uniformly to all set temperature estimate points j, and the temperature estimate points j In each case, an appropriate downstream weighting function and an upstream weighting function that meet the conditions may be selectively used.

When the temperature of the designated part i is unknown (step S251, No), the downstream weight W _1ij and the upstream weight W _2ij of the designated part i in the temperature estimation point j are _set to "0" (step S257). Then, it is determined whether the downstream weight W _1ij and the upstream weight W _2ij are calculated for all the temperature estimation points j (step S259). If there is a temperature estimation point j in which the downstream weight W _1ij and the upstream weight W _2ij are not calculated (step S259: No), the process returns to step S249 and the above-described processing is repeated. And if the weight was computed about all the temperature estimation points j (step S259: Yes), it transfers to step S261.

In step S261, it determines whether the downstream weight _W1ij and the upstream weight _W2ij are computed about all the site | parts i. If there is a part i whose downstream weight W _1ij and the upstream weight W _2ij are not calculated (step S261: No), the process returns to step S247 and the above-described processing is repeated. And if downstream weight _W1ij and upstream weight _W2ij were computed about all the site | parts i (step S261: Yes), the downstream weight _W1ij and the upstream weight for each site | part i in the calculated temperature estimation point j W _2ij is stored in the storage device (step S263), and the weight calculation process is finished.

In order to estimate and visualize the temperature distribution of the entire _fluidic system, a temperature estimating point j is set throughout the entire _{fluidic system,} and the downstream weights W _1ij and the upstream weights W _2ij are calculated for all set temperature estimate points j. It needs to be placed. In this case, it is preferable to calculate the downstream weights W _1ij and the upstream weights W _2ij for each site i in advance for all the temperature estimation points j and save them as a database in the storage device.

Next, the downstream transfer time τ _{1 ij} , the upstream transfer time τ _{2 ij} and the downstream weight W _1ij , the upstream weight W _2ij and the site i for each temperature i for the temperature estimation point j calculated as described above. The procedure of the process (temperature estimation process) for estimating the temperature of the arbitrary temperature estimation point j is demonstrated using time series temperature data of the following. 12 is a flowchart showing a processing procedure of a temperature estimation process.

In the temperature estimation process, as shown in FIG. 12, the time t ₀ at which temperature estimation is performed is first determined (step S301). Next, the temperature estimation point j which estimates temperature is specified (step S303). The processing here can be realized by increasing the value of j sequentially in the range of 1 to N at each iteration. Subsequently, the downstream weight W _1ij , the upstream weight W _2ij , the downstream propagation time τ _{1 ij} and the upstream propagation time τ _{2 ij} for each site i with respect to the designated temperature estimation point j are read and obtained from the storage device (step). S305). For example, the downstream weight W _1ij , the upstream weight W _2ij , the downstream propagation time τ _{1 ij} and the upstream propagation time τ _{2 ij} are acquired for the temperature estimation point j specified from the database described above.

According to the above formula (3), the downstream known temperature T _i (t ₀ -τ) is obtained from the time-series temperature data T _i (t) of each acquired portion i, the downstream transfer time τ _{1 ij} and the upstream transfer time τ _{2 ij} . _1ij ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ) are determined (step S307), and the downstream weight W _1ij and the upstream weight W _2ij are calculated (step S309). Subsequently, a weighted average process using the calculated downstream side weight W _1ij and upstream side weight W _2ij is performed to calculate the estimated temperature Te _j (t ₀ ) of the temperature estimation point j at time t ₀ (step S311). ). Thereafter, the estimated temperature Te _j (t ₀ ) of the calculated temperature estimation point j is stored in the storage device (step S313).

Thereafter, it is determined whether or not the estimated temperature Te _j (t ₀ ) is calculated for all the temperature estimation points j (step S315). If the estimated temperature Te _j (t ₀ ) has an uncalculated temperature estimation point j (step S315: No), the process returns to step S303 and the above-described processing is repeated. On the other hand, if the estimated temperature Te _j (t ₀ ) of all the temperature estimation points j has been calculated (step S315: Yes), the temperature estimation process ends.

As described above, in the present embodiment, the fluid component of each region i (temperature measurement region i, endothermic region i, and inflow / exit region i) with respect to the entire fluid at an arbitrary temperature estimation point j in the fluid system. forces the ratio (R _1ij, R _2ij) acquired as a, and the obtained force (R _1ij, R _2ij), based on, the weight of each portion i of the temperature estimation point j W _ij or downstream weights W _1ij and upstream weight W _2ij , the propagation time (τ _{1 ij} , τ _{2 ij} ) were calculated. The time t at which the weight W _ij for each site i is weighted to the known temperature T _i in the corresponding site i and averaged (weighted average processing) or temperature estimation is performed using time series temperature data with respect to _0, the downstream weight downstream base temperature at the portion corresponding to the i _{W 1ij t i (t 0 -τ} 1 ij) , the upstream base temperature in the area i of the weight W to the upstream side _2ij The estimated temperature of the temperature estimation point j was calculated by weighting and averaging corresponding to T _i (t ₀ + τ _{2 ij} ). Therefore, in consideration of the advection diffusion of heat under the flow field of the fluid system, the temperature of the temperature estimation point j can be estimated with high precision. According to this, even if it is difficult to actually measure temperature, such as installing the temperature measuring device 2, it becomes possible to grasp | ascertain temperature with high precision. Therefore, high-precision temperature estimation which considers heat transportation by the flow of a fluid can be realized without restricting the arrangement of the temperature measuring device 2.

If the temperature estimate point j is set throughout the fluid system, the weight W _ij for each site i for each temperature estimate point j is calculated and stored in the storage device as a weight database, for example, this weight W _ij is obtained. with also read acquisition out, simply by substituting the above equation (2) to obtain the base temperature of the parts i T _i, the estimated temperature distribution of the oil system interpolates the base temperature of the parts i T _i the instantaneous It becomes possible. In the same way for a fluid system in which temperature fluctuations in time can occur, a temperature estimating point j is set throughout the entire fluid system, and the downstream weight W _{1ij for} each site i for each temperature estimating point j, the upstream weight W _2ij and the propagation times τ _{1 ij} and τ _{2 ij} are calculated and stored in a storage device as a database, for example, and the temperature and time observed for each part i in the temperature measuring device 2 are in order. As measured and _{stored as described above} , and the temperature measured value at an arbitrary time point t can be read as time series temperature data T _i (t), the downstream weight W _1ij , the upstream weight W _2ij, and the transfer time ( τ _{1 ij} , τ _{2 ij} ) are read and acquired, and the time t is instantaneously inserted only in equation (3) using the time t ₀ and time series temperature data T _i (t) to be estimated. To estimate the temperature distribution of the fluid system at _zero Become fluent. According to this, the temperature estimating apparatus 1 of this embodiment can fully utilize for the online monitoring of the industrial process in which real-time calculation is indispensable, and can be used for operation management and a control mechanism.

Since the upstream side force R _2ij is acquired and the weight W _ij or the downstream side weight W _1ij and the upstream side weight W _2ij are calculated using this upstream side force R _2ij and the upstream transfer time τ _{2 ij} , the actual temperature is measured. The temperature at the position that becomes the upstream side of the flow as viewed from the region i can also be estimated based on the known temperature T _i (temperature actual value) of the temperature measured region i and the time series temperature data T _i (t). In addition, when the temperature of the endothermic end part i or the inflow-out part i is known, the temperature of the temperature estimation point j can be estimated further using the known temperature of these endothermic part i or the inflow-out part i. According to this, it is not necessarily necessary to arrange the temperature measuring device 2 at the most upstream position of the flow of the fluid. Therefore, the temperature of any position in the fluid system can be estimated without restricting the arrangement of the temperature measuring device 2.

In the above-described embodiment, the fluid system has been described as including the temperature measurement site i, the endothermic end region i, and the inlet / outlet region i. However, when the fluid system does not include the endothermic end region i and / or the inlet / outlet region, these are The forces (R _1ij , R _2ij ) and the propagation time (τ _{1 ij} , τ _{2 ij} ) of the excluded region i are acquired, and the weight W _ij or the downstream weight W _1ij and the upstream weight W _2ij and the downstream known temperature T _{What is} necessary is just to calculate _i (t ₀ -τ _{1 ij} ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ). For example, if the fluid system does not include the endothermic end portion i, the forces (R _1ij , R _2ij ), and the transfer time (τ) of the temperature measured portion i and the inflow and outlet portion i of the fluid at the temperature estimation point j _{1 ij} , τ _{2 ij} ) and based on the acquired forces R _1ij , R _2ij and the propagation time τ _{1 ij} , τ _{2 ij} The weight W _ij or the downstream weight W _1ij and the upstream weight W _{2ij for each i} , the downstream known temperature T _i (t ₀ τ _{1 ij} ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ) are calculated. Do it. Similarly, when the fluid system does not include the inflow and outflow region i, the forces R _1ij and R _2ij and the propagation time τ _{1 ij} of the temperature measurement region i and the endothermic region i of the fluid at the temperature estimation point j , τ _{2 ij)} for acquiring the temperature measured area i and to the weight of the heat absorbing part for each i to the temperature estimation point j W _ij or downstream weights W _1ij and upstream weights W _2ij, and a downstream-side base temperature t _i (t _0- τ _{1 ij} ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ) are calculated, and when the fluid system does not include the endothermic end portion i and the inlet / outlet portion i, the temperature at the estimated point j The weights W _ij or downstream weights W for each temperature measurement site i for the temperature estimate point j are obtained by obtaining the forces R _1ij , R _2ij and the propagation time τ _{1 ij} , τ _{2 ij} of the temperature measurement site i of the fluid. _1ij and the upstream weight W _2ij and the downstream known temperature T _i (t ₀ τ _1ij ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ) may be calculated. According to this, by providing a temperature measuring device 2 such as a thermometer or a thermocouple at an arbitrary position where at least the installation space in the fluid system can be secured, it is based on the temperature measured value measured by the temperature measuring device 2. It is possible to accurately estimate the temperature at any position within the system.

In this embodiment, the weight W _ij or the downstream weight W _1ij and the upstream weight W _2ij are calculated based on the forces R _1ij and R _2ij , and the calculated weight W _ij or the downstream weight W _1ij and the upstream weight are calculated. Although W _2ij is _supposed to be stored in the storage device, the forces R _1ij and R _2ij may be _stored , and the weights W _ij or the downstream weights W _1ij and the upstream weights W _2ij may be calculated for each temperature estimation.

(Embodiment 1)

Next, as Embodiment 1, the room is applied, and the temperature estimation of the fluid system which flows through this room, and visualization of temperature distribution are demonstrated. 13 is a schematic view showing the inside of the room 3 to be applied in the first embodiment from above.

The room 3 shown in FIG. 13 has a substantially square shape when viewed from a plane where one side is 1 (m), for example. The room 3 is provided with passages 321 and 322 having a width of about 0.2 (m) in the diagonal of the left side wall 311 and the right side wall 312 toward Fig. 13, respectively. Windows 331, 332 are attached to the ends of 321, 322, respectively. The side wall 313 of the upper side toward the FIG. 13 of the room 3 is equipped with the heat generating source which is not shown in figure, and becomes a heat generating wall which radiates heat. In four places A-D shown by attaching "x" in FIG. 13 in the room 3, thermometers 34-1 to 34-4 as a temperature measuring device for estimating temperature are provided.

In this application object, the fluid system of the estimation object flows into the room through the passage 321 from the window 331 as shown by arrow A3 in FIG. 13, specifically, the air which flows in the room 3, Air flowing out of the window 332 to the outside via the passage 322. In this application object, the installation positions A-D of the thermometers 34-1 to 34-4 are the temperature actual measurement sites i, for example, a radius of 0.05 (m) centered on the installation positions A to D which are the temperature actual measurement sites i. Let the original areas E31 to E34 be the temperature actual measurement area i, respectively. The side wall 313 provided with the heat generating source is an endothermic heat generating part i (heating part), and the wall surface area E of this side wall 313 is made into the heat generating heat generating area i, for example. Areas F and G of the windows 331 and 332, that is, longitudinal sections of the passages 321 and 322 are respectively the inflow and outflow areas i (area F of the window 331 is the inflow site and area G of the window 332 is the outflow site, respectively. ), And this region F, G is an inflow-out area i, for example. As for the heat generating source which the side wall 313 has, the temperature is controlled to 50 (degreeC), for example, and 10 (degreeC) air flows in the window 331. However, at the time of performing temperature estimation, the temperature of the heat absorption part i and the inflow-out part i shall be unknown. Below, the temperature measurement site | part i corresponding to the installation positions A-D of the thermometers 34-1-34-4 is appropriately described as temperature measurement site | parts A-D, and the foot corresponded to the wall surface area E of the side wall 313. The endothermic part i is appropriately designated as the endothermic part E, and the inflow / exit area i corresponding to the area F of the window 331 is appropriately denoted as the inflow / outflow area F, and the inflow / outflow area corresponding to the area G of the window 332. i is appropriately denoted as the inflow and outflow area G.

14 is a block diagram showing the functional configuration of the temperature estimating apparatus 10 according to the first embodiment. As shown in FIG. 14, the temperature estimating apparatus 10 is provided with the input part 11, the display part 12, the memory | storage part 13, and the control part 14, and the thermometer installed in the room 3 ( The temperature measured values from 34-1 to 34-4 are input to the control unit 14.

The input unit 11 is for the user to perform various operations such as input of information required for temperature estimation, and outputs an input signal to the control unit 14. This input unit 11 is realized by a keyboard, a mouse, a touch panel, or the like. The display unit 12 is realized by a display device such as an LCD or an EL display, and under the control of the control unit 14, displays the result of the temperature estimation and the like.

The storage unit 13 is realized by various IC memories such as ROM and RAM, such as flash memory which can be updated and recorded, an information recording medium such as a hard disk, CD-ROM, or the like and a reading device or the like connected via built-in or data communication terminals. As the temperature estimating apparatus 10 is operated, a program for realizing various functions included in the temperature estimating apparatus 10, data used during the execution of the program, and the like are recorded. The storage unit 13 has a weight database in which the weight W _ij of the temperature estimation point j set in the room 3 is registered, or a room 3 of the temperature estimation point j corresponding to the estimated temperature Te _j of the temperature estimation point j. The temperature data and the like set in association with the position in the image are stored.

The control unit 14 is realized by hardware such as a CPU. The control unit 14 transmits an instruction or data to each unit constituting the temperature estimating apparatus 10 based on an input signal input from the input unit 11, a program or data recorded in the storage unit 13, or the like. The overall operation of the temperature estimating apparatus 10 is collectively controlled. The control unit 14 includes a temperature estimating unit 141 and a temperature distribution display processing unit 143.

The temperature estimating unit 141 calculates the weight W _ij for each site i with respect to the temperature estimation point j by performing the weight calculation process in accordance with the processing procedures shown in FIGS. The weight W _ij is stored in the storage unit 13 as a weight database. Specifically, the temperature estimating unit 141 uses, for example, a finite volume method as a numerical fluid simulation, and acquires forces R _1ij and R _2ij using a standard k-ε turbulence model as a turbulence model. Calculate the weight W _ij . In the first embodiment, the flow in the height direction of the fluid system is negligible, and the two-dimensional temperature distribution of the fluid system in the room 3 is estimated. Also, in the calculation of the weight W _ij , an approximation is made by the two-dimensional model. It shall be done.

In this case, the temperature estimating unit 141 calculates the two-dimensional flow field as the actual flow field as the processing of step S3 in FIG. 4. 15 is a schematic view showing an actual flow field of the room 3 calculated here. As shown in FIG. 15, in the calculation of the flow field, the flow velocity vector V3 indicating the flow of air in the entire area of the room 3, specifically, the window 331 (as indicated by the arrow A3 in FIG. 15) ( 13 is a flow velocity vector V3 indicating the direction of the flow at each position in the room 3 of the air 3 flowing into the room from the window 332 (see FIG. 13) and out to the outside from the window 332 (see FIG. 13). Obtained.

The temperature estimation part 141 sets the temperature estimation point j to the whole area | region in the room 3 at equal intervals as a process of step S7 of FIG. The temperature estimating unit 141 uses the actual flow field shown in FIG. 15 as the processing of steps S9 to S23, and as described above, at each temperature estimating point j set in the entire area in the room 3. The downstream force R _1ij of each part i of the fluid is obtained. 16A to 16G are diagrams each showing isoline diagrams of the downstream forces R _1ij of the respective portions A to G at the respective temperature estimation points j.

Thereafter, the temperature estimating unit 141 calculates the inverted flow field in which the direction of each flow velocity vector V3 in the reverse direction of the actual flow field shown in FIG. 15 is processed as the processing of Step S25 in FIG. 5. . The temperature estimating unit 141 is a process of steps S27 to S45 in FIG. 5, and uses the inversion flow field to determine the respective portions i of the fluid at each temperature estimation point j set in the whole area in the room 3. _{Obtain the} upstream forces R _2ij . 17A to _17G are _diagrams each showing the upstream forces R _2ij of the respective parts A to G at the temperature estimation point j in an _{isometric view} .

The temperature estimating unit 141 calculates the weight W _ij using the weighting functions W (R _1ij , R _2ij ) of the formula (5) and the weight of the formula (8) as the processing of steps S47 to S59 of FIG. 6. function W (R _1ij, R _2ij), a weight W, each subjected to the calculation of the _ij, the weighting function W (R _1ij, R _2ij) to data beyiseuhwa weights W _ij for each temperature point estimation j storage unit 13 each time by using To preserve. 18A to 18G are _diagrams each showing the weights W _ij for each of the parts A to G for each temperature estimation point j calculated using the weighting functions W (R _1ij and R _2ij ) of Equation (5), respectively. to be. 19A to 19G are _equivalent diagrams showing the weights W _ij for each of the parts A to G for each temperature estimation point j calculated using the weighting functions W (R _1ij and R _2ij ) of Equation (8), respectively. to be. In Embodiment 1, the temperature of the heat absorption part E and the inflow-out part F, G is unknown, and the weight W _ij of site | part E, F, G as shown to FIG. 18E-FIG. 18G or FIG. 19E-19G. Is calculated as "0".

The temperature estimating unit 141 performs the temperature estimation process in accordance with the processing procedure shown in FIG. Based on the known temperature T _i , the estimated temperature Te _j of each temperature estimation point j is calculated by using the weight W _ij for each temperature estimation point j. Then, the temperature estimating unit 141 stores the estimated temperature Te _j of each temperature estimation point j calculated in the storage unit 13 as temperature data.

The temperature distribution display processing unit 143 refers to the temperature data estimated by the temperature estimating unit 141 and stored in the storage unit 13, and is an estimated temperature Te _j of each temperature estimation point _j , that is, a known temperature T _i . The temperature distribution of the whole fluid system in the room 3 which interpolated the temperature actual value of temperature measurement site | part i (A-D), for example, is made into the isoline diagram, and is displayed on the display part 12 as a temperature distribution monitoring screen.

In the temperature estimating apparatus 10 of the structure demonstrated above, the case where the weight W _ij calculated by the weighting function W (R _1ij , R _2ij ) of Formula (5) mentioned above is used (Experimental example 1), and Formula (8) In the case of using the weight W _ij calculated by the weighting function W (R _1ij , R _2ij ) of) (Experimental Example 2), the estimated temperature Te _j of the pollutant at each temperature estimation point j was calculated. Moreover, as a comparative example, the weight was calculated using the reverse distance weighting method which is a conventional method, and the estimated temperature Te _j of each temperature estimation point _j was computed using the obtained weight. Calculation of the weight in a comparative example is performed in the same processing procedure as the weight calculation process of FIGS. 4-6, and uses the weight W _ij 'of the inverse distance interpolation formula shown by following formula (15) in step S53. . l _ij is the linear distance of the temperature measurement site | part i and the temperature estimation point j. u calculated the weight W _ij with u = 2 as an interpolation parameter, for example.

 [Number 15]

Specifically, the temperature estimating unit 141 calculates the estimated temperature Te _j of each temperature estimation point j using the three types of weights W _ij in Experimental Examples 1, 2 and Comparative Example, respectively, and the temperature distribution display processing unit. By (143) isolinearizing the estimated temperature Te _j of each temperature estimation point j in Experimental Examples 1 and 2 and Comparative Example, an estimation result for each of Experimental Examples 1 and 2 and Comparative Example was obtained. Table 1 shows the actual temperature T _i of the temperature measurement sites A to D used for the estimation, that is, the temperature measured values of the thermometers 34-1 to 34-4 installed at the corresponding installation positions A to D in the room 3. Indicates.

(Table 1)

For comparison, the temperature distribution of true in the room 3 was computed using the numerical fluid analysis as the temperature measured value shown in Table 1 as the known temperature T _i of each temperature measuring part A-D.

20 is an isometric view of the estimation result of Experimental Example 1 according to the first embodiment, that is, the temperature distribution of the fluid system in the room 3. FIG. 21: is a figure which shows the estimation result of Experimental example 2 in Embodiment 1, and FIG. 22 is a figure which shows the estimation result of the comparative example in Embodiment 1. FIG. FIG. 23: is a figure which shows the temperature distribution of the charm in the room 3. As shown in FIG. When comparing Experimental example 1 and 2 and a comparative example, as shown in FIG. 22, in a comparative example, the temperature distribution which an isoline spreads concentrically around the temperature measurement site | part A-D is estimated. As described above, in the comparative example, the flow of air in the room 3 was not reflected in the temperature estimation, resulting in an estimation result not corresponding to the true temperature distribution in the room 3 shown in FIG. 23. On the other hand, as shown to FIG. 20, 21, in Experimental example 1 and 2 of Embodiment 1, the temperature distribution which the isoline extended long along the direction of the flow of air in the room 3 is estimated, and FIG. The estimation result corresponding to the temperature distribution of the true in the room 3 shown to was obtained. Thus, according to Embodiment 1, the temperature estimation which reflected the flow of the air in the room 3 can be implement | achieved, and the temperature distribution in the room 3 could be reproduced with high precision.

(Embodiment 2)

Next, as Embodiment 2, the water tank is applied, and the temperature estimation of the fluid system which flows in this water tank, and visualization of temperature distribution are demonstrated. 24 is a schematic view showing the inside of the water tank 4 to be applied in the second embodiment from the side. FIG. 25 is a schematic view showing the inside of the water tank 4 in FIG. 14 from above.

As for the water tank 4 shown to FIG. 24 and FIG. 25, the depth direction (up-down direction of FIG. 25) is 1 (m), and the width direction (left-right direction of FIG. 24 and FIG. 25) is 1 (m), for example. It has a rectangular shape whose depth (width in the up-down direction of FIG. 24) is 0.5 (m), and the inside of the water tank 4 is a structure which is always filled with water. That is, the upper surface of the water tank 4 is provided with two pipes 41 and 42 communicating with the internal space of the water tank 4 at both left corners toward FIG. 25, and this pipe 41, Water is injected into the water tank 4 from 42. On the other hand, in the bottom face of the water tank 4, the one pipe 43 which communicates with the internal space of the water tank 4 in the right center toward FIG. 25 is provided, and from this pipe 43 The same amount of water flows out as the total amount of water introduced from (41, 42).

In the internal space of the water tank 4, the partition plate 44 which divides half of the width direction along the width direction of the water tank 4 is provided, and the position of the depth direction of the water tank 4 is the depth of the water tank 4; It is arrange | positioned as close to the pipe 41 side by 0.2 (m) with respect to the vertical cross section S4 which passes through the direction center. In six places P41-P46 which show "x" in FIG. 24 and 25 in the water tank 4, the thermometer 45-1-45-6 as a temperature measuring apparatus for estimating temperature is provided. The position of the depth direction of the thermometer 45-1-45-6 was made into the position which becomes the center depth exactly of the water tank 4. As shown in FIG.

In this application object, the fluid system of the estimation object is water flowing in the water tank 4, and the installation positions P41-P46 of the thermometers 45-1 to 45-6 become the temperature measurement site | part i. The lower ends of the pipes 41 and 42 are inflow sites, and the upper ends of the pipes 43 are outflow sites, and these are inflow and outflow sites i. However, the entire flow path of the pipes 41 to 43 may be the inflow-out part i, or the upper ends of the pipes 41 and 42 and the lower end of the pipe 43 may be the inflow-out part i. For example, 10 (degreeC) water is injected from the pipe 41, and 50 (degreeC) water is injected from the pipe 42. As shown in FIG. However, when performing temperature estimation, the temperature of the inflow-and-exit site | part i shall be unknown. The fluid system of this application object shall not contain the heat-absorbing heat | fever site | part i, since the water surface of the water filled in the water tank 4 and the heat transfer in the inner wall surface of the water tank 4 are small enough. Although not shown in FIG. 24, FIG. 25, also in Embodiment 2, corresponding temperature measurement area | region i and inflow-out area | region i are set with respect to these temperature measurement site | parts i and the inflow-out area | region i, respectively.

FIG. 26 is a block diagram showing the functional configuration of the temperature estimating apparatus 10a according to the second embodiment. In FIG. 26, the same code | symbol is attached | subjected to the structure similar to Embodiment 1. In FIG. As shown in FIG. 26, the temperature estimating apparatus 10a is provided with the input part 11, the display part 12, the memory | storage part 13, and the control part 14a, and the thermometer installed in the water tank 4 ( The temperature measured value from 45-1 to 45-6 is input to the control part 14a.

The memory | storage part 13 is the water tank 4 of the temperature estimation point j which corresponds to the weight database which registered the weight W _ij of the temperature estimation point j set in the water tank 4, and the estimated temperature Te _j of the temperature estimation point j. The temperature data etc. which were set corresponded with the position in the inside are stored.

The control unit 14a includes a temperature estimating unit 141, a temperature data extracting unit 142a, and a temperature distribution display processing unit 143a.

The temperature estimating unit 141 calculates the weight W _ij for each site i with respect to the temperature estimation point j by performing a weight calculation process in accordance with the processing procedures shown in FIGS. 4 to 6, and calculates the calculated weight estimation point j. The stored weight W _ij is stored in the storage unit 13 as a weight database. For example, the temperature estimating unit 141 uses the finite volume method as the numerical fluid simulation, _obtains the forces R _1ij and R _2ij using the standard k-ε turbulence model as the turbulence model, and weights W _ij. To calculate. At this time, the temperature estimating unit 141 sets the temperature estimation point j throughout the water tank 4 at intervals of 0.04 (m) along the depth direction, the width direction, and the depth direction as the processing of Step S7 in FIG. 4. do. In addition, the temperature estimation part 141 calculates the weight W _ij using the weighting function W (R _1ij , R _2ij ) of Formula (5) as a process of step S47-step S59 of FIG. 6, and Formula (8) the weighting function W (R _1ij, R _2ij) the weight is performed to calculate the W _ij, respectively, using the weighting function W (R _1ij, R _2ij) each beyiseuhwa weights W _ij for each temperature point estimation j data stored by unit ( 13) to preserve. In Embodiment 2, the temperature of inflow-out part i is unknown, and the weight W _ij of this inflow-out part i is computed as "0".

The temperature estimation part 141 performs temperature estimation processing according to the processing sequence shown in FIG. 7 similarly to Embodiment 1, and is based on the temperature measured value of the temperature measurement site | part i which is known temperature T _i , and is a temperature estimation point. The estimated temperature Te _j of _j is calculated and stored in the storage unit 13 as temperature data.

The temperature data extracting unit 142a refers to the temperature data estimated by the temperature estimating unit 141 and stored in the storage unit 13, and estimates the estimated temperature Te _j in an arbitrary cross section of the water tank 4. The estimated temperature Te _j ) of the temperature estimation point j in the cross section is extracted. The cross section for extracting the estimated temperature Te _j may be configured to be fixed in advance, or may be determined according to user operation. In the case of determining according to the user operation, the temperature data extracting unit 142a receives the operation of specifying the cross section by the user through the input unit 11, and extracts the estimated temperature Te _j in the cross section designated by the user.

The temperature distribution display processing part 143a visualizes the temperature distribution in this arbitrary cross section by making isometric linearization, for example, the estimated temperature Te _j in the arbitrary cross section which the temperature data extraction part 142a extracted, and shows the temperature. It is displayed on the display part 12 as a distribution monitoring screen.

In the temperature estimating apparatus 10a of the structure demonstrated above, the case where the weight W _ij calculated by the weighting function W (R _1ij , R _2ij ) of said Formula (5) is used (Experimental example 1), and Formula (8) In the case of using the weight W _ij calculated by the weighting function W (R _1ij , R _2ij ) of) (Experimental Example 2), the estimated temperatures Te _j of the respective temperature estimation points _j were calculated. Moreover, as a comparative example, the weight was calculated using the reverse distance weighting method which is a conventional method, and the estimated temperature Te _j of each temperature estimation point _j was computed using the obtained weight. Calculation of the weight in a comparative example is performed in the same process procedure as the weight calculation process of FIGS. 4-6, and uses the weight W _ij 'shown by Formula (15) in step S53.

Specifically, the temperature estimating unit 141 calculates the estimated temperature Te _j of each temperature estimating point j using the three types of weights W _ij in Experimental Examples 1, 2 and Comparative Example, respectively, and the temperature data extracting unit. 142a extracts the estimated temperature Te _j in the horizontal section passing through the center of the depth direction of the water tank 4 from the estimated temperature Te _j of each temperature estimation point j in Experimental example 1, 2, and a comparative example, The temperature distribution display processing unit 143a obtained the estimation results for each of Experimental Examples 1, 2 and Comparative Example by isolinearizing the estimated temperature Te _j in the above-described horizontal cross section. Base temperature of temperature measured area i by using the estimated T _i, that is, indicates the temperature measured values of the corresponding installation position provided in the P41-P46 thermometers (45-1 ~ 45-6) in the tub 4 in Table 2 below.

(Table 2)

27 is an isometric view of the temperature distribution in the horizontal cross section passing through the center of the depth direction of the water tank 4, that is, the estimation result of Experiment Example 1 in Embodiment 2. FIG. FIG. 28: is a figure which shows the estimation result of the experiment example 2 in Embodiment 2, and FIG. 29 is a figure which shows the estimation result of the comparative example in Embodiment 2. FIG.

In the comparative example using the reverse distance weighting method of the related art, the weight is calculated using only geometric information such as the distance between the temperature measurement site i and the temperature estimation point j as an index, but in this method, the structure in the fluid installation is not considered. In some cases, estimation results that differ greatly from the actual temperature distribution may be obtained. That is, for example, when applied to a fluid installation having a member for blocking the flow of fluid, such as the separator plate 44 installed in the interior space of the water tank 4 of the present application, continuous beyond the separator plate 44 The temperature distribution was sometimes estimated. However, in practice, since the flow of fluid is interrupted by the separator plate 44, the temperature may be discontinuous around the separator plate 44.

In fact, comparing the experimental examples 1 and 2 with the comparative example, as shown in FIG. 27 and FIG. 28, in the estimation results of the experimental examples 1 and 2 of the second embodiment, the temperature distribution was discontinuous around the separator plate 44. The temperature distribution reflecting the influence of the flow field rectified by the separator plate 44 is estimated. On the other hand, in the estimation result of the comparative example, as shown in FIG. 29, it becomes the continuous temperature distribution beyond the partition plate 44, and the influence of the flow field rectified by the partition plate 44 is not reflected in the temperature distribution. not. Thus, in Embodiment 2, unlike the case where the conventional reverse distance weighting method was used, temperature estimation reflecting the influence of the flow field can be realized, and it was confirmed that the temperature distribution in the water tank 4 can be accurately reproduced. Therefore, according to the second embodiment, even when a fluid system having a complicated three-dimensional flow field is used as the estimation target, high-precision temperature estimation can be realized.

In order to quantitatively verify the temperature estimation accuracy in the water tank 4, the experiment was performed further by installing the thermometer in the water tank 4 further. FIG. 30: is a figure which shows the installation position of the thermometer 45-7-45-9 which were further installed in the water tank 4. FIG. As shown in FIG. 30, the thermometers 45-7-45-9 are further installed in three places P47-P49 which show "x" in the water tank 4, and measure the temperature actual value of each installation position P47-P49. In addition, the temperature estimation by Experimental example 1, 2, and the comparative example was performed using each installation position P47-P49 as temperature measurement site | part i. The position of the depth direction of the thermometers 45-7-45-9 was made into the position which becomes the center depth exactly of the water tank 4 similarly to the thermometer 45-1-45-6.

Temperature measured value in each temperature measurement site | part i which is installation positions P47-P49 of the thermometers 45-7-45-9, the estimated temperature Te _j by the experiment example 1 of Embodiment 2, and the estimated temperature by experiment example 2 Te _j and the estimated temperature Te _j by a comparative example are shown in Table 3. As shown in Table 3, the estimated temperature Te _j obtained in Experimental example 1 and 2 of Embodiment 2 is set as the value close to the temperature actual value compared with the estimated temperature Te _j obtained by the comparative example, and according to Embodiment 2 It was confirmed that the temperature estimation accuracy was excellent.

(Table 3)

(Embodiment 3)

As the third embodiment, the temperature distribution estimation and visualization when the inflow water temperature changes with time in the same water tank as the second embodiment will be described. 10 degreeC water always flows in from the pipe 41 shown in FIG. 24 and FIG. Water of a certain flow rate flows in from the pipe 42, and the water temperature initially becomes 10 degreeC and 50 degreeC from the middle. The thermometers 45-1 to 45-6 are disposed at the same positions P41, P42, P43, P44, P45 and P46 as in the second embodiment.

In the third embodiment, the downstream force R _1ij , the upstream force R _2ij , the downstream transfer time τ _{1 ij} , the upstream transfer time τ _{2 ij} , the downstream weight W _1ij and the upstream weight W _2ij Calculated using a numerical fluid simulation in the same manner as Numerical fluid simulations use the finite volume method and the standard k-ε turbulence model as the turbulence model. In the flow field calculation, water flows in from the upper end of the pipe 41 at a flow rate of 0.765 L / s, water flows in from the upper end of the pipe 42 at a flow rate of 1.531 L / s, and at a constant pressure at the lower end of the pipe 43. In the upper surface of the water tank 4, sliding conditions, side walls, and bottom walls were calculated by providing boundary conditions as wall boundary conditions using the back-face side of the wall. The delivery time was calculated as an initial temperature of water of 27 ° C, a calorific value of 2,200 kW, and a threshold temperature of 28 ° C. The temperature estimation point j was arrange | positioned at 0.04 m space | interval, and was arrange | positioned throughout the tank 4.

Formulas (14a) and (14b) were used to calculate the weighting functions W _1ij and W _2ij . The time progress of the temperature measured in each position P41, P42, P43, P44, P45, P46 in the water tank 4 with the thermometer 45-1-45-6 is shown to FIG. 31A-31F.

As time t _{0 at} which temperature estimation is performed, six time points after 1 minute, 2 minutes, 3 minutes, 4 minutes, 5 minutes, and 6 minutes after the water temperature of the pipe 42 changed to 50 ° C. think. Regarding the temperature distribution in the horizontal cross section passing through the center of the water tank 4, the downstream transfer time τ _{1 ij} calculated above from the time series temperature data T _i (t) measured at the positions P41 to P46 and the upstream The downstream known temperature T _i (t ₀ -τ _{1 ij} ) and the upstream known temperature T _i (t ₀ + τ _{2 ij} ) were extracted using the side transfer time τ _{2 ij} . And water temperature Te _j was estimated using Formula (14a) and Formula (14b).

Horizontal cross section through the center of the water tank 4 after 1 minute, 2 minutes, 3 minutes, 4 minutes, 5 minutes, and 6 minutes after the water temperature which flowed in from the pipe 42 changed to 50 degreeC Equivalence linearization was carried out about the water temperature in. Equivalence diagrams are shown in Figs. 32A to 32G. When the temperature of the water flowing from the pipe 42 changes from 10 ° C. to 50 ° C., the temperature gradually increases from the position close to the pipe 42, and even if there is a time transition of the temperature distribution, the temperature distribution is maintained. It could be estimated.

(Fourth Embodiment)

Next, as Embodiment 4, the temperature estimation and the visualization of the temperature distribution of the fluid system which flows in this hot dip galvanizing port as an application object are demonstrated. FIG. 33 is a schematic view showing the inside of the hot-dip galvanizing port 5 to be applied in Embodiment 4 from the side. In the hot dip galvanizing line, which is one of steel processes for manufacturing galvanized steel sheet used for automobiles and building materials, the steel sheet 51 is immersed in hot dip zinc in the hot dip galvanizing port 5 as illustrated in FIG. 33. The coating amount is adjusted using a coating amount control device (not shown), and predetermined post-treatment such as cooling is performed to obtain a coated steel sheet. The operating conditions are, for example, a line speed of 120 (mpm) and a sheet width of the steel sheet to 1,500 (mm).

The capacity of the molten zinc in the hot dip galvanizing port 5 of FIG. 33 to be applied in Embodiment 4 is 250 (t), for example, and the inside of the hot dip galvanizing port 5 is filled with molten zinc. This hot dip galvanizing port 5 is equipped with the induction heating apparatus 52 provided in each of the opposing inner wall surfaces parallel to the ground of FIG. The hot dip galvanizing port 5 is equipped with the ingot input part (not shown) for injecting a zinc ingot 53 into an internal space. The induction heating apparatus 52 melt | dissolves the zinc ingot 53 thrown into this ingot input part, and makes it molten zinc, and is for maintaining the temperature of this melted molten zinc at predetermined temperature.

A sink roll 54 is provided in the inner space of the hot dip galvanizing port 5, and the sink roll 54 is a plate of the steel plate 51 immersed in the hot dip zinc and conveyed in the hot dip zinc. The direction is to be redirected. Zinc consumed by adhesion to the steel plate 51 is supplied by the addition of the zinc ingot 53 to the ingot inlet (not shown). In eight places P51 to P58 denoted by "x" in FIG. 31 in the hot-dip galvanizing port 5, thermocouples 55-1 to 55-8 as a temperature measuring device for estimating the temperature are provided. Each thermocouple 55-1-55-8 is formed from one inner wall surface parallel to the ground surface of FIG. 33 of the hot dip galvanizing port 5, for example, from the inner wall surface of the front surface of the surface viewed from the sink roll 54. In the surface which becomes the distance of 300 (mm), each is provided in the positional relationship shown in FIG.

In this application object, the fluid system of the estimation object is molten zinc filled in the hot dip galvanizing port 5, and the installation positions P51-P58 of the thermocouples 55-1 to 55-8 become the temperature measurement site i. In Embodiment 4, the heating position of the induction heating apparatus 52 which is a site | part in the fluid system which heat generate | occur | produces, and two places of an ingot input part become a heat absorption part i. The temperature of the heat absorption part i corresponding to the heating position of the induction heating apparatus 52 is known, and it is specifically, 487.72 (degreeC) (refer Table 4). In addition, about the temperature of the heat absorption part i corresponding to an ingot input part, it is unknown. In the present application object, since the inflow and outflow of the fluid system of the estimation object, that is, the inflow of molten zinc into the molten zinc plating port 5 and the outflow of molten zinc out of the molten zinc plating port 5 do not exist, The fluid system of the form 4 shall not contain inflow-out part i. In Embodiment 4, the corresponding temperature measurement area | region i and heat absorption area | region i (not shown) are set with respect to these temperature measurement site | parts i and the heat absorption part i, respectively.

34 is a block diagram showing the functional configuration of the temperature estimating apparatus 10b according to the fourth embodiment. In FIG. 34, the same code | symbol is attached | subjected to the structure similar to Embodiment 1. In FIG. As shown in FIG. 34, the temperature estimating apparatus 10b is provided with the input part 11, the display part 12, the memory | storage part 13, and the control part 14b, and it is in the hot-dip galvanizing port 5 The temperature measured values from the installed thermocouples 55-1 to 55-8 are input to the control unit 14b.

The memory | storage part 13 melt | dissolves the temperature estimate point j which corresponds to the weight database which registered the weight W _ij of the temperature estimate point j set in the hot-dip galvanizing port 5, and the estimated temperature Te _j of the temperature estimate point j. The temperature data etc. which were set in correspondence with the position in the galvanizing port 5 are preserve | saved.

The control unit 14b includes a temperature estimating unit 141, a temperature data extracting unit 142b, and a temperature distribution display processing unit 143b.

The temperature estimating unit 141 calculates the weight W _ij for each site i with respect to the temperature estimation point j by performing a weight calculation process in accordance with the processing procedures shown in FIGS. 4 and 5, and calculates the calculated weight estimation point j. The stored weight W _ij is stored in the storage unit 13 as a weight database. For example, the temperature estimating unit 141 uses the finite volume method as the numerical fluid simulation, _obtains the forces R _1ij and R _2ij using the standard k-ε turbulence model as the turbulence model, and weights W _ij. To calculate. At this time, the temperature estimation part 141 sets the temperature estimation point j to the whole area | region in the hot-dip galvanizing port 5 at equal intervals as a process of step S7 of FIG. In addition, the temperature estimation part 141 calculates the weight W _ij using the weighting functions W (R _1ij , R _2ij ) of Formula (5), for example as a process of step S47-step S59 of FIG. The weight W _ij of each temperature estimation point j is converted into a database and stored in the storage unit 13. In Embodiment 4, the temperature of the heat absorption part i corresponded to an ingot input part is unknown, and the weight W _ij of the heat absorption part i corresponded to this ingot input part is calculated as "0".

The temperature estimation part 141 performs the temperature estimation process according to the processing sequence shown in FIG. 7 similarly to Embodiment 1, and the temperature actual measurement value and the induction heating apparatus 52 of the temperature measurement site | part i which are known temperature T _{i are carried out} . The estimated temperature Te _j of each temperature estimation point j is calculated on the basis of the temperature of the heat absorption part i corresponding to the heating position of and stored in the storage unit 13 as temperature data.

The temperature data extracting unit 142b refers to the temperature data estimated by the temperature estimating unit 141 and stored in the storage unit 13 to estimate the estimated temperature Te _j in any cross section of the hot dip galvanizing port 5. Extract The temperature distribution display processing part 143b visualizes the temperature distribution in this arbitrary cross section by making isometric linearization, for example, the estimated temperature Te _j in the arbitrary cross section extracted by the temperature data extraction part 142b, and temperature It is displayed on the display part 12 as a distribution monitoring screen.

In the temperature estimating apparatus 10b of the structure demonstrated above, the temperature estimation of each temperature estimation point j was performed. Specifically, the temperature estimation part 141 is, calculates the estimated temperature Te _j of temperature estimation point j, and the temperature data extraction section (142b), hot-dip galvanized port 5, for the inner side wall g The estimated temperature Te _j in the 300 mm vertical cross section was extracted, and the temperature distribution display processing part 143b obtained the estimation result by isolinearizing the estimated temperature Te _j in the vertical cross section mentioned above. Base temperature of temperature measured area i by using the estimated T _i, i.e., the actual measurement value of the thermocouple temperature (55-1 ~ 55-8) installed in the corresponding installation position P51 - P58 in the hot-dip galvanized port 5, guided It shows in Table 4 with the known temperature T _i of the heating position of the heating apparatus 52.

 (Table 4)

FIG. 35 is a diagram showing an isometric view of the estimation result in the fourth embodiment, that is, the temperature distribution in the vertical cross section of 300 mm from the inner wall surface of the hot dip galvanizing port 5. According to the fourth embodiment, temperature estimation reflecting the influence of the flow of molten zinc in the hot dip galvanizing port 5 can be realized, and the temperature distribution can be estimated with high accuracy in consideration of the flow effect of the molten zinc. In addition, since the weight W _ij for each temperature estimation point j is stored in the storage unit 13 as the weight database, the weight is given based on the known temperature T _i during the temperature estimation and the visualization of the temperature distribution. It was good just to perform an average process, and the calculation time was able to be within 1 second. Therefore, it is also possible to visualize the temperature distribution online (real time). In addition, when the temperature of the molten zinc in the hot dip galvanizing port 5 is not within a predetermined range, it is known that a surface defect occurs in the hot dip galvanized steel sheet. Therefore, the surface of the induction heating apparatus 52 is estimated by estimating the temperature of the molten zinc in the hot dip galvanizing port by the above-described processing, and controlling the temperature of the hot dip zinc in the hot dip galvanizing port to be within a predetermined range based on the estimation result. It is possible to produce a hot dip galvanized steel sheet without defects.

Specifically, as shown in FIG. 34, the control unit 14b determines whether or not the temperature of the molten zinc in the predetermined region in the hot dip galvanizing port 5 is within a predetermined threshold. ) And a temperature control unit 145b which controls the temperature of the molten zinc by operating the output of the induction heating apparatus 52 of the hot dip galvanizing port 5. In addition, the "predetermined area | region in the molten zinc plating port 5" means the place where the surface of the steel plate 51 and molten zinc contact in contact with a surface defect, the sink roll 54, and molten zinc, for example. The part which contacts, the area | region enclosed by the upper part of the sink roll 54, and the steel plate 51, etc. are said. The threshold value of the molten zinc temperature is previously input to the determination unit 144b or input by the operator through the input unit 11, and the determination unit 144b is located in a predetermined region extracted by the temperature data extraction unit 142b. It is determined whether or not the molten zinc temperature in the threshold is within. When the determination part 144b determines that the molten zinc temperature in a predetermined | prescribed area | region is out of a threshold range, the temperature control part 145b will inductively heat so that the temperature of the molten zinc in a predetermined | prescribed area may be in a threshold range. Manipulate the output of the device 52. According to the fourth embodiment, the temperature control unit 145b can control the induction heating apparatus 52 to control the molten zinc temperature in the predetermined region. Thereby, the surface defect of the steel plate 51 can be prevented.

(Embodiment 5)

Next, as Embodiment 5, the tundish for continuous casting is applied and the temperature estimation of the fluid system which flows in this tundish, and visualization of temperature distribution are demonstrated. 36 is a perspective view schematically showing a configuration of a tundish 6 to be applied in the fifth embodiment. FIG. 37: is a figure which shows the installation position of the thermocouple 64-1-64-5 provided in the tundish 6 of Embodiment 5, in the right half toward the FIG. 37 of the long side of the tundish 6, FIG. This figure schematically shows the interior thereof.

The tundish 6 shown in FIG. 36 has a rectangular shape having a depth direction of 1 (m), a width direction of 8 (m), and a height of 1 (m), and accommodates molten steel therein. In FIG. 36, the liquid surface S7 of the molten steel accommodated in the tundish 6 is shown with the broken line. The tundish 6 is provided at a nozzle 61 for injecting molten steel from a ladle and at two locations at the bottom thereof, and outlet holes 62 and 62 for introducing molten steel into a mold. ) And two plasma heating devices (63, 63) for heating the molten steel to control the temperature. The tundish 6 has a two-strand specification in which a nozzle 61 for injecting molten steel from the ladle is provided at the upper center, and outlet holes 62 and 62 to the mold are provided at both ends in the width direction.

In Embodiment 5, in the right half of the long side direction of the tundish 6 shown in FIG. 37, temperature is estimated to five places P61-P65 which show "x" in FIG. 37 in the perpendicular surface which passes through the short side direction center. The thermocouples 64-1 to 64-5 as a temperature measuring device for this purpose are provided. The thermocouples 64-1 to 64-5 are provided only at the right side of the long side direction of the tundish 6 because the tundish 6 has a symmetrical structure, but the long side direction is provided. Similarly, thermocouples 64-1 to 64-5 may be provided on the left side to be used for temperature estimation.

In this application object, the fluid system of the estimation object flows into the inside of the tundish 6 from the molten steel accommodated in the inside of the tundish 6, specifically, from the lower end of the nozzle 61, and the outflow hole 62 , 62 is molten steel flowing out of the tundish 6 to the outside (mould). In this application object, the mounting positions P61 to P65 of the thermocouples 64-1 to 64-5 become the temperature measurement site i. The heating position of the plasma heating apparatuses 63 and 63 becomes the heat absorption part i, and the lower end of the nozzle 61 and the outlet holes 62 and 62 are the inflow and outflow part i (the lower end of the nozzle 61 is the inflow part and the outflow part). The holes 62 and 62 serve as outflow sites. Since the liquid surface S7 of molten steel receives strong cooling from the outside, it is an endothermic part i. The heat absorbing portion i corresponding to the heating position of the plasma heating apparatus 63, 63 and the inlet / outlet portion i corresponding to the inlet position of the lower end of the nozzle 61 have a known temperature and correspond to the outlet holes 62, 62. The temperature of the heat absorption part i corresponding to the inflow-out part i and the liquid level S7 of molten steel is unknown.

38 is a block diagram showing the functional configuration of the temperature estimating apparatus 10c according to the fifth embodiment. In FIG. 38, the same code | symbol is attached | subjected to the structure similar to Embodiment 1. In FIG. As shown in FIG. 38, the temperature estimating apparatus 10c is provided with the input part 11, the display part 12, the memory | storage part 13, and the control part 14c, and the thermocouple provided in the tundish 6 The temperature measured value from (64-1 to 64-5) is input to the control part 14c.

The storage unit 13 is a weight database in which the weight W _ij of the temperature estimation point j set in the tundish 6 is registered, or the tundish of the temperature estimation point j corresponding to the estimated temperature Te _j of the temperature estimation point j ( 6) The temperature data etc. which were set corresponded with the position in inside are stored.

The control unit 14c includes a temperature estimating unit 141, a temperature data extracting unit 142c, a temperature distribution display processing unit 143c, a determining unit 144c, and a temperature control unit 145c.

The temperature estimating unit 141 calculates the weight W _ij for each site i with respect to the temperature estimation point j by performing the weight calculation process in accordance with the processing procedures shown in FIGS. The weight W _ij is stored in the storage unit 13 as a weight database. For example, the temperature estimating unit 141 uses the finite volume method as the numerical fluid simulation, _obtains the forces R _1ij and R _2ij using the standard k-ε turbulence model as the turbulence model, and weights W _ij. To calculate. At this time, the temperature estimating unit 141 sets the temperature estimation point j throughout the tundish 6 at equal intervals as the processing of step S7 in FIG. 4. In Embodiment 5, the temperature of the heat absorption part i corresponding to the heating position of the plasma heating apparatus 63, 63 and the inflow-out part i corresponded to the outflow holes 62, 62 are unknown, and The weight W _ij is calculated as "0".

The temperature estimation part 141 performs temperature estimation processing according to the processing sequence shown in FIG. 7 similarly to Embodiment 1, and it estimates each temperature based on the temperature measured value of the temperature measurement site | part i which is known temperature T _i . The estimated temperature Te _j of _j is calculated and stored in the storage unit 13 as temperature data.

The temperature data extracting unit 142c extracts the estimated temperature Te _j in any cross section of the tundish 6 with reference to the temperature data estimated by the temperature estimating unit 141 and stored in the storage unit 13. do. Since the molten steel in the tundish 6 is cooled at its contact surface with its bath surface or the inner wall surface of the tundish 6, the molten steel injected from the ladle flows downward to approach the outlet holes 62 and 62. As a result, the temperature decreases. As for the tundish 6, the inner wall surface is covered with the refractory, and this refractory is always in contact with hot molten steel. When the temperature of the molten steel which contacts this refractory changes rapidly, a large thermal stress will generate | occur | produce in a refractory and the problem of refractory damage will arise. Therefore, it is preferable to control the temperature of the molten steel so that the temperature of the molten steel in contact with the inner wall surface is within a predetermined threshold. So, in Embodiment 5, the temperature data extraction part 142c is a vertical cross section of the inner wall surface vicinity of the longitudinal direction of the tundish 6, for example, the inner wall surface vicinity of the front side of the paper surface of FIG. The estimated temperature Te _j is extracted.

The temperature distribution display processing unit 143c estimates the estimated temperature Te _j in an arbitrary cross section (for example, a vertical cross section near the inner wall surface in the long side direction of the tundish 6) extracted by the temperature data extraction unit 142c. For example, by isolinearizing, the temperature distribution in this arbitrary cross section is visualized and displayed on the display part 12 as a temperature distribution monitoring screen.

The determination unit 144c determines whether or not the temperature of the molten steel in the vicinity of the inner wall surface of the tundish 6, that is, the contact portion with the refractory covering the inner wall surface is within a predetermined temperature range. For example, the determination unit 144c determines whether or not the maximum or minimum temperature of the estimated temperature Te _{j in} the vertical section extracted by the temperature data extraction unit 142c is within a predetermined temperature range. The predetermined temperature range may be set to be fixed in advance, or may be determined according to user operation. In the case of determining according to the user operation, the input operation of the temperature range by the user is accepted through the input part 11, and the determination part 144c makes the above-mentioned determination according to the temperature range input by the user.

The temperature control part 145c controls the heating temperature by the plasma heating apparatus 63, 63 according to the determination result which the determination part 144c performed. Specifically, when the temperature control unit 145c determines outside the temperature range in the determination unit 144c, the temperature control unit 145c controls the temperature of the plasma heating apparatus 63, 63 so that the maximum or minimum temperature determined outside the temperature range is within the temperature range. Control the output.

In the temperature estimating apparatus 10c of the structure demonstrated above, the temperature estimation of each temperature estimation point j was performed. Specifically, a temperature estimation unit 141 calculates the estimated temperature Te _j of temperature estimation point j, and the temperature data extraction section (142c) have, for example the inner wall surface near the (wall of the tundish (6) The estimated temperature Te _j in the vertical cross section of 50 mm) was extracted, and the temperature distribution display processing part 143c obtained the estimation result by isolinearizing the estimated temperature Te _j in the above-mentioned vertical cross section. Base temperature of temperature measured area i by using the estimated T _i, i.e., the tundish 6, the temperature actually measured value of the corresponding installation position provided in the P61-P65 thermocouple (64-1 ~ 64-5), a plasma heating apparatus in a It shows in Table 5 with the known temperature T _i of the heating position of 63, and the known inflow temperature T _i of the nozzle 61.

(Table 5)

FIG. 39 is a diagram showing isothermal linearization of the temperature estimation result in the fifth embodiment, that is, the temperature distribution in the vertical section near the inner wall surface (50 mm from the wall surface) of the tundish 6. As shown in FIG. 39, according to the fifth embodiment, temperature estimation reflecting the influence of the molten steel flow in the tundish 6 can be realized, and the temperature distribution can be estimated with high accuracy in consideration of the flow effect of the molten steel. Furthermore, for it as a tundish (6) to extract the estimated temperature Te _j in the vertical cross-section of the inner wall near the long side direction, and the estimated temperature Te _j in this section, such to present the example contour line Drawing Therefore, the user can easily grasp | ascertain the temperature of the molten steel in the contact part with the refractory which covers an inner wall surface. Moreover, when the temperature of the molten steel in this contact part is outside the predetermined temperature range, since the output of the plasma heating apparatus 63, 63 can be controlled and temperature control of molten steel can be performed, it is inside of the tundish 6 Damage to the refractory covering the wall can be prevented.

In the above-mentioned embodiment, although the room, the water bath, the hot-dip galvanizing port, and the tundish for continuous casting were illustrated as the object of application of this invention, it is not limited to these, This invention applies widely if a fluid is related. can do. For example, in a steel process, it is applicable to the temperature estimation in a molten metal preservation hold | maintenance, a continuous casting mold, ladle, etc. In addition, the present invention is not limited to the steel field, but may be similarly applied to chemical processes and water treatment facilities. In addition, the present invention can be adapted to a fluid system of a wide flow state not only from a simple one-dimensional flow fluid system but also to a fluid three-dimensional flow system.

As described above, the method for estimating the temperature of the fluid system, the method for estimating the temperature distribution of the fluid system, the method for monitoring the temperature distribution of the fluid system, and the temperature estimating device provide a flow of fluid without restricting the arrangement of the temperature measuring device. It is suitable for realizing high-precision temperature estimation considering heat transfer by Moreover, according to the molten zinc temperature control method and the hot dip galvanized steel sheet in the hot dip galvanized pot of this invention, a hot dip galvanized steel plate without surface defects can be provided. Moreover, according to the molten steel temperature control method in the tundish of this invention, refractory damage of a tundish can be suppressed.

1, 10, 10a, 10b, 10c: temperature estimation device
2: temperature measuring device
11: Input unit
12: display unit
13: memory
14, 14a, 14b, 14c: control unit
141: temperature estimation unit
142a, 142b, 142c: temperature data extraction section
143, 143a, 143b, 143c: temperature distribution display processing unit
144b and 144c: Determination unit
145b, 145c: temperature control unit
3: room
4: tank
34-1 ～ 34-4, 45-1 ～ 45-9: Thermometer
5: hot dip galvanized port
52: induction heating apparatus
6: tundish
63: plasma heating device
55-1 ～ 55-8, 64-1 ～ 64-5: thermocouple

Claims (14)

As a temperature estimation method of a fluid system which estimates the temperature in the arbitrary temperature estimation point of the fluid system which has two or more temperature known areas,
A temperature generated through the temperature known area or generated within the temperature known area by using the position information of the temperature known area and information on a flow field of the fluid system indicating the flow of the fluid throughout the fluid system. A force acquisition step of acquiring, as a force in the temperature known area at the temperature estimation point, the proportion of the fluid that reaches the temperature estimated point in the fluid, up to the temperature estimated point, without passing through another temperature known area; ,
And a temperature estimating step of estimating the temperature at the temperature estimating point using information on the temperature of each temperature known region and the forces at the temperature estimating point. The method of claim 1,
The force acquiring step is based on an advective diffusion phenomenon caused by the flow field, without passing through another temperature known region of the fluid passing through the temperature known region or generated in the temperature known region. And a downstream force acquiring step of acquiring, as a downstream force of the temperature known region at the temperature estimate point, the proportion of the fluid that has reached the temperature estimate point in the total fluid of the temperature estimate point,
The temperature estimating step includes a step of estimating the temperature of the temperature estimating point using information on the downstream force of each temperature base region in the temperature estimating point. . The method according to claim 1 or 2,
The force acquiring step is generated in the temperature known region or through the temperature known region according to the advection diffusion phenomenon caused by the reverse flow field of the fluid system indicating the flow of the fluid and the reverse direction. Upstream side which acquires the ratio which occupies in the whole fluid of a temperature estimation point of the fluid which reached the temperature estimation point without passing through another temperature base region among the used fluid as an upstream force of the temperature base region in a temperature estimation point. Including power acquisition process,
The temperature estimating step includes a step of estimating the temperature of the temperature estimating point using information on the upstream force of each temperature known region in the temperature estimating point. . 4. The method according to any one of claims 1 to 3,
A time series temperature data acquisition step of acquiring time series temperature data including a temperature of a temperature known region where a temperature is known in a time series;
A transfer time acquisition step of acquiring a transfer time required when the fluid moves between the temperature known region and the temperature estimation point,
The temperature estimating step extracts the temperature of the temperature known region at the extraction time point from the time series temperature data using the time point in the past or the future as the extraction time point with respect to the time point at which temperature estimation is performed. And estimating the temperature of the temperature estimating point using the extracted temperature. 5. The method according to any one of claims 1 to 4,
A weight calculation step of calculating a weight of each temperature known area by using information on the forces in each temperature known area,
The temperature estimating step includes a step of estimating a weighted average process using the weight of each temperature known area to estimate the temperature of the temperature estimating point. The method of claim 5,
The temperature estimating step is a case in which the fluid system includes one or more endothermic sites where exothermic or endothermic occurs and / or one or more inflow / exit sites through which fluid flows into or out of the system. And estimating the temperature of the temperature estimation point by setting the value of the weight value for the portion whose temperature is unknown to 0 when the temperature is unknown. . 7. The method according to any one of claims 1 to 6,
And the fluid system is molten zinc in the hot dip galvanizing port. 7. The method according to any one of claims 1 to 6,
And the fluid system is molten steel in a tundish. A temperature distribution estimation method of a fluid system having a temperature distribution,
Using the temperature estimation method of the fluid system in any one of Claims 1-8, the temperature of the temperature estimation point set across the said fluid system is estimated,
And estimating a temperature estimated for each of the temperature estimation points as a temperature distribution of the fluid system. A temperature distribution monitoring method of a fluid system having a temperature distribution,
Visualizing and displaying the temperature distribution in any cross section of the fluid system based on the temperature distribution of the fluid system estimated using the method for estimating the temperature distribution of the fluid system according to claim 9. A method for monitoring the temperature distribution of a fluid system. A temperature estimating apparatus for estimating the temperature at an arbitrary temperature estimation point of a fluid system having two or more temperature known areas,
Among the fluids having passed through the temperature known area or generated in the temperature known area by using the position information of the temperature known area and information about the flow field of the fluid system indicating the flow of the fluid throughout the fluid system, Force acquiring means for acquiring, as a force in the temperature known area at the temperature estimated point, the proportion of the fluid reaching the temperature estimated point in the total fluid of the temperature estimated point without passing through the temperature known area;
And a temperature estimating means for estimating the temperature at the temperature estimating point using information on the temperature of each temperature known area and the force at the temperature estimating point. A hot dip zinc temperature control method in a hot dip galvanizing port,
A temperature extraction step of extracting a temperature of molten zinc in a predetermined region in the hot dip galvanizing port from the hot dip galvanizing port data estimated by the temperature estimating method of the fluid system according to claim 7; ,
A determination step of determining whether the extracted temperature is within a predetermined threshold range,
And in the determination step, a control step of operating an output of the heating means of the hot-dip galvanizing port so that the extracted temperature is within a threshold range when the extracted temperature is determined to be outside the threshold range. Method for controlling molten zinc temperature in galvanized port. It is manufactured using the molten zinc temperature control method in the hot dip galvanizing port of Claim 12, The hot dip galvanized steel plate characterized by the above-mentioned. As a molten steel temperature control method in a tundish,
A temperature extraction step of extracting a temperature of molten steel in a predetermined region in the tundish from the molten steel temperature data in the tundish estimated by the temperature estimation method of the fluid system according to claim 8;
A determination step of determining whether the extracted temperature is within a predetermined threshold range,
And in the determination step, a control step of operating an output of the heating means of the tundish so that the extracted temperature is within a threshold range when the extracted temperature is determined to be outside the threshold range. Molten steel temperature control method.

Images