CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

chtMultiRegionSimpleFoam - Volumetric Heat Source - OpenFOAMV2212 - not working

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 16, 2024, 04:12
Default chtMultiRegionSimpleFoam - Volumetric Heat Source - OpenFOAMV2212 - not working
  #1
New Member
 
Nayeem Bin Zahid
Join Date: Jun 2024
Location: Dhaka, Bangladesh
Posts: 1
Rep Power: 0
Nayeem12 is on a distinguished road
I am investigating a subchannel of Molten Salt Reactor. The design is essentially a graphite pipe, through which, heat generating molten salt is flowing.
For the heat generation, in ANSYS, I have simply added following equation in the flowing molten salt cell zone:
1.57 * 213.5e6 * sin((y*PI)/(3.96 [m])) * 1 [W m^-3]

Now for validation, I have extracted the exact mesh and imported to OpenFOAM, added exactly same boundary conditions, and for adding the heat source to the fluid, I have added an fvOptions to the fluid region with following lines:

Code:
energySource
{
    type            scalarCodedSource;
    active          true;
    name            sourceTime; 

    scalarCodedSourceCoeffs
    {
        selectionMode   cellZone;
        cellZone        fuel_region; 
        volumeMode      absolute; 
        fields          (h);          
        
        fieldNames      (h);
        name            sourceTime;
        
        codeInclude
        #{

        #};

        codeCorrect
        #{
            Pout<< "**codeCorrect**" << endl;
        #};

        codeAddSupRho
        #{
            const scalarField& V = mesh_.V();
            const scalarField& celly = mesh_.C().component(1);
            scalarField& hSource = eqn.source();

            // Initialize min and max with the first element of celly
            scalar minCelly = celly[0];
            scalar maxCelly = celly[0];

            forAll(celly, i)
            {
                hSource[i] = -(1.57*213.5e6*sin((3.14159*celly[i])/3.96))*V[i];

                // Update min and max values
                if (celly[i] < minCelly)
                {
                    minCelly = celly[i];
                }
                if (celly[i] > maxCelly)
                {
                    maxCelly = celly[i];
                }
            }

            // Print out the min and max celly values
            Pout << "Minimum celly value: " << minCelly << endl;
            Pout << "Maximum celly value: " << maxCelly << endl;

            Pout << "**codeAddSup**" << endl;
        #};

        codeConstrain
        #{
            Pout<< "**codeConstrain**" << endl;
        #};
    }
}

limitTemperature
{
    type    limitTemperature;
    active  true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        min    839;
        max    2400;
    }
}

But the result is far from what ANSYS predicted. Even the temperature profile is completely different.

To check if I have added the boundary conditions properly, I have changed the scalarCodedSource to scalarSemiImplicitSource and added a constant specific heat generation term to the fluid region with following fvOptions:

Code:
scalarSource
{
    active  true;
    type    scalarSemiImplicitSource;
    scalarSemiImplicitSourceCoeffs
    {
        volumeMode  specific;
        selectionMode   all;
        injectionRateSuSp
        {
            h   (213.5e6 0);
        }
    }
}


limitTemperature
{
    type    limitTemperature;
    active  true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        min    839;
        max    2400;
    }
}
And did the same for ANSYS, and the temperature profile came out as perfectly same.

This led me to believe I am possibly doing something wrong with the sinusoidal volumetric heat source I am using in fvOptions?
Nayeem12 is offline   Reply With Quote

Reply

Tags
chtmultiregionsimpefoam, conjugate heat transfer, openfoam 2212, volumetric heat source

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
MapFields utility for staggered grid saeed jamshidi OpenFOAM Running, Solving & CFD 17 January 23, 2024 11:36
[OpenFOAM.com] swak4foam compiling issues on a cluster saj216 OpenFOAM Installation 5 January 17, 2023 16:05
what is swap4foam ?? AB08 OpenFOAM 28 February 2, 2016 01:22
centOS 5.6 : paraFoam not working yossi OpenFOAM Installation 2 October 9, 2013 01:41
[swak4Foam] build problem swak4Foam OF 2.2.0 mcathela OpenFOAM Community Contributions 14 April 23, 2013 13:59


All times are GMT -4. The time now is 18:06.