A Numerical Methodology for Enforcing Maximum Principles and the Non-Negative Constraint for Transient Diffusion Equations
Year: 2016
Communications in Computational Physics, Vol. 19 (2016), Iss. 1 : pp. 53–93
Abstract
Transient diffusion equations arise in many branches of engineering and applied sciences (e.g., heat transfer and mass transfer), and are parabolic partial differential equations. It is well-known that these equations satisfy important mathematical properties like maximum principles and the non-negative constraint, which have implications in mathematical modeling. However, existing numerical formulations for these types of equations do not, in general, satisfy maximum principles and the non-negative constraint. In this paper, we present a methodology for enforcing maximum principles and the non-negative constraint for transient anisotropic diffusion equation. The proposed methodology is based on the method of horizontal lines in which the time is discretized first. This results in solving steady anisotropic diffusion equation with decay equation at every discrete time-level. We also present other plausible temporal discretizations, and illustrate their shortcomings in meeting maximum principles and the non-negative constraint. The proposed methodology can handle general computational grids with no additional restrictions on the time-step. We illustrate the performance and accuracy of the proposed methodology using representative numerical examples. We also perform a numerical convergence analysis of the proposed methodology. For comparison, we also present the results from the standard single-field semi-discrete formulation and the results from a popular software package, which all will violate maximum principles and the non-negative constraint.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/cicp.180615.280815a
Communications in Computational Physics, Vol. 19 (2016), Iss. 1 : pp. 53–93
Published online: 2016-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 41
-
Qualitative properties of discrete nonlinear parabolic operators
Horváth, Róbert | Faragó, István | Karátson, JánosNumerische Mathematik, Vol. 143 (2019), Iss. 3 P.529
https://doi.org/10.1007/s00211-019-01062-6 [Citations: 0] -
Large-Scale Optimization-Based Non-negative Computational Framework for Diffusion Equations: Parallel Implementation and Performance Studies
Chang, J. | Karra, S. | Nakshatrala, K. B.Journal of Scientific Computing, Vol. 70 (2017), Iss. 1 P.243
https://doi.org/10.1007/s10915-016-0250-5 [Citations: 10] -
Variational inequality approach to enforcing the non-negative constraint for advection–diffusion equations
Chang, J. | Nakshatrala, K.B.Computer Methods in Applied Mechanics and Engineering, Vol. 320 (2017), Iss. P.287
https://doi.org/10.1016/j.cma.2017.03.022 [Citations: 22] -
A Positivity-Preserving Finite Volume Scheme with Least Square Interpolation for 3D Anisotropic Diffusion Equation
Xie, Hui | Xu, Xuejun | Zhai, Chuanlei | Yong, HengJournal of Scientific Computing, Vol. 89 (2021), Iss. 3
https://doi.org/10.1007/s10915-021-01629-2 [Citations: 2] -
Physics-informed machine learning models for predicting the progress of reactive-mixing
Mudunuru, M.K. | Karra, S.Computer Methods in Applied Mechanics and Engineering, Vol. 374 (2021), Iss. P.113560
https://doi.org/10.1016/j.cma.2020.113560 [Citations: 5] -
On Numerical Stabilization in Modeling Double-Diffusive Viscous Fingering
Shabouei, M. | Nakshatrala, K. B.Transport in Porous Media, Vol. 132 (2020), Iss. 1 P.39
https://doi.org/10.1007/s11242-020-01379-z [Citations: 4] -
Chemo-thermal model and Gaussian process emulator for combustion synthesis of Ni/Al composites
Shabouei, Mohammad | Subber, Waad | Williams, Cedric W. | Matouš, Karel | Powers, Joseph M.Combustion and Flame, Vol. 207 (2019), Iss. P.153
https://doi.org/10.1016/j.combustflame.2019.05.038 [Citations: 16] -
Do Current Lattice Boltzmann Methods for Diffusion and Advection-Diffusion Equations Respect Maximum Principle and the Non-Negative Constraint?
Karimi, S. | Nakshatrala, K. B.Communications in Computational Physics, Vol. 20 (2016), Iss. 2 P.374
https://doi.org/10.4208/cicp.181015.270416a [Citations: 6] -
Optimal finite-differences discretization for the diffusion equation from the perspective of large-deviation theory
Smith, Naftali R
Journal of Statistical Mechanics: Theory and Experiment, Vol. 2024 (2024), Iss. 4 P.043201
https://doi.org/10.1088/1742-5468/ad363f [Citations: 1] -
A scalable variational inequality approach for flow through porous media models with pressure-dependent viscosity
Mapakshi, N.K. | Chang, J. | Nakshatrala, K.B.Journal of Computational Physics, Vol. 359 (2018), Iss. P.137
https://doi.org/10.1016/j.jcp.2018.01.022 [Citations: 12] -
Mechanics-Based Solution Verification for Porous Media Models
Shabouei, M. | Nakshatrala, K. B.Communications in Computational Physics, Vol. 20 (2016), Iss. 5 P.1127
https://doi.org/10.4208/cicp.OA-2016-0007 [Citations: 10] -
Local-Mean Preserving Post-Processing Step for Non-Negativity Enforcement in PET Imaging: Application to 90Y-PET
Millardet, Maël | Moussaoui, Saïd | Mateus, Diana | Idier, Jérôme | Carlier, ThomasIEEE Transactions on Medical Imaging, Vol. 39 (2020), Iss. 11 P.3725
https://doi.org/10.1109/TMI.2020.3003428 [Citations: 2]