製品情報

Core Package

Advanced Package

Two-phase Package
二相流・混相流

CVD Package

Rhino CFD

事例

必要環境

導入費用

導入実績

PHOENICSの歴史

流体解析事例コラム

お客様インタビュー

流体解析の基礎技術情報

PHOENICS機能一覧

トップページPHOENICSと熱流体解析 > 流体解析の歴史的背景(2)

◎流体解析の歴史的背景(4)

【第10回】スポルディングの生涯概要

それで、Brian SpaldingはCFDを発明しましたか?その答えは、明らかに無条件の「はい」ではありません。CFDの起源は1928年またはそれ以前にさかのぼることができます。ロスアラモスのHarlow氏、学術研究機関の他の人たちも、後にCFDと呼ばれた分野に大きく貢献していました。彼らはすでに、ICグループによって使用されたりまたは再発明されたりする主要なコンセプトの多くを検討していました。Spaldingはその名前さえ付けませんでした。CFDの名前はPat Roache [1972]によって彼の有名な本で初めて使われました。しかし、別の質問をしましょう。Spaldingなしで今日知っているようなCFDがありますか?その答えは無条件ではありません。誰よりも、エンジニアにとって興味ある問題にCFDを適用する、エンジニアリングツールとしてのCFDを作成しました。今日の市販されているCFDソフトウェアツールのほとんどは、Spaldingと彼のグループが1965-1975年に渡って行った作業に基づいています。

彼は能動的リーダーであり、効率的でエンジニアに優しいCFD方法論の開発につながった斬新なアイデアの重要な貢献者でした。この10年間の彼の重要な洞察は、しばしばブレークスルーを可能にし、重大な瞬間に焦点を絞った。彼は100人以上の大学院生を抱えていました。これらの学生の多くは、自分自身で重要な専門的貢献をしました。LaunderとWhitelawの両方がSpaldingと緊密に連携していた1965年から1972年までの期間を考えると、「IC」手法のトーチを運ぶ学生達が、一人の厳密な監督下でグループ化されたことは真に驚くべきものです。これらの学生の多くは、世界の舞台で独特の印をつくっています。今日、CFD、乱流、燃焼のよく知られた名前には、第1世代または第2世代のICとつながっています。

任意のベクトルが回転(ベクトル)および非回転(スカラー)成分として表せることは古くから知られています。その圧力はスカラー成分に関係しており、ベクトルとしての渦は古典的な流体力学の一部です。HarlowとChorinは、速度補正をベクトルと任意のスカラーステージに分解することを提案してきました。Harlowのグループは、連続方程式を圧力の方程式に変えるためにそれを使用しました。しかし、非定常流れの数値計算に必要なので、定常的な流れ(エンジニアにとって日常的な関心事)のための非効率的な仕組みでした。PatankarとSpalding(1972)の研究は、経済的且つ実用的の手法をツールに組み込んで、定常流れにも簡単に使えるようにしました。だから、圧力投影法ではSpaldingに功を奏してはいけませんが、定常流れのための実用的な代替手段として使用しています。

別の例を挙げる。有限差分は長年にわたって存在してきました。数学者はすでに離散空間を試していたとBanach、Hilbertは記述しました。Thom [1928]、von Neumann [1944]、Southwell [1946]、Courant et al[1952]、およびその他数多くの研究者達が流体解析に適用しました。Spaldingは同じツールから始まり、楕円形から双曲線に及ぶ非線形方程式を扱う際に、FDとTaylor Seriesのデメリットに即座に認識しました。彼は書かれていた数式と実際の物理学の関係を見抜くことができました。「物理量」(フラックス、熱など)が保存されたやり方で動くのです。これは最終的にFVのアプローチにつながりました。ICグループだけがFVアプローチを使用したことではないです。Harlowのグループも本質的にFVの観点から考えており、フラックスと保存された量の式で表現しました。しかし、彼らは「タンクとチューブ」のような単純な視覚的な類推によってそれを表現したことはありません。彼らは、出版物に対してもコントロールボリュームで統一することはしませんでした。 Edwards [1968]は、1960年代半ばに始まったローレンス・リヴァーモアのTRUMPというコードで、FVの考えと非構造格子アプローチを使用していました。しかし、私は大部分のCFD研究者は、今日でもエドワーズの仕事を知らないと考えています。私は確かに1978年まで彼のことを聞いていません。ちょうど「巻き上げ」のように、FVはおそらく再発明されました。Spaldingは、この技術を広く使用したことは間違いありません。


