search

pecos / tps

Torch Plasma Simulator
BSD 3-Clause "New" or "Revised" License
8 stars 2 forks source link

Multi-component reacting flow #118

Closed dreamer2368 closed 2 years ago

dreamer2368 commented 2 years ago

Main changes:

  1. Plasma condition is mainly determined by two different variables: ambipolar and twoTemperature
  2. PerfectMixture class for any plasma condition as determined above.
  3. TransportProperty computes transport property, separate from GasMixture class.
  4. Chemistry class implements mass-action law with a list of Reaction classes
  5. WorkingFluid works as preset fluid, providing a combination of GasMixture, TransportProperty and Chemistry. 5-1. supports DRY_AIR and USER_DEFINED
  6. Reaction class determines reaction rate (forward and backward) and its model and parameters can be specified from input file, reaction by reaction.
  7. WallBC: Isothermal/adiabatic wall BC extended for 2T plasma 7-1. Isothermal wall sets single temperature at the wall even for 2T plasma inside.
  8. InletBC, OutletBC: Inlet/outlet reflecting BC extended for 2T plasma 8-1. Inlet flow is set to be at single temperature.
  9. SpongeZone is extended for 2T plasma
  10. ArgonMinimalTransport supports transport property for ternary argon mixture (Ar, Ar+, e). Implemented and verified.
  11. Fluxes is extended for 2T plasma.

Minor changes:

  1. utils/l2diff: computes L2-norm difference of two solutions that are run by the same runfile.
  2. utils/beam_mesh: generates 2d rectangular mesh with various size/periodicity.
  3. utils/binaryic: generates a 2d restart file for an argon binary mixture with sinusoidal variation in mass fraction.
  4. utils/tanhic, utils/sineic: generates a 2d restart file that smoothly changes from one state to the other in one direction.

To do:

  1. Verify with a 2d manufactured solution
  2. Implement source term for electron energy equation
  3. Extend ArgonMinimalTransport to support excited states.
  4. Sheath wall BC