Opened 7 years ago

Last modified 6 years ago

#196 new defect - convergence

constant pressure model aborts due to excessive pressure convergence failures

Reported by: peter Owned by: nobody
Priority: major Milestone:
Component: pressure convergence Version: trunk
Keywords: Cc:


The following problem was reported by Alexei Baskin:


I am trying to run Cloudy for different AGN models, changing the incident continuum (nuF_nu) and hydrogen density of the slab, while holding its total pressure constant. The code crashes for log n(H)=10 and log nuFnu(1Ryd)=10 or 11 (the log of stop column is 23 in both cases). The code crashes due to convergence failures (the output is included below). Resetting the limit of failures to 100 does not help. In addition, assuming 10-times solar metallicity causes the same problem for a much wider range of assumed parameters.

Thank you, Alexei

The test case failed with c10_branch r4666. The reporter tested it against the c08.01 release.

Attachments (1) (302 bytes) - added by peter 7 years ago.
test case that fails

Download all attachments as: .zip

Change History (4)

Changed 7 years ago by peter

Attachment: added

test case that fails

comment:1 Changed 6 years ago by Ryan Porter

Component: pressure convergenceconvergence (generic)
Owner: nobody deleted
Type: defect - code abortsdefect - convergence
Version: c10_branchtrunk

Something is really ringing the code's bell here. Still running after nine hours, still waiting for the first zone to print. I'm changing this one to critical.

comment:2 Changed 6 years ago by Ryan Porter

Looking again, "print last" was issued in this sim. It's not unreasonable that a hydrostatic model iterating to convergence would take 10 or so hours. Running again to see if it finishes.

comment:3 Changed 6 years ago by Ryan Porter

Component: convergence (generic)pressure convergence
Owner: set to nobody

Still fails with multiple pressure convergence failures at r6340. It starts around zone 55 of iteration 2. But the temperature here is very close to the one in this branch of atmdat_3body.cpp:

226 else if( alogt <= 5.1139434 )

If so, this would be the same problem seen in #158, #176, #181.

Note: See TracTickets for help on using tickets.