Skip to content

Tuning the perpendicular-flap/solid-solids4foam case #893

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
MakisH wants to merge 1 commit into
base: develop
Choose a base branch
from
Draft

Tuning the perpendicular-flap/solid-solids4foam case #893

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
24 changes: 24 additions & 0 deletions perpendicular-flap/solid-solids4foam/constant/mechanicalProperties.withDefaultValues
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,24 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2512 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
arch "LSB;label=32;scalar=64";
class dictionary;
location "constant";
object mechanicalProperties.withDefaultValues;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

planeStress no;

mechanical 1 ( rubber { type linearElastic ; rho rho [ 1 -3 0 0 0 0 0 ] 3000 ; E E [ 1 -1 -2 0 0 0 0 ] 4e+06 ; nu nu [ 0 0 0 0 0 0 0 ] 0.3 ; solvePressureEqn false ; pressureSmoothingScaleFactor 100 ; } ; );


// ************************************************************************* //
12 changes: 1 addition & 11 deletions perpendicular-flap/solid-solids4foam/constant/solidProperties
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,16 +10,6 @@ solidModel linearGeometryTotalDisplacement;

linearGeometryTotalDisplacementCoeffs
{
// Maximum number of momentum correctors
nCorrectors 1000;

// Solution tolerance for displacement
solutionTolerance 1e-08;

// Alternative solution tolerance for displacement
alternativeTolerance 1e-08;

// Write frequency for the residuals
infoFrequency 100;
solutionAlgorithm PETScSNES;
}

50 changes: 50 additions & 0 deletions perpendicular-flap/solid-solids4foam/constant/solidProperties.withDefaultValues
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,50 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2512 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
arch "LSB;label=32;scalar=64";
class dictionary;
location "constant";
object solidProperties.withDefaultValues;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solidModel linearGeometryTotalDisplacement;

linearGeometryTotalDisplacementCoeffs
{
solutionAlgorithm PETScSNES;
dampingCoeff dampingCoeff [ 0 0 -1 0 0 0 0 ] 0;
solutionTolerance 1e-06;
alternativeTolerance 1e-07;
materialTolerance 1e-05;
infoFrequency 100;
nCorrectors 10000;
writeResidualField false;
relaxationMethod fixed;
restart false;
residualFile false;
stabilisation
{
momentum
{
type diffStencilLaplacian;
scaleFactor 0.1;
}
pressure
{
type diffStencilLaplacian;
scaleFactor 0.1;
}
}
}


// ************************************************************************* //
4 changes: 2 additions & 2 deletions perpendicular-flap/solid-solids4foam/system/fvSchemes
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ ddtSchemes

gradSchemes
{
default leastSquaresS4f 0;
default leastSquaresS4f;
}

divSchemes
Expand All @@ -28,7 +28,7 @@ divSchemes

laplacianSchemes
{
default none;
default Gauss linear corrected;
laplacian(DD,D) Gauss linear corrected;
laplacian(DDD,DD) Gauss linear corrected;
}
Expand Down
49 changes: 28 additions & 21 deletions perpendicular-flap/solid-solids4foam/system/fvSolution
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,3 +1,12 @@
/*--------------------------------*- C++ -*----------------------------------*\
| solids4foam: solid mechanics and fluid-solid interaction simulations |
| Version: v2.3 |
| Web: https://solids4foam.github.io |
| Disclaimer: This offering is not approved or endorsed by OpenCFD Limited, |
| producer and distributor of the OpenFOAM software via |
| www.openfoam.com, and owner of the OPENFOAM® and OpenCFD® |
| trade marks. |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
Expand All @@ -6,33 +15,31 @@ FoamFile
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
"D|DD"
D
{
solver PCG;
preconditioner FDIC;
tolerance 1e-09;
relTol 0.1;
}
}
solver petsc;

relaxationFactors
{
// Under-relaxating the D equation by 0.99-0.9999 can improve convergence in
// some cases, in particular when there are solidContact boundaries
equations
{
//D 0.999;
}
options
{
snes_type newtonls;
snes_monitor;
snes_converged_reason;
snes_mf;
snes_mf_operator;

// Under-relaxating the D field by 0.1-0.9 can improve convergence in some
// cases
fields
{
//D 0.9;
ksp_type lgmres;
ksp_gmres_restart "100";
ksp_converged_reason;

pc_type lu;
snes_lag_preconditioner_persists true;
snes_lag_preconditioner "-2";
}
}
}


// ************************************************************************* //