Ticket #105 (closed physics: worksforme)

Opened 2 years ago

Last modified 15 months ago

why does H2-H balance differ so much in coll_t3 between trunk and newmole

Reported by: gary Owned by: nobody
Priority: blocker Milestone: newmole onto trunk
Component: chemical network Version: trunk
Keywords: Cc:

Description

the H2 - H balance is very different for coll_t3 on the trunk and newmole. most other sims are fairly similar. In general we believe newmole more than the trunk for this type of calculation. this is a very simple sim - what process makes the balance different on the trunk and newmole?

Change History

Changed 2 years ago by rjrw

Controlling rates on newmole for the H/H2 balance are:

H,H2O=>OH,H2 src  1.879338 of 1.890086 [H   2.171103e+08,H2O    1914444]	        
H2,OH=>H2O,H snk  1.879281 of 1.890086 [H2   4.88918e+09,OH      170.6716]

-- have added commands to generate the diagnostics to coll_t3.in on newmole.

Changed 15 months ago by rjrw

  • status changed from new to closed
  • resolution set to worksforme

I ran coll_t3 for 10 iterations on the trunk, and 2 on newmole. For "C1", newmole is rock steady in space and between iterations; it looks like the trunk is heading its way.

For the other species reported in coll_t3.ovr, newmole is still stable, but the convergence a lot less striking: since the initial differences are greater, a similar level of relative drift for the trunk results means it doesn't get as far.

Other results suggest that this simulation is close to a threshold in physical parameters, so the results are likely to be sensitive to small parameter changes. Closing ticket assuming behaviour is now understood (and convergence results suggest behaviour of newmole is more reliable than trunk).

Changed 15 months ago by rjrw

  • milestone changed from C10 branch to newmole onto trunk

This wasn't a C10_branch milestone event...

Note: See TracTickets for help on using tickets.