Publications
Papers, posters, reports, and other publications produced by the people (past and present) of this laboratory. Conference products are listed after other documents for the given year.
2012
 E.T. Meier, A.H. Glasser, V.S. Lukin, and U. Shumlak, "Modeling open boundaries in dissipative MHD simulation," Journal of Computational Physics 231 (7), 2963 (2012).
2011
 U. Shumlak, R. Lilly, N. Reddell, E. Sousa, and B. Srinivasan, "Advanced physics calculations using a multifluid plasma model," Computer Physics Communications 182, 00104655 (2011).
 B. Srinivasan and U. Shumlak, "Analytical and computational study of the ideal full twofluid plasma model and asymptotic approximations for HallMHD," Physics of Plasmas 18 (9), 092113 (2011).
 J. Loverich, A. Hakim, and U. Shumlak, "A discontinuous Galerkin method for ideal twofluid plasma equations," Communications in Computational Physics 9, 240 (2011).
 W. Lowrie, V.S. Lukin, and U. Shumlak, "A priori mesh quality metric error analysis applied to a highorder finite element method," Journal of Computational Physics 230 (14), 5564 (2011).
 B. Srinivasan, A. Hakim, and U. Shumlak, "Numerical methods for twofluid dispersive fast MHD phenomena," Communications in Computational Physics 10, 183 (2011).
 Eric Meier, Alan Glasser, Vyacheslav Lukin, and Uri Shumlak, "Open boundary conditions for dissipative MHD," FiftyThird Annual American Physical Society Meeting of the Division of Plasma Physics, Salt Lake City, Utah, November 2011.
 Robert Lilly and Uri Shumlak, "Multidimensional Plasma Sheath Modeling Using The Three Fluid Plasma Model," FiftyThird Annual American Physical Society Meeting of the Division of Plasma Physics, Salt Lake City, Utah, November 2011. Poster
 Eder Sousa, Uri Shumlak, and Guang Lin, "Lacunabased Artificial Boundary Condition And Uncertainty Quantification of the TwoFluid Plasma Model," FiftyThird Annual American Physical Society Meeting of the Division of Plasma Physics, Salt Lake City, Utah, November 2011. Poster
 W. Lowrie, U. Shumlak, A.H. Glasser, and V.S. Lukin, "MultiBlock Development and Application to a Shear Flow ZPinch Kink Mode Stabilization," FiftyThird Annual American Physical Society Meeting of the Division of Plasma Physics, Salt Lake City, Utah, November 2011.
 Brian A. Nelson, A.H. Glasser, T.R. Jarboe, C.C. Kim, G.J. Marklin, W. Lowrie, E.T. Meier, R.D. Milroy, U. Shumlak, C.R. Sovinec, J.B. O'Bryan, E. Held, J.Y. Ji, and V.S. Lukin, "Simulations of ICC Experiments by the PSICenter," FiftyThird Annual American Physical Society Meeting of the Division of Plasma Physics, Salt Lake City, Utah, November 2011.
 Noah Reddell, Uri Shumlak, and Guoyong Fu, “Framework for advanced plasma simulations on manycore architectures and results,” FiftyThird Annual American Physical Society Meeting of the Division of Plasma Physics, Salt Lake City, Utah, November 2011.
 S. Gilliam, " A 13Moment TwoFluid Plasma Physics Model Based on a Pearson TypeIV Distribution Function"
