Manuel Quezada de Luna

Manuel Quezada de Luna, David I. Ketcheson

We study two-dimensional wave propagation in materials whose properties vary periodically in one direction only. High order homogenization is carried out to derive a dispersive effective medium approximation. One-dimensional materials with constant impedance exhibit no effective dispersion. We show that a new kind of effective dispersion may arise in two dimensions, even in materials with constant impedance. This dispersion is a macroscopic effect of microscopic diffraction caused by spatial variation in the sound speed. We analyze this dispersive effect by using high-order homogenization to derive an anisotropic, dispersive effective medium. We generalize to two dimensions a homogenization approach that has been used previously for one-dimensional problems. Pseudospectral solutions of the effective medium equations agree to high accuracy with finite volume direct numerical simulations of the variable-coefficient equations.

Manuel Quezada de Luna, J. Haydel Collins, Christopher E. Kees

We present a robust numerical method for solving incompressible, immiscible two-phase flows. The method extends the monolithic phase conservative level set method with embedded redistancing by Quezada de Luna et al. [38] and a semi-implicit high-order projection scheme for variable-density flows by Guermond and Salgado [17]. The level set method can be initialized conveniently via a simple phase indicator field, which is pre-processed to obtain an approximate signed distance function. To do this, we propose a new PDE-based redistancing method. We also improve the scheme in [38] to provide more accuracy and robustness in full two-phase flow simulations. Specifically, we perform an extra step to ensure convergence to the signed distance level set function and simplify other aspects of the original scheme. Lastly, we introduce consistent artificial viscosity to stabilize the momentum equations in the context of the projection scheme. This stabilization is algebraic, has no tunable parameters and is suitable for unstructured meshes and arbitrary refinement levels. The overall methodology includes few numerical tuning parameters; however, for the wide range of problems that we solve, we identify only one parameter that strongly affects performance of the computational model and provide a value that provides accurate results across all the benchmarks presented. The result is a robust, accurate, and efficient two-phase flow model, which is mass- and volume-conserving on unstructured meshes and has low user input requirements for real applications.