shocks in a 1-d spherically symmetric model

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

 August 19, 2008, 22:31 shocks in a 1-d spherically symmetric model #1 Alex Guest   Posts: n/a Sponsored Links 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!

 August 20, 2008, 03:27 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.

 August 20, 2008, 06:53 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post chiven OpenFOAM Bugs 7 August 23, 2011 02:52 djk1301 ANSYS 3 February 1, 2011 23:04 rystokes CFX 3 August 9, 2009 19:13 binubtharayil FLUENT 2 August 5, 2009 16:01 Margherita Cadorin CFX 0 October 29, 2008 06:24