Merge branch 'HrvojeJasak'
This commit is contained in:
commit
0b7bacfbb3
25 changed files with 110 additions and 83 deletions
|
@ -4,7 +4,14 @@
|
||||||
U = UEqn.H()/UEqn.A();
|
U = UEqn.H()/UEqn.A();
|
||||||
|
|
||||||
phi = fvc::interpolate(U) & mesh.Sf();
|
phi = fvc::interpolate(U) & mesh.Sf();
|
||||||
surfaceScalarField phiU = phi;
|
|
||||||
|
// Bug fix: must change name of phi on copy to keep objectRegistry happy
|
||||||
|
// HJ, 7/Nov/2010
|
||||||
|
surfaceScalarField phiU
|
||||||
|
(
|
||||||
|
"phiU",
|
||||||
|
phi
|
||||||
|
);
|
||||||
|
|
||||||
surfaceScalarField phip = fvc::interpolate(psiByRho)*phi;
|
surfaceScalarField phip = fvc::interpolate(psiByRho)*phi;
|
||||||
|
|
||||||
|
|
|
@ -202,7 +202,14 @@ void Foam::kineticTheoryModel::solve()
|
||||||
|
|
||||||
word scheme("div(phi,Theta)");
|
word scheme("div(phi,Theta)");
|
||||||
|
|
||||||
volScalarField alpha = alpha_;
|
// Bug fix: copying field alpha without mane change vreaks objectRegistry
|
||||||
|
// HJ, 7/Nov/2010
|
||||||
|
volScalarField alpha
|
||||||
|
(
|
||||||
|
"alphaMinMax",
|
||||||
|
alpha_
|
||||||
|
);
|
||||||
|
|
||||||
alpha.max(1.0e-6);
|
alpha.max(1.0e-6);
|
||||||
const scalar sqrtPi = sqrt(mathematicalConstant::pi);
|
const scalar sqrtPi = sqrt(mathematicalConstant::pi);
|
||||||
|
|
||||||
|
|
|
@ -51,7 +51,7 @@ Foam::BlockBiCGStabSolver<Type>::BlockBiCGStabSolver
|
||||||
BlockLduPrecon<Type>::New
|
BlockLduPrecon<Type>::New
|
||||||
(
|
(
|
||||||
matrix,
|
matrix,
|
||||||
this->dict().subDict("preconditioner")
|
this->dict()
|
||||||
)
|
)
|
||||||
)
|
)
|
||||||
{}
|
{}
|
||||||
|
|
|
@ -51,7 +51,7 @@ Foam::BlockCGSolver<Type>::BlockCGSolver
|
||||||
BlockLduPrecon<Type>::New
|
BlockLduPrecon<Type>::New
|
||||||
(
|
(
|
||||||
matrix,
|
matrix,
|
||||||
this->dict().subDict("preconditioner")
|
this->dict()
|
||||||
)
|
)
|
||||||
)
|
)
|
||||||
{}
|
{}
|
||||||
|
|
|
@ -85,7 +85,7 @@ Foam::BlockGMRESSolver<Type>::BlockGMRESSolver
|
||||||
BlockLduPrecon<Type>::New
|
BlockLduPrecon<Type>::New
|
||||||
(
|
(
|
||||||
matrix,
|
matrix,
|
||||||
this->dict().subDict("preconditioner")
|
this->dict()
|
||||||
)
|
)
|
||||||
),
|
),
|
||||||
nDirs_(readLabel(this->dict().lookup("nDirections")))
|
nDirs_(readLabel(this->dict().lookup("nDirections")))
|
||||||
|
|
|
@ -1,8 +1,8 @@
|
||||||
/*--------------------------------*- C++ -*----------------------------------*\
|
/*--------------------------------*- C++ -*----------------------------------*\
|
||||||
| ========= | |
|
| ========= | |
|
||||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
|
||||||
| \\ / O peration | Version: 1.5-dev |
|
| \\ / O peration | Version: 1.6-ext |
|
||||||
| \\ / A nd | Web: http://www.OpenFOAM.org |
|
| \\ / A nd | Web: www.extend-project.de |
|
||||||
| \\/ M anipulation | |
|
| \\/ M anipulation | |
|
||||||
\*---------------------------------------------------------------------------*/
|
\*---------------------------------------------------------------------------*/
|
||||||
FoamFile
|
FoamFile
|
||||||
|
|
|
@ -1,8 +1,8 @@
|
||||||
/*--------------------------------*- C++ -*----------------------------------*\
|
/*--------------------------------*- C++ -*----------------------------------*\
|
||||||
| ========= | |
|
| ========= | |
|
||||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
|
||||||
| \\ / O peration | Version: dev |
|
| \\ / O peration | Version: 1.6-ext |
|
||||||
| \\ / A nd | Web: http://www.OpenFOAM.org |
|
| \\ / A nd | Web: www.extend-project.de |
|
||||||
| \\/ M anipulation | |
|
| \\/ M anipulation | |
|
||||||
\*---------------------------------------------------------------------------*/
|
\*---------------------------------------------------------------------------*/
|
||||||
FoamFile
|
FoamFile
|
||||||
|
@ -21,25 +21,25 @@ FoamFile
|
||||||
{
|
{
|
||||||
type patch;
|
type patch;
|
||||||
nFaces 51;
|
nFaces 51;
|
||||||
startFace 15151;
|
startFace 23114;
|
||||||
}
|
}
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type patch;
|
type patch;
|
||||||
nFaces 51;
|
nFaces 51;
|
||||||
startFace 15202;
|
startFace 23165;
|
||||||
}
|
}
|
||||||
walls
|
walls
|
||||||
{
|
{
|
||||||
type wall;
|
type wall;
|
||||||
nFaces 436;
|
nFaces 512;
|
||||||
startFace 15253;
|
startFace 23216;
|
||||||
}
|
}
|
||||||
frontBack
|
frontBack
|
||||||
{
|
{
|
||||||
type empty;
|
type empty;
|
||||||
nFaces 15420;
|
nFaces 23340;
|
||||||
startFace 15689;
|
startFace 23728;
|
||||||
}
|
}
|
||||||
)
|
)
|
||||||
|
|
||||||
|
|
|
@ -30,6 +30,7 @@ boundaryField
|
||||||
top
|
top
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 1.0e-8;
|
inletValue uniform 1.0e-8;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
top
|
top
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform (0 0 0);
|
inletValue uniform (0 0 0);
|
||||||
value uniform ( 0 0 0 );
|
value uniform ( 0 0 0 );
|
||||||
}
|
}
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
top
|
top
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phib;
|
||||||
value uniform (0 0 0);
|
value uniform (0 0 0);
|
||||||
inletValue uniform ( 0 0 0 );
|
inletValue uniform ( 0 0 0 );
|
||||||
}
|
}
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
top
|
top
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 0;
|
inletValue uniform 0;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
top
|
top
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 10.0;
|
inletValue uniform 10.0;
|
||||||
value uniform 10.0;
|
value uniform 10.0;
|
||||||
}
|
}
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
top
|
top
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 1.0;
|
inletValue uniform 1.0;
|
||||||
value uniform 1.0;
|
value uniform 1.0;
|
||||||
}
|
}
|
||||||
|
|
|
@ -1,8 +1,8 @@
|
||||||
/*--------------------------------*- C++ -*----------------------------------*\
|
/*--------------------------------*- C++ -*----------------------------------*\
|
||||||
| ========= | |
|
| ========= | |
|
||||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
|
||||||
| \\ / O peration | Version: 1.5 |
|
| \\ / O peration | Version: 1.6-ext |
|
||||||
| \\ / A nd | Web: http://www.OpenFOAM.org |
|
| \\ / A nd | Web: www.extend-project.de |
|
||||||
| \\/ M anipulation | |
|
| \\/ M anipulation | |
|
||||||
\*---------------------------------------------------------------------------*/
|
\*---------------------------------------------------------------------------*/
|
||||||
FoamFile
|
FoamFile
|
||||||
|
@ -10,6 +10,7 @@ FoamFile
|
||||||
version 2.0;
|
version 2.0;
|
||||||
format binary;
|
format binary;
|
||||||
class polyBoundaryMesh;
|
class polyBoundaryMesh;
|
||||||
|
location "constant/polyMesh";
|
||||||
object boundary;
|
object boundary;
|
||||||
}
|
}
|
||||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 10.0;
|
inletValue uniform 10.0;
|
||||||
value uniform 10.0;
|
value uniform 10.0;
|
||||||
}
|
}
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 1.0;
|
inletValue uniform 1.0;
|
||||||
value uniform 1.0;
|
value uniform 1.0;
|
||||||
}
|
}
|
||||||
|
|
|
@ -1,8 +1,8 @@
|
||||||
/*--------------------------------*- C++ -*----------------------------------*\
|
/*--------------------------------*- C++ -*----------------------------------*\
|
||||||
| ========= | |
|
| ========= | |
|
||||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
|
||||||
| \\ / O peration | Version: 1.5 |
|
| \\ / O peration | Version: 1.6-ext |
|
||||||
| \\ / A nd | Web: http://www.OpenFOAM.org |
|
| \\ / A nd | Web: www.extend-project.de |
|
||||||
| \\/ M anipulation | |
|
| \\/ M anipulation | |
|
||||||
\*---------------------------------------------------------------------------*/
|
\*---------------------------------------------------------------------------*/
|
||||||
FoamFile
|
FoamFile
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 1.0e-7;
|
inletValue uniform 1.0e-7;
|
||||||
value uniform 1.0e-7;
|
value uniform 1.0e-7;
|
||||||
}
|
}
|
||||||
|
|
|
@ -1908,6 +1908,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform (0 0 0);
|
inletValue uniform (0 0 0);
|
||||||
value nonuniform List<vector>
|
value nonuniform List<vector>
|
||||||
25
|
25
|
||||||
|
|
|
@ -1908,6 +1908,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phib;
|
||||||
inletValue uniform (0 0 0);
|
inletValue uniform (0 0 0);
|
||||||
value nonuniform List<vector>
|
value nonuniform List<vector>
|
||||||
25
|
25
|
||||||
|
|
|
@ -1908,6 +1908,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 1;
|
inletValue uniform 1;
|
||||||
value uniform 1;
|
value uniform 1;
|
||||||
}
|
}
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 0.1;
|
inletValue uniform 0.1;
|
||||||
value uniform 0.1;
|
value uniform 0.1;
|
||||||
}
|
}
|
||||||
|
|
|
@ -29,6 +29,7 @@ boundaryField
|
||||||
outlet
|
outlet
|
||||||
{
|
{
|
||||||
type inletOutlet;
|
type inletOutlet;
|
||||||
|
phi phia;
|
||||||
inletValue uniform 1e-8;
|
inletValue uniform 1e-8;
|
||||||
value uniform 1e-8;
|
value uniform 1e-8;
|
||||||
}
|
}
|
||||||
|
|
|
@ -1,8 +1,8 @@
|
||||||
/*--------------------------------*- C++ -*----------------------------------*\
|
/*--------------------------------*- C++ -*----------------------------------*\
|
||||||
| ========= | |
|
| ========= | |
|
||||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
|
||||||
| \\ / O peration | Version: 1.5 |
|
| \\ / O peration | Version: 1.6-ext |
|
||||||
| \\ / A nd | Web: http://www.OpenFOAM.org |
|
| \\ / A nd | Web: www.extend-project.de |
|
||||||
| \\/ M anipulation | |
|
| \\/ M anipulation | |
|
||||||
\*---------------------------------------------------------------------------*/
|
\*---------------------------------------------------------------------------*/
|
||||||
FoamFile
|
FoamFile
|
||||||
|
@ -10,39 +10,37 @@ FoamFile
|
||||||
version 2.0;
|
version 2.0;
|
||||||
format ascii;
|
format ascii;
|
||||||
class polyBoundaryMesh;
|
class polyBoundaryMesh;
|
||||||
|
location "constant/polyMesh";
|
||||||
object boundary;
|
object boundary;
|
||||||
}
|
}
|
||||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
4
|
4
|
||||||
(
|
(
|
||||||
inlet
|
inlet
|
||||||
{
|
{
|
||||||
type patch;
|
type patch;
|
||||||
nFaces 25;
|
nFaces 25;
|
||||||
startFace 3650;
|
startFace 3650;
|
||||||
}
|
}
|
||||||
|
outlet
|
||||||
outlet
|
{
|
||||||
{
|
|
||||||
type patch;
|
type patch;
|
||||||
nFaces 25;
|
nFaces 25;
|
||||||
startFace 3675;
|
startFace 3675;
|
||||||
}
|
}
|
||||||
|
walls
|
||||||
walls
|
{
|
||||||
{
|
|
||||||
type wall;
|
type wall;
|
||||||
nFaces 150;
|
nFaces 150;
|
||||||
startFace 3700;
|
startFace 3700;
|
||||||
}
|
}
|
||||||
|
defaultFaces
|
||||||
defaultFaces
|
{
|
||||||
{
|
|
||||||
type empty;
|
type empty;
|
||||||
nFaces 3750;
|
nFaces 3750;
|
||||||
startFace 3850;
|
startFace 3850;
|
||||||
}
|
}
|
||||||
)
|
)
|
||||||
|
|
||||||
// ************************************************************************* //
|
// ************************************************************************* //
|
||||||
|
|
Reference in a new issue