Search

The motion of ions in the complex fluids widely appears in hydrodynamics, geodynamics and geophysics. Setting up suitable mathematical model, performing accurate and efficient numerical simulations are essential to understand the underlying physical principle in these phenomena. In this paper, we present a linear, second-order accurate in time, decoupled and unconditionally energy stable numerical scheme for solving an electrohydrodynamic power-law model with variable electrical conductivity. Based on a logarithm transformation for conductivity and the zero energy contribution property of the nonlinear coupling terms in the model, we derive an equivalent system by introducing a nonlocal auxiliary variable. The two-step backward differentiation scheme and the finite element method are used for the temporal and spatial discretizations, respectively. To decouple the system, the nonlocal splitting technique is employed to yield several linear subsystems, which can be solved very efficiently. Numerical simulations are carried out to demonstrate the second-order convergence in time and the unconditional energy stability. Moreover, the effects of Coulombic force and power-law fluid exponents are numerically investigated. • An electrohydrodynamic model with variable electrical conductivity is found satisfying an energy dissipation law. • The proposed scheme inherits the energy dissipation law at the fully discrete level. • The proposed scheme is a linear, second-order accurate in time and decoupled method. • The proposed scheme can be solved very efficiently with using a nonlocal splitting technique. [ABSTRACT FROM AUTHOR]

The stationary Darcy–Brinkman equations in the double-diffusive convection, which model the heat and mass transfer phenomena, are considered in this paper. Based on a suitable contractive operator, the existence and uniqueness of the problem are firstly proved by using the fixed point theorem. The regularities of the weak solution are also derived. Then, the Newton iterative method is studied for solving the nonlinear discrete system generated from the finite element approximation, including the stability and the optimal error estimates regarding the spatial mesh size and the iterative factor. The analysis indicates that the viscosity coefficient has more impact on the numerical algorithm than the thermal conductivity and the mass diffusivity coefficients. Finally, many numerical examples are shown to confirm the correctness of the theoretical prediction. [ABSTRACT FROM AUTHOR]