Search

We present an approach for constructing finite-volume methods for flux-divergence forms to any order of accuracy defined as the image of a smooth mapping from a rectangular discretization of an abstract coordinate space. Our approach is based on two ideas. The first is that of using higher-order quadrature rules to compute the flux averages over faces that generalize a method developed for Cartesian grids to the case of mapped grids. The second is a method for computing the averages of the metric terms on faces such that freestream preservation is automatically satisfied. We derive detailed formulas for the cases of fourth-order accurate discretizations of linear elliptic and hyperbolic partial differential equations. For the latter case, we combine the method so derived with Runge-Kutta time discretization and demonstrate how to incorporate a high-order accurate limiter with the goal of obtaining a method that is robust in the presence of discontinuities and underresolved gradients. For both elliptic and hyperbolic problems, we demonstrate that the resulting methods are fourth-order accurate for smooth solutions.

We describe an algorithm for the numerical solution of a phase-field model (PFM) of microstructure evolution in polycrystalline materials. The PFM system of equations includes a local order parameter, a quaternion representation of local orientation and a species composition parameter. The algorithm is based on the implicit integration of a semidiscretization of the PFM system using a backward difference formula (BDF) temporal discretization combined with a Newton-Krylov algorithm to solve the nonlinear system at each time step. The BDF algorithm is combined with a coordinate-projection method to maintain quaternion unit length, which is related to an important solution invariant. A key element of the Newton-Krylov algorithm is the selection of a preconditioner to accelerate the convergence of the Generalized Minimum Residual algorithm used to solve the Jacobian linear system in each Newton step. Results are presented for the application of the algorithm to 2D and 3D examples.

A coupled mode system is derived to investigate a three-wave parametric instability leading to energy transfer between co-propagating laser beams crossing in a plasma flow. The model includes beams of finite width refracting in a prescribed transverse plasma flow with spatial and temporal gradients in velocity and density. The resulting paraxial light equations are discretized spatially with a Crank-Nicholson-type scheme, and these algebraic constraints are nonlinearly coupled with ordinary differential equations in time that describe the ion acoustic response. The entire nonlinear differential-algebraic system is solved using an adaptive, backward-differencing method coupled with Newton's method. A numerical study is conducted in two dimensions that compares the intensity gain of the fully time-dependent coupled mode system with the gain computed under the further assumption of a strongly damped ion acoustic response. The results demonstrate a time-dependent gain suppression when the beam diameter is commensurate with the velocity gradient scale length. The gain suppression is shown to depend on time-dependent beam refraction and is interpreted as a time-dependent frequency shift.

We investigate the use of adaptive mesh refinement in the simulation of laser plasma filamentation. A numerical algorithm is constructed to solve model equations consisting of a fluid approximation of a quasineutral plasma combined with a paraxial light propagation model. The algorithm involves high-resolution plasma and light model discretizations on a block-structured, locally refined grid hierarchy, which is dynamically modified during the time integration to follow evolving fine-scale solution features. Comparisons of the efficiency of this approach to that of uniform grid calculations are presented.

This paper describes a numerical method for the solution of plasma fluid equations on block-structured, locally refined grids. The plasmas under consideration are typical of those used for the processing of semiconductors. The governing equations consist of a drift?diffusion model of the electrons, together with an energy equation, coupled via Poisson's equation to a system of Euler equations for each ion species augmented with electric field, collisional, and source/sink terms. A discretization previously developed for a uniform spatial grid is generalized to enable local grid refinement. This extension involves the time integration of the discrete system on a hierarchy of levels, each of which represents a degree of refinement, together with synchronization steps to ensure consistency across levels. This approach represents an advancement of methodologies developed for neutral flows using block-structured adaptive mesh refinement (AMR) to include the significant additional effect of the electrostatic forces that couple the ion and electron fluid components. Numerical results that assess the accuracy and efficiency of the method and illustrate the importance of using adequate resolution are also presented.

