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

Conjugate heat transfer high speed gas flow

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 1 Post By jcoelho5
  • 1 Post By jcoelho5

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 6, 2015, 04:20
Default Conjugate heat transfer high speed gas flow
  #1
New Member
 
Jakob Trydal
Join Date: Mar 2014
Posts: 1
Rep Power: 0
jaktry89 is on a distinguished road
Hi!
First, I must say that this forum has been of great help so far. Many of the problems I have encountered have been solved by reading different threads.

Now I hope someone can help point me in the right direction with my current problem.
I am trying to simulate pressurization of a gas pipeline. This involves high speed gas flow through an orifice. The geometry looks something like this:



The input pressure is 235 bar, and the gas temperature is 260 K. Due to the high pressure, very high flowrates occur.

I am interested in the temperature distribution in the pipe material. If the temperature drops below a certain limit, the material can become brittle and fail at much lower stresses than normal.

I have tried to use chtMultiRegionFoam, and I have managed to create a mesh with one fluid and one solid region. However, it has come to my attention that this solver might not be appropriate for the problem. As far as I have understood, even though it is a compressible solver, it assumes that rho is not a function of pressure, only of temperature.

One possible solution I can think of is replacing the fluid solver in chtMultiRegionFoam with rhoCentralFoam or sonicFoam. However, I don’t think my programming skills is sufficient to do something like that. I have tried to search the forum to see if something like this has been done, but I was unable to find anything.

I am wondering if it is possible to solve the fluid and pipe regions with different solvers (e.g. rhoCentralFoam and laplacianFoam), and somehow couple these, without creating a new solver?
Maybe this is impossible, but any help would be appreciated!

Best Regards
Jakob
jaktry89 is offline   Reply With Quote

Old   November 10, 2016, 12:06
Default
  #2
New Member
 
Alaa Bejaoui
Join Date: Oct 2016
Posts: 2
Rep Power: 0
AlaaBejaoui is on a distinguished road
Hello Jakob,

I have exactly the same Problem. Did you find a solution to your problem? I’dbeverythankful

Best Regardsfrom Germany
AlaaBejaoui is offline   Reply With Quote

Old   December 29, 2016, 04:31
Default
  #3
Member
 
Join Date: Sep 2016
Posts: 63
Rep Power: 9
sitajeje is on a distinguished road
Dear Jakob,

I am pretty new in OpenFOAM, but if the perfectGas is chosen for the equation of state, is the density still pressure independent? Your experience would be very helpful for me. Thank you very much in advance!

Sitajeje
sitajeje is offline   Reply With Quote

Old   December 29, 2016, 11:34
Default
  #4
New Member
 
Dawid Gacek
Join Date: Nov 2015
Location: Poland
Posts: 4
Rep Power: 10
DawidG is on a distinguished road
Hello everyone!

@sitajeje - If you chose perfectGas, density is dependent of pressure and temperature. Here you can read more about that
http://cfd.direct/openfoam/user-guide/thermophysical/
DawidG is offline   Reply With Quote

Old   May 18, 2021, 05:17
Default
  #5
New Member
 
Alosha Lovett
Join Date: Apr 2021
Posts: 4
Rep Power: 5
Alosha is on a distinguished road
Hi Jakob,

I am doing something very similar however have approached it using sonicFoam.

Did you ever manage to progress at all? I am a little stuck on how to implement appropriate heat transfer between my gas and bounding geometry.

Thanks,
Alex.
Alosha is offline   Reply With Quote

Old   August 30, 2021, 12:40
Default Tabulated thermophysical properties
  #6
New Member
 
Joao Coelho
Join Date: Jun 2021
Posts: 23
Rep Power: 4
jcoelho5 is on a distinguished road
For future reference, the new OpenFOAM 9 allows any option for thermophysical properties, which do not limit us in what choice to use for each property.

I was trying to simulate the flow of cryogenic fuel inside a pipe, considering all thermodynamic properties changes.
Being so, I was able to use tabulated properties for every property of interest.

transport tabulated;
thermo hTabulated;
equationOfState rhoTabulated;

thermodynamics
{
Hf 0;
Sf 0;
Hs
{
low (Plow Tlow);
high (Phigh Thigh);
values
Np Nt
(
( ... ... ... ...) // for Plow, Nt values of Hs, For Nt uniform distributed values of T between Tlow and Thigh
.
.
.
//repeat Np times (for Np uniform distributed values of P between Plow and Phigh)
.
.
.
( ... ... ... ...) // for Phigh, Nt values of Hs, For Nt uniform distributed values of T between Tlow and Thigh
);
}


.... Repeat table for Cp and Cv
}

transport
{
... Repeat table for mu and kappa
}

equationOfState
{
... Repeat table for rho
}

You can do all this by retrieving data from CoolProp or REFPROP for example.

Nevertheless, I had some problems at the beginning of the run using this setup. It turned out it was a consequence of fast changes in the properties in the first time steps. I advise you to start with constant properties (rho, Cp, etc) until you reach at least a physical solution and then upgrade the solution with these changes.

The OF page for these models (tabulated, hTabulated and rhoTabulated) did not mention the use of the table in this format. I realise it was like that looking at the source code. I believe they forgot to update their description.

Hope this can help anyone trying to use these features in the future.

Disclaimer: I am not an expert in CFD, neither in OpenFOAM. I simply tried to do this by trial and error. If you have any questions, feel free to post them here, but the most probable thing is I won't know how to answer if I did not have the same problem. But leave it here anyway, more knowledgeable people may help you.


P.S. - I had run this with buoyantPimpleFoam/chtMultiRegion (which requires psiThermo). Recently, I tried using the same feature for rhoCentralFoam (which requires rhoThermo) but it does not work, and it said I needed to use pengRobinson or perfectGas for EoS.

After some research, I found this in the commit description:
The thermophysical property package is constructed automatically from the etc/codeTemplates/dynamicCode files for the corresponding base thermo type (e.g. fluidThermo), currently these are provided only for fluidThermo but the others will be added shortly.
https://github.com/OpenFOAM/OpenFOAM...5cfb4b529bb087

I guess it was not added yet for rhoThermo.
jherb likes this.

Last edited by jcoelho5; August 31, 2021 at 07:34. Reason: new info
jcoelho5 is offline   Reply With Quote

Old   January 21, 2022, 06:26
Default
  #7
Member
 
Stanley John
Join Date: Sep 2018
Posts: 79
Rep Power: 7
sjohn2 is on a distinguished road
Can yo kindly post your thermophysical file?
It am still bit unsure.
sjohn2 is offline   Reply With Quote

Old   January 21, 2022, 09:50
Default
  #8
New Member
 
Joao Coelho
Join Date: Jun 2021
Posts: 23
Rep Power: 4
jcoelho5 is on a distinguished road
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport tabulated;
thermo hTabulated;
equationOfState rhoTabulated;

specie specie;
energy sensibleEnthalpy;
}

mixture
{
specie
{
molWeight 16.22239;
}
#include "thermoProp_tables_FUEL"
}




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


in another file (constant/thermoProp_tables_FUEL), I have the following:


(changing linearly from 50 bar to 100 bar; 75K to 925K)

HTML Code:
thermodynamics

