## Hybrid discrete (H TN) approximations to the equation of radiative transfer

 dc.contributor.advisor James A. Rossmanith dc.contributor.author Shin, Minwoo dc.contributor.department Mathematics dc.date 2019-09-22T01:07:23.000 dc.date.accessioned 2020-06-30T03:17:32Z dc.date.available 2020-06-30T03:17:32Z dc.date.copyright Wed May 01 00:00:00 UTC 2019 dc.date.embargo 2019-10-22 dc.date.issued 2019-01-01 dc.description.abstract

The linear kinetic transport equations are ubiquitous in many application areas, including as a model for neutron transport in nuclear reactors and the propagation of electromagnetic radiation in astrophysics. The main computational challenge in solving the linear transport equations is that solutions live in a high-dimensional phase space that must be sufficiently resolved for accurate simulations. The three standard computational techniques for solving the linear transport equations are the (1) implicit Monte Carlo, (2) discrete ordinate(S\$_N\$), and (3) spherical harmonic(P\$_N\$) methods. Monte Carlo methods are stochastic methods for solving time-dependent nonlinear radiative transfer problems. In a traditional Monte Carlo method when photons are absorbed, they are reemitted in a distribution which is uniform over the entire spatial cell where the temperature is assumed constant, resulting in loss of information. In implicit Monte Carlo(IMC) methods, photons are reemitted from the place where they were actually absorbed, which improves the accuracy. Overall, IMC method improves stability, flexibility, and computational efficiency \cite{fleck}. The S\$_N\$ method solves the transport equation using a quadrature rule to reconstruct the energy density. This method suffers from so-called "ray effect", which are due to the approximation of the double integral over a unit sphere by a finite number of discrete angular directions \cite{chai}. The P\$_N\$ approximation is based on expanding the part of the solution that depends on velocity direction (i.e., two angular variables) into spherical harmonics. A big challenge with the P\$_N\$ approach is that the spherical harmonics expansion does not prevent the formation of negative particle concentrations. The idea behind my research is to develop on an alternative formulation of P\$_N\$ approximations that hybridizes aspects of both P\$_N\$ and S\$_N\$. Although the basic scheme does not guarantee positivity of the solution, the new formulation allows for the introduction of local limiters that can be used to enforce positivity.

dc.format.mimetype application/pdf dc.identifier archive/lib.dr.iastate.edu/etd/17316/ dc.identifier.articleid 8323 dc.identifier.contextkey 15016628 dc.identifier.s3bucket isulib-bepress-aws-west dc.identifier.submissionpath etd/17316 dc.identifier.uri https://dr.lib.iastate.edu/handle/20.500.12876/31499 dc.language.iso en dc.source.bitstream archive/lib.dr.iastate.edu/etd/17316/Shin_iastate_0097E_17805.pdf|||Fri Jan 14 21:20:27 UTC 2022 dc.subject.disciplines Applied Mathematics dc.subject.keywords discontinuous Galerkin method dc.subject.keywords discrete ordinate dc.subject.keywords hyperbolic PDE dc.subject.keywords positivity-preserving limiters dc.subject.keywords radiative transfer dc.subject.keywords spherical harmonics dc.title Hybrid discrete (H TN) approximations to the equation of radiative transfer dc.type article dc.type.genre dissertation dspace.entity.type Publication relation.isOrgUnitOfPublication 82295b2b-0f85-4929-9659-075c93e82c48 thesis.degree.discipline Applied Mathematics thesis.degree.level dissertation thesis.degree.name Doctor of Philosophy
##### Original bundle
Now showing 1 - 1 of 1
No Thumbnail Available
Name:
Shin_iastate_0097E_17805.pdf
Size:
8.54 MB
Format: