d35773f4fb
Uses adjoint method to block regions of the domain causing total pressure loss, e.g. recirculation zones. pitzDaily tutorial case supplied.
77 lines
1.9 KiB
C++
77 lines
1.9 KiB
C++
/*--------------------------------*- C++ -*----------------------------------*\
|
|
| ========= | |
|
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
|
| \\ / O peration | Version: dev |
|
|
| \\ / A nd | Web: www.OpenFOAM.org |
|
|
| \\/ M anipulation | |
|
|
\*---------------------------------------------------------------------------*/
|
|
FoamFile
|
|
{
|
|
version 2.0;
|
|
format ascii;
|
|
class dictionary;
|
|
location "system";
|
|
object fvSchemes;
|
|
}
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
ddtSchemes
|
|
{
|
|
default steadyState;
|
|
}
|
|
|
|
gradSchemes
|
|
{
|
|
default Gauss linear;
|
|
|
|
grad(p) Gauss linear;
|
|
grad(U) Gauss linear;
|
|
}
|
|
|
|
divSchemes
|
|
{
|
|
default none;
|
|
|
|
div(phi,U) Gauss upwind;
|
|
div(phi,k) Gauss upwind;
|
|
div(phi,epsilon) Gauss upwind;
|
|
div((nuEff*dev(grad(U).T()))) Gauss linear;
|
|
|
|
div(-phi,Ua) Gauss upwind;
|
|
div((nuEff*dev(grad(Ua).T()))) Gauss linear;
|
|
}
|
|
|
|
laplacianSchemes
|
|
{
|
|
default Gauss linear corrected;
|
|
|
|
laplacian(nuEff,U) Gauss linear corrected;
|
|
laplacian(1,p) Gauss linear corrected;
|
|
laplacian(1|A(U),p) Gauss linear corrected;
|
|
laplacian(DkEff,k) Gauss linear corrected;
|
|
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
|
|
|
|
laplacian(nuEff,Ua) Gauss linear corrected;
|
|
laplacian(1|A(Ua),pa) Gauss linear corrected;
|
|
}
|
|
|
|
interpolationSchemes
|
|
{
|
|
default linear;
|
|
interpolate(U) linear;
|
|
}
|
|
|
|
snGradSchemes
|
|
{
|
|
default corrected;
|
|
}
|
|
|
|
fluxRequired
|
|
{
|
|
default no;
|
|
p;
|
|
pa;
|
|
}
|
|
|
|
// ************************************************************************* //
|