lanl / lanl/singularity-eos

Mitigation strategies for Maxwell-constructed EOS in the PTE Solvers

Open
#613 1 comment 0 reactions 1 assignee View on GitHub

@Yurlungur is already working on this.

Since Mar 12, 2026.

bug
Dominant language
C++
Stars
39
Forks
22
Avg merge
10h 22m
Merged PRs (30d)
3

Description

A user reported issues with Maxwell constructed EOS. The primary problem is that in a Maxwell constructed region, P = 0 for a wide range of densities and temperatures. We have a mitigation for densities, by having a minimum volume fraction step. But for temperatures, we do not, which can make the Jacobian ill conditioned. We should mitigate in two ways:

  1. Adding a minimum gradient size in temperature
  2. Adding a maximum temperature, to avoid temperatures diverging to NaN.

Contributor guide

No contributing guide indexed for this repository

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.