JP2012233869A - Method for estimating temperature in fluid system, method for estimating temperature distribution in fluid system, method for monitoring temperature distribution in fluid system, temperature estimation device, method for controlling temperature of molten zinc in hot-dip galvanization pot, molten zinc plated steel plate and method for controlling temperature of liquid steel in tundish - Google Patents
Method for estimating temperature in fluid system, method for estimating temperature distribution in fluid system, method for monitoring temperature distribution in fluid system, temperature estimation device, method for controlling temperature of molten zinc in hot-dip galvanization pot, molten zinc plated steel plate and method for controlling temperature of liquid steel in tundish Download PDFInfo
- Publication number
- JP2012233869A JP2012233869A JP2011210259A JP2011210259A JP2012233869A JP 2012233869 A JP2012233869 A JP 2012233869A JP 2011210259 A JP2011210259 A JP 2011210259A JP 2011210259 A JP2011210259 A JP 2011210259A JP 2012233869 A JP2012233869 A JP 2012233869A
- Authority
- JP
- Japan
- Prior art keywords
- temperature
- fluid system
- estimation
- fluid
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
- G01K7/42—Circuits effecting compensation of thermal inertia; Circuits for predicting the stationary value of a temperature
- G01K7/427—Temperature calculation based on spatial modeling, e.g. spatial inter- or extrapolation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/02—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow
-
- C—CHEMISTRY; METALLURGY
- C23—COATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; CHEMICAL SURFACE TREATMENT; DIFFUSION TREATMENT OF METALLIC MATERIAL; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL; INHIBITING CORROSION OF METALLIC MATERIAL OR INCRUSTATION IN GENERAL
- C23C—COATING METALLIC MATERIAL; COATING MATERIAL WITH METALLIC MATERIAL; SURFACE TREATMENT OF METALLIC MATERIAL BY DIFFUSION INTO THE SURFACE, BY CHEMICAL CONVERSION OR SUBSTITUTION; COATING BY VACUUM EVAPORATION, BY SPUTTERING, BY ION IMPLANTATION OR BY CHEMICAL VAPOUR DEPOSITION, IN GENERAL
- C23C2/00—Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor
- C23C2/04—Hot-dipping or immersion processes for applying the coating material in the molten state without affecting the shape; Apparatus therefor characterised by the coating material
- C23C2/06—Zinc or cadmium or alloys based thereon
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24H—FLUID HEATERS, e.g. WATER OR AIR HEATERS, HAVING HEAT-GENERATING MEANS, e.g. HEAT PUMPS, IN GENERAL
- F24H1/00—Water heaters, e.g. boilers, continuous-flow heaters or water-storage heaters
- F24H1/18—Water-storage heaters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Combustion & Propulsion (AREA)
- Thermal Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
- Heat-Pump Type And Storage Water Heaters (AREA)
Abstract
Description
本発明は、流体系の温度推定方法、流体系の温度分布推定方法、流体系の温度分布モニタリング方法、温度推定装置、溶融亜鉛めっきポット内の溶融亜鉛温度制御方法、溶融亜鉛めっき鋼板、およびタンディッシュ内の溶鋼温度制御方法に関する。 The present invention relates to a fluid system temperature estimation method, a fluid system temperature distribution estimation method, a fluid system temperature distribution monitoring method, a temperature estimation device, a hot dip galvanizing temperature control method in a hot dip galvanizing pot, a hot dip galvanized steel sheet, and a tongue. The present invention relates to a method for controlling the temperature of molten steel in a dish.
産業プロセスで用いられる流体設備や建造物内部等の壁で区切られた空間を流れる流体系の温度分布を直接知るためには、流体系内に温度分布の特徴を捉えるのに十分な数、配置で温度計測装置を設置する必要がある。しかしながら、流体設備や建造内部等は複雑な形状をしていることが多く、温度計測装置を設置できない場所も存在する。また、流体系が高温である場合や高腐食性である場合には、温度計測装置を用いた温度の計測が制限される場合があり、温度分布を知るのに十分な数、配置で温度計測装置を設置できないことが多い。このように温度計測装置を設置できないことによる温度実測部位の不足を補うためには、計測された温度実測値から流体系全体の温度分布を推定し、補間する必要がある。 In order to directly know the temperature distribution of the fluid system flowing through the space separated by walls such as fluid facilities used in industrial processes and buildings, the number and arrangement sufficient to capture the characteristics of the temperature distribution in the fluid system It is necessary to install a temperature measuring device. However, fluid facilities and building interiors often have complicated shapes, and there are places where temperature measuring devices cannot be installed. In addition, when the fluid system is hot or highly corrosive, the temperature measurement using the temperature measurement device may be limited, and the temperature measurement with a sufficient number and arrangement to know the temperature distribution. Often equipment cannot be installed. Thus, in order to compensate for the shortage of the temperature measurement part due to the inability to install the temperature measurement device, it is necessary to estimate and interpolate the temperature distribution of the entire fluid system from the measured temperature measurement value.
温度分布を推定する対象物が均質な固体である場合、幾何学的に近い位置にある温度実測部位と温度推定点との間の温度相関が大きいことから、スプライン補間等のよく知られた補間法を用いて比較的容易に温度分布を推定し、補間することが可能である。例えば、実測部位と推定点との間の距離をもとに推定点の温度等を含む各種の値を推定する逆距離加重法と呼ばれる方法が知られている(例えば、非特許文献1を参照)。非特許文献1記載の方法は、実測部位の位置と推定点との間の距離を算出し、算出された距離が大きい実測部位の値ほど重みが小さくなるような重み付けをし、重み付き平均として推定点における値を推定する。具体的には、非特許文献1記載の方法は、距離の逆数(l/li)の累乗を重みとする次式(1)を用い、例えば温度推定点の温度を推定する。ここで、Telは、温度推定点における推定温度であり、liは、温度実測部位iの位置と温度推定点との間の距離であり、Tiは、温度実測部位iにおける温度実測値である。また、uは、正の値をとる補間パラメータである。
When the object whose temperature distribution is to be estimated is a homogeneous solid, a well-known interpolation such as spline interpolation is used because of the large temperature correlation between the temperature measurement site and the temperature estimation point that are close to the geometric position. The temperature distribution can be estimated and interpolated relatively easily using the method. For example, a method called an inverse distance weighting method for estimating various values including the temperature of the estimated point based on the distance between the actually measured part and the estimated point is known (for example, see Non-Patent Document 1). ). The method described in
これに対し、温度分布を推定する対象物が流体系である場合、対流による熱輸送が発生するため、温度実測部位と温度推定点とが幾何学的に近い位置にある場合であっても温度相関が大きいとは限らない。このような流体系の温度分布を推定する方法として、例えば、特許文献1には、温水の移動速度と貯湯槽内での実測温度の履歴とをもとに貯湯槽内の温度分布を推定する温度分布推定システムが開示されている。また、特許文献2には、航法装置で得た位置データである航跡および水温計により得た水温データとともに、潮流計で得た潮流ベクトルデータを重ねて表示する水温分布表示装置が開示されている。また、特許文献3には、屋内に設けたセンサーの実測値をもとに境界条件を構築し、指定箇所の温度や湿度、二酸化炭素濃度等の環境状態を熱伝導に関する式またはナビエストークスの式を用いて推定する空調用センサーシステムが開示されている。
On the other hand, if the object whose temperature distribution is to be estimated is a fluid system, heat transfer by convection occurs, so even if the measured temperature part and the estimated temperature point are close to the geometric position, the temperature The correlation is not always great. As a method for estimating the temperature distribution of such a fluid system, for example, in
一方、流体系の流れ場に関する技術として、例えば非特許文献2、非特許文献3、および特許文献4には、換気効率の指標の1つとして用いる吹出口および吸込口の勢力範囲の概念が記載されている。この非特許文献2、非特許文献3、および特許文献4に記載されている吹出口の勢力範囲は、複数の吹出口を備えた室内のある特定の点に注目したときに、検討対象とする吹出口からの気流がその点にどれだけ到達しているかを表す。また、吸込口の勢力範囲は、検討対象とする吸込口を通して排出される空気の室内各点での分布状態を表す。また、非特許文献2、非特許文献3、および特許文献4には、数値解析による勢力範囲の算出方法についても記載されている。
On the other hand, for example, Non-Patent
しかしながら、非特許文献1記載の逆距離加重法では、温度実測部位と温度推定点との間の距離liのみに基づいた重みによる重み付け平均によって例えば温度を推定するため、得られる温度推定結果には流体の流れの影響が反映されない。このため、実際の流体プロセス等の熱輸送に対する流動の寄与が非常に大きい流体系では、流速が大きい場合と小さい場合とで温度分布が大きく異なるにも関わらず、同じ温度分布が推定されてしまう。したがって、流体の流れによる熱輸送が支配的となる流体系への適用は困難である。
However, in the inverse distance weighting method described in
また、特許文献1記載のシステムは、1次元の流体流れにしか適用できず、3次元の流体流れ場を有する流体系への適用は困難である。また、特許文献2には、水温の分布と潮流の方向およびその速さとに相関があることが記載されてはいるものの、その相関を具体的にどのように求めて水温分布を推定するかは開示されておらず、流体系が必要とする非計測箇所の推定の精度が得られない。さらに、水温データと潮流ベクトルデータとは、2次元の海面上の場合のみについて想定されているため、3次元の流体流れ場を有する流体系への適用は困難である。
Further, the system described in
また、特許文献3記載のシステムのように、熱伝導に関する式またはナビエストークスの式を用いる方法では、空気が室内に流入する窓や空調設備の吹出口といった空気の流れの最上流位置の全てにセンサーを設けることが必須となる。このため、最上流位置の全てにセンサーを設置できない場合、流体系の全域で温度分布が推定できなくなる。したがって、流体の流入位置にセンサーを設置できないような流体系には適用が困難である。
In addition, as in the system described in
また、非特許文献2および非特許文献3記載の方法は、空気の吹出口や吸込口から流入出した流体の存在位置を可視化することのみに注目しており、温度推定への適用を想定していない。特許文献4記載の方法も同様に、空気齢の空間分布を算出することに着目しており、温度推定への適用を想定していない。
In addition, the methods described in
本発明は、上記に鑑みてなされたものであって、温度計測装置の配置に制約を与えることなく流体の流れによる熱輸送を考慮した高精度な温度推定を実現することができる流体系の温度推定方法、流体系の温度分布推定方法、流体系の温度分布モニタリング方法、および温度推定装置を提供することを目的とする。 The present invention has been made in view of the above, and the temperature of a fluid system capable of realizing highly accurate temperature estimation in consideration of heat transport due to the flow of the fluid without restricting the arrangement of the temperature measuring device. It is an object to provide an estimation method, a fluid system temperature distribution estimation method, a fluid system temperature distribution monitoring method, and a temperature estimation device.
また、本発明の他の目的は、表面欠陥がない溶融亜鉛めっき鋼板を製造可能な溶融亜鉛めっきポット内の溶融亜鉛温度制御方法、およびこの溶融亜鉛めっきポット内の溶融亜鉛温度制御方法を用いて製造した溶融亜鉛めっき鋼板を提供することにある。さらに、本発明の他の目的は、タンディッシュの耐火物損傷を抑制可能なタンディッシュ内の溶鋼温度制御方法を提供することにある。 Another object of the present invention is to use a hot dip galvanizing temperature control method in a hot dip galvanizing pot capable of producing a hot dip galvanized steel sheet having no surface defects, and a hot dip galvanizing temperature control method in the hot dip galvanizing pot. The object is to provide a manufactured hot-dip galvanized steel sheet. Furthermore, the other 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.
上述した課題を解決し、目的を達成するため、本発明にかかる流体系の温度推定方法は、温度既知領域が2箇所以上ある流体系の任意の温度推定点における温度を推定する流体系の温度推定方法であって、前記温度既知領域の位置情報と流体系全域における流体の流れを表す流体系の流れ場に関する情報とを用いて、温度既知領域を通過した、または温度既知領域内で生成した流体のうち、他の温度既知領域を通過することなく前記温度推定点まで到達した流体の、温度推定点の全流体中に占める比率を温度推定点における温度既知領域の勢力として取得する勢力取得工程と、各温度既知領域の温度と前記温度推定点における勢力とに関する情報を用いて、前記温度推定点における温度を推定する温度推定工程と、を含むことを特徴とする。 In order to solve the above-described problems and achieve the object, a temperature estimation method for a fluid system according to the present invention is a temperature of a fluid system that estimates a temperature at an arbitrary temperature estimation point of a fluid system having two or more known temperature regions. An estimation method, using the position information of the temperature known area and the information on the flow field of the fluid system representing the flow of the fluid in the entire fluid system, or passed through the temperature known area or generated in the temperature known area A power acquisition step of acquiring a ratio of the fluid estimated to reach the temperature estimation point without passing through another temperature known area of the fluid in the total fluid of the temperature estimation point as the power of the temperature known area at the temperature estimation point And a temperature estimation step of estimating the temperature at the temperature estimation point using information on the temperature of each temperature known region and the power at the temperature estimation point.
上述した課題を解決し、目的を達成するため、本発明にかかる流体系の温度分布推定方法は、温度分布を有する流体系の温度分布推定方法であって、上記の発明を用いて前記流体系の全域に設定した温度推定点の温度を推定し、前記各温度推定点について推定した温度を前記流体系の温度分布として推定することを特徴とする。 In order to solve the above-described problems and achieve the object, a temperature distribution estimation method for a fluid system according to the present invention is a temperature distribution estimation method for a fluid system having a temperature distribution. The temperature estimated at the temperature estimation points set over the entire area is estimated, and the temperature estimated for each temperature estimation point is estimated as the temperature distribution of the fluid system.
上述した課題を解決し、目的を達成するため、本発明にかかる流体系の温度分布モニタリング方法は、温度分布を有する流体系の温度分布モニタリング方法であって、上記の発明を用いて推定した前記流体系の温度分布をもとに、前記流体系の任意の断面における温度分布を可視化して画面表示することを特徴とする。 In order to solve the above-described problems and achieve the object, a temperature distribution monitoring method for a fluid system according to the present invention is a temperature distribution monitoring method for a fluid system having a temperature distribution, which is estimated using the above-described invention. Based on the temperature distribution of the fluid system, the temperature distribution in an arbitrary cross section of the fluid system is visualized and displayed on the screen.
上述した課題を解決し、目的を達成するため、本発明にかかる温度推定装置は、温度既知領域が2箇所以上ある流体系の任意の温度推定点における温度を推定する温度推定装置であって、前記温度既知領域の位置情報と流体系全域における流体の流れを表す流体系の流れ場に関する情報とを用いて、温度既知領域を通過した、または温度既知領域内で生成した流体のうち、他の温度既知領域を通過することなく前記温度推定点まで到達した流体の、温度推定点の全流体中に占める比率を温度推定点における温度既知領域の勢力として取得する勢力取得手段と、各温度既知領域の温度と前記温度推定点における勢力とに関する情報を用いて、前記温度推定点における温度を推定する温度推定手段と、を備えることを特徴とする。 In order to solve the above-described problems and achieve the object, a temperature estimation device according to the present invention is a temperature estimation device that estimates a temperature at an arbitrary temperature estimation point of a fluid system having two or more temperature known regions, Using the positional information of the temperature known region and the information on the flow field of the fluid system representing the flow of the fluid in the entire fluid system, other fluids that have passed through the temperature known region or generated in the temperature known region A force acquisition means for acquiring, as a force of the temperature known area at the temperature estimated point, a ratio of the fluid that has reached the temperature estimated point without passing through the temperature known area and occupies the total fluid of the temperature estimated point; and each temperature known area Temperature estimation means for estimating the temperature at the temperature estimation point using information relating to the temperature at the temperature estimation point.
上述した課題を解決し、目的を達成するため、本発明にかかる溶融亜鉛めっきポット内の溶融亜鉛温度制御方法は、本発明にかかる流体系の温度推定方法により推定した前記溶融亜鉛めっきポット内の溶融亜鉛温度データから、前記溶融亜鉛めっきポット内の所定の領域における溶融亜鉛の温度を抽出する温度抽出ステップと、抽出した温度が、所定の閾値範囲内にあるか否かを判定する判定ステップと、前記判定ステップにおいて、前記抽出した温度が閾値範囲外と判定された場合、前記抽出した温度が閾値範囲内となるよう前記溶融亜鉛めっきポットの加熱手段の出力を操作する制御ステップと、を含むことを特徴とする。 In order to solve the above-mentioned problems and achieve the object, a hot dip galvanizing temperature control method in a hot dip galvanizing pot according to the present invention includes a hot dip galvanizing pot estimated by the fluid system temperature estimating method according to the present invention. A temperature extraction step for extracting the temperature of the molten zinc in a predetermined region in the hot dip galvanizing pot from the hot dip zinc temperature data, and a determination step for determining whether or not the extracted temperature is within a predetermined threshold range; And a control step of operating the output of the heating means of the hot dip galvanizing pot so that the extracted temperature is within the threshold range when the extracted temperature is determined to be outside the threshold range in the determination step. It is characterized by that.
上述した課題を解決し、目的を達成するため、本発明にかかる溶湯亜鉛めっき鋼板は、本発明にかかる溶融亜鉛めっきポット内の溶融亜鉛温度制御方法を用いて製造したことを特徴とする。 In order to solve the above-described problems and achieve the object, the hot dip galvanized steel sheet according to the present invention is characterized by being manufactured using the hot dip galvanizing temperature control method in the hot dip galvanizing pot 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 is based on the molten steel temperature data in the tundish estimated by the fluid system temperature estimation method according to the present invention. In the temperature extraction step of extracting the temperature of the molten steel in the predetermined region in the tundish, the determination step of determining whether the extracted temperature is within a predetermined threshold range, and the extraction in the determination step And a control step of operating the output of the tundish heating means so that the extracted temperature is within the threshold range when the temperature is determined to be outside the threshold range.
本発明によれば、温度計測装置の配置に制約を与えることなく流体の流れによる熱輸送を考慮した高精度な温度推定を実現することができる。本発明によれば、表面欠陥がない溶融亜鉛めっき鋼板を製造することができる。本発明によれば、タンディッシュの耐火物損傷を抑制することができる。 ADVANTAGE OF THE INVENTION According to this invention, the highly accurate temperature estimation which considered the heat transport by the flow of a fluid can be implement | achieved, without giving restrictions to arrangement | positioning of a temperature measuring device. According to the present invention, a hot dip galvanized steel sheet having no surface defects can be produced. According to the present invention, refractory damage of a tundish can be suppressed.
以下、図面を参照して、本発明の流体系の温度推定方法、流体系の温度分布推定方法、流体系の温度分布モニタリング方法、および温度推定装置を実施するための形態について説明する。なお、この実施の形態によって本発明が限定されるものではない。また、図面の記載において、同一部分には同一の符号を付して示している。 DESCRIPTION OF EMBODIMENTS Hereinafter, embodiments of a fluid system temperature estimation method, a fluid system temperature distribution estimation method, a fluid system temperature distribution monitoring method, and a temperature estimation device according to the present invention will be described with reference to the drawings. Note that the present invention is not limited to the embodiments. Moreover, in description of drawing, the same code | symbol is attached | subjected and shown to the same part.
〔本発明の概念〕
図1は、本発明の概念を説明するための機能ブロック図である。図1に示すように、本発明は、温度既知領域が2箇所以上ある流体中の任意の温度推定点における温度を推定するものである。詳しくは、本発明は、2箇所以上の温度既知領域の位置情報である座標と流体系全域における流体の流れを表す流体系の流れ場に関する情報とを用いて温度推定点における各温度既知領域の勢力に関する情報を取得し、各温度既知領域の温度実測値(既知温度)と温度推定点における勢力とに関する情報を用いて流体中の任意の温度推定点における温度を推定する。各温度既知領域の勢力は、温度推定点における全流体のうち、温度既知領域から流れ場又は反転流れ場による移流拡散現象に従って流れてきた流体であって、且つ、温度既知領域から他の温度既知領域を通過することなく温度推定点まで到達した流体の比率(寄与率)のことを意味している。
[Concept of the present invention]
FIG. 1 is a functional block diagram for explaining the concept of the present invention. As shown in FIG. 1, the present invention estimates the temperature at an arbitrary temperature estimation point in a fluid having two or more known temperature regions. Specifically, the present invention uses coordinates that are positional information of two or more known temperature regions and information on the flow field of the fluid system that represents the fluid flow in the entire fluid system, and describes each temperature known region at the temperature estimation point. Information on the power is acquired, and the temperature at an arbitrary temperature estimation point in the fluid is estimated using information on the temperature measurement value (known temperature) in each temperature known region and the power at the temperature estimation point. The power of each temperature known region is the fluid that has flowed from the temperature known region according to the advection diffusion phenomenon by the flow field or the reverse flow field among all the fluids at the temperature estimation point, and other temperature known from the temperature known region. It means the ratio (contribution rate) of the fluid that has reached the temperature estimation point without passing through the region.
〔温度推定装置の構成〕
図2は、本発明を実施するための装置構成の一例を示すブロック図である。図2に示す温度推定装置1は、温度を推定する推定対象の流体系内の所定の温度実測部位に設置される1つ以上の温度計測装置2と接続される。温度推定装置1は、CPU、フラッシュメモリ等のROMやRAMといった各種ICメモリ、ハードディスク、各種 記憶媒体等の記憶装置、通信装置、表示装置や印刷装置等の出力装置、入力装置等を備えた公知のハードウェア構成で実現でき、例えばワークステーションやパソコン等の汎用コンピュータを用いることができる。
[Configuration of temperature estimation device]
FIG. 2 is a block diagram showing an example of an apparatus configuration for carrying out the present invention. A
この温度推定装置1は、温度既知領域を少なくとも2箇所含む流体系を推定対象とし、温度既知領域の温度(既知温度)をもとに流体系内の所定の温度推定点における温度を推定する。温度推定点の位置および数は適宜設定できる。代表的な温度既知領域としては、温度計測装置を配置して温度を直接測定した温度実測領域がある。その際、温度実測値が既知温度となる。図2は温度計測装置を用いた場合の装置構成図である。上記説明では温度実測領域を温度既知領域としたが、温度既知領域は温度が既知であればどのような領域でもよく、温度実測領域に限定されない。
This
温度場がほぼ定常とみなせる場合は、温度推定を行う時点付近の温度の瞬時値をそのまま既知温度としてかまわない。もし温度場が時間的に変化しうる場合は、既知温度と観測された時間を逐次保存した時系列温度データが必要となるので、温度既知領域において温度と時間とを対応させて逐次時系列で保存するようにする。そして、任意の時間における温度既知領域の温度を適宜抽出もしくは補間もしくは外挿して出力できるようにすると良い。温度が例えば固定である等既知であり、記憶装置に予め保存されている場合等、何らかの手段で取得が可能な場合には、この温度も温度既知領域の温度と同等に扱い、温度と時間とを対応させて時系列で保存しておき、任意の時間における温度の値を適宜抽出もしくは補間もしくは外挿して出力できるようにすると良い。 When the temperature field can be regarded as almost steady, the instantaneous value of the temperature near the time point at which the temperature is estimated may be used as the known temperature as it is. If the temperature field can change over time, time series temperature data that sequentially stores the known temperature and the observed time is required. Try to save. Then, it is preferable that the temperature in the temperature known region at an arbitrary time can be appropriately extracted, interpolated, or extrapolated and output. If the temperature is known, for example, fixed, and stored in a storage device in advance, it can be obtained by some means. Are stored in time series so that the temperature value at an arbitrary time can be appropriately extracted, interpolated, or extrapolated for output.
温度既知領域はその領域の熱流体的な特性によって温度実測領域、発吸熱領域、および流入出領域に分類して考えることができる。温度実測領域は温度実測部位、もしくは温度実測部位とその近傍まで含めた領域を指す。発吸熱領域は発吸熱部位、もしくは発吸熱部位とその近傍まで含めた領域を指す。流入出領域は流入出部位、もしくは流入出部位とその近傍まで含めた領域を指す。温度実測部位とは、実際に温度を計測するなどの手段により温度が既知となっており、その部位またはその近傍まで含めた領域で流体の流入出や発熱・吸熱が起こっていない流体系内の点、面、または領域のことをいう。温度実測部位は必ずしも直接温度を計測している部位に限定されない。例えばモデル式等で他のパラメータから温度を換算することができる部位、制御装置等により温度を制御している部位等の温度が間接的に既知となっている部位も含まれる。流体系内の温度実測部位の位置や数、例えば温度計測装置2の設置位置や数は、適宜設定できる。発吸熱部位とは、発熱または吸熱が生じている流体系内の点、面、または領域のことをいう。流入出部位とは、系内への流体の流入または系外への流体の流出が生じている点、面、または領域のことをいう。発吸熱部位および流入出部位に関しては、温度が未知の部位も含めることができ、これにより温度推定の信頼性を向上させることも可能である。
The temperature known region can be considered by classifying it into a temperature measurement region, a heat generating / absorbing region, and an inflow / outflow region according to the thermohydrodynamic characteristics of the region. The temperature measurement region indicates a temperature measurement region or a region including the temperature measurement region and its vicinity. An endothermic / endothermic region refers to an endothermic / absorptive site, or a region that includes the endothermic / endothermic site and its vicinity. The inflow / outflow region indicates an inflow / outflow region or a region including the inflow / outflow region and its vicinity. The temperature measurement part is a part of the fluid system where the temperature is already known by means such as actually measuring the temperature, and fluid is not flowing in or out of the part or the area including the vicinity. A point, surface, or area. The temperature measurement part is not necessarily limited to the part that directly measures the temperature. For example, a part where the temperature can be converted from other parameters by a model formula or the like, a part where the temperature is indirectly known, such as a part where the temperature is controlled by a control device or the like are included. The position and number of the temperature measurement sites in the fluid system, for example, the installation position and number of the
〔温度の推定原理〕
先ず、温度推定点における温度の推定原理について説明する。なお、以下の説明では、推定対象の流体系が、K箇所の温度実測部位i(i=1〜K)と、L箇所の発熱部位または吸熱部位である発吸熱部位i(i=K+1〜K+L)と、M箇所の流入部位または流出部位である流入出部位i(i=K+L+1〜K+L+M)とを含むこととし、この流体系内に設定されるN箇所の温度推定点j(j=1〜N)の温度を推定することとする。
[Principle of temperature estimation]
First, the temperature estimation principle at the temperature estimation point will be described. In the following description, the estimation target fluid system includes K temperature measurement sites i (i = 1 to K) and L heat generation sites or heat generation sites i (i = K + 1 to K + L). ) And inflow / outflow sites i (i = K + L + 1 to K + L + M) which are M inflow sites or outflow sites, and N temperature estimation points j (j = 1 to 1) set in the fluid system The temperature of N) will be estimated.
本実施の形態の温度推定装置1は、流体系の流れ場を用い、温度推定点jにおける流体の勢力、詳細には、温度推定点jにおける流体の温度実測部位i、発吸熱部位i、および流入出部位iの勢力を取得する。そして、温度推定装置1は、この流体の勢力を指標値として用いることで、その流れ場のもとでの熱の移流拡散を考慮した温度推定点jの温度を推定する。以下では、温度推定装置1は、前述の勢力として、下流側勢力R1ijおよび上流側勢力R2ijの2種類の勢力(R1ij,R2ij)を取得し、温度推定の指標値として用いる。
The
図3−1〜図3−7を参照して、下流側勢力および上流側勢力の概念について説明する。図3−1は流体系の一例であり、流体100、容器101、および仕切り板102で構成されている。図3−1に示す流体系内には流れ場Fがあり、図中の矢印の向きに循環する流れを形成しているとする。簡単のため、図3−1の流体系には外部からの流入、外部への流出、化学反応等による流体の生成や消失は起こらないものとする。図3−1に示す流体系には3箇所の温度既知領域R1、R2、R3があり、それぞれ丸で図示している。このとき、流体系内の任意の温度推定点における温度既知領域R1の下流側勢力の定義は次の通りである。
The concept of the downstream force and the upstream force will be described with reference to FIGS. FIG. 3A is an example of a fluid system, and includes a fluid 100, a
すなわち、温度推定点における温度既知領域R1の下流側勢力は、温度推定点における全流体のうち、温度既知領域R1から流れ場Fによる移流拡散現象に従って流れてきた流体であり、温度既知領域R1から他の温度既知領域、本例では温度既知領域R2および温度既知領域R3を通過することなく温度推定点まで到達した流体の比率と定義される。この定義によって温度既知領域R1の下流側勢力は、流体系内の全ての場所に対して算出することができる。同様に、温度既知領域R2の下流側勢力および温度既知領域R3の下流側勢力も算出できる。この結果、温度既知領域R1、R2、R3に対する下流側勢力分布I11,I12、I13はそれぞれ図3−2、図3−3、および図3−4に示すようになる。図3−2、図3−3、および図3−4から明らかなように、温度既知領域R1、R2、R3に対する下流側勢力分布I11、I12、I13は、温度既知領域R1、R2、R3の各領域から流れに沿って下流側に伸びた分布となる。また、途中に別の温度既知領域があると、その領域を避けた形の分布になる。温度既知領域R1、R2、R3の下流側勢力分布で示された領域はそれぞれ温度既知領域R1、R2、R3から流れてきた流体が多く含まれている領域であるため、それぞれ温度既知領域R1、R2、R3の既知温度と強い温度相関を持つ。 That is, the downstream side force of the temperature known region R1 at the temperature estimation point is a fluid that has flowed from the temperature known region R1 according to the advection diffusion phenomenon by the flow field F out of all the fluids at the temperature estimation point, and from the temperature known region R1. It is defined as the ratio of the fluid that has reached the temperature estimation point without passing through another temperature known region, in this example, the temperature known region R2 and the temperature known region R3. With this definition, the downstream force of the temperature known region R1 can be calculated for all locations in the fluid system. Similarly, the downstream force of the temperature known region R2 and the downstream force of the temperature known region R3 can also be calculated. As a result, the downstream power distributions I11, I12, and I13 for the temperature known regions R1, R2, and R3 are as shown in FIGS. 3-2, 3-3, and 3-4, respectively. As is apparent from FIGS. 3-2, 3-3, and 3-4, the downstream power distributions I11, I12, and I13 with respect to the temperature known regions R1, R2, and R3 correspond to the temperature known regions R1, R2, and R3. The distribution extends from each region to the downstream side along the flow. Further, if there is another temperature known area in the middle, the distribution is in a form avoiding that area. The regions indicated by the downstream side power distribution of the temperature known regions R1, R2, R3 are regions containing a large amount of fluid flowing from the temperature known regions R1, R2, R3, respectively. It has a strong temperature correlation with the known temperatures of R2 and R3.
一方、温度推定点における温度既知領域R1の上流側勢力の定義は次の通りである。まず、上記流れ場Fに対して、流速ベクトルの大きさは同じで、方向のみを全て反転させた流れ場(本明細書では、反転流れ場と呼ぶ)を取得する。そして、温度推定点における全流体のうち、温度既知領域R1から反転流れ場による移流拡散現象に従って流れてきた流体であり、温度既知領域R1から他の温度既知領域、本例では温度既知領域R2および温度既知領域R3を通過することなく推定点まで到達した流体の比率を温度推定点における温度既知領域R1の上流側勢力と定義する。上記定義により温度既知領域R1の上流側勢力は流体系内の全ての場所に対して算出することができる。同様にして、温度既知領域R2の上流側勢力および温度既知領域R3の上流側勢力を算出できる。この結果、温度既知領域R1、R2、R3に対する上流側勢力分布I21,I22、I23はそれぞれ図3−5、図3−6、および図3−7に示すようになる。図3−5、図3−6、および図3−7から明らかなように、温度既知領域R1、R2、R3に対する上流側勢力分布は、温度既知領域R1、R2、R3の各領域から流れとは反対向きの上流側に伸びた分布となる。また、途中に別の温度既知領域があると、その領域を避けた形の分布になる。温度既知領域R1、R2、R3の上流側勢力分布で示された領域にある流体の多くはそれぞれ温度既知領域R1、R2、R3へ流れていくことになるため、温度既知領域R1、R2、R3の上流側勢力分布が示す領域の温度と温度既知領域R1、R2、R3の既知温度とは強い温度相関を持つ。 On the other hand, the definition of the upstream power of the temperature known region R1 at the temperature estimation point is as follows. First, with respect to the flow field F, a flow field having the same velocity vector magnitude and having only the direction reversed (referred to as a reversed flow field in this specification) is acquired. Of all the fluids at the temperature estimation point, the fluid has flowed from the temperature known region R1 according to the advection diffusion phenomenon due to the reverse flow field, and from the temperature known region R1 to another temperature known region, in this example, the temperature known region R2 and The ratio of the fluid that has reached the estimated point without passing through the known temperature region R3 is defined as the upstream side force of the known temperature region R1 at the estimated temperature point. Based on the above definition, the upstream force of the temperature known region R1 can be calculated for all locations in the fluid system. Similarly, the upstream side force of the temperature known region R2 and the upstream side force of the temperature known region R3 can be calculated. As a result, the upstream side force distributions I21, I22, and I23 with respect to the temperature known regions R1, R2, and R3 are as shown in FIGS. 3-5, 3-6, and 3-7, respectively. As is clear from FIGS. 3-5, 3-6, and 3-7, the upstream side force distribution with respect to the temperature known regions R1, R2, and R3 flows from each region of the temperature known regions R1, R2, and R3. Has a distribution extending upstream in the opposite direction. Further, if there is another temperature known area in the middle, the distribution is in a form avoiding that area. Since most of the fluid in the region indicated by the upstream side force distribution of the temperature known regions R1, R2, R3 flows to the temperature known regions R1, R2, R3, respectively, the temperature known regions R1, R2, R3 The temperature in the region indicated by the upstream side power distribution and the known temperatures in the temperature known regions R1, R2, and R3 have a strong temperature correlation.
下流側勢力R1ijおよび上流側勢力R2ijを取得するため、本実施の形態では、K箇所の温度実測部位i、L箇所の発吸熱部位i、およびM箇所の流入出部位iの各部位iに対して、各部位iをそれぞれ包含し、かつお互いが重複しないような有限の領域i(i=1〜K+L+M)を設定する。具体的には、温度実測部位iに対応する領域iとして温度実測領域iを設定し、発吸熱部位iに対応する領域iとして発吸熱領域iを設定し、流入出部位iに対応する領域iとして流入出領域iを設定する。設定する温度実測領域i、発吸熱領域i、および流入出領域iの形状は、その温度実測部位i、発吸熱部位i、または流入出部位iを包含していればどのような形状でも構わない。すなわち、各領域iは、例えば点や線、面であってもよいし、3次元の有限な体積を持つ領域としてもよい。 In order to obtain the downstream side force R 1ij and the upstream side force R 2ij , in this embodiment, each part i of the K temperature measurement part i, the L heat generation / absorption part i, and the M inflow / outflow part i. In contrast, a finite region i (i = 1 to K + L + M) that includes each part i and does not overlap each other is set. Specifically, the temperature measurement area i is set as the area i corresponding to the temperature measurement area i, the heat generation / absorption area i is set as the area i corresponding to the heat generation / absorption area i, and the area i corresponding to the inflow / outflow area i is set. As shown in FIG. The shape of the temperature measurement region i, the heat generation / absorption region i, and the inflow / outflow region i to be set may be any shape as long as it includes the temperature measurement region i, the heat generation / absorption region i, or the inflow / outflow region i. . That is, each region i may be, for example, a point, a line, or a surface, or may be a region having a three-dimensional finite volume.
幅広い対象に適用できるようにするためには、以下の方法を用いて温度実測部位i、発吸熱部位i、および流入出部位iに対応する温度実測領域i、発吸熱領域i、および流入出領域iを設定するとよい。例えば、温度実測部位iは、流体系内に設置される温度計測装置2の設置位置である。このため、温度実測領域iは、この温度実測部位iである温度計測装置2の設置位置を中心とした半径rの球領域を設定し、温度実測領域iとするとよい。設定する球領域の半径rを大きい値とすると推定される流体系の温度分布は急峻な温度分布になり、小さい値とすると平滑化された温度分布となる。具体的な半径rの値は、流体系の流動特性によって最適値が異なる。例えば1辺が1mの水槽を適用対象とし、この水槽内に6個の温度計測装置2を適所に設置して温度実測部位iを6箇所設ける場合、温度実測領域iは、例えば半径r=0.05m程度の球領域とするのが好ましい。この球領域の半径rは、温度実測領域i、発吸熱領域i、流入出領域iの各領域iが重複しない限りにおいて、どのような値としてもよいが、各温度実測部位iについてそれぞれ設定する温度実測領域iの半径は、全て同じにするのが好ましい。
In order to be applicable to a wide range of objects, the temperature measurement region i, the heat generation / absorption region i, and the temperature measurement region i corresponding to the inflow / outflow region i, the heat generation / absorption region i, and the inflow / outflow region using the following method i may be set. For example, the temperature measurement site i is an installation position of the
流体系内に加熱装置による加熱や吸熱装置による吸熱、化学反応等によって発熱または吸熱が生じる領域が含まれる場合、その領域が発吸熱部位iとなる。この場合には、この領域を発吸熱領域iとする。例えば、発熱や吸熱が流体系の端、具体的には、例えば推定対象の流体系の流れを区画する設備等の壁面や、推定対象の流体系の浴面で生じる場合は、この壁面や浴面を発吸熱領域iとする。また、例えば、流体系内に浸漬された固体が発熱し、または吸熱する場合、化学反応等によって発熱または吸熱を生じる物質が流体系内に浸漬される場合は、この固体の表面を発吸熱領域iとする。また、発熱や吸熱が流体系内の一部の領域で生じる場合、例えば誘導加熱装置が流体系内に設けられる場合には、加熱エネルギーが印加される流体系内の領域を発吸熱領域iとする。 When the fluid system includes a region where heat is generated or endothermic due to heating by a heating device, heat absorption by a heat absorption device, chemical reaction, or the like, that region becomes a heat generation / heat absorption site i. In this case, this region is defined as a heat generation / absorption region i. For example, if heat generation or heat absorption occurs at the end of the fluid system, specifically, for example, on the wall surface of a facility or the like that divides the flow of the estimation target fluid system or the bath surface of the estimation target fluid system, The surface is defined as a heat generation / absorption region i. In addition, for example, when a solid immersed in a fluid system generates heat or absorbs heat, a substance that generates heat or endotherm due to a chemical reaction or the like is immersed in the fluid system. i. Further, when heat generation or heat absorption occurs in a part of the fluid system, for example, when an induction heating device is provided in the fluid system, the region in the fluid system to which heating energy is applied is referred to as a heat generation / heat absorption region i. To do.
系内への流体の流入や系外への流体の流出が存在する場合、その流入する領域や流出する領域が流入出部位iとなる。この場合には、この領域を流入出領域iとする。例えば、推定対象の流体系の流れを区画する境界面から流体が流入し、またはこの境界面から流体が流出する場合は、該当する境界面を流入出領域iとする。 When there is an inflow of fluid into the system or an outflow of fluid outside the system, the inflow region or the outflow region is the inflow / outflow site i. In this case, this region is referred to as an inflow / outflow region i. For example, when the fluid flows in from the boundary surface that divides the flow of the estimation target fluid system, or the fluid flows out from the boundary surface, the corresponding boundary surface is set as the inflow / outflow region i.
ただし、温度実測領域i、発吸熱領域i、および流入出領域iとする領域i,i´同士(i=1〜K+L+M,i´=1〜K+L+M,i≠i´)が重複すると、下流側勢力R1ijおよび上流側勢力R2ijを取得できなくなる。このため、各領域iが必ず温度実測領域i、発吸熱領域i、および流入出領域iのいずれか1つに属するように、すなわち、各領域iが重複しないように、温度実測領域i、発吸熱領域i、および流入出領域iの形状および大きさを決定する必要がある。 However, if the temperature measurement region i, the heat generation / absorption region i, and the regions i and i ′ (i = 1 to K + L + M, i ′ = 1 to K + L + M, i ≠ i ′) as the inflow / outflow region i overlap, the downstream side The power R 1ij and the upstream power R 2ij cannot be acquired. For this reason, each of the temperature measurement areas i, the heat generation areas i, and the heat generation / absorption areas i, and the inflow / outflow areas i, so that each area i does not overlap. It is necessary to determine the shape and size of the heat absorption region i and the inflow / outflow region i.
注目する1つの領域(注目領域)iを通過した、または注目領域i内で生成した流体であって、かつ他の領域i´を通過することなく温度推定点jまで到達した流体を、温度推定点jにおける注目領域iの流体成分と定義する。そして、温度推定点jにおける全流体に対する注目領域iの流体成分の比率を、温度推定点jにおける該当する部位iの勢力と定義し、流体系の流れ場(以下、「実際の流れ場」と呼ぶ。)を用いて取得した勢力を温度推定点jにおける該当する部位iの下流側勢力R1ij、流体系の反転流れ場を用いて取得した勢力を温度推定点jにおける該当する部位iの上流側勢力R2ijと定義する。 Temperature estimation is performed on a fluid that has passed through one region of interest (region of interest) i or that has been generated within the region of interest i and that has reached the temperature estimation point j without passing through another region i ′. It is defined as the fluid component of the region of interest i at point j. Then, the ratio of the fluid component of the region of interest i to the total fluid at the temperature estimation point j is defined as the force of the corresponding part i at the temperature estimation point j, and the flow field of the fluid system (hereinafter referred to as “actual flow field”). The downstream force R 1ij of the corresponding part i at the temperature estimation point j and the power acquired using the reversal flow field of the fluid system are upstream of the corresponding part i at the temperature estimation point j. It is defined as side force R 2ij .
すなわち、下流側勢力R1ijは、流体系の流れ場(以下、「実際の流れ場」と呼ぶ。)を用いて取得する。実際の流れ場は、例えば数値シミュレーションや、実機、実機を模擬した実験装置等を用いて計算する。例えば、推定対象の流体系全域の流速ベクトル、具体的には、流体系全域を同一サイズで区画した各領域における流体の向きおよび流速を表す流速ベクトルを求め、実際の流れ場とする。 That is, the downstream side force R 1ij is acquired using a fluid flow field (hereinafter referred to as “actual flow field”). The actual flow field is calculated using, for example, a numerical simulation, an actual machine, an experimental apparatus that simulates the actual machine, and the like. For example, a flow velocity vector over the entire fluid system to be estimated, specifically, a flow velocity vector representing the direction and flow velocity of the fluid in each region where the entire fluid system is partitioned with the same size, is obtained as an actual flow field.
この実際の流れ場を用い、温度推定点jにおける流体成分の比率を全ての領域iについて算出し、下流側勢力R1ijとして取得する。具体的には、温度実測領域i、発吸熱領域i、および流入出領域iの各領域iについて、該当する領域iを通過した、またはこの領域i内で生成した流体であり、かつ他の温度実測領域i´、発吸熱領域i´、および流入出領域i´を通過することなく温度推定点jまで到達した流体(流体成分)の、温度推定点jの全流体に対する比率を算出し、下流側勢力R1ijとする。 Using this actual flow field, the ratio of the fluid component at the temperature estimation point j is calculated for all the regions i, and is acquired as the downstream force R 1ij . Specifically, each region i of the temperature measurement region i, the heat generation / absorption region i, and the inflow / outflow region i is a fluid that has passed through the corresponding region i or generated in this region i, and has other temperatures. The ratio of the fluid (fluid component) that has reached the temperature estimation point j without passing through the measurement region i ′, the heat generation / absorption region i ′, and the inflow / outflow region i ′ to the total fluid at the temperature estimation point j is calculated downstream. The side force R 1ij is assumed.
一方、上流側勢力R2ijは、流体系の反転流れ場を用いて取得する。この反転流れ場は、実際の流れ場として求めた流速ベクトルを全て反転させることで得られる。そして、この反転流れ場を用い、温度推定点jにおける流体成分の比率を全ての領域iについて算出し、上流側勢力R2ijとして取得する。具体的には、温度実測領域iおよび流入出領域iの各領域iについて、該当する領域iを通過した、またはこの領域i内で生成した流体であり、かつ他の温度実測領域i´、発吸熱領域i´、または流入出領域i´を通過することなく温度推定点jまで到達した流体(流体成分)の温度推定点jの全流体に対する比率を算出し、上流側勢力R2ijとする。また、発吸熱領域iについては、流体成分の比率を「0」とし、上流側勢力R2ijとする(上流側勢力R2ij=0とする)。発吸熱部位iとこの発吸熱部位iの下流側の位置とには温度相関があるものの、発吸熱部位iとこの発吸熱部位iの上流側の位置とには温度相関がないためである。 On the other hand, the upstream side force R 2ij is acquired by using the reverse flow field of the fluid system. This inversion flow field can be obtained by inverting all the flow velocity vectors obtained as the actual flow field. Then, using this reversal flow field, the ratio of the fluid component at the temperature estimation point j is calculated for all the regions i, and is acquired as the upstream force R 2ij . Specifically, each region i of the temperature measurement region i and the inflow / outflow region i is a fluid that has passed through the corresponding region i or generated in this region i, and other temperature measurement regions i ′, The ratio of the fluid (fluid component) reaching the temperature estimation point j without passing through the heat absorption region i ′ or the inflow / outflow region i ′ with respect to all the fluids at the temperature estimation point j is calculated as an upstream force R 2ij . For the heat generation / absorption region i, the ratio of the fluid component is set to “0”, and the upstream force R 2ij is set (upstream force R 2ij = 0). This is because although there is a temperature correlation between the endothermic part i and the position downstream of the endothermic part i, there is no temperature correlation between the endothermic part i and the upstream position of the endothermic part i.
その後、以上のようにして取得した下流側勢力R1ijおよび上流側勢力R2ijの対である勢力(R1ij,R2ij)をもとに、温度推定点jに対する部位i毎の重みWij、詳細には、各部位iの既知温度に対して重み付けを行うための重みWijを算出する。例えば、単調非減少関数となる重み関数W(R1,R2)を用い、重みWijをWij=W(R1ij,R2ij)として算出する。 Then, above manner downstream force acquired by R 1ij and upstream force R force is a pair of 2ij (R 1ij, R 2ij) based on the weight W ij of each region i for temperature estimation point j, Specifically, a weight W ij for weighting the known temperature of each part i is calculated. For example, a weight function W (R 1 , R 2 ) that is a monotonous non-decreasing function is used, and the weight W ij is calculated as W ij = W (R 1ij , R 2ij ).
なお、非特許文献2では吹出口および吸込口の勢力範囲が定義されている。これは本発明で述べている下流側勢力および上流側勢力と類似しているが、異なった概念である。すなわち、非特許文献2の吹出口および吸込口の勢力範囲は、流体系の流入部位と流出部位、つまり境界条件として設定できる部位に対してのみ適用できる方法である。従って、実測部位や発吸熱部位に対しては吹出口および吸込口の勢力範囲は定義できない。これに対して、本発明で新たに考案した下流側勢力および上流側勢力は、境界のみならず流体系内部に存在する実測部位や発吸熱部位に対しても定義できる。本発明では実測部位や発吸熱部位を考慮することが極めて重要であり、本発明で新たに考案した下流側勢力および上流側勢力の概念を使うことが必須である。
In
流れ場がほぼ定常であるとみなせる場合は、その後、以上のようにして取得した下流側勢力R1ijおよび上流側勢力R2ijの対である勢力(R1ij,R2ij)をもとに、温度推定点jに対する部位i毎の重みWij、詳細には、各部位iの既知温度に対して重み付けを行うための重みWijを算出する。例えば、単調非減少関数となる重み関数W(R1ij,R2ij)を用い、重みWijをWij=W(R1ij,R2ij)として算出する。 If the flow field can be regarded as being substantially constant, then the above manner to a pair of the obtained downstream force R 1ij and upstream force R 2ij in force (R 1ij, R 2ij) on the basis of the temperature weight W ij of each region i for estimating point j, in particular, calculates the weight W ij for weighting relative to the known temperature of each portion i. For example, a weight function W (R 1ij , R 2ij ) that is a monotonic non-decreasing function is used, and the weight W ij is calculated as W ij = W (R 1ij , R 2ij ).
そして、温度実測部位i、発吸熱部位i、および流入出部位iの各部位iの既知温度Tiと、算出した温度推定点jに対する部位i毎の重みWijとを用いた重み付き平均によって、温度推定点jにおける推定温度を算出する。各部位iの既知温度Tiと、温度推定点jの推定温度Tejとの関係は、次式(2)によって表される。温度実測部位iの既知温度Tiは、温度計測装置2によって計測される温度実測値である。発吸熱部位iおよび流入出部位iの既知温度Tiは、該当する部位iの温度が既知の場合にその値を用いる。発吸熱部位iや流入出部位iの温度が未知である場合には、該当する部位iに対する重みWijの値を「0」に置き換えた上で、次式(2)に従って温度推定点jの推定温度Tejを算出する。
Then, by a weighted average using the known temperature T i of each part i of the temperature actually measured part i, the heat generating and absorbing part i, and the inflow / outflow part i, and the weight W ij for each part i with respect to the calculated temperature estimation point j The estimated temperature at the temperature estimation point j is calculated. Relationship between the known temperature T i, the estimated temperature Te j temperature estimate point j of each part i is expressed by the following equation (2). The known temperature T i of the temperature measurement site i is a temperature measurement value measured by the
一方、本発明は流れ場がほぼ定常であるとみなせるならば、温度場が時間的に変化しうる場合でも利用できる。この場合は、上記勢力とともに後述する伝達時間を取得し、時系列温度データを用いて温度を推定する。時系列温度データは温度実測部位i、発吸熱部位i、および流入出部位iの各部位iにおいて観測された実測温度もしくは他の手段で確認された既知の温度と、観測された時間とを対応させて逐次記録したものであり、任意の時間tにおける温度Ti(t)を実測温度と時間との記録から補間、外挿して出力できるものである。計測器故障などにより時間tにおいて実測温度が観測されていない場合は温度が未知であるとして扱っても良いが、他の時間において実測温度が観測されている場合は、近傍の時間におけるデータを補間、外挿した温度をTi(t)として出力するようにしても良い。 On the other hand, the present invention can be used even when the temperature field can be changed with time if the flow field can be considered to be almost steady. In this case, a transmission time, which will be described later, is acquired together with the power, and the temperature is estimated using time-series temperature data. The time-series temperature data corresponds to the observed time and the actually measured temperature observed at each site i of the temperature measured site i, the heat-exothermic site i, and the inflow / outflow site i or other means. The temperature T i (t) at an arbitrary time t can be interpolated and extrapolated from the recorded actual temperature and time and output. If the measured temperature is not observed at time t due to instrument failure, etc., the temperature may be treated as unknown, but if the measured temperature is observed at other times, the data at a nearby time is interpolated. The extrapolated temperature may be output as T i (t).
伝達時間は、各部位iと温度推定点jとの間を流体が移流拡散によって移動するのに要する時間である。具体的には、伝達時間は下流側伝達時間と上流側伝達時間とからなり、流体が各部位iから温度推定点jへ移動するのに要する時間が下流側伝達時間τ1ij、流体が温度推定点jから各部位iへ移動するのに要する時間が上流側伝達時間τ2ijとなる。そして、時系列温度データと温度を推定したい時間t0と伝達時間(τ1ij、τ2ij)を用いて、温度を推定したい時間から伝達時間分だけ過去もしくは未来の時点における温度を算出して既知温度とする。具体的には、温度を推定したい時間t0を基準にして下流側伝達時間τ1ij分だけ過去の時点で観測された部位iの温度を、温度推定点jにおける部位iの下流側既知温度とする。すなわち、時系列温度データから時間t0−τ1ijにおける実測温度Ti(t0−τ1ij)を出力し、下流側既知温度とすればよい。同様に、温度を推定したい時間t0を基準にして上流側伝達時間τ2ij分だけ未来の時点で観測された部位iの温度を温度推定点jにおける部位iの上流側既知温度とする。すなわち、時系列温度データから時間t0+τ2ijにおける実測温度Ti(t0+τ2ij)を出力し、上流側既知温度とすればよい。 The transmission time is a time required for the fluid to move between each part i and the temperature estimation point j by advection diffusion. Specifically, the transmission time includes a downstream transmission time and an upstream transmission time. The time required for the fluid to move from each part i to the temperature estimation point j is the downstream transmission time τ 1ij , and the fluid is the temperature estimation. The time required to move from the point j to each part i is the upstream transmission time τ 2ij . Then, using the time series temperature data, the time t 0 at which the temperature is to be estimated and the transmission time (τ 1ij , τ 2ij ), the temperature at the past or future time is calculated by the transmission time from the time at which the temperature is to be estimated. Let it be temperature. Specifically, the temperature of the part i observed at the past time point by the downstream transmission time τ 1ij with respect to the time t 0 at which the temperature is to be estimated is set as the downstream known temperature of the part i at the temperature estimation point j. To do. That is, the measured temperature T i (t 0 -τ 1ij ) at the time t 0 -τ 1ij may be output from the time series temperature data to be the downstream known temperature. Similarly, the temperature of the part i observed at the future time point by the upstream transmission time τ 2ij with respect to the time t 0 at which the temperature is to be estimated is set as the upstream known temperature of the part i at the temperature estimation point j. That is, the measured temperature T i (t 0 + τ 2ij ) at the time t 0 + τ 2ij may be output from the time series temperature data to be the upstream known temperature.
下流側勢力R1ijおよび上流側勢力R2ijの対である勢力(R1ij,R2ij)をもとに、下流側重みW1ijと上流側重みW2ij、詳細には、各部位iの下流側既知温度Ti(t0−τ1ij)、上流側既知温度Ti(t0+τ2ij)それぞれに対応する重みW1ij、W2ijを算出する。下流側重みW1ijに対しては任意の上流側勢力R2ijに対して下流側勢力R1ijの単調非減少関数となる重み関数W1(R1ij,R2ij)を用いて重みW1ijをW1ij=W1(R1ij,R2ij)とし、上流側重みW2ijに対しては任意の下流側勢力R1ijに対して上流側勢力R2の単調非減少関数となる重み関数W2(R1ij,R2ij)を用いて重みW2ijをW2ij=W2(R1ij,R2ij)として算出すればよい。 Based on the force (R 1ij , R 2ij ) which is a pair of the downstream force R 1ij and the upstream force R 2ij , the downstream weight W 1ij and the upstream weight W 2ij , specifically, the downstream side of each part i Weights W 1ij and W 2ij corresponding to the known temperature T i (t 0 −τ 1ij ) and the upstream known temperature T i (t 0 + τ 2ij ) are calculated. Downstream the weight W 1ij weighting function W 1 including a monotone non-decreasing function of the downstream force R 1ij for any upstream force R 2ij for (R 1ij, R 2ij) the weight W 1ij with W 1ij = W 1 (R 1ij , R 2ij ), and for the upstream weight W 2ij , a weight function W 2 (R that is a monotonic non-decreasing function of the upstream force R 2 with respect to an arbitrary downstream force R 1ij 1 ij , R 2ij ) and the weight W 2ij may be calculated as W 2ij = W 2 (R 1ij , R 2ij ).
温度実測部位i、発吸熱部位i、および流入出部位iの各部位iの下流側既知温度Ti(t0−τ1ij)と温度推定点jに対する部位i毎の下流側重みW1ij、および上流側既知温度Ti(t0+τ2ij)と上流側重みW2ijを用いた重み付き平均によって、時間t0における温度推定点jの推定温度を算出する。時間t0における温度推定点jの推定温度Tej(t0)は次式(3)によって表される。 The downstream actual temperature T i (t 0 -τ 1ij ) and the downstream weight W 1ij for each part i with respect to the temperature estimation point j of each part i of the temperature measurement part i, the heat generation / absorption part i, and the inflow / outflow part i, and The estimated temperature of the temperature estimation point j at time t 0 is calculated by a weighted average using the upstream known temperature T i (t 0 + τ 2ij ) and the upstream weight W 2ij . The estimated temperature Te j (t 0 ) at the temperature estimation point j at time t 0 is expressed by the following equation (3).
発吸熱部位iおよび流入出部位iの下流側既知温度Ti(t0−τ1ij)、上流側既知温度Ti(t0+τ2ij)は、該当する部位iの温度が既知の場合にその値を用いる。発吸熱部位iや流入出部位iの温度が未知である場合には、該当する部位iに対する下流側重みW1ijと上流側重みW2ijの値を「0」に置き換えた上で、上記式(3)に従って温度推定点jの推定温度Tej(t0)を算出する。 The downstream side known temperature T i (t 0 -τ 1ij ) and the upstream side known temperature T i (t 0 + τ 2ij ) of the endothermic part i and the inflow / outlet part i are obtained when the temperature of the corresponding part i is known. Use the value. Originating when the temperature of the endothermic sites i and flows out site i is unknown, after replacing the value of the downstream weight W 1ij upstream weighting W 2ij for the corresponding site i to "0", the formula ( According to 3), the estimated temperature Te j (t 0 ) of the temperature estimation point j is calculated.
次に、温度推定装置1が行う処理手順について図4〜図12を参照して説明する。対象の流体系の温度場がほぼ定常とみなせる場合は図4〜図7の処理手順を用い、対象の流体系の温度場が時間変化しうる場合は図8〜図12の処理手順を用いる。温度推定装置1は、図4〜図7または図8〜図12に示す処理手順に従って処理を行うことで流体系の温度推定方法、流体系の温度分布推定方法、および流体系の温度分布モニタリング方法を実施する。ここで説明する処理は、この処理を実現するためのプログラムを例えば温度推定装置1の記憶装置に保存しておき、このプログラムを読み出して実行することで実現できる。
Next, a processing procedure performed by the
最初に対象の流体系の温度場がほぼ定常とみなせる場合の処理手順について図4〜図7を参照して説明する。先ず、温度推定装置1が上記した重みWijを算出するために行う処理(重み算出処理)の手順について説明する。図4〜図6は、重み算出処理の処理手順を示すフローチャートである。ここでは、数値流体シミュレーションを用いた温度分布解析によって勢力(R1ij,R2ij)を取得する方法を例にとって重みWijを算出する場合の処理手順を例示する。
First, a processing procedure in the case where the temperature field of the target fluid system can be regarded as almost steady will be described with reference to FIGS. First, the procedure of the process (weight calculation process) performed by the
図4および図5に示す重み算出処理では、先ず、図4に示すように、流れ場取得工程として、数値流体シミュレーションを用い、推定対象の流体系の代表的な流れ場計算条件を設定し(ステップS1)、設定した流れ場計算条件をもとに定常流れ場を計算して、実際の流れ場とする(ステップS3)。ここで、流体の流れ場の計算は、公知技術を用いて行う。具体的には、流体の流れ場と温度場とを求めることができる流体解析ソルバーであれば、市販品を含め何を用いてもよく、例えば、ANSYS FLUENT(登録商標)等を用いることで実際の流れ場を計算する。また、従来から、2次元の流れ場を算出する方法や、3次元の流れ場を算出する方法が知られているが、推定対象の流体系の特徴に応じて、2次元の流れ場または3次元の流れ場を算出する方法を適宜選択して用いることとしてもよい。 In the weight calculation process shown in FIGS. 4 and 5, first, as shown in FIG. 4, as a flow field acquisition process, a representative fluid field calculation condition of a fluid system to be estimated is set using a numerical fluid simulation ( Step S1), a steady flow field is calculated based on the set flow field calculation conditions to obtain an actual flow field (step S3). Here, the calculation of the fluid flow field is performed using a known technique. Specifically, any fluid analysis solver can be used as long as it is a fluid analysis solver capable of obtaining a fluid flow field and a temperature field. For example, ANSYS FULL (registered trademark) or the like can be used in practice. Calculate the flow field. Conventionally, a method for calculating a two-dimensional flow field and a method for calculating a three-dimensional flow field are known. Depending on the characteristics of the fluid system to be estimated, a two-dimensional flow field or 3 A method for calculating a three-dimensional flow field may be appropriately selected and used.
続いて、領域設定工程として、流体系内の温度実測部位i、発吸熱部位i、および流入出部位iの各部位iに対応する温度実測領域i、発吸熱領域i、および流入出領域iの各領域i(i=1〜K+L+M)を設定する(ステップS5)。また、推定点設定工程として、流体系内に温度推定点j(j=1〜N)を設定する(ステップS7)。続いて、温度実測部位i、発吸熱部位i、および流入出部位iの各部位iの中から、下流側勢力R1ijを取得する部位iを指定する(ステップS9)。ここでの処理は、ステップS9〜ステップS23の繰り返しの度にiの値を1〜K+L+Mの範囲で順次インクリメントしていくことで実現できる。 Subsequently, as a region setting step, the temperature measurement region i, the heat generation / absorption region i, and the inflow / outflow region i corresponding to each region i of the temperature measurement region i, the heat generation / absorption region i, and the inflow / outflow region i in the fluid system. Each region i (i = 1 to K + L + M) is set (step S5). Further, as an estimated point setting step, a temperature estimated point j (j = 1 to N) is set in the fluid system (step S7). Subsequently, the site i from which the downstream force R 1ij is acquired is designated from the site i of the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i (step S9). This process can be realized by sequentially incrementing the value of i in the range of 1 to K + L + M every time step S9 to step S23 are repeated.
続いて、温度推定点jを指定する(ステップS11)。ここでの処理は、ステップS11〜ステップS21の繰り返しの度にjの値を1〜Nの範囲で順次インクリメントしていくことで実現できる。その後、下流側勢力取得工程として、先ず、数値流体シミュレーションに必要な境界条件を与えるが、ここでは、温度分布解析によって勢力を算出するため、境界条件として、各領域iにおける温度の値を与える。具体的には、指定した部位iに対応する領域iの温度を「1」に固定する境界条件を与え、他の領域i´(i≠i´)の温度を「0」に固定する境界条件を与える(ステップS13)。その後、実際の流れ場を用い、移流拡散現象の数値流体シミュレーションを行って、与えた境界条件で温度分布解析を行う。具体的には、定常温度分布の計算を行い(ステップS15)、得られた定常温度分布に従って温度推定点jにおける温度値を取得する(ステップS17)。この温度値が、実際の流れ場での温度推定点jにおける領域iの流体成分の比率に相当する。そして、取得した温度値を、温度推定点jにおける指定した部位iの下流側勢力R1ijの値とする(ステップS19)。 Subsequently, the temperature estimation point j is designated (step S11). This process can be realized by sequentially incrementing the value of j in the range of 1 to N every time step S11 to step S21 are repeated. Thereafter, as the downstream side force acquisition step, first, boundary conditions necessary for the numerical fluid simulation are given. Here, in order to calculate the forces by temperature distribution analysis, the temperature value in each region i is given as the boundary condition. Specifically, a boundary condition for fixing the temperature of the region i corresponding to the designated part i to “1” is given, and the boundary condition for fixing the temperature of the other region i ′ (i ≠ i ′) to “0”. (Step S13). Then, using an actual flow field, a numerical fluid simulation of the advection-diffusion phenomenon is performed, and temperature distribution analysis is performed under the given boundary conditions. Specifically, the steady temperature distribution is calculated (step S15), and the temperature value at the temperature estimation point j is acquired according to the obtained steady temperature distribution (step S17). This temperature value corresponds to the ratio of the fluid component in the region i at the temperature estimation point j in the actual flow field. Then, the acquired temperature value is set as the value of the downstream force R 1ij of the designated site i at the temperature estimation point j (step S19).
続いて、温度推測値jのすべてについて下流側勢力R1ijを取得したか否かを判定する。下流側勢力R1ijが未取得の温度推測値jがある場合には(ステップS21:No)、ステップS11に戻って上記した処理を繰り返す。続いて、温度実測部位i、発吸熱部位i、および流入出部位iの全ての部位iについて下流側勢力R1ijを取得したか否かを判定する。下流側勢力R1ijが未取得の部位iがある場合には(ステップS23:No)、ステップS9に戻って上記した処理を繰り返す。全ての部位iについて下流側勢力R1ijを取得したならば(ステップS23:Yes)、続いて、図5に示すように、反転流れ場取得工程として、実際の流れ場の流速ベクトルを反転させた流れ場を反転流れ場として計算する(ステップS25)。そして、温度実測部位i、発吸熱部位i、および流入出部位iの各部位iの中から、上流側勢力R2ijを取得する部位iを指定する(ステップS27)。ステップS9と同様に、ステップS27〜ステップS45の繰り返しの度にiの値を1〜K+L+Mの範囲で順次インクリメントしていけばよい。 Subsequently, it is determined whether or not the downstream power R 1ij has been acquired for all of the estimated temperature values j. If there is an estimated temperature value j for which the downstream force R 1ij has not been acquired (step S21: No), the process returns to step S11 and the above-described processing is repeated. Subsequently, it is determined whether or not the downstream power R 1ij has been acquired for all the sites i of the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i. If there is a site i for which the downstream force R 1ij has not been acquired (step S23: No), the process returns to step S9 and the above-described processing is repeated. If the downstream forces R 1ij have been acquired for all the sites i (step S23: Yes), then, as shown in FIG. 5, the flow velocity vector of the actual flow field is inverted as an inverted flow field acquisition step. The flow field is calculated as an inverted flow field (step S25). And the site | part i which acquires upstream side force R2ij is specified from each site | part i of the temperature measurement site | part i, the heat-release / heat absorption site | part i, and the inflow / outflow site | part i (step S27). Similarly to step S9, the value of i may be sequentially incremented in the range of 1 to K + L + M each time step S27 to step S45 are repeated.
続いて、上流側勢力R2ijを取得する温度推定点jを指定する(ステップS29)。ステップS11と同様に、ステップS29〜ステップS43の繰り返しの度にjの値を1〜Nの範囲で順次インクリメントしていけばよい。続いて、指定した部位iが温度実測部位iまたは流入出部位iの場合と、発吸熱部位iの場合とで処理を分岐する。すなわち、指定した部位iが温度実測部位iまたは流入出部位iの場合には(ステップS31:Yes)、上流側勢力取得工程として、先ず、この指定した温度実測部位iまたは流入出部位iに対応する領域iの温度を「1」に固定する境界条件を与え、他の領域i´(i≠i´)の温度を「0」に固定する境界条件を与える(ステップS33)。その後、反転流れ場を用いて数値流体シミュレーションを行い、与えた境界条件で温度分布解析 を行う。具体的には、定常温度分布の計算を行い(ステップS35)、得られた定常温度分布に従って温度推定点jにおける温度値を取得する(ステップ S37)。この温度値が、反転流れ場での温度推定点jにおける領域iの流体成分の比率に相当する。そして、取得した温度値を、温度推定点jにおける指定した部位iの上流側勢力R2ijの値とし(ステップS39)、その後ステップS43に移行する。 Subsequently, the temperature estimation point j from which the upstream power R 2ij is acquired is designated (step S29). Similarly to step S11, the value of j may be sequentially incremented in the range of 1 to N each time step S29 to step S43 are repeated. Subsequently, the process branches depending on whether the designated site i is the temperature measurement site i or the inflow / outflow site i and the heat generation / absorption site i. That is, when the designated site i is the temperature measured site i or the inflow / outflow site i (step S31: Yes), first, as the upstream side force acquisition step, the designated temperature measured site i or the inflow / outflow site i corresponds. A boundary condition for fixing the temperature of the region i to be fixed to “1” is given, and a boundary condition for fixing the temperature of the other region i ′ (i ≠ i ′) to “0” is given (step S33). Then, numerical fluid simulation is performed using the reversal flow field, and temperature distribution analysis is performed under the given boundary conditions. Specifically, the steady temperature distribution is calculated (step S35), and the temperature value at the temperature estimation point j is acquired according to the obtained steady temperature distribution (step S37). This temperature value corresponds to the ratio of the fluid component in the region i at the temperature estimation point j in the reverse flow field. Then, the acquired temperature value is set as the value of the upstream force R 2ij of the designated part i at the temperature estimation point j (step S39), and then the process proceeds to step S43.
一方、指定した部位iが温度実測部位iまたは流入出部位iではなく、発吸熱部位iの場合には(ステップS31:No)、指定した部位iの上流側勢力R2ijの値を「0」とし(ステップS41)、その後ステップS43に移行する。そして、ステップS43では、全ての温度推定点jについて上流側勢力R2ijを取得したか否かを判定する。上流側勢力R2ijが未取得の温度推定点jがある場合には(ステップS43:No)、ステップS29に戻って上記した処理を繰り返す。全ての温度推定点jについて上流側勢力R2ijを取得したならば(ステップS43:Yes)、全ての部位iについて上流側勢力R2ijを取得したか否かを判定する。上流側勢力R2ijが未取得の部位iがある場合には(ステップS45:No)、ステップS27に戻って上記した処理を繰り返す。 On the other hand, when the designated part i is not the temperature actually measured part i or the inflow / outflow part i but the heat generating and absorbing part i (step S31: No), the value of the upstream side force R 2ij of the designated part i is set to “0”. (Step S41), and then the process proceeds to Step S43. In step S43, it is determined whether or not the upstream power R 2ij has been acquired for all temperature estimation points j. When there is a temperature estimation point j for which the upstream power R 2ij has not been acquired (step S43: No), the process returns to step S29 and the above-described processing is repeated. If the upstream force R 2ij has been acquired for all the temperature estimation points j (step S43: Yes), it is determined whether or not the upstream force R 2ij has been acquired for all the sites i. If there is a site i for which the upstream force R 2ij has not been acquired (step S45: No), the process returns to step S27 and the above-described processing is repeated.
上記下流側勢力R1ijおよび上流側勢力R2ijの取得法では、指定された温度推定点jに対して、下流側勢力R1ijおよび上流側勢力R2ijを取得するとしたが、本手法は容易に流体領域全体に拡張し、流体領域全体の下流側勢力分布および上流側勢力分布を取得することができる。具体的には、流体領域内を十分細かく覆うような配置、たとえば数値流体シミュレーションの全ての計算グリッドj’の位置にそれぞれ温度推定点jを配置し、全ての温度推定点jに対して下流側勢力R1ijおよび上流側勢力R2ijを取得すれば、流体領域全体の下流側勢力分布および上流側勢力分布を取得できる。 The acquisition method of the downstream force R 1ij and upstream force R 2ij, for a given temperature estimation point j, was to get the downstream force R 1ij and upstream force R 2ij, this technique is easily The entire fluid region can be expanded to obtain the downstream force distribution and the upstream force distribution of the entire fluid region. Specifically, the temperature estimation points j are arranged so as to cover the fluid region sufficiently finely, for example, at the positions of all the calculation grids j ′ of the numerical fluid simulation, and downstream of all the temperature estimation points j. If the force R 1ij and the upstream force R 2ij are acquired, the downstream force distribution and the upstream force distribution of the entire fluid region can be acquired.
全ての部位iについて上流側勢力R2ijを取得したならば(ステップS45:Yes)、続いて、図6に示すように、重み算出工程として、先ず、各部位iの中から、重みWijを算出する部位iを指定する(ステップS47)。ステップS9と同様に、ステップS47〜ステップS59の繰り返しの度にiの値を1〜K+L+Mの範囲で順次インクリメントしていけばよい。続いて、重みWijを算出する温度推定点jを指定する(ステップS49)。ステップS11と同様に、ステップS49〜ステップS57の繰り返しの度にjの値を1〜Nの範囲で順次インクリメントしていけばよい。 If the upstream forces R 2ij have been acquired for all the parts i (step S45: Yes), then, as shown in FIG. 6, as a weight calculation step, first, the weight W ij is selected from each part i. The part i to be calculated is designated (step S47). Similarly to step S9, the value of i may be sequentially incremented in the range of 1 to K + L + M each time step S47 to step S59 are repeated. Subsequently, the temperature estimation point j for calculating the weight W ij is designated (step S49). Similarly to step S11, the value of j may be sequentially incremented in the range of 1 to N every time step S49 to step S57 are repeated.
続いて、指定した部位iの温度が既知であるか否かを判定する。一般に指定した部位iが温度実測部位iの場合、温度は既知である。ただし、計測器の故障などで一時的に温度観測が不可能となる場合があるが、このような場合は実測部位iの温度を未知としても良い。一方、発吸熱部位iまたは流入出部位iについては、温度が未知の場合がある。このため、指定した部位iの温度が既知の場合には(ステップS51:Yes)、指定した部位iの勢力(R1ij,R2ij)をもとに、重み関数W(R1ij,R2ij)を用いて温度推定点jにおける指定した部位iの重みWijを算出する(ステップS53)。 Subsequently, it is determined whether or not the temperature of the designated part i is known. In general, when the designated site i is the temperature measurement site i, the temperature is known. However, there may be a case where the temperature cannot be observed temporarily due to a failure of the measuring instrument or the like. In such a case, the temperature of the actually measured part i may be unknown. On the other hand, the temperature may not be known for the heat generating and absorbing part i or the inflow / outflow part i. For this reason, when the temperature of the designated site | part i is known (step S51: Yes), weight function W ( R1ij , R2ij ) based on the influence ( R1ij , R2ij ) of the designated site | part i. Is used to calculate the weight W ij of the designated part i at the temperature estimation point j (step S53).
重み関数W(R1ij,R2ij)は、次式(4)に示すように、任意のR2ijに対してR1ijの単調非減少関数となり、任意のR1ijに対してR2ijの単調非減少関数となるような関数であって、かつ下流側勢力R1ijおよび上流側勢力R2ijがともに「0」の場合、すなわち、R1ij=R2ij=0の場合に「0」となる関数であれば、どのようなものでも適用できる。 Weighting function W (R 1ij, R 2ij), as shown in the following equation (4) becomes a monotonically non-decreasing function of R 1ij for any R 2ij, monotonous R 2ij to any R 1ij non A function that is a decreasing function, and a function that is “0” when the downstream force R 1ij and the upstream force R 2ij are both “0”, that is, when R 1ij = R 2ij = 0. Anything can be applied.
重み関数W(R1ij,R2ij)は、空間スケールや流速スケール、温度実測部位間の間隔等によって最適な関数形が変わってくるが、比較的簡単で、どのような流体系を推定対象とする場合であっても幅広く利用できる重み関数として、次式(5)に示すような下流側勢力R1ijおよび上流側勢力R2ijの部位i毎の平均値を算出する重み関数W(R1ij,R2ij)が挙げられる。例えば、この次式(5)に示す重み関数W(R1ij,R2ij)は、推定対象の流体系における発吸熱部位iの有無が把握できない場合、または、流体系が発吸熱部位iを含むものの、その正確な位置が把握できない場合等に適している。 The optimal function form of the weight function W (R 1ij , R 2ij ) varies depending on the space scale, the flow velocity scale, the interval between the temperature measurement sites, and the like. as a weight function which can be widely used even in the case where the weighting function to calculate the average value of each region i of the downstream force R 1ij and upstream force R 2ij as shown in the following equation (5) W (R 1ij, R 2ij ). For example, the weighting function W (R 1ij , R 2ij ) shown in the following equation (5) can be used when the existence or nonexistence of the endothermic part i in the fluid system to be estimated cannot be grasped, or the fluid system includes the endothermic part i. However, it is suitable when the exact position cannot be grasped.
全ての発吸熱部位iおよび流入出部位iの温度が既知である、または、全ての発吸熱部位iおよび流入出部位iの近傍に温度実測部位iが存在していて(全ての発吸熱部位iおよび流入出部位iの予め設定される所定の距離範囲内にそれぞれ温度実測部位iが存在していて)かつ全ての発吸熱部位iおよび流入出部位iが温度実測部位iから見て流れの上流側にある場合には、上流側勢力R2ijよりも下流側勢力R1ijの方が精度が高いため、重み関数W(R1ij,R2ij)として、下流側勢力R1ijのみを用いる次式(6)に示す重み関数W(R1ij,R2ij)を用いるとよい。 The temperatures of all the heat generation / absorption sites i and the inflow / outflow sites i are known, or the temperature measurement sites i exist in the vicinity of all the heat generation / absorption sites i and the inflow / outflow sites i (all the heat generation / absorption sites i). And the temperature measurement site i exists within a predetermined distance range set in advance of each of the inflow / outflow sites i) and all the heat generation / absorption sites i and the inflow / outflow sites i are upstream of the flow as viewed from the temperature measurement site i. Since the downstream side force R 1ij is more accurate than the upstream side force R 2ij, the following equation using only the downstream side force R 1ij as the weighting function W (R 1ij , R 2ij ) The weight function W (R 1ij , R 2ij ) shown in 6) may be used.
全ての発吸熱部位iおよび流入出部位iの温度が既知である、または、全ての発吸熱部位iおよび流入出部位iの近傍に温度実測部位iが存在していてかつ全ての発吸熱部位iおよび流入出部位iが温度実測部位iから見て流れの下流側にある場合には、下流側勢力R1ijよりも上流側勢力R2ijの方が精度が高いため、重み関数W(R1ij,R2ij)として、上流側勢力R2ijのみを用いる次式(7)に示す重み関数W(R1ij,R2ij)を用いるとよい。 The temperatures of all the heat generation / absorption sites i and inflow / outflow sites i are known, or the temperature measurement sites i exist in the vicinity of all the heat generation / absorption sites i and the inflow / outflow sites i, and all the heat generation / absorption sites i When the inflow / outflow site i is on the downstream side of the flow as viewed from the temperature measurement site i, the upstream side force R 2ij is more accurate than the downstream side force R 1ij, so the weight function W (R 1ij , As R 2ij ), it is preferable to use a weighting function W (R 1ij , R 2ij ) represented by the following equation (7) using only the upstream power R 2ij .
温度分布に大きく寄与する発吸熱部位iまたは流入出部位iの流体系内の位置がすべて特定できており、この発吸熱部位iおよび流入出部位iのうちの一部または全ての温度が未知である場合には、次式(8)に示す重み関数W(R1ij,R2ij)を用いるとよい。S1jは、温度が既知である温度実測部位i、発吸熱部位i、および流入出部位iについての下流側勢力R1ijの総和であり、S2jは、温度が既知である温度実測部位i、発吸熱部位i、および流入出部位iについての上流側勢力R2ijの総和である。 All the positions of the heat generation / absorption site i or the inflow / outflow region i in the fluid system that greatly contribute to the temperature distribution can be specified, and the temperature of a part or all of the heat generation / absorption site i and the inflow / outflow site i is unknown. In some cases, a weighting function W (R 1ij , R 2ij ) shown in the following equation (8) may be used. S 1j is the sum of the downstream measured force R 1ij for the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i where the temperature is known, and S 2j is the temperature measurement site i where the temperature is known, It is the sum total of the upstream side power R 2ij for the endothermic part i and the inflow / outlet part i.
温度分布に大きく寄与する発吸熱部位iまたは流入出部位iの流体系内の位置がすべて特定できており、この発吸熱部位iおよび流入出部位iのうちの一部または全ての温度が未知であり、かつ温度推定点jにおいて下流側勢力R1ijおよび上流側勢力R2ijともに小さな値となってしまうような場合は、次式(9)に示す重み関数W(R1ij,R2ij)を用いるとよい。S1jは、温度が既知である温度実測部位i、発吸熱部位i、および流入出部位iについての下流側勢力R1ijの総和であり、Savejは、温度が既知である温度実測部位i、発吸熱部位i、および流入出部位iについての下流側勢力と上流側勢力との平均値(1/2)×(R1ij+R2ij)の総和である。 All the positions of the heat generation / absorption site i or the inflow / outflow region i in the fluid system that greatly contribute to the temperature distribution can be specified, and the temperature of a part or all of the heat generation / absorption site i and the inflow / outflow site i is unknown. There, and if that becomes the downstream force R 1ij and upstream force R 2ij both small value in the temperature estimate point j uses a weighting function W in the following equation (9) (R 1ij, R 2ij) Good. S 1j is the sum of the temperature measurement site i where the temperature is known, the heat generation / absorption site i, and the downstream force R 1ij for the inflow / outflow site i, and S avej is the temperature measurement site i where the temperature is known, It is the sum total of the average value (1/2) × (R 1ij + R 2ij ) of the downstream side force and the upstream side force for the heat generating / absorbing site i and the inflow / outflow site i.
重み関数W(R1ij,R2ij)は、ステップS7で設定した全ての温度推定点jに対して同じものを一律に適用してもよいし、温度推定点j毎に、条件に合った適切な重み関数W(R1ij,R2ij)を選択的に用いることとしてもよい。 The same weighting function W (R 1ij , R 2ij ) may be applied uniformly to all the temperature estimation points j set in step S7, or an appropriate one that meets the conditions for each temperature estimation point j. The weight function W (R 1ij , R 2ij ) may be selectively used.
図6に戻り、以上のようにして重みWijを算出したならば、ステップS57に移行する。また、指定した部位iの温度が未知の場合には(ステップS51:No)、温度推定点jにおける指定した部位iの重みWijを「0」とし(ステップS55)、その後ステップS59に移行する。ステップS57では、全ての温度推定点jについて重みWijを算出したか否かを判定する。重みWijが未算出の温度推定点jがある場合には(ステップS57:No)、ステップS49に戻って上記した処理を繰り返す。そして、全ての温度推定点jについて重みWijを算出したならば(ステップS57:Yes)、ステップS59に移行する。 Returning to FIG. 6, when the weight W ij is calculated as described above, the process proceeds to step S57. If the temperature of the designated part i is unknown (step S51: No), the weight W ij of the designated part i at the temperature estimation point j is set to “0” (step S55), and then the process proceeds to step S59. . In step S57, it is determined whether or not the weights Wij have been calculated for all temperature estimation points j. When there is a temperature estimation point j for which the weight W ij has not been calculated (step S57: No), the process returns to step S49 and the above-described processing is repeated. If the weights W ij are calculated for all temperature estimation points j (step S57: Yes), the process proceeds to step S59.
ステップS59では、全ての部位iについて重みWijを算出したか否かを判定する。重みWijが未算出の部位iがある場合には(ステップS59:No)、ステップS47に戻って上記した処理を繰り返す。そして、全ての部位iについて重みWijを算出したならば(ステップS59:Yes)、算出した温度推定点jにおける部位i毎の重みWijを記憶装置に保存し(ステップS61)、重み算出処理を終える。流体系内のある温度推定点jの推定温度Tejは、この温度推定点jにおける部位i毎の重みWijを上記した手順で算出すれば推定できるが、流体系全体の温度分布を推定し、可視化するためには、流体系の全域に温度推定点jを設定し、設定した全ての温度推定点jについて重みWijを算出しておく必要がある。この場合には、事前に全ての温度推定点jについて部位i毎の重みWijを算出し、記憶装置にデータベース(重みデータベース)として保存しておくことが好ましい。 In step S59, it is determined whether or not the weights Wij have been calculated for all the parts i. If there is a part i for which the weight W ij has not been calculated (step S59: No), the process returns to step S47 and the above processing is repeated. If the weights W ij have been calculated for all the parts i (step S59: Yes), the weights W ij for each part i at the calculated temperature estimation point j are stored in the storage device (step S61), and weight calculation processing is performed. Finish. The estimated temperature Te j at a certain temperature estimation point j in the fluid system can be estimated by calculating the weight W ij for each part i at the temperature estimation point j by the above procedure, but the temperature distribution of the entire fluid system is estimated. In order to visualize the temperature, it is necessary to set the temperature estimation point j over the entire area of the fluid system and calculate the weights W ij for all the set temperature estimation points j. In this case, it is preferable to calculate weights W ij for each part i for all temperature estimation points j in advance and store them in a storage device as a database (weight database).
次に、以上のようにして算出した温度推定点jに対する部位i毎の重みWijを用いて任意の温度推定点jの温度を推定するための処理(温度推定処理)の手順について説明する。図7は、温度推定処理の処理手順を示すフローチャートである。 Next, a procedure of a process (temperature estimation process) for estimating the temperature of an arbitrary temperature estimation point j using the weight W ij for each part i with respect to the temperature estimation point j calculated as described above will be described. FIG. 7 is a flowchart showing a processing procedure of the temperature estimation processing.
温度推定処理では、図7に示すように、温度推定工程として、先ず、温度実測部位i、発吸熱部位i、および流入出部位iの既知温度Tiを取得する(ステップS71)。温度実測部位iについては、該当する温度実測部位iに設置された温度計測装置2から入力される温度実測値を既知温度Tiとして取得する。発吸熱部位iおよび流入出部位iについては、該当する部位iに温度計測装置が設置されており、温度を計測している場合や、該当する部位iの温度が例えば固定である等既知であり、記憶装置に予め保存されている場合等、何らかの手段で取得が可能な場合には、これを取得する。
In the temperature estimation process, as shown in FIG. 7, as the temperature estimation step, first, the temperature measurement part i, the heat generation / absorption part i, and the known temperature Ti of the inflow / outflow part i are acquired (step S71). For the temperature measurement part i, the temperature measurement value input from the
続いて、温度を推定する温度推定点jを指定する(ステップS73)。ここでの処理は、ステップS73〜ステップS81の繰り返しの度にjの値を1〜Nの範囲で順次インクリメントしていくことで実現できる。続いて、指定した温度推定点jに対する部位i毎の重みWijを記憶装置から読み出して取得する(ステップS75)。例えば、上記した重みデータベースから指定した温度推定点jについての重みWijを取得する。そして、上記した式(2)に従い、ステップS71で取得した各部位iの既知温度TiとステップS75で取得した重みWijとを用いた重み付き平均処理を行って、温度推定点jの推定温度Tejを算出する(ステップS77)。その後、算出した温度推定点jの推定温度Tejを記憶装置に保存する(ステップS79)。 Subsequently, a temperature estimation point j for estimating the temperature is designated (step S73). This process can be realized by sequentially incrementing the value of j in the range of 1 to N every time step S73 to step S81 are repeated. Subsequently, the weight W ij for each part i with respect to the designated temperature estimation point j is read from the storage device and acquired (step S75). For example, the weight W ij for the temperature estimation point j designated from the weight database is acquired. Then, according to the above equation (2), a weighted average process using the known temperature T i of each part i acquired in step S71 and the weight W ij acquired in step S75 is performed to estimate the temperature estimation point j. The temperature Te j is calculated (step S77). Thereafter, the estimated temperature Te j at the calculated temperature estimation point j is stored in the storage device (step S79).
その後、全ての温度推定点jについて推定温度Tejを算出したか否かを判定する。推定温度Tejが未算出の温度推定点jがある場合には(ステップS81:No)、ステップS73に戻って上記した処理を繰り返す。一方、全ての温度推定点jの推定温度Tejを算出したならば(ステップS81:Yes)、温度推定処理を終える。 Thereafter, it is determined whether or not the estimated temperature Te j has been calculated for all temperature estimated points j. If there is a temperature estimation point j for which the estimated temperature Te j has not been calculated (step S81: No), the process returns to step S73 and the above-described processing is repeated. On the other hand, if the estimated temperatures Te j of all the temperature estimation points j have been calculated (step S81: Yes), the temperature estimation process ends.
次に、対象の流体系の温度場が時間変動しうる場合の処理手順について同様に図8〜図12を参照して説明する。温度場が変動しうる場合は、前記勢力に加え、時系列の温度データの取得、伝達時間の算出、下流側重みと上流側重みとの算出、下流側既知温度と上流側既知温度との算出が必要になる。 Next, a processing procedure in the case where the temperature field of the target fluid system may fluctuate over time will be described with reference to FIGS. If the temperature field can fluctuate, in addition to the power, acquisition of time-series temperature data, calculation of transmission time, calculation of downstream weight and upstream weight, calculation of downstream known temperature and upstream known temperature Is required.
まず、伝達時間、すなわち下流側伝達時間τ1ijと上流側伝達時間τ2ijとの算出について説明する。図8は、伝達時間算出処理の処理手順を示すフローチャートである。下流側伝達時間τ1ijは移流拡散によって流体が温度実測部位i、発吸熱部位i、および流入出部位iから温度推定点jへ移動するのに要する時間を意味し、上流側伝達時間τ2ijは流体が温度推定点jから温度実測部位i、発吸熱部位i、および流入出部位iへ移動するのに要する時間を意味する。τ1ijは温度実測部位i、発吸熱部位i、および流入出部位iから見て流れの下流側の温度推定点jの方向へ流体が移動するのに要する時間なので、下流側伝達時間と呼び、同様に、τ2ijは温度実測部位i、発吸熱部位i、および流入出部位iから見て流れの上流側の温度推定点jの方向から流体が移動するのに要する時間なので、上流側伝達時間と呼ぶ。以下、下流側伝達時間と上流側伝達時間の対(τ1ij、τ2ij)を伝達時間と呼ぶ。 First, calculation of the transmission time, that is, the downstream transmission time τ 1ij and the upstream transmission time τ 2ij will be described. FIG. 8 is a flowchart showing a processing procedure of the transmission time calculation process. The downstream transmission time τ 1ij means the time required for the fluid to move from the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i to the temperature estimation point j by advection diffusion, and the upstream transmission time τ 2ij is It means the time required for the fluid to move from the temperature estimation point j to the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i. Since τ 1ij is the time required for the fluid to move in the direction of the temperature estimation point j on the downstream side of the flow as viewed from the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i, it is called the downstream transmission time, Similarly, since τ 2ij 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 viewed from the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i, the upstream transmission time Call it. Hereinafter, a pair (τ 1ij , τ 2ij ) of downstream transmission time and upstream transmission time is referred to as transmission time.
以下、伝達時間算出方法の一例として、温度の数値流体シミュレーションを用いた伝達時間(τ1ij、τ2ij)の算出方法を説明する。 Hereinafter, as an example of the transmission time calculation method, a calculation method of the transmission time (τ 1ij , τ 2ij ) using a temperature numerical fluid simulation will be described.
まず、数値流体シミュレーションを用い、流体系の代表的な境界条件を設定した後(ステップS101)、設定した境界条件に基づいて流れ場を算出する(ステップS103)。この流れ場は、温度場がほぼ定常とみなせる場合の手順(図4のステップS3)でもとめた実際の流れ場と同じものであるので、前記実際の流れ場をそのまま使ってもよい。 First, a representative 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 obtained in the procedure (step S3 in FIG. 4) when the temperature field can be regarded as almost steady, the actual flow field may be used as it is.
次に、流体系の温度実測部位i、発吸熱部位i、および流入出部位i(i=1〜K+L+M)を設定した後(ステップS105)、温度推定点j(j=1〜N) を設定する(ステップS107)。設定した温度実測部位i、発吸熱部位i、および流入出部位i(i=1〜K+L+M)と温度推定点j(j=1〜N)から、伝達時間(τ1ij、τ2ij)を算出する部位iと温度推定点jをそれぞれ指定する(ステップS109およびステップS111)。ここでの処理は、ステップS109〜ステップS131の繰り返しの度に、iの値を1からK+L+Mの範囲で順次インクリメントし、 ステップS109〜ステップS129の繰り返しの度に、jの値を1〜Nの範囲で順次インクリメントしていくことで実現できる。 Next, after setting the fluid temperature measurement site i, heat generation / absorption site i, and inflow / outflow site i (i = 1 to K + L + M) (step S105), the temperature estimation point j (j = 1 to N) is set. (Step S107). The transmission time (τ 1ij , τ 2ij ) is calculated from the set temperature measurement site i, heat generation / absorption site i, inflow / outflow site i (i = 1 to K + L + M) and temperature estimation point j (j = 1 to N). The site | part i and the temperature estimated point j are each designated (step S109 and step S111). In this processing, the value of i is sequentially incremented in the range of 1 to K + L + M every time the steps S109 to S131 are repeated, and the value of j is changed from 1 to N every time the steps S109 to S129 are repeated. This can be achieved by sequentially incrementing the range.
続いて、流体系全体に初期温度T0(単位K)を与えるとともに(ステップS113)、部位iの位置に発熱量S(単位W)の発熱条件を設定する(ステップS115)。この条件で温度分布の非定常計算を行い(ステップS117)、温度推定点jにおける温度上昇挙動を計算する。温度推定点jにおける温度が閾値温度TC(単位K)に到達したら、温度がT0からTCに到達するまでにかかった時間τ1ijを算出する(ステップS119)。τ1ijが下流側伝達時間となる。初期温度T0は、伝達時間に影響を与えない値なので、どのような値を与えても良い。発熱量Sおよび閾値温度TCに関しては、対象の流体系によって最適値が異なる。例えば、溶融亜鉛めっきポット、溶銑保持炉およびタンディッシュの一般的な場合、S=2,000KW、TC=T0+1K程度とすればよい。 Subsequently, the initial temperature T 0 (unit K) is given to the entire fluid system (step S113), and the heat generation condition of the heat generation amount S (unit W) is set at the position of the part i (step S115). Under this condition, the temperature distribution is unsteadyly calculated (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), a time τ 1ij taken until the temperature reaches T C from T 0 is calculated (step S119). τ 1ij is the downstream transmission time. Since the initial temperature T 0 is a value that does not affect the transmission time, any value may be given. For the heating value S and the threshold temperature T C, the optimum value differs by a fluid system of a subject. For example, in the general case of a hot dip galvanizing pot, a hot metal holding furnace, and a tundish, S = 2,000 KW and T C = T 0 + 1K may be set.
同様にして、流体系全体に初期温度T0を与えた後(ステップS121)、温度推定点jの位置に発熱量Sを与え(ステップS123)、温度分布の非定常計算を行い(ステップS125)、部位iの位置の温度がT0からTCに到達するまでにかかった時間τ2ijを算出する(ステップS127)。τ2ijが上流側伝達時間となる。伝達時間(τ1ij、τ2ij)は、移流拡散によって流体が温度実測部位i、発吸熱部位i、および流入出部位iから温度推定点jへ移動するのに要する時間および温度推定点jから温度実測部位i、発吸熱部位i、および流入出部位iへ流体が移動するのに要する時間と対応する指標であれば何でも良く、定義方法は特に限定されない。 Similarly, after the initial temperature T 0 is given to the entire fluid system (step S121), the heat generation amount S is given to the position of the temperature estimation point j (step S123), and the temperature distribution is unsteadyly calculated (step S125). Then, the time τ 2ij required for the temperature at the position i to reach T C from T 0 is calculated (step S127). τ 2ij is the upstream transmission time. The transmission time (τ 1ij , τ 2ij ) is the time required for the fluid to move from the temperature measurement site i, the heat-absorption / heat-absorption site i, and the inflow / outflow site i to the temperature estimation point j by the advection diffusion, Any method may be used as long as it is an index corresponding to the time required for the fluid to move to the actual measurement site i, the heat generation / absorption site i, and the inflow / outflow site i, and the definition method is not particularly limited.
ステップS129では、全ての温度推定点jについて伝達時間を算出したか否かを判定する。伝達時間(τ1ij、τ2ij)が未算出の温度推定点jがある場合には(ステップS129:No)、ステップS111に戻って上記した処理を繰り返す。そして、全ての温度推定点jについて伝達時間を算出したならば(ステップS129:Yes)、ステップS131に移行する。ステップS131では、全ての部位iについて伝達時間を算出したか否かを判定する。伝達時間(τ1ij、τ2ij)が未算出の部位iがある場合には(ステップS131:No)、ステップS109に戻って上記した処理を繰り返す。そして、全ての部位iについて伝達時間(τ1ij、τ2ij)を算出したならば(ステップS131:Yes)、算出した部位i毎の下流側伝達時間τ1ij、上流側伝達時間τ2ijを記憶装置に保存し(ステップS133)、伝達時間算出処理を終了する。流体系全体の温度分布を推定し、可視化するためには、事前に全ての温度推定点jについて部位i毎の下流側伝達時間τ1ijおよび上流側伝達時間τ2ijを算出し、記憶装置にデータベースとして保存しておくことが好ましい。 In step S129, it is determined whether or not the transmission time has been calculated for all temperature estimation points j. If there is a temperature estimation point j for which the transmission time (τ 1ij , τ 2ij ) has not been calculated (step S129: No), the process returns to step S111 and the above-described processing is repeated. If the transmission time is calculated for all temperature estimation points j (step S129: Yes), the process proceeds to step S131. In step S131, it is determined whether or not the transmission time has been calculated for all the parts i. If there is a part i for which the transmission time (τ 1ij , τ 2ij ) has not been calculated (step S131: No), the process returns to step S109 and the above-described processing is repeated. If the transmission times (τ 1ij , τ 2ij ) have been calculated for all the parts i (step S131: Yes), the downstream transmission time τ 1ij and the upstream transmission time τ 2ij for each part i are stored in the storage device. (Step S133), and the transmission time calculation process ends. In order to estimate and visualize the temperature distribution of the entire fluid system, the downstream transmission time τ 1ij and the upstream transmission time τ 2ij for each part i are calculated in advance for all temperature estimation points j, and the database is stored in the storage device. It is preferable to store as
続いて、図9および図10に示す、流れ場取得工程、領域設定工程、下流側勢力取得工程および上流側勢力取得工程を行う。図9および図10は、下流側重みW1ijおよび上流側重みW2ijの算出処理の処理手順を示すフローチャートである。これらは、対象の流体系の温度場がほぼ定常とみなせる場合について述べた前記手順(図4および図5のステップS1〜ステップS45)と同じ手順でよい。 Subsequently, a flow field acquisition process, a region setting process, a downstream power acquisition process, and an upstream power acquisition process shown in FIGS. 9 and 10 are performed. FIG. 9 and FIG. 10 are flowcharts showing a processing procedure for calculating the downstream weight W 1ij and the upstream weight W 2ij . These may be the same procedure as the above-described procedure (steps S1 to S45 in FIGS. 4 and 5) described in the case where the temperature field of the target fluid system can be regarded as almost steady.
次に、温度推定装置1が重みを算出するために行う処理(重み算出処理)の手順について説明する。温度場が時間変動しうる場合は、重みとして、下流側重みW1ijおよび/または上流側重みW2ijを算出する。全ての部位iについて下流側勢力R1ijおよび上流側勢力R2ijを取得したならば(ステップS201〜ステップS245、図8および図9参照)、続いて、図11に示すように、各部位iの中から、下流側重みW1ijおよび上流側重みW2ijを算出する部位iを指定する(ステップS247)。ここで、繰り返しの度にiの値を1〜K+L+Mの範囲で順次インクリメントしていけばよい。
Next, a procedure of processing (weight calculation processing) performed by the
続いて、下流側重みW1ijおよび上流側重みW2ijを算出する温度推定点jを指定する(ステップS249)。繰り返しの度にjの値を1〜Nの範囲で順次インクリメントしていけばよい。続いて、指定した部位iの温度が既知であるか否かを判定する(ステップS251)。一般に指定した部位iが温度実測部位iの場合、温度は既知である。ただし、計測器の故障などで一時的に温度観測が不可能となる場合、実測部位iの温度を未知としても良い。一方、発吸熱部位iまたは流 入出部位iについては、温度が未知の場合がある。このため、指定した部位iの温度が既知の場合には(ステップS251:Yes)、指定した部位iの勢力(R1ij,R2ij)をもとに、下流側重み関数W1(R1ij,R2ij)および上流側重み関数W2(R1ij,R2ij)を用いて温度推定点jにおける指定した部位iの下流側重みW1ijおよび上流側重みW2ijを算出する(ステップS253およびステップS255)。
Subsequently, the temperature estimation point j for calculating the downstream weight W 1ij and the upstream weight W 2ij is designated (step S249). What is necessary is just to increment the value of j sequentially in the range of 1-N at each repetition. Subsequently, it is determined whether or not the temperature of the designated part i is known (step S251). In general, when the designated site i is the temperature measurement site i, the temperature is known. However, when the temperature cannot be observed temporarily due to a failure of the measuring instrument or the like, the temperature of the actual measurement site i may be unknown. On the other hand, the temperature may not be known for the heat generation / absorption site i or the inflow / outflow site i. Therefore, when the temperature of the specified site i is known (step S251: Yes), power (R 1ij, R 2ij) for the specified site i on the basis of, downstream weighting function W 1 (R 1ij, R 2ij) and the
下流側重み関数W1(R1ij,R2ij)および上流側重み関数W2(R1ij,R2ij)は、次式(10)に示すように、任意のR2ijに対してR1ijの単調非減少関数となり、任意のR1ijに対してR2ijの単調非減少関数となるような関数であって、かつ下流側勢力R1ijおよび上流側勢力R2ijがともに「0」の場合、すなわち、R1ij=R2ij=0の場合に「0」となる関数であれば、どのようなものでも適用できる。 Downstream weighting function W 1 (R 1ij, R 2ij ) and the upstream weighting function W 2 (R 1ij, R 2ij ) , as shown in the following equation (10), monotonic R 1ij for any R 2ij becomes non-decreasing function, if a function such that monotonically non-decreasing function of R 2ij for any R 1ij, and downstream force R 1ij and upstream force R 2ij are both "0", i.e., Any function that is “0” when R 1ij = R 2ij = 0 can be applied.
下流側重み関数W1(R1ij,R2ij)および上流側重み関数W2(R1ij,R2ij)は、空間スケールや流速スケール、温度実測部位間の間隔等によって最適な関数形が変わってくる。比較的簡単で、どのような流体系を推定対象とする場合であっても幅広く利用できる下流側重み関数W1(R1ij,R2ij)および上流側重み関数W2(R1ij,R2ij)として、次式(11a)、(11b)に示すように、下流側勢力の0.5倍を下流側重みとし、上流側勢力の0.5倍を上流側重みとする下流側重み関数W1(R1ij,R2ij)および上流側重み関数W2(R1ij,R2ij)が挙げられる。例えば、この次式(11)は、推定対象の流体系における発吸熱部位iの有無が把握できない場合、または、流体系が発吸熱部位iを含むものの、その正確な位置が把握できない場合等に適している。 The downstream weight function W 1 (R 1ij , R 2ij ) and the upstream weight function W 2 (R 1ij , R 2ij ) have different optimal function forms depending on the spatial scale, the flow velocity scale, the interval between the temperature measurement parts, and the like. come. The downstream weight function W 1 (R 1ij , R 2ij ) and the upstream weight function W 2 (R 1ij , R 2ij ) that are relatively simple and can be used widely regardless of what fluid system is to be estimated. As shown in the following equations (11a) and (11b), a downstream weight function W 1 having a downstream weight of 0.5 times the downstream force and an upstream weight of 0.5 times the upstream force. (R 1ij , R 2ij ) and the upstream weight function W 2 (R 1ij , R 2ij ). For example, the following equation (11) is used when the presence / absence of the heat generation / absorption site i in the fluid system to be estimated cannot be grasped or when the fluid system includes the heat generation / absorption site i but the exact position thereof cannot be grasped. Is suitable.
全ての発吸熱部位iおよび流入出部位iの温度が既知である、または、全ての発吸熱部位iおよび流入出部位iの近傍に温度実測部位iが存在していて(全ての発吸熱部位iおよび流入出部位iの予め設定される所定の距離範囲内にそれぞれ温度実測部位iが存在していて)かつ全ての発吸熱部位iおよび流入出部位iが温度実測部位iから見て流れの上流側にある場合には、次式(12a)、(12b)に示すように下流側重み関数W1(R1ij,R2ij)は下流側勢力R1ijとし、上流側重み関数W2(R1ij,R2ij)は0とすると良い。 The temperatures of all the heat generation / absorption sites i and the inflow / outflow sites i are known, or the temperature measurement sites i exist in the vicinity of all the heat generation / absorption sites i and the inflow / outflow sites i (all the heat generation / absorption sites i). And the temperature measurement site i exists within a predetermined distance range set in advance of each of the inflow / outflow sites i) and all the heat generation / absorption sites i and the inflow / outflow sites i are upstream of the flow as viewed from the temperature measurement site i. If it is on the side, as shown in the following equations (12a) and (12b), the downstream weight function W 1 (R 1ij , R 2ij ) is the downstream force R 1ij and the upstream weight function W 2 (R 1ij). , R 2ij ) is preferably zero.
全ての発吸熱部位iおよび流入出部位iの温度が既知である、または、全ての発吸熱部位iおよび流入出部位iの近傍に温度実測部位iが存在していてかつ全ての発吸熱部位iおよび流入出部位iが温度実測部位iから見て流れの下流側にある場合には、次式(13a)、(13b)に示すように下流側重み関数W1(R1ij,R2ij)は0とし、上流側重み関数W2(R1ij,R2ij)は上流側勢力R2ijとすると良い。 The temperatures of all the heat generation / absorption sites i and inflow / outflow sites i are known, or the temperature measurement sites i exist in the vicinity of all the heat generation / absorption sites i and the inflow / outflow sites i, and all the heat generation / absorption sites i When the inflow / outflow part i is on the downstream side of the flow as viewed from the temperature measurement part i, the downstream weight function W 1 (R 1ij , R 2ij ) is expressed by the following equations (13a) and (13b): It is preferable that the upstream side weight function W 2 (R 1ij , R 2ij ) is the upstream side force R 2ij .
温度分布に大きく寄与する発吸熱部位iまたは流入出部位iの流体系内の位置がすべて特定できており、この発吸熱部位iおよび流入出部位iのうちの一部または全ての温度が未知である場合には、次式(14a)、(14b)に示すとおり、下流側重み関数W1(R1ij,R2ij)は下流側勢力R1ijとし、上流側重み関数W2(R1ij,R2ij)は式(14b)を用いるとよい。S1jは、温度が既知である温度実測部位i、発吸熱部位i、および流入出部位iについての下流側勢力R1ijの総和であり、Savejは、温度が既知である温度実測部位i、発吸熱部位i、および流入出部位iについての下流側勢力と上流側勢力の平均値(1/2)×(R1ij+R2ij)の総和である。 All the positions of the heat generation / absorption site i or the inflow / outflow region i in the fluid system that greatly contribute to the temperature distribution can be specified, and the temperature of a part or all of the heat generation / absorption site i and the inflow / outflow site i is unknown. In some cases, as shown in the following equations (14a) and (14b), the downstream weight function W 1 (R 1ij , R 2ij ) is the downstream force R 1ij and the upstream weight function W 2 (R 1ij , R 2ij ) may use the equation (14b). S 1j is the sum of the temperature measurement site i where the temperature is known, the heat generation / absorption site i, and the downstream force R 1ij for the inflow / outflow site i, and S avej is the temperature measurement site i where the temperature is known, This is the sum of the average value (1/2) × (R 1ij + R 2ij ) of the downstream side force and the upstream side force for the heat generating and absorbing part i and the inflow / outlet part i.
下流側重み関数W1(R1ij,R2ij)および上流側重み関数W2(R1ij,R2ij)は、設定した全ての温度推定点jに対して同じものを一律に適用してもよいし、温度推定点j毎に、条件に合った適切な下流側重み関数、上流側重み関数を選択的に用いることとしてもよい。 The downstream weight function W 1 (R 1ij , R 2ij ) and the upstream weight function W 2 (R 1ij , R 2ij ) may be uniformly applied to all set temperature estimation points j. Then, for each temperature estimation point j, an appropriate downstream weight function and upstream weight function that meet the conditions may be selectively used.
指定した部位iの温度が未知の場合には(ステップS251,No)、温度推定点jにおける指定した部位iの下流側重みW1ijおよび上流側重みW2ijを「0」とする(ステップS257)。そして、全ての温度推定点jについて下流側重みW1ijおよび上流側重みW2ijを算出したか否かを判定する(ステップS259)。下流側重みW1ijおよび上流側重みW2ijが未算出の温度推定点jがある場合には(ステップS259:No)、ステップS249に戻って上記した処理を繰り返す。そして、全ての温度推定点jについて重みを算出したならば(ステップS259:Yes)、ステップS261に移行する。 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 at the temperature estimation point j are set to “0” (step S257). . Then, it is determined whether or not the downstream weight W 1ij and the upstream weight W 2ij have been calculated for all temperature estimation points j (step S259). If there is a temperature estimation point j for which the downstream weight W 1ij and the upstream weight W 2ij have not been calculated (step S259: No), the process returns to step S249 and the above-described processing is repeated. If the weights are calculated for all temperature estimation points j (step S259: Yes), the process proceeds to step S261.
ステップS261では、全ての部位iについて下流側重みW1ijおよび上流側重みW2ijを算出したか否かを判定する。下流側重みW1ijおよび上流側重みW2ijが未算出の部位iがある場合には(ステップS261:No)、ステップS247に戻って上記した処理を繰り返す。そして、全ての部位iについて下流側重みW1ijおよび上流側重みW2ijを算出したならば(ステップS261:Yes)、算出した温度推定点jにおける部位i毎の下流側重みW1ijおよび上流側重みW2ijを記憶装置に保存し(ステップS263)、重み算出処理を終える。 In step S261, it is determined whether the downstream weight W 1ij and the upstream weight W 2ij have been calculated for all the parts i. If there is a part i for which the 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. If the downstream weight W 1ij and the upstream weight W 2ij are calculated for all the parts i (step S261: Yes), the downstream weight W 1ij and the upstream weight for each part i at the calculated temperature estimation point j. W 2ij is stored in the storage device (step S263), and the weight calculation process is terminated.
流体系全体の温度分布を推定し、可視化するためには、流体系の全域に温度推定点jを設定し、設定した全ての温度推定点jについて下流側重みW1ijおよび上流側重みW2ijを算出しておく必要がある。この場合には、事前に全ての温度推定点jについて部位i毎の下流側重みW1ijおよび上流側重みW2ijを算出し、記憶装置にデータベースとして保存しておくことが好ましい。 In order to estimate and visualize the temperature distribution of the entire fluid system, temperature estimation points j are set over the entire area of the fluid system, and the downstream weights W 1ij and the upstream weights W 2ij are set for all the estimated temperature points j. It is necessary to calculate. In this case, it is preferable that downstream weights W 1ij and upstream weights W 2ij for each part i are calculated in advance for all temperature estimation points j and stored in the storage device as a database.
次に、以上のようにして算出した温度推定点jに対する部位i毎の下流側伝達時間τ1ij、上流側伝達時間τ1ijおよび下流側重みW1ij、上流側重みW2ijおよび部位iごとの時系列温度データを用いて、任意の温度推定点jの温度を推定するための処理(温度推定処理)の手順について説明する。図12は、温度推定処理の処理手順を示すフローチャートである。 Next, the downstream transmission time τ 1ij , the upstream transmission time τ 1ij, the downstream weight W 1ij , the upstream weight W 2ij, and the region i for each part i with respect to the temperature estimation point j calculated as described above. A procedure of a process (temperature estimation process) for estimating the temperature of an arbitrary temperature estimation point j using the series temperature data will be described. FIG. 12 is a flowchart showing a processing procedure of the temperature estimation processing.
温度推定処理では、図12に示すように、まず温度推定を行う時点t0を決定する(ステップS301)。続いて、温度を推定する温度推定点jを指定する(ステップS303)。ここでの処理は、繰り返しの度にjの値を1〜Nの範囲で順次インクリメントしていくことで実現できる。続いて、指定した温度推定点jに対する部位i毎の下流側重みW1ij、上流側重みW2ij、下流側伝達時間τ1ijおよび上流側伝達時間τ2ijを記憶装置から読み出して取得する(ステップS305)。例えば、上記したデータベースから指定した温度推定点jについての下流側重みW1ij、上流側重みW2ij、下流側伝達時間τ1ij、および上流側伝達時間τ2ijを取得する。 A temperature estimation process, as shown in FIG. 12, first, determines the time t 0 of the temperature estimation (step S301). Subsequently, a temperature estimation point j for estimating the temperature is designated (step S303). This processing can be realized by sequentially incrementing the value of j in the range of 1 to N every time it is repeated. Subsequently, the downstream weight W 1ij , the upstream weight W 2ij , the downstream transmission time τ 1ij, and the upstream transmission time τ 2ij for each part i with respect to the designated temperature estimation point j are read from the storage device and acquired (step S305). ). For example, the downstream weight W 1ij , the upstream weight W 2ij , the downstream transmission time τ 1ij , and the upstream transmission time τ 2ij for the temperature estimation point j designated from the above database are acquired.
上記した式(3)に従い、取得した各部位iの時系列温度データTi(t)と下流側伝達時間τ1ijおよび上流側伝達時間τ2ijとから下流側既知温度Ti(t0−τ1ij)および上流側既知温度Ti(t0+τ2ij)を求めるとともに(ステップS307)、下流側重みW1ijおよび上流側重みW2ijとを算出する(ステップS309)。その後、算出した下流側重みW1ijおよび上流側重みW2ijを用いた重み付き平均処理を行って、時間t0における温度推定点jの推定温度Tej(t0)を算出する(ステップS311)。その後、算出した温度推定点jの推定温度Tej(t0)を記憶装置に保存する(ステップS313)。 According to the above equation (3), the downstream known temperature T i (t 0 −τ ) from the acquired time-series temperature data T i (t), the downstream transmission time τ 1ij and the upstream transmission time τ 2ij. 1ij ) and upstream known temperature T i (t 0 + τ 2ij ) are obtained (step S307), and downstream weight W 1ij and upstream weight W 2ij are calculated (step S309). Thereafter, a weighted average process using the calculated downstream weight W 1ij and upstream weight W 2ij is performed to calculate an estimated temperature Te j (t 0 ) of the temperature estimation point j at time t 0 (step S311). . Thereafter, the estimated temperature Te j (t 0 ) of the calculated temperature estimation point j is stored in the storage device (step S313).
その後、全ての温度推定点jについて推定温度Tej(t0)を算出したか否かを判定する(ステップS315)。推定温度Tej(t0)に未算出の温度推定点jがある場合には(ステップS315:No)、ステップS303へ戻って上記した処理を繰り返す。一方、全ての温度推定点jの推定温度Tej(t0)を算出したならば(ステップS315:Yes)、温度推定処理を終える。 Thereafter, it is determined whether or not the estimated temperature Te j (t 0 ) has been calculated for all temperature estimation points j (step S315). If there is an uncalculated temperature estimation point j in the estimated temperature Te j (t 0 ) (step S315: No), the process returns to step S303 and the above-described processing is repeated. On the other hand, if the estimated temperatures Te j (t 0 ) of all the temperature estimation points j have been calculated (step S315: Yes), the temperature estimation process ends.
以上説明したように、本実施の形態では、流体系内の任意の温度推定点jにおける全流体に対する領域i(温度実測領域i、発吸熱領域i、および流入出領域i)毎の流体成分の比率を勢力(R1ij,R2ij)として取得し、取得した勢力(R1ij,R2ij)をもとに、温度推定点jに対する部位i毎の重みWijまたは下流側重みW1ijおよび上流側重みW2ij、伝達時間(τ1ij,τ2ij)を算出することとした。そして、この部位i毎の重みWijを該当する部位iにおける既知温度Tiに重み付けして平均(重み付き平均処理)する、もしくは、時系列温度データを用いて、温度推定を行う時点t0に対して、下流側重みW1ijを該当する部位iにおける下流側既知温度Ti(t0−τ1ij)に、上流側重みW2ijを該当する部位iにおける上流側既知温度Ti(t0+τ2ij)に対応させて重み付けして平均することで、温度推定点jの推定温度を算出することとした。したがって、流体系の流れ場のもとでの熱の移流拡散を考慮して、温度推定点jの温度を高精度に推定することできる。これによれば、実際に温度計測装置2を設置する等して温度を実測するのが困難な場所であっても、温度を高精度に把握することが可能となる。したがって、温度計測装置2の配置に制約を与えることなく流体の流れによる熱輸送を考慮した高精度な温度推定を実現することができる。
As described above, in the present embodiment, the fluid component of each region i (temperature measurement region i, heat generation / absorption region i, and inflow / outflow region i) for all fluids at an arbitrary temperature estimation point j in the fluid system. the ratio force (R 1ij, R 2ij) obtained as the acquired power (R 1ij, R 2ij) based on the weight W ij or downstream weights W 1ij and upstream of each region i for temperature estimation point j The weight W 2ij and the transmission time (τ 1ij , τ 2ij ) are calculated. The time t 0 to an average by weighting the known temperature T i at the site i to the appropriate weight W ij for each this site i (weighted average processing), or, when using a series temperature data, the temperature estimation In contrast, the downstream side weight W 1ij is set to the downstream known temperature T i (t 0 −τ 1ij ) at the corresponding part i, and the upstream side weight W 2ij is set to the upstream side known temperature T i (t 0). + Τ 2ij ) is weighted and averaged to calculate the estimated temperature of the temperature estimation point j. Therefore, the temperature at the temperature estimation point j can be estimated with high accuracy in consideration of the advection diffusion of heat under the flow field of the fluid system. According to this, even in a place where it is difficult to actually measure the temperature by actually installing the
流体系内の全域に温度推定点jを設定し、各温度推定点jに対する部位i毎の重みWijを算出して例えば重みデータベースとして記憶装置に保存しておけば、この重みWijを読み出して取得するとともに、各部位iの既知温度Tiを取得して上記した式(2)に代入するだけで、瞬時に各部位iの既知温度Tiを補間した流体系の温度分布を推定することが可能となる。また、温度の時間変動が起こりうる流体系に対しても同様に、流体系内の全域に温度推定点jを設定し、各温度推定点jに対する部位i毎の下流側重みW1ij、上流側重みW2ij、および伝達時間(τ1ij,τ2ij)を算出して例えばデータベースとして記憶装置に保存しておき、温度計測装置2にて各部位iの温度と観測された時間とを逐次計測・保存し、任意の時点tにおける温度実測値を時系列温度データTi(t)として読み出せるようにしておけば、この下流側重みW1ij、上流側重みW2ij、および伝達時間(τ1ij,τ2ij)を読み出して取得するとともに、温度を推定したい時点t0と時系列温度データTi(t)とを用いて式(3)に代入するだけで、瞬時に時点t0における流体系の温度分布を推定することが可能となる。これによれば、本実施の形態の温度推定装置1は、リアルタイムの計算が不可欠な産業プロセスのオンラインモニタリングにも十分に活用でき、操業管理や制御機構に利用することが可能となる。
If temperature estimation points j are set in the entire region of the fluid system, weights W ij for each part i for each temperature estimation point j are calculated and stored in a storage device as a weight database, for example, and the weights W ij are read out. acquires Te, only into equation (2) described above to obtain the known temperature T i of each portion i, estimating the temperature distribution of the fluid system by interpolating the known temperature T i of each part i instantly It becomes possible. Similarly, with respect to a fluid system in which temperature fluctuations may occur, temperature estimation points j are set throughout the fluid system, and downstream weights W 1ij for each part i with respect to each temperature estimation point j, upstream side The weight W 2ij and the transmission time (τ 1ij , τ 2ij ) are calculated and stored in, for example, a storage device as a database, and the
上流側勢力R2ijを取得し、この上流側勢力R2ijや上流側伝達時間τ2ijを用いて重みWijもしくは下流側重みW1ijと上流側重みW2ijとを算出することとしたので、温度実測部位iから見て流れの上流側となる位置の温度についても、温度実測部位iの既知温度Ti(温度実測値)や時系列温度データTi(t)をもとに推定することができる。さらに、発吸熱部位iや流入出部位iの温度が既知の場合には、これら発吸熱部位iや流入出部位iの既知温度をさらに用いて温度推定点jの温度を推定することができる。これによれば、必ずしも流体の流れの最上流位置に温度計測装置2を配置する必要がない。したがって、温度計測 装置2の配置に制約を与えることなく流体系内の任意の位置の温度を推定することができる。
Since the upstream force R 2ij is acquired and the upstream force R 2ij and the upstream transmission time τ 2ij are used to calculate the weight W ij or the downstream weight W 1ij and the upstream weight W 2ij , the temperature The temperature at the position upstream of the flow when viewed from the measured site i can also be estimated based on the known temperature T i (temperature measured value) of the temperature measured site i and the time-series temperature data T i (t). it can. Further, when the temperature of the heat generating / absorbing site i or the inflow / outflow site i is known, the temperature of the temperature estimation point j can be estimated by further using the known temperature of the heat generating / absorbing site i or the inflow / outflow site i. According to this, it is not always necessary to arrange the
上記した実施の形態では、流体系が温度実測部位i、発吸熱部位i、および流入出部位iを含むこととして説明したが、発吸熱部位iおよび/または流入出部位を含まない場合には、これらを除いた部位iの勢力(R1ij,R2ij)および伝達時間(τ1ij,τ2ij)を取得し、重みWijもしくは下流側重みW1ijおよび上流側重みW2ijと、下流側既知温度Ti(t0−τ1ij)および上流側既知温度Ti(t0+τ2ij)とを算出すればよい。例えば、流体系が発吸熱部位iを含まない場合であれば、温度推定点jにおける流体の温度実測部位iおよび流入出部位iの勢力(R1ij,R2ij)、伝達時間(τ1ij,τ2ij)を取得し、取得した勢力(R1ij,R2ij)および伝達時間(τ1ij,τ2ij)をもとに、温度推定点jに対する温度実測部位iおよび流入出部位i毎の重みWijもしくは下流側重みW1ijおよび上流側重みW2ijと、下流側既知温度Ti(t0−τ1ij)および上流側既知温度Ti(t0+τ2ij)とを算出すればよい。同様に、流体系が流入出部位iを含まない場合には、温度推定点jにおける流体の温度実測部位iおよび発吸熱部位iの勢力(R1ij,R2ij)および伝達時間(τ1ij,τ2ij)を取得して温度推定点jに対する温度実測部位iおよび発吸熱部位i毎の重みWijもしくは下流側重みW1ijおよび上流側重みW2ijと、下流側既知温度Ti(t0−τ1ij)および上流側既知温度Ti(t0+τ2ij)とを算出し、流体系が発吸熱部位iおよび流入出部位iを含まない場合には、温度推定点jにおける流体の温度実測部位iの勢力(R1ij,R2ij)および伝達時間(τ1ij,τ2ij)を取得して温度推定点jに対する温度実測部位i毎の重みWijもしくは下流側重みW1ijおよび上流側重みW2ijと、下流側既知温度Ti(t0−τ1ij)および上流側既知温度Ti(t0+τ2ij)とを算出すればよい。これによれば、少なくとも流体系内の設置スペースが確保できる任意の位置に温度計や熱電対等の温度計測装置2を設置することで、この温度計測装置2が計測する温度実測値をもとに流体系内の任意の位置における温度を精度よく推定することが可能となる。
In the embodiment described above, the fluid system has been described as including the temperature measurement site i, the heat generation / absorption site i, and the inflow / outflow site i, but when the fluid system does not include the heat generation / absorption site i and / or the inflow / outflow site, The power (R 1ij , R 2ij ) and transmission time (τ 1ij , τ 2ij ) of the part i excluding these are acquired, the weight W ij or the downstream weight W 1ij and the upstream weight W 2ij, and the downstream known temperature Ti (t 0 −τ 1ij ) and upstream side known temperature Ti (t 0 + τ 2ij ) may be calculated. For example, if the fluid system does not include the heat generation / absorption part i, the fluid temperature measurement part i and the force (R 1ij , R 2ij ) of the fluid temperature measurement part i at the temperature estimation point j, the transmission time (τ 1ij , τ 2ij ), and based on the acquired power (R 1ij , R 2ij ) and transmission time (τ 1ij , τ 2ij ), the weight W ij for each temperature measurement portion i and inflow / outflow portion i for the temperature estimation point j Alternatively, the downstream weight W 1ij and the upstream weight W 2ij , the downstream known temperature Ti (t 0 −τ 1ij ), and the upstream known temperature Ti (t 0 + τ 2ij ) may be calculated. Similarly, when the fluid system does not include the inflow / outflow part i, the force (R 1ij , R 2ij ) and the transmission time (τ 1ij , τ) of the fluid temperature measurement part i and the heat generation / heat absorption part i at the temperature estimation point j and the weight W ij or downstream weights W 1ij and upstream weights W 2ij of each temperature measured sites i and outgoing endothermic sites i for temperature estimation point j to obtain the 2ij), downstream known temperature Ti (t 0 -τ 1ij ) And the upstream known temperature Ti (t 0 + τ 2ij ), and when the fluid system does not include the heat generating / absorbing part i and the inflow / outflow part i, the force of the fluid temperature measurement part i at the temperature estimation point j (R 1ij, R 2ij) and transmission time (τ 1ij, τ 2ij) weights W ij or downstream side weight W for each temperature measured site i for acquiring and estimating the temperature point j 1ij and the upstream heavy And W 2ij, downstream known temperature Ti (t 0 -τ 1ij) and the upstream known temperature Ti (t 0 + τ 2ij) and may be calculated. According to this, by installing the
本実施の形態では、勢力(R1ij,R2ij)をもとに重みWijもしくは下流側重みW1ijおよび上流側重みW2ijを算出し、算出した重みWijもしくは下流側重みW1ijおよび上流側重みW2ijを記憶装置に保存することとしたが、勢力(R1ij,R2ij)を保存しておき、重みWijもしくは下流側重みW1ijおよび上流側重みW2ijについては、温度推定の都度算出するようにしてもよい。 In the present embodiment, the weight W ij or downstream weight W 1ij and upstream weight W 2ij are calculated based on the forces (R 1ij , R 2ij ), and the calculated weight W ij or downstream weight W 1ij and upstream are calculated. Although the side weight W 2ij is stored in the storage device, the power (R 1ij , R 2ij ) is stored, and the weight W ij or the downstream weight W 1ij and the upstream weight W 2ij It may be calculated each time.
(実施の形態1)
次に、実施の形態1として、室内を適用対象とし、この室内を流れる流体系の温度推定および温度分布の可視化について説明する。図13は、実施の形態1において適用対象とする部屋3の内部を上方から示した模式図である。
(Embodiment 1)
Next, as
図13に示す部屋3は、例えば、一辺が1(m)の平面視略正方形状を有する。この部屋3は、図13に向かって左側の側壁311および右側の側壁312の対角においてそれぞれ幅が約0.2(m)の通路321,322を備え、この通路321,322の終端にそれぞれ窓331,332が取り付けられている。部屋3の図13に向かって上側の側壁313は、不図示の発熱源を備え、熱を発する発熱壁となっている。部屋3内の図13中に「×」を付して示す4箇所A〜Dにおいて、温度を推定するための温度計測装置としての温度計34−1〜34−4が設置されている。
The
本適用対象では、推定対象の流体系が、部屋3内を流れる空気、具体的には、図13中に矢印A3で示すように窓331から通路321を経て室内に流入し、通路322を経て窓332から室外へと流出する空気である。本適用対象では、温度計34−1〜34−4の設置位置A〜Dが温度実測部位iであり、例えば温度実測部位iである設置位置A〜Dを中心とした半径0.05(m)の円領域E31〜E34をそれぞれ温度実測領域iとする。発熱源を備えた側壁313が発吸熱部位i(発熱部位)であり、例えばこの側壁313の壁面領域Eを発吸熱領域iとする。窓331,332の領域F,G、すなわち、通路321,322の終端面がそれぞれ流入出部位i(窓331の領域Fが流入部位、窓332の領域Gが流出部位)であり、例えばこの領域F,Gを流入出領域iとする。側壁313が備える発熱源は、例えばその温度が50(℃)に制御され、窓331からは、10(℃)の空気が流入する。ただし、温度推定を行う際には、発吸熱部位iおよび流入出部位iの温度は未知とする。以下では、温度計34−1〜34−4の設置位置A〜Dに相当する温度実測部位iを適宜温度実測部位A〜Dと表記し、側壁313の壁面領域Eに相当する発吸熱部位iを適宜発吸熱部位Eと表記し、窓331の領域Fに相当する流入出領域iを適宜流入出領域Fと表記し、窓332の領域Gに相当する流入出領域iを適宜流入出領域Gと表記する。
In this application target, the fluid system to be estimated flows into the
図14は、実施の形態1における温度推定装置10の機能構成を示すブロック図である。図14に示すように、温度推定装置10は、入力部11と、表示部12と、記憶部13と、制御部14とを備え、部屋3内に設置された温度計34−1〜34−4からの温度実測値が制御部14に入力される構成となっている。
FIG. 14 is a block diagram illustrating a functional configuration of the
入力部11は、ユーザが温度推定に必要な情報の入力等の各種操作を行うためのものであり、入力信号を制御部14に出力する。この入力部11は、キーボードやマウス、タッチパネル等によって実現される。表示部12は、LCDやELディスプレイ等の表示装置によって実現され、制御部14の制御のもと、温度推定の結果等を画面表示する。
The
記憶部13は、更新記録可能なフラッシュメモリ等のROMやRAMといった各種ICメモリ、内蔵或いはデータ通信端子で接続されたハードディスク、CD−ROM等の情報記録媒体およびその読取装置等によって実現されるものであり、温度推定装置10を動作させ、この温度推定装置10が備える種々の機能を実現するためのプログラムや、このプログラムの実行中に使用されるデータ等が記録される。この記憶部13には、部屋3内に設定される温度推定点jの重みWijを登録した重みデータベースや、温度推定点jの推定温度Tejを該当する温度推定点jの部屋3内における位置と対応付けて設定した温度データ等が保存される。
The
制御部14は、CPU等のハードウェアによって実現される。この制御部14は、入力部11から入力される入力信号、記憶部13に記録されるプログラムやデータ等をもとに温度推定装置10を構成する各部への指示やデータの転送等を行い、温度推定装置10全体の動作を統括的に制御する。この制御部14は、温度推定部141と、温度分布表示処理部143とを含む。
The
温度推定部141は、図4〜図6に示した処理手順に従って重み算出処理を行うことで、温度推定点jに対する部位i毎の重みWijを算出し、算出した温度推定点jについての重みWijを重みデータベースとして記憶部13に保存する。具体的には、温度推定部141は、例えば、数値流体シミュレーションとして有限体積法を用い、乱流モデルとして標準k−ε乱流モデルを利用して勢力(R1ij,R2ij)を取得し、重みWijを算出する。実施の形態1では、流体系の高さ方向の流れは無視できるものとして部屋3内の流体系の2次元の温度分布を推定することとし、重みWijの算出においても、2次元モデルで近似を行うこととする。
The
この場合には、温度推定部141は、図4のステップS3の処理として、2次元の流れ場を実際の流れ場として計算する。図15は、ここで計算される部屋3の実際の流れ場を示す模式図である。図15に示すように、流れ場の計算では、部屋3内の全域における空気の流れを表す流速ベクトルV3、具体的には、図15中に矢印A3で示すように窓331(図13を参照)から室内に流入して窓332(図13を参照)から室外に流出する空気の部屋3内の各位置での流れの向きおよびその流速を表す流速ベクトルV3が得られる。
In this case, the
温度推定部141は、図4のステップS7の処理として、温度推定点jを等間隔で部屋3内の全域に設定する。温度推定部141は、ステップS9〜ステップS23の処理として、図15に示す実際の流れ場を用い、前述のように部屋3内の全域に設定した各温度推定点jにおける流体の各部位iの下流側勢力R1ijを取得する。図16−1〜図16−7は、それぞれ各温度推定点jにおける各部位A〜Gの下流側勢力R1ijを等値線図化して示した図である。
The
その後、温度推定部141は、図5のステップS25の処理として、図15に示す実際の流れ場の各流速ベクトルV3の向きを逆向きにした反転流れ場を計算する。温度推定部141は、図5のステップS27〜ステップS45の処理として、反転流れ場を用い、部屋3内の全域に設定した各温度推定点jにおける流体の各部位iの上流側勢力R2ijを取得する。図17−1〜図17−7は、それぞれ温度推定点jにおける各部位A〜Gの上流側勢力R2ijを等値線図化して示した図である。
After that, the
温度推定部141は、図6のステップS47〜ステップS59の処理として、式(5)の重み関数W(R1ij,R2ij)を用いた重みWijの算出と、式(8)の重み関数W(R1ij,R2ij)を用いた重みWijの算出とをそれぞれ行い、用いた重み関数W(R1ij,R2ij)毎に各温度推定点jの重みWijをデータベース化して記憶部13に保存する。図18−1〜図18−7は、それぞれ式(5)の重み関数W(R1ij,R2ij)を用いて算出した各温度推定点jに対する部位A〜G毎の重みWijを等値線図化して示した図である。図19−1〜図19−7は、それぞれ式(8)の重み関数W(R1ij,R2ij)を用いて算出した各温度推定点jに対する部位A〜G毎の重みWijを等値線図化して示した図である。実施の形態1では、発吸熱部位Eおよび流入出部位F,Gの温度を未知としており、図18−5〜図18−7や図19−5〜図19−7に示すように、部位E,F,Gの重みWijは「0」として算出される。
The
温度推定部141は、図7に示した処理手順に従って温度推定処理を行うことで、温度計34−1〜34−4が計測する温度実測値である温度実測部位i(A〜D)の既知温度Tiをもとに、各温度推定点jについての重みWijを用いることで各温度推定点jの推定温度Tejを算出する。そして、温度推定部141は、算出した各温度推定点jの推定温度Tejを、温度データとして記憶部13に保存する。
The
温度分布表示処理部143は、温度推定部141が推定して記憶部13に保存した温度データを参照し、各温度推定点jの推定温度Tej、すなわち、既知温度Tiである温度実測部位i(A〜D)の温度実測値を補間した部屋3内の流体系全体の温度分布を例えば等値線図化し、温度分布モニタリング画面として表示部12に表示する。
Temperature distribution display processing unit 143 refers to the temperature data
以上説明した構成の温度推定装置10において、上記した式(5)の重み関数W(R1ij,R2ij)によって算出した重みWijを用いた場合(実験例1)と、式(8)の重み関数W(R1ij,R2ij)によって算出した重みWijを用いた場合(実験例2)とで、それぞれ各温度推定点jにおける汚染物質の推定温度Tejを算出した。また、比較例として、従来法である逆距離加重法を用いた重みの算出を行い、得られた重みを用いて各温度推定点jの推定温度Tejを算出した。比較例での重みの算出は、図4〜図6の重み算出処理と同様の処理手順で行い、ステップS53において、次式(15)で示す、逆距離補間の式の重みWij´を用いる。lijは、温度実測部位iと温度推定点jとの直線距離である。uは、補間パラメータであり、例えば、u=2として重みWijを算出した。
In the
具体的には、温度推定部141が、実験例1,2および比較例での3種類の重みWijをそれぞれ用いて各温度推定点jの推定温度Tejを算出し、温度分布表示処理部143が、実験例1,2および比較例での各温度推定点jの推定温度Tejを等値線図化することで、実験例1,2および比較例のそれぞれについての推定結果を得た。推定に用いた各温度実測部位A〜Dの既知温度Ti、すなわち、部屋3内の対応する設置位置A〜Dに設置された温度計34−1〜34−4の温度実測値を表1に示す。
Specifically, the
さらに比較のため、数値流体解析を用い、各温度実測部位A〜Dの既知温度Tiを表1に示す温度実測値として、部屋3内の真の温度分布を算出した。
For further comparison, using a numerical fluid analysis, a known temperature T i of each temperature measured sites A~D as actually measured temperature values shown in Table 1, were calculated the true temperature distribution in the
図20は、実施の形態1における実験例1の推定結果、すなわち、部屋3内の流体系の温度分布を等値線図化した図である。図21は、実施の形態1における実験例2の推定結果を示す図であり、図22は、実施の形態1における比較例の推定結果を示す図である。図23は、部屋3内の真の温度分布を示す図である。実験例1,2と比較例とを比較すると、図22に示すように、比較例では、温度実測部位A〜Dの周囲に等値線が同心円状に広がる温度分布が推定されている。このように、比較例では、部屋3内の空気の流れが温度推定に反映されず、図23に示す部屋3内の真の温度分布に対応しない推定結果となった。これに対し、図20,21に示すように、実施の形態1の実験例1,2では、どちらも等値線が部屋3内の空気の流れの方向に沿って長く伸びた温度分布が推定されており、図23に示す部屋3内の真の温度分布に対応した推定結果が得られた。このように、実施の形態1によれば、部屋3内の空気の流れを反映させた温度推定が実現でき、部屋3内の温度分布を精度よく再現することができた。
FIG. 20 is an isoline diagram of the estimation result of Experimental Example 1 in the first embodiment, that is, the temperature distribution of the fluid system in the
(実施の形態2)
次に、実施の形態2として、水槽を適用対象とし、この水槽内を流れる流体系の温度推定および温度分布の可視化について説明する。図24は、実施の形態2において適用対象とする水槽4の内部を側方から示した模式図である。図25は、図14の水槽4の内部を上方から示した模式図である。
(Embodiment 2)
Next, as
図24および図25に示す水槽4は、例えば、奥行方向(図25の上下方向)が1(m)、幅方向(図24および図25の左右方向)が1(m)、深さ(図24の上下方向の幅)が0.5(m)の直方体形状を有し、水槽4内は常に水が満たされる構成となっている。すなわち、水槽4の上面には、図25に向かって左側の両角において水槽4の内部空間と連通する2本のパイプ41,42が設けられ、このパイプ41,42から水槽4内に水が注入されるようになっている。一方、水槽4の底面には、図25に向かって右側中央において水槽4の内部空間と連通する1本のパイプ43が設けられ、このパイプ43からは、パイプ41,42から流入した水の総量と同量の水が流出するようになっている。
24 and 25, for example, the depth direction (vertical direction in FIG. 25) is 1 (m), the width direction (horizontal direction in FIGS. 24 and 25) is 1 (m), and the depth (FIG. 24 has a rectangular parallelepiped shape with a vertical width of 0.5 (m), and the
水槽4の内部空間には、水槽4の幅方向に沿って幅方向の半分を仕切る仕切り板44が配設されており、その奥行方向の位置が、水槽4の奥行方向中央を通る垂直断面S4に対して0.2(m)だけパイプ41側に寄った配置になっている。水槽4内の図24および図25中に「×」を付して示す6箇所P41〜P46において、温度を推定するための温度計測装置としての温度計45−1〜45−6が設置されている。温度計45−1〜45−6の深さ方向の位置は、水槽4のちょうど中央深さとなる位置とした。
In the internal space of the
本適用対象では、推定対象の流体系が水槽4内を流れる水であり、温度計45−1〜45−6の設置位置P41〜P46が温度実測部位iとなる。パイプ41,42の下端が流入部位、パイプ43の上端が流出部位であり、これらが流入出部位iとなる。ただし、パイプ41〜43の流路全域を流入出部位iとしてもよいし、パイプ41,42の上端やパイプ43の下端を流入出部位iとしてもよい。例えば、パイプ41からは10(℃)の水が注入され、パイプ42からは50(℃)の水が注入される。ただし、温度推定を行う際には、流入出部位iの温度は未知とする。本適用対象の流体系は、水槽4内に満たされる水の水面や水槽4の内壁面における伝熱が十分に小さいものとし、発吸熱部位iを含まないこととする。図24,図25では図示しないが、実施の形態2においても、これら温度実測部位iおよび流入出部位iについて、それぞれ対応する温度実測領域iおよび流入出領域iが設定される。
In this application target, the fluid system to be estimated is water flowing in the
図26は、実施の形態2における温度推定装置10aの機能構成を示すブロック図である。図26において、実施の形態1と同様の構成には同一の符号を付する。図26に示すように、温度推定装置10aは、入力部11と、表示部12と、記憶部13と、制御部14aとを備え、水槽4内に設置された温度計45−1〜45−6からの温度実測値が制御部14aに入力される構成となっている。
FIG. 26 is a block diagram illustrating a functional configuration of the
記憶部13は、水槽4内に設定される温度推定点jの重みWijを登録した重みデータベースや、温度推定点jの推定温度Tejを該当する温度推定点jの水槽4内における位置と対応付けて設定した温度データ等を保存する。
The
制御部14aは、温度推定部141と、温度データ抽出部142aと、温度分布表示処理部143aとを含む。
The control unit 14a includes a
温度推定部141は、図4〜図6に示した処理手順に従って重み算出処理を行うことで、温度推定点jに対する部位i毎の重みWijを算出し、算出した温度推定点jについての重みWijを重みデータベースとして記憶部13に保存する。例えば、温度推定部141は、数値流体シミュレーションとして有限体積法を用い、乱流モデルとして標準k−ε乱流モデルを利用して勢力(R1ij,R2ij)を取得し、重みWijを算出する。このとき、温度推定部141は、図4のステップS7の処理として、奥行方向、幅方向、および深さ方向に沿って0.04(m)間隔で温度推定点jを水槽4内の全域に設定する。また、温度推定部141は、図6のステップS47〜ステップS59の処理として、式(5)の重み関数W(R1ij,R2ij)を用いた重みWijの算出と、式(8)の重み関数W(R1ij,R2ij)を用いた重みWijの算出とをそれぞれ行い、用いた重み関数W(R1ij,R2ij)毎に各温度推定点jの重みWijをデータベース化して記憶部13に保存する。実施の形態2では、流入出部位iの温度を未知としており、この流入出部位iの重みWijは「0」として算出される。
The
温度推定部141は、実施の形態1と同様に、図7に示した処理手順に従って温度推定処理を行うことで、既知温度Tiである温度実測部位iの温度実測値をもとに各温度推定点jの推定温度Tejを算出し、温度データとして記憶部13に保存する。
Similarly to the first embodiment, the
温度データ抽出部142aは、温度推定部141が推定して記憶部13に保存した温度データを参照し、水槽4の任意の断面における推定温度Tej(任意の断面内の温度推定点jの推定温度Tej)を抽出する。推定温度Tejを抽出する断面は、予め固定的に設定しておく構成としてもよいし、ユーザ操作に従って決定することとしてもよい。ユーザ操作に従って決定する場合には、温度データ抽出部142aは、入力部11を介してユーザによる断面の指定操作を受け付け、ユーザが指定した断面における推定温度Tejを抽出する。
The temperature data extraction unit 142 a refers to the temperature data estimated by the
温度分布表示処理部143aは、温度データ抽出部142aが抽出した任意の断面における推定温度Tejを例えば等値線図化することでこの任意の断面における温度分布を可視化し、温度分布モニタリング画面として表示部12に表示する。
The temperature distribution display processing unit 143a visualizes the temperature distribution in the arbitrary cross section by, for example, converting the estimated temperature Te j in the arbitrary cross section extracted by the temperature data extracting unit 142a into an isometric diagram, and serves as a temperature distribution monitoring screen. It is displayed on the
以上説明した構成の温度推定装置10aにおいて、上記した式(5)の重み関数W(R1ij,R2ij)によって算出した重みWijを用いた場合(実験例1)と、式(8)の重み関数W(R1ij,R2ij)によって算出した重みWijを用いた場合(実験例2)とで、それぞれ各温度推定点jの推定温度Tejを算出した。また、比較例として、従来法である逆距離加重法を用いた重みの算出を行い、得られた重みを用いて各温度推定点jの推定温度Tejを算出した。比較例での重みの算出は、図4〜図6の重み算出処理と同様の処理手順で行い、ステップS53において式(15)に示す重みWij´を用いる。
In the
具体的には、温度推定部141が、実験例1,2および比較例での3種類の重みWijをそれぞれ用いて各温度推定点jの推定温度Tejを算出し、温度データ抽出部142aが、実験例1,2および比較例での各温度推定点jの推定温度Tejから水槽4の深さ方向の中央を通る水平断面における推定温度Tejを抽出し、温度分布表示処理部143aが、前述の水平断面における推定温度Tejを等値線図化することで、実験例1,2および比較例のそれぞれについての推定結果を得た。推定に用いた各温度実測部位iの既知温度Ti、すなわち、水槽4内の対応する設置位置P41〜P46に設置された温度計45−1〜45−6の温度実測値を表2に示す。
Specifically, the
図27は、実施の形態2における実験例1の推定結果、すなわち、水槽4の深さ方向の中央を通る水平断面における温度分布を等値線図化した図である。図28は、実施の形態2における実験例2の推定結果を示す図であり、図29は、実施の形態2における比較例の推定結果を示す図である。
FIG. 27 is an isoline diagram of the estimation result of Experimental Example 1 in the second embodiment, that is, the temperature distribution in a horizontal section passing through the center of the
従来法である逆距離加重法を用いた比較例の場合、温度実測部位iと温度推定点jとの距離という幾何学的な情報のみを指標として重みを算出するが、この方法では、流体設備内の構造が考慮されないため、実際の温度分布と大きく異なる推定結果が得られる場合があった。すなわち、例えば本適用対象の水槽4の内部空間に配設された仕切り板44のような流体の流れを遮る部材を備えた流体設備に適用した場合、仕切り板44を越えて連続的な温度分布が推定されてしまう場合があった。しかしながら、実際には、仕切り板44によって流体の流れが遮られるため、この仕切り板44を境に温度が不連続となることがある。
In the case of the comparative example using the inverse distance weighting method which is the conventional method, the weight is calculated using only the geometric information of the distance between the temperature measurement site i and the temperature estimation point j as an index. Because the internal structure is not taken into account, there are cases where an estimation result that is significantly different from the actual temperature distribution is obtained. That is, for example, when applied to a fluid facility having a member that blocks the flow of fluid, such as a
実際に、実験例1,2と比較例とを比較すると、図27および図28に示すように、実施の形態2の実験例1,2の推定結果では、仕切り板44を境に温度分布が不連続となっており、仕切り板44によって整流された流れ場の影響を反映した温度分布が推定できている。一方、比較例の推定結果では、図29に示すように、仕切り板44を乗り越えて連続的な温度分布となっており、仕切り板44によって整流された流れ場の影響が温度分布に反映されていない。このように、実施の形態2では、従来法である逆距離加重法を用いた場合と異なり、流れ場の影響を反映した温度推定が実現でき、水槽4内の温度分布を精度よく再現することが確認できた。したがって、実施の形態2によれば、複雑な3次元の流れ場を有する流体系を推定対象とする場合においても、高精度な温度推定が実現できる。
Actually, 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 is separated from the
水槽4内の温度推定精度を定量的に検証するため、水槽4内にさらに温度計を追加設置して実験を行った。図30は、水槽4内に追加設置した温度計45−7〜45−9の設置位置を示す図である。図30に示すように、水槽4内に「×」を付して示す3箇所P47〜P49に温度計45−7〜45−9を追加設置し、各設置位置P47〜P49の温度実測値を取得するとともに、各設置位置P47〜P49を温度実測部位iとして実験例1,2および比較例による温度推定を行った。温度計45−7〜45−9の深さ方向の位置は、温度計45−7〜45−9と同様に、水槽4のちょうど中央深さとなる位置とした。
In order to quantitatively verify the temperature estimation accuracy in the
温度計45−7〜45−9の設置位置P47〜P49である各温度実測部位iにおける温度実測値、実施の形態2の実験例1による推定温度Tej、実験例2による推定温度Tej、および比較例による推定温度Tejを表3に示す。表3に示すように、実施の形態2の実験例1,2で得た推定温度Tejは、比較例で得た推定温度Tejと比べて温度実測値に近い値となっており、実施の形態2による温度推定精度が優れていることが確認できた。
Temperature measured value at each temperature measured portions i, which is the installation position P47~P49 thermometer 45-7~45-9, estimated temperature of the Experimental Example 1 of the second embodiment Te j, the estimated temperature Te j of the Experimental Example 2, Table 3 shows the estimated temperature Te j according to the comparative example. As shown in Table 3, the estimated temperature Te j obtained in Experimental Examples 1 and 2 of the second embodiment is closer to the actually measured temperature value than the estimated temperature Te j obtained in the comparative example. It was confirmed that the temperature estimation accuracy of
(実施の形態3)
実施の形態3として、実施の形態2と同じ水槽で、流入水温が時間変化する場合の温度分布推定および可視化について説明する。図24および図25に示すパイプ41からは常に10℃の水が流入する。パイプ42からは一定流量の水が流入し、最初は水温が10℃、途中から水温が50℃になる。温度計45−1〜45−6は、実施の形態2と同じ位置(P41、P42、P43、P44、P45、P46)に配置される。
(Embodiment 3)
As
実施の形態3において、下流側勢力R1ij、上流側勢力R2ij、下流側伝達時間τ1ij、上流側伝達時間τ2ij、下流側重みW1ij、および上流側重みW2ijを実施の形態2と同様に数値流体シミュレーションを用いて算出した。数値流体シミュレーションは有限体積法を用い、乱流モデルとして標準k−ε乱流モデルを利用した。流れ場計算では、パイプ41上端から流量0.765L/sで水が流入し、パイプ42上端から流量1.531L/sで水が流入し、パイプ43下端では圧力一定で流出することとし、水槽4の上面は滑り条件、側壁、および底壁は壁の対数則を用いた壁境界条件として境界条件を与えて計算を行った。伝達時間の算出は、水の初期温度27℃、発熱量2,200kW、閾値温度28℃として計算した。温度推定点jは0.04m間隔で配置し、水槽4内全域に配置した。
In the third embodiment, the downstream force R 1ij , the upstream force R 2ij , the downstream transmission time τ 1ij , the upstream transmission time τ 2ij , the downstream weight W 1ij , and the upstream weight W 2ij are the same as those in the second embodiment. Similarly, calculation was performed using a numerical fluid simulation. The numerical fluid simulation uses a finite volume method and a standard k-ε turbulent model as a turbulent model. In the flow field calculation, water flows from the upper end of the
重み関数W1ij、W2ijの算出には、式(14a)および式(14b)を用いた。温度計45−1〜45−6により水槽4内の各位置(P41、P42、P43、P44、P45、P46)において測定された温度の時間推移を図31−1〜図31−6に示す。
Formulas (14a) and (14b) were used for calculating the weight functions W 1ij and W 2ij . The time transition of the temperature measured at each position (P41, P42, P43, P44, P45, P46) in the
温度推定を行う時間t0として、パイプ42の水温が50℃に変わった時間から1分後、2分後、3分後、4分後、5分後、および6分後の6つの時点を考える。水槽4の中央を通る水平断面における温度分布について、上記位置P41〜P46において測定された時系列温度データTi(t)から上記で算出した下流側伝達時間τ1ijおよび上流側伝達時間τ2ijを用いて、下流側既知温度Ti(t0−τ1ij)および上流側既知温度Ti(t0+τ2ij)を抽出した。そして式(14a)と式(14b)とを用いて、水温Tejを推定した。
As the time t 0 when the temperature is estimated, six points of time after 1 minute, 2 minutes, 3 minutes, 4 minutes, 5 minutes, and 6 minutes after the time when the water temperature of the
パイプ42から流入する水温が50℃に変わった時間から1分後、2分後、3分後、4分後、5分後、および6分後の水槽4の中央を通る水平断面における水温について等値線図化した。等値線図を図32−1〜図32−6に示す。パイプ42から流入する水の温度が10℃から50℃に変わると、パイプ42に近い位置から徐々に温度が上昇していく様子がうまく現れており、温度分布の時間推移が有る場合でも温度分布を推定することができた。
About the water temperature in a horizontal section passing through the center of the
(実施の形態4)
次に、実施の形態4として、溶融亜鉛めっきポットを適用対象とし、この溶融亜鉛めっきポット内を流れる流体系の温度推定および温度分布の可視化について説明する。図33は、実施の形態4において適用対象とする溶融亜鉛めっきポット5の内部を側方から示した模式図である。自動車や建材などに利用される亜鉛めっき鋼板を製造する鉄鋼プロセスのひとつである溶融亜鉛めっきラインでは、図33に例示するような溶融亜鉛めっきポット5において鋼板51を溶融亜鉛中に浸漬させた後、不図示の付着量制御装置でそのめっき付着量を調整し、冷却等の所定の後処理を施してめっき鋼板とする。操業条件は、例えばライン速度を120(mpm)とし、鋼板の板幅を1,500(mm)とする。
(Embodiment 4)
Next, as a fourth embodiment, a hot dip galvanizing pot is an application target, and temperature estimation of a fluid system flowing in the hot dip galvanizing pot and visualization of a temperature distribution will be described. FIG. 33 is a schematic diagram showing the inside of hot
実施の形態4で適用対象とする図33の溶融亜鉛めっきポット5の溶融亜鉛の容量は、例えば250(t)であり、溶融亜鉛めっきポット5内は溶融亜鉛で満たされている。この溶融亜鉛めっきポット5は、図33の紙面と平行な対向する内側壁面のそれぞれに設置された誘導加熱装置52を備える。溶融亜鉛めっきポット5は、内部空間に亜鉛インゴット53を投入するためのインゴット投入部(不図示)を備えている。誘導加熱装置52は、このインゴット投入部に投入される亜鉛インゴット53を溶解して溶融亜鉛とし、この溶解した溶融亜鉛の温度を所定の温度に維持するためのものである。
The hot
溶融亜鉛めっきポット5の内部空間にはシンクロール54が設置されており、このシンクロール54によって、溶融亜鉛中に浸漬されてこの溶融亜鉛中で搬送される鋼板51の通板方向が方向転換されるようになっている。鋼板51への付着によって消費される亜鉛は、インゴット投入部(不図示)に対する亜鉛インゴット53の投入により補給される。溶融亜鉛めっきポット5内の図31中に「×」を付して示す8箇所P51〜P58において、温度を推定するための温度計測装置としての熱電対55−1〜55−8が設置されている。各熱電対55−1〜55−8は、溶融亜鉛めっきポット5の図33の紙面と平行な一方の内側壁面、例えば、シンクロール54から見て紙面手前側の内側壁面からの距離が300(mm)となる面内において、各々が図33に示す位置関係で設置されている。
A
本適用対象では、推定対象の流体系が溶融亜鉛めっきポット5内に満たされる溶融亜鉛であり、熱電対55−1〜55−8の設置位置P51〜P58が温度実測部位iとなる。実施の形態4では、発熱が生じる流体系内の部位である誘導加熱装置52の加熱位置およびインゴット投入部の2箇所が発吸熱部位iとなる。誘導加熱装置52の加熱位置に相当する発吸熱部位iの温度を既知とし、具体的には、487.72(℃)とする(表4を参照)。一方、インゴット投入部に相当する発吸熱部位iの温度については、未知とする。本適用対象では、推定対象の流体系の流入出、すなわち、溶融亜鉛めっきポット5内への溶融亜鉛の流入および溶融亜鉛めっきポット5外への溶融亜鉛の流出は存在しないため、実施の形態4の流体系は、流入出部位iを含まないこととする。実施の形態4では、これら温度実測部位iおよび発吸熱部位iについて、それぞれ対応する温度実測領域iおよび発吸熱領域i(不図示)が設定される。
In this application target, the fluid system to be estimated is molten zinc filled in the hot
図34は、実施の形態4における温度推定装置10bの機能構成を示すブロック図である。図34において、実施の形態1と同様の構成には同一の符号を付する。図34に示すように、温度推定装置10bは、入力部11と、表示部12と、記憶部13と、制御部14bとを備え、溶融亜鉛めっきポット5内に設置された熱電対55−1〜55−8からの温度実測値が制御部14bに入力される構成となっている。
FIG. 34 is a block diagram illustrating a functional configuration of the
記憶部13は、溶融亜鉛めっきポット5内に設定される温度推定点jの重みWijを登録した重みデータベースや、温度推定点jの推定温度Tejを該当する温度推定点jの溶融亜鉛めっきポット5内における位置と対応付けて設定した温度データ等を保存する。
The
制御部14bは、温度推定部141と、温度データ抽出部142bと、温度分布表示処理部143bとを備える。
The control unit 14b includes a
温度推定部141は、図4および図5に示した処理手順に従って重み算出処理を行うことで、温度推定点jに対する部位i毎の重みWijを算出し、算出した温度推定点jについての重みWijを重みデータベースとして記憶部13に保存する。例えば、温度推定部141は、数値流体シミュレーションとして有限体積法を用い、乱流モデルとして標準k−ε乱流モデルを利用して勢力(R1ij,R2ij)を取得し、重みWijを算出する。このとき、温度推定部141は、図4のステップS7の処理として、温度推定点jを等間隔で溶融亜鉛めっきポット5内の全域に設定する。また、温度推定部141は、図6のステップS47〜ステップS59の処理として、例えば式(5)の重み関数W(R1ij,R2ij)を用いた重みWijの算出を行い、各温度推定点jの重みWijをデータベース化して記憶部13に保存する。実施の形態4では、インゴット投入部に相当する発吸熱部位iの温度を未知としており、このインゴット投入部に相当する発吸熱部位iの重みWijは「0」として算出される。
The
温度推定部141は、実施の形態1と同様に、図7に示した処理手順に従って温度推定処理を行うことで、既知温度Tiである温度実測部位iの温度実測値と誘導加熱装置52の加熱位置に相当する発吸熱部位iの温度とをもとに各温度推定点jの推定温度Tejを算出し、温度データとして記憶部13に保存する。
Similarly to the first embodiment, the
温度データ抽出部142bは、温度推定部141が推定して記憶部13に保存した温度データを参照し、溶融亜鉛めっきポット5の任意の断面における推定温度Tejを抽出する。温度分布表示処理部143bは、温度データ抽出部142bが抽出した任意の断面における推定温度Tejを例えば等値線図化することでこの任意の断面における温度分布を可視化し、温度分布モニタリング画面として表示部12に表示する。
Temperature data extraction unit 142b refers to the temperature data
以上説明した構成の温度推定装置10bにおいて、各温度推定点jの温度推定を行った。具体的には、温度推定部141が、各温度推定点jの推定温度Tejを算出し、温度データ抽出部142bが、溶融亜鉛めっきポット5の例えば内側壁面から300(mm)の鉛直断面における推定温度Tejを抽出し、温度分布表示処理部143bが、前述の鉛直断面における推定温度Tejを等値線図化することで推定結果を得た。推定に用いた各温度実測部位iの既知温度Ti、すなわち、溶融亜鉛めっきポット5内の対応する設置位置P51〜P58に設置された熱電対55−1〜55−8の温度実測値を、誘導加熱装置52の加熱位置の既知温度Tiとともに表4に示す。
In the
図35は、実施の形態4における推定結果、すなわち、溶融亜鉛めっきポット5の内側壁面から300mmの鉛直断面における温度分布を等値線図化した図である。この実施の形態4によれば、溶融亜鉛めっきポット5内での溶融亜鉛の流動の影響を反映した温度推定が実現でき、溶融亜鉛の流動効果を考慮して温度分布を高精度に推定できる。また、各温度推定点jについての重みWijを重みデータベースとして記憶部13に保存しておくこととしたので、温度推定および温度分布の可視化の際は、既知温度Tiをもとに重み付き平均処理を行うだけでよく、計算時間を1秒以内とすることができた。よってオンライン(リアルタイム)での温度分布の可視化も可能である。なお、溶融亜鉛めっきポット5内の溶融亜鉛の温度が所定範囲内にない場合、溶融亜鉛めっき鋼板に表面欠陥が発生することが知られている。従って、上述の処理によって溶融亜鉛めっきポット内の溶融亜鉛の温度を推定し、推定結果に基づいて溶融亜鉛めっきポット内の溶融亜鉛の温度が所定範囲内となるように誘導加熱装置52を制御することによって、表面欠陥のない溶融亜鉛めっき鋼板を製造することができる。
FIG. 35 is an isoline diagram of the estimation result in the fourth embodiment, that is, the temperature distribution in a vertical cross section of 300 mm from the inner wall surface of the hot
具体的には、図34に示すように、制御部14bは、溶融亜鉛めっきポット5内の所定の領域における溶融亜鉛の温度が所定の閾値内にあるか否かを判定する判定部144bと、溶融亜鉛めっきポット5の誘導加熱装置52の出力を操作して溶融亜鉛の温度を制御する温度制御部145bとを備えている。なお、「溶融亜鉛めっきポット5内の所定の領域」とは、例えば表面欠陥に影響を与える鋼板51の表面と溶融亜鉛が接触する箇所や、シンクロール54と溶融亜鉛が接触する箇所や、シンクロール54上部と鋼板51とで囲まれた領域などをいう。溶融亜鉛温度の閾値は予め判定部144bに入力されるか、操作者により入力部11を介して入力され、判定部144bは、温度データ抽出部142bが抽出した、所定の領域における溶融亜鉛温度が閾値内にあるか否かを判定する。判定部144bが、所定の領域における溶融亜鉛温度が閾値範囲外と判定した場合、温度制御部145bは、所定の領域における溶融亜鉛の温度が閾値範囲内となるよう誘導加熱装置52の出力を操作する。本実施の形態4によれば、温度制御部145bが誘導加熱装置52を制御することにより所定の領域における溶融亜鉛温度を制御することが可能となる。これにより、鋼板51の表面欠陥を防止できる。
Specifically, as shown in FIG. 34, the control unit 14b determines whether or not the temperature of the molten zinc in a predetermined region in the hot
(実施の形態5)
次に、実施の形態5として、連続鋳造用のタンディッシュを適用対象とし、このタンディッシュ内を流れる流体系の温度推定および温度分布の可視化について説明する。図36は、実施の形態5において適用対象とするタンディッシュ6の構成を模式的に示した斜視図である。図37は、実施の形態5のタンディッシュ6に設置される熱電対64−1〜64−5の設置位置を示す図であり、タンディッシュ6の長辺側方の図37に向かって右側半分について、その内部の様子を模式的に示している。
(Embodiment 5)
Next, as a fifth embodiment, a tundish for continuous casting is applied, and temperature estimation and visualization of a temperature distribution of a fluid system flowing in the tundish will be described. FIG. 36 is a perspective view schematically showing the configuration of the
図36に示すタンディッシュ6は、奥行方向が1(m)、幅方向が8(m)、高さが1(m)の直方体形状を有し、内部に溶鋼を収容する。図36において、タンディッシュ6に収容される溶鋼の液面S7を破線で示している。このタンディッシュ6は、 取鍋からの溶鋼を注入するノズル61と、底部の2箇所に設けられ、溶鋼を鋳型へと導入するための流出孔62,62と、溶鋼を加熱して温度を制御する2つのプラズマ加熱装置63,63とを備える。タンディッシュ6は、取鍋からの溶鋼を注入するノズル61が中央上部に設けられ、鋳型への流出孔 62,62が幅方向の両端に設けられた2ストランド仕様である。
The
実施の形態5では、図37に示すタンディッシュ6の長辺方向の右側半分において、短辺方向中央を通る鉛直面内の図37中に「×」を付して示す5箇所P61〜P65に、温度を推定するための温度計測装置としての熱電対64−1〜64−5が設置されている。タンディッシュ6の長辺方向右側だけに注目し、右側の5箇所に熱電対64−1〜64−5を設置したのは、タンディッシュ6が左右対称の構造を有するためであるが、長辺方向左側にも同様に熱電対64−1〜64−5を設置し、温度推定に用いるようにしてもよい。
In the fifth embodiment, in the right half of the long side direction of the
本適用対象では、推定対象の流体系が、タンディッシュ6の内部に収容される溶鋼、具体的には、ノズル61の下端からタンディッシュ6の内部に流入し、流出孔62,62からタンディッシュ6の外部(鋳型)へと流出する溶鋼である。本適用対象では、熱電対64−1〜64−5の設置位置P61〜P65が温度実測部位iとなる。プラズマ加熱装置63,63の加熱位置が発吸熱部位iとなり、ノズル61の下端および流出孔62,62が流入出部位i(ノズル61の下端が流入部位、流出孔62,62が流出部位)となる。溶鋼の液面S7は、外部から強い冷却を受けるので発吸熱部位iとなる。プラズマ加熱装置63,63の加熱位置に相当する発吸熱部位iとノズル61の下端の流入位置に相当する流入出部位iは温度が既知であり、流出孔62,62に相当する流入出部位i、溶鋼の液面S7に相当する発吸熱部位iの温度は未知である。
In this application target, the estimation target fluid system flows into the
図38は、実施の形態5における温度推定装置10cの機能構成を示すブロック図である。図38において、実施の形態1と同様の構成には同一の符号を付する。図38に示すように、温度推定装置10cは、入力部11と、表示部12と、記憶部13と、制御部14cとを備え、タンディッシュ6内に設置された熱電対64−1〜64−5からの温度実測値が制御部14cに入力される構成となっている。
FIG. 38 is a block diagram illustrating a functional configuration of the
記憶部13は、タンディッシュ6内に設定される温度推定点jの重みWijを登録した重みデータベースや、温度推定点jの推定温度Tejを該当する温度推定点jのタンディッシュ6内における位置と対応付けて設定した温度データ等を保存する。
The
制御部14cは、温度推定部141と、温度データ抽出部142cと、温度分布表示処理部143cと、判定部144cと、温度制御部145cとを含む。
The
温度推定部141は、図4〜図6に示した処理手順に従って重み算出処理を行うことで、温度推定点jに対する部位i毎の重みWijを算出し、算出した温度推定点jについての重みWijを重みデータベースとして記憶部13に保存する。例えば、温度推定部141は、数値流体シミュレーションとして有限体積法を用い、乱流モデルとして標準k−ε乱流モデルを利用して勢力(R1ij,R2ij)を取得し、重みWijを算出する。このとき、温度推定部141は、図4のステップS7の処理として、温度推定点jを等間隔でタンディッシュ6内の全域に設定する。実施の形態5では、プラズマ加熱装置63,63の加熱位置に相当する発吸熱部位iおよび流出孔62,62に相当する流入出部位iの温度を未知としており、これらの部位iの重みWijは「0」として算出される。
The
温度推定部141は、実施の形態1と同様に、図7に示した処理手順に従って温度推定処理を行うことで、既知温度Tiである温度実測部位iの温度実測値をもとに各温度推定点jの推定温度Tejを算出し、温度データとして記憶部13に保存する。
Similarly to the first embodiment, the
温度データ抽出部142cは、温度推定部141が推定して記憶部13に保存した温度データを参照し、タンディッシュ6の任意の断面における推定温度Tejを抽出する。タンディッシュ6内の溶鋼は、その浴面やタンディッシュ6の内側壁面との接触部分において冷却されるため、取鍋から注入された溶鋼は、下方に流れて流出孔62,62に近づくに従って温度が低下していく。タンディッシュ6は、内側壁面が耐火物で覆われており、この耐火物は常に高温の溶鋼と接触している。この耐火物と接触している溶鋼の温度が急激に変化すると、耐火物に大きな熱応力が発生し、耐火物損傷の問題が起こる。したがって、内側壁面と接触する溶鋼の温度が所定の閾値内となるように溶鋼の温度を制御することが好ましい。そこで、実施の形態5では、温度データ抽出部142cは、例えばタンディッシュ6の長辺方向の内側壁面近傍、例えば図37の紙面手前側の内側壁面近傍の鉛直断面における推定温度Tejを抽出する。
The temperature data extraction unit 142 c refers to the temperature data estimated by the
温度分布表示処理部143cは、温度データ抽出部142cが抽出した任意の断面(例えばタンディッシュ6の長辺方向の内側壁面近傍の鉛直断面)における推定温度Tejを例えば等値線図化することでこの任意の断面における温度分布を可視化し、温度分布モニタリング画面として表示部12に表示する。
The temperature distribution display processing unit 143c, for example, plots the estimated temperature Te j in an arbitrary cross section (for example, a vertical cross section in the vicinity of the inner wall surface in the long side direction of the tundish 6) extracted by the temperature data extraction unit 142c. The temperature distribution in this arbitrary cross section is visualized and displayed on the
判定部144cは、タンディッシュ6の内側壁面近傍、すなわち、内側壁面を覆う耐火物との接触部分における溶鋼の温度が所定の温度範囲内であるか否かを判定する。例えば、判定部144cは、温度データ抽出部142cが抽出した鉛直断面における推定温度Tejの最大温度または最低温度が所定の温度範囲内であるか否かを判定する。所定の温度範囲は予め固定的に設定しておく構成としてもよいし、ユーザ操作に従って決定することとしてもよい。ユーザ操作に従って決定する場合には、入力部11を介してユーザによる温度範囲の入力操作を受け付け、判定部144cは、ユーザが入力した温度範囲に従って前述の判定を行う。
The determination unit 144c determines whether or not the temperature of the molten steel in the vicinity of the inner wall surface of the
温度制御部145cは、判定部144cが行った判定結果に従ってプラズマ加熱装置63,63による加熱温度を制御する。具体的には、温度制御部145cは、判定部144cにおいて温度範囲外と判定した場合には、その温度範囲外と判定した最大温度または最低温度が温度範囲内となるようにプラズマ加熱装置63,63の出力を制御する。
The
以上説明した構成の温度推定装置10cにおいて、各温度推定点jの温度推定を行った。具体的には、温度推定部141が、各温度推定点jの推定温度Tejを算出し、温度データ抽出部142cが、タンディッシュ6の例えば内側壁面近傍(壁面から50mm)の鉛直断面における推定温度Tejを抽出し、温度分布表示処理部143cが、前述の鉛直断面における推定温度Tejを等値線図化することで推定結果を得た。推定に用いた各温度実測部位iの既知温度Ti、すなわち、タンディッシュ6内の対応する設置位置P61〜P65に設置された熱電対64−1〜64−5の温度実測値を、プラズマ加熱装置63の加熱位置の既知温度Tiおよびノズル61の既知流入温度Tiとともに表5に示す。
In the
図39は、実施の形態5における温度推定結果、すなわち、タンディッシュ6の内側壁面近傍(壁面から50mm)の鉛直断面における温度分布を等値線図化した図である。図39に示すように、この実施の形態5によれば、タンディッシュ6内の溶鋼流動の影響を反映した温度推定が実現でき、溶鋼の流動効果を考慮して温度分布を高精度に推定できる。また、例えばタンディッシュ6の長辺方向の内側壁面近傍の鉛直断面の推定温度Tejを抽出し、この断面における推定温度Tejを例えば等値線図化して提示することができるので、ユーザは、内側壁面を覆う耐火物との接触部分における溶鋼の温度を容易に把握できる。また、この接触部分における溶鋼の温度が所定の温度範囲外の場合に、プラズマ加熱装置63,63の出力を制御して溶鋼の温度制御を行うことができるので、タンディッシュ6の内側壁面を覆う耐火物損傷を防止できる。
FIG. 39 is an isoline diagram of the temperature estimation result in the fifth embodiment, that is, the temperature distribution in the vertical section near the inner wall surface of the tundish 6 (50 mm from the wall surface). As shown in FIG. 39, according to the fifth embodiment, the temperature estimation reflecting the influence of the molten steel flow in the
上記した実施の形態では、本発明の適用対象として部屋、水槽、溶融亜鉛めっきポット、および連続鋳造用のタンディッシュを例示したが、これらに限定されるものではなく、本発明は、流体がかかわるものであれば幅広く適用することができる。例えば、鉄鋼プロセスでは、溶融金属保持炉や、連続鋳造鋳型、取鍋等における温度推定に適用が可能である。また、鉄鋼分野に限らず、化学プロセスや水処理設備等にも同様に適用が可能である。また、本発明は、単純な1次元流れの流体系のみならず、複雑な3次元流れとなる流体系まで幅広い流動状態の流体系に対して適応することが可能である。 In the above-described embodiment, the room, the water tank, the hot dip galvanizing pot, and the tundish for continuous casting are exemplified as the application target of the present invention. However, the present invention is not limited thereto, and the present invention involves a fluid. If it is a thing, it can apply widely. For example, the steel process can be applied to temperature estimation in a molten metal holding furnace, a continuous casting mold, a ladle or the like. Further, the present invention is not limited to the steel field, and can be similarly applied to chemical processes, water treatment facilities, and the like. Further, the present invention can be applied not only to a simple one-dimensional flow fluid system but also to a fluid system having a wide range of fluid states, including a fluid system having a complicated three-dimensional flow.
以上のように、本発明の流体系の温度推定方法、流体系の温度分布推定方法、流体系の温度分布モニタリング方法、および温度推定装置は、温度計測装置の配置に制約を与えることなく流体の流れによる熱輸送を考慮した高精度な温度推定を実現するのに適している。また、本発明の溶融亜鉛めっきポット内の溶融亜鉛温度制御方法および溶融亜鉛めっき鋼板によれば、表面欠陥がない溶融亜鉛めっき鋼板を提供することができる。また、本発明のタンディッシュ内の溶鋼温度制御方法によれば、タンディッシュの耐火物損傷を抑制することができる。 As described above, the fluid system temperature estimation method, the fluid system temperature distribution estimation method, the fluid system temperature distribution monitoring method, and the temperature estimation device of the present invention can be used without restricting the arrangement of the temperature measurement device. It is suitable for realizing highly accurate temperature estimation considering heat transport by flow. Moreover, according to the hot dip galvanization temperature control method and the hot dip galvanized steel sheet in the hot dip galvanizing pot of the present invention, a hot dip galvanized steel sheet free from surface defects can be provided. Moreover, according to the molten steel temperature control method in the tundish of this invention, the refractory damage of a tundish can be suppressed.
1,10,10a,10b,10c 温度推定装置
2 温度計測装置
11 入力部
12 表示部
13 記憶部
14,14a,14b,14c 制御部
141 温度推定部
142a,142b,142c 温度データ抽出部
143,143a,143b,143c 温度分布表示処理部
144b,144c 判定部
145b,145c 温度制御部
3 部屋
4 水槽
34−1〜34−4,45−1〜45−9 温度計
5 溶融亜鉛めっきポット
52 誘導加熱装置
6 タンディッシュ
63 プラズマ加熱装置
55−1〜55−8,64−1〜64−5 熱電対
1, 10, 10a, 10b, 10c
Claims (14)
前記温度既知領域の位置情報と流体系全域における流体の流れを表す流体系の流れ場に関する情報とを用いて、温度既知領域を通過した、または温度既知領域内で生成した流体のうち、他の温度既知領域を通過することなく前記温度推定点まで到達した流体の、温度推定点の全流体中に占める比率を温度推定点における温度既知領域の勢力として取得する勢力取得工程と、
各温度既知領域の温度と前記温度推定点における勢力とに関する情報を用いて、前記温度推定点における温度を推定する温度推定工程と、
を含むことを特徴とする流体系の温度推定方法。 A fluid system temperature estimation method for estimating a temperature at an arbitrary temperature estimation point of a fluid system having two or more known temperature regions,
Using the positional information of the temperature known region and the information on the flow field of the fluid system representing the flow of the fluid in the entire fluid system, other fluids that have passed through the temperature known region or generated in the temperature known region A force acquisition step of acquiring, as a force of the temperature known region at the temperature estimated point, a ratio of the fluid that has reached the temperature estimated point without passing through the temperature known region to the total fluid of the temperature estimated point;
A temperature estimation step of estimating the temperature at the temperature estimation point using information on the temperature of each temperature known region and the power at the temperature estimation point;
A method for estimating the temperature of a fluid system, comprising:
前記温度推定工程は、前記温度推定点における各温度既知領域の下流側勢力に関する情報を用いて前記温度推定点の温度を推定する工程を含む
ことを特徴とする請求項1に記載の流体系の温度推定方法。 According to the advection diffusion phenomenon by the flow field, the force acquisition step reaches the temperature estimation point without passing through another temperature known region among the fluids that have passed through the temperature known region or generated in the temperature known region. A downstream force acquisition step of acquiring the ratio of the estimated fluid in the total fluid at the temperature estimation point as the downstream force of the temperature known region at the temperature estimation point,
2. The fluid system according to claim 1, wherein the temperature estimation step includes a step of estimating a temperature of the temperature estimation point using information on a downstream power of each temperature known region at the temperature estimation point. Temperature estimation method.
前記温度推定工程は、前記温度推定点における各温度既知領域の上流側勢力に関する情報を用いて前記温度推定点の温度を推定する工程を含む
ことを特徴とする請求項1または2に記載の流体系の温度推定方法。 The force acquisition step includes, among the fluids that have passed through the temperature known region or generated in the temperature known region, according to the advection diffusion phenomenon caused by the reversal flow field of the fluid system that represents a flow opposite to the fluid flow. An upstream force acquisition step of acquiring, as the upstream force of the temperature known area at the temperature estimated point, the ratio of the temperature estimated point in the total fluid of the fluid that has reached the temperature estimated point without passing through another temperature known area; Including
3. The flow according to claim 1, wherein the temperature estimation step includes a step of estimating a temperature of the temperature estimation point using information on an upstream power of each temperature known region at the temperature estimation point. System temperature estimation method.
前記流体が前記温度既知領域と前記温度推定点との間で移動する際に要する伝達時間を取得する伝達時間取得工程と、
を含み、
前記温度推定工程は、温度推定を行う時点に対して前記伝達時間だけ過去、もしくは未来における時点を抽出時点とし、前記時系列温度データから前記抽出時点における温度既知領域の温度を抽出し、抽出された温度を用いて前記温度推定点の温度を推定する工程を含むことを特徴とする請求項1〜3のいずれか1つに記載の流体系の温度推定方法。 A time-series temperature data acquisition step for acquiring time-series temperature data including the temperature in a temperature known region where the temperature is known in time series;
A transmission time acquisition step of acquiring a transmission time required when the fluid moves between the temperature known region and the temperature estimation point;
Including
The temperature estimation step is performed by extracting a time point in the past or in the future for the temperature estimation time as an extraction time point and extracting the temperature of the temperature known region at the extraction time point from the time series temperature data. The method for estimating a temperature of a fluid system according to claim 1, further comprising: estimating a temperature of the temperature estimation point using a measured temperature.
前記温度推定工程は、各温度既知領域の前記重みを用いた重み付き平均処理を行って、前記温度推定点の温度を推定する工程を含むことを特徴とする請求項1〜4のいずれか1つに記載の流体系の温度推定方法。 Using information about the power of each temperature known area, including a weight calculating step for calculating the weight of each temperature known area,
The temperature estimation step includes a step of estimating a temperature of the temperature estimation point by performing a weighted average process using the weight of each temperature known region. The temperature estimation method of the fluid system as described in one.
請求項1〜8のいずれか1つに記載の流体系の温度推定方法を用いて前記流体系の全域に設定した温度推定点の温度を推定し、
前記各温度推定点について推定した温度を前記流体系の温度分布として推定することを特徴とする流体系の温度分布推定方法。 A temperature distribution estimation method for a fluid system having a temperature distribution,
Estimating the temperature of the temperature estimation point set in the entire area of the fluid system using the fluid system temperature estimation method according to any one of claims 1 to 8,
A temperature distribution estimation method for a fluid system, wherein the temperature estimated for each temperature estimation point is estimated as a temperature distribution of the fluid system.
請求項9に記載の流体系の温度分布推定方法を用いて推定した前記流体系の温度分布をもとに、前記流体系の任意の断面における温度分布を可視化して画面表示することを特徴とする流体系の温度分布モニタリング方法。 A temperature distribution monitoring method for a fluid system having a temperature distribution,
A temperature distribution in an arbitrary cross section of the fluid system is visualized and displayed on the screen based on the temperature distribution of the fluid system estimated using the temperature distribution estimation method of the fluid system according to claim 9. To monitor temperature distribution of fluid system.
前記温度既知領域の位置情報と流体系全域における流体の流れを表す流体系の流れ場に関する情報とを用いて、温度既知領域を通過した、または温度既知領域内で生成した流体のうち、他の温度既知領域を通過することなく前記温度推定点まで到達した流体の、前記温度推定点の全流体中に占める比率を前記温度推定点における温度既知領域の勢力として取得する勢力取得手段と、
各温度既知領域の温度と前記温度推定点における勢力とに関する情報を用いて、前記温度推定点における温度を推定する温度推定手段と、
を備えることを特徴とする温度推定装置。 A temperature estimation device for estimating a temperature at an arbitrary temperature estimation point of a fluid system having two or more known temperature regions,
Using the positional information of the temperature known region and the information on the flow field of the fluid system representing the flow of the fluid in the entire fluid system, other fluids that have passed through the temperature known region or generated in the temperature known region A force acquisition means for acquiring, as a force of the temperature known region at the temperature estimated point, a ratio of the fluid that has reached the temperature estimated point without passing through the temperature known region to the total fluid of the temperature estimated point;
Temperature estimation means for estimating the temperature at the temperature estimation point using information on the temperature of each temperature known region and the power at the temperature estimation point;
A temperature estimation device comprising:
請求項7に記載の流体系の温度推定方法により推定した前記溶融亜鉛めっきポット内の溶融亜鉛温度データから、前記溶融亜鉛めっきポット内の所定の領域における溶融亜鉛の温度を抽出する温度抽出ステップと、
抽出した温度が、所定の閾値範囲内にあるか否かを判定する判定ステップと、
前記判定ステップにおいて、前記抽出した温度が閾値範囲外と判定された場合、前記抽出した温度が閾値範囲内となるよう前記溶融亜鉛めっきポットの加熱手段の出力を操作する制御ステップと、
を含むことを特徴とする溶融亜鉛めっきポット内の溶融亜鉛温度制御方法。 A hot dip galvanizing temperature control method in a hot dip galvanizing pot,
A temperature extraction step of extracting a temperature of molten zinc in a predetermined region in the hot dip galvanizing pot from hot dip galvanizing temperature data in the hot dip galvanizing pot estimated by the fluid system temperature estimating method according to claim 7; ,
A determination step of determining whether the extracted temperature is within a predetermined threshold range;
In the determination step, when it is determined that the extracted temperature is outside the threshold range, a control step of operating the output of the heating means of the hot dip galvanizing pot so that the extracted temperature is within the threshold range;
A method for controlling the temperature of a hot dip galvanized hot dip galvanizing pot.
請求項8に記載の流体系の温度推定方法により推定した前記タンディッシュ内の溶鋼温度データから、前記タンディッシュ内の所定の領域における溶鋼の温度を抽出する温度抽出ステップと、
抽出した温度が、所定の閾値範囲内にあるか否かを判定する判定ステップと、
前記判定ステップにおいて、前記抽出した温度が閾値範囲外と判定された場合、前記抽出した温度が閾値範囲内となるよう前記タンディッシュの加熱手段の出力を操作する制御ステップと、
を含むことを特徴とするタンディッシュ内の溶鋼温度制御方法。 A method for controlling the temperature of molten steel in a tundish,
A temperature extraction step for extracting the temperature of the molten steel in a predetermined region in the tundish from the molten steel temperature data in the tundish estimated by the fluid system temperature estimation method according to claim 8;
A determination step of determining whether the extracted temperature is within a predetermined threshold range;
In the determination step, when the extracted temperature is determined to be outside the threshold range, a control step of operating the output of the tundish heating means so that the extracted temperature is within the threshold range;
A method for controlling the temperature of molten steel in a tundish, comprising:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011210259A JP4984000B1 (en) | 2010-09-30 | 2011-09-27 | Fluid system temperature estimation method, fluid system temperature distribution estimation method, fluid system temperature distribution monitoring method, temperature estimation device, hot dip galvanized steel temperature control method in hot dip galvanizing pot, hot dip galvanized steel sheet, and molten steel in tundish Temperature control method |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010222844 | 2010-09-30 | ||
JP2010222844 | 2010-09-30 | ||
JP2011092434 | 2011-04-18 | ||
JP2011092434 | 2011-04-18 | ||
JP2011210259A JP4984000B1 (en) | 2010-09-30 | 2011-09-27 | Fluid system temperature estimation method, fluid system temperature distribution estimation method, fluid system temperature distribution monitoring method, temperature estimation device, hot dip galvanized steel temperature control method in hot dip galvanizing pot, hot dip galvanized steel sheet, and molten steel in tundish Temperature control method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP4984000B1 JP4984000B1 (en) | 2012-07-25 |
JP2012233869A true JP2012233869A (en) | 2012-11-29 |
Family
ID=45893063
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011210259A Active JP4984000B1 (en) | 2010-09-30 | 2011-09-27 | Fluid system temperature estimation method, fluid system temperature distribution estimation method, fluid system temperature distribution monitoring method, temperature estimation device, hot dip galvanized steel temperature control method in hot dip galvanizing pot, hot dip galvanized steel sheet, and molten steel in tundish Temperature control method |
Country Status (4)
Country | Link |
---|---|
JP (1) | JP4984000B1 (en) |
KR (1) | KR101318209B1 (en) |
CN (1) | CN103124899B (en) |
WO (1) | WO2012043624A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015078916A (en) * | 2013-10-17 | 2015-04-23 | 日本電信電話株式会社 | Air current estimation method and device therefor |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5935837B2 (en) * | 2014-07-07 | 2016-06-15 | Jfeスチール株式会社 | Flow state estimation method and flow state estimation apparatus for molten steel |
JP6163590B1 (en) * | 2016-05-31 | 2017-07-12 | 株式会社日阪製作所 | Simulation method, simulation program, and simulation apparatus including storage medium incorporating the program |
CN110806235B8 (en) * | 2019-11-15 | 2021-12-17 | 北京科技大学 | Indoor environment monitoring method, device, equipment and storage medium |
CN116878693B (en) * | 2023-09-05 | 2023-11-21 | 济宁市金桥煤矿 | Mining electromechanical equipment monitoring and management method and system |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01307696A (en) * | 1988-06-06 | 1989-12-12 | Hitachi Ltd | In-pile temperature monitoring system of fast reactor |
JP2919532B2 (en) * | 1990-03-02 | 1999-07-12 | 三菱化学株式会社 | Analysis condition setting method for thermal fluid analysis in furnace |
JP3033466B2 (en) * | 1995-04-28 | 2000-04-17 | 住友金属工業株式会社 | Blast furnace operation method |
JP3952329B2 (en) * | 1997-08-27 | 2007-08-01 | 日立金属株式会社 | Airflow simulation method in a factory building |
JP2000283526A (en) * | 1999-03-25 | 2000-10-13 | Internatl Business Mach Corp <Ibm> | Air-conditioning system and method therefor |
JP4693731B2 (en) * | 2006-09-21 | 2011-06-01 | 株式会社東芝 | Air conditioning sensor system |
JP2008241139A (en) | 2007-03-27 | 2008-10-09 | Sanyo Electric Co Ltd | Sterilizing device |
JP2008248353A (en) * | 2007-03-30 | 2008-10-16 | Jfe Steel Kk | Method for producing hot dip metal-plated metal sheet, and production equipment therefor |
KR100868457B1 (en) | 2007-05-31 | 2008-11-11 | 주식회사 포스코 | Galvannealed steel sheet having superior adhesiveness of plated film and method for manufacturing the same |
JP2009241139A (en) * | 2008-03-31 | 2009-10-22 | Kobe Steel Ltd | Forecasting method for molten steel temperature within tundish, and management method |
CN101294946B (en) * | 2008-06-05 | 2011-06-15 | 武汉科技大学 | Model modeling and optimization method for mixed on-line estimating killing oxygen by adding aluminum |
JP4968388B2 (en) * | 2010-02-05 | 2012-07-04 | Jfeスチール株式会社 | Fluid system temperature estimation method, Fluid system temperature concentration estimation method, Fluid system temperature distribution estimation method, Fluid system temperature monitoring method, Molten metal equipment molten metal temperature control method, Molten metal equipment concentration and Temperature estimation method, hot dip galvanized steel sheet and fluid temperature estimation device |
-
2011
- 2011-09-27 JP JP2011210259A patent/JP4984000B1/en active Active
- 2011-09-28 KR KR1020137008653A patent/KR101318209B1/en active IP Right Grant
- 2011-09-28 CN CN201180046338.5A patent/CN103124899B/en active Active
- 2011-09-28 WO PCT/JP2011/072176 patent/WO2012043624A1/en active Application Filing
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015078916A (en) * | 2013-10-17 | 2015-04-23 | 日本電信電話株式会社 | Air current estimation method and device therefor |
Also Published As
Publication number | Publication date |
---|---|
KR101318209B1 (en) | 2013-10-15 |
CN103124899B (en) | 2014-08-13 |
CN103124899A (en) | 2013-05-29 |
KR20130042055A (en) | 2013-04-25 |
JP4984000B1 (en) | 2012-07-25 |
WO2012043624A1 (en) | 2012-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4984000B1 (en) | Fluid system temperature estimation method, fluid system temperature distribution estimation method, fluid system temperature distribution monitoring method, temperature estimation device, hot dip galvanized steel temperature control method in hot dip galvanizing pot, hot dip galvanized steel sheet, and molten steel in tundish Temperature control method | |
Lee et al. | An evaluation of empirically-based models for predicting energy performance of vapor-compression water chillers | |
Chakraborty et al. | Estimation of surface heat flux in continuous casting mould with limited measurement of temperature | |
Kang et al. | Approximate methods for uncertainty analysis of water distribution systems | |
CN109929955A (en) | A kind of detection method of blast furnace crucible corrosion situation | |
Britsch et al. | Natural circulation FLiBe loop overview | |
Kim et al. | Experimental and numerical analysis of heat transfer phenomena in a sensor tube of a mass flow controller | |
Melissari et al. | Development of a heat transfer dimensionless correlation for spheres immersed in a wide range of Prandtl number fluids | |
Kychkin | Synthesizing a system for remote energy monitoring in manufacturing | |
Zagaria et al. | Monitoring erosion and skull profile in blast furnace hearth | |
CN113272748A (en) | Analysis system and analysis method | |
JP4695376B2 (en) | Heating or cooling characteristic evaluation method and apparatus, reaction vessel operation management method and apparatus, computer program, and computer-readable recording medium | |
JP4973802B2 (en) | Method for estimating concentration and temperature of substance component in fluid system, temperature control of molten zinc in hot dip galvanizing pot, and control method of aluminum concentration in molten zinc | |
JP4968388B2 (en) | Fluid system temperature estimation method, Fluid system temperature concentration estimation method, Fluid system temperature distribution estimation method, Fluid system temperature monitoring method, Molten metal equipment molten metal temperature control method, Molten metal equipment concentration and Temperature estimation method, hot dip galvanized steel sheet and fluid temperature estimation device | |
JP5765035B2 (en) | Concentration estimation method for fluid constituent material in fluid system, concentration distribution estimation method for fluid constituent material in fluid system, concentration distribution monitoring method for fluid constituent material in fluid system, and concentration estimation device | |
JP2005315581A (en) | Thermal fatigue evaluation system for piping | |
JP4833621B2 (en) | Method, apparatus, computer program, and computer-readable recording medium for estimating temperature or heat flux of reaction vessel | |
CN103400037A (en) | Method for determining welding temperature of straight seam welded pipe | |
Umbricht et al. | Optimal estimation of thermal diffusivity in an energy transfer problem | |
JP5923854B2 (en) | Method for estimating concentration of substance component in fluid system, method for estimating concentration distribution, method for monitoring concentration, method for managing aluminum concentration in hot dip galvanized pot accommodated in hot dip galvanizing pot, and apparatus for estimating concentration of substance component in fluid system | |
Han et al. | Sensitivity analysis about transient three-dimensional IHCP with multi-parameters in an elbow pipe with thermal stratification | |
Prakash et al. | Estimating Non-Linear Heat Flux in Continuous Billet Casting Process With Limited Sensors | |
Gómez et al. | Modeling and optimization of triple tube heat exchangers. Theoretical formulation, CFD model and experimental contrast | |
Ward | Gravity stratified mixed convection in liquid metal pools: implications and experimental interpretation for LMFRs | |
CN116306144A (en) | Thermal fluid-solid coupling simulation modeling and experimental method for 1045 steel quenching process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120409 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4984000 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150511 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |