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

Making Openfoam radiation solver

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 24, 2017, 03:01
Post Only radiation validation problem
  #1
Member
 
Jung da yoon
Join Date: Nov 2017
Posts: 53
Rep Power: 4
kane is on a distinguished road
HI all
I am trying to make only radiation solver, so far with little success.
Here are the codes.

1. radiation.C
Code:
\*---------------------------------------------------------------------------*/

#include "fvCFD.H"
#include "radiationModel.H"
#include "simpleControl.H"
#include "fvOptions.H"

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

int main(int argc, char *argv[])
{
    #include "setRootCase.H"
    #include "createTime.H"
    #include "createMesh.H"
    #include "createFvOptions.H"
    #include "initContinuityErrs.H"
    #include "createControl.H"


        Info<< "Reading Temperature field\n" << endl;
        volScalarField T
        (
                IOobject
                (
                        "T",
                        runTime.timeName(),
                        mesh,
                        IOobject::MUST_READ,
                        IOobject::AUTO_WRITE
                ),
        mesh
        );

        autoPtr<radiation::radiationModel> radiation
        (
                radiation::radiationModel::New(T)
        );


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

    Info<< "\nStarting time loop\n" << endl;

    while (simple.loop())
    {
        Info<< "Time = " << runTime.timeName() << nl << endl;
        {
                radiation->correct();
        }

        runTime.write();

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
            << nl << endl;
    }

    Info<< "End\n" << endl;

    return 0;
}
2. Make/files

Code:
radiationFoam.C

EXE = $(FOAM_USER_APPBIN)/radiationFoam
3. Make/options
Code:
EXE_INC = \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/sampling/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
    -I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \

EXE_LIBS = \
    -lfiniteVolume \
    -lfvOptions \
    -lsampling \
    -lmeshTools \
    -lfluidThermophysicalModels \
    -lspecie \
    -lradiationModels \
    -lmeshTools
Also I'm trying to solve a simple radiation case with radiation solver.
I have tried Fvdom radiation model and finally I want to know heat flux on walls.

I have attached the following case named Enclosure.

Problem :
1. I don't know if only radiation solver is correct. So plz tell me errors or bugs

2. I want to know heat flux on walls. As a result of paraview program, I think 'qin' means 'heat flux on walls'. Is it correct?

3. After running this case with paraview, I checked that qin of front and back walls is different from other walls. Why it happens ??

I'm very confused with only radiation solver and case.
Plz help me !!

Happly christmas
Attached Files
File Type: zip Enclosure.zip (4.9 KB, 3 views)

Last edited by kane; December 25, 2017 at 00:38. Reason: Adding information+
kane is offline   Reply With Quote

Old   February 3, 2018, 20:25
Default
  #2
Member
 
Jaydeep
Join Date: Jun 2015
Posts: 34
Rep Power: 6
jaydeep is on a distinguished road
Hi Kane,

If you want to solve radiation only -- I would recommend using any of the heat transfer solvers. For e.g. you can pick up buoyantSimpleFoam, and turn "FrozenFlow" on in fvSolution. It will stop mass momentum equations and will only solve energy.

p.s. I used the same approach with chtMultiRegionSimpleFoam
jaydeep is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Sharing links for two phase solver packages developed by openfoam community swap_9068 OpenFOAM Programming & Development 1 April 2, 2017 06:43
[swak4Foam] funkyDoCalc with OF2.3 massflow NiFl OpenFOAM Community Contributions 11 November 1, 2016 07:43
How to rewrite a standard OpenFOAM solver as a C++ class cfbaptista OpenFOAM Programming & Development 7 March 23, 2016 05:50
Frequently Asked Questions about Installing OpenFOAM wyldckat OpenFOAM Installation 0 January 1, 2014 20:21
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 18:51


All times are GMT -4. The time now is 14:47.