View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0003581 | OpenFOAM | Bug | public | 2020-10-22 16:12 | 2020-11-23 11:31 |
| Reporter | vitor.geraldes@tecnico.ulisboa.pt | Assigned To | henry | ||
| Priority | normal | Severity | major | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Platform | Unix | OS | Other | OS Version | (please specify) |
| Product Version | 8 | ||||
| Fixed in Version | 8 | ||||
| Summary | 0003581: bug in externalWallHeatFluxTemperature BC | ||||
| Description | The radiation effect is not well implemented. When emissivity > 0, the temperature distribution is not well computed. In the attached testing case is based on the following case, OpenFOAM-8/tutorials/multiphase/compressibleInterFoam/laminar/depthCharge2D/ The scale of the mesh was reduced by a factor of 0.01, and the bubble was removed. The external temperature was increased to 400 K. As we can see by the simulation, before one 1 s, the internal temperature increases beyond 400 K, which is a violation of the energy balance. | ||||
| Steps To Reproduce | Run the scipt "Allrun" in the attached case, which uses the solver "compressibleInterFoam". | ||||
| Tags | No tags attached. | ||||
|
|
emissivityErrorTestCase.tar.gz (593,018 bytes) |
|
|
How do you think radiation should be implemented in externalWallHeatFluxTemperature? Can you provide a patch? |
|
|
Have you tried relaxing the radiation contribution and using more outer iterations to converge the heat transfer? |
|
|
externalWallHeatFluxTemperatureFvPatchScalarField.C (14,295 bytes) |
|
|
The point of the current implementation is to avoid the "/(Ta-Tp+VSMALL)" which can cause serious problem, particularly if (Ta-Tp) = -VSMALL. The stabilisation will have to be done more carefully than this to avoid crashes. |
|
|
If we factor (Ta^4 - Tp^4) as (Ta - Tp)((Ta^3 + Tp^3) + TaTp(Ta + Tp)) we can avoid the division by (Ta -Tp) |
|
|
Resolved in OpenFOAM-8 by commit 132a11d3674a0bd71761bedd5d3e75d52b373b82 Resolved in OpenFOAM-dev by commit bdbb2ae0a83a5aa06ed55c472ef7dac602bcb077 |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2020-10-22 16:12 | vitor.geraldes@tecnico.ulisboa.pt | New Issue | |
| 2020-10-22 16:12 | vitor.geraldes@tecnico.ulisboa.pt | File Added: emissivityErrorTestCase.tar.gz | |
| 2020-11-21 19:50 | henry | Note Added: 0011703 | |
| 2020-11-21 23:16 | henry | Note Added: 0011728 | |
| 2020-11-21 23:44 | vitor.geraldes@tecnico.ulisboa.pt | File Added: externalWallHeatFluxTemperatureFvPatchScalarField.C | |
| 2020-11-22 01:13 | henry | Note Added: 0011730 | |
| 2020-11-22 16:46 | henry | Note Added: 0011732 | |
| 2020-11-23 11:31 | henry | Assigned To | => henry |
| 2020-11-23 11:31 | henry | Status | new => resolved |
| 2020-11-23 11:31 | henry | Resolution | open => fixed |
| 2020-11-23 11:31 | henry | Fixed in Version | => 8 |
| 2020-11-23 11:31 | henry | Note Added: 0011736 |
