
[Sponsors] 
October 15, 2013, 09:39 
multispecies plasma test case

#1 
Member
Tom Jentink
Join Date: Jan 2013
Posts: 48
Rep Power: 5 
I've just installed the newest 2.0.8 from github.
I'm trying to run the multispecies reactargon case, and am getting a lot of errors. I'm using the config file from the newest release. see below:  SU2_CFD plasma_argon_3D.cfg WARNING: unrecognized option in the config. file: MIXTURE_MOLAR_MASS. WARNING: unrecognized option in the config. file: BLOCK_DIAGONAL_JACOBIAN. WARNING: unrecognized option in the config. file: MG_RESTART_CYCLE. WARNING: unrecognized option in the config. file: NUMBER_OF_SPECIES_IN_PLASMA. WARNING: unrecognized option in the config. file: NUMBER_OF_CHEMICAL_REACTIONS. WARNING: unrecognized option in the config. file: PARTICLE_MASS. WARNING: unrecognized option in the config. file: MAGNETIC_DIPOLE. WARNING: unrecognized option in the config. file: SOUR_NUM_METHOD_ELEC. WARNING: unrecognized option in the config. file: VISC_NUM_METHOD_ELEC. WARNING: unrecognized option in the config. file: CYCLE_ADAPT. WARNING: unrecognized option in the config. file: HORIZONTAL_PLANE_POSITION. WARNING: unrecognized option in the config. file: HORIZONTAL_PLANE_MARKER. WARNING: unrecognized option in the config. file: ONESHOT_CAUCHY_EPS. WARNING: No species CFL numbers specified for plasma multitimestepping, using mean flow CFL parameters WARNING: No species CFL ramp iteration specified for plasma multitimestepping, using mean flow CFL parameters WARNING: No maximum CFL specified for plasma multitimestepping, using mean flow CFL parameters   _____ _ _ ___   / ____     __ \ Web: su2.stanford.edu    (___     )  Twitter: @su2code   \___ \     / / Forum: www.cfdonline.com/Forums/su2/   ____)   __  / /_   _____/ \____/ ____ Suite (Computational Fluid Dynamic Code)   Release 2.0.8    Stanford University Unstructured (SU2).   Copyright (C) 20122013 Aerospace Design Laboratory (ADL).   SU2 is distributed in the hope that it will be useful,   but WITHOUT ANY WARRANTY; without even the implied warranty of   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU   Lesser General Public License (version 2.1) for more details.    Physical case definition  Plasma equations (with viscosity). Using 3 species Argon gas model. Mach number: 4.6. Angle of attack (AoA): 0 deg, and angle of sideslip (AoS): 0 deg. No restart solution, use the values at infinity (freestream). The reference length/area (force coefficient) is 1. The reference length (moment computation) is 1. Reference origin (moment computation) is (0.25, 0, 0). Surface(s) where the force coefficients are evaluated: BulletHead. Surface(s) where the objective function is evaluated: Surface(s) plotted in the output file: BulletHead. Surface(s) affected by the design variables: Input mesh file name: mesh_plasma_argon_3D.su2  Space numerical integration  1st order Roe solver for the plasma inviscid terms. Average of gradients (1st order) for computation of viscous flow terms. Piecewise constant integration of the plasma source terms. Gradient computation using GreenGauss theorem.  Time numerical integration  Local time stepping (steady state simulation). Euler implicit method for the plasma equations. CFL ramp definition. factor: 1.1, every 100 iterations, with a limit of 8. CourantFriedrichsLewy number: 0.5  Convergence criteria  Maximum number of iterations: 5. Reduce the density residual 6 orders of magnitude. The minimum bound for the density residual is 10^(22). Start convergence criteria at iteration 10.  Output information  Writing a flow solution every 250 iterations. Writing the convergence history every 1 iterations. The output file format is Tecplot ASCII (.dat). Convergence history file name: history. Surface flow coefficients file name: surface. Flow variables file name: solution. Restart flow file name: restart.dat.  Config file boundary information  Farfield boundary marker(s): Inlet, Exit. Symmetry plane boundary marker(s): Sidewall1, Sidewall2. Constant heat flux wall boundary marker(s): BulletHead.  Flow & Nondimensionalization information  Viscous flow: Computing pressure using the ideal gas law based on the freestream temperature and a density computed from the Reynolds number. Input conditions: Grid conversion factor to meters: 1 Ratio of specific heats: 1.66667 Specific gas constant (J/(kg.K)): 208 Freestream pressure (N/m^2): 0 Freestream temperature (K): 810 Freestream density (kg/m^3): 0 Freestream velocity (m/s): (2437.57,0,0) Freestream velocity magnitude (m/s):2437.57 Freestream energy (kg.m/s^2): nan Freestream viscosity (N.s/m^2): 3.66516e05 Reference values: Reference pressure (N/m^2): 1 Reference temperature (K): 1 Reference energy (kg.m/s^2): nan Reference density (kg/m^3): 1 Reference velocity (m/s): 1 Reference viscosity (N.s/m^2): 1 Resulting nondimensional state: Mach number (nondimensional): 4.6 Reynolds number (nondimensional): 0 Reynolds length (m): 1 Froude number (nondimensional): 778.388 Lenght of the baseline wave (nondimensional): 3.8069e+06 Negative pressure, temperature or density is not allowed! Specific gas constant (nondimensional): 208 Freestream temperature (nondimensional): 810 Freestream pressure (nondimensional): 0 Freestream density (nondimensional): 0 Freestream velocity (nondimensional): (2437.57,0,0) Freestream velocity magnitude (nondimensional): 2437.57 Freestream energy (nondimensional): nan Freestream viscosity (nondimensional): 3.66516e05 Force coefficients computed using freestream values.  Read grid file information  Three dimensional problem. 46327 interior elements. 31600 tetrahedra. 13104 hexahedra. 1623 pyramids. 22557 points. 6 surface markers. 315 boundary elements in index 0 (Marker = BulletHead). 238 boundary elements in index 1 (Marker = Cylinder). The configuration file doesn't have any definition for marker Cylinder!! 

