<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Sep 29, 2015 at 3:02 AM, Justin Chang <span dir="ltr"><<a href="mailto:jychang48@gmail.com" target="_blank">jychang48@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi all,<br>
<br>
I think I have asked a previous question related to this, and have seen other threads about related issues, but say I have this anisotropic dispersion tensor:<br>
<br>
D = (D_m + alpha_L*norm(v))*Identity + (alpha_L - alpha_T)*outer(v,v)/norm(v)<br>
<br>
The anisotropy of this tensor comes from the ratio of alpha_L/alpha_T and velocity v, and is mesh independent. That said, can tuning the multi-grid preconditioners parameters of GAMG/HYPRE/etc help prevent things like negative concentrations and violations of the discrete maximum principles?<br></blockquote><div><br></div><div>In my view (Barry and Jed can tell me I am wrong), this is the job of the interpolator. You already guarantee the bound</div><div>in your smoothers, so this is where is gets messed up. You can certainly use your own interpolators in GMG. GAMG</div><div>is harder since it makes the aggregates, and I do not understand Hypre.</div><div><br></div><div>The way I imagine doing this is to setup the interpolation as a saddle point problem, where the thing you want to enforce</div><div>is a constraint, so the interpolation is a solve instead of jsut a matrix. This is exactly what I did for my Phd when we transferred</div><div>fluid-structure solutions from one mesh to another.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks,</blockquote></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>