{

	Hf	0;

	Sf	0;

	Hs

	{

		low		(5.000e+06 75.0);

		high	(1.000e+07 925.0);

		values

		11 18

		(

			(-1.178803e+05 5.104859e+04 2.422008e+05 6.487439e+05 7.939995e+05 9.227508e+05 1.051353e+06 1.184809e+06 1.328738e+06 1.477549e+06 1.634881e+06 1.800852e+06 1.975394e+06 2.158334e+06 2.349433e+06 2.548414e+06 2.754973e+06 2.968791e+06 )

			(-1.170362e+05 5.171680e+04 2.413779e+05 6.356783e+05 7.875740e+05 9.184711e+05 1.048230e+06 1.182433e+06 1.326888e+06 1.476098e+06 1.633741e+06 1.799962e+06 1.974708e+06 2.157819e+06 2.349063e+06 2.548168e+06 2.754834e+06 2.968747e+06 )

			(-1.161919e+05 5.238880e+04 2.406817e+05 6.218159e+05 7.810957e+05 9.142002e+05 1.045127e+06 1.180078e+06 1.325056e+06 1.474663e+06 1.632616e+06 1.799084e+06 1.974033e+06 2.157313e+06 2.348701e+06 2.547929e+06 2.754702e+06 2.968709e+06 )

			(-1.153474e+05 5.306449e+04 2.400937e+05 6.071368e+05 7.745729e+05 9.099414e+05 1.042045e+06 1.177743e+06 1.323243e+06 1.473245e+06 1.631504e+06 1.798218e+06 1.973369e+06 2.156818e+06 2.348348e+06 2.547698e+06 2.754577e+06 2.968676e+06 )

			(-1.145027e+05 5.374375e+04 2.395995e+05 5.916930e+05 7.680153e+05 9.056980e+05 1.038984e+06 1.175429e+06 1.321448e+06 1.471843e+06 1.630407e+06 1.797365e+06 1.972716e+06 2.156331e+06 2.348003e+06 2.547474e+06 2.754458e+06 2.968649e+06 )

			(-1.136579e+05 5.442646e+04 2.391877e+05 5.756641e+05 7.614343e+05 9.014736e+05 1.035947e+06 1.173138e+06 1.319672e+06 1.470458e+06 1.629325e+06 1.796525e+06 1.972074e+06 2.155854e+06 2.347666e+06 2.547258e+06 2.754346e+06 2.968629e+06 )

			(-1.128129e+05 5.511252e+04 2.388487e+05 5.594160e+05 7.548431e+05 8.972718e+05 1.032936e+06 1.170869e+06 1.317916e+06 1.469090e+06 1.628257e+06 1.795697e+06 1.971443e+06 2.155387e+06 2.347338e+06 2.547049e+06 2.754241e+06 2.968614e+06 )

			(-1.119677e+05 5.580184e+04 2.385747e+05 5.435082e+05 7.482566e+05 8.930967e+05 1.029950e+06 1.168623e+06 1.316180e+06 1.467739e+06 1.627204e+06 1.794882e+06 1.970822e+06 2.154930e+06 2.347018e+06 2.546847e+06 2.754142e+06 2.968604e+06 )

			(-1.111224e+05 5.649432e+04 2.383590e+05 5.285828e+05 7.416911e+05 8.889520e+05 1.026993e+06 1.166402e+06 1.314464e+06 1.466405e+06 1.626165e+06 1.794079e+06 1.970213e+06 2.154481e+06 2.346706e+06 2.546653e+06 2.754049e+06 2.968601e+06 )

			(-1.102770e+05 5.718987e+04 2.381960e+05 5.151530e+05 7.351646e+05 8.848417e+05 1.024065e+06 1.164204e+06 1.312769e+06 1.465088e+06 1.625141e+06 1.793289e+06 1.969615e+06 2.154043e+06 2.346403e+06 2.546466e+06 2.753963e+06 2.968603e+06 )

			(-1.094314e+05 5.788839e+04 2.380811e+05 5.034469e+05 7.286961e+05 8.807701e+05 1.021168e+06 1.162032e+06 1.311094e+06 1.463789e+06 1.624132e+06 1.792512e+06 1.969027e+06 2.153614e+06 2.346108e+06 2.546286e+06 2.753884e+06 2.968611e+06 )

		);

	}

	Cp

	{

		low		(5.000e+06 75.0);

		high	(1.000e+07 925.0);

		values

		11 18

		(

			(3.363368e+03 3.500509e+03 4.680558e+03 3.432670e+03 2.642175e+03 2.548970e+03 2.610251e+03 2.734956e+03 2.893184e+03 3.060625e+03 3.233022e+03 3.405569e+03 3.575525e+03 3.741261e+03 3.901774e+03 4.056434e+03 4.204853e+03 4.346815e+03 )

			(3.362161e+03 3.493151e+03 4.571544e+03 3.733175e+03 2.704999e+03 2.578750e+03 2.628339e+03 2.747325e+03 2.902316e+03 3.067615e+03 3.238557e+03 3.410064e+03 3.579248e+03 3.744396e+03 3.904448e+03 4.058741e+03 4.206862e+03 4.348580e+03 )

			(3.360974e+03 3.485994e+03 4.477441e+03 4.092451e+03 2.770596e+03 2.608976e+03 2.646494e+03 2.759675e+03 2.911413e+03 3.074564e+03 3.244054e+03 3.414526e+03 3.582944e+03 3.747506e+03 3.907101e+03 4.061029e+03 4.208856e+03 4.350332e+03 )

			(3.359806e+03 3.479030e+03 4.395113e+03 4.516157e+03 2.838809e+03 2.639591e+03 2.664695e+03 2.771996e+03 2.920470e+03 3.081471e+03 3.249512e+03 3.418954e+03 3.586610e+03 3.750591e+03 3.909733e+03 4.063300e+03 4.210834e+03 4.352069e+03 )

			(3.358656e+03 3.472251e+03 4.322282e+03 4.999728e+03 2.909405e+03 2.670531e+03 2.682920e+03 2.784277e+03 2.929482e+03 3.088332e+03 3.254930e+03 3.423348e+03 3.590247e+03 3.753652e+03 3.912343e+03 4.065552e+03 4.212796e+03 4.353794e+03 )

			(3.357525e+03 3.465649e+03 4.257248e+03 5.516967e+03 2.982058e+03 2.701726e+03 2.701143e+03 2.796509e+03 2.938444e+03 3.095146e+03 3.260307e+03 3.427706e+03 3.593854e+03 3.756687e+03 3.914932e+03 4.067786e+03 4.214743e+03 4.355504e+03 )

			(3.356412e+03 3.459218e+03 4.198709e+03 6.009435e+03 3.056345e+03 2.733097e+03 2.719340e+03 2.808683e+03 2.947351e+03 3.101910e+03 3.265641e+03 3.432028e+03 3.597431e+03 3.759697e+03 3.917500e+03 4.070001e+03 4.216673e+03 4.357201e+03 )

			(3.355317e+03 3.452949e+03 4.145653e+03 6.393163e+03 3.131741e+03 2.764560e+03 2.737486e+03 2.820787e+03 2.956198e+03 3.108622e+03 3.270931e+03 3.436314e+03 3.600977e+03 3.762681e+03 3.920046e+03 4.072199e+03 4.218588e+03 4.358885e+03 )

			(3.354240e+03 3.446838e+03 4.097274e+03 6.597493e+03 3.207623e+03 2.796026e+03 2.755554e+03 2.832812e+03 2.964980e+03 3.115279e+03 3.276176e+03 3.440563e+03 3.604493e+03 3.765640e+03 3.922570e+03 4.074378e+03 4.220488e+03 4.360554e+03 )

			(3.353180e+03 3.440877e+03 4.052926e+03 6.611357e+03 3.283271e+03 2.827400e+03 2.773518e+03 2.844748e+03 2.973691e+03 3.121878e+03 3.281374e+03 3.444775e+03 3.607978e+03 3.768572e+03 3.925073e+03 4.076538e+03 4.222371e+03 4.362211e+03 )

			(3.352138e+03 3.435061e+03 4.012081e+03 6.483608e+03 3.357892e+03 2.858585e+03 2.791352e+03 2.856585e+03 2.982327e+03 3.128419e+03 3.286526e+03 3.448948e+03 3.611431e+03 3.771479e+03 3.927554e+03 4.078680e+03 4.224240e+03 4.363854e+03 )

		);

	}

	Cv

	{

		low		(5.000e+06 75.0);

		high	(1.000e+07 925.0);

		values

		11 18

		(

			(2.337519e+03 2.005374e+03 1.880687e+03 1.760242e+03 1.722873e+03 1.809590e+03 1.949279e+03 2.116728e+03 2.298235e+03 2.483265e+03 2.667948e+03 2.849435e+03 3.026106e+03 3.197016e+03 3.361599e+03 3.519517e+03 3.670584e+03 3.814727e+03 )

			(2.338665e+03 2.006754e+03 1.877775e+03 1.784391e+03 1.730230e+03 1.813340e+03 1.951619e+03 2.118341e+03 2.299420e+03 2.484160e+03 2.668640e+03 2.849980e+03 3.026540e+03 3.197365e+03 3.361881e+03 3.519745e+03 3.670768e+03 3.814875e+03 )

			(2.339805e+03 2.008129e+03 1.875521e+03 1.809042e+03 1.737548e+03 1.817046e+03 1.953930e+03 2.119936e+03 2.300593e+03 2.485047e+03 2.669328e+03 2.850523e+03 3.026974e+03 3.197714e+03 3.362163e+03 3.519973e+03 3.670953e+03 3.815024e+03 )

			(2.340938e+03 2.009497e+03 1.873790e+03 1.833482e+03 1.744795e+03 1.820704e+03 1.956210e+03 2.121511e+03 2.301755e+03 2.485927e+03 2.670012e+03 2.851063e+03 3.027406e+03 3.198062e+03 3.362445e+03 3.520202e+03 3.671138e+03 3.815173e+03 )

			(2.342064e+03 2.010860e+03 1.872481e+03 1.856547e+03 1.751937e+03 1.824307e+03 1.958459e+03 2.123068e+03 2.302905e+03 2.486800e+03 2.670691e+03 2.851600e+03 3.027836e+03 3.198409e+03 3.362726e+03 3.520431e+03 3.671324e+03 3.815323e+03 )

			(2.343184e+03 2.012217e+03 1.871518e+03 1.876574e+03 1.758939e+03 1.827851e+03 1.960675e+03 2.124606e+03 2.304045e+03 2.487666e+03 2.671366e+03 2.852135e+03 3.028265e+03 3.198756e+03 3.363008e+03 3.520660e+03 3.671510e+03 3.815473e+03 )

			(2.344297e+03 2.013568e+03 1.870842e+03 1.891655e+03 1.765763e+03 1.831331e+03 1.962859e+03 2.126125e+03 2.305173e+03 2.488526e+03 2.672037e+03 2.852668e+03 3.028693e+03 3.199102e+03 3.363289e+03 3.520889e+03 3.671696e+03 3.815625e+03 )

			(2.345404e+03 2.014912e+03 1.870406e+03 1.900366e+03 1.772371e+03 1.834743e+03 1.965008e+03 2.127625e+03 2.306290e+03 2.489378e+03 2.672704e+03 2.853198e+03 3.029120e+03 3.199448e+03 3.363571e+03 3.521118e+03 3.671883e+03 3.815776e+03 )

			(2.346504e+03 2.016251e+03 1.870172e+03 1.902632e+03 1.778724e+03 1.838081e+03 1.967123e+03 2.129105e+03 2.307395e+03 2.490224e+03 2.673367e+03 2.853726e+03 3.029545e+03 3.199793e+03 3.363852e+03 3.521348e+03 3.672071e+03 3.815929e+03 )

			(2.347599e+03 2.017583e+03 1.870110e+03 1.899848e+03 1.784788e+03 1.841343e+03 1.969203e+03 2.130566e+03 2.308489e+03 2.491063e+03 2.674026e+03 2.854252e+03 3.029969e+03 3.200138e+03 3.364134e+03 3.521578e+03 3.672259e+03 3.816082e+03 )

			(2.348687e+03 2.018909e+03 1.870194e+03 1.894030e+03 1.790528e+03 1.844523e+03 1.971246e+03 2.132008e+03 2.309572e+03 2.491896e+03 2.674681e+03 2.854776e+03 3.030392e+03 3.200482e+03 3.364415e+03 3.521809e+03 3.672447e+03 3.816236e+03 )

		);

	}

}