October 15, 2013, 16:15 

#2 
Member
Sean R. Copeland
Join Date: Jan 2013
Posts: 40
Rep Power: 5 
Hi Tom,
The reactinggas solver is currently undergoing a complete overhaul and, unfortunately, the test case in the repository and in the user guide is outofdate. I apologize for the inconvenience, I've temporarily removed it from both locations until I can appropriately document the changes and set up a proper testcase for the community. Thanks for bringing it to my attention. The new implementation will be finished in a month or two and I encourage you to stay tuned for updates. Thanks. Regards, Sean 

October 15, 2013, 16:34 

#3 
Member
Tom Jentink
Join Date: Jan 2013
Posts: 48
Rep Power: 5 
I actually got the first part of it to run by checking and rechecking the input items and config_structure.cpp and taking out the wrong and discontinued items.
I started it euler, and it runs and gives me a solution that looks proper. However, no matter what I try, I get nan's when restarting. Whether I run 25 iterations and restart or whether I run 15,000 and restart. 

October 15, 2013, 16:46 

#4 
Member
Sean R. Copeland
Join Date: Jan 2013
Posts: 40
Rep Power: 5 
Hi Tom,
I'm glad to hear that you've seen some life in the reactinggas solver. There are likely settings in the configuration file that are referencing defunct and obsolete subroutines within the old reacting gas solver that are contributing to your difficulties. The old code will be removed for the next interim release. This may also be contributing to your restartfile issues. Francisco and I sent you a mesh and configuration file the other day, I would encourage you to try running that case as a starting point. Be sure you are referencing the TNE2_EULER solver (NOT the PLASMA_EULER solver), and set the GAS_MODEL=N2. Again, I want to stress that the reactinggas solver is not supported at this time. Many of the settings are not welldocumented. If you continue to have difficulty, please send me an email and I will do my best to assist. Thanks, Sean 

October 15, 2013, 16:52 

#5 
Member
Tom Jentink
Join Date: Jan 2013
Posts: 48
Rep Power: 5 
Yes, I planned on continuing with that tne case you sent earlier... thanks for that.
I'll hold off on the plasma case for now. 

May 6, 2016, 22:17 

#6 
New Member
California
Join Date: Apr 2016
Posts: 2
Rep Power: 0 
Removing this reply to post elsewhere.
Last edited by abmason2; May 13, 2016 at 00:18. Reason: It should be converted to a new post. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Error when run test case turb_ONERAM  mistuning  SU2  1  October 3, 2013 01:17 
Case (VOF + Solidification&Melting) will not work with Species Transport? ERROR ?  andy.eco  FLUENT  4  September 14, 2010 17:27 
OF 1.6  Ubuntu 9.10 (64bit)  GLIBCXX_3.4.11 not found  piprus  OpenFOAM Installation  22  February 25, 2010 14:43 
Porous Media test case  Alex  FLUENT  0  April 9, 2006 08:23 
c1 body test case  Eric Lenormand  Main CFD Forum  0  March 2, 2000 07:54 