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

Engine simulation with topological changes

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 3, 2009, 04:58
Default Hi All, my task is to simul
  #1
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 8
marico is on a distinguished road
Hi All,

my task is to simulate an reciprocating engine with one piston and one valve. At least the opening of the valve needs topological changes rather layer addition. Probably the best model for this task should be linearValveLayersFvMesh.
The mesh solver would be displacementLaplacian as I have to define the displacements according to an analytical function of time (piston movement...).
Now I have got a few questions:

1. Is this possible as it is described here? Has anyone experience (testcase...) with that?

2. For what is the variable "pistonVelocity" of linearValveLayersFvMesh? Is it just an approach for calculation? cause velocity changes a lot with time.

3. Whats the easiest way to disable CFD-calculation. I first want to test mesh-deformation... Is there a dummy-solver or a command to disable solving of fluid equations?

I have no experience in C++ so my goal is to avoid dealing much with make... ;)
Thanks for your answers

Marco
marico is offline   Reply With Quote

Old   March 3, 2009, 11:58
Default Hi Marco Have you tried the
  #2
Senior Member
 
Alexandre Pereira
Join Date: Mar 2009
Posts: 155
Rep Power: 8
alexandrepereira is on a distinguished road
Hi Marco

Have you tried the simpleEngineTutorial from Zagreb University OpenFoam page...?

http://powerlab.fsb.hr/ped/kturbo/Op...neTutorial.tgz


I think that using this tutorial and understanding how valve motion is done is a step... but so far I have only managed to get a core dump using engineFoam with this tutorial...

Tried to run it in OpenFOAM-1.5...

Hve you tried it already...?

Best regards

Alex
alexandrepereira is offline   Reply With Quote

Old   March 5, 2009, 04:13
Default Hi Alexandre, thanks for yo
  #3
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 8
marico is on a distinguished road
Hi Alexandre,

thanks for your reply. I looked into it and firstly I'm not able to display the case in paraview. (I'm using 1.5-dev, everything built from sources) Error message when pressing "Apply" to display patches sounds like:

*** glibc detected *** /home.../paraview-real: munmap_chunk(): invalid pointer: 0x0900ab20 ***
...

(other cases work for me)

After looking into the case files I think that the tutorial is out of date in relation to the use of dynamic meshes...?!

Any hints concerning "pistonVelocity" and disabling CFD-calculation?

Best regards
Marco
marico is offline   Reply With Quote

Old   March 5, 2009, 08:00
Default Hi Marco Unfortunately I ha
  #4
Senior Member
 
Alexandre Pereira
Join Date: Mar 2009
Posts: 155
Rep Power: 8
alexandrepereira is on a distinguished road
Hi Marco

Unfortunately I have no clue...

A workable tutorial on moving mesh and dynamicmeshDict is the icoDyMFoam tutotial of movingCone...

