Skip to content

Instantly share code, notes, and snippets.

@TristanWright
Last active April 15, 2016 18:41
Show Gist options
  • Save TristanWright/a695539b1a7d7d03e92aa7a5aaff556b to your computer and use it in GitHub Desktop.
Save TristanWright/a695539b1a7d7d03e92aa7a5aaff556b to your computer and use it in GitHub Desktop.

Simput - PyFr - euler-vortex-2d

Instructions to recreate euler-vortex-2d.ini in Simput. Unless mentioned below do not change default values.

Constants

  • delete cpTref
  • add three custom constants,
    • S = 13.5
    • M = 0.4
    • R = 1.5

Solver

Settings

  • order: 3
  • set Viscosity Correction to nothing
  • set Shock Capturing to nothing

Time Integrator

  • initial time = 0.0
  • final time = 50.0
  • time step = 0.005

Interfaces

  • delete LDG Beta
  • delete LDG Tau

Solver Interfaces

  • change Type to Linear

Solver Elements

  • Add one, change nothing

##Solution For each of these sections add an item and change the type to the section's name. Rename the property in the sidebar if desired.

Plugin NaN check

  • nsteps = 50

Plugin Writer

  • disk write time interval = 10.0
  • basedir = .
  • basename = euler_vortex_2d-{t:.1f}

ics

  • Initial Density = pow(1 - S*S*M*M*(gamma - 1)*exp(2*%(f)s)/(8*pi*pi), 1/(gamma - 1))
  • Initial X velocity = S*y*exp(%(f)s)/(2*pi*R)
  • Initial Y velocity = 1 - S*x*exp(%(f)s)/(2*pi*R)
  • Initial static pressure dist =
    1/(gamma*M*M)*pow(1 - S*S*M*M*(gamma - 1)*exp(2*%(f)s)/(8*pi*pi), gamma/(gamma - 1))
  • add a helper function: f = ((1 - x*x - y*y)/(2*R*R))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment