interFoamでの粘度式の改造

[aki]

OpenFoam初心者のため、知識不足が多々あるかと思います。

またこのグループで初めて質問を投稿するため、失礼があるかもしれません。
大学でinterFoamで計算をさせて研究をしております。
元の粘度の式v=α1​v1​+(1−α1​)v2​

を界面付近で粘度が大きく上昇していくような関数に改造したいのですがやり方がよくわからないのでご教授お願いします。OpenFoamのバージョンはv2206です。

変更したい式↓

ν=αν1​+(1−α)ν2​+νint×​exp(−(α−0.5)^2/2w^2​)

vintとwは自分で決める定数です

[aki]

計算というファイルがinterFoamで計算させたいファイルです。
もう一つのがcatGPTで調べながらやってみたのですがうまくいきませんでした。
myincompressibleの myincompressibletwoPhaseMixtureの中を変えてコンパイルさせたのですができませんでした。

2026年1月23日金曜日 0:43:16 UTC+9 aki: