JP2003028819A - Method for calculating heat transfer by piston effect at supercritical fluid - Google Patents
Method for calculating heat transfer by piston effect at supercritical fluidInfo
- Publication number
- JP2003028819A JP2003028819A JP2001217306A JP2001217306A JP2003028819A JP 2003028819 A JP2003028819 A JP 2003028819A JP 2001217306 A JP2001217306 A JP 2001217306A JP 2001217306 A JP2001217306 A JP 2001217306A JP 2003028819 A JP2003028819 A JP 2003028819A
- Authority
- JP
- Japan
- Prior art keywords
- temperature
- physical property
- calculating
- calculated
- heat transfer
- 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.)
- Pending
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P20/00—Technologies relating to chemical industry
- Y02P20/50—Improvements relating to the production of bulk chemicals
- Y02P20/54—Improvements relating to the production of bulk chemicals using solvents, e.g. supercritical solvents or ionic liquids
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Complex Calculations (AREA)
Abstract
Description
【0001】[0001]
【発明の属する技術分野】本発明は、超臨界流体におけ
るピストン効果による熱移動の算出方法に関する。TECHNICAL FIELD The present invention relates to a method for calculating heat transfer due to a piston effect in a supercritical fluid.
【0002】[0002]
【従来の技術】近年、超臨界流体の多方面での利用が期
待されている。超臨界流体は、粘度、拡散速度等の性質
が気体と液体の中間に当たる状態にあるため、気体のよ
うにどこにでも拡散する性質と液体のように成分を溶解
する性質とを併せ持ち、例えば、食品や医療品等の分野
では、有機媒体の代替物として利用することができ、環
境に優しい技術として注目されている。ところで、超臨
界流体を工業的に利用するためには、超臨界流体におけ
る熱に関する、例えば、超臨界流体内での熱移動を正確
に把握することが不可欠とされている。2. Description of the Related Art In recent years, the use of supercritical fluids in various fields is expected. Since supercritical fluids have properties such as viscosity and diffusion speed that are in the middle of gas and liquid, they have the property of diffusing anywhere like gas and the property of dissolving components like liquid. It can be used as a substitute for an organic medium in the fields of medicine, medical products, and the like, and is attracting attention as an environment-friendly technology. By the way, in order to industrially use the supercritical fluid, it is indispensable to accurately grasp the heat transfer in the supercritical fluid, for example, the heat transfer in the supercritical fluid.
【0003】従来、超臨界流体中では、熱拡散率は微小
であるため、自然対流が発生しない場合には、熱移動は
ほとんど無いものと考えられてきた。ところが、ドイツ
のNitscheとStraub等の1986年の報告
によれば、宇宙実験において、対流の発生しない微小重
力下において超臨界流体中で高速な熱移動現象が見られ
ることが明らかになった。この超臨界流体における熱移
動現象は、フランスのZappoli等によって、19
90年にピストン効果と名付けられた。Conventionally, since the thermal diffusivity is very small in a supercritical fluid, it has been considered that there is almost no heat transfer when natural convection does not occur. However, according to a report by Nitsche and Straub et al. In Germany in 1986, it was revealed in space experiments that a high-speed heat transfer phenomenon was observed in a supercritical fluid under microgravity in which convection did not occur. The heat transfer phenomenon in this supercritical fluid was reported by Zappoli et al.
It was named the piston effect in 1990.
【0004】ピストン効果とは、上述のごとく、超臨界
流体の臨界点近傍における高速な熱移動現象を指し、従
来の熱移動現象である、熱伝導、対流、輻射とは異なる
伝熱機構とも言われているが、理論的には未だ不明な点
も多い。このピストン効果は、前出のNitsche、
Straub、および日本の小貫等によってさらに研究
が進められ、特に、小貫等はピストン効果に関係ある物
性値を一定とした温度分布に関する基礎式を基に熱移動
現象を解析している。As described above, the piston effect refers to a high-speed heat transfer phenomenon in the vicinity of the critical point of a supercritical fluid, and is also called a heat transfer mechanism different from the conventional heat transfer phenomena of heat conduction, convection, and radiation. However, there are many unclear points in theory. This piston effect is due to Nitsche,
Further research has been conducted by Straub and Konuki of Japan. In particular, Konuki et al. Analyze the heat transfer phenomenon based on a basic equation relating to temperature distribution with constant physical properties related to the piston effect.
【0005】[0005]
【発明が解決しようとする課題】しかし、上記従来の超
臨界流体のピストン効果による熱移動の解析は、未だ研
究段階にとどまるものであり、工業化レベルでの利用が
有効な算出方法を提供するものではない。また、上記従
来の研究では、ピストン効果による熱移動の算出に必要
な物性値を一定値として計算しているために、実際の超
臨界流体の熱移動現象に即した精度良い算出結果が得ら
れないという問題もあった。However, the conventional analysis of heat transfer due to the piston effect of the supercritical fluid is still in the research stage, and it provides a calculation method that is effectively used at the industrial level. is not. Further, in the above-mentioned conventional research, since the physical property value necessary for calculating the heat transfer due to the piston effect is calculated as a constant value, an accurate calculation result according to the actual heat transfer phenomenon of the supercritical fluid can be obtained. There was also the problem of not having it.
【0006】本発明の目的は、上記の従来の問題点に鑑
み、温度、圧力による物性値変化を考慮すると共に、コ
ンピュータ等の情報処理装置を用いて高精度でかつ容易
に算出することを可能にした超臨界流体におけるピスト
ン効果による熱移動の算出方法を提供することにある。In view of the above-mentioned conventional problems, an object of the present invention is to consider changes in physical properties due to temperature and pressure, and to calculate with high accuracy and easily using an information processing device such as a computer. Another object of the present invention is to provide a method for calculating heat transfer due to the piston effect in a supercritical fluid described above.
【0007】[0007]
【課題を解決するための手段】本発明は、上記の課題を
解決するために、次のような手段を採用した。The present invention employs the following means in order to solve the above problems.
【0008】第1の手段は、超臨界流体におけるピスト
ン効果による熱移動の算出方法において、算出に必要な
解析領域、解析時間、初期条件、および境界条件の各情
報を入力する第1のステップと、初期温度T0および初
期圧力p0に基づいて物性値を取得する第2のステップ
と、前記第2のステップまたは後述する第6のステップ
で取得された物性値を用いて、∂T/∂t=∇(λ
(T、p)・∇T)/ρ(T、p)・cp(T、p)式
に基づいて熱伝導による温度Tを算出する第3のステッ
プと、前記第2のステップまたは後述する第6のステッ
プで取得された物性値を用いて、∂p/∂t=∫Vρ
(T、p)・αp(T、p)・[∂T/∂t]・dV/
∫Vρ(T、p)・χT(T、p)・dV式に基づいて
圧力pを算出する第4のステップと、前記第2のステッ
プまたは後述する第6のステップで取得された物性値を
用いて、Q=[1−cV(T、p)/cp(T、p)]
・[∂T/∂p]ρ(T、p)・∂p/∂t式に基づい
て発熱量Qを算出する第5のステップと、前記第3のス
テップで算出された温度Tおよび前記第4のステップで
算出された圧力pに基づいて物性値を取得する第6のス
テップと、第5のステップで算出された発熱量Qおよび
第6のステップで取得された物性値を用いて、∂T/∂
t=∇(λ(T、p)・∇T)/ρ(T、p)・c
p(T、p)+[1−c V(T、p)/cp(T、
p)]・[∂T/∂p]ρ(T、p)・∂p/∂t式に
基づいてピストン効果を考慮した温度Tを算出する第7
のステップとからなることを特徴とする。ただし、上記
の各物性値は、熱伝導率λ(T、p)、定圧比熱c
p(T、p)、密度ρ(T、p)、定積比熱cV(T、
p)、体膨張率αp(T、p)、等温圧縮率χT(T、
p)、[∂T/∂p]ρ(T、p)。The first means is a fix in a supercritical fluid.
In the method of calculating heat transfer due to
Analysis area, analysis time, initial condition, and boundary condition information
First step of inputting information and initial temperature T0And first
Period pressure p0Second step of obtaining physical property values based on
And the second step or the sixth step described later
∂T / ∂t = ∇ (λ
(T, p) ・ ∇T) / ρ (T, p) ・ cp(T, p) formula
The third step for calculating the temperature T due to heat conduction based on
And a second step or a sixth step described later.
∂p / ∂t = ∫ using the physical property values acquired inVρ
(T, p) / αp(T, p) ・ [∂T / ∂t] ・ dV /
∫Vρ (T, p) · χTBased on (T, p) · dV formula
The fourth step of calculating the pressure p and the second step
Or the physical property values obtained in the sixth step described below
Using Q = [1-cV(T, p) / cp(T, p)]
・ [∂T / ∂p]ρBased on (T, p) ・ ∂p / ∂t formula
The fifth step of calculating the heat generation amount Q by
In the temperature T calculated in step and the fourth step
A sixth step for obtaining a physical property value based on the calculated pressure p
Step and the calorific value Q calculated in the fifth step and
∂T / ∂ using the physical property values obtained in the sixth step
t = ∇ (λ (T, p) · ∇T) / ρ (T, p) · c
p(T, p) + [1-c V(T, p) / cp(T,
p)] ・ [∂T / ∂p]ρ(T, p) · ∂p / ∂t formula
7th which calculates the temperature T which considered the piston effect based on
And the steps of. However, above
The physical properties of are the thermal conductivity λ (T, p) and the constant pressure specific heat c
p(T, p), density ρ (T, p), constant volume specific heat cV(T,
p), body expansion rate αp(T, p), isothermal compressibility χT(T,
p), [∂T / ∂p]ρ(T, p).
【0009】第2の手段は、第1の手段において、前記
第3のステップおよび前記第7のステップで算出された
温度Tについて、経時的に算出された最新の温度T
n+1と経時的に前記最新の温度Tn+1の直前に算出
された温度Tnとに基づいて両者間の収束を判定し、収
束していないと判定されたときは前記第4のステップか
らの処理を繰り返させ、収束していると判定されたとき
は後述する第9のステップの処理を行わせる第8のステ
ップと、各ステップの算出結果を記憶ないし表示する第
9のステップとからなることを特徴とする。A second means is the latest temperature T calculated with time with respect to the temperature T calculated in the first step and the third step and the seventh step.
n + 1 and based on the time to the latest temperature T n + 1 of the temperature T n calculated immediately before to determine the convergence between the two, when it is determined not converged processing from the fourth step Is repeated, and when it is determined that the convergence has occurred, an eighth step for performing the processing of the ninth step described later and a ninth step for storing or displaying the calculation result of each step. Characterize.
【0010】第3の手段は、第2の手段において、前記
第9のステップの処理後、所定のタイムステップ数を処
理したかを判定し、処理していないと判定されたときは
前記第3のステップからの処理を行わせる第10のステ
ップとからなることを特徴とする。The third means determines, in the second means, whether or not a predetermined number of time steps have been processed after the processing of the ninth step, and if it is determined that the predetermined number of time steps have not been processed, the third means. And a tenth step for performing the processing from step.
【0011】第4の手段は、算出に必要な解析領域、解
析時間、初期条件、および境界条件の各情報を入力する
第1のステップと、初期温度T0および初期圧力p0に
基づいて物性値を取得する第2のステップと、前記第2
のステップで取得された物性値を用いて、∂T/∂t=
∇(λ(T、p)・∇T)/ρ(T、p)・cp(T、
p)式に基づいて熱伝導による温度Tを算出する第3の
ステップと、前記第2のステップまたは後述する第6の
ステップで取得された物性値を用いて、∂p/∂t=∫
Vρ(T、p)・αp(T、p)・[∂T/∂t]・d
V/∫Vρ(T、p)・χT(T、p)・dV式に基づ
いて圧力pを算出する第4のステップと、前記第2のス
テップまたは後述する第6のステップで取得された物性
値を用いて、Q=[1−cV(T、p)/cp(T、
p)]・[∂T/∂p]ρ(T、p)・∂p/∂t式に
基づいて発熱量Qを算出する第5のステップと、前記第
3のステップで算出された温度Tおよび前記第4のステ
ップで算出された圧力pに基づいて物性値を取得する第
6のステップと、第5のステップで算出された発熱量Q
および第6のステップで取得された物性値を用いて、∂
T/∂t=∇(λ(T、p)・∇T)/ρ(T、p)・
cp(T、p)+[1−cV(T、p)/cp(T、
p)]・[∂T/∂p]ρ(T、p)・∂p/∂t式に
基づいてピストン効果を考慮した温度Tを算出する第7
のステップと、前記第7のステップで算出された温度T
について、経時的に算出された最新の温度Tn+1と経
時的に前記最新の温度Tn+1の直前に算出された温度
Tnとに基づいて両者間の収束を判定し、収束していな
いと判定されたときは前記第4のステップからの処理を
繰り返させ、収束していると判定されたときは後述する
第9のステップの処理を行わせる第8のステップと、各
ステップの算出結果を記憶ないし表示する第9のステッ
プと、前記第9のステップの処理後、所定のタイムステ
ップ数を処理したかを判定し、処理していないと判定さ
れたときは前記第3のステップからの処理を行わせる第
10のステップと、からなることを特徴とする超臨界流
体におけるピストン効果による熱移動を算出するプログ
ラム。ただし、上記の各物性値は、熱伝導率λ(T、
p)、定圧比熱cp(T、p)、密度ρ(T、p)、定
積比熱cV(T、p)、体膨張率αp(T、p)、等温
圧縮率χT(T、p)、[∂T/∂p]ρ(T、p)。A fourth means is to input the information of the analysis region, the analysis time, the initial condition and the boundary condition necessary for the calculation, and the physical properties based on the initial temperature T 0 and the initial pressure p 0. A second step of obtaining a value, and the second step
∂T / ∂t = using the physical property values obtained in step
∇ (λ (T, p) · ∇T) / ρ (T, p) · c p (T,
p) ∂p / ∂t = ∫ using the third step of calculating the temperature T due to heat conduction based on the equation and the physical property values obtained in the second step or the sixth step described later.
V ρ (T, p) ・ α p (T, p) ・ [∂T / ∂t] ・ d
V / ∫ V ρ (T, p) · χ T (T, p) · dV obtained in the fourth step of calculating the pressure p based on the equation and the second step or the sixth step described later. by using the physical property values, Q = [1-c V (T, p) / c p (T,
p)] · [∂T / ∂p] ρ (T, p) · ∂p / ∂t based on the fifth step of calculating the calorific value Q, and the temperature T calculated in the third step. And a sixth step of obtaining a physical property value based on the pressure p calculated in the fourth step, and a calorific value Q calculated in the fifth step.
And using the physical property values obtained in the sixth step, ∂
T / ∂t = ∇ (λ (T, p) ・ ∇T) / ρ (T, p) ・
c p (T, p) + [1-c V (T, p) / c p (T,
p)] · [∂T / ∂p] ρ (T, p) · ∂p / ∂t based on the formula, the temperature T considering the piston effect is calculated.
And the temperature T calculated in the seventh step
For determines convergence between them on the basis of the temperature T n calculated immediately before the time the calculated latest temperature T n + 1 and over time the latest temperature T n + 1, and does not converge determination When it is determined that the process from the fourth step is repeated, and when it is determined that the process is converged, an eighth step for performing the process of a ninth step described later and a calculation result of each step are stored. Or the ninth step of displaying and after the processing of the ninth step, it is determined whether or not a predetermined number of time steps have been processed, and when it is determined that the predetermined time step has not been processed, the processing from the third step is performed. A program for calculating heat transfer due to a piston effect in a supercritical fluid, characterized by comprising: a tenth step to be performed. However, the above-mentioned physical property values are the thermal conductivity λ (T,
p), pressure specific heat c p (T, p), the density ρ (T, p), specific heat at constant volume c V (T, p), volume expansion ratio α p (T, p), isothermal compressibility chi T (T , P), [∂T / ∂p] ρ (T, p).
【0012】第5の手段は、コンピュータ読み取り可能
な記録媒体として、第4の手段に記載の超臨界流体にお
けるピストン効果による熱移動を算出するプログラムを
記録したことを特徴とする。A fifth means is characterized in that, as a computer-readable recording medium, a program for calculating the heat transfer due to the piston effect in the supercritical fluid described in the fourth means is recorded.
【0013】[0013]
【発明の実施の形態】従来、Straub等によって提
案され、超臨界流体におけるピストン効果による熱移動
現象を解析するための計算式としては下式のものが知ら
れている。DESCRIPTION OF THE PREFERRED EMBODIMENTS Conventionally, the following formula has been known as a calculation formula proposed by Straub et al. For analyzing a heat transfer phenomenon due to a piston effect in a supercritical fluid.
【0014】
∂T/∂t=∇(λ・∇T)/ρcp+[1−cV/cp]・
[∂T/∂p]ρ・∂p/∂t・・・ (1)’
∂p/∂t=∫Vραp[∂T/∂t]dV/∫VρχTdV・・・
・・(2)’
ここで、上記各式において、λは熱伝導率、cpは定圧
比熱、ρは密度、cVは定積比熱(比熱比γ)、αは体
膨張率、cpは定圧比熱、χTは等温圧縮率からなる物
性値である。[0014] ∂T / ∂t = ∇ (λ · ∇T) / ρc p + [1-c V / c p] · [∂T / ∂p] ρ · ∂p / ∂t ··· (1) '∂p / ∂t = ∫ V ρα p [∂T / ∂t] dV / ∫ V ρχ T dV ... (2)' where λ is thermal conductivity and c p the specific heat at constant pressure, [rho is the density, c V is specific heat at constant volume (specific heat ratio gamma), alpha is the body expansion, the c p the specific heat at constant pressure, the chi T is a physical property value consisting of isothermal compressibility.
【0015】従来においても、これらの偏微分方程式を
適当な計算領域、境界条件、初期条件のもとに解くこと
により、超臨界流体におけるピストン効果による熱移動
現象を解明することが可能とされている。しかし、これ
らの偏微分方程式を解析的に解くのは一般的に困難なこ
とであり、また、本件発明者の知見によれば、超臨界流
体では温度、圧力の変化、特に、臨界点近傍では温度、
圧力のわずかな変化でも、これらの物性値に大きな変化
を与えることが分かり、各物性値が適正な値に設定され
なければ、実際の現象に即したピストン効果による熱移
動を算出することができない。Conventionally, it has been possible to solve the heat transfer phenomenon due to the piston effect in a supercritical fluid by solving these partial differential equations under appropriate calculation regions, boundary conditions and initial conditions. There is. However, it is generally difficult to solve these partial differential equations analytically, and according to the knowledge of the inventor of the present invention, in supercritical fluids, changes in temperature and pressure, especially in the vicinity of the critical point, temperature,
It has been found that even a slight change in pressure causes a large change in these physical property values, and unless each physical property value is set to an appropriate value, it is not possible to calculate the heat transfer due to the piston effect that matches the actual phenomenon. .
【0016】そこで、本件発明においては、温度、圧力
の変化に伴って変化する物性値の変化を有効な手段によ
って逐次算出過程に反映させ、かつコンピュータ処理に
よって高精度でかつ容易に算出することを可能にした、
超臨界流体におけるピストン効果による熱移動の算出方
法を提案するものである。以下に、本発明の一実施形態
を図1ないし図6を用いて説明する。本発明では、上述
のごとく、実際の超臨界流体におけるピストン効果によ
る熱移動現象においては、各物性値が温度T、圧力pに
よって大きな影響を受けることから、各物性値を温度T
および圧力pの関数として表す。従って、各物性値であ
る、熱伝導率λはλ(T、p)、定圧比熱cpはc
p(T、p)、密度ρはρ(T、p)、比熱比γはγ
(T、p)(定積比熱cVはc V(T、p))、体膨張
率αはαp(T、p)、等温圧縮率χTはχT(T、
p)、そして、[∂T/∂p]ρは[∂T/∂p]
ρ(T、p)として表わされ、上式(1)’式および
(2)’は下式のように表される。Therefore, in the present invention, temperature, pressure
The change of physical property value that changes with the change of
Is reflected in the sequential calculation process, and computer processing is performed.
Therefore, it is possible to calculate with high accuracy and easily.
Calculation method of heat transfer due to piston effect in supercritical fluid
It proposes a law. Below, one embodiment of the present invention
Will be described with reference to FIGS. In the present invention, the above
As a result, due to the piston effect in an actual supercritical fluid,
In the heat transfer phenomenon, each physical property value changes to temperature T and pressure p.
Therefore, the physical properties are greatly affected.
And as a function of pressure p. Therefore, in each physical property value
The thermal conductivity λ is λ (T, p) and the constant pressure specific heat cpIs c
p(T, p), density ρ is ρ (T, p), specific heat ratio γ is γ
(T, p) (constant volume specific heat cVIs c V(T, p)), body expansion
The rate α is αp(T, p), isothermal compressibility χTIs χT(T,
p), and [∂T / ∂p]ρIs [∂T / ∂p]
ρIs expressed as (T, p) and is expressed by the above equation (1) ′ and
(2) 'is expressed by the following equation.
【0017】 ∂T/∂t=∇(λ(T、p)・∇T)/ρ(T、p)・cp(T、p)+ [1−cV(T、p)/cp(T、p)]・ [∂T/∂p]ρ(T、p)・∂p/∂t・・・(1) ∂p/∂t=∫Vρ(T、p)・αp(T、p)[∂T/∂t]dV/ ∫Vρ(T、p)・χT(T、p)・dV・・・(2)。[0017] ∂T / ∂t = ∇ (λ ( T, p) · ∇T) / ρ (T, p) · c p (T, p) + [1-c V (T, p) / c p (T, p)] ・ [∂T / ∂p] ρ (T, p) ・ ∂p / ∂t ... (1) ∂p / ∂t = ∫ V ρ (T, p) ・ α p ( T, p) [∂T / ∂t ] dV / ∫ V ρ (T, p) · χ T (T, p) · dV ··· (2).
【0018】以下に、上記の(1)式および(2)式に
基づいて、本件発明に係る超臨界流体におけるピストン
効果による熱移動の算出方法を図1を用いて説明する。A method for calculating heat transfer due to the piston effect in the supercritical fluid according to the present invention will be described below with reference to the above equations (1) and (2) with reference to FIG.
【0019】図1は、本件発明に係る超臨界流体におけ
るピストン効果による熱移動の算出方法を示すフローチ
ャートである。FIG. 1 is a flow chart showing a method for calculating heat transfer due to the piston effect in a supercritical fluid according to the present invention.
【0020】まず、ステップ1において、算出に必要な
解析領域、解析時間、初期条件、境界条件等の各種の情
報を設定し入力する。本実施形態では、超臨界流体とし
ては水を用いた。なお、一般に用いられる水は温度64
7.096K、圧力22.064MPaで臨界点に達し
超臨界流体となる。算出には、図2に示すように、1次
元の物理モデルを採用した。First, in step 1, various kinds of information such as analysis area, analysis time, initial condition and boundary condition necessary for calculation are set and input. In this embodiment, water is used as the supercritical fluid. In addition, generally used water has a temperature of 64.
At 7.096 K and a pressure of 22.064 MPa, it reaches a critical point and becomes a supercritical fluid. For the calculation, as shown in FIG. 2, a one-dimensional physical model was adopted.
【0021】計算領域は、同図に示すように、1次元で
長さは0.01mであり、計算領域は153個のコント
ロールボリューム(なお、境界における厚さ0のコント
ロールボリュームも含む)で分割され、各コントロール
ボリュームの中心点に153点の格子点を設定した。な
お、本実施形態の物理モデルでは、一端付近の温度変化
が大きいので、この一端付近に多くの格子点を設けるた
めに、一端から長さ0.001mの範囲に格子点の半数
を置いた。下式は、n番目のコントロールボリュームの
厚さを表す式である。d=0.001×(1.0−e
−α)e−nα/(1.0−e−77α)ここで、αは
格子点をどのくらいの密度で配置するかの係数であり、
本実施形態では、α=0.1とし、前記0.001mの
位置から他端の0.01mまではコントロールボリュー
ムを均等に配置した。As shown in the figure, the calculation area is one-dimensional and has a length of 0.01 m, and the calculation area is divided by 153 control volumes (including a control volume having a thickness of 0 at the boundary). Then, 153 grid points were set at the center points of the respective control volumes. In the physical model of the present embodiment, since the temperature change near one end is large, half of the grid points are placed within a range of 0.001 m from one end in order to provide many grid points near this one end. The following formula is a formula expressing the thickness of the n-th control volume. d = 0.001 * (1.0-e
−α ) e −nα /(1.0−e −77α ) where α is a coefficient of how densely the lattice points are arranged,
In this embodiment, α = 0.1, and the control volume is evenly arranged from the position of 0.001 m to the other end of 0.01 m.
【0022】問題としては、時間t=0において、前記
一端の温度が急に上昇し、定常状態において、領域全体
が前記一端の温度と同じになるまでのピストン効果によ
る熱移動を算出するものとし、境界条件としては、同図
に示すように、一端が温度一定(Th)、他端が断熱
(∂T/∂x=0)という条件に設定し、初期条件は、
初期圧力が22.1MPa、前記一端における初期温度
が651K、その他の領域は650Kとし、解析時間と
して、タイムステップは、0.1秒とし、タイムステッ
プ数は任意数設定した。The problem is that the heat transfer due to the piston effect is calculated until the temperature at the one end suddenly rises at time t = 0 and the entire region becomes the same as the temperature at the one end in the steady state. As the boundary condition, as shown in the figure, one end is set to a constant temperature (Th) and the other end is adiabatic (∂T / ∂x = 0), and the initial condition is
The initial pressure was 22.1 MPa, the initial temperature at the one end was 651 K, and the other region was 650 K. As the analysis time, the time step was set to 0.1 seconds, and the number of time steps was set to an arbitrary number.
【0023】次に、ステップ2において、初期温度と初
期圧力から必要とされる物性値を取得する。ここで、初
期温度をT0、初期圧力p0とすると、上記の各物性値
は、熱伝導率λ(T0、p0)、定圧比熱cp(T0、
p0)、密度ρ(T0、p0)、比熱比γ(T0、
p0)(定積比熱cV(T0、p0))、体膨張率αp
(T0、p0)、等温圧縮率χT(T0、p0)、[∂
T/∂p]ρ(T0、p0)で表せる。なお、本実施形
態では、これらの物性値は、1997に作成されたIA
PWS実用国際状態式に基づいて作成された日本機械学
会蒸気表を参照して取得するように設定した。Next, in step 2, required physical property values are acquired from the initial temperature and the initial pressure. Here, the initial temperature T 0, when the initial pressure p 0, the physical property values described above, the thermal conductivity λ (T 0, p 0) , the constant-pressure specific heat c p (T 0,
p 0 ), density ρ (T 0 , p 0 ), specific heat ratio γ (T 0 ,
p 0 ) (constant volume specific heat c V (T 0 , p 0 )), body expansion coefficient α p
(T 0 , p 0 ), isothermal compressibility χ T (T 0 , p 0 ), [∂
T / ∂p] ρ (T 0 , p 0 ). In addition, in the present embodiment, these physical property values are IA created in 1997.
It was set so as to be acquired by referring to the steam table of the Japan Society of Mechanical Engineers created based on the PWS practical international state equation.
【0024】次に、ステップ3において、上記(1)式
において、前段の熱伝導項と後段の発熱項のうち、前段
の熱伝導項、∂T/∂t=∇(λ(T、p)・∇T)/
ρ(T、p)・cp(T、p)に基づいて、領域全体の
温度T1を求める。ここで、物性値として、熱伝導率λ
(T0、p0)、密度ρ(T0、p0)、定圧比熱cp
(T0、p0)を用いると、上式は下式のように表され
る。
∂T/∂t=∇(λ(T0、p0)・∇T)/ρ(T0、p0)・
cp(T0、p0)・・・(3−1)
温度T1の算出に当たって、図3に示すように、コント
ロールボリュームの任意のコントロールボリュームにお
いて、任意の厚さ△xのコントロールボリュームおよび
その格子点S、およびそのコントロールボリュームに隣
接するコントロールボリュームの格子点W、E、格子点
Wと格子点S間および格子点Sと格子点E間の差分を
(δx)w、(δx)e、各コントロールボリューム間
の界面をwおよびeとし、格子点Sにおける過去(既
知)の温度、例えば、初期温度をTS 0とし、格子点
W、S、Eにおける温度をそれぞれ、TW 1、TS 1、
TE 1とすると、(3−1)式を離散化して表すと、下
式のように表される。
aSTS 1=aETE 1+aWTW 1+b・・・(4−1)
ここで、aE=λe(T0、p0)/(δx)e、aW
=λw(T0、p0)/(δx)w、aS 0=ρ
(T0、p0)・cp(T0、p0)・△x/△t、b
=aS・TS 0、aS=aE+aW+aS 0 、△tは
タイムステップ間隔であり、密度ρ(T0、p0)およ
び定圧比熱cp(T0、p0)はコントロールボリュー
ムSの物性値を使用する。また熱伝導率λe(T0、p
0)、λw(T 0、p0)については、コントロールボ
リュームでの界面w、eでの値が必要であり、実際に求
められているのは、コントロールボリュームW、Eの値
であるので、補間計数をfe≡(δx)e+/(δx)
eと定義すると、図4に示すように、界面eにおける熱
伝導率λe(T0、p0)は、λe(T0、p0)=
((1−fe)/λS(T0、p0)+fe/λE(T
0、p0))で表される。従って、(4−1)式を利用
して、詳細は省略するが、領域全体の各格子点の温度を
未知数として連立方程式(マトリックス)を作成し、さ
らに、これに境界条件を入れて連立方程式を解くと、初
期温度T0から熱伝導のみを考慮した、最初のタイムス
テップにおける領域全体の温度T1が求められる。Next, in step 3, the above equation (1) is used.
Of the heat conduction term of the first stage and the heat generation term of the second stage,
Heat conduction term of ∂T / ∂t = ∇ (λ (T, p) ・ ∇T) /
ρ (T, p) · cpBased on (T, p)
Temperature T1Ask for. Here, as the physical property value, the thermal conductivity λ
(T0, P0), Density ρ (T0, P0), Constant pressure specific heat cp
(T0, P0), The above equation is expressed as
It
∂T / ∂t = ∇ (λ (T0, P0) ・ ∇T) / ρ (T0, P0) ・
cp(T0, P0) ・ ・ ・ (3-1)
Temperature T1As shown in Fig. 3, when calculating
For any control volume of the roll volume
And a control volume of arbitrary thickness Δx and
Next to the grid point S and its control volume
Grid points W, E, grid points of the control volume that touches
Differences between W and grid point S and between grid point S and grid point E
(Δx)w, (Δx)eBetween each control volume
Let w and e be the interfaces of the
Knowledge) temperature, for example, the initial temperature is TS 0And grid points
The temperatures at W, S, and E are TW 1, TS 1,
TE 1Then, when the equation (3-1) is discretized and expressed,
It is expressed as an expression.
aSTS 1= AETE 1+ AWTW 1+ B ... (4-1)
Where aE= Λe(T0, P0) / (Δx)e, AW
= Λw(T0, P0) / (Δx)w, AS 0= Ρ
(T0, P0) ・ Cp(T0, P0) ・ Δx / Δt, b
= AS・ TS 0, AS= AE+ AW+ AS 0 , Δt
Is the time step interval, and the density ρ (T0, P0) And
And constant pressure specific heat cp(T0, P0) Is the control volume
Use the physical property value of S. And the thermal conductivity λe(T0, P
0), Λw(T 0, P0) For control
The values at the interfaces w and e in the ryme are necessary, and are actually calculated.
What is being set is the value of the control volume W, E
Therefore, the interpolation count is fe≡ (δx)e +/ (Δx)
e, The heat at the interface e is defined as shown in FIG.
Conductivity λe(T0, P0) Is λe(T0, P0) =
((1-fe) / ΛS(T0, P0) + Fe/ ΛE(T
0, P0)). Therefore, use equation (4-1)
Then, although details are omitted, the temperature of each grid point in the entire area is
Create simultaneous equations (matrix) as unknowns,
In addition, if the simultaneous equations are solved by adding boundary conditions to this,
Period temperature T0From the first time
Temperature T of the whole area at the step1Is required.
【0025】次に、ステップ4にて、上記(2)式、∂
p/∂t=∫Vρ(T、p)・αp(T、p)・[∂T
/∂t]・dV/∫Vρ(T、p)・χT(T、p)・
dVを用いて圧力pを求める。(2)式は、物性値とし
て、各コントロールボリュームおける密度ρ(T0、p
0)、体膨張率αp(T0、p0)、等温圧縮率χ
T(T0、p0)、および各格子点における熱伝導のみ
を考慮した温度T1に基づく[∂T1/∂t]を用いる
と、下式のように表される。
∂p/∂t=∫Vρ(T0、p0)αp(T0、p0)・[∂T1/∂t]
・dV/∫Vρ(T0、p0)・χT(T0、p0)・dV・・
・・(5−1)
この(5−1)式を数値積分して圧力p1を求める。Next, in step 4, the above equation (2), ∂
p / ∂t = ∫ V ρ ( T, p) · α p (T, p) · [∂T
/ ∂t] · dV / ∫ V ρ (T, p) · χ T (T, p) ·
The pressure p is determined using dV. Equation (2) is the physical property value of the density ρ (T 0 , p in each control volume).
0 ), body expansion rate α p (T 0 , p 0 ), isothermal compression rate χ
Using T (T 0 , p 0 ) and [∂T 1 / ∂t] based on the temperature T 1 in which only the heat conduction at each lattice point is taken into consideration, the following expression is obtained. ∂p / ∂t = ∫ V ρ ( T 0, p 0) α p (T 0, p 0) · [∂T 1 / ∂t] · dV / ∫ V ρ (T 0, p 0) · χ T (T 0 , p 0 ) · dV ···· (5-1) This equation (5-1) is numerically integrated to obtain the pressure p 1 .
【0026】この式で特に注目すべき点は、圧力の時間
変化であり、領域内で温度が上がったことにより、圧力
がどれだけ上昇したかを把握するものである。ここで数
値積分は、ピストン効果を精度良く反映させるため、精
度の良いものでなければならない。また、温度が求めら
れている格子点の間隔が等間隔ではないことから、不等
間隔を扱えなければならない。従って、この計算におい
て不等間隔を扱うことのできるシンプソン(simps
on)の公式を用いて、領域全体を積分して圧力を求め
る。The point of particular interest in this equation is the change in pressure over time, and it is to grasp how much the pressure has risen due to the rise in temperature within the region. Here, the numerical integration must be accurate because it accurately reflects the piston effect. In addition, since the lattice points for which the temperature is required are not evenly spaced, it is necessary to handle unequal spacing. Therefore, in this calculation, simpson (simps) that can handle unequal intervals
on) to calculate the pressure by integrating over the entire area.
【0027】次に、ステップ5にて、上記(1)式の後
段の発熱項、Q=[1−cV/cp]・[∂T/∂p]
ρ・∂p/∂tから発熱量Qを求める。ここで、物性値
として、定積比熱cV(T0、p0)、定圧比熱c
p(T0、p0)、密度ρ(T0、p0)、比熱比γ
(T0、p0)(定積比熱cV(T0、p0))、体膨
張率αp(T0、p0)、等温圧縮率χT(T0、
p0)、[∂T/∂p]ρ(T0、p0)を用いると、
発熱量Q1は、下式で表される。
Q1=[1−cV(T0、p0)/cp(T0、p0)]・
[∂T/∂p]ρ(T0、p0)・∂p/∂t・・・(6−1)
ここで、圧力pは重力による影響の部分を除き、領域全
体で一定とし、本実施形態でも、圧力pは領域全体で一
定としている。Next, in step 5, above (1) in the subsequent stage of the heating section, Q = [1-c V / c p] · [∂T / ∂p]
The calorific value Q is calculated from ρ · ∂p / ∂t. Here, as the physical property values, constant volume specific heat c V (T 0 , p 0 ), constant pressure specific heat c
p (T 0 , p 0 ), density ρ (T 0 , p 0 ), specific heat ratio γ
(T 0 , p 0 ) (constant volume specific heat c V (T 0 , p 0 )), body expansion rate α p (T 0 , p 0 ), isothermal compression rate χ T (T 0 ,
p 0 ), [∂T / ∂p] ρ (T 0 , p 0 )
The calorific value Q 1 is represented by the following formula. Q 1 = [1-c V (T 0, p 0) / c p (T 0, p 0)] · [∂T / ∂p] ρ (T 0, p 0) · ∂p / ∂t ·· (6-1) Here, the pressure p is constant over the entire region except for the portion affected by gravity, and in this embodiment also, the pressure p is constant over the entire region.
【0028】次に、ステップ6では、ステップ3および
ステップ4にてそれぞれ求められた温度T1および圧力
p1を用いて、各物性値を求める。これらの物性値はス
テップ2において求めたと同様に、IAPWS実用国際
状態式に基づいて作成された日本機械学会蒸気表から、
熱伝導率λ(T1、p1)、定圧比熱cp(T1、
p1)、密度ρ(T1、p1)、比熱比γ(T1、
p1)(定積比熱cV(T1、p1)、体膨張率α
p(T1、p1)、等温圧縮率χ T(T1、p1)、
[∂T/∂p]ρ(T1、p1)を取得する。Next, in step 6, step 3 and
Temperature T obtained in step 41And pressure
p1Each physical property value is calculated using. These physical properties are
Similar to the one obtained in Step 2, IAPWS Practical International
From the Japan Society of Mechanical Engineers steam table created based on the state equation,
Thermal conductivity λ (T1, P1), Constant pressure specific heat cp(T1,
p1), Density ρ (T1, P1), Specific heat ratio γ (T1,
p1) (Constant volume specific heat cV(T1, P1), Body expansion rate α
p(T1, P1), Isothermal compressibility χ T(T1, P1),
[∂T / ∂p]ρ(T1, P1) To get.
【0029】次に、ステップ7においてピストン効果を
考慮した温度Tを求める。熱伝導のみを考慮して求めら
れた温度T1に対して、ピストン効果をも考慮した温度
をT2とし、(1)式から温度T2を求める。算出法
は、ステップ3で求めたと同様の手法によるが、ステッ
プ3で求めた方法と比べて、発熱量Q1が加算される点
で相違する。(1)式を離散化した式は、下式のように
表される。
aSTS 2=aETE 2+aWTW 2+b・・・(7−1)
ここで、aE=λe(T1、p1)/(δx)e、aW
=λw(T1、p1)/(δx)w、aS 1=ρ
(T1、p1)・cp(T1、p1)・△x/△t、b
=aS 1・TS 1+Q1Δx、aS=aE+aW+aS
1 、△tはタイムステップ間隔である。ここでも詳細
は省略するが、(7−1)式を利用して、領域全体につ
いて連立方程式(マトリックス)を作成し、さらに、こ
れに境界条件を入れて連立方程式を解くと、温度T1か
らピストン効果をも考慮した温度T2が求められる。Next, in step 7, the temperature T considering the piston effect is obtained. With respect to the temperature T 1 obtained by considering only the heat conduction, the temperature T 2 which also takes the piston effect into consideration is set as T 2, and the temperature T 2 is obtained from the equation (1). The calculation method is the same as that obtained in step 3, but differs from the method obtained in step 3 in that the calorific value Q 1 is added. The equation obtained by discretizing the equation (1) is expressed as the following equation. a S T S 2 = a E T E 2 + a W T W 2 + b (7-1) Here, a E = λ e (T 1 , p 1 ) / (δx) e , a W
= Λ w (T 1 , p 1 ) / (δx) w , a S 1 = ρ
(T 1, p 1) · c p (T 1, p 1) · △ x / △ t, b
= A S 1 · T S 1 + Q 1 Δx, a S = a E + a W + a S
1 and Δt are time step intervals. Although details are omitted here, if equations (7-1) are used to create simultaneous equations (matrix) for the entire region, and boundary equations are added to these to solve the simultaneous equations, the temperature T 1 is calculated. The temperature T 2 that takes into consideration the piston effect is also required.
【0030】ピストン効果は、領域内のどこかでわずか
でも温度変化(上昇、下降)のあることが必要である
が、本発明の算出方法では熱伝導とピストン効果とによ
る熱移動現象を想定しており、この熱伝導による温度変
化がピストン効果を引き起こし、発熱項Qが領域全体に
影響を及ぼすものであるからピストン効果を表している
といえる。The piston effect requires a slight temperature change (rise, fall) anywhere in the region, but the calculation method of the present invention assumes a heat transfer phenomenon due to heat conduction and the piston effect. Therefore, it can be said that the temperature change due to the heat conduction causes the piston effect, and the heat generation term Q affects the entire region, so that it represents the piston effect.
【0031】次に、ステップ8にて、n+1回目、n回
目の計算でそれぞれ求められた温度Tn+1、Tnとす
るとき、所定値εに対して全ての格子点で|(Tn+1
−Tn)/Tn|<εを判定する。温度T2と温度T1
間では収束することはあり得ないので、温度T2と温度
T1しか求められていない段階では、収束判定を行うこ
となく、ステップ4からの処理を行うようにしてもよ
い。その場合は、ステップ8での収束判定は|(T3−
T2)/T2|<εからとなる。なお、収束していると
判定されたときはステップ9の処理を行い、収束してい
ると判定されないときはステップ4からの処理を繰り返
す。なお、εとしては0.01%(10−5)を用い
た。Next, in step 8, assuming that the temperatures T n + 1 and T n are respectively obtained by the n + 1-th and n-th calculations, | (T n + 1) at all lattice points with respect to a predetermined value ε.
-T n ) / T n | <ε is determined. Temperature T 2 and temperature T 1
Since there is no possibility of convergence during the period, the processing from step 4 may be performed without performing the convergence determination at the stage where only the temperatures T 2 and T 1 are obtained. In that case, the convergence determination in step 8 is | (T 3 −
T 2 ) / T 2 | <ε. If it is determined that the convergence has occurred, the process of step 9 is performed, and if it is not determined that the convergence has occurred, the processes from step 4 are repeated. As ε, 0.01% (10 −5 ) was used.
【0032】ここで、収束していないと判定されたとき
のステップ4での処理は、先のステップ4での処理と同
様に、各コントロールボリュームおける密度ρ(T1、
p1)、体膨張率αp(T1、p1)、等温圧縮率χT
(T1、p1)、および各格子点における温度T2に基
づく[∂T2/∂t]を用い、下式を数値積分して、圧
力p2を求める。
∂p/∂t=∫Vρ(T1、p1)αp(T1、p1)・[∂T2/∂t]・
dV/∫Vρ(T1、p1)・χT(T1、p1)・dV・・
・・(5−2)。Here, the processing in step 4 when it is determined that they have not converged is similar to the processing in step 4 above, the density ρ (T 1 ,
p 1 ), body expansion rate α p (T 1 , p 1 ), isothermal compression rate χ T
Using (T 1 , p 1 ) and [∂T 2 / ∂t] based on the temperature T 2 at each grid point, the following formula is numerically integrated to obtain the pressure p 2 . ∂p / ∂t = ∫ V ρ (T 1 , p 1 ) α p (T 1 , p 1 ) ・ [∂T 2 / ∂t] ・ dV / ∫ V ρ (T 1 , p 1 ) ・ χ T (T 1 , p 1 ) · dV ···· (5-2).
【0033】次いで、ステップ5にて、先と同様にし
て、下式から発熱量Q2を求める。Q2=[1−c
V(T1、p1)/cp(T1、p1)]・
[∂T/∂p]ρ(T1、p1)・∂p/∂t・・・(6−2)。Then, in step 5, the calorific value Q 2 is calculated from the following equation in the same manner as above. Q 2 = [1-c
V (T 1, p 1) / c p (T 1, p 1)] · [∂T / ∂p] ρ (T 1, p 1) · ∂p / ∂t ··· (6-2).
【0034】次いで、ステップ6にて、先と同様にし
て、ステップ3およびステップ4にてそれぞれ求められ
た温度T2および圧力p2を用いて、各物性値を求め
る。これらの物性値はステップ2において求めたと同様
の手法により、熱伝導率λ(T2、p2)、定圧比熱c
p(T2、p2)、密度ρ(T2、p2)、比熱比γ
(T2、p2)(定積比熱cV(T2、p2)、体膨張
率αp(T2、p2)、等温圧縮率χT(T2、
p2)、[∂T/∂p]ρ(T2、p2)を求める。Then, in step 6, each physical property value is obtained using the temperature T 2 and the pressure p 2 obtained in step 3 and step 4, respectively, in the same manner as above. These physical property values are calculated by the same method as that obtained in step 2, the thermal conductivity λ (T 2 , p 2 ), the constant pressure specific heat c
p (T 2 , p 2 ), density ρ (T 2 , p 2 ), specific heat ratio γ
(T 2 , p 2 ) (constant volume specific heat c V (T 2 , p 2 ), body expansion rate α p (T 2 , p 2 ), isothermal compression rate χ T (T 2 ,
p 2 ) and [∂T / ∂p] ρ (T 2 , p 2 ) are obtained.
【0035】さらに、ステップ7にて、先と同様に、下
式からピストン効果を考慮した温度T3を求める。
aSTS 3=aETE 3+aWTW 3+b・・・(7−2)
ここで、aE=λe(T2、p2)/(δx)e、aW
=λw(T2、p2)/(δx)w、aS 2=ρ
(T2、p2)・cp(T2、p2)・△x/△t、b
=aS 2TS 2+Q2Δx、aS=aE+aW+aS 2
となり、詳細は省略するが、(7−2)式を利用し
て、領域全体について連立方程式(マトリックス)を作
成し、さらに、これに境界条件を入れて連立方程式を解
くと、温度T2からピストン効果をも考慮した温度T3
が求められる。Further, in step 7, similarly to the above, the temperature T 3 considering the piston effect is obtained from the following equation. a S T S 3 = a E T E 3 + a W T W 3 + b (7-2) where a E = λ e (T 2 , p 2 ) / (δx) e , a W
= Λ w (T 2 , p 2 ) / (δx) w , a S 2 = ρ
(T 2, p 2) · c p (T 2, p 2) · △ x / △ t, b
= A S 2 T S 2 + Q 2 Δx, a S = a E + a W + a S 2
Although the details are omitted, if equations (7-2) are used to create simultaneous equations (matrix) for the entire region, and then the simultaneous equations are solved by including boundary conditions in these, the temperature T 2 is calculated. Temperature T 3 that also considers the piston effect
Is required.
【0036】次いで、ステップ8で再び、所定値εに対
して、全ての格子点で|(T3−T2)/T2|<εを
判定することにより、温度T3と温度T2間の収束を判
定する。収束していないと判定されたときは、上記に述
べたように、再度ステップ4からの処理を繰り返すが、
このような処理を繰り返して、ステップ8での収束判定
において、全ての格子点で|(Tn+1−Tn)/Tn
|<εを満足していると判定された時は、ステップ9に
移行する。Then, in step 8, again, for a predetermined value ε, by determining | (T 3 −T 2 ) / T 2 | <ε at all the lattice points, the temperature between the temperatures T 3 and T 2 is determined. Determine the convergence of. If it is determined that they have not converged, as described above, the process from step 4 is repeated,
By repeating such processing, in the convergence determination in step 8, | (T n + 1 −T n ) / T n at all grid points
If it is determined that | <ε is satisfied, the process proceeds to step 9.
【0037】ステップ9では、各ステップ、例えば、ス
テップ8での処理の結果算出された各コントロールボリ
ュームにおける全ての格子点の温度Tn+1を記憶し、
また必要に応じて、任意のタイムステップ毎にグラフ表
示する。In step 9, the temperatures T n + 1 of all the lattice points in each control volume calculated as a result of the processing in each step, for example, step 8 are stored,
If necessary, a graph is displayed for each arbitrary time step.
【0038】次に、ステップ10では、次のタイムステ
ップの処理を行うか否かを判定し、次のタイムステップ
での処理を行う場合は、それまで算出された算出結果を
初期値として、ステップ3において、再び熱伝導のみの
新たな温度T1を算出するための処理を開始する。設定
された所定のタイムステップ数の処理を終えたとき全て
の処理を終了する。Next, at step 10, it is determined whether or not the processing at the next time step is to be performed. When the processing at the next time step is to be performed, the calculation result calculated up to that point is used as an initial value to perform the step. In 3, the process for calculating the new temperature T 1 only for heat conduction is started again. When the processing of the set predetermined number of time steps is completed, all the processing is ended.
【0039】図5は、本発明の超臨界流体のピストン効
果による熱移動の算出方法によって算出された結果を示
す図である。同図において、横軸は領域(0.000〜
0.010m)の位置x(m)、縦軸は温度T(K)、
t(秒)=0.1,1.0,10,100はそれぞれ算
出処理の経過時間を表す。FIG. 5 is a diagram showing the results calculated by the method for calculating the heat transfer by the piston effect of the supercritical fluid of the present invention. In the figure, the horizontal axis represents the region (0.000 to 0).
0.010 m) position x (m), the vertical axis represents temperature T (K),
t (seconds) = 0.1, 1.0, 10, 100 respectively represent the elapsed time of the calculation process.
【0040】これは、超臨界水(臨界点温度Tc=64
7.096K、圧力pc=22.064MPa)では、
初期温度T0=650K、初期圧力p0=22.1MP
a、t=0において一端(x=0.000)の温度が急
に651Kになった場合の温度分布の時間変化を示すも
のであり、ピストン効果は熱の伝わりが速く、領域全体
が同じ温度になるという特性があるが、この結果は、そ
の状況を良く表している。特に、t=0.1秒後でも、
t=1.0秒後でも他端(x=0.010)で温度上昇
が見られ、これは熱伝導だけでは考えられない現象であ
り、また、境界から徐々に熱が伝わるのではなくて、領
域全体に温度が伝わるという傾向が見られ、これは微小
重力実験で得られた結果とよく一致し、本発明の超臨界
流体の熱移動の算出方法が、ピストン効果による熱移動
の極めて有効な算出方法であることが解る。This is due to supercritical water (critical point temperature Tc = 64
At 7.096 K and pressure pc = 22.064 MPa),
Initial temperature T 0 = 650K, initial pressure p 0 = 22.1MP
It shows the change over time in the temperature distribution when the temperature at one end (x = 0.000) suddenly reaches 651K at a and t = 0. The piston effect shows that heat transfer is fast and the entire region has the same temperature. However, this result shows the situation well. In particular, even after t = 0.1 seconds,
Even after t = 1.0 seconds, a temperature rise was observed at the other end (x = 0.010), which is a phenomenon that cannot be considered only by heat conduction, and heat is not gradually transferred from the boundary. , There is a tendency that the temperature is transmitted to the entire region, which is in good agreement with the result obtained in the microgravity experiment, and the calculation method of the heat transfer of the supercritical fluid of the present invention is extremely effective for the heat transfer by the piston effect. It turns out that this is a simple calculation method.
【0041】なお、図6は本発明の超臨界流体のピスト
ン効果による熱移動の算出方法を実現するためのハード
ウェア構成の一例を示す図である。同図において、1は
コンピュータプログラムに従って各種の演算を実行する
中央演算処理装置、2は演算処理に必要なデータ、例え
ば、熱移動の算出に必要な解析領域、解析時間、初期条
件、境界条件等の各種の情報を入力する入力装置、3
は、超臨界流体のピストン効果による熱移動の算出を実
行するプログラム等を記憶するメインメモリ、4は超臨
界流体のピストン効果による熱移動の算出において必要
なデータ、例えば、算出された温度Tや圧力pに基づい
て物性値を取得するために参照される日本機械学会蒸気
表等のデータや、超臨界流体の熱移動算出の結果が記憶
されるファイル装置、5は超臨界流体のピストン効果に
よる熱移動の算出結果等を画像表示したりプリントアウ
トする出力装置である。FIG. 6 is a diagram showing an example of a hardware configuration for realizing the method for calculating the heat transfer by the piston effect of the supercritical fluid according to the present invention. In the figure, 1 is a central processing unit for executing various calculations according to a computer program, 2 is data required for the calculation processing, for example, an analysis region required for calculation of heat transfer, analysis time, initial conditions, boundary conditions, etc. Input device for inputting various information of 3
Is a main memory for storing a program or the like for calculating the heat transfer due to the piston effect of the supercritical fluid, and 4 is data necessary for calculating the heat transfer due to the piston effect of the supercritical fluid, such as the calculated temperature T or File device 5 in which data such as the steam table of the Japan Society of Mechanical Engineers referred to for obtaining the physical property value based on the pressure p and the result of heat transfer calculation of the supercritical fluid are stored This is an output device that displays an image or prints out the calculation result of heat transfer and the like.
【0042】[0042]
【発明の効果】請求項1ないし請求項3に記載の発明に
よれば、はじめに熱伝導を算出し、それに基づいてピス
トン効果を算出するようにしているので、ピストン効果
が、わずかでも温度上昇(ないし温度下降)しなければ
発生しない性質を見極めた算出方法が得られる。また、
温度、圧力の変化に伴って変化する物性値の変化を逐次
算出過程に反映させ、かつコンピュータ処理によって高
精度でかつ容易に算出するようにしたので、工業化レベ
ルでの超臨界流体におけるピストン効果による熱移動の
算出が容易となり、各種の工業分野での応用が可能とな
る。According to the invention described in claims 1 to 3, the heat transfer is first calculated and the piston effect is calculated based on the calculated heat transfer. It is possible to obtain a calculation method that identifies the properties that do not occur unless the temperature decreases). Also,
Changes in physical properties that change with changes in temperature and pressure are sequentially reflected in the calculation process, and are calculated with high precision and ease by computer processing. Calculation of heat transfer becomes easy, and application in various industrial fields becomes possible.
【0043】請求項2に記載の発明によれば、収束判定
を行いながら、ピストン効果による熱移動を算出してい
るので、正確な算出結果が得られる。According to the second aspect of the present invention, since the heat transfer due to the piston effect is calculated while performing the convergence determination, an accurate calculation result can be obtained.
【0044】請求項3に記載の発明によれば、タイムス
テップ数を適正な値に設定することにより、各時間毎の
詳細な超臨界流体におけるピストン効果による熱移動に
よる詳細な温度分布を得ることができる。According to the third aspect of the present invention, by setting the number of time steps to an appropriate value, it is possible to obtain a detailed temperature distribution by heat transfer due to the piston effect in the detailed supercritical fluid at each time. You can
【0045】請求項4記載の発明によれば、はじめに熱
伝導を算出し、それに基づいてピストン効果を算出する
ようにしているので、ピストン効果が、わずかでも温度
上昇(ないし温度下降)しなければ発生しない性質を見
極めた算出プログラムが得られる。また、温度、圧力の
変化に伴って変化する物性値の変化を逐次算出過程に反
映させ、また収束判定を行いながら、ピストン効果によ
る熱移動を算出しているので、ピストン効果による熱移
動の正確な算出結果が得られる。その結果、工業化レベ
ルでの超臨界流体におけるピストン効果による熱移動の
算出が容易となり、各種の工業分野での応用が可能とな
る。According to the fourth aspect of the present invention, the heat conduction is first calculated and the piston effect is calculated based on the heat conduction. Therefore, if the piston effect does not increase in temperature at all (or decrease in temperature). It is possible to obtain a calculation program that identifies properties that do not occur. In addition, changes in physical properties that change with changes in temperature and pressure are sequentially reflected in the calculation process, and heat transfer due to the piston effect is calculated while performing convergence determination, so the heat transfer due to the piston effect is accurate. It is possible to obtain various calculation results. As a result, it becomes easy to calculate the heat transfer due to the piston effect in the supercritical fluid at the industrial level, and the application in various industrial fields becomes possible.
【0046】請求項5記載の発明によれば、記録媒体に
本発明に係るプログラムを記録して頒布することによ
り、広く本発明の利用が可能となる。According to the invention of claim 5, the program according to the present invention is recorded in a recording medium and distributed, whereby the present invention can be widely used.
【図1】本件発明に係る超臨界流体におけるピストン効
果による熱移動の算出方法を示すフローチャートであ
る。FIG. 1 is a flowchart showing a method for calculating heat transfer due to a piston effect in a supercritical fluid according to the present invention.
【図2】格子点、コントロールボリューム、タイムステ
ップを示した超臨界流体のピストン効果による熱移動を
算出するための物理モデルを示す図である。FIG. 2 is a diagram showing a physical model for calculating heat transfer due to a piston effect of a supercritical fluid showing lattice points, control volumes, and time steps.
【図3】格子点W、S、Eとコントロールボリュームと
の関係を示す図である。FIG. 3 is a diagram showing a relationship between grid points W, S, E and a control volume.
【図4】格子点S、Eとコントロールボリュームの界面
eとの関係を示す図である。FIG. 4 is a diagram showing a relationship between lattice points S and E and an interface e of a control volume.
【図5】本発明に係る超臨界流体のピストン効果による
熱移動の算出方法によって算出された温度分布の時間変
化を示す図である。FIG. 5 is a diagram showing a time change of a temperature distribution calculated by a method of calculating heat transfer due to a piston effect of a supercritical fluid according to the present invention.
【図6】本発明の超臨界流体におけるピストン効果によ
る熱移動の算出方法を実現するためのハードウェア構成
の一例を示す図である。FIG. 6 is a diagram showing an example of a hardware configuration for realizing a method of calculating heat transfer due to a piston effect in a supercritical fluid of the present invention.
1 中央演算処理装置 2 入力装置 3 メインメモリ 4 ファイル装置 5 出力装置 1 Central processing unit 2 input devices 3 main memory 4 file device 5 Output device
───────────────────────────────────────────────────── フロントページの続き (72)発明者 相澤 崇史 宮城県仙台市宮城野区苦竹4−2−1 独 立行政法人産業技術総合研究所 東北セン ター内 (72)発明者 金久保 光央 宮城県仙台市宮城野区苦竹4−2−1 独 立行政法人産業技術総合研究所 東北セン ター内 (72)発明者 斎藤 功夫 宮城県仙台市宮城野区苦竹4−2−1 独 立行政法人産業技術総合研究所 東北セン ター内 Fターム(参考) 2G040 AA01 AB08 AB12 BA22 HA15 HA16 5B056 BB00 HH00 ─────────────────────────────────────────────────── ─── Continued front page (72) Inventor Takashi Aizawa 4-2-1 Maitake, Miyagino-ku, Sendai City, Miyagi Prefecture, Germany National Institute of Advanced Industrial Science and Technology Tohoku Sen Inside (72) Inventor Mitsuo Kanakubo 4-2-1 Maitake, Miyagino-ku, Sendai City, Miyagi Prefecture, Germany National Institute of Advanced Industrial Science and Technology Tohoku Sen Inside (72) Inventor Isao Saito 4-2-1 Maitake, Miyagino-ku, Sendai City, Miyagi Prefecture, Germany National Institute of Advanced Industrial Science and Technology Tohoku Sen Inside F-term (reference) 2G040 AA01 AB08 AB12 BA22 HA15 HA16 5B056 BB00 HH00
Claims (5)
条件、および境界条件の各情報を入力する第1のステッ
プと、 初期温度T0および初期圧力p0に基づいて物性値を取
得する第2のステップと、 前記第2のステップで取得された物性値を用いて、 ∂T/∂t=∇(λ(T、p)・∇T)/ρ(T、p)
・cp(T、p) に基づいて熱伝導による温度Tを算出する第3のステッ
プと、 前記第2のステップまたは後述する第6のステップで取
得された物性値を用いて、 ∂p/∂t=∫Vρ(T、p)・αp(T、p)・[∂
T/∂t]・dV/∫Vρ(T、p)・χT(T、p)
・dV に基づいて圧力pを算出する第4のステップと、 前記第2のステップまたは後述する第6のステップで取
得された物性値を用いて、 Q=[1−cV(T、p)/cp(T、p)]・[∂T
/∂p]ρ(T、p)・∂p/∂t に基づいて発熱量Qを算出する第5のステップと、 前記第3のステップで算出された温度Tおよび前記第4
のステップで算出された圧力pに基づいて物性値を取得
する第6のステップと、 第5のステップで算出された発熱量Qおよび第6のステ
ップで取得された物性値を用いて、 ∂T/∂t=∇(λ(T、p)・∇T)/ρ(T、p)
・cp(T、p)+[1−cV(T、p)/cp(T、
p)]・[∂T/∂p]ρ(T、p)・∂p/∂t に基づいてピストン効果を考慮した温度Tを算出する第
7のステップと、 からなることを特徴とする超臨界流体におけるピストン
効果による熱移動の算出方法。 ただし、上記の各物性値は、熱伝導率λ(T、p)、定
圧比熱cp(T、p)、密度ρ(T、p)、定積比熱c
V(T、p)、体膨張率αp(T、p)、等温圧縮率χ
T(T、p)、[∂T/∂p]ρ(T、p)。1. A physical property value is obtained based on a first step of inputting each information of an analysis region, an analysis time, an initial condition, and a boundary condition necessary for calculation, and an initial temperature T 0 and an initial pressure p 0. Using the second step and the physical property values obtained in the second step, ∂T / ∂t = ∇ (λ (T, p) · ∇T) / ρ (T, p)
∂p / using the third step of calculating the temperature T due to heat conduction based on c p (T, p) and the physical property values obtained in the second step or the sixth step described later. ∂t = ∫ V ρ (T, p) · α p (T, p) · [∂
T / ∂t] ・ dV / ∫ V ρ (T, p) ・ χ T (T, p)
- a fourth step of calculating the pressure p on the basis of the dV, using the sixth to physical property values obtained in the step of the second step or later, Q = [1-c V (T, p) / c p (T, p) ] · [∂T
/ ∂p] ρ (T, p) · ∂p / ∂t based on the fifth step, the temperature T calculated in the third step and the fourth step
∂T using the sixth step of acquiring the physical property value based on the pressure p calculated in the step, the calorific value Q calculated in the fifth step, and the physical property value acquired in the sixth step. / ∂t = ∇ (λ (T, p) ・ ∇T) / ρ (T, p)
· C p (T, p) + [1-c V (T, p) / c p (T,
p)] · [∂T / ∂p] ρ (T, p) · ∂p / ∂t based on the seventh step of calculating the temperature T considering the piston effect, and Calculation method of heat transfer due to piston effect in critical fluids. However, the physical properties of the above, thermal conductivity lambda (T, p), pressure specific heat c p (T, p), the density ρ (T, p), specific heat at constant volume c
V (T, p), body expansion rate α p (T, p), isothermal compression rate χ
T (T, p), [∂T / ∂p] ρ (T, p).
について、経時的に算出された最新の温度Tn+1と経
時的に前記最新の温度Tn+1の直前に算出された温度
Tnとに基づいて両者間の収束を判定し、収束していな
いと判定されたときは前記第4のステップからの処理を
繰り返させ、収束していると判定されたときは後述する
第9のステップの処理を行わせる第8のステップと、 各ステップの算出結果を記憶ないし表示する第9のステ
ップと、 からなることを特徴とする請求項1に記載の超臨界流体
の熱移動算出方法。2. The temperature T calculated in the seventh step
For determines convergence between them on the basis of the temperature T n calculated immediately before the time the calculated latest temperature T n + 1 and over time the latest temperature T n + 1, and does not converge determination When it is determined that the process from the fourth step is repeated, and when it is determined that the process is converged, an eighth step for performing the process of a ninth step described later and a calculation result of each step are stored. Or a ninth step of displaying, and the method for calculating heat transfer of a supercritical fluid according to claim 1, wherein:
イムステップ数を処理したかを判定し、処理していない
と判定されたときは前記第3のステップからの処理を行
わせる第10のステップと、 からなることを特徴とする請求項2に記載の超臨界流体
におけるピストン効果による熱移動の算出方法。3. After the processing of the ninth step, it is determined whether or not a predetermined number of time steps have been processed, and when it is determined that the predetermined number of time steps have not been processed, the processing from the third step is performed. The method for calculating heat transfer due to the piston effect in a supercritical fluid according to claim 2, characterized in that:
条件、および境界条件の各情報を入力する第1のステッ
プと、 初期温度T0および初期圧力p0に基づいて物性値を取
得する第2のステップと、 前記第2のステップで取得された物性値を用いて、 ∂T/∂t=∇(λ(T、p)・∇T)/ρ(T、p)
・cp(T、p) に基づいて熱伝導による温度Tを算出する第3のステッ
プと、 前記第2のステップまたは後述する第6のステップで取
得された物性値を用いて、 ∂p/∂t=∫Vρ(T、p)・αp(T、p)・[∂
T/∂t]・dV/∫Vρ(T、p)・χT(T、p)
・dV に基づいて圧力pを算出する第4のステップと、 前記第2のステップまたは後述する第6のステップで取
得された物性値を用いて、 Q=[1−cV(T、p)/cp(T、p)]・[∂T
/∂p]ρ(T、p)・∂p/∂t に基づいて発熱量Qを算出する第5のステップと、 前記第3のステップで算出された温度Tおよび前記第4
のステップで算出された圧力pに基づいて物性値を取得
する第6のステップと、 第5のステップで算出された発熱量Qおよび第6のステ
ップで取得された物性値を用いて、 ∂T/∂t=∇(λ(T、p)・∇T)/ρ(T、p)
・cp(T、p)+[1−cV(T、p)/cp(T、
p)]・[∂T/∂p]ρ(T、p)・∂p/∂t に基づいてピストン効果を考慮した温度Tを算出する第
7のステップと、 前記第7のステップで算出された温度Tについて、経時
的に算出された最新の温度Tn+1と経時的に前記最新
の温度Tn+1の直前に算出された温度Tnとに基づい
て両者間の収束を判定し、収束していないと判定された
ときは前記第4のステップからの処理を繰り返させ、収
束していると判定されたときは後述する第9のステップ
の処理を行わせる第8のステップと、 各ステップの算出結果を記憶ないし表示する第9のステ
ップと、 前記第9のステップの処理後、所定のタイムステップ数
を処理したかを判定し、処理していないと判定されたと
きは前記第3のステップからの処理を行わせる第10の
ステップと、 からなることを特徴とする超臨界流体におけるピストン
効果による熱移動を算出するプログラム。ただし、上記
の各物性値は、熱伝導率λ(T、p)、定圧比熱c
p(T、p)、密度ρ(T、p)、定積比熱cV(T、
p)、体膨張率αp(T、p)、等温圧縮率χT(T、
p)、[∂T/∂p]ρ(T、p)。4. A physical property value is acquired based on a first step of inputting each information of an analysis region, an analysis time, an initial condition, and a boundary condition necessary for calculation, and an initial temperature T 0 and an initial pressure p 0. Using the second step and the physical property values obtained in the second step, ∂T / ∂t = ∇ (λ (T, p) · ∇T) / ρ (T, p)
∂p / using the third step of calculating the temperature T due to heat conduction based on c p (T, p) and the physical property values obtained in the second step or the sixth step described later. ∂t = ∫ V ρ (T, p) · α p (T, p) · [∂
T / ∂t] ・ dV / ∫ V ρ (T, p) ・ χ T (T, p)
- a fourth step of calculating the pressure p on the basis of the dV, using the sixth to physical property values obtained in the step of the second step or later, Q = [1-c V (T, p) / c p (T, p) ] · [∂T
/ ∂p] ρ (T, p) · ∂p / ∂t based on the fifth step, the temperature T calculated in the third step and the fourth step
∂T using the sixth step of acquiring the physical property value based on the pressure p calculated in the step, the calorific value Q calculated in the fifth step, and the physical property value acquired in the sixth step. / ∂t = ∇ (λ (T, p) ・ ∇T) / ρ (T, p)
· C p (T, p) + [1-c V (T, p) / c p (T,
p)] · [∂T / ∂p] ρ (T, p) · ∂p / ∂t based on the seventh step of calculating the temperature T considering the piston effect, and the seventh step. and the temperature T, to determine the convergence between the two on the basis of the temperature T n calculated immediately before the time the calculated latest temperature T n + 1 and over time the latest temperature T n + 1, have converged If it is determined that there is not, the process from the fourth step is repeated, and if it is determined that the process has converged, an eighth step that causes the process of the ninth step described later to be performed, and calculation of each step. A ninth step of storing or displaying the result, and after the processing of the ninth step, it is determined whether or not a predetermined number of time steps have been processed, and when it is determined that the processing has not been performed, the third step is executed. The 10th step When a program for calculating the heat transfer by the piston effect in the supercritical fluid, characterized in that it consists of. However, the above-mentioned physical property values are thermal conductivity λ (T, p), constant pressure specific heat c
p (T, p), the density ρ (T, p), specific heat at constant volume c V (T,
p), the body expansion rate α p (T, p), the isothermal compression rate χ T (T,
p), [∂T / ∂p] ρ (T, p).
コンピュータ読み取り可能な記録媒体。5. A computer-readable recording medium in which the program according to claim 4 is recorded.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001217306A JP2003028819A (en) | 2001-07-17 | 2001-07-17 | Method for calculating heat transfer by piston effect at supercritical fluid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001217306A JP2003028819A (en) | 2001-07-17 | 2001-07-17 | Method for calculating heat transfer by piston effect at supercritical fluid |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2003028819A true JP2003028819A (en) | 2003-01-29 |
Family
ID=19051654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001217306A Pending JP2003028819A (en) | 2001-07-17 | 2001-07-17 | Method for calculating heat transfer by piston effect at supercritical fluid |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2003028819A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103308551A (en) * | 2013-06-05 | 2013-09-18 | 中国石油大学(华东) | Experimental device and method for measuring surface coefficient of heat transfer of supercritical carbon dioxide |
CN103439356A (en) * | 2013-06-05 | 2013-12-11 | 中国石油大学(华东) | Experimental apparatus and method for constant-pressure specific heat measurement of supercritical carbon dioxide |
-
2001
- 2001-07-17 JP JP2001217306A patent/JP2003028819A/en active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103308551A (en) * | 2013-06-05 | 2013-09-18 | 中国石油大学(华东) | Experimental device and method for measuring surface coefficient of heat transfer of supercritical carbon dioxide |
CN103439356A (en) * | 2013-06-05 | 2013-12-11 | 中国石油大学(华东) | Experimental apparatus and method for constant-pressure specific heat measurement of supercritical carbon dioxide |
CN103439356B (en) * | 2013-06-05 | 2016-08-10 | 中国石油大学(华东) | The experimental provision measured for supercritical carbon dioxide specific heat at constant pressure and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Turek et al. | Numerical simulation and benchmarking of a monolithic multigrid solver for fluid-structure interaction problems with application to hemodynamics | |
Ge et al. | Numerical simulation of flow in mechanical heart valves: grid resolution and the assumption of flow symmetry | |
Arina | Numerical simulation of near-critical fluids | |
Borg et al. | Controllers for imposing continuum-to-molecular boundary conditions in arbitrary fluid flow geometries | |
Prates et al. | Inverse strategies for identifying the parameters of constitutive laws of metal sheets | |
Miguélez et al. | Residual stresses in orthogonal cutting of metals: the effect of thermomechanical coupling parameters and of friction | |
US20170011147A1 (en) | Simulation device, simulation program and simulation method for liquid metal | |
Kumar et al. | Micro-structural impact of different strut shapes and porosity on hydraulic properties of Kelvin-like metal foams | |
Mirzaeifar et al. | Coupled thermo-mechanical analysis of shape memory alloy circular bars in pure torsion | |
Mohasseb et al. | A novel approach to solve inverse heat conduction problems: Coupling scaled boundary finite element method to a hybrid optimization algorithm | |
Scheerlinck et al. | Finite element computation of unsteady phase change heat transfer during freezing or thawing of food using a combined enthalpy and Kirchhoff transform method | |
Mihoubi et al. | Transfer phenomena during the drying of a shrinkable product: modeling and simulation | |
Montoya et al. | Sensitivity analysis in thermoelastic problems using the complex finite element method | |
Aly et al. | Modelling of non-Darcy flows through porous media using extended incompressible smoothed particle hydrodynamics | |
Mohammadiun et al. | Application of sequential function specification method in heat flux monitoring of receding solid surfaces | |
Wu et al. | Numerical simulation of subcooled nucleate boiling by coupling level-set method with moving-mesh method | |
Yang et al. | A random field model based on nodal integration domain for stochastic analysis of heat transfer problems | |
Banaszkiewicz | Online determination of transient thermal stresses in critical steam turbine components using a two-step algorithm | |
Zagrouba et al. | Drying of clay. II Rheological modelisation and simulation of physical phenomena | |
Le Page et al. | Development of an approximate empirical-CFD model estimating coupled heat and water transfers of stacked food products placed in airflow | |
JP2003028819A (en) | Method for calculating heat transfer by piston effect at supercritical fluid | |
Singh et al. | Study of combined free convection and surface radiation in closed cavities partially heated from below | |
Duda | Simplification of 3D transient heat conduction by reducing it to an axisymmetric heat conduction problem and a new inverse method of the problem solution | |
Nejat et al. | A Newton–Krylov finite volume algorithm for the power-law non-Newtonian fluid flow using pseudo-compressibility technique | |
CN112084723B (en) | Optical fiber preform one-time stretching process simulation method and device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20041217 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060110 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20060516 |