-
Notifications
You must be signed in to change notification settings - Fork 41
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
k-epsilon model #57
Comments
What's y'all intended purpose here for the k-e model? (Natural circulation or convective flows?) |
Thanks for your interest @lcarasik ! Ideally, this tool should as generally as possible support turbulent or near-turbulent flow regimes in liquid fueled MSRs. This k-epsilon feature is not intended for a specific problem, but rather to support realistic simulation of the broad class of liquid fueled MSRs. Natural convection is a significant component of the flow in MSRs. But for ordinary operation, the flow is driven primarily by the pumps. To us, k-epsilon seemed the simplest appropriate model for this scenario, but any advice is welcome! |
I like the approach but I would warn you against trying to use it in transitional regimes. It might be worth looking into transitional RANS modeling if you get close. Is this 2-D? Or 3-D? Also, for finer tooth simulations where buoyancy is a big contributor to turbulence, y'all will want to look into how you handle the contributions of it in the k and e equations. |
Hey Lane,
I don’t think there’s a particular application of the TH in mind at the moment. Most of the prior work has been trying to reproduce MSRE results, which had a tiny Re=900 or so.
Going up to full-power liquid-fueled MSRs, the Reynolds number seems to be very turbulent, so transitional stuff probably isn’t a huge worry at the moment.
What’s the difference in 2D and 3D?
… On Dec 13, 2017, at 12:56 AM, Lane Carasik ***@***.***> wrote:
I like the approach but I would warn you against trying to use it in transitional regimes. It might be worth looking into transitional RANS modeling if you get close.
Is this 2-D? Or 3-D?
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub, or mute the thread.
|
Hey Gavin, I assume y'all are just using N-S (no RANS) for the MSRE work? It is still worth pointing out that if you get into transitional regimes (anything below 3000 for internal flows), you're going to see issues with standard RANS approaches. Turbulence and transitional behavior is inherently 3-D phenomena. You won't capture a lot in 2-D that you will in 3-D. Though it may not be important for y'all's work. |
Just for the record, this is false:
:) |
Which in particular then? Sorry, not that I was aware of, haha.
Gavin Ridley
… On Dec 13, 2017, at 9:52 AM, Katy Huff ***@***.***> wrote:
Just for the record, this is false:
I don’t think there’s a particular application of the TH in mind at the moment.
:)
—
You are receiving this because you commented.
Reply to this email directly, view it on GitHub, or mute the thread.
|
Depends on which of my proposals get funded. :) In seriousness though, a number of designs in the open moltres issues (some of which my current graduate students are pursuing) as well as MSRE itself (as noted in the review of our paper) could be expected to reach near turbulent regimes during design basis accident transients. |
It really sounds like you'd want someone to take a look at transition RANS after implementation of k-e! |
Oh haha, nice! Out of curiosity, if it's not classified info, how could the MSRE go turbulent? The only thing I can imagine would be somehow blocking half of the channels. Aside from that, the reviewer's claim of near turbulence didn't agree with @lindsayad's calculation at all (which I trust more). The pump ran at 80% power since the heat transfer calcs were slightly off, which led them to not run at the full 10 MW(th). So Re could max out around 1200, neglecting half of the channels being blocked. |
Yes, @lcarasik , particularly someone with your expertise! @gridley Not classified. I said 'near turbulent regimes during design basis accident transients' not 'turbulent in a steady state scenario.' In any case, there are plenty of examples of potentially laminar-to-turbulent even in not-far-off-normal MSRE experience. To quote the following about the MSRE experience.... "The regular pattern of the graphite support grid is discontinued near the centerline of the core where the control rods and the surveillance specimin holder are located. This allows greater salt velocities past these components for cooling.... The reynolds number of these channels is in the order of 3000 and is therefore in the transitional region." (http://moltensalt.org/references/static/downloads/pdf/ORNL-TM-3229.pdf) I'm not sure our published model took into account these larger central channels, as many don't. Additionally, not to let any cats out of any bags, but some atypical dynamic operation of some of these reactors (intentional or otherwise) might be constrained by the impact of rapid power oscillations on the stability of the velocity field. All of that is to say I don't think we can ignore near-turbulent behavior for interesting transients. Anyway, the geometry varies all over the place among the proposed designs. Some have pipes, some grids, some plates, some pools. The tool shouldn't be constrained to laminar regimes. |
I agree! Laminar is boring (comparatively). This should be exciting.
… On Dec 13, 2017, at 12:52 PM, Katy Huff ***@***.***> wrote:
Yes, @lcarasik , particularly someone with your expertise!
I said 'near turbulent regimes during design basis accident transients' not 'turbulent in a steady state scenario.' In any case, there are plenty of examples of potentially laminar-to-turbulent even in not-far-off-normal MSRE experience. To quote the following about the MSRE experience.... "The regular pattern of the graphite support grid is discontinued near the centerline of the core where the control rods and the surveillance specimin holder are located. This allows greater salt velocities past these components for cooling.... The reynolds number of these channels is in the order of 3000 and is therefore in the transitional region." I'm not sure our published model took into account these larger central channels, as many don't.
Additionally, not to let any cats out of any bags, but some atypical dynamic operation of some of these reactors (intentional or otherwise) might be constrained by the impact of rapid power oscillations on the stability of the velocity field. All of that is to say I don't think we can ignore near-turbulent behavior for interesting transients.
Anyway, the geometry varies all over the place among the proposed designs. Some have pipes, some grids, some plates, some pools. The tool shouldn't be constrained to laminar regimes.
—
You are receiving this because you commented.
Reply to this email directly, view it on GitHub, or mute the thread.
|
I'm working towards implementing and demonstrating this model for MSFRs modeling. I noticed that @lindsayad had already coded much of the work in the |
Just a quick drop in comment, please make sure the flow is fully turbulent if you apply it! If it isn't fully turbulent, you'll applying the turbulence model outside of it's applicability. |
@smpark7 let's chat in person about your question:
|
was accidentally previously filed in the wrong repo (arfc/publications#44).
From @lindsayad :
The k-epsilon turbulence model was never finished since it really required SUPG stabilization to be relevant. Now, SUPG and PSPG are implemented in the navier-stokes module. So the k-epsilon model could be finished if someone wanted to work on it. It's very close to completion; the major piece remaining are the wall functions.
While Dr. Lindsay is not able to work on this directly, consult with @lindsayad to get appropriate initial direction.
The text was updated successfully, but these errors were encountered: