diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/TEqn.H b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/TEqn.H index 708d968d6c..65c92dab51 100644 --- a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/TEqn.H +++ b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/TEqn.H @@ -1,14 +1,14 @@ { - kappat = turbulence->nut()/Prt; - kappat.correctBoundaryConditions(); + alphat = turbulence->nut()/Prt; + alphat.correctBoundaryConditions(); - volScalarField kappaEff("kappaEff", turbulence->nu()/Pr + kappat); + volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat); fvScalarMatrix TEqn ( fvm::ddt(T) + fvm::div(phi, T) - - fvm::laplacian(kappaEff, T) + - fvm::laplacian(alphaEff, T) == radiation->ST(rhoCpRef, T) ); diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H index 342af079d8..3609fe66fa 100644 --- a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H +++ b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H @@ -65,12 +65,12 @@ ); // kinematic turbulent thermal thermal conductivity m2/s - Info<< "Reading field kappat\n" << endl; - volScalarField kappat + Info<< "Reading field alphat\n" << endl; + volScalarField alphat ( IOobject ( - "kappat", + "alphat", runTime.timeName(), mesh, IOobject::MUST_READ, @@ -116,4 +116,3 @@ pRefValue - getRefCellValue(p, pRefCell) ); } - diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H index 76adf004cb..a0c5d24b2b 100644 --- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H +++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/TEqn.H @@ -1,13 +1,13 @@ { - kappat = turbulence->nut()/Prt; - kappat.correctBoundaryConditions(); + alphat = turbulence->nut()/Prt; + alphat.correctBoundaryConditions(); - volScalarField kappaEff("kappaEff", turbulence->nu()/Pr + kappat); + volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat); fvScalarMatrix TEqn ( fvm::div(phi, T) - - fvm::laplacian(kappaEff, T) + - fvm::laplacian(alphaEff, T) ); TEqn.relax(); diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H index 477a322833..6d7fd8d823 100644 --- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H +++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H @@ -65,12 +65,12 @@ ); // kinematic turbulent thermal thermal conductivity m2/s - Info<< "Reading field kappat\n" << endl; - volScalarField kappat + Info<< "Reading field alphat\n" << endl; + volScalarField alphat ( IOobject ( - "kappat", + "alphat", runTime.timeName(), mesh, IOobject::MUST_READ, diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/setRegionSolidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/setRegionSolidFields.H index a36fa7e98f..37eea84e92 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/setRegionSolidFields.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/solid/setRegionSolidFields.H @@ -12,7 +12,7 @@ if (!thermo.isotropic()) { - tmp tkappaCp = thermo.Kappa()/cp; + tmp tkappaByCp = thermo.Kappa()/cp; const coordinateSystem& coodSys = coordinates[i]; tAnialpha = @@ -29,13 +29,13 @@ IOobject::NO_WRITE ), mesh, - tkappaCp().dimensions(), + tkappaByCp().dimensions(), zeroGradientFvPatchVectorField::typeName ) ); volSymmTensorField& Anialpha = tAnialpha(); - Anialpha.internalField() = coodSys.R().transformVector(tkappaCp()); + Anialpha.internalField() = coodSys.R().transformVector(tkappaByCp()); Anialpha.correctBoundaryConditions(); } diff --git a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/UEqn.H b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/UEqn.H index 430e60c848..eaaa3e60c9 100644 --- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/UEqn.H +++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/UEqn.H @@ -5,7 +5,7 @@ tmp UEqn fvm::ddt(U) + fvm::div(phi, U) + turbulence->divDevReff(U) - == + == fvOptions(U) ); @@ -17,7 +17,7 @@ rAU = 1.0/UEqn().A(); if (pimple.momentumPredictor()) { - solve(UEqn() == -fvc::grad(p) + fvOptions(U)); + solve(UEqn() == -fvc::grad(p)); fvOptions.correct(U); } diff --git a/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C b/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C index 2406cb667c..14bf0ee9eb 100644 --- a/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C +++ b/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -47,6 +47,8 @@ Description #include "snapParameters.H" #include "layerParameters.H" #include "vtkSetWriter.H" +#include "faceSet.H" +#include "motionSmoother.H" using namespace Foam; @@ -675,6 +677,31 @@ int main(int argc, char *argv[]) } + { + // Check final mesh + Info<< "Checking final mesh ..." << endl; + faceSet wrongFaces(mesh, "wrongFaces", mesh.nFaces()/100); + motionSmoother::checkMesh(false, mesh, motionDict, wrongFaces); + const label nErrors = returnReduce + ( + wrongFaces.size(), + sumOp