Show simple item record

dc.contributor.advisorGanapol, Barry D.en_US
dc.contributor.advisorMuscat, Anthony J.en_US
dc.contributor.authorFord, Kristopher
dc.creatorFord, Kristopheren_US
dc.date.accessioned2014-06-13T17:58:37Z
dc.date.available2014-06-13T17:58:37Z
dc.date.issued2014
dc.identifier.urihttp://hdl.handle.net/10150/321611
dc.description.abstractThe one dimensional transient heat conduction equation was numerically modeled through matrix diagonalization and three time-discretization schemes. The discrete methods were first-order backward, second-order backward, and implicit finite difference schemes. All simulations used the central difference formula in the space dimension. Two relevant physical systems were considered: a uniformly conducting slab and a melting block of ice. The latter lead to a moving boundary system, or Stefan problem. The multiphysics of melting was numerically modeled through alternating updates of temperature and melt front profiles. Iterative simulations were run with regularly refined discretization meshes in both systems. In the case of the conducting slab, temperature at a fixed point in space and time was considered. For the Stefan problem, the melt front movement after a set time was the physical solution of interest. The accuracy of the convergent results was increased using Richardson acceleration and the Wynn's epsilon algorithm. Accuracy was improved for the moving boundary problem as well, but to a significantly lesser degree. The relative errors improved by five and two orders of magnitude for the conducting block and melting ice simulations, respectively. These relative errors were used to determine that matrix diagonalization is the most accurate numerical solution among the four considered. In both simulation convergence and acceleration potential, matrix diagonalization was superior to the implicit and explicit discretization solutions. However, matrix diagonalization required significantly more computational time. With the enhancement of convergence acceleration, the finite difference schemes obtained similar relative errors to the diagonalization model. This demonstrated the value of convergence acceleration in the classic dilemma for every programmer. There is always a balance struck between model sophistication, accuracy, and computational time. Convergence acceleration allows for a simpler numerical model to achieve comparable accuracy, and in less time than that required for sophisticated numerical models. The numerical models were also compared for stability through parameters that arose in each simulation. These parameters were the Courant-Friederichs-Lewy (CFL) condition and diagonalized eigenvalues. Though diagonalization was found to be the most accurate, it was determined that the backwards finite difference solutions are the easiest to evaluate for stability. In these solution methods, the CFL value allows the stability to be determined prior to running the simulation.
dc.language.isoen_USen
dc.publisherThe University of Arizona.en_US
dc.rightsCopyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author.en_US
dc.subjectConvergence Accelerationen_US
dc.subjectHeat Transferen_US
dc.subjectRichardson Accelerationen_US
dc.subjectStefan Problemen_US
dc.subjectWynn's Epsilon Algorithmen_US
dc.subjectChemical Engineeringen_US
dc.subjectComputationen_US
dc.titleComparison of Convergence Acceleration Algorithms Across Several Numerical Models of 1-Dimensional Heat Conductionen_US
dc.typetexten
dc.typeElectronic Thesisen
thesis.degree.grantorUniversity of Arizonaen_US
thesis.degree.levelmastersen_US
dc.contributor.committeememberGanapol, Barry D.en_US
dc.contributor.committeememberMuscat, Anthony J.en_US
dc.contributor.committeememberSaez, Avelino E.en_US
thesis.degree.disciplineGraduate Collegeen_US
thesis.degree.disciplineChemical Engineeringen_US
thesis.degree.nameM.S.en_US
refterms.dateFOA2018-06-13T04:00:45Z
html.description.abstractThe one dimensional transient heat conduction equation was numerically modeled through matrix diagonalization and three time-discretization schemes. The discrete methods were first-order backward, second-order backward, and implicit finite difference schemes. All simulations used the central difference formula in the space dimension. Two relevant physical systems were considered: a uniformly conducting slab and a melting block of ice. The latter lead to a moving boundary system, or Stefan problem. The multiphysics of melting was numerically modeled through alternating updates of temperature and melt front profiles. Iterative simulations were run with regularly refined discretization meshes in both systems. In the case of the conducting slab, temperature at a fixed point in space and time was considered. For the Stefan problem, the melt front movement after a set time was the physical solution of interest. The accuracy of the convergent results was increased using Richardson acceleration and the Wynn's epsilon algorithm. Accuracy was improved for the moving boundary problem as well, but to a significantly lesser degree. The relative errors improved by five and two orders of magnitude for the conducting block and melting ice simulations, respectively. These relative errors were used to determine that matrix diagonalization is the most accurate numerical solution among the four considered. In both simulation convergence and acceleration potential, matrix diagonalization was superior to the implicit and explicit discretization solutions. However, matrix diagonalization required significantly more computational time. With the enhancement of convergence acceleration, the finite difference schemes obtained similar relative errors to the diagonalization model. This demonstrated the value of convergence acceleration in the classic dilemma for every programmer. There is always a balance struck between model sophistication, accuracy, and computational time. Convergence acceleration allows for a simpler numerical model to achieve comparable accuracy, and in less time than that required for sophisticated numerical models. The numerical models were also compared for stability through parameters that arose in each simulation. These parameters were the Courant-Friederichs-Lewy (CFL) condition and diagonalized eigenvalues. Though diagonalization was found to be the most accurate, it was determined that the backwards finite difference solutions are the easiest to evaluate for stability. In these solution methods, the CFL value allows the stability to be determined prior to running the simulation.


Files in this item

Thumbnail
Name:
azu_etd_13343_sip1_m.pdf
Size:
1.036Mb
Format:
PDF

This item appears in the following Collection(s)

Show simple item record