Skip to content

Commit

Permalink
Further developments.
Browse files Browse the repository at this point in the history
  • Loading branch information
@nikola-m
nikola-m committed Sep 29, 2020
1 parent 9e90100 commit 6f4d0b9
Show file tree
Hide file tree
Showing 244 changed files with 30,649 additions and 821 deletions.
Empty file modified Makefile
100644 → 100755
View file
Empty file.
19 changes: 19 additions & 0 deletions examples/RBCells/0/T
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
internalField
uniform
300.0
boundaryField
floor
Dirichlet
uniform
301.0
ceiling
Dirichlet
uniform
300.0
sideWalls
Neumann
zeroGradient
frontAndBack
Neumann
zeroGradient

21 changes: 21 additions & 0 deletions examples/RBCells/0/U
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,21 @@
internalField
uniform
1e-4 0.0 0.0
boundaryField
floor
Dirichlet
uniform
0. 0. 0.
ceiling
Dirichlet
uniform
0. 0. 0.
sideWalls
Dirichlet
uniform
0. 0. 0.
frontAndBack
Neumann
zeroGradient


19 changes: 19 additions & 0 deletions examples/RBCells/0/epsilon
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
internalField
uniform
1e-5
boundaryField
floor
Dirichlet
uniform
0.
ceiling
Dirichlet
uniform
0.
sideWalls
Dirichlet
uniform
0.
frontAndBack
Neumann
zeroGradient
20 changes: 20 additions & 0 deletions examples/RBCells/0/k
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@
internalField
uniform
1e-5
boundaryField
floor
Dirichlet
uniform
0.
ceiling
Dirichlet
uniform
0.
sideWalls
Dirichlet
uniform
0.
frontAndBack
Neumann
zeroGradient

15 changes: 15 additions & 0 deletions examples/RBCells/0/omega
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,15 @@
internalField
uniform
31.94
boundaryField
ulaz
Dirichlet
uniform
31.94
zid
Dirichlet
uniform
31.94
izlaz
Neumann
zeroGradient
5 changes: 5 additions & 0 deletions examples/RBCells/README.md
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
The example is taken from OpenFOAM (www.openfoam.com) examples of heat transfer, more specifically from the examples for buoyantBoussinesqPimpleFoam solver.

OpenFOAM/OpenFOAM-v2006/tutorials/heatTransfer/buoyantBoussinesqPimpleFoam/BenardCells


