CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

forward step: treating singularity pt for WENO FV

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

LinkBack Thread Tools Display Modes
Old   January 16, 2007, 11:20
Default forward step: treating singularity pt for WENO FV
Posts: n/a

I have implemented the WENO FV scheme. When I try to simulate the Forward facing step (Mach 3), I found that the code blows after certain number of iterations.

I checked that the pressure goes to negative values at some time near the corner of the step.

The paper by Shu mentions about applying an extra boundary condition at the singularity point (corner of the step). But, I am not sure if that is the problem with my current simulation.

I apply the wall boundary condition by reflecting the values on all the ghost cells and also at only one ghost cell, while maintaining others at a near infinite value. Both gives similar solutions. But, which one is best?

If anyone had tried this problem, kindly share your experience.

Thank you very much for your help.
  Reply With Quote

Old   January 18, 2007, 03:15
Default Re: forward step: treating singularity pt for WENO
Posts: n/a
I suggest you to try to fix the singularity, then run again. If problem persists, there has to be a bug somewhere in the code. That is because, I implemented that case with a multiblock grid and it worked well. There shouldn't be a need for an order reduction for Mach 3.

Fixing the singularity which is at the corner of the step is another issue. Try to reset the values according to the the paper of Woodward and Colella. Like, resetting density variable in the first four cells in the first row beginning from the corner, and resetting magnitudes of velocities not their directions in the first two cells, which are located in the second row just above the step.
  Reply With Quote


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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
Transient simulation not converging skabilan OpenFOAM Running, Solving & CFD 12 September 17, 2007 17:48
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 18:35.