From 950d4a078fa77ef328a5f0268c3e56839a29eb51 Mon Sep 17 00:00:00 2001 From: Henry Weller <http://cfd.direct> Date: Thu, 9 Jun 2016 18:09:27 +0100 Subject: [PATCH] stressComponents: utility replaced by the 'R' and 'components' functionObject used with the '-postProcess' option --- .../stressField/stressComponents/Make/files | 3 - .../stressField/stressComponents/Make/options | 10 - .../stressComponents/stressComponents.C | 223 ------------------ 3 files changed, 236 deletions(-) delete mode 100644 applications/utilities/postProcessing/stressField/stressComponents/Make/files delete mode 100644 applications/utilities/postProcessing/stressField/stressComponents/Make/options delete mode 100644 applications/utilities/postProcessing/stressField/stressComponents/stressComponents.C diff --git a/applications/utilities/postProcessing/stressField/stressComponents/Make/files b/applications/utilities/postProcessing/stressField/stressComponents/Make/files deleted file mode 100644 index db34ff3cac..0000000000 --- a/applications/utilities/postProcessing/stressField/stressComponents/Make/files +++ /dev/null @@ -1,3 +0,0 @@ -stressComponents.C - -EXE = $(FOAM_APPBIN)/stressComponents diff --git a/applications/utilities/postProcessing/stressField/stressComponents/Make/options b/applications/utilities/postProcessing/stressField/stressComponents/Make/options deleted file mode 100644 index e5d2760627..0000000000 --- a/applications/utilities/postProcessing/stressField/stressComponents/Make/options +++ /dev/null @@ -1,10 +0,0 @@ -EXE_INC = \ - -I$(LIB_SRC)/transportModels \ - -I$(LIB_SRC)/finiteVolume/lnInclude \ - -I$(LIB_SRC)/meshTools/lnInclude - -EXE_LIBS = \ - -lgenericPatchFields \ - -lincompressibleTransportModels \ - -lfiniteVolume \ - -lmeshTools diff --git a/applications/utilities/postProcessing/stressField/stressComponents/stressComponents.C b/applications/utilities/postProcessing/stressField/stressComponents/stressComponents.C deleted file mode 100644 index 35161c0dff..0000000000 --- a/applications/utilities/postProcessing/stressField/stressComponents/stressComponents.C +++ /dev/null @@ -1,223 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | - \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation - \\/ M anipulation | -------------------------------------------------------------------------------- -License - This file is part of OpenFOAM. - - OpenFOAM is free software: you can redistribute it and/or modify it - under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - OpenFOAM is distributed in the hope that it will be useful, but WITHOUT - ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or - FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License - for more details. - - You should have received a copy of the GNU General Public License - along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>. - -Application - stressComponents - -Description - Calculates and writes the scalar fields of the six components of the stress - tensor sigma for each time. - -\*---------------------------------------------------------------------------*/ - -#include "fvCFD.H" -#include "incompressible/singlePhaseTransportModel/singlePhaseTransportModel.H" -#include "zeroGradientFvPatchFields.H" - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -int main(int argc, char *argv[]) -{ - timeSelector::addOptions(); - - #include "setRootCase.H" - #include "createTime.H" - - instantList timeDirs = timeSelector::select0(runTime, args); - - #include "createMesh.H" - - forAll(timeDirs, timeI) - { - runTime.setTime(timeDirs[timeI], timeI); - - Info<< "Time = " << runTime.timeName() << endl; - - IOobject Uheader - ( - "U", - runTime.timeName(), - mesh, - IOobject::MUST_READ - ); - - // Check U exists - if (Uheader.headerOk()) - { - mesh.readUpdate(); - - Info<< " Reading U" << endl; - volVectorField U(Uheader, mesh); - - #include "createPhi.H" - - singlePhaseTransportModel laminarTransport(U, phi); - - volSymmTensorField sigma - ( - IOobject - ( - "sigma", - runTime.timeName(), - mesh, - IOobject::NO_READ, - IOobject::AUTO_WRITE - ), - laminarTransport.nu()*2*dev(symm(fvc::grad(U))) - ); - - - volScalarField sigmaxx - ( - IOobject - ( - "sigmaxx", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - sigma.component(symmTensor::XX) - ); - sigmaxx.write(); - - volScalarField sigmayy - ( - IOobject - ( - "sigmayy", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - sigma.component(symmTensor::YY) - ); - sigmayy.write(); - - volScalarField sigmazz - ( - IOobject - ( - "sigmazz", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - sigma.component(symmTensor::ZZ) - ); - sigmazz.write(); - - volScalarField sigmaxy - ( - IOobject - ( - "sigmaxy", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - sigma.component(symmTensor::XY) - ); - sigmaxy.write(); - - volScalarField sigmaxz - ( - IOobject - ( - "sigmaxz", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - sigma.component(symmTensor::XZ) - ); - sigmaxz.write(); - - volScalarField sigmayz - ( - IOobject - ( - "sigmayz", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - sigma.component(symmTensor::YZ) - ); - sigmayz.write(); - - volVectorField Ub - ( - IOobject - ( - "Ub", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - U, - zeroGradientFvPatchVectorField::typeName - ); - Ub.correctBoundaryConditions(); - Ub.write(); - - volScalarField sigmaUn - ( - IOobject - ( - "sigmaUn", - runTime.timeName(), - mesh, - IOobject::NO_READ - ), - 0.0*sigma.component(symmTensor::YZ) - ); - - volScalarField::Boundary& sigmaUnBf = - sigmaUn.boundaryFieldRef(); - - forAll(sigmaUn.boundaryField(), patchi) - { - sigmaUnBf[patchi] = - ( - mesh.boundary()[patchi].nf() - & sigma.boundaryField()[patchi] - )().component(vector::X); - } - - sigmaUn.write(); - } - else - { - Info<< " No U" << endl; - } - - Info<< endl; - } - - Info<< "End" << endl; - - return 0; -} - - -// ************************************************************************* // -- GitLab