Ticket #117 (new defect - convergence)

Opened 2 years ago

initial solvers fail badly in atomic region with SNR abundances

Reported by: gary Owned by: nobody
Priority: major Milestone: C12 branch
Component: chemical network Version: trunk
Keywords: Cc:

Description

attached sim fails with non-conservation in C08 and trunk r3667. Robin reports does fine in newmole. This might be an interesting sim to add to the test suite after newmole comes to trunk

CoolHeatError?: Te: 3.0133e+03 C: 2.3546e-16 H: 1.8927e-17 (C-H)/H: 1.1440e+01 CoolHeatError?: Te: 2.9530e+03 C: 2.2965e-16 H: 1.8987e-17 (C-H)/H: 1.1095e+01 CoolHeatError?: Te: 2.8939e+03 C: 2.2428e-16 H: 1.9051e-17 (C-H)/H: 1.0773e+01

PROBLEM DISASTER PressureTotal?: the chemical species are not conserved. The sum of the densities of the atoms and ions for Carbon is 8.146e+02 which is greater than the gas-phase abundance of the element, 1.245e+02. Something that cannot happen, has happened. This is TotalInsanity?, I live in ../service.cpp.

Attachments

savin.in (8.3 kB) - added by gary 2 years ago.
input script jon posted on yahoo site

Change History

Changed 2 years ago by gary

input script jon posted on yahoo site

Note: See TracTickets for help on using tickets.