[Djordje]

Dear foamers,

I have some troubles with comment parts of this code:

PHP Code:

label myinlet = mesh_.boundaryMesh().findPatchID("inlet");
const polyPatch& cPatch = mesh_.boundaryMesh()[myinlet];
const fvPatchVectorField& Uin = U().boundaryField()[myinlet];
//const fvPatchScalarField& UPhiIN = UPhi_()->boundaryField()[myinlet];
const surfaceScalarField& magSf = mesh_.magSf();

scalar IFR = 0.0;
scalar IEFzDir = 0.0;
//scalar IEFPhiDir = 0.0; 
    
forAll(cPatch, facei) 
{ 
    IFR += (Uin[facei].z())*magSf.boundaryField()[myinlet][facei];
    IEFzDir += pow(Uin[facei].z(),3)*magSf.boundaryField()[myinlet][facei];
    //IEFPhiDir += sqr(UPhiIN[facei])*(Uin[facei].z())*magSf.boundaryField()[myinlet][facei];
}   

//SwS_=IEFPhiDir/(0.001+IEFzDir); 

Info<<endl<<l_<<endl;
Info<<Ctau_<<endl;
Info<<"InFlowRate = "<<IFR<<" m^3/s"<<endl;
Info<<"InEnFluxZdir = "<<IEFzDir<<" m^5/s^3"<<endl;
Info<<"InEnFluxPhiDir = "<<IEFPhiDir<<" m^5/s^3"<<endl;
//Info<<"SwS = "<<SwS_<<endl<<endl<<endl; 


I want to calculate the swirl strenght (SwS) on inlet section. The code without comments works fine up to 200 iterations approximately, and then it breaks with next message:

PHP Code:

l [0 0 0 0 0 0 0] 0.09
Ctau [0 0 0 0 0 0 0] 0.182
InFlowRate = 0.647061086243 m^3/s
InEnFluxZdir = 307.942132036 m^5/s^3
InEnFluxPhiDir = 37.1113979524 m^5/s^3
SwS = 0.12051380301


ExecutionTime = 96.61 s  ClockTime = 97 s

Time = 222

smoothSolver:  Solving for Ux, Initial residual = 0.00173071772014, Final residual = 0.000106084441356, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.00173017679547, Final residual = 0.000105942528607, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.000569924142604, Final residual = 3.36294823708e-05, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.0077950722647, Final residual = 7.5300467066e-05, No Iterations 158
time step continuity errors : sum local = 0.000354384134085, global = 7.11056780688e-05, cumulative = 5.24229577773e-05
smoothSolver:  Solving for epsilon, Initial residual = 0.00130095491961, Final residual = 8.68883144141e-05, No Iterations 4
smoothSolver:  Solving for k, Initial residual = 0.00253765686116, Final residual = 0.000191029833772, No Iterations 4
[HP-Z6-G4:14979] *** Process received signal ***
[HP-Z6-G4:14979] Signal: Segmentation fault (11)
[HP-Z6-G4:14979] Signal code:  (-6)
[HP-Z6-G4:14979] Failing at address: 0x3e900003a83
[HP-Z6-G4:14979] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f895e59af20]
[HP-Z6-G4:14979] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f895e59ae97]
[HP-Z6-G4:14979] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3ef20)[0x7f895e59af20]
[HP-Z6-G4:14979] [ 3] /home/djnovkovic/foam/foam-extend-4.0/lib/linux64GccDPOpt/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels10kEpsHamNov7correctEv+0x2e62)[0x7f8962696d22]
[HP-Z6-G4:14979] [ 4] simpleFoam(+0x1caf7)[0x561b420b9af7]
[HP-Z6-G4:14979] [ 5] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f895e57db97]
[HP-Z6-G4:14979] [ 6] simpleFoam(+0x1e1ba)[0x561b420bb1ba]
[HP-Z6-G4:14979] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node HP-Z6-G4 exited on signal 11 (Segmentation fault).
-------------------------------------------------------------------------- 


I think that I have problem with method of acces Uphi_() values on the inlet section. I have tried to use:
const fvPatchScalarField& UPhiIN = UPhi_().boundaryField()[myinlet];
but then I can't compile code and I obtain next error:

PHP Code:

error: ‘class Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >’ has no member named ‘boundaryField’
 const fvPatchScalarField& UPhiIN = UPhi_().boundaryField()[myinlet]; 


What is wrong in this part of my code?

[Djordje]

I find solution. The problem was in Uphi_() velocitiy definition. Previously I defined it with next set of commands:

...

 X_                                               
    (
          mesh_.C().component(vector::X)
    ),

    Y_
    (
          mesh_.C().component(vector::Y)
    ),

    Z_
    (
          mesh_.C().component(vector::Z)
    ),

    r_
    (

          sqrt(sqr(X_)+sqr(Y_))+rSmall_
    ),

    SinPhi_
    (
          Y_/r_
    ),

    CosPhi_
    (
          X_/r_
    ),

...

tmp<volScalarField> kEpsHamNov::Ux_() const
{
    return tmp<volScalarField>
    (
        new volScalarField
        (
            IOobject
            (
                "Ux_",
                runTime_.timeName(),
                U_.db(),
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
           U_.component(vector::X)
        )
    );
}

tmp<volScalarField> kEpsHamNov::Uy_() const
{
    return tmp<volScalarField>
    (
        new volScalarField
        (
            IOobject
            (
                "Uy_",
                runTime_.timeName(),
                U_.db(),
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
           U_.component(vector::Y)
        )
    );

}

tmp<volScalarField> kEpsHamNov::UPhi_() const
{
    return tmp<volScalarField>
    (
        new volScalarField
        (
            IOobject
            (
                "UPhi_",
                runTime_.timeName(),
                U_.db(),
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
           -Ux_()*SinPhi_+Uy_()*CosPhi_
        )
    );

}

...

Later I made modification and I calculate Uphi_ velocity in problematic part of code on this manner:

label myinlet = mesh_.boundaryMesh().findPatchID("inlet");
const polyPatch& cPatch = mesh_.boundaryMesh()[myinlet];
const fvPatchVectorField& Uin = U().boundaryField()[myinlet];
const fvPatchScalarField& SinIN = SinPhi_.boundaryField()[myinlet];
const fvPatchScalarField& CosIN = CosPhi_.boundaryField()[myinlet];
const surfaceScalarField& magSf = mesh_.magSf();

scalar IFR = 0.0;
scalar IEFzDir = 0.0;
scalar IEFPhiDir = 0.0;

forAll(cPatch, facei)
{
    IFR += (Uin[facei].z())*magSf.boundaryField()[myinlet][facei];
    IEFzDir += pow(Uin[facei].z(),3)*magSf.boundaryField()[myinlet][facei];
    IEFPhiDir += sqr(-Uin[facei].x()*SinIN[facei]+Uin[facei].y()*CosIN[facei])*
                 (Uin[facei].z())*magSf.boundaryField()[myinlet][facei];
}

SwS_=IEFPhiDir/(0.001+IEFzDir);

Info<<endl<<l_<<endl;
Info<<Ctau_<<endl;
Info<<"InFlowRate = "<<IFR<<" m^3/s"<<endl;
Info<<"InEnFluxZdir = "<<IEFzDir<<" m^5/s^3"<<endl;
Info<<"InEnFluxPhiDir = "<<IEFPhiDir<<" m^5/s^3"<<endl;
Info<<"SwS = "<<SwS_<<endl<<endl<<endl;

It converge fine now.