2010
 E. T. Meier, V. S. Lukin, and U. Shumlak, "Spectral element spatial discretization error in solving highly anisotropic heat conduction equation," Computer Physics Communications 181, 837 (2010).
 N. Reddell, R. Lilly, U. Shumlak, E. Sousa, and B. Srinivasan, "Computational Algorithm for the MultiFluid Plasma Model  WARPX," FiftySecond Annual American Physical Society Meeting of the Division of Plasma Physics, Chicago, Illinois, November 2010.
 R. Lilly and U. Shumlak, "Plasma Sheath Modeling Using The Three Fluid Plasma Model," FiftySecond Annual American Physical Society Meeting of the Division of Plasma Physics, Chicago, Illinois, November 2010.
 E.T. Meier, A.H. Glasser, V.S. Lukin, and U. Shumlak, "Lacunaebased open boundary conditions for dissipative MHD," FiftySecond Annual American Physical Society Meeting of the Division of Plasma Physics, Chicago, Illinois, November 2010.
 W. Lowrie, U. Shumlak, V.S. Lukin, and A.H. Glasser, "MultiBlock Development and Applications of HiFi," FiftySecond Annual American Physical Society Meeting of the Division of Plasma Physics, Chicago, Illinois, November 2010.
 B.A. Nelson, A.H. Glasser, T.R. Jarboe, C.C. Kim, G.J. Marklin, W. Lowrie, E.T. Meier, R.D. Milroy, U. Shumlak, C.R. Sovinec, J.B. O’Bryan, E. Held, J.Y. Ji, and V.S. Lukin, "Recent Results of ICC Experiment Simulations by the PSICenter," FiftySecond Annual American Physical Society Meeting of the Division of Plasma Physics, Chicago, Illinois, November 2010.
 B. Srinivasan, G.L. Delzanno, X. Tang, and U. Shumlak, "Applications of the TwoFluid Plasma Model to Simulate Instabilities with and without Background Magnetic Fields," FiftySecond Annual American Physical Society Meeting of the Division of Plasma Physics, Chicago, Illinois, November 2010.
 U. Shumlak, R. Lilly, E. Meier, N. Reddell, E. Sousa, and B. Srinivasan, "Advanced Physics Calculations using a MultiFluid Plasma Model," Conference on Computational Physics, Trondheim, Norway, June 2010.
2009
 B. Srinivasan and U. Shumlak, "A Study of 3Dimensional Plasma Configurations using the TwoFluid Plasma Model," ThirtySixth IEEE International Conference on Plasma Science, San Diego, California, June 2009. Poster

B. Srinivasan and U. Shumlak, "Implicit TimeIntegration Schemes for Finite Element TwoFluid Code," FiftyFirst Annual American Physical Society Meeting of the Division of Plasma Physics, Atlanta, Georgia, November 2009. Poster

A. Susanto, "Development of Electromagnetic Solvers for Use with the TwoFluid Plasma Algorithm"
2008
 B. Srinivasan, U. Shumlak, and A. Hakim, "Comparison and Applications of TwoFluid Plasma Algorithms," Thirty Ninth American Institute of Aeronautics & Astronautics Plasmadynamics and Lasers Conference, Seattle, Washington, June 2008. AIAA20083787
 B. Srinivasan and U. Shumlak, "Comparison of TwoFluid Plasma Models," High Performance Computing Users Group Conference, Seattle, Washington, July 2008. Paper
 B. Srinivasan and U. Shumlak, "A Comparison between the TwoFluid Plasma Model and HallMHD for Captured Physics and Computational Effort," Fiftieth Annual American Physical Society Meeting of the Division of Plasma Physics, Dallas, Texas, November 2008. Poster
2007
 A. Hakim and U. Shumlak, "TwoFluid Physics and FieldReversed Configurations," Physics of Plasmas 14, 055911 (2007).
Abstract. In this paper, algorithms for the solution of twofluid plasma equations are presented and applied to the study of fieldreversed configurations (FRCs). The twofluid model is more general than the often used magnetohydrodynamic (MHD) model. The model takes into account electron inertia, charge separation, and the full electromagnetic field equations, and it allows for separate electron and ion motion. The algorithm presented is the highresolution wave propagation scheme. The wave propagation method is based on solutions to the Riemann problem at cell interfaces. Operator splitting is used to incorporate the Lorentz and electromagnetic source terms. The algorithms are benchmarked against the Geospace Environmental Modeling Reconnection Challenge problem. Equilibrium of FRC is studied. It is shown that starting from a MHD equilibrium produces a relaxed twofluid equilibrium with strong flows at the FRC edges due to diamagnetic drift. The azimuthal electron flow causes lowerhybrid drift instabilities (LHDI), which can be captured if the ion gyroradius is well resolved. The LHDI is known to be a possible source of anomalous resistivity in many plasma configurations. LHDI simulations are performed in slab geometries and are compared to recent experimental results.
 B. Srinivasan and U. Shumlak, "Applying Asymptotic Approximations to the Full TwoFluid Plasma System to Study Reduced Fluid Models," ThirtyFourth Annual IEEE International Conference on Plasma Sciences, Albuquerque, New Mexico, June 2007. Poster
 B. Srinivasan and U. Shumlak, "Comparing Reduced Fluid Models to the Full TwoFluid Plasma System," FortyNinth Annual American Physical Society Meeting of the Division of Plasma Physics, Orlando, Florida, November 2007. Poster
