電力中央研究所

報告書「電力中央研究所報告」は当研究所の研究成果を取りまとめた刊行物として、昭和28年より発行されております。 一部の報告書はPDF形式で全文をダウンロードすることができます。

※ PDFのファイルサイズが大きい場合には、ダウンロードに時間がかかる場合がございます。 ダウンロードは1回のクリックで開始しますので、ダウンロードが完了するまで、複数回のクリックはなさらないようご注意願います。

電力中央研究所 報告書(電力中央研究所報告)

報告書データベース 詳細情報


報告書番号

U02042

タイトル(和文)

複雑形状流れに対する非定常乱流解析コードの開発(その1) --有限要素法LES解析コードSMART-FEMの開発--

タイトル(英文)

Development of computer codes for an unsteady turbulent flow with complex geometry. Part 1: Development of SMART-FEM code based on finite element method for large-eddy simulation

概要 (図表や脚注は「報告書全文」に掲載しております)

大渦解析(LES)の基礎方程式を厳密に導出し、空間平均化と空間偏微分の可換性の仮定が壁近傍では誤差の原因となることを明かにした。この誤差には空間平均を施さない速度と圧力が含まれるため、誤差低減のスキームをLES解析に導入することは困難と考えられる。この誤差を認識した上で、有限要素法を用いた非定常乱流解析コードをLES方程式に基づき開発した。近似関数としては複雑形状データの作成が容易な低次要素を用いたが、この要素固有の数値振動の可能性があるため、その抑制手法を導入した。本コードを平行平板内の3次元乱流に適用したところ、要素平均の流速に基づく平均流速と乱れ量は既往の結果と良く一致することが示され、本コードの妥当性を検証した。さらに、複雑形状流れである2次元十字構造物周りの流れを解析し、構造物に課せられる抗力・揚力・回転モーメントを算出し、その変動周期がほぼ妥当であることを確認した。

概要 (英文)

The exact form of the space-filtered basic equations for the large-eddy simulation (LES) was derived, which revealed that the commutativity between space-filtering and spatial partial differentiation breaks down near flow boundaries such as walls and inlets/outlets. The extra terms due to the commutativity errors involve unfiltered velocity and pressure, which could not be handled in the LES computation and then must be neglected. Understanding these unavoidable errors of the LES formulation, a finite element LES code (SMART-FEM) was developed to simulate an unsteady turbulent flow with complex geometry. A low-order element (multi-linear interpolation for velocity, and piecewise constant for pressure) was used here to have flexibility for mesh generation of complex geometries. Since this element tends to invoke spatial instability, some countermeasures are provided to suppress the numerical oscillation. The present code was verified for a 3-D parallel channel flow by comparing the averaged-velocity profile and magnitude of the velocity fluctuation with those of a direct simulation result. Furthermore, the code was applied to 2-D turbulent flow around a cruciform structure to compute the drag force, lift force and rotating moment. The primary fluctuating period turned out to be consistent with the expected values.

報告書年度

2002

発行年月

2003/04

報告者

担当氏名所属

江口 譲

我孫子研究所流体科学部

服部 康男

我孫子研究所流体科学部

西原 崇

我孫子研究所流体科学部

キーワード

和文英文
乱流 Turbulent flow
高レイノルズ数 High Reynolds number
大渦解析 Large-eddy simulation
スマゴリンスキー・モデル Smagorinsky model
有限要素法 Finite element method
Copyright (C) Central Research Institute of Electric Power Industry