CFD Online URL
[Sponsors]
Home > Forums > CFX

Stirling Motor

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 2, 2013, 10:02
Default Stirling Motor
  #1
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
hallo Everybody

here you can find a simple simulation of stirlig alpha engine.
With the same condition I've made a theorical calculation and I expect
to gain positive work.
With CFX simulation there isn't any positive work

Any advice ?
Attached Files
File Type: zip stirling.zip (4.2 KB, 11 views)
bertozzi_marco is offline   Reply With Quote

Old   September 2, 2013, 19:27
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Please post an image of what you are simulating.

Why do you have gravity enabled? What does it do?

This sounds like an FAQ: http://www.cfd-online.com/Wiki/Ansys..._inaccurate.3F
ghorrocks is offline   Reply With Quote

Old   September 3, 2013, 03:58
Default
  #3
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
Gravity does nothing, it's only for work buoyant model.

in this tread you see the model I l'like to simulate
Simulation of a regenerator in stirling engine

I don't think to find something useful in FAQ, the model give resultrs
very different from theoric calculation, I expect a power of 1.5 / 2 Kw
and the model give -0.2 Kw negative work

If you review results you see the fluid doesn't gain enough heat from bonduary,
and pressure doesn't increase during expansion.

Is known this type of engine work for sure with condition you see in the model
bertozzi_marco is offline   Reply With Quote

Old   September 3, 2013, 20:00
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
What do you mean by "work buoyant model"?

Can you post an image of your geometry, your mesh and the flows you are getting (particularly if you think some part of the flow looks wrong).
ghorrocks is offline   Reply With Quote

Old   September 4, 2013, 06:32
Default
  #5
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
Pressure: http://youtu.be/l-2CPJg-llU

Temperature: http://youtu.be/edIgvUl11rc

and images of geometry and mesh
Attached Images
File Type: jpg mesh.jpg (94.4 KB, 21 views)
File Type: jpg geometry.jpg (83.2 KB, 18 views)
bertozzi_marco is offline   Reply With Quote

Old   September 4, 2013, 07:01
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Your mesh is too coarse. Your heat transfer will be miles off with a mesh that coarse - which is exactly what you have found. If you had gone through the points described in the FAQ you would have worked that out for yourself.
ghorrocks is offline   Reply With Quote

Old   September 4, 2013, 07:14
Default
  #7
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
I will do a simple model to understand only the influence of the coarse mesh
and I will post the results......
thanks for the moment....
bertozzi_marco is offline   Reply With Quote

Old   September 4, 2013, 19:06
Default
  #8
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
here you find simulation with poor mesh: http://youtu.be/ln_CnugS1mI

and with fine mesh: http://youtu.be/xtPTV8xm6Fc

there are some difference of about 100 k at the end of simulation, but the question now is how fine should be mesh?
very fine mesh need a lot of CPU time, poor mesh need less CPU time
but is inacurate, so there are some guidelines to deside appropriate mesh
refinemet ?
bertozzi_marco is offline   Reply With Quote

Old   September 4, 2013, 19:47
Default
  #9
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
This is all discussed in the FAQ I posted right at the start

In short: You need to perform a sensitivity analysis on mesh size. Choose a parameter of importance to you (maybe the work generated) and generate a series of meshes, each with different mesh size. And don't make it different by 10%, the difference needs to be large, like halving the edge length. Do simulations of all these meshes and plot your parameter versus mesh size. Hopefully it will start converging on a value as the mesh is refined. Then you choose an accuracy you are happy to live with and read off the mesh size required to give it.

There are more sophisticated methods of doing this (grid refinement index etc) which are highly recommended and mean you can do this quicker with less simulations - but the concept is the same.

And yes, this will inevitably lead to big simulations, and probably lead to simulations too big to run on a desktop PC. That is why people make supercomputers to run CFD simulations. It is because they have to use a computer that big to run it.
ghorrocks is offline   Reply With Quote

Old   September 5, 2013, 04:26
Default
  #10
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
I understand,
About the model I'd like to explore, I think the heat transfer be the core
parameter for the simulation. Because the simulation of the engine will
take in account a lot of variables, I can use a simple model like you have
seen in the videos before, to explore only the heat transfer phenomena.
Once found the appropriate mesh size I can transfer it to the engine model
and obtain a valid solution with one run. What is your opinion about ?
bertozzi_marco is offline   Reply With Quote

Old   September 5, 2013, 04:28
Default
  #11
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Yes, that is a good plan - Providing the heat transfer is the key parameter in the simulation. Even if that is not the case it is a good start.

Don't forget fluid flow along the pipe and in the entry/exit will also contribute.
ghorrocks is offline   Reply With Quote

