Skip to content

Conversation

@JRChreim
Copy link

No description provided.

henryleberre and others added 30 commits May 23, 2022 21:42
Adding cylindrical coordinates to the restart of a simulation.
m_phase_change: removing m_palpha_eps and ptgalpha_eps from global parameters, and adding it to user-input level

m_start_up: adding palpha_eps and ptgalpha_eps to the startup dictionary

m_global_parametes: setting palpha_eps and ptgalpha_eps as default and setting theis default values.

m_mpi_proxy: broadcasting palpha_eps and ptgalpha_eps
Changes in both pT-equilibrium and pTg-equilibrium. Note that, while the pT-equilibrium algorithm is globally convergent given the physical constraints are satisfied, this is yet not true for the pTg-equilibrium. Currently debugging the case in which there is a pressure jump.
- printing some variables of interest for debugging
- using Newton's solver instead of fixed-point iteration
- updating state variables for the pTg-equilibrium
- correcting Jacobian matrix and Residual array for the Newton Solver
- cleaned unused variables
- pT-equilibrium as IC for the pTg-equilibrium, as it improves stability of the method
- removed unnecessary subroutines as they have already been called in th pT-equilibrium (as part of the new IC)
- renamed a few variables for consistency
- removed many repeated printing statements and substituted for a new subroutine that prints whatever is needed (temporary)
- Added Mass depletion test to check whether phase-change can happen at the current thermodynamic state and mass conditions:
test zero liquid phase
test zero vapor phase
decide if either one is a solution
if yes, continue with this solution. If not, start phase-change
- added function to print output variables so error can be tracked easily
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

5 participants