OnWork: initial approximation via potentialFoam

This commit is contained in:
L-Nafaryus 2021-02-24 15:36:17 +05:00
parent b8d46e3a93
commit a9959b1f68
12 changed files with 110819 additions and 363029 deletions

16
run.sh
View File

@ -1,14 +1,24 @@
#!/usr/bin/env bash
path=build/simple-cubic/0.1
path=.
#build/simple-cubic/0.1
# python src/genmesh.py
# python src/prefoam.py
rm -r postProcessing processor* *.log logs constant/polymesh
ideasUnvFoam -case $path mesh.unv
checkMesh -case $path
polyDualMesh 70 -overwrite
checkMesh -case $path -allGeometry -allTopology
foamDictionary -case $path constant/polyMesh/boundary -entry entry0.wall.type -set wall
#potentialFoam -case $path
potentialFoam -case $path
decomposePar -case $path
foamDictionary -case $path 0/U -entry boundary.inlet.type -set pressureInletVelocity
foamDictionary -case $path 0/U -entry boundary.inlet.value -set 'uniform (0 0 0)'
mpirun -np 4 --oversubscribe simpleFoam -parallel -case $path > $path/simpleFoam.log

View File

@ -10,7 +10,7 @@ if not os.path.exists(build):
###
alpha = [0.1, 0.15, 0.2]
alpha = [0.1] #, 0.15, 0.2]
simpleCubic = os.path.join(src, "simple-cubic/main.py")
# Body-centered cubic

File diff suppressed because it is too large Load Diff

View File

@ -10,32 +10,33 @@ FoamFile
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
internalField uniform 0;
boundaryField
{
inlet
{
type fixedValue;
value uniform 0.001;
type fixedValue;
value uniform 0.001;
}
outlet
{
type fixedValue;
value uniform 0;
type fixedValue;
value uniform 0;
}
wall
{
type zeroGradient;
type zeroGradient;
}
}
// ************************************************************************* //

View File

@ -1,40 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
3
(
inlet
{
type patch;
nFaces 168;
startFace 109912;
}
outlet
{
type patch;
nFaces 164;
startFace 110080;
}
wall
{
type patch;
nFaces 5834;
startFace 110244;
}
)
// ************************************************************************* //

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

File diff suppressed because it is too large Load Diff

View File

@ -1,249 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class faceSet;
location "constant/polyMesh/sets";
object nonOrthoFaces;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
225
(
86025
83989
86038
83995
104489
26674
81986
100435
84052
82006
108643
90219
108654
26734
108656
26764
80016
41104
88215
41119
77985
94381
77997
41134
84163
26824
41164
84187
61659
108785
82170
100604
96543
100642
84263
88360
100661
96574
92496
24924
88416
100713
104810
12664
104828
24969
100749
104850
24999
100787
100814
88533
104932
98797
12784
78333
16909
80413
80419
80435
92729
86606
92796
10879
92810
109203
10909
25269
58039
103102
8899
25284
99043
25329
88831
25344
78595
88839
8969
78631
105277
105279
25409
99138
9029
82763
80734
17264
109434
80767
101262
86927
9109
27549
88990
68519
84904
82868
68539
105417
101338
101340
107494
107496
95245
42004
97319
97321
44114
91222
44124
78953
109678
89198
44144
99485
101549
101552
44214
97473
87244
5324
105694
5344
97506
95465
109815
44284
83197
105732
5384
101641
42249
101653
25884
5404
101674
87342
25904
5424
27974
101706
107861
25944
25964
83311
81271
81272
89468
87420
97675
26004
97702
40359
105900
103856
60854
79287
81356
83406
26064
95704
79322
60894
42484
95766
93733
91686
99882
93751
93753
104010
99927
13914
99933
81501
83552
99940
91751
99944
91753
85612
104059
95868
85646
83610
83627
108209
26304
106191
77533
77536
95991
102141
100113
77585
83744
79660
104247
96058
106312
42824
16214
94043
102240
106340
16234
81790
16254
77708
61324
89997
42894
83858
100248
83867
85923
83877
85926
83878
16294
90028
83907
96214
96215
)
// ************************************************************************* //

View File

@ -1,458 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class faceSet;
location "constant/polyMesh/sets";
object skewFaces;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
434
(
24668
24672
24676
24680
24683
24688
24692
24696
57468
24700
57472
57476
57480
57483
57488
57492
57496
57500
57503
57508
51495
51498
51501
51508
51512
51516
51520
51523
12628
12632
51548
51552
51556
12648
12652
12656
12660
12663
12668
12672
12676
12680
57738
12683
57742
12688
57746
12692
12696
12700
57758
12703
12708
12712
12716
12720
12723
12728
57788
12732
57792
12736
57796
12740
12743
57800
57803
12748
12752
12756
12760
57818
12763
57822
12768
57826
12772
57830
12776
57833
12780
12783
12788
57848
12792
57852
12796
57856
12800
57860
57863
57878
57882
57886
57890
57893
8953
8957
8973
8977
8981
8985
8988
8993
8997
9001
9005
9008
9013
9017
9021
9025
9028
9033
9037
9041
9045
9048
9053
9057
9061
54163
25503
54183
54187
54191
54195
25523
54198
25527
54203
25531
54207
25535
25538
54211
25543
25547
25551
25555
25558
64508
64512
64516
64520
64523
64528
64532
64536
64540
64543
5268
9367
5272
9371
5276
9375
5280
5283
9383
5288
9387
5292
9391
5296
9395
5300
9398
5303
5308
5312
5316
9413
5320
9417
40138
5323
9421
25808
5328
9425
9428
5332
5336
5340
40158
5343
40162
9443
25828
5348
40166
9447
25832
5352
40170
9451
5356
40173
9455
5360
40178
9458
5363
40182
5368
5372
5376
9473
25858
5380
9477
25862
5383
9481
25866
5388
9485
25870
9488
5392
5396
5400
5403
9503
25888
5408
9507
25892
5412
9511
25896
5416
9515
25900
5420
9518
25903
5423
5428
5432
5436
25918
5440
25922
25926
25930
25933
25948
25952
25956
25960
25963
25978
25982
25986
25990
25993
26008
26012
26016
26020
26023
26038
26042
26046
50623
26050
50627
26053
50631
50635
50638
50643
26068
50647
26072
50651
26076
50655
26080
50658
26083
26098
26102
26106
26110
26113
26128
26132
26136
26140
26143
26158
26162
26166
26170
26173
26188
26192
26196
26200
26203
26218
26222
26226
26230
26233
26248
26252
26256
26260
26263
26278
26282
26286
26290
26293
26308
26312
26316
26320
26323
50903
50907
50911
26338
50915
50918
26342
26346
50923
40683
26350
50927
40687
26353
50931
40691
50935
40695
40698
50943
40703
26368
50947
40707
26372
40711
26376
40715
26380
40718
26383
40723
40727
40731
26398
26402
26406
40743
26410
40747
26413
40751
40755
16198
61258
61262
61266
61270
61273
16218
61278
16222
61282
16226
61286
16230
61290
61293
16238
16242
16246
16250
61308
16253
61312
16258
61316
16262
61320
16266
61323
16270
16273
16278
61338
16282
61342
16286
61346
16290
61350
61353
16298
16302
16306
16310
61368
61372
16318
61376
16322
61380
16326
61383
16330
16333
16338
61398
16342
61402
16346
61406
16350
16353
61410
61413
16358
16362
16366
16370
16373
)
// ************************************************************************* //

View File

@ -58,6 +58,10 @@ solvers
potentialFlow
{
nNonOrthogonalCorrectors 10;
PhiRefCell 0;
PhiRefPoint 0;
PhiRefValue 0;
Phi 0;
}
cache