This paper describes a numerical method for the solution of a system of plasma fluid equations. The fluid model is similar to those employed in the simulation of high-density, low-pressure plasmas used in semiconductor processing. The governing equations consist of a drift-diffusion model of the electrons, together with an internal energy equation, coupled via Poisson's equation to a system of Euler equations for each ion species augmented with electrostatic force, collisional, and source/sink terms. The time integration of the full system is performed using an operator splitting that conserves space charge and avoids dielectric relaxation timestep restrictions. The integration of the individual ion species and electrons within the time-split advancement is achieved using a second-order Godunov discretization of the hyperbolic terms, modified to account for the significant role of the electric field in the propagation of acoustic waves, combined with a backward Euler discretization of the parabolic terms. Discrete boundary conditions are employed to accommodate the plasma sheath boundary layer on underresolved grids. The algorithm is described for the case of a single Cartesian grid as the first step toward an implementation on a locally refined grid hierarchy in which the method presented here may be applied on each refinement level.

A preconditioner for the conjugate gradient solution of the linear system arising from a standard finite element discretization of a self-adjoint, second-order elliptic problem is proposed. To achieve an efficient implementation on multiprocessors, the preconditioner employs a domain decomposition strategy designed to substantially reduce the amount of interdomain (and therefore interprocessor) communication required to invert the preconditioner against the residual vector in each conjugate gradient iteration. A spectral equivalence result and some numerical experiments are presented.

The goal of this project was to investigate the utility of parallel adaptive mesh refinement (AMR) in the simulation of laser plasma interaction (LPI). The scope of work included the development of new numerical methods and parallel implementation strategies. The primary deliverables were (1) parallel adaptive algorithms to solve a system of equations combining plasma fluid and light propagation models, (2) a research code implementing these algorithms, and (3) an analysis of the performance of parallel AMR on LPI problems. The project accomplished these objectives. New algorithms were developed for the solution of a system of equations describing LPI. These algorithms were implemented in a new research code named ALPS (Adaptive Laser Plasma Simulator) that was used to test the effectiveness of the AMR algorithms on the Laboratory's large-scale computer platforms. The details of the algorithm and the results of the numerical tests were documented in an article published in the Journal of Computational Physics [2]. A principal conclusion of this investigation is that AMR is most effective for LPI systems that are ''hydrodynamically large'', i.e., problems requiring the simulation of a large plasma volume relative to the volume occupied by the laser light. Since the plasma-only regions require less more » resolution than the laser light, AMR enables the use of efficient meshes for such problems. In contrast, AMR is less effective for, say, a single highly filamented beam propagating through a phase plate, since the resulting speckle pattern may be too dense to adequately separate scales with a locally refined mesh. Ultimately, the gain to be expected from the use of AMR is highly problem-dependent. One class of problems investigated in this project involved a pair of laser beams crossing in a plasma flow. Under certain conditions, energy can be transferred from one beam to the other via a resonant interaction with an ion acoustic wave in the crossing region. AMR provides an effective means of achieving adequate resolution in the crossing region while avoiding the expense of using the same fine grid everywhere, including the region between the beams where no LPI occurs. We applied ALPS to a suite of problems modeling crossed beam experiments performed on the Omega laser at the University of Rochester. Our simulations contributed to the theoretical interpretation of these experiments, which was recently published in Physical Review Letters [4]. This project has advanced the Laboratory's computational capabilities in the area of AMR algorithms and their application to LPI problems. The knowledge gained and software developed will contribute to the computational tools available for use in the design and interpretation of experiments to be performed at the National Ignition Facility (NIF) in support of Laboratory missions in stockpile stewardship, energy research and high energy density science. « less

A new method is presented for calculating the scattering of an arbitrary electromagnetic wave by a bounded, perfectly conducting body of general shape. The strategy is to replace the corresponding exterior boundary-value problem by an approximate problem on the boundary of the scattering body. This involves the introduction of a certain bilinear form and non-local boundary operator, together with the use of a special class of known solutions of the reduced Maxwell's equations satisfying the Sommerfeld radiation boundary conditions at infinity. Two computer programs implementing this method are described and numerical results showing the successful application of this method to some model problems are presented.

In this paper we consider the numerical solution of the Helmholtz equation in an exterior domain in $\mathbb{R}^3 $. We introduce a new variational formulation of the problem involving a nonlocal boundary condition. We show that the discrete bilinear form associated with our variational formulation satisfies the “inf-sup” condition, and we give a quasi-optimal error estimate for the approximate solution. Moreover, we show a certain density property of our subspaces which plays a crucial role in the justification of the method.