I also wish i could learn more on cylinder mesh topological changes in engineFoam... :-(

BRGDS

Alex
alexandrepereira is offline   Reply With Quote

Old   March 5, 2009, 08:26
Default Hello Marco, In order to avoi
  #5
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 647
Rep Power: 11
dmoroian is on a distinguished road
Hello Marco,
In order to avoid solving the flow, you have to remove the part of the code that performs this task. Just open the engineFoam.C file and remove lines from 95 to 119 (this is at least for version 1.4.1 that I have access for the moment). Save the modifications and close the file. Then in the same directory type wclean; wmake.
This should be all!

I hope this is helpful,
Dragos
dmoroian is offline   Reply With Quote

Old   March 6, 2009, 05:53
Default Hi, thank you Dragos for ad
  #6
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 8
marico is on a distinguished road
Hi,

thank you Dragos for advice, I just created "dummyFoam" that way out of icoDyMFoam.
I just tried to modify the simpleEngine Tutorial to work with the actual OF (1.5-dev). When solving I get the message:

Selecting dynamicFvMesh linearValveLayersFvMesh
Time = 0
Adding zones and modifiers to the mesh
Adding point and face zones


Not all zones and patches needed in the definition have been found. Please check your mesh definition.
Error code: 001111

Here's my dynamicMeshDict:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
object dynamicMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMeshLib "libtopoChangerFvMesh.so";

dynamicFvMesh linearValveLayersFvMesh;

solver displacementLaplacian;

diffusivity uniform;

linearValveLayersFvMeshCoeffs
{
pistonVelocity (0 0 1);

slider
{
inside valveCurtainPort1 valveCurtainPort2;
outside valveCurtainCyl1 valveCurtainCyl2;
}

layer
{
patch piston valve2 valve4 cylinderHead;
}

minThickness 0.001;
maxThickness 0.003;
}


Maybe someone is able to find the error... I also attach my case:



Note: BC's and solver settings are probably not well, the problem occured seems to be independent from that.


Thank You
Marco
marico is offline   Reply With Quote

Old   March 6, 2009, 06:13
Default http://www.cfd-online.com/Ope
  #7
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 8
marico is on a distinguished road
engine.tar
marico is offline   Reply With Quote

Old   March 9, 2009, 02:56
Default Hi again, one problem I fou
  #8
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 8
marico is on a distinguished road
Hi again,

one problem I found out is probably that the patches defined in the dynamicMeshDict above to the keywords "slider" and "layer" are not all recognized. Only the first one, e.g. "piston", is recognized. Is there a way to define multiple patches to these keywords? Another way would be to combine the patches in mesh setup but that would lead to another problem: I'm not any longer able to define different boundary conditions (which I need to) to them.
Any ideas?

Thank You
Marco
marico is offline   Reply With Quote

Old   March 9, 2009, 10:48
Default Hi once more, I created a s
  #9
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 8
marico is on a distinguished road
Hi once more,

I created a simple testcase with BlockMesh geometry. The geometry has some interfaces. Two of them are for the slider. Others result out of my "block choice". To use the mesh I have to merge these internal patches?!
After doing this OF creates patchZones and pointZones files.
Then I try to run the dynamic mesh solver and it complains about meshmodifiers. (a meshModifiers file isn't created, XXXZones files aren't updated due to the dynamicMeshDict)
May it be a bug (sorry for that formulation ;) that the dynamic mesh forgets to create the XXXZones files it needs, because of the fact, that in polyMesh-directory exist similar files resulting from the merge-process???
Or am I just to stupid to use the right tools?

Thank You
Marco
marico is offline   Reply With Quote

Old   March 10, 2009, 06:50
Default With silent cries and hope for
  #10
Member
 
Marco Müller
Join Date: Mar 2009
Location: Germany
Posts: 94
Rep Power: 8
marico is on a distinguished road
With silent cries and hope for attention I'm waiting for someone who knows how to deal with the linearValve(Layers)FvMesh...
So far I've created the meshModifiers and zones files manually and they seem to work.
When trying to run mesh deformation it appears:

...
Time = 0
Sliding interfaces decoupled
Coupling sliding interfaces
segmentation faultī

ARRGH.... 8-/

As described before my mesh consists of 3 regions... 2 of them form the "surrounding mesh" and one cylinder inside has to move...
I furthermore tried two other examples from the message board using the linearValveThings and they run.

And what is the *** pistonVelocity for???

Thank You
Marco
marico is offline   Reply With Quote

Old   March 10, 2009, 13:10
Default Dear Marco, I am following
  #11
New Member
 
Stefan Vogel
Join Date: Mar 2009
Posts: 1
Rep Power: 0
svogel is on a distinguished road
Dear Marco,

I am following the thread. I just know that
Tommaso is using this kind of technique.

tommaso.lucchini--at--polimi.it

Where are you located in Germany, I also want to be able to generate engine grids. (Tried it already with simple deformation using snappyhexmesh).

You could write me an email to zwitschervogel--at--web.de

kind regards

Stefan
svogel is offline   Reply With Quote

Old   March 11, 2009, 09:33
Default Dear Marco, If you need mor
  #12
Member
 
Tommaso Lucchini
Join Date: Mar 2009
Posts: 82
Rep Power: 8
lucchini is on a distinguished road
Dear Marco,

If you need more information or support about simulations of IC engines in OpenFOAM with valve and piston motion + topological changes, please contact me at my e-mail address.

Regards

Tommaso
lucchini is offline   Reply With Quote

Old   April 14, 2013, 07:18
Default
  #13
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Iran - Tehran
Posts: 185
Rep Power: 4
sasanghomi is on a distinguished road
Dear Dr.Luchinni,

I am simulating an engine with sonicTurbDyMEngineFoam and simpleEngineTopoFvMesh library ..
I imported the mesh from GAMBIT to openfoam But unfortunately when I ran the case I had an error Motion continuity errors ...

Starting time loop Courant Number mean: 0 max: 0 velocity magnitude: 0 deltaT = 2.6455e-06 Crank angle = 0.0238095 CA-deg CG: Solving for motionUx, Initial residual = 0, Final residual = 0, No Iterations 0 CG: Solving for motionUy, Initial residual = 0, Final residual = 0, No Iterations 0 CG: Solving for motionUz, Initial residual = 0.213596, Final residual = 0.000179652, No Iterations 20 CG: Solving for motionUx, Initial residual = 4.89248e-09, Final residual = 4.42683e-10, No Iterations 1 CG: Solving for motionUy, Initial residual = 0, Final residual = 0, No Iterations 0 CG: Solving for motionUz, Initial residual = 0.000100573, Final residual = 7.10585e-08, No Iterations 21 Volume: new = 48492.4 old = 48492.3 change = 0.0719435 Motion continuity errors : sum local = 1.22866e-15, maximum = 2.88108e-14 Floating point exception (core dumped)

How can I overcome this error?
I need your help...I sent you an email but unfortunately you didn't answer me.
I appreciate any help from you..

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

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
NEED HELP on IC engine simulation with AVL FIRE! Mohsen khalafi Main CFD Forum 2 February 6, 2014 08:47
IC Engine Simulation Anubhav Sinha FLUENT 6 July 28, 2011 02:58
simulation oil in combustion engine ? Berg FLUENT 0 April 23, 2006 09:25
IC Engine simulation in Fluent Ravi Kiran FLUENT 0 March 6, 2006 00:48
IC Engine Flow simulation Apurv Main CFD Forum 2 February 8, 2000 01:12


All times are GMT -4. The time now is 21:43.