DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Interface- and discontinuity-aware numerical schemes for plasma 3-T radiation diffusion in two and three dimensions

Abstract

A set of numerical schemes is developed for two- and three-dimensional time-dependent 3-T radiation diffusion equations in systems involving multi-materials. To resolve sub-cell structure, interface reconstruction is implemented within any cell that has more than one material. Therefore, the system of 3-T radiation diffusion equations is solved on two- and three-dimensional polyhedral meshes. The focus of the development is on the fully coupling between radiation and material, the treatment of nonlinearity in the equations, i.e., in the diffusion terms and source terms, treatment of the discontinuity across cell interfaces in material properties, the formulations for both transient and steady states, the property for large time steps, and second order accuracy in both space and time. The discontinuity of material properties between different materials is correctly treated based on the governing physics principle for general polyhedral meshes and full nonlinearity. The treatment is exact for arbitrarily strong discontinuity. The scheme is fully nonlinear for the full nonlinearity in the 3-T diffusion equations. Three temperatures are fully coupled and are updated simultaneously. The scheme is general in two and three dimensions on general polyhedral meshes. The features of the scheme are demonstrated through numerical examples for transient problems and steady states. Finally,more » the effects of some simplifications of numerical schemes are also shown through numerical examples, such as linearization, simple average of diffusion coefficient, and approximate treatment for the coupling between radiation and material.« less

Authors:
ORCiD logo [1];  [1]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
OSTI Identifier:
1458921
Alternate Identifier(s):
OSTI ID: 1250021
Report Number(s):
LA-UR-14-29261
Journal ID: ISSN 0021-9991; TRN: US1901521
Grant/Contract Number:  
AC52-06NA25396
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 300; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
73 NUCLEAR PHYSICS AND RADIATION PHYSICS; 97 MATHEMATICS AND COMPUTING; diffusion; radiation; 3-T; coupling; interface; implicit

Citation Formats

Dai, William W., and Scannapieco, Anthony J. Interface- and discontinuity-aware numerical schemes for plasma 3-T radiation diffusion in two and three dimensions. United States: N. p., 2015. Web. doi:10.1016/j.jcp.2015.07.041.
Dai, William W., & Scannapieco, Anthony J. Interface- and discontinuity-aware numerical schemes for plasma 3-T radiation diffusion in two and three dimensions. United States. https://doi.org/10.1016/j.jcp.2015.07.041
Dai, William W., and Scannapieco, Anthony J. Thu . "Interface- and discontinuity-aware numerical schemes for plasma 3-T radiation diffusion in two and three dimensions". United States. https://doi.org/10.1016/j.jcp.2015.07.041. https://www.osti.gov/servlets/purl/1458921.
@article{osti_1458921,
title = {Interface- and discontinuity-aware numerical schemes for plasma 3-T radiation diffusion in two and three dimensions},
author = {Dai, William W. and Scannapieco, Anthony J.},
abstractNote = {A set of numerical schemes is developed for two- and three-dimensional time-dependent 3-T radiation diffusion equations in systems involving multi-materials. To resolve sub-cell structure, interface reconstruction is implemented within any cell that has more than one material. Therefore, the system of 3-T radiation diffusion equations is solved on two- and three-dimensional polyhedral meshes. The focus of the development is on the fully coupling between radiation and material, the treatment of nonlinearity in the equations, i.e., in the diffusion terms and source terms, treatment of the discontinuity across cell interfaces in material properties, the formulations for both transient and steady states, the property for large time steps, and second order accuracy in both space and time. The discontinuity of material properties between different materials is correctly treated based on the governing physics principle for general polyhedral meshes and full nonlinearity. The treatment is exact for arbitrarily strong discontinuity. The scheme is fully nonlinear for the full nonlinearity in the 3-T diffusion equations. Three temperatures are fully coupled and are updated simultaneously. The scheme is general in two and three dimensions on general polyhedral meshes. The features of the scheme are demonstrated through numerical examples for transient problems and steady states. Finally, the effects of some simplifications of numerical schemes are also shown through numerical examples, such as linearization, simple average of diffusion coefficient, and approximate treatment for the coupling between radiation and material.},
doi = {10.1016/j.jcp.2015.07.041},
journal = {Journal of Computational Physics},
number = C,
volume = 300,
place = {United States},
year = {Thu Jul 30 00:00:00 EDT 2015},
month = {Thu Jul 30 00:00:00 EDT 2015}
}

Journal Article:

Citation Metrics:
Cited by: 8 works
Citation information provided by
Web of Science

Save / Share: