View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0002617 | OpenFOAM | Contribution | public | 2017-07-17 14:25 | 2017-07-17 16:49 |
| Reporter | huberta | Assigned To | henry | ||
| Priority | normal | Severity | minor | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| Platform | GNU/Linux | OS | Ubuntu | OS Version | 15.04 |
| Product Version | dev | ||||
| Fixed in Version | dev | ||||
| Summary | 0002617: BuoyantPimpleFoam crashes when no pressure is fixed | ||||
| Description | I tried to run the hotroom tutorial with buoyantPimpleFOAM just changing all the BC to walls but using Boussinesq linear approximation to calculate the density. This case would then look similar to the hotroom tutorial available for the solver buoyantBoussinesqPimpleFOAM with BC changed to wall everywhere. However, the simulation crashes when running it with buoyantPimpleFOAM while it does not with buoyantBoussinesqPimpleFOAM | ||||
| Steps To Reproduce | A way to solve the problem is to set a reference pressure in a fixed point for the buoyantPimpleFOAM solver in the createFields.H file, the same way as it is done for the buoyantBoussinesqPimpleFOAM i.e.: label pRefCell = 0; scalar pRefValue = 0.0; setRefCell ( p, p_rgh, pimple.dict(), pRefCell, pRefValue ); if (p_rgh.needReference()) { p += dimensionedScalar ( "p", p.dimensions(), pRefValue - getRefCellValue(p, pRefCell) ); } | ||||
| Additional Information | By doing so, in the case the pressure in not fixed via a pressure inlet BC or outlet, it would help to stabilize the calculation and prevent it from crashing | ||||
| Tags | No tags attached. | ||||
|
|
Setting a reference pressure is only correct if the flow is incompressible as is the case with theBoussinesq equation of state. Have you tried running with your "fix" with the perfect gas or other compressible equations of state? |
|
|
you are right I was more focused on the incompressible case. When run with any other model (except polynomial) it does not need the "fix" suggested and the fix even makes the calculation crash in these cases. So to sum up, the solver works fine with the set-up above without any modification. It is just when dealing with polynomial properties or Boussinesq approximation the "fix" suggested is needed |
|
|
Can you provide a patch will work for any equation of state? |
|
|
Resolved by commit 40ddc7204530d7724d31e04edf660707965e4c0e |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2017-07-17 14:25 | huberta | New Issue | |
| 2017-07-17 14:36 | henry | Note Added: 0008403 | |
| 2017-07-17 15:17 | huberta | Note Added: 0008404 | |
| 2017-07-17 15:22 | henry | Note Added: 0008405 | |
| 2017-07-17 16:49 | henry | Assigned To | => henry |
| 2017-07-17 16:49 | henry | Status | new => resolved |
| 2017-07-17 16:49 | henry | Resolution | open => fixed |
| 2017-07-17 16:49 | henry | Fixed in Version | => dev |
| 2017-07-17 16:49 | henry | Note Added: 0008406 |
