Hello Bill,
I have attached my geometry and desired boundary conditions.
Thank you for your time,
Ohi
________________________________
From: overture-bounce@xxxxxxxxxxxxx <overture-bounce@xxxxxxxxxxxxx> on behalf
of Bill Henshaw <wdhenshaw@xxxxxxxxx>
Sent: Tuesday, July 18, 2017 5:23:08 PM
To: overture@xxxxxxxxxxxxx
Subject: [overture] Re: Matrix has a row with all zero coefficients
I'm afraid I don't follow what you are saying. Can you maybe make a sketch of
the geometry and boundary conditions you want to use?
On Tue, Jul 18, 2017 at 2:46 PM, Ohiremen L Dibua
<odibua@xxxxxxxxxxxx<mailto:odibua@xxxxxxxxxxxx>> wrote:
Hello Bill,
I tested it on tcm3 as well and saw that it looked good. I have been going
through my code and saw that it looks fine as well. In those cases though, I
think that the boundary conditions are the same for every grid, where as for me
they are mixed. I think that one issue is that the current geometry is assuming
a boundary layer problem which causes the artificial discontinuity I am
observing whenever the corner boundary intersects with other grids. I am
thinking of testing this by attempting to use a fillet to create a similar
geometry/mesh. Do you have any recommendations of other mappings, besides the
smooth polygon, that I can use to create this grid?
Thank you for your time,
Ohi
________________________________
From: overture-bounce@xxxxxxxxxxxxx<mailto:overture-bounce@xxxxxxxxxxxxx>
<overture-bounce@xxxxxxxxxxxxx<mailto:overture-bounce@xxxxxxxxxxxxx>> on behalf
of Bill Henshaw <wdhenshaw@xxxxxxxxx<mailto:wdhenshaw@xxxxxxxxx>>
Sent: Monday, July 17, 2017 1:36:33 PM
To: overture@xxxxxxxxxxxxx<mailto:overture@xxxxxxxxxxxxx>
Subject: [overture] Re: Matrix has a row with all zero coefficients
The solution using the test code tcm3 on your grid looks fine to me.
On Mon, Jul 17, 2017 at 2:07 PM, Ohiremen L Dibua
<odibua@xxxxxxxxxxxx<mailto:odibua@xxxxxxxxxxxx>> wrote:
Hello Bill,
I apologize. I sent the wrong script previously. The above corresponds to the
grid I created. I do not use "prevent hole cutting," and it appears to properly
interpolate. However, my issue is that in the middle of the a column of zero
appears to be the solution where the grids interpolate, rather than the linear
gradient that occurs at either side of the boundary. Based on my intuition, I
do not think this should be the case. Is there a way to combat this issue?
Bestm
Ohi
________________________________
From: overture-bounce@xxxxxxxxxxxxx<mailto:overture-bounce@xxxxxxxxxxxxx>
<overture-bounce@xxxxxxxxxxxxx<mailto:overture-bounce@xxxxxxxxxxxxx>> on behalf
of Bill Henshaw <wdhenshaw@xxxxxxxxx<mailto:wdhenshaw@xxxxxxxxx>>
Sent: Monday, July 17, 2017 10:58:54 AM
To: overture@xxxxxxxxxxxxx<mailto:overture@xxxxxxxxxxxxx>
Subject: [overture] Re: Matrix has a row with all zero coefficients
hi Ohi,
The script you sent does not seem to correspond to the grid in the figure.
Where the grids join in the middle they should be interpolating from one
another (bc=0).
You should not need to use "prevent hole cutting"
...Bill
On Mon, Jul 17, 2017 at 1:18 PM, Ohiremen L Dibua
<odibua@xxxxxxxxxxxx<mailto:odibua@xxxxxxxxxxxx>> wrote:
Hello Bill,
I wanted to follow up on the grid that you sent me. To make it more like the
domain I want, I removed the inlet and added two more main channels and
corners. I did not have to use hole cutting here, but there appears to be an
issue when I solve the Laplace equation within this domain. At the points where
the left and right sets of grid are interpolated, there appears to be a
discontinuity in the solution. I have attached the Ogen file, a picture of my
mesh grid, and a picture of the solution. I have played a lot with the number
of lines in my grid and their overlap length, but it appears to have no effect.
I wanted to know if you have any tips for me as to how I could fix/debug this
issue?
Thank you for your time,
Ohi
Attachment:
IMG_20170719_081632.jpg
Description: IMG_20170719_081632.jpg