多かれ少なかれ同じ話がk-εモデルで繰り返される。基礎方程式はKolmogorov [1942]とPrandlt [1945]によって書かれました。Rotta [1951]は1950年代の知識に加えて、長さの尺度と頻度の概念について詳述した。Davidov(1961)は、非常に複雑な3次モーメント閉鎖の一環として正式にεの閉形式を導出し、提案したようです。Harlow and Nakayama [1968]はε方程式を開発し、いくつかの流れに対して2方程式渦粘性法でその有効性を実証しました。Launderの博士コースの学生であったHanjalicは、k-lとk-klベースの変数を使って壁付近と壁から遠く離れた両方の流れを予測する際に、ICの同僚が経験している困難を知っていました。彼はDavidovのε方程式の係数を適合させ、初めてk-ε渦粘度式で試し、同様な式が壁と自由流れの両方を模擬できることを見出した。その後、それを3方程式モデル(k-ε-)に組み込み、せん断応力とひずみが異なる場所で消失する流れを予測した(Hanjalic, 1970)。Bill Jones(Jones&Launder,1972)は、k-εモデルを拡張して、サブレイヤを壁の右に渡って適用できるようにしました(サブレイヤが普遍的でない状況を可能にします)。k-l、k-kl、k-ωモデルを探索するまでのSpaldingは、「標準」モデルの利点と将来の研究にもたらす順序を見抜きました。彼は他のモデルでの作業を中止し、ICグループは選択モデルとしてk-εを採用しました。これはLaunder and Spalding[1972]とLaunder and Spalding [1974]の画期的結果であり、工学応用のための一般的なツールとしてのk-εモデルの公式化、「標準化」と受け入れにつながりました。

Spaldingの洞察力と実用的なツールに焦点を当てたため、これらの技術を確立させ、且つ技術者にとって有用になりました。彼は注目点を数学と現象の物理特性から現場の実課題へと変更しました。私の見解では、ブライアンはCFDの「科学」を発明しなかったが、彼は依然としてエンジニアリング設計と実践のためのCFDの「芸術」と「技術」を磨いた人物であると言えるでしょう。彼でなければ、CFDは学術機関や研究機関で実践された難解な科学にとどまっているかもしれません。

ブライアンは、彼のキャリアの上で工学と科学の様々な分野で活躍されました。彼は、燃焼、乱流、熱と物質移動、CFDにおいて、先駆的かつ永続的な貢献をしてきました。特にCFDだけを見ると、ブライアンの貢献(発明または再発明)は少なくとも以下が言えます:

1.有限体積法
2.一般化された輸送方程式により、すべての二次対流拡散系を統一
3.風上的数値処理
4.圧力プロジェクション
5.スタッガード格子
6.K-ε乱流モデル
7.壁熱伝達においての乱流エネルギーに注目
8.乱流と乱流エネルギの相互作用を考慮した渦崩壊モデル
9.多相流に対するIPSA
10.放射線に対する6-FluxモデルとIMMERSOL法

彼が初めて数値解法の文脈において、すべての二次輸送方程式を一つの一般化輸送方程式として解くことを提案しました。したがって、彼は一貫して単一の「普遍的な」方法によって流体の流れ、熱および物質輸送の処理、および機械的な応力を「統一」する彼の生涯にわたる目標をある程度達成しました。乱流、多相、固液相互作用および壁距離計算に関する彼の他の研究は、これらと同じように一般的な適応性はまだ見ていないが、将来的に実を結ぶ重要な進路研究を指摘しています。

ニュートンの言葉を借りれば、Spaldingは巨人の肩に立っていて、仲間よりも遠くが見えました。彼は、統一的な流れ、熱および物質輸送が技術者に実用的な道具につながることを予見しました。彼は、数学ではなく物理学を見れば、CFDツールの受け入れが広がることを予見していました。彼は設計ツールに変わったCFDが工学に革命を起こすと予見していました。彼はCFDに商業的機会があることを予見しました。

「統一」と「大胆な洞察」を伴う他人のアイデアの「採用」によるパターンは、彼の生涯の共通するテーマでああります。Spaldingの思考の中で共通する特徴的なスレッドを特徴づけるなら、「数学」ではなく「物理学」であることでしょう。彼の主な貢献は複雑な数学的概念を伴うものですが、彼はその概念の背後にある物理学を見ていたので、彼のブレークスルーが起こりました。彼は数学方程式の森で迷子になるのを拒否し、これらの方程式が表現しようとしていたプロセスを見ていたので、「木」を見ることができました。彼はエンジニアのように思いました。彼は繰り返し「私はこの問題を解決するのに役立つのか」という質問を繰り返しました。これは何よりも彼の重要な貢献です。彼は、2相、乱流、境界条件などに関する私たちの知識の中で「ギャップ」を抱いていました。彼が研究開発の上で「障害」が出会った時は、いつもエンジニアの解決策を見つけました。それがうまくいったらたとえ理論的には多少厳密でないところが有ったとしても、それは仕事を終わらせるのに十分でした。CFDで燃焼を試みることができると思っていた人はいませんでしたので、私たちの知る限りでは多くのギャップがあったとき、彼は燃焼のための渦巻き破砕のようなものを提案しました!

