CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   OpenFOAM 6, chtMultiRegionFoam residualcontrol for steady state not working (https://www.cfd-online.com/Forums/openfoam-solving/206032-openfoam-6-chtmultiregionfoam-residualcontrol-steady-state-not-working.html)

qtruong August 28, 2018 10:45

OpenFOAM 6, chtMultiRegionFoam residualcontrol for steady state not working
 
Hi,

Since chtMultiRegionSimpleFoam was removed in OpenFOAM 6, I tried running a steady state case with chtMultiRegionFoam, with a residualcontrol (shown below). But it does not stop when converged. Has anyone tested if residualcontrol works in the new chtMultiRegionFoam?

Code:

PIMPLE
{
    residualcontrol
    {
          U    1e-4;
    }
}


Robin.Kamenicky September 2, 2018 07:46

Hi qtruong,

could you try
Code:

residualControl
instead of
Code:

residualcontrol
Kind regards,
Robin

qtruong September 3, 2018 11:48

Quote:

Originally Posted by Robin.Kamenicky (Post 705011)
Hi qtruong,

could you try
Code:

residualControl
instead of
Code:

residualcontrol
Kind regards,
Robin

thank you for the response.

Sorry, it was my typo while writing the post. It was actually "residualControl" in the code, and it's confirmed by the log file:

Code:

PIMPLE: Region gas
PIMPLE: Convergence criteria found
            U:    tolerance  0.0001

Otherwise, with lower "c", there would be "no convergence criteria found".

Nonetheless, the convergence is reached, but chtMultiRegionFoam keeps running until my endTime.

Robin.Kamenicky September 3, 2018 14:04

Hi qtruong,

Quote:

Originally Posted by qtruong (Post 705165)
Nonetheless, the convergence is reached, but chtMultiRegionFoam keeps running until my endTime.

the chtMultiRegionFoam should run till the endTime. It is transient solver not steady-state as chtMultiRegionSimpleFoam.

The residualControl allows to exit PIMPLE outer loops when the residual criteria are fulfilled. Then it continues with next time step.

Kind regards,
Robin

qtruong September 4, 2018 12:47

Quote:

Originally Posted by Robin.Kamenicky (Post 705176)
Hi qtruong,


the chtMultiRegionFoam should run till the endTime. It is transient solver not steady-state as chtMultiRegionSimpleFoam.

The residualControl allows to exit PIMPLE outer loops when the residual criteria are fulfilled. Then it continues with next time step.

Kind regards,
Robin

Is there then a way to solve steady state conjugate heat transfer in OpenFoam 6 ?

Robin.Kamenicky September 4, 2018 13:12

Hi qtruong,
Quote:

Originally Posted by qtruong (Post 705288)
Is there then a way to solve steady state conjugate heat transfer in OpenFoam 6 ?

There is no solver for steady state conjugated heat transfer in OpenFoam6. However, if your problem is steady state problem, the transient solver will converge to steady state.

Eventually you can try to use local time stepping (localEuler time derivative approximation).

Another possibility is to use OpenFoam5 and chtMultiRegionSimpleFoam. If I understand right, the difference between OpenFoam5 and OpenFoam6 is only in usability for conjugated heat transfer solvers.

Hope it helps,
Robin

julieng March 2, 2019 10:37

Hello,


I try to result a CHT steady state case with openfoam 6 and chtMultiRegionFoam. I am not able to solve it, it crashes after few time steps.

chtMultiRegionSimpleFoam done the job with Openfoam 5, I obtained good accuracy with analytical solutions.



I read your post Robin.Kamenicky, you wrote :


Code:

There is no solver for steady state conjugated heat transfer in  OpenFoam6. However, if your problem is steady state problem, the  transient solver will converge to steady state.

What is the reason to have give up the steady state solver in openfoam 6 ?


Best regards

Robin.Kamenicky March 2, 2019 11:53

Hi julieng,

Actually, I have been mistaken. The official documentation of OF6 for chtMultiRegionFoam tells:
Code:

Solver for steady or transient fluid flow and solid heat conduction, with
    conjugate heat transfer between regions, buoyancy effects, turbulence,
    reactions and radiation modelling.

What was the error?

Kind regrads,
Robin

Robin.Kamenicky March 2, 2019 12:00

Hi julieng,

according to tutorials. You can setup scheme for time derivative to
Code:

steadyState
Kind Regards,
Robin

julieng March 4, 2019 08:15

Hello Robin,

Yes it works for stationary cases, I have test it. I have same results than openfoam v5.

I try the new functionnality
"wallHeatTransferCoeff"

I read on the Openfoam 6 release:
Data Processing: function objects for individual regions in a multi-region simulation [ commit a5a034 ]; wallHeatTransferCoeff function object to calculate the wall heat transfer coefficient [ commit 99841e ]; in wallHeatFlux, improved efficiency of heat flux calculation [ commit 6584fa ] and corrected contribution of radiative heat flux [ commit 396259 ].

But I see that it works only for incompressible fluids… cht is for compressible fluids. Maybe someone knows how to modify the file.


I see the post of bloerb
https://www.cfd-online.com/Forums/op...ing-print.html

But I Don't know how to do. I am a total beginner in function object.

Best regards

hokhay December 18, 2019 03:24

Quote:

Originally Posted by julieng (Post 726688)
Hello Robin,

Yes it works for stationary cases, I have test it. I have same results than openfoam v5.

I try the new functionnality
"wallHeatTransferCoeff"

I read on the Openfoam 6 release:
Data Processing: function objects for individual regions in a multi-region simulation [ commit a5a034 ]; wallHeatTransferCoeff function object to calculate the wall heat transfer coefficient [ commit 99841e ]; in wallHeatFlux, improved efficiency of heat flux calculation [ commit 6584fa ] and corrected contribution of radiative heat flux [ commit 396259 ].

But I see that it works only for incompressible fluids… cht is for compressible fluids. Maybe someone knows how to modify the file.


I see the post of bloerb
https://www.cfd-online.com/Forums/op...ing-print.html

But I Don't know how to do. I am a total beginner in function object.

Best regards

Do you need to change the anything in fvSolution and controlDict too? Such as relaxationFactors or Courant number.

Thank you


All times are GMT -4. The time now is 00:25.