Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Nlte and p1 model rebase #239

Draft
wants to merge 18 commits into
base: main
Choose a base branch
from
Draft

Nlte and p1 model rebase #239

wants to merge 18 commits into from

Commits on Oct 19, 2023

  1. Configuration menu
    Copy the full SHA
    3c9e4c0 View commit details
    Browse the repository at this point in the history
  2. Set up the nlte case.

    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    b7adbbd View commit details
    Browse the repository at this point in the history
  3. Configuration menu
    Copy the full SHA
    d952f01 View commit details
    Browse the repository at this point in the history
  4. Fix issue with negative unphisical number densities. Negative values …

    …odduc due to interpolation error so we force this value to zero in face_integrator.cpp and BCintegrator.cpp. We also check for negative values at the end of the iteration with Check_undershoot. undershoot can occur due to discretisation error in the solution of the transport equation. The use of a TVD scheme could mitigate this effect.
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    ba1057d View commit details
    Browse the repository at this point in the history
  5. An attempt to run the two-temperature plasma model but it failed. The…

    … electron number density approaches zero at remote regions and Te goes to inf. We need to think how to address this one.
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    8c3441d View commit details
    Browse the repository at this point in the history
  6. Added a input-file flag to enable saving the temperature field. This …

    …field is to be read by the P1 stand-alone code.
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    072dbeb View commit details
    Browse the repository at this point in the history
  7. The coupling of the P1 radiation model is complete. The input file ha…

    …s been modified to add options for the P1 model. some h5 files related to P1 have been added to test/rad-data
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    6d8a75a View commit details
    Browse the repository at this point in the history
  8. I have added a class RadiativeDecay to deal with radiative decay reac…

    …tions. It is a low level class. Currently it can huddle specific reactions for our specific chemistry model like Ar(4p) -> Ar(r), Ar(4p) -> Ar(m), and Ar(r) -> Ar(g). It calculates the effective Einstein A coefficients for these reactions based on the A coef of the transitions involved with these lumped species. For the number densities of the internal states, we just assume that they are bolzmann distributed at a given Te evaluated by the model. The method also calculates the escape factors and consider them in the evaluation of the effective A coef.
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    de27f14 View commit details
    Browse the repository at this point in the history
  9. Configuration menu
    Copy the full SHA
    5e96715 View commit details
    Browse the repository at this point in the history
  10. Configuration menu
    Copy the full SHA
    a9585d0 View commit details
    Browse the repository at this point in the history
  11. Fixed issue with the order of the species. I have added the functions…

    … SetStateIndices and SetSpeciesIndices in class GasMixture to set an idenx system for the sate variables and species variables. These indeces are used in GasMixture and derived classes.
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    274a56b View commit details
    Browse the repository at this point in the history
  12. I made some changes to test-table.cpp to make it compatible with the …

    …changes in h5ReadTable. h5ReadTable now gets groupName as an extra argument.
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    3c8bc24 View commit details
    Browse the repository at this point in the history
  13. Configuration menu
    Copy the full SHA
    5de5aa6 View commit details
    Browse the repository at this point in the history
  14. Fixed an if statement in M2ulPhyS.cpp related to the selection of the…

    … radiation_ class that was causing some tests to fail.
    malamast committed Oct 19, 2023
    Configuration menu
    Copy the full SHA
    2895528 View commit details
    Browse the repository at this point in the history

Commits on Oct 20, 2023

  1. fixed some issues in source_term related with the selection of the ra…

    …diation model that caused some tests to fail.
    malamast committed Oct 20, 2023
    Configuration menu
    Copy the full SHA
    f5b0c59 View commit details
    Browse the repository at this point in the history
  2. Configuration menu
    Copy the full SHA
    00f038b View commit details
    Browse the repository at this point in the history
  3. moved the rad-data for nec in rad-data/NEC/ and chenged the paths in …

    …the plasma.lte1/2d.ini input files.
    malamast committed Oct 20, 2023
    Configuration menu
    Copy the full SHA
    5208f37 View commit details
    Browse the repository at this point in the history
  4. Configuration menu
    Copy the full SHA
    19f57c4 View commit details
    Browse the repository at this point in the history