dear all!!
i would like to
dear all!!
i would like to implement the P1 radiation model (used in the buoyantSimplRadiationFoam solver) into the chtMultiRegionFoam solver. i m using the 1.5 version and started first by implementing the P1model into the transient buoyantFoam solver which i think works well as a parameter study with respect to the absorption/emission coefficients showed plausible results. Now i m about to realize the implementation into the chtMultiRegionFoam solver, but the same parameter study indicates that the energy balance is decoupled from the radiative transfer equation (plausible behavior of the radiative flux but no changes in temperature distribution). my implementation looks like this: 1.) create a radiationModel "rad" in creatFluidFields.H with the temperature field thermof[i].T(). 2.) the object "rad" is given to the function solveEnthalpyEquation, done at the top of hEqn.H and solveEnthalpyEquation.C. 3.) add the radiation source term to the enthalpy equation in solveEnthalpyEquation.C in the form rad.Sh(thermo). i ll keep working on it and report, and would greatly appreciate any comments!! thx in advance!! aram 
hi!!
now i face another fun
hi!!
now i face another fundamental problem with this new solver: i simulated a 10x5x10 m3 room with a small solid heater (2x05x2 m3) in the middle of the room. when holding the heater temperature constant over the simulation time, results seem to be ok; but at the moment i prescribe the conjugated heat transfer BC (solidWallHeatFluxTemperatureCoupled/solidWallTemperatureCoupled) at the solid/fluid interface i get surface temperatures higher than the initial heater temperature. as a result i had a closer look at the solidWallHeatFluxTemperatureCoupled/solidWallTemperatureCoupled BC and have two remarks: 1) why is the flux not divided by the patch surface for the calculation of the gradient, as the flux, calculated in solidWallTemperatureCoupledFvPatchScalarField.C is determined in the units of Watt? _solidWallHeatFluxTemperatureCoupledFvPatchScalarF ield.C: gradient() = refCast<const>(neighbourField).flux()/K; _solidWallTemperatureCoupledFvPatchScalarField.C: Foam::tmp<foam::scalarfield> Foam::solidWallTemperatureCoupledFvPatchScalarFiel d::flux() const { const fvPatchScalarField& Kw = patch().lookupPatchField<volscalarfield,>(KName_); const fvPatchScalarField& Tw = *this; return Tw.snGrad()*patch().magSf()*Kw; } 2) how and what kind of information about the Tfield is given from solidWallHeatFluxTemperatureCoupled to solidWallTemperatureCoupled? thx for any comments!! aram 
1) You are better off starting
1) You are better off starting from 1.5.x. It fixes the flux issue and introduces a mixed boundary conditions which makes setup much easier. See the tutorial case.
2) from memory: K and neighbour internal field. 
1) oh yeah, actually i m using
1) oh yeah, actually i m using 1.5.x (updated with git from 1.5). sorry for the confusion!
2)ok, but when i m looking into solidWallTemperatureCoupledFvPatchScalarField.C i can basically find the member functions flux() and updateCoeffs() including the operation operator==(neighbourField); but what is this operator doing or where is K and neighbour internal field given over? i searched the codes from which this BC type is derived and was not able to figure out what operator==() means. 3) additionally i checked out again the chtMultiRegionFoam and simulated a simplified case of the tutorial case (heater without "wings", no left/right solid and no inlet in the upper domain). First i used convertToMeters = 1 like the tutorial, and then convertToMeters = 10 (with additional changes in makeCellSets.setSet). when i compare the temperature plots, i see plausible results for convertToMeters = 1, but the problem of a higher heater surface temperature than the initial one (508 K vs 500 K) for convertToMeters = 10. i m about to investigate further in how the data flow at the fluidsolid interface is solved and greatly appreciate comments on that!! thx a lot! aram 
hi!
ad 3) i refined the mes
hi!
ad 3) i refined the mesh (twice) for the simplified and scaled chtMultiRegionFoam tutorial case, but still have the problem of a higher heater surface temperature than the initial one (now 502 K vs 500 K). i really have no idea what i m doing wrong here. thx in advance for any help! aram 
Hi Aram,
I would like to mo
Hi Aram,
I would like to model the heating of a solid in air. I think that I need to do what you are doing. I found on a forum the tutorial "multiRegionHeater" but I find it very complicated for a start. Would you mind sending me the simplified case that you made? Maybe would could work together from this point? I will need to add surface radiation and other terms in the solid in the future. Thanks, Jean PS: I have no idea if this can help, but did you check the stability if an explicit scheme is used? (value of the Courant number) 
hi jean!
give me your mail
hi jean!
give me your mail address so i can send you the setup. concerning the courant number, the timestep is adjusted to keep Co=0.3. aram 
Hi Aram,
That would be very

