View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0001989 | OpenFOAM | Bug | public | 2016-02-04 21:49 | 2016-02-09 17:03 |
| Reporter | lifeinhand | Assigned To | henry | ||
| Priority | normal | Severity | major | Reproducibility | always |
| Status | closed | Resolution | suspended | ||
| Platform | GNU/Linux | OS | Ubuntu | OS Version | 12.04 |
| Summary | 0001989: createTurbulenceFields or R utility not compatible with timeVaryingMappedFixedValueFvPatchField | ||||
| Description | simulation completed with RAS turbulence model ( have tested kEpsilon, realisableKE), with inlet boundary condition, timeVaryingMappedFixedValueFvPatchField, and specified (U, k and epsilon). When using R utility to calculate Reynolds Stress, it will fail. | ||||
| Steps To Reproduce | take Tutorial pitzDailyExpInlet run cd tutorial/incompressible/simpleFoam/pitzDailyExpInlet blockMesh simpleFoam # after solved, use either createTurbulenceFields or R it will give error | ||||
| Additional Information | Writing turbulence field omega --> FOAM FATAL ERROR: object of type N4Foam9DataEntryIdEE is not allocated From function Foam::autoPtr<T>::operator->() in file /home/openfoam/OpenFOAM/OpenFOAM-2.4.0/src/OpenFOAM/lnInclude/autoPtrI.H at line 176. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::timeVaryingMappedFixedValueFvPatchField<double>::write(Foam::Ostream&) const at ??:? #3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::writeEntry(Foam::word const&, Foam::Ostream&) const at ??:? #4 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::writeData(Foam::Ostream&) const at ??:? #5 Foam::regIOobject::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const at ??:? #6 ? at ??:? #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #8 ? at ??:? Aborted | ||||
| Tags | No tags attached. | ||||
|
|
There is no automatic way to convert the timeVaryingMappedFixedValue BC for k into the equivalent for R. It would be possible to simply evaluate R for the current state of all the BCs and output them all as type "calculated" which would be OK for post-processing but not for changing from say a k-epsilon to an R-epsilon turbulence model. |
|
|
Could you clarify if you need R for post-processing purposes only or are you running k-epsilon as a start-up simulation for running Reynolds stress? |
|
|
Thank you for looking into this issue.The createTurbulenceFields and R utility, I only used for post-processing purpose. The method you indicated earlier works fine for me, for now. I have not tested running k-epsilon as a start-up simulation for running Reynolds stress yet. |
|
|
OK, I think it would be good if the R utility (and the R option of createTurbulenceFields) had the option to generate calculated BC for R for post-processing but this is not high on my priority list. It would be good if someone would either fund this work or contribute a patch for this functionality. In the meantime I will close this report as you are able to work around the issue. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2016-02-04 21:49 | lifeinhand | New Issue | |
| 2016-02-04 22:05 | henry | Note Added: 0005923 | |
| 2016-02-09 13:38 | henry | Note Added: 0005927 | |
| 2016-02-09 16:31 | lifeinhand | Note Added: 0005928 | |
| 2016-02-09 17:03 | henry | Note Added: 0005929 | |
| 2016-02-09 17:03 | henry | Status | new => closed |
| 2016-02-09 17:03 | henry | Assigned To | => henry |
| 2016-02-09 17:03 | henry | Resolution | open => suspended |
