<div dir="ltr"><div>Since your using staggered grids, the physical boundary of your domain will be approximate by a "stair-case" type of boundary. (Correct me if this is not what you were thinking to do) Thus, imposing traction boundary conditions on the stair-case boundary should be no more complex than it was in your standard cube domain with staggered grids. The only exception is that you have to implement, in a cell-wise manner, the imposition of the traction condition. The functionality should already exist in your original cube staggered grid implementation, but possible the implementation of this boundary condition was done "wall-wise" rather than cell-wise.<br>
</div><div><div><br><br><br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
> 3. I'm trying to put both domains in a single matrix to avoid the<br>
> difficulty I would have if I want to consider only the domain A. In this<br>
> case I would need a traction free boundary condition on the irregular<br>
> boundary of domain A, and it seems a bit too challenging for me to<br>
> incorporate it with the staggered grid. If there is an idea to implement<br>
> this and if you think this could be more suitable than the approach in 2<br>
> above, I would like to learn about that too!<br>
<br>
</div>Complexity of implementing boundary conditions on staggered grids is one<br>
reason some people turn to other discretization technology, such as<br>
finite elements.<br>
</blockquote></div><br></div></div></div></div>