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

WENO for traveling/steepening acoustic wave

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 19, 2009, 00:58
Default WENO for traveling/steepening acoustic wave
  #1
New Member
 
Nathan Grube
Join Date: Mar 2009
Posts: 13
Rep Power: 8
ngrube is on a distinguished road
Hello,

I am puzzled by the failure of a well debugged and validated WENO finite difference code in the case of a traveling acoustic wave propagating over a long distance. I expected the originally sinusoidal wave to steepen into a sawtooth or N-wave pattern. What I see instead is the development of shocks (where they should be) plus somewhat weaker but still significant shocks approximately halfway in between, so I get twice as many shocks as I should.

I might faintly recall seeing a discussion in a journal article about the inability of WENO to handle this case. Or maybe I just dreamed that. I have been searching with all the keywords I can think of and can find nothing. Any help would be appreciated.

Thanks,
Nathan
ngrube is offline   Reply With Quote

Old   March 19, 2009, 11:09
Default
  #2
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 188
Rep Power: 8
harishg is on a distinguished road
Did you check if the method worked with lower-order shock capturing schemes ? I recall a discussion on this topic in the book of Leveque for long distance propagation.
harishg is offline   Reply With Quote

Old   March 19, 2009, 11:16
Default Leveque
  #3
New Member
 
Nathan Grube
Join Date: Mar 2009
Posts: 13
Rep Power: 8
ngrube is on a distinguished road
Thanks for your reply.

No, I have not tried lower order shock capturing on it. I need high order since this code will be used for turbulence. This case is just a test case I wanted to use to see how a new turbulence model behaves in the presence of these shocks. (But this test case is 1-dimensional and hence laminar of course.)

I will look for something relevant from Leveque. I have requested his finite difference book from another library. Will post again if I learn anything.

Last edited by ngrube; March 19, 2009 at 15:27.
ngrube is offline   Reply With Quote

Old   March 19, 2009, 15:25
Default Time integration
  #4
New Member
 
Nathan Grube
Join Date: Mar 2009
Posts: 13
Rep Power: 8
ngrube is on a distinguished road
It just occurred to me that if I end up with twice as many shocks as I should have, then the total variation is not diminishing. I have been using a Williamson 3rd order low-storage Runge-Kutta scheme for time advancement, and now that I think about it, I am not certain that that is TVD. I know many people use RK with WENO methods, but maybe the low storage RK that I am using is not suitable?

Last edited by ngrube; March 19, 2009 at 15:32. Reason: typo
ngrube is offline   Reply With Quote

Old   March 19, 2009, 15:36
Default
  #5
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 188
Rep Power: 8
harishg is on a distinguished road
There is a version of the low-storage TVD R-K scheme which was implemented by Shu and Gottelieb. If you google for it, you might find the report. The book of leveque is " FVM for hyperbolic problems". He talks about the possible failure of linearized Riemann solvers such as Roe one in that. I am not exactly sure though if thats the exact thing that you might be looking for.
harishg is offline   Reply With Quote

Old   March 19, 2009, 17:00
Default TVD low-storage Runge-Kutta
  #6
New Member
 
Nathan Grube
Join Date: Mar 2009
Posts: 13
Rep Power: 8
ngrube is on a distinguished road
Got it. Thanks. For anyone else who is interested, the reference is

Total Variation Diminishing Runge-Kutta Schemes
by Sigal Gottlieb and Chi-Wang Shu
Mathematics of Computation vol 67, number 221, January 1998, pp 73-85

Also, the Leveque info is in section 15.3.7, "Failure of Linearized Solvers." Not sure yet whether that is pertinent. I plan to try a TVD RK scheme first. Will post with the results when done.
ngrube is offline   Reply With Quote

Old   March 20, 2009, 06:51
Default Better reference
  #7
jed
Member
 
Jed Brown
Join Date: Mar 2009
Posts: 56
Rep Power: 10
jed is on a distinguished road
This field (now known as "strong stability preserving" time discretization, WENO is not strictly TVD by the way) has evolved significantly in the last decade. A good current reference is

Code:
@article{gottlieb2009hos,
  title={{High order strong stability preserving time discretizations}},
  author={Gottlieb, S. and Ketcheson, D.I. and Shu, C.W.},
  journal={Journal of Scientific Computing},
  volume={38},
  number={3},
  pages={251--289},
  year={2009},
  publisher={Springer}
}
a preprint is at

http://www.amath.washington.edu/~ketch/David_Ketcheson/Publications_files/SSPreviewFINAL.PDF

if you don't have access to the journal.

A selection of optimal SSP Runge-Kutta and multistep (explicit and implicit) are tabulated at
http://www.amath.washington.edu/~ketch/SSP/
jed is offline   Reply With Quote

Old   March 20, 2009, 09:28
Default
  #8
New Member
 
Nathan Grube
Join Date: Mar 2009
Posts: 13
Rep Power: 8
ngrube is on a distinguished road
Thanks, Jed, for those very helpful references and for cluing me in to the term "SSP." As you say, it is true that WENO isn't TVD of course. I was just hoping that the time integration was responsible for the non-TVD behavior since that is much easier to fix than moving away from the WENO scheme, which I am pretty much obligated to use.
ngrube is offline   Reply With Quote

Old   March 20, 2009, 17:11
Default Problem seems to be with WENO, not RK
  #9
New Member
 
Nathan Grube
Join Date: Mar 2009
Posts: 13
Rep Power: 8
ngrube is on a distinguished road
Well, I've coded a TVD RK scheme and get almost identical results. I would say that my problem is with trying to use WENO for this.
ngrube is offline   Reply With Quote

Old   March 20, 2009, 20:31
Default
  #10
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 188
Rep Power: 8
harishg is on a distinguished road
Start with first order upwind, Godunov schemes and then try higher order scheme. I think its also possible to reduce order of ENO scheme by using a smaller stencil size.
harishg is offline   Reply With Quote

Old   March 9, 2015, 04:59
Unhappy toothy turbulent structures
  #11
Member
 
Hassan Nemati
Join Date: Apr 2012
Posts: 30
Rep Power: 5
hnemati is on a distinguished road
Hi all
I am simulating DNS of developing pipe flow with highly temperature dependent property. I saw-tooth this structures in my enthalpy and velocity fluctuations. My simulations is highly resolved and I know that they are physical. I could not find litretures to get help from it in this issue. Can somebody help me in this case? or introduce some refs.?
Good day guys
Attached Images
File Type: jpg Screenshot from 2015-03-09 09:51:58.jpg (66.9 KB, 7 views)
hnemati is offline   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
ENO & WENO jha Main CFD Forum 5 March 4, 2008 08:04
WENO-1D Hyperbolic Bob Main CFD Forum 0 January 16, 2007 16:21
acoustic wave in water Viktor Bessonov Main CFD Forum 0 February 3, 2005 05:56
Acoustic wave Viktor Bessonov Main CFD Forum 2 January 28, 2005 05:16
Tue, Apr 2000 ?, traveling backward in time? John C. Chien Main CFD Forum 1 May 26, 2000 14:07


All times are GMT -4. The time now is 06:38.