transport

{

	mu

	{

		low		(5.000e+06 75.0);

		high	(1.000e+07 925.0);

		values

		11 18

		(

			(4.581839e-04 9.894276e-05 4.098806e-05 1.037836e-05 1.148123e-05 1.286222e-05 1.424161e-05 1.558284e-05 1.688270e-05 1.813810e-05 1.935452e-05 2.053499e-05 2.168239e-05 2.279934e-05 2.388818e-05 2.495100e-05 2.598965e-05 2.700579e-05 )

			(4.610169e-04 9.950840e-05 4.171352e-05 1.076275e-05 1.164671e-05 1.297667e-05 1.433268e-05 1.566000e-05 1.695053e-05 1.819876e-05 1.940961e-05 2.058557e-05 2.172921e-05 2.284297e-05 2.392907e-05 2.498949e-05 2.602602e-05 2.704027e-05 )

			(4.638521e-04 1.000730e-04 4.240857e-05 1.121877e-05 1.182464e-05 1.309640e-05 1.442668e-05 1.573901e-05 1.701969e-05 1.826038e-05 1.946542e-05 2.063671e-05 2.177649e-05 2.288698e-05 2.397027e-05 2.502824e-05 2.606261e-05 2.707494e-05 )

			(4.666897e-04 1.006367e-04 4.307752e-05 1.175900e-05 1.201539e-05 1.322144e-05 1.452360e-05 1.581987e-05 1.709017e-05 1.832294e-05 1.952194e-05 2.068840e-05 2.182420e-05 2.293135e-05 2.401177e-05 2.506724e-05 2.609942e-05 2.710980e-05 )

			(4.695298e-04 1.011994e-04 4.372380e-05 1.239602e-05 1.221930e-05 1.335181e-05 1.462343e-05 1.590255e-05 1.716196e-05 1.838645e-05 1.957916e-05 2.074064e-05 2.187236e-05 2.297608e-05 2.405357e-05 2.510650e-05 2.613645e-05 2.714486e-05 )

			(4.723725e-04 1.017611e-04 4.435016e-05 1.313933e-05 1.243665e-05 1.348751e-05 1.472615e-05 1.598705e-05 1.723505e-05 1.845087e-05 1.963707e-05 2.079342e-05 2.192095e-05 2.302116e-05 2.409566e-05 2.514601e-05 2.617370e-05 2.718010e-05 )

			(4.752179e-04 1.023219e-04 4.495888e-05 1.399022e-05 1.266762e-05 1.362854e-05 1.483173e-05 1.607334e-05 1.730942e-05 1.851622e-05 1.969568e-05 2.084673e-05 2.196997e-05 2.306660e-05 2.413805e-05 2.518577e-05 2.621116e-05 2.721553e-05 )

			(4.780661e-04 1.028818e-04 4.555186e-05 1.493606e-05 1.291232e-05 1.377487e-05 1.494016e-05 1.616140e-05 1.738507e-05 1.858248e-05 1.975497e-05 2.090058e-05 2.201941e-05 2.311238e-05 2.418073e-05 2.522578e-05 2.624883e-05 2.725114e-05 )

			(4.809172e-04 1.034408e-04 4.613072e-05 1.594826e-05 1.317078e-05 1.392647e-05 1.505140e-05 1.625122e-05 1.746197e-05 1.864963e-05 1.981493e-05 2.095494e-05 2.206927e-05 2.315850e-05 2.422369e-05 2.526602e-05 2.628671e-05 2.728693e-05 )

			(4.837714e-04 1.039989e-04 4.669681e-05 1.698867e-05 1.344287e-05 1.408329e-05 1.516542e-05 1.634277e-05 1.754013e-05 1.871767e-05 1.987555e-05 2.100983e-05 2.211954e-05 2.320496e-05 2.426694e-05 2.530651e-05 2.632479e-05 2.732290e-05 )

			(4.866288e-04 1.045561e-04 4.725135e-05 1.802188e-05 1.372837e-05 1.424527e-05 1.528219e-05 1.643602e-05 1.761952e-05 1.878658e-05 1.993683e-05 2.106522e-05 2.217022e-05 2.325176e-05 2.431046e-05 2.534723e-05 2.636308e-05 2.735905e-05 )

		);

	}

	kappa

	{

		low		(5.000e+06 75.0);

		high	(1.000e+07 925.0);

		values

		11 18

		(

			(2.303868e-01 1.716501e-01 1.037495e-01 3.273543e-02 3.567220e-02 4.129326e-02 4.843069e-02 5.659053e-02 6.554126e-02 7.510357e-02 8.514228e-02 9.555263e-02 1.062521e-01 1.171734e-01 1.282594e-01 1.394611e-01 1.507356e-01 1.620455e-01 )

			(2.307379e-01 1.722519e-01 1.051569e-01 3.442327e-02 3.616692e-02 4.172952e-02 4.876216e-02 5.686126e-02 6.577181e-02 7.530486e-02 8.532135e-02 9.571415e-02 1.063994e-01 1.173088e-01 1.283849e-01 1.395780e-01 1.508450e-01 1.621484e-01 )

			(2.310882e-01 1.728494e-01 1.064978e-01 3.642646e-02 3.686556e-02 4.218282e-02 4.910239e-02 5.713733e-02 6.600601e-02 7.550876e-02 8.550237e-02 9.587718e-02 1.065479e-01 1.174452e-01 1.285111e-01 1.396955e-01 1.509550e-01 1.622518e-01 )

			(2.314378e-01 1.734427e-01 1.077814e-01 3.879661e-02 3.760880e-02 4.265350e-02 4.945143e-02 5.741875e-02 6.624384e-02 7.571525e-02 8.568532e-02 9.604172e-02 1.066975e-01 1.175826e-01 1.286382e-01 1.398138e-01 1.510656e-01 1.623557e-01 )

			(2.317866e-01 1.740319e-01 1.090149e-01 4.156309e-02 3.839834e-02 4.314178e-02 4.980929e-02 5.770551e-02 6.648530e-02 7.592432e-02 8.587020e-02 9.620775e-02 1.068484e-01 1.177210e-01 1.287660e-01 1.399327e-01 1.511768e-01 1.624601e-01 )

			(2.321347e-01 1.746172e-01 1.102039e-01 4.470003e-02 3.923538e-02 4.364783e-02 5.017597e-02 5.799757e-02 6.673036e-02 7.613595e-02 8.605699e-02 9.637525e-02 1.070004e-01 1.178603e-01 1.288947e-01 1.400522e-01 1.512885e-01 1.625650e-01 )

			(2.324821e-01 1.751985e-01 1.113532e-01 4.809115e-02 4.012053e-02 4.417168e-02 5.055143e-02 5.829489e-02 6.697898e-02 7.635010e-02 8.624567e-02 9.654422e-02 1.071536e-01 1.180006e-01 1.290241e-01 1.401725e-01 1.514008e-01 1.626703e-01 )

			(2.328288e-01 1.757761e-01 1.124670e-01 5.152980e-02 5.722869e-02 4.471324e-02 5.093562e-02 5.859744e-02 6.723114e-02 7.656677e-02 8.643622e-02 9.671464e-02 1.073080e-01 1.181418e-01 1.291544e-01 1.402934e-01 1.515137e-01 1.627762e-01 )

			(2.331747e-01 1.763500e-01 1.135485e-01 5.479075e-02 4.203361e-02 4.527229e-02 5.132844e-02 5.890516e-02 6.748680e-02 7.678591e-02 8.662862e-02 9.688648e-02 1.074635e-01 1.182839e-01 1.292853e-01 1.404149e-01 1.516271e-01 1.628826e-01 )

			(2.335200e-01 1.769204e-01 1.146006e-01 5.773374e-02 4.305846e-02 4.584847e-02 5.172977e-02 5.921797e-02 6.774591e-02 7.700750e-02 8.682285e-02 9.705975e-02 1.076201e-01 1.184270e-01 1.294171e-01 1.405371e-01 1.517411e-01 1.629894e-01 )

			(2.338645e-01 1.774872e-01 1.156257e-01 6.033915e-02 4.412509e-02 4.644127e-02 5.213946e-02 5.953581e-02 6.800842e-02 7.723150e-02 8.701887e-02 9.723440e-02 1.077778e-01 1.185709e-01 1.295496e-01 1.406599e-01 1.518556e-01 1.630967e-01 )

		);

	}

}



