k, omega の発散は抗力係数にどの程度影響するのかについて
166 views
Skip to first unread message
Daichry
Oct 25, 2020, 9:11:57 AM10/25/20
to OpenFOAM
初めまして、大学の研究で橋梁の桁についての数値解析を行っている谷本と申します。
OpenFOAMを使用し始めて3カ月程度の初心者です。
現在、橋梁のモデルケースを解析する前段階で長方形断面のモデルを使用していますが、境界条件やスキームなどの設定が上手くいかず抗力係数が過去の実験結果（辺長比が同じ長方形断面だと2.1程度）とかけ離れたものとなってしまいます。
原因として考えられるのはｋとomegaの発散ですが、それらをピンポイントで抑える方法が分からないためご教授いただけましたら幸いです。
稚拙な質問で申し訳ありませんが、抗力係数がおかしい原因が分からず困っております。是非よろしくお願いします。
以下モデルケース、チェックメッシュの結果、初期条件、スキーム、ソルバー、logです。
#Mesh作成にはANSYSのstudent.verを使っておりポリヘドラルメッシュ（50万）をfluent3DMeshToFoamで変換して利用しております。
#境界メッシュは10層で奥行きが26㎜、縦78㎜、横292㎜の長方形断面です。
#ライセンスの関係で奥行きが十分に取れず薄くなっています。
添付画像のMeshdata1~3に詳細を添付しました。
#非圧縮でpimpleFoam, RANSのkomegaSSTモデルを使用して計算を行っております。
#計算修了後、yPlusの計算を行った結果はこちらです。
*******************************************************
yPlus
Create time
Create mesh for time = 30
PIMPLE: No convergence criteria found
PIMPLE: No corrector convergence criteria found
Calculations will do 5 corrections
PIMPLE: Operating solver in transient mode with 5 outer correctors
Time = 30
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
model kOmegaSST;
turbulence on;
printCoeffs on;
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}
No MRF models present
No finite volume options present
yPlus yPlus write:
writing object yPlus
patch deck y+ : min = 0.411235, max = 149.339, average = 52.2314
***********************************************************************
checkMesh
Create time
Create polyMesh for time = 0
Time = 0
Mesh stats
points: 2285408
faces: 3059556
internal faces: 2954573
cells: 506149
faces per cell: 11.8821
boundary patches: 6
point zones: 0
face zones: 1
cell zones: 1
Overall number of cells of each type:
hexahedra: 336
prisms: 0
wedges: 1
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 505812
Breakdown of polyhedra by number of faces:
faces number of cells
6 6
7 25907
8 95146
9 13715
10 21574
11 46592
12 54474
13 71842
14 77388
15 57095
16 28689
17 10218
18 2598
19 490
20 67
21 10
22 1
Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
inlet 175 300 ok (non-closed singly connected)
upperwall 375 417 ok (non-closed singly connected)
lowerwall 373 430 ok (non-closed singly connected)
outlet 160 303 ok (non-closed singly connected)
frontandback 91001 167611 ok (non-closed singly connected)
deck 12899 25697 ok (non-closed singly connected)
Checking geometry...
Overall domain bounding box (-1170 -1014 0) (2242 936 26)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (-2.77701e-19 4.55947e-19 1.03612e-14) OK.
Max cell openness = 3.32536e-16 OK.
Max aspect ratio = 57.8045 OK.
Minimum face area = 0.000180675. Maximum face area = 10457.6. Face area magnitudes OK.
Min volume = 0.00263102. Max volume = 135819. Total volume = 1.72396e+08. Cell volumes OK.
Mesh non-orthogonality Max: 66.0569 average: 4.94713
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 2.50466 OK.
Coupled point location match (average 0) OK.
Mesh OK.
*******************************************************************
log の一部
Build : 8-9b73cf21a682
Exec : pimpleFoam
Date : Oct 23 2020
Time : 14:14:44
Host : "fujical-VirtualBox"
PID : 20196
I/O : uncollated
Case : /home/fuji-cal/OpenFOAM/fuji-cal-8/run/2ddeckmodel.102301
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
PIMPLE: No convergence criteria found
PIMPLE: No corrector convergence criteria found
Calculations will do 5 corrections
PIMPLE: Operating solver in transient mode with 5 outer correctors
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
RAS
{
model kOmegaSST;
turbulence on;
printCoeffs on;
alphaK1 0.85;
alphaK2 1;
alphaOmega1 0.5;
alphaOmega2 0.856;
gamma1 0.555556;
gamma2 0.44;
beta1 0.075;
beta2 0.0828;
betaStar 0.09;
a1 0.31;
b1 1;
c1 10;
F3 false;
}
No MRF models present
No finite volume options present
Courant Number mean: 6.19131e-05 max: 0.0785238
--> FOAM Warning :
From function void Foam::timeControl::read(const Foam::dictionary&)
in file db/functionObjects/timeControl/timeControl.C at line 89
Reading "/home/fuji-cal/OpenFOAM/fuji-cal-8/run/2ddeckmodel.102301/system/controlDict/functions/forceCoeffs1" from line 56 to line 73
Using deprecated 'outputControl'
Please use 'writeControl' with 'writeInterval'
forces forceCoeffs1:
Not including porosity effects
forceCoeffs forceCoeffs1:
Not including porosity effects
Starting time loop
forceCoeffs forceCoeffs1 write:
Cm = 7.84381e+06
Cd = 15698.7
Cl = -6.07253e-17
Cl(f) = 7.84381e+06
Cl(r) = -7.84381e+06
Courant Number mean: 6.19131e-05 max: 0.0785238
deltaT = 0.0012
Time = 0.0012
PIMPLE: Iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.61679e-05, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 1.81678e-05, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 5.50009e-06, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.0098103, No Iterations 401
DICPCG: Solving for p, Initial residual = 0.379531, Final residual = 0.00378355, No Iterations 147
time step continuity errors : sum local = 1.72564e-09, global = -1.5697e-10, cumulative = -1.5697e-10
DICPCG: Solving for p, Initial residual = 0.0512197, Final residual = 0.000466348, No Iterations 278
DICPCG: Solving for p, Initial residual = 0.0195146, Final residual = 9.53883e-07, No Iterations 433
time step continuity errors : sum local = 7.6857e-13, global = -7.29039e-15, cumulative = -1.56977e-10
PIMPLE: Iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 0.0691074, Final residual = 9.6541e-09, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.211492, Final residual = 1.03328e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.214156, Final residual = 4.31805e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.0256359, Final residual = 0.000252309, No Iterations 390
DICPCG: Solving for p, Initial residual = 0.097952, Final residual = 0.00097534, No Iterations 142
time step continuity errors : sum local = 6.96188e-10, global = -5.69477e-11, cumulative = -2.13925e-10
DICPCG: Solving for p, Initial residual = 0.013071, Final residual = 0.000128193, No Iterations 260
DICPCG: Solving for p, Initial residual = 0.00827721, Final residual = 9.07004e-07, No Iterations 400
time step continuity errors : sum local = 7.57534e-13, global = 1.23749e-15, cumulative = -2.13923e-10
PIMPLE: Iteration 3
DILUPBiCG: Solving for Ux, Initial residual = 2.97769e-07, Final residual = 2.97769e-07, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 3.6411e-07, Final residual = 3.6411e-07, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 6.36898e-07, Final residual = 6.36898e-07, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00353934, Final residual = 3.51338e-05, No Iterations 107
DICPCG: Solving for p, Initial residual = 0.00263034, Final residual = 2.47212e-05, No Iterations 30
time step continuity errors : sum local = 2.09063e-11, global = -3.34571e-12, cumulative = -2.17269e-10
DICPCG: Solving for p, Initial residual = 0.00199849, Final residual = 1.93237e-05, No Iterations 30
DICPCG: Solving for p, Initial residual = 0.00174241, Final residual = 9.81544e-07, No Iterations 366
time step continuity errors : sum local = 8.30641e-13, global = 3.50788e-15, cumulative = -2.17266e-10
PIMPLE: Iteration 4
DILUPBiCG: Solving for Ux, Initial residual = 1.67031e-08, Final residual = 1.67031e-08, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 2.50486e-08, Final residual = 2.50486e-08, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 7.05902e-08, Final residual = 7.05902e-08, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00162569, Final residual = 1.56526e-05, No Iterations 28
DICPCG: Solving for p, Initial residual = 0.00146115, Final residual = 1.36897e-05, No Iterations 24
time step continuity errors : sum local = 1.16128e-11, global = -8.55153e-13, cumulative = -2.18121e-10
DICPCG: Solving for p, Initial residual = 0.00138129, Final residual = 1.38059e-05, No Iterations 22
DICPCG: Solving for p, Initial residual = 0.0013202, Final residual = 9.68584e-07, No Iterations 256
time step continuity errors : sum local = 8.21722e-13, global = -1.67549e-14, cumulative = -2.18138e-10
PIMPLE: Iteration 5
DILUPBiCG: Solving for Ux, Initial residual = 5.58491e-09, Final residual = 5.58491e-09, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 8.7926e-09, Final residual = 8.7926e-09, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 1.3724e-08, Final residual = 1.3724e-08, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00129284, Final residual = 1.25477e-05, No Iterations 28
DICPCG: Solving for p, Initial residual = 0.00124621, Final residual = 1.17323e-05, No Iterations 23
time step continuity errors : sum local = 9.95841e-12, global = -2.15339e-13, cumulative = -2.18353e-10
DICPCG: Solving for p, Initial residual = 0.00121898, Final residual = 1.18152e-05, No Iterations 21
DICPCG: Solving for p, Initial residual = 0.00119244, Final residual = 9.5548e-07, No Iterations 89
time step continuity errors : sum local = 8.1105e-13, global = -1.83956e-13, cumulative = -2.18537e-10
DILUPBiCG: Solving for omega, Initial residual = 0.000366817, Final residual = 7.56339e-10, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 6.55929e-08, No Iterations 1
ExecutionTime = 261.91 s ClockTime = 262 s
forceCoeffs forceCoeffs1 write:
Cm = 1.50761e+14
Cd = 3.01864e+11
Cl = -1.72679e+07
Cl(f) = 1.50761e+14
Cl(r) = -1.50761e+14
Courant Number mean: 7.45941e-05 max: 0.0651621
deltaT = 0.00143654
Time = 0.00263654
PIMPLE: Iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 0.264943, Final residual = 4.6607e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.445641, Final residual = 2.99096e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.447424, Final residual = 1.13014e-08, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.0742494, Final residual = 0.00074037, No Iterations 387
DICPCG: Solving for p, Initial residual = 0.666121, Final residual = 0.0063918, No Iterations 142
time step continuity errors : sum local = 3.59176e-09, global = 3.0614e-10, cumulative = 8.76032e-11
DICPCG: Solving for p, Initial residual = 0.135759, Final residual = 0.00131444, No Iterations 278
DICPCG: Solving for p, Initial residual = 0.0573419, Final residual = 9.63783e-07, No Iterations 461
time step continuity errors : sum local = 3.95298e-13, global = -3.63801e-16, cumulative = 8.76028e-11
PIMPLE: Iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 1.7419e-06, Final residual = 1.7419e-06, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 1.91935e-06, Final residual = 1.91935e-06, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 2.77208e-06, Final residual = 2.77208e-06, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.018042, Final residual = 0.00017486, No Iterations 281
DICPCG: Solving for p, Initial residual = 0.011211, Final residual = 0.000109483, No Iterations 34
time step continuity errors : sum local = 5.1197e-11, global = 5.92234e-12, cumulative = 9.35252e-11
DICPCG: Solving for p, Initial residual = 0.00733238, Final residual = 7.26042e-05, No Iterations 89
DICPCG: Solving for p, Initial residual = 0.00596898, Final residual = 9.67796e-07, No Iterations 352
time step continuity errors : sum local = 4.54828e-13, global = -5.05972e-15, cumulative = 9.35201e-11
PIMPLE: Iteration 3
DILUPBiCG: Solving for Ux, Initial residual = 2.64362e-08, Final residual = 2.64362e-08, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 4.84712e-08, Final residual = 4.84712e-08, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 2.65505e-07, Final residual = 2.65505e-07, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00518144, Final residual = 4.94749e-05, No Iterations 8
DICPCG: Solving for p, Initial residual = 0.00474364, Final residual = 4.20451e-05, No Iterations 6
time step continuity errors : sum local = 1.98797e-11, global = 1.07356e-12, cumulative = 9.45937e-11
DICPCG: Solving for p, Initial residual = 0.00448853, Final residual = 3.93608e-05, No Iterations 5
DICPCG: Solving for p, Initial residual = 0.00436799, Final residual = 9.92001e-07, No Iterations 336
time step continuity errors : sum local = 4.69085e-13, global = 5.03017e-15, cumulative = 9.45987e-11
PIMPLE: Iteration 4
DILUPBiCG: Solving for Ux, Initial residual = 4.52698e-09, Final residual = 4.52698e-09, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 5.14559e-09, Final residual = 5.14559e-09, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 1.95722e-08, Final residual = 1.95722e-08, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00431949, Final residual = 4.03306e-05, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.00422228, Final residual = 2.72179e-05, No Iterations 3
time step continuity errors : sum local = 1.28875e-11, global = 1.07822e-13, cumulative = 9.47065e-11
DICPCG: Solving for p, Initial residual = 0.00415975, Final residual = 2.71514e-05, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.00412515, Final residual = 9.19757e-07, No Iterations 245
time step continuity errors : sum local = 4.35502e-13, global = 6.41049e-15, cumulative = 9.47129e-11
PIMPLE: Iteration 5
DILUPBiCG: Solving for Ux, Initial residual = 1.52748e-09, Final residual = 1.52748e-09, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 1.98898e-09, Final residual = 1.98898e-09, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 5.13226e-09, Final residual = 5.13226e-09, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00410072, Final residual = 2.65004e-05, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.00407441, Final residual = 3.87458e-05, No Iterations 2
time step continuity errors : sum local = 1.83515e-11, global = 2.20863e-14, cumulative = 9.4735e-11
DICPCG: Solving for p, Initial residual = 0.00405929, Final residual = 2.54138e-05, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.00403699, Final residual = 9.95953e-07, No Iterations 46
time step continuity errors : sum local = 4.71722e-13, global = 8.27465e-14, cumulative = 9.48178e-11
DILUPBiCG: Solving for omega, Initial residual = 0.000437939, Final residual = 1.82893e-09, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.306679, Final residual = 1.33175e-08, No Iterations 1
ExecutionTime = 445.64 s ClockTime = 446 s
*************************************************************
LastTimeに飛びます
Courant Number mean: 0.00372625 max: 1.42495
deltaT = 0.06
Time = 30
PIMPLE: Iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 1.34724e-05, Final residual = 1.12682e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.000173067, Final residual = 3.7028e-07, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.00093385, Final residual = 5.20196e-07, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.00827369, Final residual = 8.27173e-05, No Iterations 15
DICPCG: Solving for p, Initial residual = 0.00407076, Final residual = 3.88498e-05, No Iterations 9
time step continuity errors : sum local = 6.84859e-12, global = -2.62363e-14, cumulative = 7.11667e-11
DICPCG: Solving for p, Initial residual = 0.00362324, Final residual = 3.58469e-05, No Iterations 8
DICPCG: Solving for p, Initial residual = 0.00301529, Final residual = 9.93042e-07, No Iterations 125
time step continuity errors : sum local = 1.74653e-13, global = -2.0413e-14, cumulative = 7.11463e-11
PIMPLE: Iteration 2
DILUPBiCG: Solving for Ux, Initial residual = 4.84091e-08, Final residual = 4.84091e-08, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 3.36764e-07, Final residual = 3.36764e-07, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 2.71462e-06, Final residual = 2.71462e-06, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00385197, Final residual = 3.78654e-05, No Iterations 22
DICPCG: Solving for p, Initial residual = 0.00234348, Final residual = 2.22723e-05, No Iterations 7
time step continuity errors : sum local = 3.91108e-12, global = -2.07317e-14, cumulative = 7.11256e-11
DICPCG: Solving for p, Initial residual = 0.00194461, Final residual = 1.79149e-05, No Iterations 9
DICPCG: Solving for p, Initial residual = 0.00171581, Final residual = 9.98594e-07, No Iterations 84
time step continuity errors : sum local = 1.7531e-13, global = -1.83146e-14, cumulative = 7.11072e-11
PIMPLE: Iteration 3
DILUPBiCG: Solving for Ux, Initial residual = 7.88687e-10, Final residual = 7.88687e-10, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 6.44574e-09, Final residual = 6.44574e-09, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 5.11874e-08, Final residual = 5.11874e-08, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00163541, Final residual = 1.60644e-05, No Iterations 5
DICPCG: Solving for p, Initial residual = 0.00144347, Final residual = 1.27454e-05, No Iterations 5
time step continuity errors : sum local = 2.23842e-12, global = -1.84831e-14, cumulative = 7.10888e-11
DICPCG: Solving for p, Initial residual = 0.00139004, Final residual = 1.34224e-05, No Iterations 4
DICPCG: Solving for p, Initial residual = 0.00128373, Final residual = 9.4911e-07, No Iterations 43
time step continuity errors : sum local = 1.66765e-13, global = -1.77811e-14, cumulative = 7.1071e-11
PIMPLE: Iteration 4
DILUPBiCG: Solving for Ux, Initial residual = 1.77377e-10, Final residual = 1.77377e-10, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 1.12403e-09, Final residual = 1.12403e-09, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 9.6197e-09, Final residual = 9.6197e-09, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00124774, Final residual = 1.23973e-05, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.00116837, Final residual = 1.0015e-05, No Iterations 3
time step continuity errors : sum local = 1.76037e-12, global = -1.78131e-14, cumulative = 7.10532e-11
DICPCG: Solving for p, Initial residual = 0.0011364, Final residual = 9.25199e-06, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.00108061, Final residual = 9.93739e-07, No Iterations 21
time step continuity errors : sum local = 1.74732e-13, global = -1.77744e-14, cumulative = 7.10354e-11
PIMPLE: Iteration 5
DILUPBiCG: Solving for Ux, Initial residual = 1.05469e-10, Final residual = 1.05469e-10, No Iterations 0
DILUPBiCG: Solving for Uy, Initial residual = 6.4617e-10, Final residual = 6.4617e-10, No Iterations 0
DILUPBiCG: Solving for Uz, Initial residual = 6.46747e-09, Final residual = 6.46747e-09, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.00106605, Final residual = 8.96248e-06, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.00101868, Final residual = 8.87481e-06, No Iterations 2
time step continuity errors : sum local = 1.56081e-12, global = -1.77728e-14, cumulative = 7.10176e-11
DICPCG: Solving for p, Initial residual = 0.00100858, Final residual = 8.98269e-06, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.000981306, Final residual = 9.92853e-07, No Iterations 20
time step continuity errors : sum local = 1.74645e-13, global = -1.76776e-14, cumulative = 7.09999e-11
DILUPBiCG: Solving for omega, Initial residual = 0.000175908, Final residual = 2.18351e-07, No Iterations 1
bounding omega, min: -0.365801 max: 145.896 average: 1.23775
DILUPBiCG: Solving for k, Initial residual = 3.61166e-05, Final residual = 8.45969e-09, No Iterations 1
bounding k, min: -0.000372912 max: 0.0798642 average: 0.00103804
ExecutionTime = 31028.7 s ClockTime = 31043 s
forceCoeffs forceCoeffs1 write:
Cm = 9.95693e+08
Cd = 1.97131e+06
Cl = 4456.16
Cl(f) = 9.95696e+08
Cl(r) = -9.95691e+08
End
以上、宜しくお願いいたします
Daichry
Oct 25, 2020, 9:20:48 AM10/25/20
to OpenFOAM
ファイルの添付ができなかったためこちらに添付します。
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
ｋ
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.00015;
boundaryField
{
inlet
{
type fixedValue;
value uniform 0.00015;
}
outlet
{
type zeroGradient;
}
upperwall
{
type symmetryPlane;
}
lowerwall
{
type symmetryPlane;
}
frontandback
{
type slip;
}
deck
{
type kqRWallFunction;
value uniform 0.00015;
}
}
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
nut
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type zeroGradient;
}
upperwall
{
type symmetryPlane;
}
lowerwall
{
type symmetryPlane;
}
frontandback
{
type slip;
}
deck
{
type nutkWallFunction;
value uniform 0;
}
}
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
omega
dimensions [0 0 -1 0 0 0 0];
internalField uniform 0.764;
boundaryField
{
inlet
{
type fixedValue;
value uniform 0.764;
}
outlet
{
type zeroGradient;
}
upperwall
{
type symmetryPlane;
}
lowerwall
{
type symmetryPlane;
}
frontandback
{
type slip;
}
deck
{
type omegaWallFunction;
value uniform 0.764;
}
}
****************************************************
p
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
upperwall
{
type symmetryPlane;
}
lowerwall
{
type symmetryPlane;
}
frontandback
{
type slip;
}
deck
{
type zeroGradient;
}
}
**************************************************
U
dimensions [0 1 -1 0 0 0 0];
internalField uniform (1 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (1 0 0);
}
outlet
{
type zeroGradient;
}
upperwall
{
type symmetryPlane;
}
lowerwall
{
type symmetryPlane;
}
frontandback
{
type slip;
}
deck
{
type fixedValue;
value uniform (0 0 0);
}
}
**********************************************
boundary
6
(
inlet
{
type patch;
nFaces 175;
startFace 2954573;
}
upperwall
{
type symmetryPlane;
inGroups List<word> 1(symmetryPlane);
nFaces 375;
startFace 2954748;
}
lowerwall
{
type symmetryPlane;
inGroups List<word> 1(symmetryPlane);
nFaces 373;
startFace 2955123;
}
outlet
{
type patch;
nFaces 160;
startFace 2955496;
}
frontandback
{
type slip;
inGroups List<word> 1(slip);
nFaces 91001;
startFace 2955656;
}
deck
{
type wall;
inGroups List<word> 1(wall);
nFaces 12899;
startFace 3046657;
}
)
********************************************
controlDict
application pimpleFoam;
startFrom latestTime;
startTime 0;
stopAt endTime;
endTime 30;
deltaT 0.001;
writeControl adjustableRunTime;
writeInterval 0.3;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
adjustTimeStep true;
maxCo 5;
functions
{
forceCoeffs1
{
type forceCoeffs; //Mandatory entries
functionObjectLibs ("libforces.so");
patches (deck); //change to your patch name
outputControl timeStep; //optional entries
timeInterval 0.1;
log yes;
p p; //Field names & reference
U U;
rho rhoInf;
log true;
rhoInf 1; //Fake
CofR (0.146 -0.039 0.013); //origin of the moment calculation
liftDir (0 1 0); //direction on Cl
dragDir (1 0 0); //direction on Cd
pitchAxis (0 0 1); //pitch axis
magUInf 1; //Freestream velocity magnitude[m/s]
lRef 0.078; //Reference length i.e. h(0.078)[m]
Aref 0.002028; //Reference area B(0.026)*h(0.078)[m2]
}
}
*******************************************
fvScheme
ddtSchemes
{
default backward;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss linearUpwind grad(U);
div(phi,k) Gauss linearUpwind grad(k);
div(phi,omega) Gauss linearUpwind grad(omega);
div(phi,R) Gauss linearUpwind;
div(R) Gauss linearUpwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited 0.5;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited 0.5;
}
fluxRequired
{
default no;
p ;
}
wallDist
{
method meshWave;
}
*************************************************
fvSolution
solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-6;
relTol 0.01;
}
pFinal
{
$p;
tolerance 1e-6;
relTol 0;
}
"(U|k|omega)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
"(U|k|omega)Final"
{
$U;
tolerance 1e-05;
relTol 0;
}
}
PIMPLE
{
nNonOrthogonalCorrectors 1;
nOuterCorrectors 5;
nCorrectors 2;
}
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
長くなりましたが、宜しくお願いいたします。
Reply all
Reply to author
Forward
0 new messages