# Implementing Radiation Model into buoyantBoussinesqPisoFoam compiling error

 Register Blogs Members List Search Today's Posts Mark Forums Read

 December 16, 2009, 09:53 #2 New Member   David Huckaby Join Date: Jul 2009 Posts: 21 Rep Power: 9 Fabian, 2) "thermo" is generally declared in "createFields.H", see buoyantPisoFoam. 4) I don't think you need to start from scratch. To build a transient solver with bouyancy, I think you could add radiation to buoyantPisoFoam or transient terms to bouyantRadiationSimpleFoam. This would avoid the errors from the hEqn.H 5) The previous workshop (Basic & Advanced training) has some good tutorial material: http://www.openfoamworkshop.org/2009/4th_Workshop/ Hope this helps. Dave

 December 17, 2009, 13:42 #3 New Member   Fabian Hampp Join Date: Dec 2009 Location: Abu Dhabi Posts: 8 Rep Power: 8 Thanks Dave, your hints helped a lot. I follow your recommendation and made bouyantRadiationSimpleFoam transient. Unfortunately I run into a different error as follows. I can simulate a few time steps but three things are concerning myself. 1) why do i calculate rho in the beginning and with the pressure, is this correct? 2) The solver does not need any iteration for the calculation of rho. Is this correct or could it be possible that I made the solver anywhere incompressible? and 3) what does this error mean? Thanks a lot in advance, Best regards Fabian !!! Error Message: Courant Number mean: 0.0148611 max: 0.138988 Time = 0.28 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.224731, Final residual = 1.45725e-09, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.899876, Final residual = 1.18918e-09, No Iterations 2 DILUPBiCG: Solving for h, Initial residual = 0.177634, Final residual = 6.19952e-09, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.48412, Final residual = 0.0440553, No Iterations 8 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 0.00124261, global = 2.08827e-05, cumulative = -3.15106e-05 DICPCG: Solving for p, Initial residual = 0.0285576, Final residual = 8.32981e-07, No Iterations 36 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors : sum local = 4.33048e-08, global = 2.76174e-09, cumulative = -3.15079e-05 DILUPBiCG: Solving for epsilon, Initial residual = 0.161125, Final residual = 1.08227e-06, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.552436, Final residual = 1.88081e-11, No Iterations 2 ExecutionTime = 7.11 s ClockTime = 7 s Courant Number mean: 0.0468625 max: 0.219 Time = 0.3 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Ux, Initial residual = 0.466304, Final residual = 5.59278e-11, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.744887, Final residual = 4.29536e-09, No Iterations 2 DILUPBiCG: Solving for h, Initial residual = 0.245595, Final residual = 7.89469e-07, No Iterations 3 #0 Foam::error:rintStack(Foam::Ostream&) in "/home/fhampp/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/fhampp/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib/libc.so.6" #3 Foam::hRhoThermo > > > >::calculate() in "/home/fhampp/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #4 Foam::hRhoThermo > > > >::correct() in "/home/fhampp/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libbasicThermophysicalModels.so" #5 Uninterpreted: ./buoyantPisoRadiationFoam #6 __libc_start_main in "/lib/libc.so.6" #7 Uninterpreted: ./buoyantPisoRadiationFoam Floating point exception

 December 18, 2009, 15:24 #4 New Member   David Huckaby Join Date: Jul 2009 Posts: 21 Rep Power: 9 Fabian, 1) I think there are two different "rho" fields. A transported "rho" calculated by the mass conservation equation and thermodynamic "rho" calculated from "p" and "T". 2) Yes, the solver does not need to iterate to solve for rho, since this is a fully explicit equation. fvc (explicit) vs. fvm (implicit) Looking through the code there are some subtle difference between between the implementation of the PISO and SIMPLE algorthms, thus it would be better to use a PISO solver and add radiation as opposed to the opposite. Have you tried decreasing the time-step at least during the initial transient ? Dave

 February 8, 2010, 09:59 How to Include Radition heat sources in BuoyantPISOFoam #5 Member   Maruthamuthu Venkatraman Join Date: Mar 2009 Location: Norway Posts: 80 Rep Power: 9 Have you succeeded in implimenting Radiation sources in BuoyantPisoFoam ? If so, then could you give me the instructions to follow the same. Thanks

 September 5, 2010, 22:28 #6 Member   Robert Ong Join Date: Aug 2010 Posts: 60 Rep Power: 8 Hi David and Fabian, I have a question which may sound kind of stupid.... I think buoyantBoussinesqPisoFoam is only meant to solve incompressible flow, if this is the case then why did you need to bother about hEqn.h and thermo stuffs? Isn't there any other way to implement incompressible flow with radiative heat transfer apart from this? Thank you for your time and attention. Robert

 October 26, 2010, 16:46 #7 New Member   Join Date: Jul 2010 Posts: 25 Rep Power: 8 Hi Fabian and Dhuckaby, Greetings, I guess from your conservation, u have included transient terms into simpleradiationFoam solver. I am also trying to add conduction and convection heat transfer into simpleradiationFoam solver. For this, i am thinking to use transport equation to the radiation solver. Is this the right way to do this? Would you please tell me how you included transient terms into radiation solver with steps to follow? Any help will be greatly appreciated, Thanks a lot,

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post padmanathan OpenFOAM 9 October 13, 2009 05:17 waynezw0618 OpenFOAM Installation 1 February 18, 2009 11:12 jango Open Source Meshers: Gmsh, Netgen, CGNS, ... 3 November 9, 2007 14:29 fanzhong Meng FLUENT 4 May 15, 2006 11:40 CFD user CFX 3 November 25, 2002 16:16

All times are GMT -4. The time now is 09:02.