Issue with wedge domains

[hinman-stf]

There seems to be an issue regarding pimpleCentralFoam with wedge domains. To replicate, please run the attached case. This is a standard benchmark case for supersonic compressible solvers. The grid is uniform and is wedge type. Using rhoCentralFoam the solution is correct. Using pimpleCentralFoam, there is an odd spurious jet type phenomenon along the central axis.

The issue is independent of flux scheme (have tried VanLeer and VanAlbada) as well as local time stepping, Euler or backward time stepping, or significant changes in CFL. This indicates to me it is a bug in the flux formulation when dealing with the wedge (empty) central axis boundary.
CDNozzleValidation.tar.gz

I am running this suite of solvers with OpenFOAM 6.

[vkorchagova]

Dear hinman-stf,

blockMeshDict seems uncorrect. When I run the blockMesh command in your case, I see topology errors.

Best regards, Victoria

[hinman-stf]

The mesh was made using pointwise. You do not need to rerun the mesh. Just the solver.

…

On Wed., Oct. 9, 2019, 6:15 a.m. vkorchagova ***@***.***> wrote: Dear hinman-stf, blockMeshDict seems uncorrect. When I run the blockMesh command in your case, I see topology errors. Best regards, Victoria — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#6?email_source=notifications&email_token=ANNJTZQCFJHWXX3WGUTHLRLQNXDPLA5CNFSM4I6HBVL2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEAXVVEA#issuecomment-539974288>, or mute the thread <https://github.com/notifications/unsubscribe-auth/ANNJTZWZVQFVLFLRZXL55ZDQNXDPLANCNFSM4I6HBVLQ> .

[hinman-stf]

This is the same case: https://curiosityfluids.com/2016/07/14/converging-diverging-nozzle-v2-rhocentralfoam/ On Wed., Oct. 9, 2019, 10:10 p.m. Schuyler Hinman <hinmanws@gmail.com> wrote:

…

The mesh was made using pointwise. You do not need to rerun the mesh. Just the solver. On Wed., Oct. 9, 2019, 6:15 a.m. vkorchagova ***@***.***> wrote: > Dear hinman-stf, > > blockMeshDict seems uncorrect. When I run the blockMesh command in your > case, I see topology errors. > > Best regards, Victoria > > — > You are receiving this because you authored the thread. > Reply to this email directly, view it on GitHub > <#6?email_source=notifications&email_token=ANNJTZQCFJHWXX3WGUTHLRLQNXDPLA5CNFSM4I6HBVL2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEAXVVEA#issuecomment-539974288>, > or mute the thread > <https://github.com/notifications/unsubscribe-auth/ANNJTZWZVQFVLFLRZXL55ZDQNXDPLANCNFSM4I6HBVLQ> > . >

[vkorchagova]

Yes, I see. Good.
Next, one more correction: please change "sensibleInternalEnergy" to "sensibleEnthalpy" in the thermophysicalProperties dictionary.

[hinman-stf]

Yes. Sorry. The case was originally set up for rhoCentralFoam. I quickly swapped it over to test the wedge domain issue brought to my attention by a student. This should not make a difference for the issue.

…

On Thu, Oct 10, 2019 at 6:35 AM vkorchagova ***@***.***> wrote: Yes, I see. Good. Next, one more correction: please change "sensibleInternalEnergy" to "sensibleEnthalpy" in the thermophysicalProperties dictionary. — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#6?email_source=notifications&email_token=ANNJTZUYTMH4Q7CSWUJAAF3QN4OQHA5CNFSM4I6HBVL2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEA4CUFI#issuecomment-540551701>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/ANNJTZQWYVNAXIZ2E22CXNDQN4OQHANCNFSM4I6HBVLQ> .

-- *Schuyler Hinman, PhD* Director Specialized Thermo Fluids Solutions Ltd. Ph: 403-829-2796 E: wshinman@stfsol.com W: www.stfsol.com