|
[Sponsors] |
sliding mesh for valves & sonicTurbDyMEngineFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 6, 2013, 18:16 |
sliding mesh for valves & sonicTurbDyMEngineFoam
|
#1 |
Senior Member
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15 |
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 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. |
|
March 8, 2013, 05:35 |
|
#2 |
Senior Member
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15 |
I appreciate any help from you....
|
|
March 9, 2013, 08:11 |
|
#3 |
Member
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 13 |
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? |
|
March 9, 2013, 09:28 |
|
#4 |
Senior Member
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15 |
Hi Martin ,
Can you explain what does stitchMesh do exactly? Is it necessary for our simulation??? thanks, Sasan. |
|
March 9, 2013, 13:13 |
|
#5 |
Member
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 13 |
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.
|
|
March 9, 2013, 18:42 |
|
#6 |
Senior Member
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15 |
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. |
|
March 11, 2013, 11:15 |
|
#7 |
Senior Member
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15 |
please help me for solving this problem......How can I remove this error???
|
|
March 12, 2013, 06:50 |
|
#8 |
Senior Member
Karl-Johan Nogenmyr
Join Date: Mar 2009
Location: Linköping
Posts: 279
Rep Power: 21 |
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 |
|
March 28, 2013, 07:53 |
|
#9 |
Member
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 13 |
Sorry, I have been very busy recently. Does the Kalle's solution help?
|
|
March 28, 2013, 11:22 |
|
#10 |
Senior Member
Sasan Ghomi
Join Date: Sep 2012
Location: Denmark
Posts: 292
Rep Power: 15 |
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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |