Numerical resolution of the hyperbolic heat equation using smoothed mathematical functions instead of Heaviside and Dirac delta distributions
Autor: | M. J. Rivera, J. A. López Molina, Enrique Berjano, Vicent Romero-García, Macarena Trujillo |
---|---|
Rok vydání: | 2013 |
Předmět: |
Heaviside distribution
Heaviside step function General Chemical Engineering Multiphysics Numerical analysis Dirac delta function Numerical method Condensed Matter Physics Atomic and Molecular Physics and Optics Polynomial interpolation TECNOLOGIA ELECTRONICA symbols.namesake Distribution (mathematics) Dirac distribution Gaussian function symbols Hyperbolic bioheat equation Applied mathematics Non-Fourier heat MATEMATICA APLICADA Wave heat Mathematics Second derivative |
Zdroj: | RiuNet. Repositorio Institucional de la Universitat Politécnica de Valéncia instname |
ISSN: | 0735-1933 |
DOI: | 10.1016/j.icheatmasstransfer.2013.05.017 |
Popis: | The hyperbolic bioheat equation (HBE) has been used to model heating applications involving very short power pulses. This equation includes two mathematical distributions (Heaviside and Delta) which have to be necessarily substituted for smoothed mathematical functions when the HBE is solved by numerical methods. This study focuses on which type of smoothed functions would be suitable for this purpose, i.e. those which would provide solutions similar to those obtained analytically from the original Heaviside and Delta distributions. The logistic function was considered as a substitute for the Heaviside function, while its derivative and the probabilistic Gaussian function were considered as substitutes for the Delta distribution. We also considered polynomial interpolation functions, in particular, the families of smoothed functions with continuous second derivative without overshoot used by COMSOL Multiphysics. All the smoothed functions were used to solve the HBE by the Finite Element Method (COMSOL Multiphysics), and the solutions were compared to those obtained analytically from the original Heaviside and Delta distributions. The results showed that only the COMSOL smoothed functions provide a numerical solution almost identical to the analytical one. Finally, we demonstrated mathematically that in order to find a suitable smoothed function (f) that must adequately substitute any mathematical distribution (D) in the HBE, the difference D - f must have compact support. (c) 2013 Elsevier Ltd. All rights reserved. This work received financial support from the Spanish "Plan Nacional de I + D + I del Ministerio de Ciencia e Innovacion" Grant No. TEC2011-27133-C02-01 and from Universitat Politenica de Valencia (PAID-06-11 Ref. 1988). V. Romero Garcia is grateful for the support of "Programa de Contratos Post-Doctorales con Movilidad UPV del Campus de Excelencia (CEI-01-11)" and FEDER Project MAT2009-09438. |
Databáze: | OpenAIRE |
Externí odkaz: |