uniform heating of fluid region
Hi,
I am using bousinesqBuoyantSimpleFoam to model a heat exchager. I would like to give uniform heating to whole fluid region as a source term (volumetric heating in W/m^3). Can anybody suggest me how to do this in OpenFOAM. Regards Santhosh. 
Doubt regd TEqn.H
I am looking at the TEqn.H in boussinesqBuoyantSimepleFoam. following is snippet copied from the above file.
Code:
fvScalarMatrix TEqn Code:
fvm::Sp(fvc::div(phi), T) Code:
div(sV) = sdiv(V)+V.grad(s) Combination of first two terms will results in the following, Code:
phi.gradT Regards Santhosh. 
I have added the source term in the TEqn.H as below.
Code:
fvScalarMatrix TEqn I will post back if there is any other alternative. Regards Santhosh 
That sounds interesting.
Please let us know, if you have some results. Regards. 
Hi Oli,
That seems to be working, I am solving very big heat exchanged problem, As for as residuals and priliminary results I am not facing any major problem (Compared to the one without source terms). The only thing I am worried about is Eugene's following post.. http://www.cfdonline.com/Forums/ope...newpost.html I am also finding many problems with these set of solvers, especially convergence of Temperature is very very slow. forget about second order results (Remember I am a novice user still, and I think I am continue to be novice as long as openfoam does not provide documentation) I am doing lot of parametric study (Thanks to my organization for proving me the infinite compute power on world's faster super computer, EKA) Regards Santhosh 
Sensitivity of Turbulent prandtl number
1 Attachment(s)
Hi,
Please look at the residual file attached here, Unknowingly I initially used the turbulent Prandtl number (Prt) equals 5. After reading the documentation I came to know Prt is used as a contribution to production term in turbulent kinetic energy equation. Also for ke model it has to be aroung 0.85. So I changed the value of Prt to 0.8 from 5. I observed lot of variation in residue plots. In the attached plot, sensitivity with viscosity (which I changed matching experimental) and Prt. The results are similar in case of variation of viscosity (although very little change). But with variation of Prt I found completely different residual. I am using upwind numerical schemes for interpolation. Please can anybody through a light on the use of Prt. Is it okey to use Prt=5 if not can you suggest chages to get residual down. I have observed simular variation in residual for other parameteers also. Thanks Santhosh 
Hi Santhos,
I am new on heat transfer problems using cfd and I would like to model a constant volume heat generation, following the approach you proposed, since it seems to me that there is no already implemented function to model it in OF. Before starting, I would like to know it you have some results showing that this is a good approach, i.e. cfd matches with experimental results. regards, mad 
Hi,
I have added source term as I have mentioned above. I have even completed one test case with the modified solver. But unfortunately, I have not compared the results with the any benchmark case. Qualitatively It seemed OK to me. If you find any benchmark cases to test it, please let me know, and I am happy to test it for you. Regards, Santhosh 
Hi Santhosh,
finally I included the heat source contribution on chtMRFoam, as you suggested above. However, as explained here, the new solver gives acceptable values only if I use a single region. In the case of multiple region, of which only one is heated, the solver does not perform well. I am thinking that this may be due to the coupling between the regions. Have you any experience in that? Any suggestion is welcome. Regards, mad 
adding term in solidWallMixedTemperatureCoupled
Hi all,
trying to find an explanation on the strange results of my cht simulation including heat source, I ended up with the following: solidWallMixedTemperatureCoupled in OF 1.6.x and 1.7 implements a coupling that differs from OF 1.6 and previous. In OF 1.6 (code): Code:
forAll(*this, i) Code:
this>refValue() = nbrIntFld; My doubts raise when applying a heat source. The 1D steady state conduction equation + heat source says that: dT2/dx2 + H/K = 0. Integrating once gives dT/dx  H/K*x = 0. It means that, in the case of an heat source refGrad is not equal to zero, but it is proportional to the heat souce itself. For this reason, I am thinking to modify the solidWallMixedTemperatureCoupled to include this contribution as well. K is already read into solidWallMixedTemperatureCoupled.C, and I can easily modify the solver to read H as well. What I am missing is how to get the size of the cell normally to the coupling interface. Is there anyone that can confirm my approach? anyone that can help me to understand how to get the cell size? thanks in advance from any suggestion I may get. cheers, mad 
Hey Santhosh,
Quote:
I think I am now able to answer the question. Basically we have the term d/dxj (rho * u_j * T) which includes the density. This term goes to d/dxj (rho * u_j * T) = rho * d/d_xj ( u_j * T) + u_j* T * d rho /dxj The first term equals to HTML Code:
fvm::div(phi, T) From the continuity equation we know d/dxj (rho * uj)=0 According to product rule d/dxj (rho * uj ) = uj * d rho/dxj + rho * duj/dxj =0 and therefore uj * d rho / dxj =  rho * duj/dxj Therefore the second term becomes u_j* T * d rho /dxj =  T * rho * duj/dxj Which is in OF the term HTML Code:
fvm::SuSp(fvc::div(phi), T) The additional term is zero for incompressible flows. 
Thanks for your effort in neatly explaining. Actually a while ago, My professor cleared my doubt about post. Sorry I forgot to post it back to forum. The explanation was similar to the one you explained.
Thanks again Santhosh 
some error
Quote:
I used user idea for building the source term in my solver and i specify Q as volumetric heat source dimension my solver is compiling properly but when i am using the solver it is giving me some kind of dimension error i am not able to debug it if you can help. Code:
Different dimensions for = 
Can any of you tell me How to implement the volumetric source term in 'buoyantSimpleFoam'?

Hello,
I am reviving an old thread but I was wondering if anyone has any insights on this one. I am trying to do exactly what the above posters intended, adding a volumetric heat source term to buoyantboussinesq solver. Unfortunately, adding the source term in the Teqn file does not work, as the solver does not recognize the Q term (I also had to change rho to rhok, hope that is correct). I realize that the field Q is not recognized, do I also need to add it in the createFields dictionary? Forgive my ignorance, new at adjusting solvers :( Thank you so much in advance. Kind regards, Theodore.  Hello, Reviving an old thread in the hope of people still being around. i'm trying to add the heat source term in the same solver without much success. Anyone has a more detailed approach as to how this was implemented (i.e. the Q/(rho*cp) term). Thanks in advance. Theodore. 
Dear Theodore,
you need to declare it in createFields.H. You should do it accordingly to the other fields which are declared in there. It would be easier, though, to understand your problem if you could post your code and the error message. Kind Regards Anne 
All times are GMT 4. The time now is 10:58. 