CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > Siemens > STAR-CCM+

Segregated flow not converging properly - DARPA Suboff submarine

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 18, 2017, 14:28
Default Segregated flow not converging properly - DARPA Suboff submarine
  #1
New Member
 
Join Date: Jan 2017
Posts: 13
Rep Power: 9
Jon Faried is on a distinguished road
Hi, everyone!

I'm currently trying to validate the benchmark DARPA Suboff (submarine) case using STAR-CCM+, and am at a loss as to why I haven't been able to get the residuals to converge properly and for the results to come out correctly. Any insight would be appreciated.

And now for the details...

The mesh is a simple axi-symmetric one (only the bare hull of the sub is currently being tested, no fins and rudders) with a y+ equaling roughly 1, with the axes removed and replaced with H-meshes. I'm afraid I don't have any photos currently handy, but the completed mesh looks like a dome (or half a sphere).

I've managed to successfully run the exact same mesh on ANSYS CFX and obtained excellent results, so I don't think the mesh is the problem, unless there's something about meshes that requires special treatment in STAR-CCM+, in which case I'd value your input.

The submarine has no angle of attack, and is exposed to a constant flow with a Reynold's number of 14,000,000. The mesh is split into four main boundaries: the hull (no-slip wall), the symmetry plane (only half of the axisymmetric mesh is necessary), an inflow section (half of the dome, where the fluid comes in with a constant velocity), and an pressure outlet section (backflow direction is extrapolated; pressure is static, though 'environmental' was tried as well).

I'm currently using the segregated flow model (RANS, K-Omega SST, Liquid). The turbulence residuals are fine (quickly reaching E-8) whereas the mass and momentum residuals drop quite slowly (down to E-3, couple thousand iterations) before eventually (approx. 70000 iters) jumping back up to E-1 and staying put. One way or another, for a simple case like this one I'd expect a fairly quick convergence, so something seems definitely off here.

The pressure and velocity results seem to be OK throughout the entire process, whereas the turbulent viscosity eventually goes haywire and produces wrong results. The coefficient of drag completely misses the mark as well.

I've tried running the laminar case in order to see if turbulent viscosity was the problem, but the case completely diverges. Trying the Coupled Flow solver doesn't work either - it too diverges (even for relatively low Courant numbers).

I get the feeling that I'm missing something trivial here... any help would be most appreciated, and I'd be happy to give more details if necessary. Thanks in advance!
Jon Faried is offline   Reply With Quote

Old   January 19, 2017, 16:25
Default
  #2
Senior Member
 
Join Date: Nov 2010
Location: USA
Posts: 1,232
Rep Power: 25
me3840 is on a distinguished road
You have a dome but you are using velocity and pressure conditions? That's probably the issue. You will have velocity traveling tangent to its face at some point connecting to a cell with zero gradient and a pressure specified. That's just a nightmare. I would switch your domain to be a cylinder instead.
me3840 is offline   Reply With Quote

Old   January 19, 2017, 17:07
Default
  #3
New Member
 
Join Date: Jan 2017
Posts: 13
Rep Power: 9
Jon Faried is on a distinguished road
Quote:
Originally Posted by me3840 View Post
You have a dome but you are using velocity and pressure conditions? That's probably the issue. You will have velocity traveling tangent to its face at some point connecting to a cell with zero gradient and a pressure specified. That's just a nightmare. I would switch your domain to be a cylinder instead.
Thanks for the reply!

While what you're saying certainly makes sense, I get the feeling that this isn't the issue at hand, since:

A. The same mesh (boundary conditions included) was tried using CFX and another solver, and the results in both cases were accurate. Is there a fundamental difference that you know about between the way CFX and STAR-CCM+ treat boundary conditions? (or in general, for that matter?)

B. What's more, the velocity and pressure solutions seemed to be quite accurate, even when the residuals jumped from E-3 to E-1. There certainly wasn't anything amiss at the interface between the two dome shaped boundaries (inlet-outlet), I checked that area.

Have you ever encountered a particular mesh which one solver tackled effectively whereas STAR-CCM+ had difficulty with it? If so, I'd love to hear about it, perhaps I can glean a bit of info that'll be practical in this case.

Again, thank you for providing your two cents, much appreciated.

Sent from my A0001 using CFD Online Forum mobile app
Jon Faried is offline   Reply With Quote

Old   January 19, 2017, 18:54
Default
  #4
Senior Member
 
Join Date: Nov 2010
Location: USA
Posts: 1,232
Rep Power: 25
me3840 is on a distinguished road
Quote:
Originally Posted by Jon Faried View Post
Thanks for the reply!

A. The same mesh (boundary conditions included) was tried using CFX and another solver, and the results in both cases were accurate. Is there a fundamental difference that you know about between the way CFX and STAR-CCM+ treat boundary conditions? (or in general, for that matter?)

B. What's more, the velocity and pressure solutions seemed to be quite accurate, even when the residuals jumped from E-3 to E-1. There certainly wasn't anything amiss at the interface between the two dome shaped boundaries (inlet-outlet), I checked that area.
A. Yes, I've seen this many times. There's lots of "stuff" that happens behind every commercial code which makes BCs and discretization schemes that look the same act differently. If something is theoretically a poor idea but "just works" in a commercial code, it's often because the code is designed to help people who don't follow best practices.

B. Define accurate here? Remember the residual itself is only a measure of relative error. Perhaps your initial solution is very good? Perhaps your mesh is coarse enough to diffuse the physics in a location to give a stable but poorly converged solution? Have you looked at the unnormalized residuals? They could be quite low while the relative residual does not converge well.

I'm surprised you don't see any issues in that area, perhaps that means I don't understand how the domain is set up after all, a picture may help. In STAR-CCM+ you can use the 'temporary storage' option on the solver to enable you to view the residuals as a function of space, I would recommend you do that here to see where the largest residuals are.
me3840 is offline   Reply With Quote

Reply

Tags
coupled solver, darpa suboff, residuals, segregated flow, star ccm+

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
Re- 3.4x10^7- flow wont run properly shafee Main CFD Forum 12 March 25, 2016 19:05
Coupled flow vs. segregated ? cfdblender Main CFD Forum 3 December 22, 2015 12:09
stability of converging flow Bo Busk Jensen Main CFD Forum 2 October 17, 2001 03:38
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 22:31
Flow visualization vs. Calculated flow patterns Francisco Saldarriaga Main CFD Forum 1 August 3, 1999 00:18


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