<div dir="ltr"><div dir="ltr">On Mon, Jan 17, 2022 at 10:49 AM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Thibault Bridel-Bertomeu <<a href="mailto:thibault.bridelbertomeu@gmail.com" target="_blank">thibault.bridelbertomeu@gmail.com</a>> writes:<br>
<br>
> Sorry I wasn't clear I think.<br>
> I was thinking that for finite volume codes, you have two options: either<br>
> you use ghost cells that you fill up with the appropriate primitive state<br>
> for the boundary condition, or you directly enforce the boundary condition<br>
> at the face.<br>
> In the PetscFV examples, the DM_BC_NATURAL_RIEMANN type of boundary<br>
> condition with ghost cells is always used, but what about _not_ generating<br>
> the ghost cells and using a DM_BC_ESSENTIAL to enforce the boundary<br>
> condition directly at the face (or at the quadrature points on the face if<br>
> you will) ?<br>
<br>
Essential boundary conditions change the function space in which you seek a solution. There is no dof at the face in an FV method, so how would you propose modifying the function space.<br>
<br>
Also, FV methods with higher than first order require a reconstruction (generally nonlinear) using neighbor cells to create a richer (usually piecewise linear) representation inside cells.<br>
</blockquote></div><br clear="all"><div>What you suggest (flux boundary conditions) would require short-circuiting the Riemann solver loop to stick in the flux. It is doable, but we opted against it in the initial</div><div>design because it seemed more complex than ghost cells and did not seem to have any advantages. For multi-domain things, wouldn't you require that the cell average of</div><div>the adjacent cell be used as the ghost value for consistency? We would need a new interface to mark existing cells as ghosts, but that is quite simple compared to a flux BC</div><div>addition.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>