CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > SU2

Gas flow through a channel with sudden radius expansion

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 16, 2014, 00:21
Default Gas flow through a channel with sudden radius expansion
  #1
New Member
 
Brazil
Join Date: May 2014
Posts: 8
Rep Power: 11
Blooper is on a distinguished road
(a brief introduction before the actual problem)

Hello fellow members!

I'm currently studying mechanical engineering, and one of my greatest areas of interest lies on fluid dynamics (and I'm still a complete rookie at it, mind you). My teacher asked me to try out/learn more about SU2 on his behalf, so here I am.

(the actual problem)

I went throught some of the tutorials and found them to be very good (only problem I found was that the configuration text files which describes the problems themselves weren't on a proper format on my notepad, I had to rearrange all the text myself, does anyone know how I could avoid such work?).

After this, I started trying to make a configuration file for my own problem, using one of the tutorials' configurations as a base file, but I'm finding some issues. The problem itself is one developed by Patankar, and goes as the title says: it's a gas flow inside a channel with radius r, and it suddenly expands to a radius R > r. It's also important to state that the tubes on first part of the flow (with radius r) are isolated, whereas the tubes on the second part, with radius R, are maintained on a Temperature T.

Before trying the mesh specific to this problem, I'm trying to make a configuration file to the tutorial on the channel with a bump and applying some changes (mostly, I'm trying to add the Temperature to the lower wall)
Link for the channel bump tutorial:
(http://adl-public.stanford.edu/docs/...p+in+a+Channel)

I tried applying the Navier Stokes equations to the physical problem (the tutorial uses Euler) and everything went OK. But when I try to change the wall conditions from MARKER_EULER to MARKER_ISOTHERMAL (in the proper format, with entity name and temperature), the solver crashes.

Does anyone have any idea what could be the problem?

I'll include the little changes I've made along with the original commands:

Original:

PHYSICAL_PROBLEM= EULER

FREESTREAM_PRESSURE= 101300.0

FREESTREAM_TEMPERATURE= 288.0

MARKER_EULER= ( upper_wall, lower_wall )

New:

PHYSICAL_PROBLEM= NAVIER_STOKES

KIND_TURB_MODEL= SA

FREESTREAM_TEMPERATURE= 300.0

REYNOLDS_NUMBER= 5000000.0

MARKER_ISOTHERMAL= ( upper_wall, 500.0 )

MARKER_ISOTHERMAL= ( lower_wall, 500.0 )

I'm inclined to believe that i didn't describe the boundary counditions properly, so if anyone could shed any light on the matter I'd be really thankful!
Blooper is offline   Reply With Quote

Old   May 18, 2014, 22:43
Default
  #2
Super Moderator
 
Francisco Palacios
Join Date: Jan 2013
Location: Long Beach, CA
Posts: 404
Rep Power: 15
fpalacios is on a distinguished road
Quote:
Originally Posted by Blooper View Post
(a brief introduction before the actual problem)

Hello fellow members!

I'm currently studying mechanical engineering, and one of my greatest areas of interest lies on fluid dynamics (and I'm still a complete rookie at it, mind you). My teacher asked me to try out/learn more about SU2 on his behalf, so here I am.

(the actual problem)

