diff --git a/applications/solvers/combustion/dieselEngineFoam/Make/files b/applications/solvers/combustion/dieselEngineFoam/Make/files deleted file mode 100644 index f643bc9144..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/Make/files +++ /dev/null @@ -1,3 +0,0 @@ -dieselEngineFoam.C - -EXE = $(FOAM_APPBIN)/dieselEngineFoam diff --git a/applications/solvers/combustion/dieselEngineFoam/UEqn.H b/applications/solvers/combustion/dieselEngineFoam/UEqn.H deleted file mode 100644 index c6caf5989c..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/UEqn.H +++ /dev/null @@ -1,14 +0,0 @@ - fvVectorMatrix UEqn - ( - fvm::ddt(rho, U) - + fvm::div(phi, U) - + turbulence->divDevRhoReff(U) - == - rho*g - + dieselSpray.momentumSource() - ); - - if (pimple.momentumPredictor()) - { - solve(UEqn == -fvc::grad(p)); - } diff --git a/applications/solvers/combustion/dieselEngineFoam/YEqn.H b/applications/solvers/combustion/dieselEngineFoam/YEqn.H deleted file mode 100644 index 46d0f364e6..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/YEqn.H +++ /dev/null @@ -1,48 +0,0 @@ -tmp > mvConvection -( - fv::convectionScheme::New - ( - mesh, - fields, - phi, - mesh.divScheme("div(phi,Yi_h)") - ) -); - -{ - combustion->correct(); - dQ = combustion->dQ(); - label inertIndex = -1; - volScalarField Yt(0.0*Y[0]); - - forAll(Y, i) - { - if (Y[i].name() != inertSpecie) - { - volScalarField& Yi = Y[i]; - - fvScalarMatrix YiEqn - ( - fvm::ddt(rho, Yi) - + mvConvection->fvmDiv(phi, Yi) - - fvm::laplacian(turbulence->muEff(), Yi) - == - dieselSpray.evaporationSource(i) - + combustion->R(Yi) - ); - - YiEqn.relax(); - YiEqn.solve(mesh.solver("Yi")); - - Yi.max(0.0); - Yt += Yi; - } - else - { - inertIndex = i; - } - } - - Y[inertIndex] = scalar(1) - Yt; - Y[inertIndex].max(0.0); -} diff --git a/applications/solvers/combustion/dieselEngineFoam/createFields.H b/applications/solvers/combustion/dieselEngineFoam/createFields.H deleted file mode 100644 index bbb9000943..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/createFields.H +++ /dev/null @@ -1,104 +0,0 @@ -Info<< "Creating combustion model\n" << endl; - -autoPtr combustion -( - combustionModels::psiChemistryCombustionModel::New - ( - mesh - ) -); - -psiChemistryModel& chemistry = combustion->pChemistry(); - -hsCombustionThermo& thermo = chemistry.thermo(); - -basicMultiComponentMixture& composition = thermo.composition(); -PtrList& Y = composition.Y(); - -word inertSpecie(thermo.lookup("inertSpecie")); - -if (!composition.contains(inertSpecie)) -{ - FatalErrorIn(args.executable()) - << "Specified inert specie '" << inertSpecie << "' not found in " - << "species list. Available species:" << composition.species() - << exit(FatalError); -} - -volScalarField rho -( - IOobject - ( - "rho", - runTime.timeName(), - mesh - ), - thermo.rho() -); - -Info<< "Reading field U\n" << endl; -volVectorField U -( - IOobject - ( - "U", - runTime.timeName(), - mesh, - IOobject::MUST_READ, - IOobject::AUTO_WRITE - ), - mesh -); - - -volScalarField& p = thermo.p(); -const volScalarField& psi = thermo.psi(); -const volScalarField& T = thermo.T(); -volScalarField& hs = thermo.hs(); - - -#include "compressibleCreatePhi.H" - -Info << "Creating turbulence model.\n" << nl; -autoPtr turbulence -( - compressible::turbulenceModel::New - ( - rho, - U, - phi, - thermo - ) -); - -// Set the turbulence into the combustion model -combustion->setTurbulence(turbulence()); - -Info<< "Creating field DpDt\n" << endl; -volScalarField DpDt -( - fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p) -); - - -multivariateSurfaceInterpolationScheme::fieldTable fields; - -forAll(Y, i) -{ - fields.add(Y[i]); -} -fields.add(hs); - -volScalarField dQ -( - IOobject - ( - "dQ", - runTime.timeName(), - mesh, - IOobject::NO_READ, - IOobject::AUTO_WRITE - ), - mesh, - dimensionedScalar("dQ", dimEnergy/dimTime, 0.0) -); diff --git a/applications/solvers/combustion/dieselEngineFoam/createSpray.H b/applications/solvers/combustion/dieselEngineFoam/createSpray.H deleted file mode 100644 index 8c6706e0ba..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/createSpray.H +++ /dev/null @@ -1,37 +0,0 @@ -Info << "Constructing Spray" << endl; - -PtrList gasProperties(Y.size()); -forAll(gasProperties, i) -{ - gasProperties.set - ( - i, - new gasThermoPhysics - ( - dynamic_cast&> - (thermo).speciesData()[i] - ) - ); -} - -spray dieselSpray -( - U, - rho, - p, - T, - composition, - gasProperties, - thermo, - g -); - -scalar gasMass0 = fvc::domainIntegrate(rho).value(); - -if (dieselSpray.twoD()) -{ - gasMass0 *= constant::mathematical::twoPi/dieselSpray.angleOfWedge(); -} - -gasMass0 -= - dieselSpray.injectedMass(runTime.value()) - dieselSpray.liquidMass(); diff --git a/applications/solvers/combustion/dieselEngineFoam/hsEqn.H b/applications/solvers/combustion/dieselEngineFoam/hsEqn.H deleted file mode 100644 index 7c0b7f9dc3..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/hsEqn.H +++ /dev/null @@ -1,21 +0,0 @@ -{ - fvScalarMatrix hsEqn - ( - fvm::ddt(rho, hs) - + mvConvection->fvmDiv(phi, hs) - - fvm::laplacian(turbulence->alphaEff(), hs) - == - DpDt - + combustion->Sh() - + dieselSpray.heatTransferSource()() - - ); - - hsEqn.relax(); - hsEqn.solve(); - - thermo.correct(); - - Info<< "min/max(T) = " - << min(T).value() << ", " << max(T).value() << endl; -} diff --git a/applications/solvers/combustion/dieselEngineFoam/pEqn.H b/applications/solvers/combustion/dieselEngineFoam/pEqn.H deleted file mode 100644 index 0493a1ff6c..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/pEqn.H +++ /dev/null @@ -1,65 +0,0 @@ -rho = thermo.rho(); - -volScalarField A(UEqn.A()); -U = UEqn.H()/A; - -if (pimple.transonic()) -{ - surfaceScalarField phid - ( - "phid", - fvc::interpolate(psi) - *((fvc::interpolate(U) & mesh.Sf()) - fvc::meshPhi(rho, U)) - ); - - while (pimple.correctNonOrthogonal()) - { - fvScalarMatrix pEqn - ( - fvm::ddt(psi, p) - + fvm::div(phid, p) - - fvm::laplacian(rho/A, p) - == - Sevap - ); - - pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter()))); - - if (pimple.finalNonOrthogonalIter()) - { - phi == pEqn.flux(); - } - } -} -else -{ - phi = fvc::interpolate(rho) - *((fvc::interpolate(U) & mesh.Sf()) - fvc::meshPhi(rho, U)); - - while (pimple.correctNonOrthogonal()) - { - fvScalarMatrix pEqn - ( - fvm::ddt(psi, p) - + fvc::div(phi) - - fvm::laplacian(rho/A, p) - == - Sevap - ); - - pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter()))); - - if (pimple.finalNonOrthogonalIter()) - { - phi += pEqn.flux(); - } - } -} - -#include "rhoEqn.H" -#include "compressibleContinuityErrs.H" - -U -= fvc::grad(p)/A; -U.correctBoundaryConditions(); - -DpDt = fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p); diff --git a/applications/solvers/combustion/dieselEngineFoam/rhoEqn.H b/applications/solvers/combustion/dieselEngineFoam/rhoEqn.H deleted file mode 100644 index ff594d40b4..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/rhoEqn.H +++ /dev/null @@ -1,64 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | - \\ / A nd | Copyright (C) 2011 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 . - -Global - rhoEqn - -Description - Solve the continuity for density. - -\*---------------------------------------------------------------------------*/ - -volScalarField Sevap -( - IOobject - ( - "Sevap", - runTime.timeName(), - mesh, - IOobject::NO_READ, - IOobject::NO_WRITE - ), - mesh, - dimensionedScalar("zero", dimensionSet(1, -3, -1, 0, 0), 0.0) -); - -forAll(Y, i) -{ - if (dieselSpray.isLiquidFuel()[i]) - { - Sevap += dieselSpray.evaporationSource(i); - } -} - -{ - solve - ( - fvm::ddt(rho) - + fvc::div(phi) - == - Sevap - ); -} - -// ************************************************************************* // diff --git a/applications/solvers/combustion/dieselEngineFoam/spraySummary.H b/applications/solvers/combustion/dieselEngineFoam/spraySummary.H deleted file mode 100644 index 5c1e43dc99..0000000000 --- a/applications/solvers/combustion/dieselEngineFoam/spraySummary.H +++ /dev/null @@ -1,30 +0,0 @@ - label Nparcels = dieselSpray.size(); - reduce(Nparcels, sumOp