GIT: revert: revert bad commit b6c4e144c4

This commit is contained in:
mattijs
2012-09-20 15:19:25 +01:00
parent 6835300c3c
commit cbdba6b471
10 changed files with 21 additions and 73 deletions

View File

@ -81,14 +81,8 @@ Foam::displacementSBRStressFvMotionSolver::displacementSBRStressFvMotionSolver
diffusivityPtr_ diffusivityPtr_
( (
motionDiffusivity::New(fvMesh_, coeffDict().lookup("diffusivity")) motionDiffusivity::New(fvMesh_, coeffDict().lookup("diffusivity"))
), )
solveOnPoints0_(coeffDict().lookupOrDefault("solveOnPoints0", false)) {}
{
if (solveOnPoints0_)
{
Info<< type() << " : solving on points0" << endl;
}
}
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
@ -129,18 +123,11 @@ void Foam::displacementSBRStressFvMotionSolver::solve()
diffusivityPtr_->correct(); diffusivityPtr_->correct();
pointDisplacement_.boundaryField().updateCoeffs(); pointDisplacement_.boundaryField().updateCoeffs();
pointField oldPoints;
if (solveOnPoints0_)
{
oldPoints = fvMesh_.points();
movePoints(points0());
}
surfaceScalarField Df(diffusivityPtr_->operator()()); surfaceScalarField Df(diffusivityPtr_->operator()());
volTensorField gradCd(fvc::grad(cellDisplacement_)); volTensorField gradCd(fvc::grad(cellDisplacement_));
tmp<fvMatrix<vector> > laplacianDf Foam::solve
( (
fvm::laplacian fvm::laplacian
( (
@ -148,16 +135,7 @@ void Foam::displacementSBRStressFvMotionSolver::solve()
cellDisplacement_, cellDisplacement_,
"laplacian(diffusivity,cellDisplacement)" "laplacian(diffusivity,cellDisplacement)"
) )
);
// if (solveOnPoints0_)
// {
// movePoints(oldPoints);
// }
Foam::solve
(
laplacianDf
+ fvc::div + fvc::div
( (
Df Df
@ -195,12 +173,6 @@ void Foam::displacementSBRStressFvMotionSolver::solve()
) )
*/ */
); );
if (solveOnPoints0_)
{
movePoints(points0());
}
} }

View File

@ -64,8 +64,6 @@ class displacementSBRStressFvMotionSolver
//- Diffusivity used to control the motion //- Diffusivity used to control the motion
autoPtr<motionDiffusivity> diffusivityPtr_; autoPtr<motionDiffusivity> diffusivityPtr_;
//- Solve on base mesh
Switch solveOnPoints0_;
// Private Member Functions // Private Member Functions

View File

@ -17,7 +17,7 @@ FoamFile
application rhoPimplecFoam; application rhoPimplecFoam;
startFrom latestTime; //startTime; startFrom startTime;
startTime 0; startTime 0;

View File

@ -17,7 +17,7 @@ FoamFile
application rhoSimplecFoam; application rhoSimplecFoam;
startFrom latestTime; //startTime; startFrom startTime;
startTime 0; startTime 0;

View File

