Is PyFR suitable for supersonic cases?

[Ray]

Dear all,

I've tried a few of supersonic cases and failed each time. And there's no supersonic cases in published posts. I'm wondering if PyFR is suitable for supersonic cases

in normal working condition?

Best regards,

Ray

[Vincent, Peter E]

Hi Ray,

How did it fail? What Mach numbers were you trying?

We have run cases with relatively low Mach numbers (<2) with PyFR using artificial viscosity (see the User Guide):

• shock-capturing — shock capturing scheme:

  none | artificial-viscosity

  Cheers
  
  
  Peter
  
  
  Dr Peter Vincent MSci ARCS DIC PhD
  Reader in Aeronautics and EPSRC Fellow
  Department of Aeronautics
  Imperial College London
  South Kensington
  London
  SW7 2AZ
  UK
  
  web: www.imperial.ac.uk/aeronautics/research/vincentlab
  twitter: @Vincent_Lab
  
  
  
  
  
  
  

• --
  You received this message because you are subscribed to the Google Groups "PyFR Mailing List" group.
  To unsubscribe from this group and stop receiving emails from it, send an email to pyfrmailingli...@googlegroups.com.
  To post to this group, send email to pyfrmai...@googlegroups.com.
  Visit this group at https://groups.google.com/group/pyfrmailinglist.
  For more options, visit https://groups.google.com/d/optout.
  

  
  

[Ray]

Hi Peter,

I've tried this case, Mach 3 Wind Tunnel With a Step. It's in the book, AN ANALYSIS OF STABILITY OF THE FLUX RECONSTRUCTION FORMULATION WITH APPLICATIONS TO SHOCK CAPTURING. I got RuntimeError: NaNs detected at t = 0.05... Is there any problem in my configure file? 

Thanks!

Ray

[Park, Jin Seok]

Hi Ray,

I think Forward facing step problem with M=3 is one of the tough benchmark problems of shock-capturing methods.

It may yield negative value issues around the expansion corner and you may need some special cures (ex. refining mesh around the corner).

Furthermore, you didn't turn on the artificial viscosity option in solver setion as mentioned the previous Peter's mail.

The current shock capturing method (artificial viscoisty) of PyFR runs quite well for the flow with weak shock or oblique shock.

Please find the attached sample mesh and ini files. I think you can run similar supersonic cases.

Regards,

Jin Seok

Dr. Jin Seok Park PhD

---

From: pyfrmai...@googlegroups.com <pyfrmai...@googlegroups.com> on behalf of Ray <hera...@gmail.com>
Sent: 08 December 2016 10:31:12
To: PyFR Mailing List
Subject: [pyfrmailinglist] Re: Is PyFR suitable for supersonic cases?

 

--

[Ray]

Hi Jin Seok,

Thanks a lot! It works.

Regards,

Ray