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

engine simulation with mesh motion and topological changes

Register Blogs Community New Posts Updated Threads Search

Like Tree6Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 18, 2011, 07:00
Default
  #81
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
ok I must be missing something, I keep getting the same error, do I need the dynamicFvMeshDict? it seems like perhaps I need to remove '# include "dynamicFvMesh.H" ' from my solver, but then that brings other errors, like, createDynamicFvMesh.H does not work any more :S
abminternet is offline   Reply With Quote

Old   April 18, 2011, 07:50
Default
  #82
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
wait, I modified the solver, I think I answered my own question
abminternet is offline   Reply With Quote

Old   April 18, 2011, 07:59
Default
  #83
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
hey Peter, so I modified the solver to be able to use engineTopoChangerMesh, and seems to compile alright, do you maybe have a test case I could run? there are many keywords I need to add and I'm not familiar with :S Would appreciate it
abminternet is offline   Reply With Quote

Old   April 18, 2011, 08:13
Default
  #84
Senior Member
 
Join Date: Oct 2009
Posts: 140
Rep Power: 17
Peter_600 is on a distinguished road
Try to replace # include "createDynamicFvMesh.H" with
# include "createEngineDynamicMesh.H"
in the sonicTurbDyMEngineFoam. C file and recompile.
Peter_600 is offline   Reply With Quote

Old   April 18, 2011, 08:31
Default
  #85
Senior Member
 
Join Date: Oct 2009
Posts: 140
Rep Power: 17
Peter_600 is on a distinguished road
I don t have a test case. But feel free to post any question you have. We may debug it together here.
Peter_600 is offline   Reply With Quote

Old   April 18, 2011, 08:54
Default
  #86
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
great I have one well I am testing verticalValves now but perhaps you know, it asks for valveHeadPatch, but also for poppetPatch, I thought they were the same, so I know what is the poppetPatch from my simpleEngine test case, but do u know what is the valveHeadPatch?
abminternet is offline   Reply With Quote

Old   April 18, 2011, 09:38
Default
  #87
Senior Member
 
Join Date: Oct 2009
Posts: 140
Rep Power: 17
Peter_600 is on a distinguished road
http://powerlab.fsb.hr/ped/kturbo/Op...Zagreb2007.pdf
Maybe this can be useful for you.
Peter_600 is offline   Reply With Quote

Old   April 18, 2011, 09:54
Default
  #88
Senior Member
 
Join Date: Oct 2009
Posts: 140
Rep Power: 17
Peter_600 is on a distinguished road
Try for valveHeadPatch -> none;
Peter_600 is offline   Reply With Quote

Old   April 18, 2011, 11:34
Default
  #89
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
thx for the tutorial so one quick question, what is the difference between valveTopTol and deformationLift? they are from the class verticalValves and engineVerticalValves respectively, but from the definition, it looks like they both do the same thing, disable top layering for the valves when lift < valveTopTol (or deformationLift).
abminternet is offline   Reply With Quote

Old   May 20, 2011, 15:30
Default
  #90
New Member
 
Dalibor
Join Date: Aug 2010
Posts: 6
Rep Power: 16
jajcevic is on a distinguished road
Hi guys,
I have started the tutorial "simpleEngineStem" posted by abm, but after only several time steps following error occurs:

"Floating point exception"

I am using 1.6-ext and "sonicTurbDyMEngineFoam".


Do you have any idea why I can not successfully run this tutorial?
Attached Files
File Type: gz simpleEngineStem.tar.gz (10.5 KB, 70 views)
Berstuleski and enginejob like this.
jajcevic is offline   Reply With Quote

Old   May 23, 2011, 11:55
Default
  #91
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
hi Peter,

could you send me the engineGeometry file you used for your accordionEngineMesh class and for the verticalValves too? I would really appreciate it
abminternet is offline   Reply With Quote

Old   May 23, 2011, 15:57
Default
  #92