ブライアンは私にとってリチャード・ファインマンの有名な言葉の最高の例です:「あなたは何を気にしていますか?」彼は、彼が最も顕著な特徴は「彼は彼の仕事をした。何がもっと良い方法か知らなかったため」と感じています。彼はVirgil:「Possunt quia possunt videntur」を引用しています(彼らはできると思うのでそうすることができます)。彼はそれを「自分ができないことを知らないのでできます」と変えたがっていました。他の人が既に問題を宣言していて、解決できない、またはあまりにも嫌なことを知らされていないことを単に知らなかった(あるいは、おそらく気にしなかった)人は、数学的および物理的な議論によってそれを証明しようとしていました。彼は単に問題の本質に対処するための実用的な解決法を考案するだけでした。

ブライアンは、問題の本質を理解し、解決策を提案する巨大な知的能力を持っています。私が苦労している問題を彼に相談すると、問題の本質を見ることができることを何度も何度も経験しました。ブライアンは新しいアイデアを生み出すだけでなく、彼と一緒に働く他の人たちに創造性を引き出す素晴らしい能力を持っています。彼の優勢な知性と鋭い洞察力のほかに、眼の前の利益に混乱されない、「平常心」を持っているので、彼は即座の関心事に焦点を当て、彼の巨大な業績に繋がりました。ブライアンがある問題に集中したとき、彼はその目的に直接関係のないすべてのものを無視することができます。私は、Brianが世界を二つのグループに分けて見ていると言うことを敢然と言います:彼と共に働く人と、彼の目的に関係のない人に分けます。 ここで彼は彼の85歳の誕生日を迎います。私は彼と一緒にいた彼の姿を思い出します:ガラスを手にし、素敵な女の子の話を注意深く聞いている姿を。。。
ハッピーバースデー、ブライアン。



◎参考文献

