JP2006258391A - 浮力を伴う乱流の流体的及び熱的諸特性の推定方法及び推定プログラム - Google Patents
浮力を伴う乱流の流体的及び熱的諸特性の推定方法及び推定プログラム Download PDFInfo
- Publication number
- JP2006258391A JP2006258391A JP2005078548A JP2005078548A JP2006258391A JP 2006258391 A JP2006258391 A JP 2006258391A JP 2005078548 A JP2005078548 A JP 2005078548A JP 2005078548 A JP2005078548 A JP 2005078548A JP 2006258391 A JP2006258391 A JP 2006258391A
- Authority
- JP
- Japan
- Prior art keywords
- turbulent
- heat flux
- term
- calculated
- velocity
- 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
Abstract
【課題】熱移動によって発生した浮力を伴う乱流の流体的及び熱的諸特性を推定する方法およびコンピュータによって実行可能な，浮力を伴う乱流の流体的及び熱的諸特性の推定プログラムを提供する．
【解決手段】熱移動を伴う流れ場についての乱流の2方程式モデルにおいて，レイノルズ応力を
，強制対流に関する項をF，浮力に関する項をG
とすると
の式を用いてレイノルズ応力を算出して，速度場の計算をし，乱流熱流束を
，乱流熱流束の強制対流に関する項をFt，乱流熱流束の浮力に関する項をGtとすると，
の式を用いて乱流熱流束を算出して，上記速度場における温度場の計算をすることを特徴とする，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定方法
【選択図】 図１
【解決手段】熱移動を伴う流れ場についての乱流の2方程式モデルにおいて，レイノルズ応力を
，強制対流に関する項をF，浮力に関する項をG
とすると
の式を用いてレイノルズ応力を算出して，速度場の計算をし，乱流熱流束を
，乱流熱流束の強制対流に関する項をFt，乱流熱流束の浮力に関する項をGtとすると，
の式を用いて乱流熱流束を算出して，上記速度場における温度場の計算をすることを特徴とする，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定方法
【選択図】 図１
Description
本発明は，各種機械，電子機器等の流体による放熱設計や，冷暖房時の室内気流の解析等において必要となる，熱移動によって発生した浮力を伴う乱流の流体的及び熱的諸特性を推定する方法およびコンピュータによって実行可能な，浮力を伴う乱流の流体的及び熱的諸特性の推定プログラムに関するものである．
熱移動を伴う流れ場を有する場合の，放熱特性や流体の圧力損失を算出することは，たとえば電気機器等の電子部品の実装密度を高め，省スペースに設計する場合などに必要である．これらを実験的に求めることが一般に行われているが，試作等を要するため多くの時間とコストがかかるという問題がある．
熱移動を伴う流れ場を，シュミレーションし解析的に解く為には，連続の方程式，運動の方程式であるナビエストークス方程式及びエネルギー式を解く必要がある． 流れが層流である場合は，流れの現象をあらわす運動方程式等は乱流の場合に比較して簡単であり，方程式を解くことにより計算できる場合が多いが，実際の機器においては，多くの場合流れは乱流となっており，解析的に解を求めることはできないのが実情である．
乱流の場合にこれらの方程式を解くことを困難にしているのは，レイノルズ平均化されたナビエストークス方程式のなかにあるレイノルズ応力項と，同じくレイノルズ平均化されたエネルギー式の中にある乱流熱流束項である．レイノルズ応力項は速度変動部分の積の平均であり，乱流熱流束項は速度変動と温度変動の積の平均であるから，乱流現象を厳密に小変動にいたるまで計算しなければ算出できない．厳密に解析する方法として，ＤＮＳ（Direct Numerical Simulation）があるが，計算できるのはレイノルズ数の小さな（つまり寸法と，速度の積が小さな）場合に限られ，機器設計等での実際の流れ場においては，
計算すべき格子点数が膨大になりコンピューターの計算容量から計算は実質的には不可能である．
そこで，このレイノルズ応力項や乱流熱流束項を簡単な式に代表させ，方程式を解くことが行われている．この簡単な式をモデル式と言っている．従来モデルの例としてMKCOモデル式がある．以下にMKCOモデル式を示す．レイノルズ応力項は次のように表される．
熱移動を伴う流れ場を，シュミレーションし解析的に解く為には，連続の方程式，運動の方程式であるナビエストークス方程式及びエネルギー式を解く必要がある． 流れが層流である場合は，流れの現象をあらわす運動方程式等は乱流の場合に比較して簡単であり，方程式を解くことにより計算できる場合が多いが，実際の機器においては，多くの場合流れは乱流となっており，解析的に解を求めることはできないのが実情である．
乱流の場合にこれらの方程式を解くことを困難にしているのは，レイノルズ平均化されたナビエストークス方程式のなかにあるレイノルズ応力項と，同じくレイノルズ平均化されたエネルギー式の中にある乱流熱流束項である．レイノルズ応力項は速度変動部分の積の平均であり，乱流熱流束項は速度変動と温度変動の積の平均であるから，乱流現象を厳密に小変動にいたるまで計算しなければ算出できない．厳密に解析する方法として，ＤＮＳ（Direct Numerical Simulation）があるが，計算できるのはレイノルズ数の小さな（つまり寸法と，速度の積が小さな）場合に限られ，機器設計等での実際の流れ場においては，
計算すべき格子点数が膨大になりコンピューターの計算容量から計算は実質的には不可能である．
そこで，このレイノルズ応力項や乱流熱流束項を簡単な式に代表させ，方程式を解くことが行われている．この簡単な式をモデル式と言っている．従来モデルの例としてMKCOモデル式がある．以下にMKCOモデル式を示す．レイノルズ応力項は次のように表される．
ここで，
は，速度場における浮力のダンピング関数であって，次式で表される．
乱流熱流束項は次式で表される．
ここで，
は，速度場における浮力のダンピング関数であって，次式で表される．
以上のMKCOモデル式の記号の意味は以下のとおりである．
はそれぞれ乱流エネルギーkの速度平均勾配における生成項，浮力による生成項，散逸率である．大文字Cで表される項は係数であり，式または値で与えられる．その他の記号は他の式の場合とあわせ，この後に記載する．
このモデル式は浮力を伴う乱流を計算するのに優れたモデル式とされているが，次の欠点がある．
1)乱流熱流束項
このモデル式は浮力を伴う乱流を計算するのに優れたモデル式とされているが，次の欠点がある．
1)乱流熱流束項
が平均温度勾配
に比例するため，平均温度勾配がほとんどない場合に乱流熱流束が計算上ゼロになるが，実際はゼロでない．つまり，現象を正しく表していないため，この項が影響する場合の計算精度が悪くなる．
2)式が場合分けによって選択する形になっているため，座標依存性があり，座標軸と重力方向によって，同一の式で計算ができない．つまり重力方向に対し流路が斜め形状の場合などでは計算が困難である．
村上周三著 「CFDによる環境設計工学」東京大学出版 2000年153〜165―ページ
2)式が場合分けによって選択する形になっているため，座標依存性があり，座標軸と重力方向によって，同一の式で計算ができない．つまり重力方向に対し流路が斜め形状の場合などでは計算が困難である．
村上周三著 「CFDによる環境設計工学」東京大学出版 2000年153〜165―ページ
本発明の目的は，機械の放熱設計等に必要な，ヌッセルト数や摩擦係数などを求めるための，熱移動を伴う流れ場の乱流の流体的かつ熱的性質を，高精度で，重力方向がいかなる場合でも，同一の式で予測できる推定方法及び推定プログラムを提供することである．
請求項１記載の発明は，熱移動を伴う流れ場についての乱流の2方程式モデルにおいて，レイノルズ応力
を，強制対流に関する項をF，浮力に関する項をGとすると
の式を用いてレイノルズ応力を算出して，速度場の計算をし，乱流熱流束を
，乱流熱流束の強制対流に関する項をFt，乱流熱流束の浮力に関する項をGtとすると，
の式を用いて乱流熱流束を算出して，上記速度場における温度場の計算をすることを特徴とする，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定方法である． 本発明に従えば，レイノルズ応力を強制対流に関する項Fと浮力に関する項Gの和で表し，乱流熱流束を強制対流に関する項Ftと浮力に関する項Gtの和で表しているので，乱流の諸量を正確に予測することができ，かつこれらの式が重力の方向にかかわらず統一された式なのでいかなる重力方向であっても，同一式で計算し予測することが可能となる．
請求項２に記載の本発明は熱移動を伴う流れ場についての乱流の2方程式モデルを用いて，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定するコンピュータを，流れ場の温度，速度等を定める条件となる必要なデータを入力するための入力手段，入力手段によって入力されたデータに基づいて
請求項２に記載の本発明は熱移動を伴う流れ場についての乱流の2方程式モデルを用いて，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定するコンピュータを，流れ場の温度，速度等を定める条件となる必要なデータを入力するための入力手段，入力手段によって入力されたデータに基づいて
の式を用いてレイノルズ応力を算出して，速度場の計算を行い，
の式を用いて乱流熱流束を算出して，上記速度場における温度場の計算をすることを特徴とする，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値を求める演算手段，および演算手段によって求められた速度分布，温度分布，摩擦係数，ヌッセルト数，レイノルズ応力，乱流熱流束等の特性値を表示する出力手段，として機能させるコンピュータによって実行可能な熱移動を伴う流れ場の速度分布，温度分布，摩擦係数，ヌッセルト数，レイノルズ応力，乱流熱流束等の特性値の推定プログラムである．本発明に従えば，コンピュータにより上記乱流諸元の推定プログラムが実行され、入力手段により，乱流諸元を生成するために必要な，たとえば，流れ場の形状，流体的及び熱的条件を特定するための諸元，格子点などの各種データが入力されると，演算手段は前記入力手段から入力されたデータに基づいて
の式を使って，速度場支配方程式と温度場支配方程式を解き，速度及び温度等の乱流諸元を演算し，出力手段にその演算結果を表示することができる．このようにして，出力手段に乱流の速度及び温度に関するの諸元が表示されるので，乱流の予測ツールをコンピュータ上に実現することができる．
請求項１の本発明に従えば，レイノルズ応力を強制対流に関する項Fと浮力に関する項Gの和で表し乱流熱流速を強制対流に関する項Ftと浮力に関する項Gtの和で表しているので乱流の諸量を正確に予測することができ，かつこれらの式が重力の方向にかかわらず統一された式なのでいかなる重力方向であっても、同一式で計算し予測することが可能となる．
請求項２の本発明に従えば，コンピュータにより上記乱流諸元の推定プログラムが実行され，手段により，乱流諸元を生成するために必要な，例えば，流れ場の形状，流体的及び熱的条件を特定するための諸元，格子点などの各種データが入力されると，演算手段は前記入力手段から入力されたデータに基づいて
請求項２の本発明に従えば，コンピュータにより上記乱流諸元の推定プログラムが実行され，手段により，乱流諸元を生成するために必要な，例えば，流れ場の形状，流体的及び熱的条件を特定するための諸元，格子点などの各種データが入力されると，演算手段は前記入力手段から入力されたデータに基づいて
の式を使って，速度場支配方程式と温度場支配方程式を解き，速度及び温度等の乱流諸元を演算し，出力手段にその演算結果を表示することができる．このようにして，出力手段に乱流の速度及び温度に関するの諸元が表示されるので，乱流の予測ツールをコンピュータ上に実現することができる．
熱移動を伴う乱流の流体的及び熱的現象は周知のレイノルズ平均化された以下の式
連続の式
連続の式
ナビエストークスの式
及びエネルギー式
で表される．
ただし，式には，周知のBoussinesq近似が施されている．以上の式で，繰り返される添え字はテンソルの総和規約に従う．実際のケースでは，主流の温度
ただし，式には，周知のBoussinesq近似が施されている．以上の式で，繰り返される添え字はテンソルの総和規約に従う．実際のケースでは，主流の温度
を与条件として与えると，重力加速度
，や流体の物性値である流体の密度
，粘性係数
，熱拡散係数
は既知であり，レイノルズ応力
と乱流熱流束
の値を算出できると，未知数が平均圧力
(1つ)，平均速度ベクトル
(3方向で3つ)，平均温度
(１つ)の合計6つであるのに対し，式の数が連続の式（１つ），ナビエストークス式（3方向で3つ），エネルギー式(１つ）で合計6つとなり式を解くことができる．
本発明の特色とするところはこのレイノルズ応力と乱流熱流束を算出する式を物理現象を最もよく表すと見られる次の式としたところにある．レイノルズ応力を
本発明の特色とするところはこのレイノルズ応力と乱流熱流束を算出する式を物理現象を最もよく表すと見られる次の式としたところにある．レイノルズ応力を
，強制対流に関する項をF，浮力に関する項をGとすると
の式を用いてレイノルズ応力
を算出する．また，乱流熱流束を
，乱流熱流束の強制対流に関する項をFt，乱流熱流束の浮力に関する項をGtとすると，
の式を用いて乱流熱流束
を算出する．
以上の式の計算に必要な乱流エネルギー
以上の式の計算に必要な乱流エネルギー
，乱流エネルギー散逸率
，温度乱れ強度
，温度乱れ強度散逸率
の輸送方程式はそれぞれ次式で与えられる．
ここで，
は乱流エネルギーの生成項，
は浮力項，
は温度乱れ強度の生成項である．
の乱流拡散項については，一般化されたこう配拡散近似に基づき，以下のようにモデル化されている．
式中のモデル関数，およびモデル定数を以下にまとめて示す．
ここで，
は無次元距離
（局所レイノルズ数でもある）と乱流レイノルズ数
の調和平均形で定義された修正レイノルズ数である．次に式中におけるモデル関数，およびモデル定数を以下にまとめて示す．
ここで用いられている記号の意味するところは以下の通りである．
以下に具体的方法について説明する．
図１は本発明の実施の一形態である乱流の熱的・流体的諸元を求めるフローチャートである．図２は乱流の熱的・流体的諸元の推定プログラムを実行して，シュミレーションするために用いられるコンピュータ１の電気的構成を示すブロック図である．ステップa1で，コンピュータ１により当該プログラムが実行され，計算すべき流体，形状，温度等の与条件が入力手段２により入力されて初期設定される．乱流の熱的・流体的諸元の推定プログラムは記憶手段３に記憶され演算処理装置４からの実行指令によって読み出されて実行し，出力手段である表示手段５に計算結果を表示させ，乱流の熱的・流体的諸元の推定システムを構築する．
計算は流路を多数の格子状の微小体積に分割し，各微小体積で上記の各式が成り立つように未知の項を求める方式で行う．上記の各式は，実質微分D/Dtや偏微分(
図１は本発明の実施の一形態である乱流の熱的・流体的諸元を求めるフローチャートである．図２は乱流の熱的・流体的諸元の推定プログラムを実行して，シュミレーションするために用いられるコンピュータ１の電気的構成を示すブロック図である．ステップa1で，コンピュータ１により当該プログラムが実行され，計算すべき流体，形状，温度等の与条件が入力手段２により入力されて初期設定される．乱流の熱的・流体的諸元の推定プログラムは記憶手段３に記憶され演算処理装置４からの実行指令によって読み出されて実行し，出力手段である表示手段５に計算結果を表示させ，乱流の熱的・流体的諸元の推定システムを構築する．
計算は流路を多数の格子状の微小体積に分割し，各微小体積で上記の各式が成り立つように未知の項を求める方式で行う．上記の各式は，実質微分D/Dtや偏微分(
など)を含んでいるのでこのままでは解けないのでこれらを何らかの形で近似化する必要があり，このような近似化の手法として差分法，有限要素法，有限体積法などがある．たとえば，有限体積法では各微小部分についての積分と直線近似を用いて離散化し，多元一次方程式に置き換えて計算する．これらの手法は周知である．平均圧力
は連続の式とナビエストークスの式から周知のＳＩＭＰＬＥ解法等で求める．この手法も周知である．
DNSによる厳密な解との比較のために実施した例の場合について説明する．
図3に示す本形態は，一方が高温で，他方が低温である２枚の垂直平板間の乱流流れについてその熱的・流体的諸元を推定する方法について説明する．計算式は三次元の場合を記載するが，計算はＤＮＳ計算のできる，二次元の完全発達流について行うので，実質的には一次元である．計算格子は流れに直角方向に145点とした．計算手法は有限体積法によった．熱的流体的予条件はレイノルズ数150，グラスホッフ数9.6×10■である．
ステップａ１で初期設定すると，ステップａ２で，入力された諸条件に基づいて，以上に述べた各方程式を解き乱流の解を求める計算が開始される．
まず，ステップａ２１でレイノルズ応力
DNSによる厳密な解との比較のために実施した例の場合について説明する．
図3に示す本形態は，一方が高温で，他方が低温である２枚の垂直平板間の乱流流れについてその熱的・流体的諸元を推定する方法について説明する．計算式は三次元の場合を記載するが，計算はＤＮＳ計算のできる，二次元の完全発達流について行うので，実質的には一次元である．計算格子は流れに直角方向に145点とした．計算手法は有限体積法によった．熱的流体的予条件はレイノルズ数150，グラスホッフ数9.6×10■である．
ステップａ１で初期設定すると，ステップａ２で，入力された諸条件に基づいて，以上に述べた各方程式を解き乱流の解を求める計算が開始される．
まず，ステップａ２１でレイノルズ応力
を計算する．計算に必要な，式の中に現れる各項は，繰り返し算の１回目は，適当に仮定した値を用い，繰り返し算の２回目以降はそれ以前に求めた値を使うことにより計算される．計算は格子の各点について行われる．レイノルズ応力
が求まると，次に，ステップａ２２でナビエストークス式を計算し，平均速度ベクトル
を求める．次にステップａ２３で乱流エネルギー
，乱流エネルギー散逸率
を前述の輸送方程式を同様に有限体積法で解くことにより求める．さらに以上の計算の結果を使い，平均圧力を求め，その他の以降の計算のために必要な項の値を求める．ステップａ２４で乱流熱流束
の計算を既述の式を使い有限体積法で解き，その結果を用いてステップａ２５エネルギー式を同じく有限体積法で解き平均温度
を得る．次にステップａ２６で温度乱れ強度
，温度乱れ強度散逸率
の輸送方程式を解きこれらの値を求める．ステップａ３で以上のようにして求めた平均速度ベクトル
，乱流エネルギー
，乱流エネルギー散逸率
，平均温度
，温度乱れ強度
，温度乱れ強度散逸率
らの値が収束するまで，ステップａ２を繰り返す．収束すれば，乱流の熱的・流体的状態が求められる．以上の結果を用いて，工学的に必要な数値をステップａ４で求めることができる．たとえば流体の圧力損失係数は壁近傍の平均速度分布から求めることができる．壁近傍の温度分布からヌッセルト数を求めることができる．
図４は本計算によって得られたレイノルズ応力の一成分
図４は本計算によって得られたレイノルズ応力の一成分
を実線で示し，厳密解であるＤＮＳ値を丸印で示し，従来のＭＫＣＯ法によるものを破線で比較している．図５は乱流熱流束
について本計算結果を実線で示し，厳密解であるＤＮＳ値を丸印で示し，従来のＭＫＣＯ法によるものを破線で比較している．
レイノルズ応力について，差は少ないが，乱流熱流束について，従来は算出できなかったのが，できるようになり，乱流の構造の推定が，より詳細にできるようになったことがわかる．
レイノルズ応力について，差は少ないが，乱流熱流束について，従来は算出できなかったのが，できるようになり，乱流の構造の推定が，より詳細にできるようになったことがわかる．
以上説明したように，本発明によるレイノルズ応力と、乱流熱流速のモデル式は，乱流の物理現象を従来のモデルより正確に表現し，しかも流れと重力の方向にかかわらず同じ式で計算できるので，機械の放熱設計などにおいて，実験の手間を削減し，より予測精度を上げ，摩擦係数や，ヌッセルト数を計算し，圧力損失や放熱量の正確な予測を可能とし，速やかにコンパクトで効率のよい機械を設計する手段を提供するものである．
１ コンピュータ
２ 入力手段
３ 記憶手段
４ 演算処理装置
５ 表示手段
２ 入力手段
３ 記憶手段
４ 演算処理装置
５ 表示手段
Claims (2)
- 熱移動を伴う流れ場についての乱流の2方程式モデルにおいて，
レイノルズ応力を
，強制対流に関する項をF，浮力に関する項をGとすると
の式を用いてレイノルズ応力を算出して，速度場の計算をし，乱流熱流束
を，乱流熱流束の強制対流に関する項をFt，乱流熱流束の浮力に関する項をGtとすると，
の式を用いて乱流熱流束を算出して，上記速度場における温度場の計算をすることを特徴とする，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定方法 - 熱移動を伴う流れ場についての乱流の2方程式モデルを用いて，
乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定するコンピュータを，流れ場の温度，速度等を定める条件となる必要なデータを入力するための入力手段，入力手段によって入力されたデータに基づいてレイノルズ応力を
，強制対流に関する項をF，浮力に関する項をGとすると
の式を用いてレイノルズ応力を算出して，速度場の計算を行い,乱流熱流束を
，乱流熱流束の強制対流に関する項をFt，乱流熱流束の浮力に関する項をGtとすると，
の式を用いて乱流熱流束を算出して，上記速度場における温度場の計算をすることを特徴とする，乱流流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値を求める演算手段，および演算手段によって求められた速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値を表示する出力手段，として機能させるコンピュータによって実行可能な熱移動を伴う流れ場の速度分布，温度分布，レイノルズ応力，乱流熱流束等の特性値の推定プログラム
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005078548A JP4797157B2 (ja) | 2005-03-18 | 2005-03-18 | 浮力を伴う乱流の流体的及び熱的諸特性の推定方法及び推定プログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005078548A JP4797157B2 (ja) | 2005-03-18 | 2005-03-18 | 浮力を伴う乱流の流体的及び熱的諸特性の推定方法及び推定プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2006258391A true JP2006258391A (ja) | 2006-09-28 |
| JP4797157B2 JP4797157B2 (ja) | 2011-10-19 |
Family
ID=37097847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005078548A Active JP4797157B2 (ja) | 2005-03-18 | 2005-03-18 | 浮力を伴う乱流の流体的及び熱的諸特性の推定方法及び推定プログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4797157B2 (ja) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008129953A (ja) * | 2006-11-22 | 2008-06-05 | Denso Techno Corp | 管路内流れの圧力損失の評価方法及び装置 |
| CN102621094A (zh) * | 2012-03-07 | 2012-08-01 | 北京大学 | 一种热对流斑图测试仪及其测量方法 |
| CN106897537A (zh) * | 2017-03-14 | 2017-06-27 | 清华大学 | 含三维或曲面外形结构体的温度场与热流同时重构方法 |
| CN109948170A (zh) * | 2017-12-20 | 2019-06-28 | 中国科学院长春光学精密机械与物理研究所 | 一种光学遥感相机热稳定性的仿真方法 |
| CN111222293A (zh) * | 2020-01-10 | 2020-06-02 | 大连交通大学 | 基于sngr方法的高速列车气动噪声计算方法 |
| CN113899657A (zh) * | 2021-12-06 | 2022-01-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于复合式纹影技术的可压缩雷诺应力测量方法 |
| CN113899525A (zh) * | 2021-12-06 | 2022-01-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于复合式纹影技术的可压缩雷诺应力测量系统 |
| CN114754717A (zh) * | 2022-03-21 | 2022-07-15 | 天津大学 | 一种基于热学原理测量冰层厚度的方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07146194A (ja) * | 1993-09-28 | 1995-06-06 | Nippon Sheet Glass Co Ltd | 流体中の圧力分布を求める方法 |
-
2005
- 2005-03-18 JP JP2005078548A patent/JP4797157B2/ja active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07146194A (ja) * | 1993-09-28 | 1995-06-06 | Nippon Sheet Glass Co Ltd | 流体中の圧力分布を求める方法 |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008129953A (ja) * | 2006-11-22 | 2008-06-05 | Denso Techno Corp | 管路内流れの圧力損失の評価方法及び装置 |
| CN102621094A (zh) * | 2012-03-07 | 2012-08-01 | 北京大学 | 一种热对流斑图测试仪及其测量方法 |
| CN102621094B (zh) * | 2012-03-07 | 2013-11-27 | 北京大学 | 一种热对流斑图测试仪及其测量方法 |
| CN106897537A (zh) * | 2017-03-14 | 2017-06-27 | 清华大学 | 含三维或曲面外形结构体的温度场与热流同时重构方法 |
| CN109948170A (zh) * | 2017-12-20 | 2019-06-28 | 中国科学院长春光学精密机械与物理研究所 | 一种光学遥感相机热稳定性的仿真方法 |
| CN111222293A (zh) * | 2020-01-10 | 2020-06-02 | 大连交通大学 | 基于sngr方法的高速列车气动噪声计算方法 |
| CN111222293B (zh) * | 2020-01-10 | 2023-07-25 | 大连交通大学 | 基于sngr方法的高速列车气动噪声计算方法 |
| CN113899657A (zh) * | 2021-12-06 | 2022-01-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于复合式纹影技术的可压缩雷诺应力测量方法 |
| CN113899525A (zh) * | 2021-12-06 | 2022-01-07 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种基于复合式纹影技术的可压缩雷诺应力测量系统 |
| CN113899657B (zh) * | 2021-12-06 | 2022-02-18 | 中国空气动力研究与发展中心高速空气动力研究所 | 基于复合式纹影技术的可压缩雷诺应力测量方法 |
| CN114754717A (zh) * | 2022-03-21 | 2022-07-15 | 天津大学 | 一种基于热学原理测量冰层厚度的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4797157B2 (ja) | 2011-10-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4797157B2 (ja) | 浮力を伴う乱流の流体的及び熱的諸特性の推定方法及び推定プログラム | |
| Ozden et al. | Shell side CFD analysis of a small shell-and-tube heat exchanger | |
| Francois et al. | Computations of drop dynamics with the immersed boundary method, part 1: numerical algorithm and buoyancy-induced effect | |
| JP2020027658A (ja) | 安定した陽的拡散の性能及び精度の向上 | |
| Luo et al. | Flow and forced-convection characteristics of turbulent flow through parallel plates with periodic transverse ribs | |
| Niceno et al. | Unstructured large eddy and conjugate heat transfer simulations of wall-bounded flows | |
| Rozza et al. | Reduced basis methods and a posteriori error estimators for heat transfer problems | |
| Hernández-López et al. | Analysis of laminar and turbulent natural, mixed and forced convection in cavities by heatlines | |
| Ferreira et al. | A combination of implicit and adaptative upwind tools for the numerical solution of incompressible free surface flows | |
| Lacanette et al. | Numerical simulation of gas-jet wiping in steel strip galvanizing process | |
| Shahane | Numerical simulations of die casting with uncertainty quantification and optimization using neural networks | |
| Dhinsa et al. | Investigation into the performance of turbulence models for fluid flow and heat transfer phenomena in electronic applications | |
| Selimefendigil | Numerical investigation and POD-based interpolation of natural convection cooling of two heating blocks in a square cavity | |
| Belharizi et al. | RANS and LES computations of natural convection in a square cavity | |
| Zhao | Stable computation of turbulent flows with a low-Reynolds-number k-ϵ turbulence model and explicit solver | |
| Kasturi Rangan et al. | A face-based immersed boundary method for compressible flows using a uniform interpolation stencil | |
| Shahane et al. | Simulations of die casting with uncertainty quantification | |
| Zając et al. | Novel method for fast FEM simulation of chips with integrated microchannel cooling | |
| Ferreira et al. | Incompressible Turbulent Flow Simulation Using the κ‐ɛ Model and Upwind Schemes | |
| Hu et al. | Validation and Evaluation of Turbulence Models in Thermal Predictions of a Small Data Center Facility | |
| Ravishankar et al. | Enhanced cooling of electronic components using fluid flow under high adverse pressure gradient | |
| JP2005250855A (ja) | 筐体表面の温度予測装置及び方法、電子機器の製造装置及び方法、筐体表面の温度予測システム、サーバ、データベースサーバ、クライアント、筐体表面の温度予測プログラム及びこれを記録した記録媒体 | |
| Martinez Carrascal | Implementation in OpenFoam of a thermal-fluid analysis for thermal internal flows | |
| Fiala et al. | Determination of the Velocity Field for the Calculation of Wall Heat Transfer Coefficients | |
| Gontijo et al. | The numerical modeling of thermal turbulent wall flows with the classical κ-ε model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080118 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101116 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110112 |
|
| 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 Effective date: 20110705 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |