Skip to content

Commit 06fc36e

Browse files
committed
source
1 parent 5f4a972 commit 06fc36e

File tree

443 files changed

+64109
-0
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

443 files changed

+64109
-0
lines changed

examples/advection/advection.bcd

Lines changed: 12 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,12 @@
1+
$TITL OpenFVM
2+
$FILE Boundary conditions file
3+
$BOUND 1 3 Description of $CODE
4+
10180 1 Adiabatic wall - Physical surfaces (ID u v w p T s)
5+
32 0.0 0.0 0.0 0.0 0.0 0.0
6+
10020 2 Cyclic - Physical surfaces (ID u v w p T s)
7+
30 0.0 0.0 0.0 0.0 0.0 0.0
8+
31 0.0 0.0 0.0 0.0 0.0 0.0
9+
10250 1 Volume - Physical volumes (ID u v w p T s)
10+
35 0.0 0.0 1.0 0.0 exp(-0.5*((z-0.05)/0.08)*(z-0.05)/0.08)) 0.0
11+
$ENDF
12+

examples/advection/advection.geo

Lines changed: 43 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,43 @@
1+
Mesh.MshFileVersion=1;
2+
3+
lc = 0.5;
4+
5+
dx = 0.05;
6+
dy = 0.05;
7+
dz = 1.0;
8+
9+
nz1 = 100;
10+
11+
Point(1) = {0.0,0.0,0.0,lc};
12+
Point(2) = {dx,0.0,0.0,lc};
13+
Point(3) = {dx,dy,0.0,lc};
14+
Point(4) = {0.0,dy,0.0,lc};
15+
16+
Line(1) = {1,2};
17+
Line(2) = {2,3};
18+
Line(3) = {3,4};
19+
Line(4) = {4,1};
20+
21+
Line Loop(5) = {1,2,3,4};
22+
23+
Plane Surface(6) = {5};
24+
25+
Extrude Surface {6, {0,0,dz} } {
26+
Layers { {nz1}, {9000}, {1.0} };
27+
Recombine;
28+
};
29+
30+
Recombine Surface {6};
31+
32+
// Hot surface
33+
Physical Surface(30) = {6};
34+
35+
// Cold surface
36+
Physical Surface(31) = {28};
37+
38+
// Side walls
39+
Physical Surface(32) = {15,19,23,27};
40+
41+
// Box volume
42+
Physical Volume (35) = {9000};
43+

examples/advection/advection.mtl

Lines changed: 29 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,29 @@
1+
$TITL OpenFVM
2+
$FILE Material file
3+
$MTRL 1 12 Description of $CODE
4+
20012 1 Compressibility of fluid 0
5+
0.0
6+
20015 1 Density of fluid 0
7+
1.0
8+
20021 1 Viscosity of fuild 0
9+
1E-5
10+
20022 1 Specific heat of fuild 0
11+
1.0
12+
20023 1 Thermal conductivity of fluid 0
13+
0.02
14+
21012 1 Compressibilty of fluid 1
15+
0.0
16+
21015 1 Density of fluid 1
17+
1.0
18+
21021 1 Viscosity of fluid 1
19+
1E-5
20+
21022 1 Specific heat of fuild 1
21+
1.0
22+
21023 1 Thermal conductivity of fluid 1
23+
0.18
24+
26050 1 Constant surface tension (fluid 0 / fluid 1)
25+
0.0
26+
26060 1 Thermal conductivity (boundary)
27+
10.0
28+
$ENDF
29+

examples/advection/advection.par