equationOfState

{

	rho

	{

		low		(5.000e+06 75.0);

		high	(1.000e+07 925.0);

		values

		11 18

		(

			(4.769768e+02 4.103729e+02 3.160783e+02 5.848828e+01 4.018529e+01 3.192277e+01 2.682432e+01 2.326460e+01 2.062766e+01 1.853533e+01 1.684881e+01 1.545587e+01 1.428342e+01 1.328143e+01 1.241429e+01 1.165589e+01 1.098657e+01 1.039123e+01 )

			(4.772255e+02 4.109196e+02 3.184894e+02 6.706308e+01 4.477428e+01 3.531370e+01 2.958190e+01 2.561472e+01 2.269333e+01 2.037882e+01 1.851737e+01 1.698225e+01 1.569149e+01 1.458921e+01 1.363581e+01 1.280231e+01 1.206693e+01 1.141299e+01 )

			(4.774733e+02 4.114603e+02 3.207435e+02 7.646878e+01 4.947280e+01 3.873684e+01 3.234989e+01 2.796698e+01 2.475814e+01 2.221956e+01 2.018234e+01 1.850474e+01 1.709560e+01 1.589312e+01 1.485361e+01 1.394518e+01 1.314394e+01 1.243158e+01 )

			(4.777202e+02 4.119950e+02 3.228635e+02 8.681263e+01 5.427992e+01 4.219065e+01 3.512735e+01 3.032079e+01 2.682173e+01 2.405731e+01 2.184359e+01 2.002322e+01 1.849570e+01 1.719311e+01 1.606764e+01 1.508448e+01 1.421757e+01 1.344700e+01 )

			(4.779661e+02 4.125239e+02 3.248671e+02 9.817165e+01 5.919342e+01 4.567335e+01 3.791325e+01 3.267556e+01 2.888373e+01 2.589185e+01 2.350095e+01 2.153762e+01 1.989171e+01 1.848915e+01 1.727790e+01 1.622020e+01 1.528782e+01 1.445923e+01 )

			(4.782111e+02 4.130472e+02 3.267686e+02 1.105394e+02 6.420969e+01 4.918296e+01 4.070653e+01 3.503070e+01 3.094377e+01 2.772294e+01 2.515429e+01 2.304783e+01 2.128359e+01 1.978119e+01 1.848436e+01 1.735233e+01 1.635469e+01 1.546829e+01 )

			(4.784552e+02 4.135650e+02 3.285797e+02 1.237496e+02 6.932350e+01 5.271729e+01 4.350611e+01 3.738559e+01 3.300146e+01 2.955035e+01 2.680345e+01 2.455376e+01 2.267126e+01 2.106919e+01 1.968699e+01 1.848084e+01 1.741817e+01 1.647416e+01 )

			(4.786984e+02 4.140775e+02 3.303102e+02 1.374162e+02 7.452788e+01 5.627395e+01 4.631083e+01 3.973962e+01 3.505643e+01 3.137384e+01 2.844829e+01 2.605532e+01 2.405468e+01 2.235313e+01 2.088577e+01 1.960574e+01 1.847825e+01 1.747685e+01 )

			(4.789408e+02 4.145848e+02 3.319681e+02 1.509713e+02 7.981396e+01 5.985034e+01 4.911954e+01 4.209218e+01 3.710827e+01 3.319319e+01 3.008867e+01 2.755241e+01 2.543378e+01 2.363297e+01 2.208068e+01 2.072700e+01 1.953493e+01 1.847635e+01 )

			(4.791822e+02 4.150870e+02 3.335603e+02 1.638402e+02 8.517098e+01 6.344365e+01 5.193102e+01 4.444263e+01 3.915661e+01 3.500816e+01 3.172445e+01 2.904496e+01 2.680851e+01 2.490866e+01 2.327169e+01 2.184462e+01 2.058820e+01 1.947266e+01 )

			(4.794228e+02 4.155842e+02 3.350927e+02 1.756357e+02 9.058626e+01 6.705092e+01 5.474404e+01 4.679036e+01 4.120105e+01 3.681853e+01 3.335549e+01 3.053287e+01 2.817881e+01 2.618018e+01 2.445880e+01 2.295858e+01 2.163805e+01 2.046579e+01 )

		);

	}

}
sjohn2 likes this.
jcoelho5 is offline   Reply With Quote

