CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

shocks in a 1-d spherically symmetric model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 19, 2008, 22:31
Default shocks in a 1-d spherically symmetric model
  #1
Alex
Guest
 
Posts: n/a
Hi all,

I am a complete novice in the field of CFD and I am having trouble capturing shocks in the following situation. I would greatly appreciate any help or advise you can give me.

I am modeling a spherically symmetric collapse of a gas bubble in liquid. The bubble wall's motion is a prescribed function of time. The equations inside are compressible Navier-Stokes + thermal conduction in Lagrange coordinate formulation (lagrange coordinates are very convenient here since I have a boundary that moves with the gas). I am using 2nd order accurate explicit finite difference numerical scheme and it works just fine for relatively slow collapses that do not develop shocks. What is the best way to capture shocks in this situation? Also, the shock seems to be forming in a region within the bubble where density is low, and thus grid points are sparse (thanks to the Lagrange coordinates).

Thanks!
  Reply With Quote

Old   August 20, 2008, 03:27
Default Re: shocks in a 1-d spherically symmetric model
  #2
Ananda Himansu
Guest
 
Posts: n/a
A great deal is known about shock capturing and flux limiting in Eulerian coordinates. A lot of this should translate to Lagrangian coordinates, though I am unfamiliar with the literature for the latter. Check the wiki here at cfd-online (click on the link for the free CFD reference at top right of this webpage) or check Wikipedia.

Regarding the grid sparsity, should this not be addressable by suitably clustering gridpoints in your initial condition? You could examine your final solution or your solution before it goes haywire, and prescribe what would be a good gridpoint distribution for shock-capturing at that time. Then map this better grid back to the initial condition using the Lagrange coordinates. Now use this better initial grid to recompute the bubble collapse.
  Reply With Quote

Old   August 20, 2008, 06:53
Default Re: shocks in a 1-d spherically symmetric model
  #3
jinwon park
Guest
 
Posts: n/a
I wonder how to treat the interface between the bubble and liquid(?water). The nonphysical error can be made from incorrect interface treatment. If the flow is compressible, you could insert any flux limiter or slope limiter.

  Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
multiphaseInterFoam for RAS turbulence model chiven OpenFOAM Bugs 7 August 23, 2011 02:52
Problem importing watertight model into CFX djk1301 ANSYS 3 February 1, 2011 23:04
multi fluid mixture model issue rystokes CFX 3 August 9, 2009 19:13
Error in Axi symmetric model. binubtharayil FLUENT 2 August 5, 2009 16:01
LES and combustion model Margherita Cadorin CFX 0 October 29, 2008 06:24


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