OpenFOAMで、通常のtutorialsケースに対して、potentialFoam を実行すると、
Calculating potential flow
keyword laplacian(1,p) is undefined in dictionary "/home/et/OpenFOAM/et-1.5- dev/run/pitzDaily/system/fvSchemes::laplacianSchemes"
file: /home/et/OpenFOAM/et-1.5-dev/run/pitzDaily/system/fvSchemes::laplacianSchemes from line 43 to line 49.
From function dictionary::lookupEntry(const word& keyword) const
in file db/dictionary/dictionary.C at line 213.
FOAM exiting
となって計算できない。
system/fvSchemes にて、
laplacian(1,p)
などと、定義してやれば実行してくれるのだが、一々書き換えるのも面倒でした。
そういう時に、pyFoamPotentialRunner.pyが使えます。
et@ubuntu-vm:~/OpenFOAM/et-1.5-dev/run/pitzDaily$ pyFoamPotentialRunner.py .
Setting system-directory for potentialFoam
/*—————————————————————————*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5-dev |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*—————————————————————————*/
Exec : potentialFoam -case . -writep
Date : Feb 07 2009
Time : 11:37:20
Host : ubuntu-vm
PID : 9914
Case : ./.
nProcs : 1
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p
Reading field U
Calculating potential flow
PCG: Solving for p, Initial residual = 1, Final residual = 0.00943312, No Iterations 172
continuity error = 0.16506
Interpolated U error = 0.000116284
mag(U): max: 27.0928 min: 0.127263
ExecutionTime = 0.51 s ClockTime = 1 s
End
その結果、
et@ubuntu-vm:~/OpenFOAM/et-1.5-dev/run/pitzDaily$ ls 0 -l
\u5408\u8a08 820
-rwxr-xr-x 1 et et 1324 2009-02-07 11:32 R
-rwxr-xr-x 1 et et 397279 2009-02-07 11:37 U
-rw-r–r– 1 et et 1387 2009-02-07 11:37 U.prepotential
-rwxr-xr-x 1 et et 1312 2009-02-07 11:32 epsilon
-rwxr-xr-x 1 et et 1304 2009-02-07 11:32 k
-rwxr-xr-x 1 et et 1302 2009-02-07 11:32 nuTilda
-rwxr-xr-x 1 et et 102503 2009-02-07 11:37 p
-rw-r–r– 1 et et 1296 2009-02-07 11:37 p.prepotential
-rw-r–r– 1 et et 298674 2009-02-07 11:37 phi
U、pにポテンシャル流れ場の結果が入って、
元々の初期条件は、U.prepotential、p.prepotentialに保存される。。。ということです。
ああ、便利。。。⇒これもDEXCS2009で採用しよう！