staneuski / dualFuelEngine

OpenFOAM solver based on phenomenological compression model for dual-fuel marine engines
GNU Affero General Public License v3.0
13 stars 3 forks source link

alpha_exh written in the wrong way (or completely not written) #2

Closed staneuski closed 5 years ago

staneuski commented 5 years ago

Не могу реализовать нормальный расчёт третьего поля. Удалось заставить его записывать файл, но он писал в него просто граничные условия, а поле не создавал.

multiCompression:

// Exhaust gas concentration
volScalarField alpha_exh
(
       dimensionedScalar("1", dimless, 1)
     - alpha_gas
     - alpha_air
);

createFields:

volScalarField alpha_exh
(
    IOobject
    (
        "alpha_exh",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedScalar("alpha_exh", dimless, 0)
);

Вот пример записанного при решении файла:

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0.4";
    object      alpha_exh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    injectionXMinus
    {
        type            calculated;
        value           uniform 0;
    }
    injectionXPlus
    {
        type            calculated;
        value           uniform 0;
    }
    injectionYMinus
    {
        type            calculated;
        value           uniform 0;
    }
    injectionYPlus
    {
        type            calculated;
        value           uniform 0;
    }
    inletXMinus
    {
        type            calculated;
        value           uniform 0;
    }
    inletXPlus
    {
        type            calculated;
        value           uniform 0;
    }
    inletYMinus
    {
        type            calculated;
        value           uniform 0;
    }
    inletYPlus
    {
        type            calculated;
        value           uniform 0;
    }
    outlet
    {
        type            calculated;
        value           uniform 0;
    }
    valve
    {
        type            calculated;
        value           uniform 0;
    }
    walls
    {
        type            calculated;
        value           uniform 0;
    }
}

// ************************************************************************* //
staneuski commented 5 years ago

Solved in commit 6cf9b35