Old   February 8, 2022, 02:43
Default
  #9
New Member
 
lee wei
Join Date: Dec 2020
Posts: 11
Rep Power: 5
KKlee is on a distinguished road
Hi, Joao Coelho,
Thanks for your post, but I want to know what
Quote:
value 11 18
mean, and how can I decide the two values.
KKlee is offline   Reply With Quote

Old   February 8, 2022, 07:48
Default
  #10
New Member
 
Joao Coelho
Join Date: Jun 2021
Posts: 23
Rep Power: 4
jcoelho5 is on a distinguished road
Quote:
Originally Posted by KKlee View Post
Hi, Joao Coelho,
Thanks for your post, but I want to know what mean, and how can I decide the two values.
It's the number of "rows" and "columns" of the table.
For simplicity, let's assume it was "values 3 2"


----- | T1 | T2 |
-----------------------
p1 | v(1,1) | v(1,2) |
-----------------------
p2 | v(2,1) | v(2,2) |
-----------------------
p3 | v(3,1) | v(3,2) |
-----------------------

v(1,1) is the property (for example Cp) @ pressure p1 and Temperature T1
...
v(3,2) is the property (for example Cp) @ pressure p3 and Temperature T2

So, you have 3 values of P, varying from Phigh and Plow (100bar and 50bar in the example I gave --- would result in p1=50, p2=75, p3=100 bar using the 3x2 table), and 2 values of T between Thigh and Tlow (925K and 75K in the example I gave; would result in T1=75K and T2=925K).