Hi Aram,
Nice meeting you here again. Actually I am doing the same job as you with OpenFOAM1.6, adding radiation model to the solver chtMultiRegionFoam. The implementation is almost the same, 1.) ... 2.) ... 3.) ... 4.) rad.correct() after h equation is solved. I used a case to calculate with this new solver which I called chtMultiRegionRadiationFoam but can't get reasonable results. Here are some problems I want to ask you, 1.) Can chtMultiRegionFoam be used to calculate turbulent flow? I found that the multiRegionHeater case is using laminar model, so I have such a wonder. 2.) If I want to calculate radiation in nonparticipating media, what kind of radiation model should I choose? P1 or fvDOM? 3.) If one of the solid region is semitransparent, what should I do to deal with this situation? Is there any relationship with the radiation model? It seems that fvDOM is more flexible than P1. I am looking forward to your answer. Thanks. Xinyuan 
add Raditon model to buoyantBoussinesqPisoFoam
Quote:
Hi Aram, I am simulating a HVAC application in OF v1.6 with buoyantBoussinesqPisoFoam. I intend to run this again, but this time around integrating the P1 radiation model into the buoyantBoussinesqPisoFoam solver. It seems that you implemented a similar process as mentioned in your previous post "implementing the P1model into the transient buoyantFoam solver ". I am fairly new to OpenFOAM use, could you please educate me on how to implement the P1 radiation model in buoyantBoussinesqPisoFoam? Thank you very much for your assistance. 
Hi!!
Sorry for the late response but it seems that I haven t got a message about your posts. @Xinyuan 1.) yes; check constant/<fluidRegion>/RASProperties and constant/<fluidRegion>/turbulenceProperties 2.) I optained results with both models, but those with fvDOM were more accurate 3.) I would try to use a fluidRegion and eliminate convection (g=(0 0 0)) @samulu I have no experience with buoyantBoussinesqPisoFoam but I would say the approach is the same. Study how the radiation model is implemented in buoyantSimpleRadiationFoam (creat a radiation model in creatFields.H; add the radiation source term Sh() to the energy equation hEqn.H). All the best, Aram 
1 Attachment(s)
Dear all,
I added the radiation model to chtMultiRegionFoam and modified the solidWallMixedTemperatureCoupled BC (couples solidfluid regions) to take the radiative wall heat flux into account (for now only with the fvDOM radiation model). I attached the solver (chtMultiRegionRadFoam.tar.gz) maybe its useful for you; appreciate your comments on it. Thanks in advance! Regards, Aram 
Hi Aram,
would you mind reposting your solver. The tar.gz file is invalid. Thank you, Mirko 
Hi Aram,
I am trying to set up a simple case and check the solver chtMultiregionFoam. I looked into the tutorial case and got some idea about setting up the cell sets and the region. Actually I have been working on a different topic using interFoam. But now i need to set up a simple case in chtMultiregionFOam, for demonstrating the capabilitites of chtMultiregionFoam. I need a simple case like the one you have specified in your post to understand the basic set up first. So could you please send me the set up files of the case with a heater in the middle of the room. My email I.D is kumar.kannan@uni.lu. thankyou regards K.Suresh kumar 
1 Attachment(s)
please check the following .gz file

chtMultiRegionRadFoam
Hello,
I'm using OF1.6.x and would like to try our new solver. A newbiequestion: How do I implement it? Regards Marco 
hi!
extract the files and copy the chtMultiRegionRadFoam folder to the <user>1.6.x/applications directory. go to <user>1.6.x/applications/chtMultiRegionRadFoam and execute wmake. regards, aram 
Hi,
is it just the same procedure to implement radiation in chtMultiRegionSimpleFoam? Best Regards, tH3f0rC3 
Hi,
I want to use the solver chtMultiRegionRadFoam. Does someone have a tutorial case for me so that I can study the boundary conditions and the fvshemes,... That would be very nice, because there is no tutorial of this solver. So if someone has a little case solved by his own, I would me glad to see how he changed the data (fvShemes, fvSolution, ...). Best Regards, tH3f0rC3 
All times are GMT 4. The time now is 22:37. 