Old   September 9, 2013, 05:16
Default
  #12
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
Hallo
here you can find results of simple analisys of a can heated for 3 seconds.
here you can find the immages of the most fine and coarse mesh of
the sensitive analisys, also the table of the sensitive analisys with
the target parameter, the target parameter is the temperature at center
of the can

The strange thing is that the temperature doesn't vary a lot varying the
quality of the mesh, why ?
Attached Images
File Type: jpg Immagine1.jpg (57.3 KB, 9 views)
File Type: jpg Immagine2.jpg (72.8 KB, 11 views)
File Type: jpg Immagine3.jpg (59.4 KB, 10 views)
bertozzi_marco is offline   Reply With Quote

Old   September 9, 2013, 06:53
Default
  #13
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
This is showing that for the case you have modelled you are reasonably accurate. But this case has no flow so the only thing causing heat transfer is conduction - and conduction does not include the tricky non-linear terms which are difficult to get mehs convergence for. This is why you have easily got mesh convergence in this case.

But I think you will find your actual engine case has the majority of the heat transfer from convection effects, and that will be much more sensitive to mesh density. So it is a matter of finding a simple analogy which contains all the important physics but is not too complex to model properly.

I would add an inlet and exhaust port with constant flow to your cylinder. Make the flow rate of the same magnitude as the flows you expect in the engine. Then repeat the mesh sensitivity study on this case and I suspect you will find you will require a considerably finer mesh.
ghorrocks is offline   Reply With Quote

Old   September 9, 2013, 13:01
Default
  #14
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
Here the simulation according your advices,
There is a notable difference of the temperature contour.
But at the outlet the mean temperature isn't so different.

I'd expect a big variation of temperature in the outlet tube,
I don't know if this will vary a lot the work of the engine, is clear the
distribution of temperature is better.

I attach also ccl for your kind check
Attached Images
File Type: jpg Immagine4.jpg (57.8 KB, 6 views)
File Type: jpg Immagine5.jpg (97.2 KB, 5 views)
File Type: jpg Immagine6.jpg (58.4 KB, 6 views)
File Type: jpg Immagine7.jpg (64.7 KB, 6 views)
Attached Files
File Type: zip xxx.zip (3.5 KB, 1 views)
bertozzi_marco is offline   Reply With Quote

Old   September 9, 2013, 18:49
Default
  #15
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
The temperature change is limited to just the boundary layer in the cylinder, so not much has happened yet. You need to run this for more time so the temperature effects convects into the main flow. Better still, run it as a steady state simulation and get the steady state result.
ghorrocks is offline   Reply With Quote

Old   September 11, 2013, 05:48
Default
  #16
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
I have a question about buoyant model, this model
Enable different density of fluid to float over when it is The case,
But if I have fluid speed enough to make the different density negligible
In front of the turbolence due to the speed, should I Enable buoyancy?
bertozzi_marco is offline   Reply With Quote

Old   September 11, 2013, 07:11
Default
  #17
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
The question to ask is actually - given the flow velocities in my model, will buoyancy make a significant difference? If the answer is no (which your last post seems to imply) then you can turn buoyancy off.
ghorrocks is offline   Reply With Quote

Old   September 11, 2013, 09:11
Default
  #18
Member
 
Marco
Join Date: Jul 2013
Location: Italy
Posts: 32
Rep Power: 3
bertozzi_marco is on a distinguished road
yes it is, but I'd like to understand the concept in general,
I'd think the buoyant model will be usefull when the speed of
fluid is very low or comparable with the buoyant turbolence.
So I think that where is fluid driven from high speed volume variation
(like piston engine) the buoyant model could be disabled.
bertozzi_marco is offline   Reply With Quote

Old   September 11, 2013, 09:18
Default
  #19
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 9,880
Rep Power: 78
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
If you are in doubt then activate buoyancy. In forced convection flows it does not add much (if anything) to simulation time. It only adds to simulation time in low speed simulations - but that is because it is doing what it is meant to be doing and generating flows.
ghorrocks 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
Simulation of stirling cryocooler stirling Main CFD Forum 0 March 21, 2012 02:54
Squirrel Cage Induction Motor, Thermal Model?? Manas Agarwal OpenFOAM 0 April 14, 2011 16:40
snappyMesh Motor Bike case for simpleFoam solver s.sivakumar OpenFOAM Pre-Processing 0 July 23, 2009 01:51
snappyMesh Motor Bike case for simpleFoam solver s.sivakumar OpenFOAM Running, Solving & CFD 0 July 22, 2009 09:00
about the simulation of gears air motor brose NUMECA 1 October 20, 2005 12:02


All times are GMT -4. The time now is 17:30.