<div class="gmail_quote">On Wed, Sep 21, 2011 at 05:51, Gong Ding <span dir="ltr"><<a href="mailto:gdiso@ustc.edu">gdiso@ustc.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>The background of this problem is the single event upset simulation of CMOS unit circuit, i.e. SRAM, latch, flip-flop.<br>They all very big, ~ 1M mesh nodes. Direct method is too slow and memory consuming. <br>Some stupid commercial code needs more than one week to simulate one particle event.<br>
People usually need 100-1000 events to determine the circuit behavior! <br><br><span>Fortunately，the problem is in time domain. the transport equation of carrier has a d/dt term which helps stable.<br>The most difficult part comes from metal connection. Just use above inverter, when both NMOS and PMOS are closed, <br>
the metal connection region has only displacement current from semiconductor region as eps*dE/dt, while metal region<br>the current can be written as sigma*E. The current conservation here is the boundary condition at the metal-semiconductor interface. <br>
(</span><span>Since E is \frac{\partical phi}{\partical n}, here phi is the independent variable</span><span>)<br></span></div></blockquote><div><br></div><div>How much does the time step help? Given the description, I'd expect not much unless it's unreasonably short.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><span><br>Obviously, sigma >> eps/T, the governing equation of phi has a large jump parameter at the metal/semiconductor </span><span>interface</span><span>.<br>
The eigen value analysis confirmed that the number of samllest eigen value equals to the number of floating metal region.<br>And the </span><span>eigen vector is exactly the floating region.</span><span> <br></span></div>
</blockquote><div><br></div><div>Hmm, with a single low-energy mode, could you just start the Krylov method with a nonzero initial guess and make the initial guess exactly this mode. That would ensure that the Krylov space captures it exactly, and further iterations should all come from picking up the rest of the modes. Alternatively, you could formulate a modified system for which this mode was removed analytically.</div>
<div><br></div><div>Just tossing out ideas here.</div></div>