Lines changed: 56 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,56 @@
1+
$TITL OpenFVM
2+
$FILE Parameter file
3+
$PRMT 1 25 Description of $CODE
4+
30005 1 Convection interpolation scheme
5+
1 1 1 1 1 1
6+
30020 1 Binary output
7+
1
8+
30040 1 Calculate variable (u v w p T s)
9+
0 0 0 0 1 0
10+
30100 1 Steady state
11+
0
12+
30105 1 Convergence for steady state solutions
13+
1E-6
14+
30200 1 Adjust time interval
15+
1
16+
30201 1 Maxmimum Courant number
17+
0.9
18+
30400 1 Number of saves
19+
1
20+
30450 1 Write face scalars (u v w p T s)
21+
0 0 0 0 0 0
22+
30455 1 Write face vectors (uvw)
23+
0
24+
30460 1 Write element scalars (u v w p T s)
25+
0 0 0 0 1 0
26+
30465 1 Write element vectors (uvw)
27+
0
28+
30470 1 Write vorticity (x y z)
29+
0 0 0
30+
30485 1 Probe (u v w p T s)
31+
0 0 0 0 1 0
32+
30550 1 Maximum number of non-othorgonal corrections
33+
0
34+
30600 1 Convergence criterion (matrix solution)
35+
1E-8 1E-8 1E-8 1E-8 1E-8 1E-8
36+
30601 1 Maximum number of iterations (matrix solution)
37+
5000 5000 5000 5000 5000 5000
38+
30650 1 Matrix solver (u v w p T s) (0-Jacobi, 1-SOR, 2-CGN, 3-GMRES, 4-BiCG, 5-QMR, 6-CGS, 7-BiCGStab, 8-BiCGStabM)
39+
3 3 3 3 3 3
40+
30651 1 Matrix preconditioner (0-Null, 1-Jacobi, 2-SOR, 3-ILU)
41+
3 3 3 3 3 3
42+
30800 1 Interface scheme factor - CICSAM
43+
1.0
44+
30900 1 Maximum number of CICSAM corrections
45+
2
46+
32000 1 Start time
47+
0.0
48+
32001 1 End time
49+
1.0
50+
32002 1 Time interval
51+
0.001
52+
34000 1 Gravity vector
53+
0.0 0.0 0.0
54+
$ENDF
55+
56+

examples/backstep3D/backstep3D.bcd

Lines changed: 16 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,16 @@
1+
$TITL OpenFVM
2+
$FILE Boundary conditions file
3+
$BOUND 1 5 Description of $CODE
4+
10100 1 Inlet - Physical surfaces (ID u v w p T s)
5+
84 0.125 0.0 0.0 0.0 0.0 0.0
6+
10150 1 Outlet - Physical surfaces (ID u v w p T s)
7+
85 0.0 0.0 0.0 0.0 0.0 0.0
8+
10160 1 Slip - Physical surfaces (ID u v w p T s)
9+
82 0.0 0.0 0.0 0.0 0.0 0.0
10+
10170 2 Wall - Physical surfaces (ID u v w p T s)
11+
81 0.0 0.0 0.0 0.0 0.0 0.0
12+
83 0.0 0.0 0.0 0.0 0.0 0.0
13+
10250 1 Volume - Physical volumes (ID u v w p T s)
14+
86 0.0 0.0 0.0 0.0 0.0 0.0
15+
$ENDF
16+

examples/backstep3D/backstep3D.geo