Senior Member
 
Join Date: Oct 2009
Posts: 140
Rep Power: 17
Peter_600 is on a distinguished road
Hi,

no problem. Here is my engineGeometry file for the accordionEngineMesh class. U will find a lot of things which are not used by the class, but the class requires an entry for it in the script. I didn t check which entry is really needed, thus I used it always like that. I will highlight the parts which are important for you.

I don t have a file for the verticalValves class.

Peter
Quote:
conRodLength conRodLength [0 1 0 0 0 0 0] 0.161;

bore bore [0 1 0 0 0 0 0] 0.084;

stroke stroke [0 1 0 0 0 0 0] 0.09;

clearance clearance [0 1 0 0 0 0 0] 0.0014;

rpm rpm [0 0 -1 0 0 0 0] 2000;



cylinderHeadName cylinderHead;
linerName liner;
pistonAuxPoints pistonPoints;


moveDetach on;
correctPointsMotion yes;

delta 0.0001;
offSet 0.0005;

valveTopTol 1e-3;
bottomValveTol 5e-3;
poppetValveTol 5e-3;
checkMesh on;
cylinderSetName cylinder;


engineTopoChangerMesh accordionEngineMesh;

deformAngle 34;

piston
{
patch piston;
// patch none;

coordinateSystem
{
type cartesian;
origin (0 0 0);
axis (0 0 1);
direction (0 1 0);
}

minLayer 0.0005;
maxLayer 0.0015;
pistonCellSetName emptyfaces;
pistonFaceSetName emptyfaces;
pistonPointSetName emptyPoints;

bowlInPistonPatchName off;
bowlInCylinderPatchName off;

}

accordionEngineMesh
(
valve1
{
coordinateSystem
{
type cylindrical;
origin (-0.0203 -0.015518978 0.097944030);
axis (0 -0.258819 0.965926);
direction (1 0 0);

}
bottomPatch valve1_intake_bottom;
poppetPatch valve1_intake_top;
stemPatch valve1_intake_stem;

sidePatch off;
minLift 0.00025;
curtainInPortPatch off;
curtainInCylinderPatch off;
detachInCylinderPatch valveDetachInCyl1;
detachInPortPatch valveDetachInPort1;

valveHeadPatch off;
detachFaces detach1;
detachDistance 0.95e-8;
detachTol 0.005e-8;
stemEdge (0 0 0);
diameter -10;
deformationLift 0.005;
liftProfileFile valve1.lift;
staticPoints valve_fixed_points_1;
movingPoints valve_moving_points_1;
movingInternalPoints valve_moving_inernal_points_1;
staticCells valve_fixed_cells_1;
movingCells valve_moving_cells_1;

}

valve2......
.......
);

Last edited by Peter_600; May 23, 2011 at 16:17.
Peter_600 is offline   Reply With Quote

Old   May 24, 2011, 14:41
Default
  #93
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
Hi Peter, thanks a lot!! it is very helpful. I've got a couple of questions though:

When you set a variable to off, does that actually turn it off in the motion solver, so that it won't use it?
Also, what is the difference between the movingPoints and movingInternalPoints for the valve?

and the pistonCellset, pointSet and faceSet, what should they contain? the cells in the piston bowl but not the cylinder maybe?

thanks a lot for your help
abminternet is offline   Reply With Quote

Old   May 24, 2011, 18:36
Default
  #94
Senior Member
 
Join Date: Oct 2009
Posts: 140
Rep Power: 17
Peter_600 is on a distinguished road
Hi
First of all, do you have a mesh which has a similar designed like this?

http://www.engines.polimi.it/images/...eMotionGDI.gif

Which tools did you use to create the mesh? It is crucial to have a folder (named sets) in the polyMesh folder which conists all the different regions. For instance, I create the meshes in Gambit and name the boundaries directly there. Than I use
Code:
FluentMeshToFoam -writeSets
and the sets folder is automatically created with the required files. Do you follow a similiar approach?

