# 153 -1 #print version 1 WELCOME TO WIND TUNNEL CALCULATIONS FronTier Gas version 20020808 DATE OF RUN Wed Apr 30 21:00:00 2003 # 39921 0 ##MACHINE PARAMETERS # Hostname = a15333 # Operating System = Linux # OS Release = 2.4.18 # OS Version = #1 SMP Mon Sep 9 15:17:50 EDT 2002 # CPU Type = i686 # Byte Ordering = Little Endian # Floating Point Word Size = 8 Type 'p' to obtain prompting on input: Type debug to Turn On Debugging: debug Specify Debugging Destination, screen or file: FILE List Functions to Debug - or all - end ends list Names following all will not be debugged - end ends :fill_hdf_values2d :N_half_step :T_half_step :fatalbad : bad_state : end Enter the spatial dimension of the computation: 2 Enter the remapping (Jacobian) converting the physical coordinate system to the cartesian computational coordinate system. The choices are Identity Remap (default) (i) Cylindrical Geometry (c) Type Choice Here: c Enter the computational limits in the radial direction, RL, RU: 0 0.15 Enter the computational limits in the vertical direction, ZL, ZU: 0 0.15 Enter the numbers of grid intervals of the computational grid in the r & z directions: 50 50 The topological grid is a grid used for the construction of the tracked front topology. It is constrained to be a square grid. You specify the grid in one of two ways. If you enter a single number, it will be used as a coarseness factor for the topological grid relative to the computational grid entered above. In this case the length of a topological grid block cell side is the nearest allowable multiple of the shortest side of the computational grid by the coarseness factor. Otherwise the code will read the two integers input for the number of grid cells in each coordinate direction of the topological grid. If your input values do not yield a square grid they will be corrected to produce a square grid. This correction will attempt to produce values close to those input, but if the input values are highly rectangular, the resulting values may differ considerably from those entered. The default for this input option is the nearest square grid that matches the computational grid. Generally the topological grid is coarser than the computational grid. Larger coarseness factors yield coarser grids, a value one gives the nearest square grid to the computational grid. Enter your choice (cor_fac, two integers, or return) (defaults are 200 200): 3 The topological mesh used is 67 67 General run termination/pause conditions Enter limits on real time (max time), mesh time (max timesteps), an optional initial time, and an optional stop time mode (exact or constant), (dflt = inf 2147483647 0 constant): 0.15 2 Specify the pause time mode [exact, constant(dflt), mesh]: Enter the first Pause Time (dflt = inf): Enter maximum number of time step modifications allowed during a propagation step (default = 50): 25 Select triangulation option, exact [e], fast [f], or pcs [p, dflt]: p Printing Control Request main output format(s). Options are front_plots only (F) front_plots plus interior_states (Restart format) (R) front_plots plus tri_plots (T) front_plots, interior_states, and tri_plots (A) HDF raster plots (H) SDS files (S) PROSTAR plots (P) suppress output (dflt) Enter the choices as a space separated list: R H Prompt for front_plots printing control. Specify the interval type for printing [exact, constant, mesh (default)]: mesh Enter the step interval and first step for printing (default = 1 0): 100 0 Request binary/non-binary output [b(dflt),n]: Prompt for interior_states printing control. Specify the interval type for printing [exact, constant, mesh (default)]: mesh Enter the step interval and first step for printing (default = 2 0): 100 0 Request binary/non-binary output [b(dflt),n]: Prompt for HDF_plots printing control. Specify the interval type for printing [exact, constant, mesh (default)]: const Enter the time interval and first time for printing (default = 1 0): 0.01 Request binary/non-binary output [b(dflt),n]: n The user can request that restart dumps be printed at a specified wall time interval. These dumps will be named out/Implosion/impf16RM.lastdump0 and out/Implosion/impf16RM.lastdump1 and will be alternately overwritten as the run proceeds The wall time dump frequency can be given in units of seconds, minutes (default), or hours. Indicate the units in the obvious way such as 30 minutes, 2 hours, etc. To request this option enter the wall time print frequency: 20 hours Print wall time dumps in binary (default = no): Request composition type of materials. Available types are PURE_NON_REACTIVE (PNR, default) MULTI_COMP_NON_REACTIVE (MCNR) Enter choice here: HDF plotting initialization Specify variables to be plotted (HDF). The choices are -- r-velocity (VR), z-velocity (VZ), flow speed (U), kinetic energy density (K), momentum (M), density (D), energy density (E), vorticity (VORTICITY), divergence of velocity (W), pressure (P), specific internal energy (I), temperature (T), specific enthalpy (H), entropy (S), radial component of velocity (RV), tangential component of xy velocity (TV), self similar Mach number (L), or Mach number (N). Enter choices as a space separated list, using multiple lines if needed. Terminate all lines EXCEPT THE LAST with a backslash '\' Enter choices: D Enter the coordinates of lower corner of the initial view box (dflt = 0 0): Enter the coordinates of upper corner of the initial view box (dflt = 0.25 0.25): Enter the velocity of the view box (dflt = 0 0): Specify the number of pixels in the x and y directions (dflt 600 600): Enter an optional plotting filter for DENSITY, choices are Identity filter (none, default) Log filter (log) Log1p filter (log1p) Exp filter (exp) Expm1 filter (expm1) Atan filter (atan) Tan filter (tan) Enter choice: log Enter optional scaling factors for LOG_DENSITY: Color palelette data for LOG_DENSITY can be entered as either an HDF color palette file, or as a raw color palette consisting of a binary list of unsigned chars in the form of 256 red values, followed by 256 green values, followed by 256 blue values. A raw format file is indicated by appending a blank and the indicator raw after the file name. Otherwise any file entered will be assumes to be in HDF palette format. Enter an optional color palette file: You will now be prompted for a base file name and optional directory for the HDF output. Output for each variable is to a separate file whose name contains the base name, and the prompt string for that variable. Enter a file name for the output for HDF data (default = out/Implosion/hdf/impf16RM/impf16RM): Enter the compression type, choices are None (N) Run length encoding, no data loss (R, default) JPEG, some data loss (J) Enter choice: N End HDF plotting initialization Specify additional GAS printing variables. The choices are, VELOCITY (V), RADIAL_COMPONENT_OF_VELOCITY (RV), PRESSURE (P), SPECIFIC_ENTROPY (SPECIFIC_ENTROPY), SOUND_SPEED (SOUND_SPEED), TEMPERATURE (TEMPERATURE), TANGENTIAL_COMPONENT_OF_XY_VELOCITY (TV). Enter choices as a space separated list, using multiple lines if needed. Terminate all lines EXCEPT THE LAST with a backslash '\' Enter choices: Specify interior hyperbolic difference method. Select the difference method for solving the hyperbolic system in the interior regions away from the fronts. Choices are Split Lax-Wendroff (LWS) Split Lax-Friedrichs (LFS) Split first order Godunov (G) Vectorized split Lax-Wendroff (VLS) Five point Vectorized split MUSCL (VM) Colella Piecewise Linear Method (PLM) Vectorized pseudo unsplit Lax-Wendroff (PUSLW) Vectorized pseudo unsplit MUSCL (PUSM) Colella pseudo unsplit Piecewise Linear Method (PUSPLM) Five point Vectorized 2D Unsplit MUSCL (UNVM) Enter choice here: UNVM Current values for MUSCL parameters State reconstruction = Reconstruct density, energy, velocity Riemann flux solver = Exact Riemann solver No method of characteristic solver used Irregular hyp stencil method = MUSCL Tangential sweep method = MUSCL Don't test for negative density and energies at half step Enforce monotone reconstructions at cell edges = no Link reconstructions (zero slope in one field implies zero slope in all fields) = no End List of current values for MUSCL parameters Use all defaults for MUSCL code (dflt = y): n Choose the desired type of linear reconstruction, Choices are Reconstruct density, energy, velocity (d, default) Reconstruct eigen coordinates (e) Bell-Colella-Trangenstein reconstruction (b) First order Godunov reconstruction (zero slopes) (f) Enter choice: d Test for negative density and energies at half step (dflt=no): no Choose the desired Riemann solver, Choices are Exact Riemann solver (e, default) Linear approximate Riemann solver (l) Colella-Glaz's approximate Riemann solver (c) Linear US/UP fit (Dukowicz) (d) Gamma Law fit (g) Enter choice: C Current values of the Colella-Glaz Riemann solver parameters Maximum number of iterations for Riemann solver 4 Minimum pressure jump below which mass flux 1e-06 is replaced by acoustic impedance Velocity convergence tolerance 1e-06 The strong wave tolerance above which approximate 100 Riemann solver is replaced by an exact solver Minimum allowed mass flux 1e-12 End List of Colella-Glaz Riemann solver parameters Use all defaults for Colella-Glaz Riemann solver (dflt = y): n Enter the number of iterations for the Riemann Solver (dflt = 4): Enter the minimum pressure jump below which the mass flux is replaced by the acoustic impedance (dflt = 1e-06): Enter the velocity convergence factor (dflt = 1e-06): Enter the strong wave tolerance to turn on exact solver (dflt = 100): Enter the minimum allowed mass flux (dflt = 1e-12): Choose the irregular stencil method, Choices are Lax_Wendroff (l) MUSCL (m, default) Enter choice: M Choose the tangential sweep method, Choices are Lax_Wendroff (l) Lax-Friedrichs (lf) First order godunov (g) MUSCL (m, default) Enter choice: G Current defaults for artificial viscosity parameters--- Artificial Viscosities and Heat Conductions Use nonlinear artificial viscosity = no Coefficient of nonlinear artificial viscosity = 0 Use linear artificial viscosity = no Coefficient of linear artificial viscosity = 0 Use upwind artificial viscosity = no Coefficient of upwind artificial viscosity = 0 Use MUSCL slope flattening = yes Muscl slope flattening parameter eta inverse = 2 Muscl slope flattening minimum shock strength = 0.25 Muscl slope flattening minimum shock specific volume jump = 1e-06 Coefficient of artificial heat conduction = 0 Muscl slope flattening charateristic speed cutoff = 0 Coefficient of dynamic surface tension = 0 Artificial viscosity stability coefficient = 1 Use current defaults for artificial viscosity parameters (dflt = y, type p to print defaults): n Do you wish to use slope flattening at strong waves (dflt yes, d = global defaults)): To have a global default coefficient of the wave speed weight eta, (0(max limiting) < eta < 1) enter the coefficient (default = 0.5, type d to use 0.5): To have a global default coefficient of the minimum shock jump tolerance enter the coefficient (default = 0.25, type d to use 0.25): To have a global default coefficient of the minimum specific volume jump tolerance enter the coefficient (default = 1e-06, type d to use 1e-06): To have a global default coefficient of the scaled characteristic speed cutoff enter the coefficient (default = 0, type d to use 0): To have a global default coefficient of Lapidus nonlinear artificial viscosity enter the coefficient (default = 0, type d to use 0.5): 0.5 To have a global default coefficient of linear artificial viscosity enter the coefficient (default = 0, type d to use 0.05): 0.05 To have a global default coefficient of upwind artificial viscosity enter the coefficient (default = 0, type d to use 0.1): To have a global default coefficient of artificial heat conduction enter the coefficient (default = 0, type d to use 0.05): Dynamic surface tension is used to stabilize contacts against shear instabilities. The value for this coefficient should be of the same magnitude as the wavelength (in zones) of the instabilities you wish to suppress, i.e. 2-4 zones. To have a global default coefficient of dynamic surface tension enter the coefficient (default = 0, type d to use 0): The following choices are available for a gravitational acceleration No gravity (N or default) Constant gravity (C or Y) Time dependent gravity (T) Astrophysical (central force) gravity (A) Radial gravity with constant magnitude (R) Enter choice: Dynamic tracking decision variables Tracking decisions on dynamically produced waves are based on a floating point cutoff on the wave strength. Strengths are normalized to zero for weak waves, so a tolerance of 0.0 will always signal tracking, and a very large tolerance will always signal not to track. For each interaction type, you will be asked to enter the cutoff tolerance and a tolerance type for determining whether scattered waves of the indicated type should be tracked when produced by a specific bifurcation type. The currently supported tolerance types are Never track (Never) Always track (Always) Pressure ratio minus one across the wave (Pressure) Absolute value of the Atwood number across the wave (Atwood) Mach number minus one for the state ahead of the wave (Mach) Wave strength tolerance for tracking reflected shocks at regular reflections = always track Wave strength tolerance for tracking reflected shocks at attached boundary reflection nodes = always track Wave strength tolerance for tracking the slip line at Mach reflections = always track Wave strength tolerance for tracking reflected shocks at Mach reflections = always track Wave strength tolerance for tracking the Mach stem at Mach reflections = always track Wave strength tolerance for tracking slip lines produced by shock crossings = always track Wave strength tolerance for tracking reflected shocks at shock crossings = always track Wave strength tolerance for tracking slip lines at shock overtakes = always track Wave strength tolerance for tracking reflected shocks at shock overtakes = always track Wave strength tolerance for tracking transmitted shocks at shock overtakes = always track Wave strength tolerance for tracking reflected rarefaction leading edges at shock overtakes = always track Wave strength tolerance for tracking reflected rarefaction trailing edges at shock overtakes = always track Wave strength tolerance for tracking material interfaces at shock-contact diffractions = always track Wave strength tolerance for tracking reflected shocks at shock-contact diffractions = always track Wave strength tolerance for tracking transmitted shocks at shock-contact diffractions = always track Wave strength tolerance for tracking reflected rarefaction leading edges at shock-contact diffractions = always track Wave strength tolerance for tracking reflected rarefaction trailing edges at shock-contact diffractions = always track Wave strength tolerance for tracking material interfaces at shock-contact transmission nodes = always track Wave strength tolerance for tracking transmitted shocks at shock-contact transmission nodes = always track Don't Turn off tracking of regular reflection node if node propagation fails Don't Turn off tracking of Mach node if node propagation fails Don't Turn off tracking of overtake node if node propagation fails Don't Turn off tracking of precursor rr diffraction (cluster) if node propagation fails End Dynamic tracking decision variables Use default settings for dynamic tracking (default = y): y Type yes to request automatic wave capture: Point propagation parameters have the current default values Operator split normal/tangential update Current values for options for g_npt_w_speed A wave is defined to be strong if |1 - (1/(rho*c)*|dp/du|| > Mach_tol or |rhol - rhor|/(rhol+rhor) > A_tol Mach_tol = 0.25 A_tol = 0.25 Neumann boundary states are computed by an average of a reflection symmetry contact propagation and a method of characterics calculation. The weight of the symmetry contact result is proportional to the flow gradient. The wall limiter value gives this proportionality constant. Wall_limiter = 1 vector_moc = MOC_PLUS_RH scalar_moc = RIEMANN without filtering of outgoing waves vector_ahead_state_moc = 0x81a32a0 shock_ahead_state_riem_inv_moc neumann_moc = 0x81a2bf4 neumann_riem_inv_moc geom_source_method = BACKWARD_EULER End current values for options for g_npt_w_speed Use defaults for point propagation operators (dflt = y): y The current defaults for the linear interpolation options are Linear interpolation based on conserved variables Volume interpolation coefficients Use current defaults for linear interpolation options (default = y): y Enter an upper bound for the number of components (default = 100): Specify initial interface of tracked curves Choices are Input interface by hand (type `screen') Input interface from a file (restart option - enter filename) Request default option(s) (hit `return') Enter choice: 2D front redistribution control Enter tracking algorithm, choices are: Grid free tracking(F), Grid based tracking (G), Enter choice: G front spacing control Enter the spacing for general curves in dimensionless length/mesh units (dflt = 0.75): Enter the spacing for vector type curves in dimensionless length/mesh units (dflt = 0.75): Small loop control Reflect small loop shocks (dflt = no): time step size control The current defaults for the front time step control are Time step factor = 0.75 Apply CFL at nodes = yes Maximum node separation at untangle = 2 CFL increase factor = 1.1 Minimum time step modification factor = 0.75 Maximum time step modification factor = 1.25 Use defaults for front time step control (y dflt): n Enter the time step factor (fraction of CFL condition - default 0.75): 0.45 Use node velocity to restrict CFL condition (default YES): Enter the maximum node separation at tangles (default 1.5): Enter the CFL increase factor (default 1.1): Enter the minimum time step modification factor (default 0.75): Enter the maximum time step modification factor (default 1.25): flow specified state enforcement at fronts Enforce flow specified states at fronts (dflt=yes): Type yes to propagate front in moving frame (dflt = no): Select tangent computation algorithm, choices are Linear centered SECANT vectors (SECANT, default) Fourth order LANGRANGIAN interpolation (LANGRANGIAN) Cubic SPLINE fit (SPLINE) Enter choice: To use curvature dependent limiting at scalar fronts enter the curvature factor (dflt = 0): Request problem type. Current choices are an ambient state test (AM), a trigrid test (TRI), a plane front (P), a bowshock (BO), a Meshkov instability problem (M), a shock diffraction (D), a shock transmission (T), a ramp reflection problem (RR), a contact-contact interaction (CC), a Richtmyer linear theory (RL), an astrophysical jet (AJ), an injection inlet jet (IJ), a gas injection jet (FJ), a neutrino booster colapse (NB), a Supernova simulation (SN), an imploding elliptical shock (IMP), a shock running over an expanding ramp (X), a random surface instability problem (RS), a shocked thermal layer (STL), a Richtmyer-Meshkov instability problem (RM), a Rayleigh-Taylor instability problem (RT), a bubbles and drops problem (BD), an expanding shells (ES), shock jet interaction (SJ), a Radial Rayleigh-Taylor instability problem (Radial Rayleigh Taylor), or a Kelvin-Helmholtz instability problem (KH). Enter choice here: IMP Request Neumann Wall to cut off zero radius [y,n(dflt)]: no Enter the number (>= 1) of elliptical waves (dflt = 1): 2 The inner ellipse can be chosen to as closed (full) (C, F, default), half (H), quarter (Q), or angle sector specified (A). Enter choice: Q This ellipse has been oriented in the counter clockwise direction so that the forward moving normal directions point towards the outside of the ellipse. Enter the 2 coordinates of the center of the inner ellipse: 0 0 Enter the common radius or the 2 mean radii of the inner ellipse: 0.1 0.1 Enter a scaling factor for the radius (dflt = 1): Type 'y' to rotate the inner ellipse: Enter the rotation angle (in degrees) about the axis of rotation: Type 'y' if the inner ellipse is to be fourier polynomial perturbed: y Should the perturbation be (R)andom (default), (I)nput, or (M)ixed input and random?: R Enter the number of modes or the minimum and maximum mode numbers for this random ellipse (dflt = 1): 3 6 Enter the amplitude standard deviation: 0.001 Use sine weighting on amplitudes (dflt n): n Enter the average phase: 90 Enter the bubble phase standard deviation: 0 Enter an optional three short integers for the amplitude random number generator seed: 22 31 16 Enter an optional three short integers for the phase random number generator seed: Type 'y' if the inner ellipse is to be Legendre polynomial perturbed: n Possible choices for wave types are Forward shock wave (f), Forward sound wave leading edge (fl), Forward sound wave trailing edge (ft), Backward shock wave (b), Backward sound wave leading edge (bl), Backward sound wave trailing edge (bt), Contact (c), Neumann boundary (n), Dirichlet boundary (d), Periodic boundary (p), Passive boundary (pa), Unspecified wave type (u, default), Enter the wave type for the inner ellipse (normal/forward direction = outward): c wave type = CONTACT Enter the surface tension for the inner elliptical curve (dflt = 0): 0 Type y to turn off tracking for the inner ellipse: n Enter the component number for the region inside the inner ellipse (default = 2): The 1-st ellipse can be chosen to as closed (full) (C, F, default), half (H), quarter (Q), or angle sector specified (A). Enter choice: Q This ellipse has been oriented in the counter clockwise direction so that the forward moving normal directions point towards the outside of the ellipse. Enter the 2 coordinates of the center of the 1-st ellipse (dflt = ( 0, 0)): 0 0 Enter the common radius or the 2 mean radii of the 1-st ellipse (dflt = ( 0.080625, 0.080625)): 0.11 0.11 Enter a scaling factor for the radius (dflt = 1): Type 'y' to rotate the 1-st ellipse: Type 'y' if the 1-st ellipse is to be fourier polynomial perturbed: n Type 'y' if the 1-st ellipse is to be Legendre polynomial perturbed: n Possible choices for wave types are Forward shock wave (f), Forward sound wave leading edge (fl), Forward sound wave trailing edge (ft), Backward shock wave (b), Backward sound wave leading edge (bl), Backward sound wave trailing edge (bt), Contact (c), Neumann boundary (n), Dirichlet boundary (d), Periodic boundary (p), Passive boundary (pa), Unspecified wave type (u, default), Enter the wave type for the 1-st ellipse (normal/forward direction = outward): b wave type = BACKWARD_SHOCK_WAVE Type y to turn off tracking for the 1-st ellipse: y Enter the component number for the region between the 0-th and the 1-st ellipses (default = 3): Enter the component number for the region outside the 1-st ellipse (default = 4): To have upper plane wave, enter the height above L[1]: To have lower plane wave, enter the height above L[1]: You will now be prompted for the number of different equations of state models, followed by prompts for the parameters of each EOS. The various equations of state will then be referred to by the integer that corresponds to the order in which they are prompted. Enter the number of EOS models to be used: 1 Enter the equation of state type for the material with index 0. Current choices are Obstacle (behind reflecting wall) (O), Polytropic (gamma law) gas (P), Stiffened polytropic gas (SP), Multiple component polytropic gas (MP), Sesame table lookup (SE), JWL Equation of state (J), Mie Gruneisen (M), Isentropic two phase eos (S2PH), or Generic Test (minimal SPOLY for GENERIC testing) (GT). Enter choice here (dflt = P): P Enter the ratio of specific heats (gamma), the ideal gas constant (R, PV = RT, default for R = 1) the shear viscosity coefficient (default = 0), the the bulk viscosity coefficient (default = 0), and the the thermal conductivity (default = 0) for the gas with index 0: 1.666 Use current defaults for artificial viscosity parameters (dflt = y, type p to print defaults): Use defaults for thermodynamic restrictions (dflt = yes): Initialized params Gas_param = 1 Equation of state = 11 POLYTROPIC gamma = 1.666, R = 1 bulk viscosity = 0, shear viscosity = 0 Artificial Viscosities and Heat Conductions Use nonlinear artificial viscosity = yes Coefficient of nonlinear artificial viscosity = 0.5 Use linear artificial viscosity = yes Coefficient of linear artificial viscosity = 0.05 Use upwind artificial viscosity = no Coefficient of upwind artificial viscosity = 0 Use MUSCL slope flattening = yes Muscl slope flattening parameter eta inverse = 2 Muscl slope flattening minimum shock strength = 0.25 Muscl slope flattening minimum shock specific volume jump = 1e-06 Coefficient of artificial heat conduction = 0 Muscl slope flattening charateristic speed cutoff = 0 Coefficient of dynamic surface tension = 0 Artificial viscosity stability coefficient = 1.280776406404415 composition_type = 0 PURE_NON_REACTIVE min_energy = 2.220446049250313e-16 min_pressure = 2.220446049250313e-16 vacuum_dens = 2.220446049250313e-16 raref_press = 0.9999999999999998 Init states from inner to outer (y or n, default = y): y Initialize the elliptical state inside the inner ellipse - The choices for flow initialization types are Ambient region (constant initial conditions) (A) Tabulated region (initial conditions read from a file) (TR) Random perturbation of an ambient region (RA) Elliptical region (constant thermodynamics, radial velocity) (E, default) Random perturbation of elliptical region (RE) Overlay of one dimensional flow (O) Enter choice: E Enter the gas state inside the inner ellipse rho, pr, radial velocity: 1.37 0.000156 0 Initialized elliptical state: address 0x85d82f8 density = 1.37 specific internal energy = 0.0001709738936016308 pressure = 0.000156 sound speed = 0.01377334779525483 temperature = 0.0001138686131386861 specific entropy = -13.94914300794029 mx = 0 vx = 0 my = 0 vy = 0 total energy = 0.0002342342342342342 Mach number = 0 velocity angle = 0 State type = TGAS_STATE Params state = 2 Gas_param = 2 Type 'y' to enter an arbitrary state for the region between the 0-th and the 1-st ellipses: n The choices for flow initialization types are Ambient region (constant initial conditions) (A) Tabulated region (initial conditions read from a file) (TR) Random perturbation of an ambient region (RA) Elliptical region (constant thermodynamics, radial velocity) (E, default) Random perturbation of elliptical region (RE) Overlay of one dimensional flow (O) Enter choice: E Type 'y' to input an independent behind state: n In addition to the ahead state, two more parameters are needed to specify the contact configuration. Enter the density behind the contact: 0.1 Enter the velocity jump (shear) across the contact (dflt = 0): Initialized elliptical state: address 0x85d85c0 density = 0.1 specific internal energy = 0.002342342342342342 pressure = 0.000156 sound speed = 0.05097999607689274 temperature = 0.00156 specific entropy = -7.401723402082069 mx = 0 vx = 0 my = 0 vy = 0 total energy = 0.0002342342342342342 Mach number = 0 velocity angle = 0 State type = TGAS_STATE Params state = 1 Gas_param = 1 Type 'y' to enter an arbitrary state for the region outside the outer ellipse: n The choices for flow initialization types are Ambient region (constant initial conditions) (A) Tabulated region (initial conditions read from a file) (TR) Random perturbation of an ambient region (RA) Elliptical region (constant thermodynamics, radial velocity) (E, default) Random perturbation of elliptical region (RE) Overlay of one dimensional flow (O) Enter choice: E Type 'y' to input an independent behind state: n In addition to the ahead state, one more parameter is needed to specify the shock wave configuration. The choices are The pressure behind the shock (P) The magnitude of the normal component of the gas velocity behind the shock (V) The normal component of the shock speed (S) The normal shock Mach number (M) Enter choice here: M Enter the normal shock Mach number of the wave: 10 Initialized elliptical state: address 0x85d8930 density = 0.3886297376093297 specific internal energy = 0.07517791057423524 pressure = 0.01945810352588147 sound speed = 0.2888149957067779 temperature = 0.05006848844244067 specific entropy = -3.550912984196986 mx = -0.1471434289099821 vx = -0.3786211261524668 my = 0 vy = 0 total energy = 0.05707217704039912 Mach number = 1.310946909892675 velocity angle = 180.0000000000002 State type = TGAS_STATE Params state = 1 Gas_param = 1 Enter the boundary type -- Unknown, Periodic, Reflecting, Mixed, Neumann, No Slip Neumann, Dirichlet, or Passive -- for the left boundary in the x direction: Reflecting Enter the boundary type -- Unknown, Periodic, Reflecting, Mixed, Neumann, No Slip Neumann, Dirichlet, or Passive -- for the right boundary in the x direction: Dirichlet Specify the Dirichlet boundary state time-independent boundary state set by ambient state (A, default), flow-through boundary conditions (FT), flow-through boundary with constant pressure (FP), or flow-through boundary with time dependent pressure (FD), or random velocity inlet (R), preset boundary state (P), or time-independent state specified by the user (U). Enter choice here: FT Enter the boundary type -- Unknown, Periodic, Reflecting, Mixed, Neumann, No Slip Neumann, Dirichlet, or Passive -- for the lower boundary in the y direction: Reflecting Enter the boundary type -- Unknown, Periodic, Reflecting, Mixed, Neumann, No Slip Neumann, Dirichlet, or Passive -- for the upper boundary in the y direction: Dirichlet Specify the Dirichlet boundary state time-independent boundary state set by ambient state (A, default), flow-through boundary conditions (FT), flow-through boundary with constant pressure (FP), or flow-through boundary with time dependent pressure (FD), or random velocity inlet (R), preset boundary state (P), or time-independent state specified by the user (U). Enter choice here: FT Specify parabolic steps Type 'y' to have the Navier-Stokes terms computed for several eos models, and this will turn on parabolic driver parab_driver (y, n(dflt)): N Type 'y' to have local mesh refinement : N Statistics Control Type 'y' to initiate prompting for statistics/diagnostics: y Type 'y' to request grid statistics for conserved variables: n Type 'y' to request a periodic glimpse of the solution via a plot of the component regions: n Type 'y' to obtain cross sectional plots: n Type 'y' to request interface extrema data: y Specify the interval type for printing [exact (default), constant, mesh]: exact Enter the time interval and first time for printing (default = 0.02 0): 0.02 0 Request binary/non-binary output [b,n(dflt)]: n Compute interface extrema for planar ('p', default) or radial geometry ('r'): r Enter the coordinates of the origin (default = ( 0, 0)): 0 0 Enter a sub-grid refinement factor for the averaging of the ambient state at the interface extrema (default = 2): 2 Current gas param list Number of params = 2 Param[0] Gas_param = 1 Equation of state = 11 POLYTROPIC gamma = 1.666, R = 1 bulk viscosity = 0, shear viscosity = 0 Artificial Viscosities and Heat Conductions Use nonlinear artificial viscosity = yes Coefficient of nonlinear artificial viscosity = 0.5 Use linear artificial viscosity = yes Coefficient of linear artificial viscosity = 0.05 Use upwind artificial viscosity = no Coefficient of upwind artificial viscosity = 0 Use MUSCL slope flattening = yes Muscl slope flattening parameter eta inverse = 2 Muscl slope flattening minimum shock strength = 0.25 Muscl slope flattening minimum shock specific volume jump = 1e-06 Coefficient of artificial heat conduction = 0 Muscl slope flattening charateristic speed cutoff = 0 Coefficient of dynamic surface tension = 0 Artificial viscosity stability coefficient = 1.280776406404415 composition_type = 0 PURE_NON_REACTIVE min_energy = 2.220446049250313e-16 min_pressure = 2.220446049250313e-16 vacuum_dens = 2.220446049250313e-16 raref_press = 0.9999999999999998 Param[1] Gas_param = 2 Equation of state = 11 POLYTROPIC gamma = 1.666, R = 1 bulk viscosity = 0, shear viscosity = 0 Artificial Viscosities and Heat Conductions Use nonlinear artificial viscosity = yes Coefficient of nonlinear artificial viscosity = 0.5 Use linear artificial viscosity = yes Coefficient of linear artificial viscosity = 0.05 Use upwind artificial viscosity = no Coefficient of upwind artificial viscosity = 0 Use MUSCL slope flattening = yes Muscl slope flattening parameter eta inverse = 2 Muscl slope flattening minimum shock strength = 0.25 Muscl slope flattening minimum shock specific volume jump = 1e-06 Coefficient of artificial heat conduction = 0 Muscl slope flattening charateristic speed cutoff = 0 Coefficient of dynamic surface tension = 0 Artificial viscosity stability coefficient = 1.280776406404415 composition_type = 0 PURE_NON_REACTIVE min_energy = 2.220446049250313e-16 min_pressure = 2.220446049250313e-16 vacuum_dens = 2.220446049250313e-16 raref_press = 0.9999999999999998 Enter the EOS indices of the inner and outer materials, respectively: 1 0 Enter a file name for the output for interface minimum (default = out/Implosion/intfc_extrema/impf16RM.min): Enter a file name for the output for interface maximum (default = out/Implosion/intfc_extrema/impf16RM.max): Enter a file name for the output for interface amplitude (default = out/Implosion/intfc_extrema/impf16RM.amp): Type 'y' to get data for 1%-99% levels: n Type 'y' to get data for 5%-95% levels: n Type 'y' to request layer statistics: n Type 'y' to request interface statistics: n Type 'y' to request rect state statistics: n Type 'y' to request printing of the radial amplitude and velocity data: n Type 'y' to request printing of front states along the contacts: n Type 'y' to request printing of states for all fronts: n --- End of Input ---