An AMR capable finite element diffusion solver for ALE hydrocodes [An AMR capable diffusion solver for ALEAMR]
Abstract
Here, we present a novel method for the solution of the diffusion equation on a composite AMR mesh. This approach is suitable for including diffusion based physics modules to hydrocodes that support ALE and AMR capabilities. To illustrate, we proffer our implementations of diffusion based radiation transport and heat conduction in a hydrocode called ALEAMR. Numerical experiments conducted with the diffusion solver and associated physics packages yield 2nd order convergence in the L _{2} norm.
 Authors:
 Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
 Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
 Publication Date:
 Research Org.:
 Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 1226207
 Report Number(s):
 LLNLJRNL558861
Journal ID: ISSN 10090630
 DOE Contract Number:
 AC5207NA27344
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Plasma Science and Technology; Journal Volume: 17; Journal Issue: 2
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 42 ENGINEERING; 97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; hydrodynamic simulation; heat conduction; thermal radiation; adaptive mesh refinement; finite element method
Citation Formats
Fisher, A. C., Bailey, D. S., Kaiser, T. B., Eder, D. C., Gunney, B. T. N., Masters, N. D., Koniges, A. E., and Anderson, R. W. An AMR capable finite element diffusion solver for ALE hydrocodes [An AMR capable diffusion solver for ALEAMR]. United States: N. p., 2015.
Web. doi:10.1088/10090630/17/2/04.
Fisher, A. C., Bailey, D. S., Kaiser, T. B., Eder, D. C., Gunney, B. T. N., Masters, N. D., Koniges, A. E., & Anderson, R. W. An AMR capable finite element diffusion solver for ALE hydrocodes [An AMR capable diffusion solver for ALEAMR]. United States. doi:10.1088/10090630/17/2/04.
Fisher, A. C., Bailey, D. S., Kaiser, T. B., Eder, D. C., Gunney, B. T. N., Masters, N. D., Koniges, A. E., and Anderson, R. W. 2015.
"An AMR capable finite element diffusion solver for ALE hydrocodes [An AMR capable diffusion solver for ALEAMR]". United States.
doi:10.1088/10090630/17/2/04. https://www.osti.gov/servlets/purl/1226207.
@article{osti_1226207,
title = {An AMR capable finite element diffusion solver for ALE hydrocodes [An AMR capable diffusion solver for ALEAMR]},
author = {Fisher, A. C. and Bailey, D. S. and Kaiser, T. B. and Eder, D. C. and Gunney, B. T. N. and Masters, N. D. and Koniges, A. E. and Anderson, R. W.},
abstractNote = {Here, we present a novel method for the solution of the diffusion equation on a composite AMR mesh. This approach is suitable for including diffusion based physics modules to hydrocodes that support ALE and AMR capabilities. To illustrate, we proffer our implementations of diffusion based radiation transport and heat conduction in a hydrocode called ALEAMR. Numerical experiments conducted with the diffusion solver and associated physics packages yield 2nd order convergence in the L2 norm.},
doi = {10.1088/10090630/17/2/04},
journal = {Plasma Science and Technology},
number = 2,
volume = 17,
place = {United States},
year = 2015,
month = 2
}
Other availability
Save to My Library
You must Sign In or Create an Account in order to save documents to your library.

A finite element Poisson solver for gyrokinetic particle simulations in a global field aligned mesh
A new finite element Poisson solver is developed and applied to a global gyrokinetic toroidal code (GTC) which employs the field aligned mesh and thus a logically nonrectangular grid in a general geometry. Employing test cases where the analytical solutions are known, the finite element solver has been verified. The CPU time scaling versus the matrix size employing portable, extensible toolkit for scientific computation (PETSc) to solve the sparse matrix is promising. Taking the ion temperature gradient modes (ITG) as an example, the solution from the new finite element solver has been compared to the solution from the original GTC'smore » 
A 2D finite element wave equation solver based on triangular base elements
A finite element method based on the subdivision of the physical domain in triangular subdomains in which simple local 'areale' coordinates are adopted is explored. The advantage of the method is that it straightforwardly allows grid refinement in regions where higher precision is required. The plasma model was kept simple for this 'proofofprinciple' exercise. Rather than accounting for the actual differential or integrodifferential dielectric tensor, its locally uniform plasma equivalent was adopted for 3 possible choices: the cold plasma response, the full hot Stix/Swanson plasma tensor retaining all orders in finite Larmor radius (FLR) and the more common hot tensor,more » 
Plasma wave simulation based on a versatile finite element method solver
The possibility of plasma wave simulation based on the finite element method (FEM) was explored, and a new full wave simulation code of the lower hybrid (LH) wave, Lower Hybrid wavE Analysis based on FEM (LHEAF), was developed. The efficient iterative solver for the LH wave [O. Meneghini, S. Shiraiwa, and R. Parker, Phys. Plasmas 16, 090701 (2009)] was coupled to onedimensional FokkerPlanck calculation to obtain the selfconsistent distribution function. In LHEAF, seamless handling of the core, scrape off layer, and antenna regions were realized and the wave launching structure were naturally included into the simulation. Simulation results of amore »