Peter
Peter_600 is offline   Reply With Quote

Old   May 27, 2011, 07:19
Default
  #95
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
Hi Peter,

I am following a similar approach, but I am using Ansa for meshing. The thing is, I haven't really done it before and am not sure if I can write out the sets from ansa. So the moving points would be the points moving together with the valve, and the movingInternal points would be the ones were layering is being done on the picture? Is there some kind of documentation on this?
abminternet is offline   Reply With Quote

Old   May 28, 2011, 06:07
Default
  #96
Senior Member
 
Vangelis Skaperdas
Join Date: Mar 2009
Location: Thessaloniki, Greece
Posts: 287
Rep Power: 21
vangelis is on a distinguished road
Hi there,

I would like to help on the ANSA part of the issue.
I assume you are using ANSA v13.1.3
and that you are starting in CFD mode.
If you want to create SETs of nodes or elements
in ANSA for OpenFOAM you should do the following.

Activate the SET button at the top main toolbar
Press NEW
The selection mode window opens
you can mark if you want to select Nodes or Elements
or Properties
Select from the screen the entities you want to load
in the SET and confirm with middle mouse button.
The SETs card opens
Type in the name of the SET and press OK
The SET is created and displayed in the list

Now you can output in OpenFOAM format
In the output options window
activate the Output Sets flag.

ANSA will create the sets folder inside polymesh
as needed

Hope this helps

Vangelis
vangelis is offline   Reply With Quote

Old   May 30, 2011, 07:46
Default
  #97
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
wow great tipp, thanks a lot
abminternet is offline   Reply With Quote

Old   May 30, 2011, 08:10
Default
  #98
Senior Member
 
Vangelis Skaperdas
Join Date: Mar 2009
Location: Thessaloniki, Greece
Posts: 287
Rep Power: 21
vangelis is on a distinguished road
Let me know if you need anything about ANSA
vangelis is offline   Reply With Quote

Old   May 31, 2011, 05:47
Default
  #99
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
Thanks Vangelis, I do have a question, how do I create facesets in ansa? I need a set of internal faces, and I tried with the volToshell function to transform them to a shell and then selecting them in the sets option, but when I tried to write the mesh out, I didnt see the "write sets" option, and indeed the set was not written any ideas? Would appreciate it
abminternet is offline   Reply With Quote

Old   May 31, 2011, 06:39
Default
  #100
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 16
abminternet is on a distinguished road
Hi Peter, I've also got a question for you, hope I'm not causing too much problem. Which faces would you set as the attach/detach faces? I mean, in the test case I did, the faces were at the very entrance of the port, but this faces were not moving. So, is it possible to select faces which are in the valve curtain for example? or faces in the layering zone, etc? Say in the picture you posted, which ones would be the attach/detach faces?
abminternet is offline   Reply With Quote

Reply


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


LinkBacks (?)
LinkBack to this Thread: https://www.cfd-online.com/Forums/openfoam-solving/83177-engine-simulation-mesh-motion-topological-changes.html
Posted By For Type Date
Untitled document This thread Refback February 4, 2014 12:36

Similar Threads
Thread Thread Starter Forum Replies Last Post
Dynamic moving mesh Pei-Ying Hsieh (Hsieh) OpenFOAM Running, Solving & CFD 64 June 7, 2012 11:04
engine simulation with mesh motion and topological changes abminternet OpenFOAM 0 December 16, 2010 12:47
[Commercial meshers] Good mesh for pistoncylinder application Serkan Cetin OpenFOAM Meshing & Mesh Conversion 4 November 3, 2010 08:36
Radiation and miscellaneous enhancements vtk_fan OpenFOAM Running, Solving & CFD 6 February 18, 2008 00:49
Valve action Hrvoje Jasak (Hjasak) OpenFOAM Running, Solving & CFD 0 January 13, 2005 14:23


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