@ -35,18 +35,6 @@ geometry
radius 4; radius 4;
} }
// // Define floor just so we have feature edges behave correctly
// floor
// {
// type searchablePlane;
// planeType pointAndNormal;
// pointAndNormalDict
// {
// basePoint (0 0 0);
// normalVector (0 0 1);
// }
// }
box1 box1
{ {
type searchableBox; type searchableBox;
@ -199,11 +187,6 @@ castellatedMeshControls
// Surface-wise min and max refinement level // Surface-wise min and max refinement level
level (1 1); level (1 1);
} }
// floor
// {
// // Surface-wise min and max refinement level
// level (1 1);
// }
} }
// Resolve sharp angles on fridges // Resolve sharp angles on fridges
@ -252,10 +235,7 @@ snapControls
{ {
//- Number of patch smoothing iterations before finding correspondence //- Number of patch smoothing iterations before finding correspondence
// to surface // to surface
// nSmoothPatch 3; nSmoothPatch 3;
//nSmoothPatch 0; //! Only since blockMesh aligned with surfaces it is
// actually beneficial to run with 0 smoothing.
nSmoothPatch 1;
//- Relative distance for points to be attracted by surface feature point //- Relative distance for points to be attracted by surface feature point
// or edge. True distance is this factor times local // or edge. True distance is this factor times local
@ -314,7 +294,7 @@ addLayersControls
// Relative to undistorted size of cell outside layer. // Relative to undistorted size of cell outside layer.
// is the thickness of the layer furthest away from the wall. // is the thickness of the layer furthest away from the wall.
// See relativeSizes parameter. // See relativeSizes parameter.
finalLayerThickness 1; finalLayerThickness 0.5;
//- Minimum thickness of cell layer. If for any reason layer //- Minimum thickness of cell layer. If for any reason layer
// cannot be above minThickness do not add layer. // cannot be above minThickness do not add layer.
@ -393,7 +373,7 @@ meshQualityControls
// 1 : write intermediate meshes // 1 : write intermediate meshes
// 2 : write volScalarField with cellLevel for postprocessing // 2 : write volScalarField with cellLevel for postprocessing
// 4 : write current intersections as .obj files // 4 : write current intersections as .obj files
debug 7; debug 0;
// Merge tolerance. Is fraction of overall bounding box of initial mesh. // Merge tolerance. Is fraction of overall bounding box of initial mesh.

View File

@ -46,7 +46,6 @@ boundary
( (
(3 7 6 2) (3 7 6 2)
); );
inGroups (groupWall);
} }
fixedWalls fixedWalls
{ {
@ -57,7 +56,6 @@ boundary
(2 6 5 1) (2 6 5 1)
(1 5 4 0) (1 5 4 0)
); );
inGroups (groupWall);
} }
frontAndBack frontAndBack
{ {

View File

@ -17,7 +17,7 @@ FoamFile
application porousSimpleFoam; application porousSimpleFoam;
startFrom latestTime; //startTime; startFrom startTime;
startTime 0; startTime 0;

View File

@ -17,23 +17,23 @@ FoamFile
application moveDynamicMesh; application moveDynamicMesh;
startFrom latestTime; //startTime; startFrom startTime;
startTime 0; startTime 0;
stopAt endTime; stopAt endTime;
endTime 30; endTime 25;
deltaT 1; deltaT 1;
writeControl timeStep; writeControl timeStep;
writeInterval 1; writeInterval 5;
purgeWrite 0; purgeWrite 0;
writeFormat ascii;//binary; writeFormat binary;
writePrecision 6; writePrecision 6;

View File

@ -17,7 +17,7 @@ FoamFile
application icoFoam; application icoFoam;
startFrom latestTime; startFrom startTime;
startTime 0; startTime 0;

View File

@ -16,7 +16,7 @@ FoamFile
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Which of the steps to run // Which of the steps to run
castellatedMesh false; castellatedMesh true;
snap true; snap true;
addLayers false; addLayers false;
@ -85,10 +85,10 @@ castellatedMeshControls
// This is a featureEdgeMesh, read from constant/triSurface for now. // This is a featureEdgeMesh, read from constant/triSurface for now.
features features
( (
// { {
// file "flange.eMesh"; file "flange.eMesh";
// level 0; level 0;
// } }
); );
@ -182,10 +182,10 @@ snapControls
nFeatureSnapIter 10; nFeatureSnapIter 10;
//- Detect (geometric) features by sampling the surface //- Detect (geometric) features by sampling the surface
implicitFeatureSnap true; implicitFeatureSnap false;
//- Use castellatedMeshControls::features //- Use castellatedMeshControls::features
explicitFeatureSnap false; explicitFeatureSnap true;
} }