[swak4Foam] Manipulate thermophysical properties with swak4foam?

TobM · December 9, 2015, 08:34

Hi

is it possible to manipulate thermophysical properties with swak4foam?

I calculated a temperature dependent viscosity and used the manipulateField function to write the result to "thermo:mu".

The problem is that the solver does not use the new thermo:mu for its calculations. Is there a way to do that?

fabian_roesler · December 11, 2015, 05:41

Why don't you just use polynomial transport properties?
Equation 7.2 in http://cfd.direct/openfoam/user-guide/thermophysical/ should be what you are looking for.

Cheers

Fabian

TobM · December 11, 2015, 06:03

Hi Fabian,

sadly no. Polynomial would only work with a piecewise definition, what is, as far as I know, not possible with existing OF functionality.
What I need is a piecewise definition or alternatively a table where the values are interpolated. Both is possible in swak4foam, what I don't know is, if it's possible to provide the calculated value to the solver.

Greetings

Tobi

fabian_roesler · December 11, 2015, 07:05

You could try to implement a new transport model. The polynomial transport lies here:

src/thermophysicalModels/specie/transport/polynomial/polynomialTransport.H
src/thermophysicalModels/specie/transport/polynomial/polynomialTransport.C

Probably you can change the code to your needs.

Fabian

fabian_roesler · December 11, 2015, 07:08

You could try to merge it with code snippets from the uniformFixedValue boudary condition to create a new tabulated transport class.
I'd say this is the real OpenFOAM masters way

hcl734 · January 5, 2016, 03:29

Did you made any Progress on this issue, since it is interesting for me as well.
I am trying to implement an exponential dependency on temperature for viscosity
http://www.cfd-online.com/Forums/ope...ependance.html

December 9, 2015, 08:34	Manipulate thermophysical properties with swak4foam?	#1
TobM Member Join Date: Sep 2014 Location: Germany Posts: 88 Rep Power: 12	Hi is it possible to manipulate thermophysical properties with swak4foam? I calculated a temperature dependent viscosity and used the manipulateField function to write the result to "thermo:mu". The problem is that the solver does not use the new thermo:mu for its calculations. Is there a way to do that?

December 11, 2015, 07:05		#4
fabian_roesler Senior Member Fabian Roesler Join Date: Mar 2009 Location: Germany Posts: 213 Rep Power: 18	You could try to implement a new transport model. The polynomial transport lies here: src/thermophysicalModels/specie/transport/polynomial/polynomialTransport.H src/thermophysicalModels/specie/transport/polynomial/polynomialTransport.C Probably you can change the code to your needs. Fabian

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Custom Thermophysical Properties	wsmith02	OpenFOAM	4	June 1, 2023 15:30
Transport and Thermophysical Properties	dancfd	OpenFOAM Pre-Processing	0	September 12, 2014 20:52
confused with units in solidThermophysical properties and Thermophysical properties	zahraa	OpenFOAM Running, Solving & CFD	3	August 14, 2014 09:11
table properties for thermophysical properties	romant	OpenFOAM Running, Solving & CFD	1	August 12, 2014 09:41
Thermophysical properties for natural convection	Ciefdi	OpenFOAM Running, Solving & CFD	0	November 7, 2013 12:44

December 11, 2015, 05:41		#2
fabian_roesler Senior Member Fabian Roesler Join Date: Mar 2009 Location: Germany Posts: 213 Rep Power: 18	Why don't you just use polynomial transport properties? Equation 7.2 in http://cfd.direct/openfoam/user-guide/thermophysical/ should be what you are looking for. Cheers Fabian

December 11, 2015, 06:03		#3
TobM Member Join Date: Sep 2014 Location: Germany Posts: 88 Rep Power: 12	Hi Fabian, sadly no. Polynomial would only work with a piecewise definition, what is, as far as I know, not possible with existing OF functionality. What I need is a piecewise definition or alternatively a table where the values are interpolated. Both is possible in swak4foam, what I don't know is, if it's possible to provide the calculated value to the solver. Greetings Tobi

December 11, 2015, 07:08		#5
fabian_roesler Senior Member Fabian Roesler Join Date: Mar 2009 Location: Germany Posts: 213 Rep Power: 18	You could try to merge it with code snippets from the uniformFixedValue boudary condition to create a new tabulated transport class. I'd say this is the real OpenFOAM masters way

January 5, 2016, 03:29		#6
hcl734 Member Join Date: May 2015 Posts: 68 Rep Power: 11	Did you made any Progress on this issue, since it is interesting for me as well. I am trying to implement an exponential dependency on temperature for viscosity http://www.cfd-online.com/Forums/ope...ependance.html