Search

In this work, we develop recent research on the fully mixed virtual element method (mixed-VEM) based on the Banach space for the stationary Boussinesq equation to suggest and analyze a new mixed-VEM for the stationary two-dimensional Boussinesq equation with temperature-dependent parameters in terms of the pseudostress, vorticity, velocity, pseudoheat vector and temperature fields. The well-posedness of the continuous formulation is analyzed utilizing a fixed-point strategy, a smallness assumption on the data, and some additional regularities on the solution. The discretization for the mentioned variables is based on the coupling $\mathbb{H}(\mathbf{div}_{6/5})$ -- and $\mathbf{H}(\mathrm{div}_{6/5})$ -- conforming virtual element techniques. The proposed scheme is rewritten as an equivalent fixed point operator equation, so that its existence and stability estimates have been proven. In addition, an a priori convergence analysis is established by utilizing the C\'ea estimate and a suitable assumption on data for all variables in their natural norms showing an optimal rate of convergence. Finally, several numerical examples are presented to illustrate the performance of the proposed method.

The present article introduces, mathematically analyzes, and numerically validates a new weak Galerkin (WG) mixed-FEM based on Banach spaces for the stationary Navier--Stokes equation in pseudostress-velocity formulation. More precisely, a modified pseudostress tensor, called $ \boldsymbol{\sigma} $, depending on the pressure, and the diffusive and convective terms has been introduced in the proposed technique, and a dual-mixed variational formulation has been derived where the aforementioned pseudostress tensor and the velocity, are the main unknowns of the system, whereas the pressure is computed via a post-processing formula. Thus, it is sufficient to provide a WG space for the tensor variable and a space of piecewise polynomial vectors of total degree at most 'k' for the velocity. Moreover, in order to define the weak discrete bilinear form, whose continuous version involves the classical divergence operator, the weak divergence operator as a well-known alternative for the classical divergence operator in a suitable discrete subspace is proposed. The well-posedness of the numerical solution is proven using a fixed-point approach and the discrete versions of the Babu\v{s}ka-Brezzi theory and the Banach-Ne\v{c}as-Babu\v{s}ka theorem. Additionally, an a priori error estimate is derived for the proposed method. Finally, several numerical results illustrating the method's good performance and confirming the theoretical rates of convergence are presented.

In this article, we propose and analyze a fully coupled, nonlinear, and energy-stable virtual element method (VEM) for solving the coupled Poisson-Nernst-Planck (PNP) and Navier--Stokes (NS) equations modeling microfluidic and electrochemical systems (diffuse transport of charged species within incompressible fluids coupled through electrostatic forces). A mixed VEM is employed to discretize the NS equations whereas classical VEM in primal form is used to discretize the PNP equations. The stability, existence and uniqueness of solution of the associated VEM are proved by fixed point theory. Global mass conservation and electric energy decay of the scheme are also proved. Also, we obtain unconditionally optimal error estimates for both the electrostatic potential and ionic concentrations of PNP equations in the $H^{1}$-norm, as well as for the velocity and pressure of NS equations in the $\mathbf{H}^{1}$- and $L^{2}$-norms, respectively. Finally, several numerical experiments are presented to support the theoretical analysis of convergence and to illustrate the satisfactory performance of the method in simulating the onset of electrokinetic instabilities in ionic fluids, and studying how they are influenced by different values of ion concentration and applied voltage. These tests relate to applications in the desalination of water.