<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Jan 16, 2014 at 3:04 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On Thu, Jan 16, 2014 at 2:54 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>

> On Thu, Jan 16, 2014 at 2:36 PM, Dharmendar Reddy <<a href="mailto:dharmareddy84@gmail.com">dharmareddy84@gmail.com</a>><br>
> wrote:<br>
>><br>
>> Hello,<br>
>>          I have a request related to cell geometry for FVM. I am<br>
>> working on implementing a solver for coupled electron hole transport<br>
>> problem in a semiconductor.<br>
>><br>
>> For the descretization (see the link [1] below), i want to use the FVM<br>
>> which uses voronoi box around each node as the control volume. Current<br>
>> implementation of the ComputeCellGeometry, provide the centroid of the<br>
>> cell.<br>
>><br>
>> Can i request for functionality similar to ComputeCellGeometry, which<br>
>> on call returns, the orthocenter and the area perpendiclar to each<br>
>> edge for flux computation ?<br>
><br>
><br>
> 1) I had to lookup orthocenter and I got this:<br>
> <a href="http://www.mathopenref.com/triangleorthocenter.html" target="_blank">http://www.mathopenref.com/triangleorthocenter.html</a><br>
>      which only defines it in 2D, and says it can be outside the triangle.<br>
> This does not sound promising.<br>
>      If you explain the calculation to me, I can help you implement it.<br>
><br>
</div>    I have Fortran implementation, i have not fully tested it but has<br>
the core required formulas. I willtry to clean it and send it to you.<br>
<div class="im"><br>
> 2) I do give back the edge area if you pass in the edge. Have you looked at<br>
> TS ex12? It shows me<br>
>     using this routine.<br>
><br>
<br>
</div>I will take a look at this one.<br>
<div class="im"><br>
> I will take a look at the discretization and see if I can understand it (I<br>
> see they say there is a problem with<br>
> obtuse elements).<br>
><br>
<br>
</div>Yes there will be a problem, if the element is obtuse which results in<br>
the orthocenter being outside the element.  The flux computation step<br>
will be incorrect in that scenario. My other post on fixing non<br>
Delaunay elements is related to this issue.</blockquote><div><br></div><div>Then I would at least consider the FVM we have implemented in TS ex11. It</div><div>does not have the sensitive mesh dependence, is 2nd order TVD, and works</div>
<div>in both 2D and 3D. You can see a discussion of the code in my slides from Paris</div><div>on the website.</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">
<div class="HOEnZb"><div class="h5">
<br>
>    Matt<br>
><br>
>><br>
>> Or any suggested references to look for geometry related routines.<br>
>><br>
>> [1] <a href="http://www.iue.tuwien.ac.at/phd/triebl/node30.html" target="_blank">http://www.iue.tuwien.ac.at/phd/triebl/node30.html</a><br>
>><br>
>><br>
>> Thanks<br>
>> Reddy<br>
><br>
><br>
><br>
><br>
> --<br>
> What most experimenters take for granted before they begin their experiments<br>
> is infinitely more interesting than any results to which their experiments<br>
> lead.<br>
> -- Norbert Wiener<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>