[delMonster]

Hi everyone, I'd much appreciate some advice / guidance. I've used codedMixed quite a lot to vary boundary conditions based on time. I was hoping someone has perhaps varied inlet velocity based on a temperature located at a predefined location in the system volume and could perhaps share how they achieved this. Alternatively anyone knowledgeable that would be so kind as to point me in the right direction. Thanks

[Alczem]

Hey!


I think inputValueMapper can probably achieve what you do. Here is an example by Tobias Holzmann:



https://www.linkedin.com/pulse/new-f...obias-holzmann


I tried to use it several times, but I find it very confusing you basically need a functionObjet that stores some quantities that inputValueMapper can then read and use.


Also, I could not make it work with precompiled OpenFoam for Windows (v2212). I had to recompile OpenFoam with slight modifications (cant remember what, but a thread was discussing this) and run my tests on Linux. Maybe it was fixed though.



If you manage to make it work, please share the setup, I am curious to see more examples of this feature that could be really useful!


EDIT:


This thread shows the changes to the code, but it was because of one specific boundary condition. Other boundary conditions worked even with precompiled OpenFoam.

[delMonster]

Originally Posted by Alczem

Hey!


I think inputValueMapper can probably achieve what you do. Here is an example by Tobias Holzmann:



https://www.linkedin.com/pulse/new-f...obias-holzmann


I tried to use it several times, but I find it very confusing you basically need a functionObjet that stores some quantities that inputValueMapper can then read and use.


Also, I could not make it work with precompiled OpenFoam for Windows (v2212). I had to recompile OpenFoam with slight modifications (cant remember what, but a thread was discussing this) and run my tests on Linux. Maybe it was fixed though.



If you manage to make it work, please share the setup, I am curious to see more examples of this feature that could be really useful!


EDIT:


This thread shows the changes to the code.

Thanks, really appreciate that. I'm on wsl at the moment. I'll definitely share if I can get it to work.

[sharonyue]

This can be easily done by codedFixedValue. A template can be seen as:


inlet
{
type codedFixedValue;
value uniform (10 0 0);
name linearTBC1;
code
#{
const fvMesh& mesh = this->patch().boundaryMesh().mesh();
volVectorField& U = mesh.lookupObjectRef<volVectorField>("U");
#};
}

What you need is just access your field (U in the above code, T for your code) in the boundary condition.

[delMonster]

Thanks, I'll give that a try! Much appreciated.

[delMonster]

Hi sharonyue,


I tried your proposed method below:

Originally Posted by sharonyue

This can be easily done by codedFixedValue. A template can be seen as:


inlet
{
type codedFixedValue;
value uniform (10 0 0);
name linearTBC1;
code
#{
const fvMesh& mesh = this->patch().boundaryMesh().mesh();
volVectorField& U = mesh.lookupObjectRef<volVectorField>("U");
#};
}

What you need is just access your field (U in the above code, T for your code) in the boundary condition.

Since I want to access the temperature field at the outlet boundary condition I modified it as follows:

const fvMesh& mesh = this->patch().boundaryMesh().mesh();
label bottomoutlet = mesh.boundaryMesh().findPatchID("bottom_outlet");
volScalarField& T = mesh.lookupObjectRef<volScalarField>("T");
const fvPatchScalarField& T_p = T.boundaryField()[bottomoutlet];

Info << "T at outlet = " << T_p << " K" << endl;

However, I'm not getting back a value, I'm getting the initial boundary condition type:

T at outlet = type inletOutlet;
inletValue nonuniform List<scalar> 0();
value nonuniform List<scalar> 0(); K

I'm probably missing something silly. Have you got any idea where I'm going wrong?

[delMonster]

Originally Posted by Alczem

If you manage to make it work, please share the setup, I am curious to see more examples of this feature that could be really useful!


EDIT:


This thread shows the changes to the code, but it was because of one specific boundary condition. Other boundary conditions worked even with precompiled OpenFoam.

Looking at that thread, I was not able to get this working for that boundary condition