An upwind mixed finite volume element-fractional step method and convergence analysis for three-dimensional compressible contamination treatment from nuclear waste (Q6089351)

The authors discuss a numerical simulation problem of three-dimensional compressible contamination treatment from nuclear waste. The model is a nonlinear convection-diffusion system of four partial differential equations and it determines four major physical unknowns: the pressure, the concentrations of brine and radionuclide, and the temperature. The pressure is approximated using a conservative mixed finite volume element method, and the computational accuracy is improved for Darcy velocity. Other unknowns are approximated using a composite scheme of upwind approximation and mixed finite volume element. Numerical dispersion and nonphysical oscillation are eliminated, and the convection-dominated diffusion problems are solved well with high order computational accuracy. The mixed finite volume element is conservative locally, and get the objective functions and their adjoint vector functions simultaneously. Fractional step difference is introduced to solve the concentrations of radionuclide factors, and the computational work is shortened significantly by decomposing a three-dimensional problem into three successive one-dimensional problems. An optimal order estimates in the \(L^2\)-norm is derived and proved. Some numerical examples are presented to show the effectiveness of the numerical method.