2006
 J. Loverich and U. Shumlak, "NonLinear TwoFluid Study of m=0 Sausage Instabilities in an Axisymmetric ZPinch," Physics of Plasmas 13, 082310 (2006).
Abstract. A nonlinear full fivemoment twofluid model is used to study axisymmetric instabilities in a Z pinch. When the electron velocity due to the current J is greater than the ion acoustic speed, high wavenumber sausage instabilities develop that initiate shock waves in the ion fluid. This condition corresponds to a pinch radius on the order of a few ion Larmor radii.
 A. Hakim, J. Loverich, and U. Shumlak, "High Resolution Wave Propagation Scheme for Ideal TwoFluid Plasma Equations," Journal of Computational Physics 219 (1), 418442 (2006).
Abstract. Algorithms for the solution of the fivemoment ideal TwoFluid equations are presented. The ideal TwoFluid model is more general than the often used magnetohydrodynamic (MHD) model. The model takes into account electron inertia effects, charge separation and the full electromagnetic field equations and allows for separate electron and ion motion. The algorithm presented is the high resolution wave propagation method. The wave propagation method is based on solutions to the Riemann problem at cell interfaces. Operator splitting is used to incorporate the Lorentz and electromagnetic source terms. To preserve the divergence constraints on the electric and magnetic fields two different approaches are used. In the first approach Maxwell equations are rewritten in their mixedpotential form. In the second approach the socalled perfectly hyperbolic form of Maxwell equations are used which explicitly incorporate the divergence equations into the time stepping scheme. The algorithm is applied to a onedimensional Riemann problem, ionacoustic soliton propagation and magnetic reconnection. In each case TwoFluid physics described by the ideal TwoFluid model is highlighted.
 A. Hakim, "High Resolution Wave Propagation Schemes for TwoFluid Plasma Simulations" PhD Dissertation
 R. C. Lilly, "Study on the FlowThrough ZPinch Fusion Concept" MS Thesis
2005
 J. Loverich and U. Shumlak, "A Discontinuous Galerkin Method for the Full TwoFluid Plasma Model," Computer Physics Communications 169 (3), 251 (2005).
Abstract. A discontinuous Galerkin method for the full twofluid plasma model is described. The plasma model includes complete electron and ion fluids, which allows charge separation, separate electron and ion temperatures and velocities. Complete Maxwell's equations are used including displacement current. The algorithm is validated by benchmarking against existing plasma algorithms on the GEM Challenge magnetic reconnection problem. The algorithm can be easily extended to three dimensions, higher order accuracy, general geometries and parallel platforms.
 J. Loverich and U. Shumlak, "The Discontinuous Galerkin Method for Simulation of Plasma Dynamics: The Full TwoFluid System," FortySeventh Annual American Physical Society Meeting of the Division of Plasma Physics, Denver, Colorado, November 2005.
 A. Hakim and U. Shumlak, "A High Resolution Wave Propagation Scheme for TwoFluid Plasma Equations with applications to Field Reversed Configurations," FortySeventh Annual American Physical Society Meeting of the Division of Plasma Physics, Denver, Colorado, November 2005.
 S. Vadlamani, G. Marklin, U. Shumlak, T. Jarboe, and R. Lionello, "Development of MH4D for Simulation of Emerging Concept Fusion Experiments," FortySeventh Annual American Physical Society Meeting of the Division of Plasma Physics, Denver, Colorado, November 2005.
 T.R. Jarboe, C.C. Kim, G. Marklin, A.I.D. Macnab, R.D. Milroy, B.A. Nelson, U. Shumlak, S. Vadlamani, S. Woodruff, R.A. Bayliss, C.R. Sovinec, E. Held, and JY. Ji, "Overview of The Plasma Science and Innovation Center (PSICenter)," FortySeventh Annual American Physical Society Meeting of the Division of Plasma Physics, Denver, Colorado, November 2005.
 J. Loverich, "A Discontinuous Galerkin Method for the TwoFluid Plasma System and Its Application to the ZPinch" PhD Dissertation
 B. Srinivasan, "A Comparison between the Discontinuous Galerkin Algorithm and the High Resolution Wave Propagation Algorithm for the Full TwoFluid Plasma Model" MS Thesis
