constant pressure model aborts due to excessive pressure convergence failures (trac #196)

cloudy-bot commented 13 years ago

reported by: peter

The following problem was reported by Alexei Baskin:

Hello,

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.

Migrated from https://www.nublado.org/ticket/196

{
    "status": "new",
    "changetime": "2019-02-04T13:07:06Z",
    "_ts": "1549285626580696",
    "description": "The following problem was reported by Alexei Baskin:\n\n----------------\n\nHello,\n\nI 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.\n\nThank you,\nAlexei\n\n---------------\n\nThe test case failed with c10_branch r4666. The reporter tested it against the c08.01 release.",
    "reporter": "peter",
    "cc": "",
    "resolution": "",
    "time": "2011-02-15T10:57:31Z",
    "component": "pressure convergence",
    "summary": "constant pressure model aborts due to excessive pressure convergence failures",
    "priority": "major",
    "keywords": "",
    "version": "trunk",
    "milestone": "no milestone",
    "owner": "nobody",
    "type": "defect - convergence"
}

cloudy-bot commented 13 years ago

test case that fails Attachment: alexei.in

cloudy-bot commented 13 years ago

@ryan-porter-noaccount commented:

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.

@ryan-porter-noaccount changed component from "pressure convergence" to "convergence (generic)"
@ryan-porter-noaccount changed owner from "nobody" to ""
@ryan-porter-noaccount changed type from "defect - code aborts" to "defect - convergence"
@ryan-porter-noaccount changed version from "c10_branch" to "trunk"

cloudy-bot commented 13 years ago

@ryan-porter-noaccount commented:

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.

cloudy-bot commented 13 years ago

@ryan-porter-noaccount commented:

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.

@ryan-porter-noaccount changed component from "convergence (generic)" to "pressure convergence"
@ryan-porter-noaccount changed owner from "" to "nobody"

cloudy-bot commented 13 years ago

@peter-van-hoof-noaccount changed milestone from "" to "no milestone"

cloudy-astrophysics / bug-tracker-migration-test

constant pressure model aborts due to excessive pressure convergence failures (trac #196) #198