I will make a few comments on this. - The periodic boundaries are indeed connected cell-to-cell as you speculate. - BCs are initialized by the function setupUpscalingConditions() in setupBoundaryConditions.hpp, starting on line 97. - For the explicit solver, bcs are handled in opm/porsol/euler/EulerUpstreamResidual_impl.hpp, lines 113-130. This is in the method that computes the new explicit values, the operator() of the residual class (which is done in this way to enable ttb parallelization). - For the implicit solver, bcs are handled in opm/porsol/euler/EulerUpstreamImplicit_impl.hpp, in the initObj() method.
Atgeirr Den 16. mai 2013 kl. 15:38 skrev Lars Vingli Odsæter: Hi list, We are trying to understand how boundary conditions (BCs) are implemented in the steady-state code (steadystate_test_implicit). Particularly, the BCs for the transport solver is of interest. Can someone please point out where in the code this is done? Does it exist some kind of documentation? Based on the results for some of our models, it might seems like the boundaries perpendicular to the pressure drop are connected point-wise (or cell-to-cell) for periodic BCs, i.e. so that what flows out of one boundary face, flows into the corresponding boundary face (periodic partner) on the opposite boundary. Is this correct? Best regards, Carl Fredrik and Lars _______________________________________________ Opm mailing list [email protected]<mailto:[email protected]> http://www.opm-project.org/mailman/listinfo/opm
_______________________________________________ Opm mailing list [email protected] http://www.opm-project.org/mailman/listinfo/opm
