Recursos
Proyectos/Publicaciones
This article is concerned with the mathematical and numerical analysis of a steady phase change problem for non-isothermal incompressible viscous flow. The system is formulated in terms of pseudostress, strain rate and velocity for the Navier–Stokes–Brinkman equation, whereas temperature, normal heat flux on the boundary, and an auxiliary unknown are introduced for the energy conservation equation. In addition, and as one of the novelties of our approach, the symmetry of the pseudostress is imposed in an ultra-weak sense, thanks to which the usual introduction of the vorticity as an additional unknown is no longer needed. Then, for the mathematical analysis two variational formulations are proposed, namely mixed-primal and fully-mixed approaches, and the solvability of the resulting coupled formulations is established by combining fixed-point arguments, Sobolev embedding theorems and certain regularity assumptions. We then construct corresponding Galerkin discretizations based on adequate finite element spaces, and derive optimal a priori error estimates. Finally, numerical experiments in 2D and 3D illustrate the interest of this scheme and validate the theory.
M. Álvarez, G. Gatica, B. Gómez-Vargas, and R. Ruiz-Baier. New mixed finite element methods for natural convection with phase-change in porous media. Journal of Scientific Computing, vol. 80, pp. 141-174, (2019). DOI: https://link.springer.
In this paper we study a phase change problem for non-isothermal incompressible viscous flows. The underlying continuum is modelled as a viscous Newtonian fluid where the change of phase is either encoded in the viscosity itself, or in the Brinkman–Boussinesq approximation where the solidification process influences the drag directly. We address these and other modelling assumptions and their consequences in the simulation of differentially heated cavity flows of diverse type. A second order finite element method for the primal formulation of the problem in terms of velocity, temperature, and pressure is constructed, and we provide conditions for its stability. We finally present several numerical tests in 2D and 3D, corroborating the accuracy of the numerical scheme as well as illustrating key properties of the model.
J. Woodfield, M. Álvarez, B. Gómez-Vargas, and R. Ruiz-Baier. Stability and finite element approximation of phase change models for natural convection in porous media. Journal of Computational and Applied Mathematics, vol. 360, pp.117-137, (2019).DOI: http://dx.doi.org/10.1016/j.
In this paper we develop the a posteriori error analysis of an augmented mixed-primal finite element method for the 2D and 3D versions of a stationary flow and transport coupled system, typically encountered in sedimentation–consolidation processes. The governing equations consist in the Brinkman problem with concentration-dependent viscosity, written in terms of Cauchy pseudo-stresses and bulk velocity of the mixture; coupled with a nonlinear advection – nonlinear diffusion equation describing the transport of the solids volume fraction. We derive two efficient and reliable residual-based a posteriori error estimators for a finite element scheme using Raviart–Thomas spaces of order k for the stress approximation, and continuous piecewise polynomials of degree ≤k+1 for both velocity and concentration. For the first estimator we make use of suitable ellipticity and inf–sup conditions together with a Helmholtz decomposition and the local approximation properties of the Clément interpolant and Raviart–Thomas operator to show its reliability, whereas the efficiency follows from inverse inequalities and localisation arguments based on triangle-bubble and edge-bubble functions. Next, we analyse an alternative error estimator, whose reliability can be proved without resorting to Helmholtz decompositions. Finally, we provide some numerical results confirming the reliability and efficiency of the estimators and illustrating the good performance of the associated adaptive algorithm for the augmented mixed-primal finite element method.
M. Álvarez, G.N. Gatica and R. Ruiz-Baier. A posteriori error estimation for an augmented mixed-primal method applied to sedimentation-consolidation systems. Journal of Computational Physics, vol. 367, pp. 322- 346, (2018). DOI: http://dx.doi.org/10.
We analyse the stability of a second-order finite element scheme for the primal formulation of a Brinkman-Boussinesq model where the solidification process influences the drag and the viscosity. The problem is written in terms of velocity, temperature, and pressure, and we produce numerical approximations to the flow observed in heated cavities and near ice sheets.
M. Álvarez, B. Gómez-Vargas, R. Ruiz-Baier and J. Woodfield. Stability of a second-order method for phase change in porous media flow. Proceedings in Applied Mathematics and Mechanics, 25 April, (2018). DOI: http://dx.doi.org/10.