CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

sliding mesh for valves & sonicTurbDyMEngineFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By kalle

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 6, 2013, 18:16
Default sliding mesh for valves & sonicTurbDyMEngineFoam
  #1
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15
sasanghomi is on a distinguished road
Hi foamers ,

I am simulating an engine with sonicTurbDyMEngineFoam and simpleEngineTopoFVmesh library.
I have a question about using sliding mesh for valves .

In a 2D case mesh generated with blockMesh and two blocks generated in the valve curtain and they are separated by using different vertices . My problem comes when I want to use gambit for mesh generation (Actually I'm testing gambit for 3D simulation)
In gambit I splite the face in valvecurtain and I created two edges at valvecurtain.But when I imported mesh in openfoam and I ran the case I had an error :

Code:
--> FOAM FATAL ERROR: 
Face 2229 reduced to less than 3 points.  Topological/cutting error B.
Old face: 2(22 1432) new face: 2(22 1432)
old points: 2((-0.01 -0.003 0.00107703) (-0.01 -0.003 -0.00107703))

    From function void slidingInterface::coupleInterface(polyTopoChange& ref) const
    in file polyMeshModifiers/slidingInterface/coupleSlidingInterface.C at line 1765.

FOAM aborting
I think this error comes from sliding mesh!!! I confused....
please anyone help me...any idea??How can I solve this problem??

Thanks and best regards,
Sasan.

Last edited by sasanghomi; March 8, 2013 at 05:00.
sasanghomi is offline   Reply With Quote

Old   March 8, 2013, 05:35
Default
  #2
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15
sasanghomi is on a distinguished road
I appreciate any help from you....
sasanghomi is offline   Reply With Quote

Old   March 9, 2013, 08:11
Default
  #3
Member
 
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 13
novakm is on a distinguished road
Hi sasan

I have reached the same problem as you. The problem is in the mesh generation.
I suppose that you have generated whole geometry at once in gambit. The problem is truly in the sliding interface. In the case of whole geometry generation in some "mesher" (I use Ansa) the stitchMesh (basically static sliding interface) needs the -perfect flag since the points on Cyl and Port are same (within some tolerance). But this is impossible in the case of Curtain in ICE, because of the valve inside. I have evaded this problem by generating all the parts separately (ports, cylinder, intakePorts, exHaustPorts). Then I used mergeMeshes to put them together. In this case the topological errors disappeared. But I must warn you. I have not been able to force the valves to move.
The motion of piston works fine, but when the valve should move the solver just stops with any error message. I have tracked the error and it show up, that is something wrong in the line

phi = fvc::interpolate(rho)*((fvc::interpolate(U) & mesh.Sf()) - fvc::meshPhi(rho, U));

but there my knowledge of OF ends.
Does anyone know, where is the dog buried?
novakm is offline   Reply With Quote

Old   March 9, 2013, 09:28
Default
  #4
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15
sasanghomi is on a distinguished road
Hi Martin ,

Can you explain what does stitchMesh do exactly?
Is it necessary for our simulation???

thanks,
Sasan.
sasanghomi is offline   Reply With Quote

Old   March 9, 2013, 13:13
Default
  #5
Member
 
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 13
novakm is on a distinguished road
Quote:
Originally Posted by sasanghomi View Post
Hi Martin ,

Can you explain what does stitchMesh do exactly?
Is it necessary for our simulation???

thanks,
Sasan.
StitchMesh stitch the mesh. It is not necessary for the simulation, but it is a good pointer on the validity of the mesh. If stitchMesh works correctly, then the mesh is valid for the sliding interface, i.e the sliding interface is possible.
novakm is offline   Reply With Quote

Old   March 9, 2013, 18:42
Default
  #6
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15
sasanghomi is on a distinguished road
Hi Martin

I have an error when I use stitchMesh...

--> FOAM FATAL ERROR:
Face 2229 reduced to less than 3 points. Topological/cutting error B.
Old face: 2(22 1432) new face: 2(22 1432)
old points: 2((-0.01 -0.003 0.00107703) (-0.01 -0.003 -0.00107703))

From function void slidingInterface::coupleInterface(polyTopoChange& ref) const
in file polyMeshModifiers/slidingInterface/coupleSlidingInterface.C at line 1765.

FOAM aborting

Aborted (core dumped)

It means that the slave patch and master patch have same vertices????
please correct me if I am wrong...
In Gambit I splite my geometry with face and
connect selector switch was deactive and I created two edges in the interface (2D case)...where is problem????

thanks,
Sasan.

sasanghomi is offline   Reply With Quote

Old   March 11, 2013, 11:15
Default
  #7
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15
sasanghomi is on a distinguished road
please help me for solving this problem......How can I remove this error???
sasanghomi is offline   Reply With Quote

Old   March 12, 2013, 06:50
Default
  #8
Senior Member
 
Karl-Johan Nogenmyr
Join Date: Mar 2009
Location: Linköping
Posts: 279
Rep Power: 21
kalle is on a distinguished road
Make sure your sliding interface patches are absolutely flat and that they lie in the same plane. The mesher may not do this sufficiently well. Check the cells located at coordinate (-0.01 -0.003 0.00107703). If they are awkward, remesh the whole thing and get rid of ugly cells where you use sliding interface.

K
sasanghomi likes this.
kalle is offline   Reply With Quote

Old   March 28, 2013, 07:53
Default
  #9
Member
 
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 13
novakm is on a distinguished road
Sorry, I have been very busy recently. Does the Kalle's solution help?
novakm is offline   Reply With Quote

Old   March 28, 2013, 11:22
Default
  #10
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15
sasanghomi is on a distinguished road
Hi Martin ,

I created the mesh again and I could overcome that error(thanks for karl's suggestion)...But I have a new problem in the 2D simulation that I described it in the below thread :
http://www.cfd-online.com/Forums/ope...imulation.html

And unfortunately I have another error for 3D simulation :
http://www.cfd-online.com/Forums/ope...ion-error.html

Can you help me? please help me..
I appreciate any help from you.

Thanks and best regards,
Sasan.
sasanghomi is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 07:42
[snappyHexMesh] Layers:problem with curvature giulio.topazio OpenFOAM Meshing & Mesh Conversion 10 August 22, 2012 10:03
[snappyHexMesh] external flow with snappyHexMesh chelvistero OpenFOAM Meshing & Mesh Conversion 11 January 15, 2010 20:43
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 15:09


All times are GMT -4. The time now is 20:07.