Basically, you decide how "refined" you want your table. It depends on your particular case.
jcoelho5 is offline   Reply With Quote

Old   February 8, 2022, 08:05
Default
  #11
New Member
 
lee wei
Join Date: Dec 2020
Posts: 11
Rep Power: 5
KKlee is on a distinguished road
Hi, Joao,

I got it, your reply is quite clear and really helpful, thanks.

lee
KKlee is offline   Reply With Quote

Old   February 8, 2022, 09:22
Default
  #12
New Member
 
lee wei
Join Date: Dec 2020
Posts: 11
Rep Power: 5
KKlee is on a distinguished road
Hi, Joao,
Sorry to bother, but I have another related question, you mentioned in the former post that you had used "rhoTabulated" with chtMultiRegionFoam, actually I want to use the "icoTabulated" with chtMultiRegionFoam, in which rho is only dependent with temperature. I want to know have you try the "icoTabulated" method as well? my simulation stucked by the error below when I set icoTabulated method into the thermalProperties file.
Quote:
--> FOAM FATAL ERROR:
Attempt to return primitive entry ITstream : /home/cunchang/OpenFOAM/cunchang-9/run/transient/constant/coolant/thermophysicalProperties/mixture/equationOfState/rho, line 40, IOstream: Version 2.0, format ASCII, line 0, OPENED, GOOD
primitiveEntry 'rho' comprises
on line 40 the punctuation token '('
on line 40 the punctuation token '('
on line 40 the label 300
on line 40 the doubleScalar 3.2
on line 40 the punctuation token ')'
on line 40 the punctuation token '('
on line 40 the label 400
on line 40 the doubleScalar 2.3
on line 40 the punctuation token ')'
on line 40 the punctuation token '(' ...
as a sub-dictionary

From function virtual const Foam::dictionary& Foam:rimitiveEntry::dict() const
in file db/dictionary/primitiveEntry/primitiveEntry.C at line 190.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam:rimitiveEntry::dict() const at primitiveEntry.C:?
#3 Foam::icoTabulated<Foam::specie>::icoTabulated(Foa m::dictionary const&) at ??:?
#4 Foam::hTabulatedThermo<Foam::icoTabulated<Foam::sp ecie> >::hTabulatedThermo(Foam::dictionary const&) at ??:?
#5 Foam::tabulatedTransport<Foam::species::thermo<Foa m::hTabulatedThermo<Foam::icoTabulated<Foam::speci e> >, Foam::sensibleInternalEnergy> >::tabulatedTransport(Foam::dictionary const&) at ??:?
#6 Foam::heThermo<Foam::rhoReactionThermo::composite, Foam::SpecieMixture<Foam::singleComponentMixture<F oam::tabulatedTransport<Foam::species::thermo<Foam ::hTabulatedThermo<Foam::icoTabulated<Foam::specie > >, Foam::sensibleInternalEnergy> > > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#7 Foam::heRhoThermo<Foam::rhoReactionThermo::composi te, Foam::SpecieMixture<Foam::singleComponentMixture<F oam::tabulatedTransport<Foam::species::thermo<Foam ::hTabulatedThermo<Foam::icoTabulated<Foam::specie > >, Foam::sensibleInternalEnergy> > > > >::heRhoThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#8 Foam::fluidReactionThermo::addfvMeshConstructorToT able<Foam::heRhoThermo<Foam::rhoReactionThermo::co mposite, Foam::SpecieMixture<Foam::singleComponentMixture<F oam::tabulatedTransport<Foam::species::thermo<Foam ::hTabulatedThermo<Foam::icoTabulated<Foam::specie > >, Foam::sensibleInternalEnergy> > > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#9 Foam::autoPtr<Foam::fluidReactionThermo> Foam::basicThermo::New<Foam::fluidReactionThermo>( Foam::fvMesh const&, Foam::word const&) at ??:?
#10 Foam::fluidReactionThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#11 ? in "/home/cunchang/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? in "/home/cunchang/OpenFOAM/OpenFOAM-9/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
Here is my thermalProperties file:
Quote:
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport tabulated;
thermo hTabulated;
equationOfState icoTabulated;
specie specie;
energy sensibleInternalEnergy;
}
dpdt off;