27 changes: 27 additions & 0 deletions examples/RBCells/input-piso
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,27 @@
Rayleigh-Bernard-Cells laminar TITLE
F T F LREAD,LWRITE,LTEST
T T T T F F T F F F F (LCAL(I),I=1,NPHI)
1 1 0 1 monCell,pRefCell,iPrefProcessMPoints
1e+30 1e-6 SLARGE,SORMAX
1.0 1e-3 DENSIT,VISCOS
1.0 300.0 1e-3 PRANL,TREF,BETA
T 0.0 -9.81 0.0 T LBUOY,GRAVX,GRAVY,GRAVZ,BOUSSINESQ
f 30. 0. roughWall,EROUGH,ZZERO
1.0 1.0 FACNAP,FACFLX
T T F F F LTRANSIENT,BDF,BDF2,BDF3,CN
F f f f LEVM,LASM,LLES,LDES
f f f LSGDH,LGGDH,LAFM
0 TurbModel
1e-4 0. 0. 0. 0. 0. 0. 0. UIN,VIN,WIN,TEIN,EDIN,TIN,VARTIN,CONIN
linear-f Convective scheme
Barth-Jespersen Gradient limiter
1.0 1.0 1.0 0.0 1.0 1.0 1.0 0.0 0.0 0.0 0.0 (GDS(I),I=1,NPHI)
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 0.6 0.9 0.8 (URF(I),I=1,NPHI)
3*1e-2 1e-6 7*1e-2 (SOR(I),I=1,NPHI)
20 20 20 100 10 10 10 5 5 5 5 (NSW(I),I=1,NPHI)
1000 1.0 50 1 NUMSTEP,TIMESTEP,NZAPIS,MAXIT
f f f t lstsq, lstsq_qr, lstsq_dm, gauss
1 1 NPCOR, NIGRAD
f t 2 SIMPLE,PISO,ncorr
f const_mflux
f 1.0 CoNumFix, CoNumFixValue
27 changes: 27 additions & 0 deletions examples/RBCells/input-simple
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,27 @@
Rayleigh-Bernard-Cells laminar TITLE
f t f LREAD,LWRITE,LTEST
t t t t F F T t f f F (LCAL(I),I=1,NPHI)
1 1 0 1 monCell,pRefCell,iPrefProcessMPoints
1e+30 1e-6 SLARGE,SORMAX
1.0 1e-3 DENSIT,VISCOS
1.0 300.0 1e-3 PRANL,TREF,BETA
T 0.0 -9.81 0.0 T LBUOY,GRAVX,GRAVY,GRAVZ,BOUSSINESQ
f 30. 0. roughWall,EROUGH,ZZERO
1.0 1.0 FACNAP,FACFLX
T T f f f LTRANSIENT,BDF,BDF2,BDF3,CN
F f f f LEVM,LASM,LLES,LDES
f f f LSGDH,LGGDH,LAFM
0 TurbModel
1e-4 0. 0. 0.375 18.455 0. 0. 0. UIN,VIN,WIN,TEIN,EDIN,TIN,VARTIN,CONIN
linear-f Convective scheme
Barth-Jespersen Gradient limiter
1.0 1.0 1.0 0.0 1.0 1.0 1.0 0.0 0.0 0.0 0.0 (GDS(I),I=1,NPHI)
0.7 0.7 0.7 0.3 1.0 1.0 0.7 1.0 0.6 0.9 0.8 (URF(I),I=1,NPHI)
11*1e-2 (SOR(I),I=1,NPHI)
20 20 20 100 10 10 10 5 5 5 5 (NSW(I),I=1,NPHI)
1000 1.0 50 10 NUMSTEP,TIMESTEP,NZAPIS,MAXIT
f f f t lstsq, lstsq_qr, lstsq_dm, gauss
1 1 NPCOR, NIGRAD
t f 1 SIMPLE,PISO,ncorr
f const_mflux
f 1.0 CoNumFix, CoNumFixValue
14 changes: 14 additions & 0 deletions examples/RBCells/plotResiduals
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,14 @@
set logscale y
set title "Residuals"
set ylabel 'Residual'
set xlabel 'Iteration'
set yrange[1e-5:1e+2]
plot "< cat monitor | grep 'Solving for U,' | cut -d' ' -f12 | tr -d ','" title 'Ux' with lines,\
"< cat monitor | grep 'Solving for V,' | cut -d' ' -f12 | tr -d ','" title 'Uy' with lines,\
"< cat monitor | grep 'Solving for W,' | cut -d' ' -f12 | tr -d ','" title 'Uz' with lines,\
"< cat monitor | grep 'Solving for p' | cut -d' ' -f12 | tr -d ','" title 'p' with lines, \
"< cat monitor | grep 'Solving for Temp' | cut -d' ' -f12 | tr -d ','" title 'T' with lines, \
"< cat monitor | grep 'Solving for k' | cut -d' ' -f12 | tr -d ','" title 'k' with lines, \
"< cat monitor | grep 'Solving for Omega' | cut -d' ' -f12 | tr -d ','" title 'omega' with lines
pause 5
reread
5 changes: 5 additions & 0 deletions examples/RBCells/polyMesh/boundary
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,5 @@
# bcName bcType nFaces startFace
floor wall 90 1700
ceiling wall 90 1790
sideWalls wall 20 1880
frontAndBack symmetry 1800 1900
51 changes: 51 additions & 0 deletions examples/RBCells/polyMesh/boundary-OF
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,51 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1912 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format binary;
class polyBoundaryMesh;
arch "LSB;label=32;scalar=64";
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

4
(
floor
{
type wall;
inGroups 1(wall);
nFaces 90;
startFace 1700;
}
ceiling
{
type wall;
inGroups 1(wall);
nFaces 90;
startFace 1790;
}
sideWalls
{
type wall;
inGroups 1(wall);
nFaces 20;
startFace 1880;
}
frontAndBack
{
type empty;
inGroups 1(empty);
nFaces 1800;
startFace 1900;
}
)

// ************************************************************************* //
Loading

0 comments on commit 6f4d0b9

Please sign in to comment.