[anijdon]

hello Martin;
excuse me; I did some change in the solver's code and compiled it eventually, but it seems has a main problem when running a case; It doesn't do well and gives this error:

nano.zip

I don't know what to do; would you guide me?
thanks;

[MartinB]

Hi Maryam ,

just add this line to constant/transportProperties:
nu nu [0 2 -1 0 0 0 0] 1e-5; // value for kinematic viscosity

Change the value to the correct one for your main fluid.

Martin

[anijdon]

hello;
thank you for reply;I want it read the nu from that formula which determined in the createFields not from transportProperties; If I change it's name there would be no problem?
thanks

[MartinB]

Hi,

it will not work the way you have implemented it right now.

In your createFields.H this transport model and a turbulence model (although it might be a laminar one) is defined:

    singlePhaseTransportModel laminarTransport(U, phi);
    autoPtr<incompressible::RASModel> turbulence
    (
        incompressible::RASModel::New(U, phi, laminarTransport)
    );

During the setup of these models nu is read from transportProperties. Later on in your UEqn.H the viscosity from these models is used via this line:

      + turbulence->divDevReff(U)

So your own viscosity is not used so far. You must replace the divDevReff by your own formulation, if you want to continue with your hardwired approach. You can see what has to be done by looking into OpenFOAM-xxx/src/turbulenceModels/incompressible/turbulenceModel/laminar/laminar.C.

The more elegant solution would be to define your own viscosity model. These are located in OpenFOAM-xxx/src/transportModels/incompressible/viscosityModels/. There are examples how to make a new one in the forum.

Martin

[anijdon]

hello Martin;
I am extremely thankful for your kindness but I am beginner in openFoam and not familiar to the way that you said;could you explain more and tell me what shoud I do step by step; I realy apologize to you for my asking;
thanks;

[MartinB]

Hi,

in your UEqn.H change:

    tmp<fvVectorMatrix> UEqn
    (
        fvm::div(phi, U)
      // + turbulence->divDevReff(U)               // removed
      - fvm::laplacian(nuEff, U)                         // nuEff is your effective viscosity
      - fvc::div(nuEff*dev(fvc::grad(U)().T())) // nuEff is your effective viscosity
    );

In your createFields.H remove:

    // singlePhaseTransportModel laminarTransport(U, phi);

    // autoPtr<incompressible::RASModel> turbulence
    //(
    //    incompressible::RASModel::New(U, phi, laminarTransport)
    //);

In your nanoalaki.C remove:

        // turbulence->correct();

And check the line:

    volTensorField gradU = fvc::grad(U); // gradU is not used in this context

If compilation errors occur, read them very carefully and fix them, beginning from the first one.

Martin

[anijdon]

Hello;
thank you a lot; I hope it'll resolve in this way; This problem is just for nu and not for other properties like rho or k ; Am I right?