Abstract
In the numerical solution of elliptic equations, multiscale methods typically involve two steps: the solution of families of local solutions or multiscale basis functions (an embarrassingly parallel task) associated with subdomains of a domain decomposition of the original domain, followed by the solution of a global problem. In the solution of multiphase flow problems approximated by an operator splitting method one has to solve an elliptic equation every time step of a simulation, that would require that all multiscale basis functions be recomputed.In this work, we focus on the development of a novel method that replaces a full update of local solutions by reusing multiscale basis functions that are computed at an earlier time of a simulation. The procedure is based on classical perturbation theory. It can take advantage of both an offline stage (where multiscale basis functions are computed at the initial time of a simulation) as well as of a good initial guess for velocity and pressure. The formulation of the method is carefully explained and several numerical studies are presented and discussed. They provide an indication that the proposed procedure can be of value in speeding-up the solution of multiphase flow problems by multiscale methods.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.