View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0000849 | OpenFOAM | Bug | public | 2013-05-15 03:00 | 2013-09-12 21:07 |
| Reporter | Assigned To | henry | |||
| Priority | high | Severity | major | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Platform | Linux | OS | Ubuntu | OS Version | 10.04 |
| Summary | 0000849: Violation of geometric conservation law(GCL) in OpenFOAM | ||||
| Description | Just did a GCL test in OpenFOAM 2.1.1, and found that the ALE implementation in OF fails to conserve a uniform solution field. The attached pictures show my result. As can be seen, after some deformation, the originally uniform velocity field got distorted and no longer uniform and has none zero vertical component, which is not right. Refine the mesh will only make the situation worse! Has anyone here studied this before? It's been confirmed that this violation of GCL can introduce accuracy as well as stability problems in ALE solvers. Any comments please? | ||||
| Steps To Reproduce | extract the source code and case files, recompile and run. | ||||
| Tags | No tags attached. | ||||
|
2013-05-15 03:00
|
|
|
|
Thanks for the bug-report. This is an issue with ddtPhiCorr rather than a limitation of the ALE implementation in OpenFOAM. If you comment out the ddtPhiCorr line in surfaceScalarField phiHbyA ( "phiHbyA", (fvc::interpolate(HbyA) & mesh.Sf()) + fvc::ddtPhiCorr(rAU, U, phiAbs) ); your case runs correctly producing a uniform field. The issue with the ddtPhiCorr relates to the way phiAbs is maintained on a mesh in which the faces rotate and change size. I will work on a general solution to this problem. |
|
|
This is now resolved and the updates to ddtPhiCorr will be included in the next release. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2013-05-15 03:00 |
|
New Issue | |
| 2013-05-15 03:00 |
|
File Added: testGCL.tar.gz | |
| 2013-09-02 15:44 | henry | Note Added: 0002449 | |
| 2013-09-12 21:07 | henry | Note Added: 0002488 | |
| 2013-09-12 21:07 | henry | Status | new => resolved |
| 2013-09-12 21:07 | henry | Resolution | open => fixed |
| 2013-09-12 21:07 | henry | Assigned To | => henry |
