雷諾方程中湍流渦黏性的Richardson勢算子本構模型

摘要： 證明三維湍流Richardson擴散方程中的變系數拉普拉斯算子是基本解為r-7/3的一個勢算子，并稱其為Richardson勢算子，其中，r是空間中兩點的歐幾里得距離.基于Kolmogorov標度律和湍流快擴散（superdiffusion）理論，運用隱式微積分方程建模方法提出雷諾方程中的Richardson勢算子湍流渦黏性本構方程.

關鍵詞： 湍流； 渦黏性本構； Richardson擴散方程； 基本解； Richardson勢算子； 隱式微積分方程建模

中圖分類號： O39； O241.8文獻標志碼： A

Abstract： It is verified that the Laplacian operator with varying coefficient in the 3D turbulence Richardson diffusion equation is a potential operator having the fundamental solution r-7/3， which is called the Richardson potential operator and where r is the Euclidean distance between two points. Based on the Kolmogorov scaling law and turbulence superdiffusion theory， a Richardson potential operator equation of turbulence eddy viscosity constitutive relationship in the Reynolds equation is proposed by the implicit calculus equation modeling approach.

Key words： turbulence； eddy viscosity constitutive relationship； Richardson diffusion equation； fundamental solution； Richardson potential operator； implicit calculus equation modeling

3討論

隱式微積分建模方法將微積分建模與統(tǒng)計模型深刻緊密地結合起來，揭示確定性模型與隨機模型的內在聯系.本文利用隱式微積分建模方法發(fā)展的雷諾湍流模型，是一個統(tǒng)計意義清晰的確定性物理模型.

Richardson勢算子基本解r-7/3與湍流的Kolmogorov標度率有數學、力學上的清晰聯系.另一方面，其也許與湍流的分形結構有關，深入分析這個問題是下一步工作的重點.

類似于文獻[2]的方法，也可以用Richardson勢算子構造間歇性湍流的統(tǒng)計方程Pt-γΔRP-υΔP=0（11）式中：P為概率密度函數.式（11）包含渦尺度和分子尺度的黏性擴散行為，可刻畫湍流的多尺度行為.

文獻[2]提出的分數階拉普拉斯算子渦黏性本構模型，本質上假設湍流渦的擴散服從Lévy穩(wěn)態(tài)分布[3]，但是湍流實驗數據更接近伸展高斯分布[6].本文引入的Richardson勢算子湍流渦黏性本構模型在統(tǒng)計上反映湍流的伸展高斯分布特征.此外，Lévy穩(wěn)態(tài)分布的2階矩無窮大[3]，而伸展高斯分布沒有這個問題.有關這兩個模型的數值驗證是一個非常重要的工作，可以考慮從槽道湍流問題入手.

參考文獻：

[1]BATCHELOR G K. Theory of homogeneous turbulence[M]. Cambridge： Cambridge University Press， 1953.

[2]CHEN W. A speculative study of 2/3order fractional Laplacian modeling of turbulence： Some thoughts and conjectures[J]. Chaos： An Interdisciplinary Journal of Nonlinear Science， 2006， 16（2）： 023126. DOI： 10.1063/1.2208452.

[3]SAICHEV A， ZASLAVSKY G M. Fractional kinetic equations： Solutions and applications[J]. Chaos： An Interdisciplinary Journal of Nonlinear Science， 1997， 7（4）： 753764. DOI： 10.1063/1.166272.

[4]JULLIEN M C， PARET J， TABELING P. Richardson pair dispersion in twodimensional turbulence[J]. Physical Review Letters， 1999， 82（14）： 28722875. DOI： 10.1103/PhysRevLett.82.2872.

[5]MAJDA A J， KRAMER P R. Simplified models for turbulent diffusion： Theory， numerical modeling， and physical phenomena[J]. Physics Reports， 1999， 314（45）： 237574.

[6]PORTA A L， VOTH G A， CRAWFORD A M， et al. Fluid particle accelerations in fully developed turbulence[J]. Nature， 2001， 409： 10171019.

[7]SOKOLOV I M， KLAFTER J， BLUMEN A. Ballistic versus diffusive pair dispersion in the Richardson regime[J]. Physical Review E， 2000， 61（3）： 27172722. DOI：10.1103/PhysRevE.61.2717.

[8]陳文.復雜科學與工程問題仿真的隱式微積分建模[J]. 計算機輔助工程， 2014， 23（5）： 16.

CHEN W. Implicit calculus modeling for simulating complex scientific and engineering problems[J]. Computer Aided Engineering， 2014， 23（5）： 16. DOI： 10.13340/j.cae.2014.05.001.

（編輯于杰）