Lines changed: 77 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,77 @@
1+
Mesh.MshFileVersion=1;
2+
3+
h = 0.04;
4+
5+
y0 = 5*h;
6+
7+
z = 7.5*h;
8+
9+
l1 = 10*h;
10+
l2 = 40*h;
11+
12+
n = 10;
13+
r = 5;
14+
15+
Point(1) = {-l1,h,0,0.1};
16+
Point(2) = {0,h,0,0.1};
17+
Point(3) = {0,0,0,0.1};
18+
Point(4) = {l2,0,0,0.1};
19+
Point(5) = {l2,h,0,0.1};
20+
Point(6) = {l2,y0+h,0,0.1};
21+
Point(7) = {0,y0+h,0,0.1};
22+
Point(8) = {-l1,y0+h,0,0.1};
23+
24+
Line(1) = {8,7};
25+
Line(2) = {6,7};
26+
Line(3) = {6,5};
27+
Line(4) = {5,2};
28+
Line(5) = {7,2};
29+
Line(6) = {8,1};
30+
Line(7) = {1,2};
31+
Line(8) = {3,2};
32+
Line(9) = {4,3};
33+
Line(10) = {4,5};
34+
35+
Line Loop(11) = {1,5,-7,-6};
36+
Plane Surface(12) = {11};
37+
Line Loop(13) = {2,5,-4,-3};
38+
Plane Surface(14) = {13};
39+
Line Loop(15) = {4,-8,-9,10};
40+
Plane Surface(16) = {15};
41+
42+
Transfinite Line {1,7} = l1/h*n+1 Using Progression 0.96;
43+
Transfinite Line {2,4,9} = l2/h*n+1 Using Progression 0.99;
44+
Transfinite Line {6,5,3} = y0/h*n+1 Using Progression 0.92;
45+
Transfinite Line {8,10} = r*n+1 Using Progression 1.0;
46+
47+
Transfinite Surface {14} = {7,6,5,2};
48+
Transfinite Surface {16} = {2,5,4,3};
49+
Transfinite Surface {12} = {7,2,1,8};
50+
Recombine Surface {14,16,12};
51+
52+
Extrude Surface {12, {0,0,z} } {
53+
Layers { {z/h*n}, {9000}, {1.0} };
54+
Recombine;
55+
};
56+
57+
Extrude Surface {14, {0,0,z} } {
58+
Layers { {z/h*n}, {9000}, {1.0} };
59+
Recombine;
60+
};
61+
62+
Extrude Surface {16, {0,0,z} } {
63+
Layers { {z/h*n}, {9000}, {1.0} };
64+
Recombine;
65+
};
66+
67+
// walls - slip
68+
Physical Surface(82) = {25,47};
69+
// walls - no slip
70+
Physical Surface(81) = {82,60,16,14,38,12};
71+
Physical Surface(83) = {77,73,33};
72+
// inlet
73+
Physical Surface(84) = {37};
74+
// outlet
75+
Physical Surface(85) = {59,81};
76+
77+
Physical Volume(86) = {9000};

examples/backstep3D/backstep3D.mtl

Lines changed: 29 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,29 @@
1+
$TITL OpenFVM
2+
$FILE Material file
3+
$MTRL 1 12 Description of $CODE
4+
20012 1 Compressibility of fluid 0
5+
0.0
6+
20015 1 Density of fluid 0
7+
1000.0
8+
20021 1 Viscosity of fuild 0
9+
1.0E-3
10+
20022 1 Specific heat of fuild 0
11+
1.0
12+
20023 1 Thermal conductivity of fluid 0
13+
1.0
14+
21012 1 Compressibilty of fluid 1
15+
0.0
16+
21015 1 Density of fluid 1
17+
1000.0
18+
21021 1 Viscosity of fluid 1
19+
1.0E-3
20+
21022 1 Specific heat of fuild 1
21+
1.0
22+
21023 1 Thermal conductivity of fluid 1
23+
1.0
24+
26050 1 Constant surface tension (fluid 0 / fluid 1)
25+
0.0
26+
26060 1 Thermal conductivity (boundary)
27+
1.0
28+
$ENDF
29+

examples/backstep3D/backstep3D.par

