transportProperties as a distribution instead of a scalar

saurcerer · June 7, 2022, 13:11

Hi all,

I want to use heat conductivity as a user input distribution (instead of a scalar) in transportProperties on a node by node basis. Is there a way out to handle this elegantly? setFields sets initial fields as a distribution very nicely, I want to do something like that.

Tobermory · June 9, 2022, 10:30

That's a pretty open question I think - presumably you are interested in a RAS model turbulece closure? The effective thermal diffusivity is calculated as the sum of the laminar and thermal diffusivities, with the former being related to the thermal conductivity, and the latter being calcuilated from the eddy viscosity using the turbulent Prandtl number.

If you want to override all or part of this relationship, then I imagine that you'll need to write your own version of the eddyDiffusivity class and read in the thermal conductivity & Prandtl # as fields ...

Try track down the class logic, starting with the divergence term in the energy equation for the solver that you are working on, and you'll begin to see what you need to duplicate/change.

PS - let us know how you get on!

June 7, 2022, 13:11	transportProperties as a distribution instead of a scalar	#1
saurcerer New Member Aniruddha Saha Join Date: Mar 2022 Posts: 9 Rep Power: 4	Hi all, I want to use heat conductivity as a user input distribution (instead of a scalar) in transportProperties on a node by node basis. Is there a way out to handle this elegantly? setFields sets initial fields as a distribution very nicely, I want to do something like that.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
codeFixedValue - define and use function in code block?	GregNordin	OpenFOAM Running, Solving & CFD	6	November 9, 2023 18:43
[Solved] Passive Scalar Source Term "Bug"	taalf	OpenFOAM Running, Solving & CFD	1	July 13, 2016 09:29
using a Fortran library of thermodynamics inside sutherland transport model	Mehdi3031	OpenFOAM Programming & Development	0	April 7, 2016 10:34
access to member function variable	Rahand	OpenFOAM Programming & Development	3	March 29, 2016 16:46
Specifying nonuniform boundary condition	maka	OpenFOAM Running, Solving & CFD	59	October 22, 2014 15:52

June 9, 2022, 10:30		#2
Tobermory Senior Member Join Date: Apr 2020 Location: UK Posts: 746 Rep Power: 14	That's a pretty open question I think - presumably you are interested in a RAS model turbulece closure? The effective thermal diffusivity is calculated as the sum of the laminar and thermal diffusivities, with the former being related to the thermal conductivity, and the latter being calcuilated from the eddy viscosity using the turbulent Prandtl number. If you want to override all or part of this relationship, then I imagine that you'll need to write your own version of the eddyDiffusivity class and read in the thermal conductivity & Prandtl # as fields ... Try track down the class logic, starting with the divergence term in the energy equation for the solver that you are working on, and you'll begin to see what you need to duplicate/change. PS - let us know how you get on!