|
[Sponsors] |
April 24, 2013, 04:38 |
regionCouple problem
|
#1 |
Member
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13 |
Hi all
I am trying to merge conjugate heat foam with a fluid structure interaction solver in order to calculate the temperature field in a fluid and a solid region, however when i run my case i get the following error: Create time Create dynamic mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: laplaceFaceDecomposition --> FOAM FATAL ERROR: Unknown faceTetPolyPatchFaceDecomp type regionCouple. Valid faceTetPolyPatchFaceDecomp types are : 7 ( wall processor wedge patch empty ggi symmetryPlane ) From function faceTetPolyPatchFaceDecomp::New(const polyPatch&, const tetPolyBoundaryMeshFaceDecomp&) in file tetPolyMeshFaceDecomp/tetPolyPatches/faceTetPolyPatch/newFaceTetPolyPatchFaceDecomp.C at line 59. FOAM exiting Anyone have any hints please? Best regards |
|
April 26, 2013, 21:35 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128 |
Greetings Samer,
Without knowing what steps you've taken and which OpenFOAM version/variant you're using, it's very difficult to help you. My best bet would be that you haven't properly analysed/studied each one of the solvers you're basing yourself on and that you are not taking into account the need to:
Bruno
__________________
|
|
April 27, 2013, 05:01 |
|
#3 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
Hi Samer,
The reason that you are running into problems are that you are using the tetFem-mesh motion solver from 1.6-ext. This solver does not have the a boundary type "regionCouple". You should be able begin with the use of other mesh solvers, which solves the motion on the cell centres. Since they rely on the finite volume discretisation and as do your region-coupled model, then these type of mesh motion solvers should work in your case. What I mean by other mesh solvers are those that do not originate from the tetFem framework. Kind regards Niels |
|
April 27, 2013, 18:23 |
|
#4 |
Member
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13 |
Hi Niels
I am sorry that i am new to openfoam, in the fluid/system file there is a tetFemsolution and it sets the mesh motion using the Fe, how could i change this file to set it to a fvsolver? in the movefluidmesh.H file there is: if (feMotionSolver){ ..... } else if (fvMotionSolver) { } so how can i change this to use a fv motion solver?do you have an idea about a fv solver that supports the regionCouple condition? Thank you in advance |
|
April 27, 2013, 18:29 |
|
#5 |
Member
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13 |
hello Bruno
I am using OpenFoam 1.6 ext, and i am implementing that temperature coupling in a fluid structure interaction that simulations a fluid flow in a channel with an elastic Flap downstream of the flow, and i need to write the temperature coupling at the interface fluid-solid, but i don't know what to do since i am new to openFoam. I tried so far to implement the solveEnergy.H in the FSI loop: coupledFvScalarMatrix TEqns(2); // Add fluid equation TEqns.set ( 0, new fvScalarMatrix ( fvm::ddt(T) + fvm::div(phi, T) - fvm::laplacian(DT, T) ) ); // Add solid equation TEqns.set ( 1, new fvScalarMatrix ( fvm::ddt(Tsolid) - fvm::laplacian(DTsolid, Tsolid) ) ); TEqns.solve(); and I specified T initial condition as regionCoupling, and the boundary for the flap as regionCouple, and i got the following error. I hope this helps, anymore requested information i will tell you about it. Best regards |
|
April 28, 2013, 06:23 |
|
#6 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128 |
Greetings to all!
@Samer: There are two types of solvers you can study to figure out what you need:
Best regards, Bruno
__________________
|
|
April 29, 2013, 05:08 |
|
#7 |
Member
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13 |
Greetings to all
I merged the MRconjugateHeatFoam solver with the fluid structure interaction solver, and i set the neighbour region as fluid for the solid boundary conditions and the neighbour region to solid in the fluid boundary conditions. However i get this error whenever i run the simulation: --> FOAM FATAL ERROR: request for fvMesh fluid from objectRegistry fluid failed available objects of type fvMesh are 2 ( region0 solid ) why it doesn't see the fluid neighbour region?? |
|
April 29, 2013, 05:32 |
|
#8 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
Hi Samer,
One of the meshes are (always?) called region0, at least if you partly rely on the standard creation files for one of the meshes. This means that your fluid mesh is called region0 in your simulation. Kind regards, Niels |
|
April 29, 2013, 05:37 |
|
#9 |
Member
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13 |
Hi Niels
The thing is conjugate heat foam case has a different file structure than my case, as there is the directory 0 and inside it there is subdirectories fluid and solid, whereas my fluid structure interaction case, has a fluid and solid directory and inside the fluid directory there's time step subdirectories, could this have an influence on the result or the simulation? |
|
April 29, 2013, 05:49 |
|
#10 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
It is impossible to say without the source code. However, I suspect that the fluid mesh has been created by a call to createFvMesh.H or createDynamicFvMesh.H.
These files will create a fluid mesh, which is called region0, why your file structure must obey to this. / Niels |
|
April 29, 2013, 05:57 |
|
#11 |
Member
Samer
Join Date: Jan 2013
Posts: 31
Rep Power: 13 |
hi Niels
would you like to send the code with the test case?, because the temperature residual are fine, but when it comes to compute the heat flux residuals it always say in the simulation that it is equal to 0 even there is some difference between the heat on the fluid patch and the heat on the solid patch. |
|
April 29, 2013, 06:03 |
|
#12 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,900
Rep Power: 37 |
Hi Samer,
I do unfortunately not have the time to go through your source code. Kind regards and good luck Niels |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
UDF compiling problem | Wouter | Fluent UDF and Scheme Programming | 6 | June 6, 2012 05:43 |
Gambit - meshing over airfoil wrapping (?) problem | JFDC | FLUENT | 1 | July 11, 2011 06:59 |
natural convection problem for a CHT problem | Se-Hee | CFX | 2 | June 10, 2007 07:29 |
Adiabatic and Rotating wall (Convection problem) | ParodDav | CFX | 5 | April 29, 2007 20:13 |
Is this problem well posed? | Thomas P. Abraham | Main CFD Forum | 5 | September 8, 1999 15:52 |