I went throught some of the tutorials and found them to be very good (only problem I found was that the configuration text files which describes the problems themselves weren't on a proper format on my notepad, I had to rearrange all the text myself, does anyone know how I could avoid such work?).

After this, I started trying to make a configuration file for my own problem, using one of the tutorials' configurations as a base file, but I'm finding some issues. The problem itself is one developed by Patankar, and goes as the title says: it's a gas flow inside a channel with radius r, and it suddenly expands to a radius R > r. It's also important to state that the tubes on first part of the flow (with radius r) are isolated, whereas the tubes on the second part, with radius R, are maintained on a Temperature T.

Before trying the mesh specific to this problem, I'm trying to make a configuration file to the tutorial on the channel with a bump and applying some changes (mostly, I'm trying to add the Temperature to the lower wall)
Link for the channel bump tutorial:
(http://adl-public.stanford.edu/docs/...p+in+a+Channel)

I tried applying the Navier Stokes equations to the physical problem (the tutorial uses Euler) and everything went OK. But when I try to change the wall conditions from MARKER_EULER to MARKER_ISOTHERMAL (in the proper format, with entity name and temperature), the solver crashes.

Does anyone have any idea what could be the problem?

I'll include the little changes I've made along with the original commands:

Original:

PHYSICAL_PROBLEM= EULER

FREESTREAM_PRESSURE= 101300.0

FREESTREAM_TEMPERATURE= 288.0

MARKER_EULER= ( upper_wall, lower_wall )

New:

PHYSICAL_PROBLEM= NAVIER_STOKES

KIND_TURB_MODEL= SA

FREESTREAM_TEMPERATURE= 300.0

REYNOLDS_NUMBER= 5000000.0

MARKER_ISOTHERMAL= ( upper_wall, 500.0 )

MARKER_ISOTHERMAL= ( lower_wall, 500.0 )

I'm inclined to believe that i didn't describe the boundary counditions properly, so if anyone could shed any light on the matter I'd be really thankful!

Thanks a lot for your interest in SU2,

Instead of

MARKER_ISOTHERMAL= ( upper_wall, 500.0 )

MARKER_ISOTHERMAL= ( lower_wall, 500.0 )

could you please try with

MARKER_ISOTHERMAL= ( upper_wall, 500.0, lower_wall, 500.0 )

Best Regards,

Francisco

PS.- If you are a rookie in CFD, my recommendation is to look for a similar problem in the TestCases folder and play with it before proposing something new.
fpalacios is offline   Reply With Quote

Old   May 19, 2014, 18:42
Default
  #3
New Member
 
Brazil
Join Date: May 2014
Posts: 8
Rep Power: 11
Blooper is on a distinguished road
Hey Francisco, thanks for the reply!

Indeed, when correcting the command format as you recommended, the simulation starts, unlike before. The problem now is that the script crashes after a few iterations (very few, 48 iterations). The last message is the following:

CSysSolve:modGramSchmidt: w[i+1] = NaN

I'm currently trying to work around this, but if you know what it's about then it would be of great help.

Thanks again!

Last edited by Blooper; May 19, 2014 at 18:47. Reason: typo
Blooper is offline   Reply With Quote

Old   June 5, 2014, 03:31
Default
  #4
hlk
Senior Member
 
Heather Kline
Join Date: Jun 2013
Posts: 309
Rep Power: 13
hlk is on a distinguished road
What you are describing is a symptom of the solution diverging.

Since you are using a config file which was originally for Euler it is likely that the CFL number and multigrid settings likely need to be changed. I recommend lowering the CFL number (try something a little less than one ramped up to 1 to start, if that works you can try increasing the CFL number to reduce the number of iterations) and turning of multigrid (MGLEVEL = 0).
hlk is offline   Reply With Quote

Old   June 5, 2014, 16:40
Default
  #5
New Member
 
Brazil
Join Date: May 2014
Posts: 8
Rep Power: 11
Blooper is on a distinguished road
Indeed hlk, I had already tried lowering the CFL number and it is now working again (I didn't even turn off multigrid - what would that change?). Anyway, all is running smoothly now, thanks a bunch for the reply.
Blooper is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
rhoCentralFoam for channel flow fportela OpenFOAM Running, Solving & CFD 22 June 10, 2014 21:14
Gravitational water flow in closed channel. Szymon85 CFX 7 September 3, 2013 17:28
Inlet boundary condition for gas dispersion in low Re channel flow vainilreb OpenFOAM Pre-Processing 0 February 22, 2013 04:03
Pressure loss due to sudden expansion in pipe flow Ahmed FLUENT 0 January 2, 2006 11:01
Transient natural gas flow description Leila FLUENT 0 November 29, 2003 17:06


All times are GMT -4. The time now is 13:22.