FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

motionDirectionalDiffusivity.C

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software: you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by
00013     the Free Software Foundation, either version 3 of the License, or
00014     (at your option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
00023 
00024 \*---------------------------------------------------------------------------*/
00025 
00026 #include "motionDirectionalDiffusivity.H"
00027 #include <finiteVolume/surfaceInterpolate.H>
00028 #include <finiteVolume/zeroGradientFvPatchFields.H>
00029 #include <OpenFOAM/addToRunTimeSelectionTable.H>
00030 
00031 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
00032 
00033 namespace Foam
00034 {
00035     defineTypeNameAndDebug(motionDirectionalDiffusivity, 0);
00036 
00037     addToRunTimeSelectionTable
00038     (
00039         motionDiffusivity,
00040         motionDirectionalDiffusivity,
00041         Istream
00042     );
00043 }
00044 
00045 
00046 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
00047 
00048 Foam::motionDirectionalDiffusivity::motionDirectionalDiffusivity
00049 (
00050     const fvMotionSolver& mSolver,
00051     Istream& mdData
00052 )
00053 :
00054     uniformDiffusivity(mSolver, mdData),
00055     diffusivityVector_(mdData)
00056 {}
00057 
00058 
00059 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
00060 
00061 Foam::motionDirectionalDiffusivity::~motionDirectionalDiffusivity()
00062 {}
00063 
00064 
00065 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
00066 
00067 void Foam::motionDirectionalDiffusivity::correct()
00068 {
00069     const fvMesh& mesh = mSolver().mesh();
00070 
00071     static bool first = true;
00072 
00073     if (!first)
00074     {
00075         const volVectorField& cellMotionU = 
00076             mesh.lookupObject<volVectorField>("cellMotionU");
00077 
00078         volVectorField D
00079         (
00080             IOobject
00081             (
00082                 "D",
00083                 mesh.time().timeName(),
00084                 mesh
00085             ),
00086             diffusivityVector_.y()*vector::one
00087           + (diffusivityVector_.x() - diffusivityVector_.y())*cellMotionU
00088            /(mag(cellMotionU) + dimensionedScalar("small", dimVelocity, SMALL)),
00089             zeroGradientFvPatchVectorField::typeName
00090         );
00091         D.correctBoundaryConditions();
00092 
00093         surfaceVectorField n = mesh.Sf()/mesh.magSf();
00094         faceDiffusivity_ == (n & cmptMultiply(fvc::interpolate(D), n));
00095     }
00096     else
00097     {
00098         first = false;
00099         const_cast<fvMotionSolver&>(mSolver()).solve();
00100         correct();
00101     }
00102 }
00103 
00104 
00105 // ************************ vim: set sw=4 sts=4 et: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines