TCNEQ Version - Revision history

Conrad54418: /* Boundary Conditions */

2024-09-06T17:07:44Z

‎Boundary Conditions

Conrad54418: /* Initial Conditions */

2024-09-06T17:07:34Z

‎Initial Conditions

Conrad54418: /* Boundary Conditions */

2024-09-06T17:02:52Z

‎Boundary Conditions

Conrad54418: /* Example of solver.inp file: */

2024-08-29T14:07:33Z

‎Example of solver.inp file:

Conrad54418: /* Example of solver.inp file: */

2024-08-29T14:06:47Z

‎Example of solver.inp file:

Conrad54418: /* Post-Processing */

2024-08-29T14:05:27Z

‎Post-Processing

Conrad54418: /* Example of flow.pht file */

2024-08-29T14:04:30Z

‎Example of flow.pht file

Conrad54418: /* Simulation Inputs */

2024-08-29T14:01:51Z

‎Simulation Inputs

Conrad54418: /* Example of solver.inp file: */

2024-08-29T14:01:27Z

‎Example of solver.inp file:

Conrad54418: Undo revision 2084 by Conrad54418 (talk)

2024-08-29T14:01:13Z

Undo revision 2084 by Conrad54418 (talk)

@@ Line 26: / Line 26: @@
 * pressure - Specification of static pressure over a surface, [Pa]
 ** Used to compute mole fractions of each species of the gas with Dalton's Law of partial pressures in conjunction with reference mole fractions specified in solver.inp
 * heat flux - set to zero for adiabatic wall boundary condition

@@ Line 32: / Line 32: @@
 * initial velocity - Components and magnitude of flow velocity, [m/s]
 ** If a supersonic outlet condition is used, set such that flow is initialized Mach > 1
-* initial temperature - Value used to set translational-rotational temperature, [K]
+* initial temperature - Value used to set translational-rotational temperature (related to initial vibrational temperature, see input.config in source code), [K]
-* initial scalar_1 - Initial value of turbulent eddy viscosity
+* initial scalar_1 - Initial value of turbulent eddy viscosity, [Pa-s]
 * initial pressure - Static pressure of the gas, [Pa]
 ** For multi-species flows, this value is used in combination with solver.inp values to compute the mole fraction of each species

@@ Line 21: / Line 21: @@
 Below are the recognized boundary conditions that can be applied for the current version:
 * comp1/comp2/comp3 - Specification of one/two/three components of velocity, [m/s]
-* temperature - Specification of translational-rotational temperature, [K]. By default, vibrational temperature is held in equilibrium with this value and nonequilibrium is controlled through simulation inputs.
+* temperature - Specification of translational-rotational temperature, [K]. By default, vibrational temperature is held in equilibrium with this value and nonequilibrium is controlled through simulation inputs (see input.config in source code)
 * surfID - When value is set to 702, the boundary is treated as a slip wall. If using this option, include a boundary layer mesh along the surface to ensure the wall normal direction is accurately computed.
-* scalar_1 - Turbulent eddy viscosity
+* scalar_1 - Turbulent eddy viscosity, [Pa-s]
 * pressure - Specification of static pressure over a surface, [Pa]
 ** Used to compute mole fractions of each species of the gas with Dalton's Law of partial pressures in conjunction with reference mole fractions specified in solver.inp

@@ Line 135: / Line 135: @@
 #STEP SEQUENCE
 #{
-        Step Construction  : 0 1 0 1 # 10 11 10 11
+        Step Construction  : 0 1 0 1 # laminar
 #}</nowiki>

@@ Line 50: / Line 50: @@
       Number of Timesteps: 1000
       Time Step Size: 1e-8
 #       Limit instructions : switch min max -- change switch from zero to activate
@@ Line 134: / Line 135: @@
 #STEP SEQUENCE
 #{
-        Step Construction  : 0 1 0 1
+        Step Construction  : 0 1 0 1 # 10 11 10 11
 #}</nowiki>

@@ Line 139: / Line 139: @@
 == Post-Processing ==
-An example of the ''flow.pht'' file is provided to demonstrate the ordering of the variables that can be viewed in Paraview. Note that both the ''errors'' and ''DCqpt'' fields would need to be output to the restart file for the example below to work. Saving these two fields is controlled through the solver.inp options.
+An example of the ''flow.pht'' file is provided to demonstrate the ordering of the variables that can be viewed in Paraview. Note that both the ''evisc'' and ''dwal'' fields would require turbulence modeling to be turned on for the example below to work. This is controlled through the solver.inp options.
 === Example of flow.pht file ===

@@ Line 144: / Line 144: @@
   <nowiki><?xml version="1.0" ?>
-<PhastaMetaFile number_of_pieces="24">
+<PhastaMetaFile number_of_pieces="40">
-    <GeometryFileNamePattern pattern="24-procs_case/geombc.dat.%d"
+    <GeometryFileNamePattern pattern="ent_from_NAS/geombc.dat.%d"
                              has_piece_entry="1"
                              has_time_entry="0"/>
-    <FieldFileNamePattern pattern="24-procs_case/restart.%d.%d"
+    <FieldFileNamePattern pattern="ent_from_NAS/restart.%d.%d"
                           has_piece_entry="1"
                           has_time_entry="1"/>
-    <TimeSteps number_of_steps="1"
+    <TimeSteps number_of_steps="9"
 	      auto_generate_indices="1"
-               start_index="2010"
+               start_index="65800"
-	      increment_index_by="50"
+	      increment_index_by="1000"
                start_value="0"
                increment_value_by="1">
     </TimeSteps>
-    <Fields number_of_fields="9">
+    <Fields number_of_fields="10">
       <Field phasta_field_tag="solution"
-             paraview_field_tag="rho_N2"
+             paraview_field_tag="rho_1"
              start_index_in_phasta_array="0"
              number_of_components="1"
              data_dependency="0"
              data_type="double"/>
       <Field phasta_field_tag="solution"
-             paraview_field_tag="rho_O2"
+             paraview_field_tag="rho_5"
-             start_index_in_phasta_array="1"
+             start_index_in_phasta_array="4"
              number_of_components="1"
              data_dependency="0"
@@ Line 173: / Line 191: @@
       <Field phasta_field_tag="solution"
              paraview_field_tag="velocity"
-             start_index_in_phasta_array="2"
+             start_index_in_phasta_array="5"
              number_of_components="3"
              data_dependency="0"
@@ Line 179: / Line 197: @@
       <Field phasta_field_tag="solution"
              paraview_field_tag="temp_vib"
-             start_index_in_phasta_array="5"
+             start_index_in_phasta_array="8"
              number_of_components="1"
              data_dependency="0"
              data_type="double"/>
       <Field phasta_field_tag="solution"
-             paraview_field_tag="temperature"
+             paraview_field_tag="temp"
-             start_index_in_phasta_array="6"
+             start_index_in_phasta_array="9"
              number_of_components="1"
              data_dependency="0"
-            data_type="double"/>
+	    data_type="double"/>
-      <Field phasta_field_tag="errors"
+      <Field phasta_field_tag="solution"
-             paraview_field_tag="nu"
+             paraview_field_tag="evisc"
-             start_index_in_phasta_array="9"
+             start_index_in_phasta_array="10"
              number_of_components="1"
              data_dependency="0"
              data_type="double"/>
-      <Field phasta_field_tag="DCqpt"
+      <Field phasta_field_tag="dwal"
-             paraview_field_tag="DCqpt1"
+             paraview_field_tag="dwal"
              start_index_in_phasta_array="0"
              number_of_components="1"
-             data_dependency="1"
+             data_dependency="0"
              data_type="double"/>
    </Fields>
 </PhastaMetaFile></nowiki>

@@ Line 39: / Line 39: @@
 == Simulation Inputs ==
-Below is an example of the input script for the current version of the code. Capability is included for handling multi-species flows up to number of species ''nsp'' equal to 5.
+Below is an example of the input script for the current version of the code. Capability is included for handling multi-species flows up to number of species (''nsp'') equal to 5.
 === Example of solver.inp file: ===

@@ Line 43: / Line 43: @@
 === Example of solver.inp file: ===
-  <nowiki># PHASTA_HYP Version 1 Input File
+  <nowiki>
 #SOLUTION CONTROL

← Older revision		Revision as of 14:01, 29 August 2024
Line 42:		Line 42:

	=== Example of solver.inp file: ===		=== Example of solver.inp file: ===
		+
		+	<nowiki># PHASTA_HYP Version 1 Input File

	#SOLUTION CONTROL		#SOLUTION CONTROL