mixture
{
// coefficients for Hydrogen

specie
{
molWeight 2;
}
equationOfState
{
rho ((300 3.2) (400 2.3) (500 1.89) (600 1.57) (700 1.36) (800 1.19) (900 1.06));
}

thermodynamics
{
Hf 0;
Sf 0;
Cp ((300 14900)(400 14657)(500 14463)(600 14361)(700 14436)(800 14622)(900 14816));
}
transport
{
mu ((300 8.8e-6)(400 9.4e-6)(500 1.05e-5)(600 1.17e-5)(700 1.27e-5)(800 1.39e-5)(900 1.49e-5));
kappa ((300 0.16)(400 0.186)(500 0.212)(600 0.245)(700 0.284)(800 0.33)(900 0.37));
}
}
I don't know whether I set it incorrectly or chtmultiregionfoam doesn't support this setting. Please don't hesitate to comment if you have any idea, thanks.
Lee
KKlee is offline   Reply With Quote

Old   June 14, 2022, 04:37
Default
  #13
New Member
 
Joao Coelho
Join Date: Jun 2021
Posts: 23
Rep Power: 4
jcoelho5 is on a distinguished road
Quote:
Originally Posted by KKlee View Post
Hi, Joao,
Sorry to bother, but I have another related question, you mentioned in the former post that you had used "rhoTabulated" with chtMultiRegionFoam, actually I want to use the "icoTabulated" with chtMultiRegionFoam, in which rho is only dependent with temperature. I want to know have you try the "icoTabulated" method as well? my simulation stucked by the error below when I set icoTabulated method into the thermalProperties file.

Here is my thermalProperties file:

I don't know whether I set it incorrectly or chtmultiregionfoam doesn't support this setting. Please don't hesitate to comment if you have any idea, thanks.
Lee
Sorry for the late reply KKlee, I haven't checked this post in a while. I dont know if I ever used icoTabulated. but should be similar (with only 1 column table, I believe). However, you may provide the same properties through rhoTabulated anyways: you just have to use a low and high value of Plow and Phigh, and use the same properties on those two. The interpolation will result in rho varying only with temperature since it would be the same for different pressures.
Hope it helps.
jcoelho5 is offline   Reply With Quote

Old   October 25, 2022, 15:56
Default
  #14
New Member
 
Thomas M
Join Date: Aug 2018
Posts: 20
Rep Power: 7
tmik is on a distinguished road
Also looking for something like chtMultiRegionSonicFoam.

Sonic conjugate heat transfer model
tmik is offline   Reply With Quote

Old   March 17, 2023, 07:53
Default
  #15
New Member
 
Sreehari Perumanath
Join Date: Jun 2022
Posts: 28
Rep Power: 3
sreehahaha is on a distinguished road
Quote:
Originally Posted by jcoelho5 View Post
thermoType
{
type heRhoThermo;
mixture pureMixture;
transport tabulated;
thermo hTabulated;
equationOfState rhoTabulated;

specie specie;
energy sensibleEnthalpy;
}

mixture
{
specie
{
molWeight 16.22239;
}
#include "thermoProp_tables_FUEL"
}




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


in another file (constant/thermoProp_tables_FUEL), I have the following:


(changing linearly from 50 bar to 100 bar; 75K to 925K)

HTML Code:
thermodynamics

{

	Hf	0;

	Sf	0;

	Hs

	{...
	}

	Cp

	{...
	}

	Cv

	{...
	}

}



transport

{

	mu

	{...
	}

	kappa

	{...
	}

}



equationOfState

{

	rho

	{ ...
	}

}
Hi Joao Coelho,

Could you please tell me how you evaluated Hs for your system? Also, do you know how to evaluate Es of a system?

I have a system for which I know Cv as a function of p,T, and I'd like to use eTabulated for my "thermo". But I don't know how to specify Es(p,T) within "thermodynamics". Please share if you have any insights.

Many thanks,
SP
sreehahaha is offline   Reply With Quote

Reply


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
conjugate heat transfer in OpenFOAM skuznet OpenFOAM Running, Solving & CFD 99 March 16, 2017 05:07
Quenching simulation : how to set up a conjugate heat transfer between solid&liquid Rockda FLUENT 24 August 30, 2016 06:33
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
high speed gas flow the surface of powder bed? lily CFX 3 November 12, 2005 02:17
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12


All times are GMT -4. The time now is 10:12.