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

Temperature depending density

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 23, 2016, 05:20
Post Temperature depending density
  #1
New Member
 
Lukas
Join Date: May 2016
Location: Germany
Posts: 19
Rep Power: 7
Birnensuppe is on a distinguished road
Hey,
I am struggling at the moment on the implementation of temperature depended density. I created a user function with 100 K temperature steps. (My temperature range is from 300 K to 30 000K).
The problem is, that the solver has fatal linear overflow after 5 iterations.
I also tried to implement a function which matches the density, but with 10 K temperature steps, this also creates a fatal linear overflow.
Then I modified the function so that for a temperature below 1000 K the density is constant (I know this is not physical valid, it was just a test), the simulations runs fine until the temperature exceeds 1000 K.

ERROR #004100018 has occurred in subroutine FINMES.
Message: Fatal overflow in linear solver.



Some also faced this problem:
http://www.cfd-online.com/Forums/cfx...efinition.html


Rearranging the table doesn’t works for me and also my mathematical function doesn’t work.
I you need further information about my setup, please ask.
Thank in advance for your response
Best Greetings
Lukas
Birnensuppe is offline   Reply With Quote

Old   June 23, 2016, 05:27
Default
  #2
Senior Member
 
Lance
Join Date: Mar 2009
Posts: 669
Rep Power: 19
Lance is on a distinguished road
um.... 30 000 K is pretty much, having in mind that the surface temperature of the sun is roughly 5000 K. What is your application? What physics are you trying to model? What are the actual values of the density of your fluid above 1000 K? Are you sure the continuum assumption is valid?
It seems that CFX is having a hard time dealing with your density function, and you probably need to model whatever you are modeling in another way.
Lance is offline   Reply With Quote

Old   June 23, 2016, 05:52
Default
  #3
Member
 
Join Date: Dec 2009
Posts: 44
Rep Power: 13
cfdgremlin is on a distinguished road
You need to be very careful when developing your own functions/expressions for material properties, such as density. Strictly speaking, they need to be thermodynamically consistent in order for them to work properly and produce correct solutions.

My advice would be to, if possible, use one of the existing material models to define your fluid, as these are tried and tested. If none of these fit your requirement, and you need to develop your own expressions, then there is some documentation at the end of the section on Material Properties which discusses the criteria needed for thermodynamic consistency.
cfdgremlin is offline   Reply With Quote

Old   June 23, 2016, 06:30
Default
  #4
New Member
 
Lukas
Join Date: May 2016
Location: Germany
Posts: 19
Rep Power: 7
Birnensuppe is on a distinguished road
@Lance: Yes I know. I am simulation eletric arcs with special look at temperatur radiation. I hope the maximum temperature is around 15 000 K, but the density function is implemented in a range from 300K to 30 000K.
I take the density values from Gleizes' publication: http://iopscience.iop.org/article/10...9/5/055013/pdf

For the testfile I set the density below 1000 K to 3.4692E-01 [kg/m^3] (This is the density at 1000K) after that I followed publication. Are you interessted that I post all missing 29 000 values?
For 10 000 K the density is 1.7218E-02 [kg/m^3] for example.

Yes, I am sure the continuum assumption is valid. Thank you for this thaught.

@cfdgremlin: Thank you for your advice, since I take the values from an trustworthy publication, I hope that the values are correct.

To specify, I startet with ideal gas provided by CFX, then I add temperatur depended conductivity, dynamic viscosity, specific heat capacity and thermal conductivity, it's all working fine, but not the density.

In the solver theory guide it is often mentioned, that the density can be variable. Is it possible that CFX calculated the density by itself? (I tried to find it, but wasn't sucessful)

Kind regards
Lukas
Birnensuppe is offline   Reply With Quote

Old   June 23, 2016, 07:18
Default
  #5
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
The original question is an FAQ: http://www.cfd-online.com/Wiki/Ansys...do_about_it.3F

And as the others have stated it is almost certainly linked to your extreme conditions and/or the material properties linked to it.

Note that beyond a few thousand Kelvin most materials are plasma and CFX does not have a plasma model. CFX's capabilities in the plasma regime are quite limited. So CFX may well be incapable of modelling what you are intending.
ghorrocks is offline   Reply With Quote

Old   June 24, 2016, 03:17
Default
  #6
New Member
 
Lukas
Join Date: May 2016
Location: Germany
Posts: 19
Rep Power: 7
Birnensuppe is on a distinguished road
Thank you for your reply ghorrocks.
My question is, whether it is possible to implement a temperature dependent density.
The fact, that the solver runs fine with a constant density at a value for 300 K or for 15 000 K and with a modified density which is constant. As soon as the setup exceeds the temperature range where the density is hold constant by the modified function, the solver is diverging. The difference of the value for the density from 1000 K to 1100 K is 3.4692E-01 - 3.1538E-01 = 0.3154E-01 [kg/m^3], which shouldn't be such a big step, which could disturb the solver. This matter astonish me.
In addition I don't think CFX should have a problem to simulate the setup with a temperature around 300 K. I hope you agree with me in this point. That makes me think that either I made a basic mistake or CFX has a basic failure. (I hope it's me )
Birnensuppe is offline   Reply With Quote

Old   June 24, 2016, 06:35
Default
  #7
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 16,840
Rep Power: 132
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
Yes, you can define a temperature dependant density.

Are you using ideal gas or some other compressible EOS? If not, why not? And if not, what EOS are you using?
ghorrocks is offline   Reply With Quote

Old   June 24, 2016, 08:47
Default
  #8
New Member
 
Lukas
Join Date: May 2016
Location: Germany
Posts: 19
Rep Power: 7
Birnensuppe is on a distinguished road
I am using ideal gas EOS with the option value, because I couldn't find the necessary material parameters for ideal gas EOS, which are accurate for my temperature range.
If I change the EOS to ideal gas, CFX calculates the density with the ideal gas law, but that not exact enough for my temperature range. There I want to implement an density, which is more accurate for my setu
I have attached my material setup for the density.
Attached Images
File Type: jpg CFD_BASIC_SETTINGS.JPG (27.6 KB, 32 views)
File Type: jpg CFD_Material_Properties.JPG (73.0 KB, 33 views)
Birnensuppe is offline   Reply With Quote

Reply

Tags
density, variable density

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
Ansys CFX problem: unexpected very high temperatures in premix laminar combustion faizan_habib7 CFX 4 February 1, 2016 17:00
Difficulty In Setting Boundary Conditions Moinul Haque CFX 4 November 25, 2014 17:30
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Is there a way to set a CEL expression to the density of a material at a temperature? audrey CFX 11 July 19, 2012 07:51
the material parameters depending on temperature Jason CFX 4 December 22, 2006 19:54


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