D. Agonafer, Liao Gan-Li, D.B. Spalding [1996], The LVEL Turbulence Model For Conjugate Heat Transfer At Low Reynolds Numbers, ASME International Mechanical Congress and Exposition, Atlanta.
H.Z. Barakat and J.A. Clark [1965], Transient Natural Convection in Flows in Closed Containers, Univ. of Michigan, Mech. Eng. Dept., Heat Transfer Lab, Tech Rep No. 2.
J.P. Botha and D.B. Spalding [1954], The Laminar Flame Speed Of Propane-Air ixtures With Heat-Extraction From The Flame. Proc Roy Soc, A, Vol. 225, pp 71-96.
O.R. Burggraf [1966], Analytical and Numerical Studies of Steady Separated Flows, J. Fluid Mechanics, Vol. 24, 1, pp 113-151.
L.S. Caretto, R.M. Curr and D.B. Spalding [1972], Two Numerical Methods for Three-Dimensional Boundary Layers, Computer Methods in Applied Mechanics & Engg.,, Vol. 1, 1, pp. 39-57, June.
A.J. Chorin [1967], The Numerical Solution of the Navier-Stokes Equations for an Incompressible Fluid, Bull. Am. Math. Soc., Vol. 73, 6, pp 928-931.
P.Y. Chou [1945], On the Velocity Correlations and the Solution of Equations of Turbulence Fluctuations, Quart. Appl. Math., Vol. 3, p 38.
R. Courant, E. Isaacson and M. Rees [1952], On the Solution of Non-Linear Hyperbolic Differential Equations by Finite-Differences, Comm. Pure Applied Math, Vol. 5, pp 243-255.
B.I. Davidov.[1961], On the Statistical Dynamics of an Incompressible Turbulent Fluid, Dokl. AN SSR, 136 #1, 47-50.
E.R.G. Eckert [1949], V Lieblein, Forschung Vol. 16, p 33.
A.L. Edwards [1968] TRUMP: A Computer Program for Transient and Steady-State Temperature Distributions in Multi-Dimensional Systems, Lawrence Radiation Laboratory, Univ. of California, Livermore, California, UCRL 14754.
M.P. Escudier [1967], The Turbulent Incompressible Hydrodynamic Boundary Layer, Ph.D. Thesis, Imperial College, London University, U.K.
M.P. Escudier and W.B. Nicoll [1966], The Entrainment Function in Turbulent- Boundary-Layer and Wall-Jet Calculations, J. Fluid Mechanics, Vol. 25, pp. 337-366.
A.D. Gosman, W.M. Pun, A.K. Runchal, D.B. Spalding, M. Wolfshtein [1969], Heat and mass Transfer in Recirculatimg Flows, Academic Press, London.
K. Hanjalic [1970], Two Dimensional Asymmetrical Turbulent Flow in Ducts, Ph.D. Thesis, Imperial College, London University, U.K.
F.H. Harlow and J.E. Welch [1965], Numerical Calculation of Time-dependent Viscous Incompressible Flow with Free-Surface, Physics of Fluids, Vol. 8, No. 12, pp 2182- 2189.
F. H. Harlow and P.I. Nakayama [1968], Transport of Turbulence Energy Decay Rate, LA-3854, Los Alamos Scientific Laboratory, University of California.
C.W. Hirt [1968], Heuristic Stability Theory for Finite-Difference Equations, J. Computational Physics, Vol. 2, 4, pp. 339-355.
C.L.V. Jayatilllaka [1966], The Influence of Prandtl Number and Surface Roughness on the Resistance of the Laminar Sub-Layer to Momentum and Heat Transfer, Ph.D. Thesis, Imperial
College of Science and Technology, London, U.K.
W. P. Jones and B.E. Launder [1972], The Prediction of Laminarization with a Two-Equation Model of Turbulence, Int J Heat & Mass Transfer, Vol 15, 301-314.
T. von Karman [1921], Z Angew. Math. Mech, Vol. 1, p 233.
A.N. Kolmogorov [1942], Equations of Motion of an Incompressible Turbulent Fluid, Izv Akad Nauk SSSR Ser Phys, VI No 1-2, p56.
G. Kruzhilin [1936], J Tekh Phys (USSR) Col 3, p 183.
B.E. Launder and D.B. Spalding [1972], Mathematical Models of Turbulence, Academic Press.
B.E. Launder and D.B. Spalding [1974], Computational Methods for Turbulent Flows, Comp. Meth in Appl Mech & Eng'g, Vol 3, pp 269-289.
B.R. Morton, G.I. Taylor and J.S. Turner [1956], Turbulent Gravitational Convection From Maintained And Instantaneous Sources, Proc Roy Soc, A, Vol. 234, pp 1-23.
S.V. Patankar [1967], Heat and Mass Transfer In Turbulent Boundary Layers, Ph.D. Thesis, Imperial College, London University, U.K.
S.V. Patankar, D.B. Spalding [1967], Heat and Mass Transfer In Boundary Layers, Morgan- Grampian, London.
S.V. Patankar, D.B. Spalding [1972], A Calculation Procedure for Heat, Mass And Momentum Transfer In Three-Dimensional Parabolic Flows; Int. J Heat Mass Transfer, Vol. 15, p 1787..
L. Prandtl [1945], Ueber Ein Neues Formelsystem fuer die Ausgebildete Turbulenz, Nachr. Akad. Wiss. Goettingen
G.D. Raithby [1976], Skew-Upstream Differencing Schemes for Nearly-Steady Problems Involving Fluid Flow, Computer Methods in Applied Mechanics and Engineering, Vol. 9, pp. 153-164.
C.M. Rhie and W.L. Chow [1983], Numerical Study of the Turbulent Flow Past an Airfoil with trailing Edge Separation, AIAA Journal, Vol. 21, No. 11, pp. 1525-1532.
P.J. Roache [1972], Computational Fluid Dynamics, Hermosa Publishers.
J.C. Rotta [1951], Statistische Theorie Nichthomogener Turbulenz, Zeitschrift fur Physik, Vol. 129, pp 547-572.
A.K. Runchal [1967], Three Finite-Difference Methods for Navier-Stokes Equations. Proc. 2nd Applied Mechanics Conference, Univ. of Strathclyde, Glasgow.
A.K. Runchal [1969], Transport Processes in Steady Two-Dimensional Separated Flows, Ph.D. Thesis, Imperial College of Science and Technology, London, U.K.
A. K. Runchal [1971], A Non-Staggered Finite-Difference Procedure for 3D Primitive Variable Navier-Stokes Equations and a Concise Notation for Difference Schemes, Technical Note: 3D-Duct Flow Group, 15 July, 1971, Mech. Eng. Dept., Imperial College, London.
A.K. Runchal [1972], Convergence and Accuracy of Three Finite-Difference Schemes for a Two- Dimensional Conduction and Convection Problem. Int'l J. Num. Methods in Engineering, Vol. 4, p. 540-550.
A.K. Runchal and M. Wolfshtein [1969], Numerical Integration Procedure for the Steady-State Navier-Stokes Equations. J. Mech. Eng. Sci., II, Vol. 5, p. 445-453.
P.G. Saffman [1976], Development of a Complete Model for the Calculation of Turbulent Shear Flows, Symposium on Turbulence and Dynamical Systems, Duke Univ., Durham, NC, April.
N.N. Semenov [1940], Progress Phys. Sci. (USSR), Vol. 24, p 433.
D.B. Spalding [1951], The Combustion of Liquid Fuels, Ph.D. Thesis, Cambridge University.
D.B. Spalding [1955], Some Fundamentals of Combustion, Butterworths, London.
D B Spalding [1957], Predicting the Laminar Flame Speed In Gases With Temperature-Explicit Reaction Rates, Combustion & Flame, Vol. 1, No 3, pp 296-307.
D.B. Spalding [1957], Analogue For High-Intensity Steady-Flow Combustion Phenomena, Proc. I. Mech. E., London, Vol. 171, No 10, pp 383-411.
D.B. Spalding [1961], A Single Formula For The Law Of The Wall, Trans ASME Series A, J. Appl. Mech., Vol. 28, No 3, pp 444-458.
D.B. Spalding [1963], Convective Mass Transfer: An Introduction, Edward Arnold (Publishers) Ltd., London
D.B. Spalding [1964], A Unified Theory Of Friction, Heat Transfer And Mass Transfer In The Turbulent Boundary Layer And Wall Jet, UK Aeronautical Research Council Report, ARC 25 925, March.
D.B. Spalding, [1971], Mixing and Chemical Reaction in Confined Turbulent Flames, 13th International Symposium on Combustion, The Combustion Institute, pp 649-657.
D.B. Spalding, [1972], A Novel Finite-Difference Formulation for Differential Expressions Involving Both First and Second Derivatives, Int. J Num, Methods Eng., Vol. 4, pp 551-559.
D.B. Spalding [1977], GENMIX: A General Computer Program For Two-Dimensional Parabolic Phenomena, Pergamon Press, Oxford.
D.B. Spalding [1994], A Simple Method of Calculating Distances from and Gaps between the Walls of Solid Bodies, Poster session. International Heat Transfer Conference, Inst. Chem. Eng. London.
D.B. Spalding [1996], PHOENICS Encyclopedia, Article on Radiation Models in PHOENICS, Section on IMMERSOL.
D.B. Spalding [1996], Multi-fluid Models of Turbulent Combustion, CTAC95 Biennial Conference, Melbourne, Australia, published by World Scientific Publishing Co, pp 59-81.
D.B. Spalding [2005], Unifying two fields of Computational Mechanics: Solid and Fluid, Invited Lecture at International Symposium on Science and Society, St Petersburg, Russia.
R.V. Southwell [1946], Relaxation Methods in Theoretical Physics, Oxford University press, New York, New York.
A. Thom [1928], An Investigation of Fluid Flow in Two Dimensions, Aerospace Research Council, United Kingdom, R & M No. 1194.
J. von Neumenn [1944], Proposal for Analysis of Numerical Method for the Treatment of Hydrodynamical Shock Problems, National Defense and Research Com., Report AM-551, March.
P. Vanka [1986], Block-Implicit Multigrid Solution of Navier-Stokes Equations in Primitive Variables, J. Computational Physics, Vol. 65, 1, pp. 138-158.
D.C. Wilcox [2006], Turbulence Modeling for CFD, 3rd Edition, DCW Industries, Inc., La Canada, CA 91011.
M. Wolfshtein [1968], Numerical Smearing in One Sided Difference Approximations to the Equations of Non Viscous Flow. EF/TN/3, Dept. of Mech. Eng., Imperial College, London.
M. Wolfshtien [1968], Turbulent Convection in Impinging Jets, Ph.D. Thesis, Imperial College, London University, London, U.K.
Y.B. Zeldovich and V.A. Frank-Kamenetskiy [1938], J. Phys. Chem. (USSR), Vol. 12, p 100.

 
©2005 Copyright Concentration Heat and Momentum Ltd. Welcome to this Web site ! since 8/Jan/1999