unicfdlab / hybridCentralSolvers

United collection of hybrid Central solvers - one-phase, two-phase and multicomponent versions
GNU General Public License v3.0
92 stars 59 forks source link

Issue with wedge domains #6

Open hinman-stf opened 5 years ago

hinman-stf commented 5 years ago

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 commented 5 years ago

Dear hinman-stf,

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

Best regards, Victoria

hinman-stf commented 5 years ago

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 notifications@github.com 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 https://github.com/unicfdlab/hybridCentralSolvers/issues/6?email_source=notifications&email_token=ANNJTZQCFJHWXX3WGUTHLRLQNXDPLA5CNFSM4I6HBVL2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEAXVVEA#issuecomment-539974288, or mute the thread https://github.com/notifications/unsubscribe-auth/ANNJTZWZVQFVLFLRZXL55ZDQNXDPLANCNFSM4I6HBVLQ .

hinman-stf commented 5 years ago

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 notifications@github.com 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 https://github.com/unicfdlab/hybridCentralSolvers/issues/6?email_source=notifications&email_token=ANNJTZQCFJHWXX3WGUTHLRLQNXDPLA5CNFSM4I6HBVL2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEAXVVEA#issuecomment-539974288, or mute the thread https://github.com/notifications/unsubscribe-auth/ANNJTZWZVQFVLFLRZXL55ZDQNXDPLANCNFSM4I6HBVLQ .

vkorchagova commented 5 years ago

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

hinman-stf commented 5 years ago

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 notifications@github.com 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 https://github.com/unicfdlab/hybridCentralSolvers/issues/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