2004
 U. Shumlak, C. Aberle, A. Hakim, and J. Loverich, "Plasma Simulation Algorithm for the TwoFluid Plasma Model," European Physical Society /American Physical Society Conference on Computational Physics, Division of Computational Physics, Genoa, Italy, September 2004. Poster
 J. J. Loverich and U. Shumlak, "A Discontinuous Galerkin Method for the Full TwoFluid Plasma Model," European Physical Society /American Physical Society Conference on Computational Physics, Division of Computational Physics, Genoa, Italy, September 2004. Poster
 A. Hakim, U. Shumlak, and J. Loverich, "A High Resolution Scheme for MultiDimensional TwoFluid Plasma Equations," FortySixth Annual American Physical Society Meeting of the Division of Plasma Physics, Savannah, Georgia, November 2004. Poster
 C. Aberle, "Algorithm for Solving Colocated Electromagnetic Fields with Sources" MS Thesis
2003
 U. Shumlak and J. Loverich, "Approximate Riemann Solver for the TwoFluid Plasma Model," Journal of Computational Physics 187 (2), 620 (2003).
Abstract. An algorithm is presented for the simulation of plasma dynamics using the twofluid plasma model. The twofluid plasma model is more general than the magnetohydrodynamic (MHD) model often used for plasma dynamic simulations. The twofluid equations are derived in divergence form and an approximate Riemann solver is developed to compute the fluxes of the electron and ion fluids at the computational cell interfaces and an upwind characteristicbased solver to compute the electromagnetic fields. The source terms that couple the fluids and fields are treated implicitly to relax the stiffness. The algorithm is validated with the coplanar Riemann problem, Langmuir plasma oscillations, and the electromagnetic shock problem that has been simulated with the MHD plasma model. A numerical dispersion relation is also presented that demonstrates agreement with analytical plasma waves.
 J. Loverich, U. Shumlak, C. Aberle, and A. Hakim, "A Finite Volume Algorithm for the TwoFluid Plasma System," Sixteenth American Institute of Aeronautics & Astronautics Computational Fluid Dynamics Conference, Orlando, Florida, June 2003. AIAA20034238
 A. Hakim, U. Shumlak, C. Aberle, and J. Loverich, "Maxwell Equation Solver for Plasma Simulations Based on a Mixed Potential Formulation," Sixteenth American Institute of Aeronautics & Astronautics Computational Fluid Dynamics Conference, Orlando, Florida, June 2003. AIAA20033829
 J. Loverich, "A Finite Volume Algorithm for the Twofluid Plasma System in One Dimension" MS Thesis
2002
 J. J. Loverich and U. Shumlak, "A Finite Volume Scheme for the TwoFluid Plasma System," TwentyNinth Annual IEEE International Conference on Plasma Sciences, Banff, Alberta, Canada, May 2002. Poster
 U. Shumlak, C. S. Aberle, A. Hakim, and J. J. Loverich, "Plasma Simulation Algorithm for the TwoFluid Plasma Model," American Physical Society Meeting of the Division of Computational Physics, San Diego, California, August 2002. Poster
 C. S. Aberle, A. Hakim, and U. Shumlak, "Development of a Maxwell Equation Solver for Application to TwoFluid Plasma Models," American Physical Society Meeting of the Division of Computational Physics, San Diego, California, August 2002. Poster
2001
 U. Shumlak and B. Udrea, "An Approximate Riemann Solver for MHD Computations on Parallel Architectures," Fifteenth American Institute of Aeronautics & Astronautics Computational Fluid Dynamics Conference, Anaheim, California, June 2001. AIAA20012591
 C. Aberle and U. Shumlak, "Application of Analytical Methods to Computing Numerical Flux Jacobians," Fifteenth American Institute of Aeronautics & Astronautics Computational Fluid Dynamics Conference, Anaheim, California, June 2001. AIAA20012610
2000
 W. Vuillemot, U. Shumlak, "MultiTemperature, Thermal & Ion Fraction Effects over Wedge and Bluff Body Shapes in Hypervelocity Flow" TwentySeventh Annual IEEE International Conference on Plasma Sciences, June 2000. Poster
1999
 W.W. Vuillemot and U. Shumlak, "Plasma Effects on Hypersonic Flows," FortyFirst Annual American Physical Society Meeting of the Division of Plasma Physics, Seattle, Washington, November 1999.
 B. Udrea, "An Advanced Implicit Solver for MHD" PhD Dissertation
1998
 J. T. Becerra Sagredo, "Semiimplicit Treatment of the Hall Term in Finite Volume, MHD Computations" MS Thesis