[Luis F]

Dear all,

I am trying to modify my own solver based in simpleFoam. I managed to make a copy and paste it in my own folder, but that's it.

I am working in the pitzDaily tutorial and my intention is to introduce Reynolds stresses calculated by me and simulate it. Basically I don't need to solve the transport equations for k and epsilon (Im using the k-psilon model), because I already have the Reynolds stresses.

I assume I have to modify the dicrectory createFields.H, but my knowledge in C++ is 0.

I appreciate any hint in that regard.

[RobertHB]

There is a bunch of tutorials on how to modify/create your own solver. For example on youtube [1] [2], the OpenFoam wiki [3], sourceflux [4] and available from Chalmers courses [5]. All of those are a good place to start.

[Luis F]

Dear Robert,

Thank you very much for your answer, I will take a look

[RobertHB]

Not to forget this great collection of slides from C. Soulaine on researchgate.

[Luis F]

Dear all,

I have checked all the links mentioned in this post and I reached the conclusion that the best way to handle it, is to add a new term in the momentum equation:

// Momentum predictor

tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divDevReff(U)
+ fvc::div(RN) // Added by me: influence of Rij
==
fvOptions(U)
);

UEqn().relax();

fvOptions.constrain(UEqn());

solve(UEqn() == -fvc::grad(p));

fvOptions.correct(U);

Being RN an increment of the reynolds stresses calculated by me.

Nevertheless I get an error when I try to build the solver (.dep file).

I created the RN in the createFields.H like this:

Info<< "Reading field RN\n" << endl;
volSymmTensorField RN
(
IOobject
(
"RN",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

Best regards