Lines changed: 58 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,58 @@
1+
$TITL OpenFVM
2+
$FILE Parameter file
3+
$PRMT 1 26 Description of $CODE
4+
30005 1 Convection interpolation scheme
5+
1
6+
30020 1 Binary output
7+
1
8+
30040 1 Calculate variable (u v w p T s)
9+
1 1 0 1 0 0
10+
30100 1 Steady state
11+
0
12+
30105 1 Convergence for steady state solutions
13+
1E-6
14+
30200 1 Adjust time interval
15+
1
16+
30201 1 Maxmimum Courant number
17+
0.8
18+
30400 1 Number of saves
19+
1
20+
30450 1 Write face scalars (u v w p T s)
21+
0 0 0 0 0 0
22+
30455 1 Write face vectors (uvw)
23+
0
24+
30460 1 Write element scalars (u v w p T s)
25+
0 0 0 0 0 0
26+
30465 1 Write element vectors (uvw)
27+
1
28+
30470 1 Write vorticity (x y z)
29+
0 0 0
30+
30475 1 Write stream function (xy)
31+
1
32+
30485 1 Probe (u v w p T s)
33+
0 0 0 0 0 0
34+
30550 1 Maximum number of non-othorgonal corrections
35+
0
36+
30600 1 Convergence criterion (matrix solution)
37+
1E-8
38+
30601 1 Maximum number of iterations (matrix solution)
39+
5000
40+
30650 1 Matrix solver (u v w p T s) (0-Jacobi, 1-SOR, 2-CGN, 3-GMRES, 4-BiCG, 5-QMR, 6-CGS, 7-BiCGStab, 8-BiCGStabM)
41+
7 7 7 7 7 3
42+
30651 1 Matrix preconditioner (0-Null, 1-Jacobi, 2-SOR, 3-ILU, 4-ASM)
43+
4 4 4 4 4 4
44+
30800 1 Interface scheme factor - CICSAM
45+
1.0
46+
30900 1 Maximum number of CICSAM corrections
47+
2
48+
32000 1 Start time
49+
0.0
50+
32001 1 End time
51+
1.0
52+
32002 1 Time interval
53+
0.01
54+
34000 1 Gravity vector
55+
0.0 0.0 0.0
56+
$ENDF
57+
58+

examples/channel/channel.bcd

Lines changed: 11 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,11 @@
1+
$TITL OpenFVM
2+
$FILE Boundary conditions file
3+
$BOUND 1 3 Description of $CODE
4+
10100 1 Inlet - Physical surfaces (ID u v w p T s)
5+
73 0.0 0.0 -0.1 0.0 0.0 0.0
6+
10150 1 Outlet - Physical surfaces (ID u v w p T s)
7+
74 0.0 0.0 0.0 0.0 0.0 0.0
8+
10170 1 Wall - Physical surfaces (ID u v w p T s)
9+
75 0.0 0.0 0.0 0.0 0.0 0.0
10+
$ENDF
11+

examples/channel/channel.geo

Lines changed: 49 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,49 @@
1+
Mesh.MshFileVersion=1;
2+
nT=25;
3+
nC=15;
4+
nt=8;
5+
T=1;
6+
lc=0.1;
7+
lc2=0.1;
8+
9+
Point(1) = {0.0,0.0,0.0,lc};
10+
Point(2) = {0.0,1,0.0,lc};
11+
Point(3) = {1,0,0.0,lc2};
12+
Point(4) = {-1,0,0.0,lc};
13+
Point(5) = {0,-1,0.0,lc};
14+
15+
Line(1) = {2,3};
16+
Line(2) = {3,5};
17+
Line(3) = {5,4};
18+
Line(4) = {4,2};
19+
Line Loop(5) = {1,2,3,4};
20+
Plane Surface(6) = {5};
21+
22+
Transfinite Line {2,3,4,1} = 8 Using Progression 1.0;
23+
Transfinite Surface {6} = {2,3,5,4};
24+
Recombine Surface {6};
25+
26+
Extrude Surface {6, {0.0,0.0,5.0}}
27+
{ Recombine; Layers { {nT}, {9000}, {T} }; };
28+
29+
Extrude Surface {6, {1.0,1.0,0.0}, {2.0,0.0,0.0}, -Pi/2}
30+
{ Recombine; Layers { {nC}, {9001}, {T} }; };
31+
32+
Extrude Surface {50, {1.0,-1.0,0.0}}
33+
{ Recombine; Layers { {nt}, {9002}, {T} }; };
34+
35+
36+
// Inlet
37+
Physical Surface(73) = {28};
38+
39+
// Outlet
40+
Physical Surface(74) = {72};
41+
//Physical Surface(74) = {6};
42+
//Physical Surface(74) = {50};
43+
44+
// Walls
45+
Physical Surface(75) = {15,27,23,19,37,45,49,41,59,67,71,63};
46+
47+
Physical Volume(76) = {9000};
48+
Physical Volume(77) = {9001};
49+
Physical Volume(78) = {9002};

0 commit comments

Comments
 (0)