View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0001697 | OpenFOAM | Bug | public | 2015-05-18 20:05 | 2015-05-25 10:35 |
| Reporter | Assigned To | henry | |||
| Priority | normal | Severity | minor | Reproducibility | always |
| Status | closed | Resolution | no change required | ||
| Platform | GNU/Linux | OS | Ubuntu | OS Version | 14.04 |
| Summary | 0001697: Instability of transported quantity in scalarTransportFoam | ||||
| Description | Using the solver scalarTransportFoam, the use of a diffusion coefficient lower than 0.01 gives temperatures out of bounds (lower than 0 or bigger than 1). | ||||
| Steps To Reproduce | Use a diffusion coefficient lower than 0.001 in the tutorial test case: tutorials/basic/scalarTransportFoam/pitzDaily/ | ||||
| Tags | No tags attached. | ||||
|
|
Hello! div(phi,T) bounded Gauss linearUpwind grad(T);//Gauss linearUpwind grad(T); I tried with this, it makes it bounded perfectly between 0 and 1. Dt=0.001 If u use this: div(phi,T) Gauss linearUpwind grad(T); The upper limit T is 1.05 around. |
|
|
If you need the transported property to be rigorously bounded you will need to use a bounded scheme. NVD or TVD schemes guarantee boundedness in 1D but not in 2 or 3D although they provide better boundedness than unlimited higher-order schemes. If you use linearUpwind and want to get a more bounded solution you will need to limit the gradient. If these options still do not provided the level of boundedness there are two options left: upwind or re-implement the transport equation using MULES. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2015-05-18 20:05 |
|
New Issue | |
| 2015-05-25 10:09 | sharonyue | Note Added: 0004800 | |
| 2015-05-25 10:35 | henry | Note Added: 0004801 | |
| 2015-05-25 10:35 | henry | Status | new => closed |
| 2015-05-25 10:35 | henry | Assigned To | => henry |
| 2015-05-25 10:35 | henry | Resolution | open => no change required |
