CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   Outlet fields as inlet fields, how to. (

Villo January 22, 2014 06:43

Outlet fields as inlet fields, how to.
3 Attachment(s)
i`m running a unsteady RANS simulation (k-epsilon) with twoLiquidMixingFoam.
In order to reduce the computational domain showed in the top of the attached figure1 i need to use the outlet from the bottom left geometry (at -latestTime) as inlet condition for the bottom right geometry.
I tried two ways (to take fields from the outlet of the bottom left geometry):
1) mapFields (with mapFieldsDict), the result of the mapped field (for U) is show in figure2
2) sample (with sampleDict), the result of the mapped field (for U) is show in figure3
In both cases k, epsilon, U, nut are well mapped: I have set these fields as inlet condition for the bottom left geometry in 0/U 0/k etc etc.
BUT seems to be the solver does not "feel" this nonuniform List of inlet condition and solved fields are wherever 0.
p.s. If to answer you need some dictionaries just ask!

Villo January 22, 2014 11:00

Making it more interesting...
3 Attachment(s)
Just to make it more interesting (:rolleyes: or more (un)clear...?) here you some more information:
- attached the mapFieldsDict.txt and sampleDict.txt used to map fields at outlet patch
- attached the example of 0/U obtained after run mapFields and sample utility (for the last one I have just copy the sample result after value nonuniform List<vector>)

Villo January 22, 2014 11:22

Making it more interesting... part2
Browsing for similar problems i prefer show you the constant/polymesh/boundary. The mapping/sampling target is the patch int_1B (and the mapping/sampling source is the patch int_1A)
Could be that my problem is just due from the definition of boundary condition on int_1B? So (eg) chancing it from patch to mappedPatch (and similarly for all the 0/fields)

type wall;
nFaces 84456;
startFace 3005866;
type patch;
nFaces 1637;
startFace 3090322;
type patch;
nFaces 1637;
startFace 3091959;

romant February 4, 2014 11:06

Did you solve the problem? And if yes, could you share your procedure step by step?

alexeym February 4, 2014 11:44

Did you try something like this:


        type            outletMappedUniformInlet;
        outletPatchName outlet;
        phi            phi;
        value          uniform 0;

        type            zeroGradient;

Or this is not what you want? As from the pictures I wasn't able to see any difference in geometries.

Sorry forgot that you need to rotate vectors. My solution is wrong in this case.

romant February 4, 2014 15:23

3 Attachment(s)

the way I solved something similar applied to your problem:

1. Sample the outlet patch int_1A -> with the result that you get a folder postProcessing/surfaces/latestTime/sampleValues/vectorField with contains the velocity values and then under sampleValues, which is just the name of the surface set you have faceCentres, faces and points
2. Next I renamed the faceCentres to points and discarded faces and points completely.
3. Put the new points file into a folder which has the same name as the boundary patch that you would like to use (in your case int_1B) and also create a 0 folder in there.
4. In the 0 folder you put all the sampled fields that you would like to have an inlet condition for (U, nu, p, ...)
5. The files are missing headers. I attached the necessary headers in this post, the points one gets the points.txt header, the vector fields get vector headers and the scalar field get the scalar header.
5. The full folder int_1B , containing the 0 folder and the points file are now copied to constant/boundaryData
6. The inlet boundary condition gets the following type

        type            timeVaryingMappedFixedValue;
        setAverage      true;
        offset          (0 0 0);

where the offset could be used to move it to the right location, but as your sampling and the mapping are most likely in the same location, you don't have to offset it. The setAverage is set to true in the above example, however, one must then set also an average value in the field files under boundaryData, otherwise it will set it to 0, as this is the average that I have in the header files at the moment. Set to false if no averaging should be done.
7. run the case

Villo February 5, 2014 10:21

1 Attachment(s)
Hi Alexey and Roman, thanks for your interest.

The problem (solved for RANS computation)was due to my mistake because i was using the outlet from the left geometry as inlet for the right geometry... independently from the followed procedure (mapFields, sample... or just copying the field from the time folder of interest as i did at the end) the direction of the flow is maintained... OpenFoam was feeling my flow as going out and not coming in :D
So, what i did can be summarized as:
1) On the the whole geometry create an AMI in the location where i want to split the geometry (let`s say left_AMI and right_AMI)
2) From the solved whole geometry copy the fields (from left_AMI patch) that have to be mapped
3) On the new case reduced geometry (in my case the one on the right) paste this fields in the correspondent 0/fields for the right_AMI patch (the new inlet)

This will keep the flow direction.

@Alexey i will try on one other case your suggestion "outletMappedUniformInlet"... could works nice, bypassing the need to simulate the whole geometry as i did. I will let you know if it works.
@Roman "timeVaryingMappedFixedValue" is nice but i prefer to use it just for transient inlet. By the way you show me how to use it ;)

BUT... :rolleyes:
I said "I solve it for RANS computation."... because trying to apply the procedure laid before and more as sample and mapFields, on LES computation... well, look at the new attached picture.
I thought it was due to the enumeration of cells, but if it was so I should have had the same problem in the RANS computation.
Guys, any idea?

__________________________________________________ ______________________
@Alexey I read the edit to your post only after answering. The origin of my problem was exactly the no rotation of the vector field (U). So let me understand, i cannot apply "outletMappedUniformInlet" because does not reverse the vector field? So what is the purpose of "outletMappedUniformInlet"?

romant February 5, 2014 11:46

Due to the nature of how the timeVaryingMappedFixedValue works, you can even create the sample surface in matlab, meaning you could also of course turn around the velocity vectors in matlab. When you sample in raw instead of foamFile, you get something that can easily be read into matlab, turn around the vectors and put out the 2 files (U and points) in a OpenFoam compatible version.

All times are GMT -4. The time now is 16:45.