|
[Sponsors] |
June 16, 2008, 04:30 |
Dear All,
I'm challenging a
|
#1 |
New Member
David Palko
Join Date: Mar 2009
Location: Stockholm
Posts: 17
Rep Power: 17 |
Dear All,
I'm challenging a task of performing DNS simulation of strongly heated single-phase flow in a pipe. I'm wondering if anyone has an experience with DNS simulation in a pipe geometry with heat transfer using OpenFoam? Any kind on help is highly appreciated. Regards, David |
|
June 25, 2011, 04:10 |
|
#2 | |
New Member
Join Date: Oct 2010
Posts: 24
Rep Power: 15 |
Quote:
|
||
October 21, 2017, 04:15 |
|
#3 |
Member
Alex Jarosch
Join Date: Dec 2015
Location: Austria
Posts: 32
Rep Power: 10 |
Hi,
did anybody follow up on this and performed DNS simulations of heated single phase pipe flow in OpenFOAM? Would love to hear about that. Regards, Alex |
|
July 6, 2019, 05:17 |
|
#4 |
Senior Member
Jianrui Zeng
Join Date: May 2018
Location: China
Posts: 157
Rep Power: 7 |
I am doing DNS of heated pipe flow. But I do not know how to recognize the current state of simulation. I am not sure if the simulation is going in the right way and when it will end. It costs a lot of time. Any suggestion?
|
|
July 6, 2019, 08:27 |
|
#5 |
Senior Member
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15 |
the same way you check any numerical experiment. You do some structural (statistically speaking) analysis of your results, and compare them with experiments. you have azimuthal homogeneity in pipe flows. if theres no stratification, a fully developed pipe flow behaves similar to a channel flow, for first order staristics. you can start verifying a "baseline" case with no buoyancy. Hint: Stable startification "relaminarizes" the flow.
|
|
July 6, 2019, 08:30 |
|
#6 |
Senior Member
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15 |
incidental to this discussion but... Have in mind PISO is not an algorithm suitable for DNS. it is over-dissipative in the inertial range of turbulence, even for time scales below kolmogorov's.
... Not to say you cannot use it for such purpose though. |
|
September 4, 2019, 04:03 |
|
#7 | |
Senior Member
Jianrui Zeng
Join Date: May 2018
Location: China
Posts: 157
Rep Power: 7 |
Quote:
Is there any paper or book that verifies PISO is not an algorithm suitable for DNS? Except for PISO, are there other suitable algorithms for DNS? PIMPLE? Thank you. |
||
September 5, 2019, 09:12 |
Run a DNS Re_tau=300 channels flow, driven by a constant dp/dx
|
#8 |
Senior Member
Santiago Lopez Castano
Join Date: Nov 2012
Posts: 354
Rep Power: 15 |
Hi,
The first order accuracy reduction of PISO for low (<1) CFL was acknowledged by ISSA, the proposer of PISO. You can easily check this issue by running a classic channel flow, driven by a CONSTANT pressure gradient (not by adjusting the bodyforce to maintain a mean velocity), with a higher-than-180 Re_tau. You'll see spurious wiggles in the pressure RMS profiles near the walls. Articles? First, the one of Issa (1986?). Then you can check the one by Tukovic & Jasak: "Consistent second-order time-accurate non-iterative PISO-algorithm" That skims through the inconsistency problem of PISO. Then, you can read the following: "Direct numerical simulation and wall-resolved Large Eddy Simulation in nuclear thermal hydraulics" Where they do some DNS and LES using OpenFOAM. The problem of over-dissipation in numerical codes can be tackled by running with time scales of the order of Kolmogorov's time scales, preventing such dissipation to be produced on the first place. But then there's the inconsistency problem (1st order accuracy) of PISO... As I mentioned, I'm not saying you **should not** use PISO for DNS; instead, you should prove that the solutions obtained therein are not "contaminated" by the aforementioned issues. One can simply avoid this problems either by using a code that linearizes the advection term in a consistent way for implicit codes, or project the term in time explicitely, or just use one of the many spectral/DG/FR codes that exist in the wild. (SIDE NOTE/RANT: Validating any solver/turbulence model using a channel flow with Re_tau = 180 is "cheating"; The inertial range of turbulence is almost absent, hence the flow is almost purely dissipative which, in general, masks any (unexpected) dissipative errors that the code may have. Furthermore, the turbulence models are just there for the free ride, specially for low-order (or FV) codes.) |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Heated vertical channel flow RANS | ep4 | OpenFOAM Running, Solving & CFD | 0 | November 24, 2008 09:37 |
Fluent-Flow over heated cylinder | Khyati | Main CFD Forum | 0 | July 12, 2006 00:17 |
flow in pipe | co2 | FLUENT | 5 | December 21, 2003 18:39 |
Pipe Flow | Kunal Jain | Main CFD Forum | 2 | August 27, 2001 10:47 |
buoyancy in a heated pipe | Volker Pawlik | FLUENT | 5 | February 8, 2001 10:53 |