CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   mapFields doesn't change fields (http://www.cfd-online.com/Forums/openfoam/113456-mapfields-doesnt-change-fields.html)

Werne February 20, 2013 04:42

mapFields doesn't change fields
 
Hello,
I have a problem with the command mapFields. I'm trying to map the data of a big mesh onto a smaller mesh that is completely surrounded by the bigger mesh. mapFields runs smoothely and doesn't produce any errors, but it also doesn't produce any results as the fields in the target folder aren't changed. I get the following output:
Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.7.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.com                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 1.7.x-d08d3c27ccda
Exec  : mapFields ../source_case/
Date  : Feb 20 2013
Time  : 09:29:09
Host  : ocean03
PID    : 10092
Case  : /home/werne/OpenFOAM1.7.x/werne-1.7.x/run/target_case
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Source: "../source_case" ""
Target: "/home/werne/OpenFOAM1.7.x/werne-1.7.x/run/Project" "target case"
Create databases as time
Source time: 0
Target time: 0
Create meshes
Source mesh size: 6557642      Target mesh size: 230809
 
Mapping fields for time 0
    interpolating hTotal
    interpolating h0
    interpolating manning
    interpolating h
    interpolating U
End

I have tried different approaches:
- I left both patchMap and cuttingPatches empty
- I left patchMap empty and inserted all boundaries of the smaller mesh in cuttingPatches (and vice versa)

I get always the same output and no results. I'm mapping from the time 0 to time 0 and have "dummy files" of the fields I want to map in the target folders.

Any suggestions what I could do, or how I should build up the mapFieldsDict?

nimasam February 20, 2013 05:02

mapFields has some extra options for example -consistent or -time,
you can use these options to map from different time step :) or for different conditions
if you want to see mapFields options you can easily type mapFields xxxx
then all options will be shown in terminal

Werne February 20, 2013 05:12

I already tried mapping from different time folders, but that doesn't help either.
And the consistent option shouldn't work for my case, as my meshs are not consistent, if I understand it correctly.
The bigger mesh has no common boundary with the smaller mesh. It exceeds the size of the smaller mesh in all directions.

edit:
After I've changed the geometry because of an error, I have a new problem: mapFields gets stuck at this point:
Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.7.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.com                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 1.7.x-d08d3c27ccda
Exec  : mapFields ../source_case/
Date  : Feb 20 2013
Time  : 17:35:11
Host  : ocean03
PID    : 2237
Case  : /home/werne/OpenFOAM1.7.x/werne-1.7.x/run/Project/target_case                                                                                                                                                           
nProcs : 1                                                                                                                                                                                                               
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).                                                                                                                                                       
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //                                                                                                                                           
Source: "../source_case" ""                                                                                                                                                                                         
Target: "/home/werne/OpenFOAM1.7.x/werne-1.7.x/run/Project" "target_case"                                                                                                                                                       
 
Create databases as time                                                                                                                                                                                                 
 
Source time: 0                                                                                                                                                                                                           
Target time: 0                                                                                                                                                                                                           
Create meshes                                                                                                                                                                                                             
 
Source mesh size: 3585480      Target mesh size: 230809

It is stuck at this point for about 2 hours already, but it still uses RAM and CPU. As before I used an empty mapFieldsDict and several modifications. Neither gets further than this point.

anishtain4 February 20, 2013 16:57

Have you had used any utility like Mach? when you use the mapFields it tries to map all the fields in the source to the fields in the target, this means you should have similar fields and no additional in anyone.
Though I had the same problem but -consistent solved my case

Werne February 21, 2013 04:40

No, I've not used any other utility. And I have the same fields in the source and in the target case.
And -consistent won't work for me, as my cases are not consisten. It always fails as I have only 3 boundries in my source case but 9 in my target case.

anishtain4 February 23, 2013 13:39

Your case is strange, I can't think of a specific reason but only guessing and throwing some solutions.
what are your initial conditions for the target case? try testing them as uniform


All times are GMT -4. The time now is 18:46.