Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Inertial Response of Wind Power Plants: A Comparison of Frequency-Based Inertial Control and Stepwise Inertial Control

The frequency regulation capability of a wind power plant plays an important role in enhancing frequency reliability especially in an isolated power system with high wind power penetration levels. A comparison of two types of inertial control methods, namely frequency-based inertial control (FBIC) and stepwise inertial control (SIC), is presented in this paper. Comprehensive case studies are carried out to reveal features of the different inertial control methods, simulated in a modified Western System Coordination Council (WSCC) nine-bus power grid using real-time digital simulator (RTDS) platform. The simulation results provide an insight into the inertial control methods under various scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2001

Segmented Domain Decomposition Multigrid For 3-D Turbomachinery Flows

A Segmented Domain Decomposition Multigrid (SDDMG) procedure was developed for three-dimensional viscous flow problems as they apply to turbomachinery flows. The procedure divides the computational domain into a coarse mesh comprised of uniformly spaced cells. To resolve smaller length scales such as the viscous layer near a surface, segments of the coarse mesh are subdivided into a finer mesh. This is repeated until adequate resolution of the smallest relevant length scale is obtained. Multigrid is used to communicate information between the different grid levels. To test the procedure, simulation results will be presented for a compressor and turbine cascade. These simulations are intended to show the ability of the present method to generate grid independent solutions. Comparisons with data will also be presented. These comparisons will further demonstrate the usefulness of the present work for they allow an estimate of the accuracy of the flow modeling equations independent of error attributed to numerical discretization.

Celestina, M. L.↗
NASA NTRS2004

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗
DOE OSTI2018

The Role of Dredge-up in Double White Dwarf Mergers

We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs (WDs) with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two WDs is $\lesssim$1.0 M ⊙ . Mergers involving a CO and a He WD have been suggested as a possible formation channel for R Coronae Borealis (RCB)–type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of 18 O, the merger must avoid dredging up large amounts of 16 O, or else it will be difficult to produce sufficient 18 O to explain the oxygen ratio observed to be of order unity. Here, we performed a total of nine simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, >10 -2 M ⊙ of 16 O is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO WD with a ~0.1 M ⊙ layer of helium on top, we find that no 16 O is being dredged up, while in the q = 0.8 simulation <10 -4 M ⊙ of 16 O has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Model Predictive Control of A Matrix-Converter Based Solid State Transformer for Utility Grid Interaction

The matrix converter solid state transformer (MC-SST), formed from the back-to-back connection of two three-to-single-phase matrix converters, is studied for use in the interconnection of two ac grids. The matrix converter topology provides a light weight and low volume single-stage bidirectional ac-ac power conversion without the need for a dc link. Thus, the lifetime limitations of dc-bus storage capacitors are avoided. However, space vector modulation of this type of MC-SST requires to compute vectors for each of the two MCs, which must be carefully coordinated to avoid commutation failure. An additional controller is also required to control power exchange between the two ac grids. In this paper, model predictive control (MPC) is proposed for an MC-SST connecting two different ac power grids. The proposed MPC predicts the circuit variables based on the discrete model of MC-SST system and the cost function is formulated so that the optimal switch vector for the next sample period is selected, thereby generating the required grid currents for the SST. Simulation and experimental studies are carried out to demonstrate the effectiveness and simplicity of the proposed MPC for such MC-SST-based grid interfacing systems.

Xue, Yaosuo↗
DOE OSTI2022

Decentralized Distribution System Restoration with Grid-Forming/Following Inverter-Based Resources

The high penetration of distributed energy resources (DERs) in active distribution systems has posed challenges to the centralized distribution system restoration (DSR) strategies in current practice. On the other hand, the advancement in smart inverter technologies enables the bottom-up restoration capability. This paper is motivated to develop a 3-layered hierarchical framework for decentralized DSR, based on the grid-forming (GFM) and grid-following (GFL) grid-edge inverters. The first layer presents the tertiary control, which determines the load pickup schedule and generation dispatch of DERs, using the alternating direction method of the multipliers algorithm. The second layer consists of two control functions: GFM control, which regulates voltage and frequency, establishing a stable grid for GFL inverters to follow; and GFL control, which regulates the real and reactive power. In the third layer, the primary control is proposed to regulate the inverter voltage and current, which is developed based on the virtual oscillator control (VOC). Furthermore, the developed framework is tested in the modified IEEE 13-node test feeder. Two scenarios of grid-connected and islanded operating modes are designed, and simulation results demonstrate the effectiveness of decentralized DSR strategies for controlling grid-edge inverters to enhance the distribution system resilience.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

ASGarD: Adaptive Sparse Grid Discretization

Many areas of science exhibit physical processes that are described by high dimensional partial differential equations (PDEs), e.g., the 4D, 5D and 6D models describing magnetized fusion plasmas, models describing quantum chemistry, or derivatives pricing. Such problems are affected by the so-called “curse of dimensionality” where the number of degrees of freedom (or unknowns) required to be solved for scales as N D where N is the number of grid points in any given dimension D. A simple, albeit naive, 6D example is demonstrated in the left panel of Figure 1. With N = 1000 grid points in each dimension, the memory required just to store the solution vector, not to mention forming the matrix required to advance such a system in time, would exceed an exabyte - and also the available memory on the largest of supercomputers available today. The right panel of Figure 1 demonstrates potential savings for a range of problem dimensionalities and grid resolution. While there are methods to simulate such high-dimensional systems, they are mostly based on Monte-Carlo methods, which rely on a statistical sampling such that the resulting solutions include noise. Since the noise in such methods can only be reduced at a rate proportional to $\sqrt{N_p}$ where N p is the number of Monte-Carlo samples, there is a need for continuum, or grid/mesh-based methods for high-dimensional problems, which both do not suffer from noise and bypass the curse of dimensionality. We present a simulation framework that provides such a method using adaptive sparse grids.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗
DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Integrated State Estimation and Contingency Analysis Software Implementation using High Performance Computing Techniques

Power system simulation tools are traditionally developed in sequential mode and codes are optimized for single core computing only. However, the increasing complexity in the power grid models requires more intensive computation. The traditional simulation tools will soon not be able to meet the grid operation requirements. Therefore, power system simulation tools need to evolve accordingly to provide faster and better results for grid operations. This paper presents an integrated state estimation and contingency analysis software implementation using high performance computing techniques. The software is able to solve large size state estimation problems within one second and achieve a near-linear speedup of 9,800 with 10,000 cores for contingency analysis application. The performance evaluation is presented to show its effectiveness.

Chen, Yousu↗
DOE OSTI2023

Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

In the last decades, the Arctic climate has changed dramatically, with the loss of multiyear sea ice one of the clearest consequences. These changes have occurred on relatively rapid timescales, and both accurate short-term Arctic prediction (e.g., 10 days to three months) and climate projection of future Arctic scenarios present ongoing challenges. Here we describe a representation of the Arctic ocean and sea ice in a ultrahigh resolution simulation in which the horizontal grid mesh reduces from 8 km at the equator to 2 km at the poles (UH8to2) for the years 2017-2020. We find the simulation reproduces observed distributions of seasonal sea-ice thickness and concentration realistically, although concentration is biased low in the spring and summer and low biases in thickness are found in the central and eastern basins in the fall. Volume, fresh water, and heat transports through key passages are realistic, lying within observationally determined ranges. Climatological comparisons reveal that the UH8to2 Atlantic Water is shallower, warmer, and saltier than the World Ocean Atlas 2018 climatology for 2005-2017 in the eastern basin. Our analysis suggests that these biases, combined with a lack of stratification in the upper 100 m of the simulated ocean, contribute to the winter biases in modeled sea ice thickness. This relationship between biases in the sea ice and ocean points to a potential positive feedback within the model, illuminating challenges for long term model predictive power in a changing Arctic climate.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Medium-Fidelity CFD Modeling of Multicopter Wakes for Airborne Sensor Measurements

A steady-state multicopter simulation capability leveraging automated Cartesian grid generation and a blade element source term rotor model was used to investigate multi-rotor aerodynamics in and out of ground effect. Simulations demonstrated that this model is able to predict reasonably accurate thrust in ground effect for both single and multiple propeller cases. The method yielded accurate lift and drag predictions for a commercial quadcopter when compared to recent wind tunnel test data. This aircraft was then simulated in a variety of flight conditions, including both hover and edgewise forward flight, to determine if generalized guidelines for airborne sensor placement could be developed. Velocity perturbation magnitudes showed large regions of affected air upstream of the vehicle at low speeds, which contracted as the vehicle’s speed was increased. Placing the sensors more than one characteristic length ahead of or above the aircraft reduces errors from rotor-induced flow at higher speeds. For low-speed forward flight, the entrainment of flow into the propeller wakes introduces significant flow turning above the aircraft, suggesting that a forward location, with more moderate perturbations, could be advantageous for these conditions.

ARMD↗
DOE OSTI2014

A Grid Intensity Map-Based Source Model for Monte Carlo Modeling of IMRT Field Irradiation - Paper 5

An accurate source model is important for TPS (Treatment Planning System) dose calculation and verification, especially for IMRT (Intensity Modulated Radiotherapy) irregular and inhomogeneous fields. Monte Carlo simulation is the most accurate way to calculate dose distributions in the heterogeneous media like human bodies. A full radiotherapy modeling generally includes Linac modeling and phantom dose calculation. Full Linac modeling is very time-consuming and requires highly detailed knowledge of the accelerator composition and structure, however, which is very difficult to obtain for commercial interests. Often all relevant parameters are not known with sufficient accuracy, and thus a time-consuming tuning of e.g. the electron beam spot size and energy must be performed. A multiple source model (MSM) regards particles coming from different parts of an accelerator as if they are coming from different virtual sources. MSM uses a set of virtual sources representing the treatment beam (e.g. primary photons, head scatter, electron contamination) by a set of analytic functions or histograms with parameters derived from dosimetric measurements or a full MC simulation. Tillikainen et al built a three sub-source model by matching with measurement data using an automatic process to optimize the free model parameters for a Varian 2100 C/D accelerator. Lin et al have built a convenient five-photon-sub-source model for Varian 2300C 6 MV-X beam based on Monte Carlo simulated Linac for ARTS (Accurate Radiotherapy System). The particle generation of MSM is very fast compared to full Linac simulation, and its sampling efficiency is not limited by the size of the phase space file (PSF). However, generally MSM is generated based on regular field PSF and needs an irregular component to form the beam shape meanwhile apt to bringing error Linac modeling and MSM producing from PSF analysis actually both are multi-freedom optimization process by adjusting several parameters (e.g. electron energy and angular divergence). Linac modeling adjusts its parameters from the beginning of electron hitting target through all components to fit the simulation to the measurement. MSM tunes its PSF-derived parameters by optimizing to minimize the deviations between the simulation and the measurement. Since the full Linac modeling and the MSM model both need to fine tune their parameters to match the simulation to the measurement, and sometimes maybe offsets the real accelerator for its multi-parameter optimization process. So, why not fine tune directly all source model parameters just behind all Linac components to fit the simulation to the measurement? A grid intensity-based surface source model was suggested here, which assumed the particle emitting from the bottom surface of Linac last beam trimmer MLC (Multi Leaf Collimator). The particle emitting position is sampled on the grid intensity map optimized by TPS, and the particle weight is adjusted by its original grid intensity. This source model not only takes advantage of the Monte Carlo method in the heterogeneous media, but also avoids effectively the transportation speed shortcoming of Monte Carlo simulation through Linac components. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2019

Towards Improved Real-Time Observability of Behind-Meter PhotoVoltaic Systems: A Data-Driven Approach

The stochasticity in energy supply introduced by increased PV penetration in grid distribution systems can lead to frequent and sudden supply-demand imbalances, thereby causing voltage and frequency fluctuations. Enhanced observability, in real-time, of grid attributes such as real and reactive power flows, power injections/consumption and voltages is becoming increasingly essential to ensure smooth grid operation under these new large-scale PV integration scenarios. Two significant challenges prevent grid operators from achieving this objective: 1) Computational complexity of quasi-static timeseries simulations, and 2) Partial observability of behind-meter grid assets. In this work, we take a step towards this objective by enhancing the observability of smart inverter equipped PV systems located behind AMI meters. Smart inverters convert the DC output of PV systems into real and reactive power, the ratio of which is governed by a power factor setting. This in turn is governed by the real-time voltage at its terminal, as per a pre-set volt/var curve. While it is possible to infer the power factor setting - thereby, real and reactive power injection of the PV system - at each time interval by quasi-static time series simulations, such methods are not scalable to large grid sizes. We propose to reduce the reliance on simulations by developing a deep neural network based data-driven model which predicts the time dependent power factors of the PV systems in the grid. The model uses a small amount of simulated data to learn the dependence of power factor on 1) load and apparent PV power injection (root mean square of real and reactive power) data, and 2) AMI net-load data (for BTM PV assets). Our results show that our model can predict power factors with mean absolute error of less than 6.2% under full observability scenario and under 8.6% under the partial observability scenario. Thus, our data-driven model can reduce the amount of simulations required to compute real and reactive injections from PV systems in the distribution system with or without complete information.

14 SOLAR ENERGY↗
DOE OSTI2023

Simulating Sea-Ice Deformation in Viscous-Plastic Sea-Ice Models With CD-Grids

Linear kinematic features (LKFs) are found everywhere in the Arctic sea-ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. In viscous-plastic (VP) sea-ice models, the simulation of LKFs depends on several factors such as the grid resolution, the numerical solver convergence, and the placement of the variables on the mesh. In this study, we compare two recently proposed discretization with a CD-grid placement with respect to their ability to reproduce LKFs. The first (CD1) is based on a nonconforming finite element discretization, whereas the second (CD2) uses a conforming subgrid discretization. To analyze their resolution properties, we evaluate runs from different models (e.g., FESOM, MPAS) on a benchmark problem using quadrilateral, hexagonal and triangular meshes. Our findings show that the CD1 setup simulates more deformation structure than the CD2 setup. This highlights the importance of the type of spatial discretization for the simulation of LKFs. Due to the higher number of degrees of freedom, both CD-grids resolve more LKFs than traditional A, B, and C-grids at fixed mesh level. This is an advantage of the CD-grid approach, as high spatial mesh resolution is needed in VP sea-ice models to simulate LKFs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Efficient Simulation of Cascading Outages Using an Energy Function-Embedded Quasi-Steady-State Model

Here, this paper proposed an energy function-embedded quasi-steady-state model for efficient simulation of cascading outages on a power grid while addressing transient stability concerns. Compared to quasi-steady-state models, the proposed model incorporates short-term dynamic simulation and an energy function method to efficiently evaluate the transient stability of a power grid together with outage propagation without transient stability simulation. Cascading outage simulation using the proposed model conducts three steps for each disturbance such as a line outage. First, it performs time-domain simulation for a short term to obtain a post-disturbance trajectory. Second, along the trajectory, the system state with the local maximum potential energy is found and used as the initial point to search for a relevant unstable equilibrium by Newton's method. Third, the transient energy margin is estimated based on this unstable equilibrium to predict an out-of-step condition with generators. The proposed energy function-embedded quasi-steady-state model is tested in terms of its accuracy and time performance on an NPCC 140-bus power system and compared to a quasi-steady-state model embedding transient stability simulation.

Guo, Zhenping [Univ. of Tennessee, Knoxville, TN (↗
DOE OSTI2020

Real-Time Federated Cyber-Transmission-Distribution Testbed Architecture for the Resiliency Analysis

With the ongoing automation driven by the push toward the smart electric grid and the advancement in associated cyber infrastructure, the interaction between physical [electric transmission and distribution (T&D) systems], cyber (communication, automation, and control), and human (grid operators and decision-makers) is increasingly becoming more complex. This creates the requirement of analyzing the effect of the transmission system on the distribution system and vice versa with consideration of the additional complexity of the cyber infrastructure. Such an integrated testbed will also help with resiliency analysis, where resiliency refers to the ability of the system to continue serving energy to the critical loads even with limited extreme contingencies. Interaction of the physical power grid with the cyber layer can be effectively modeled using real-time (RT) simulator for developing and validating various operational and control algorithms. Testbeds using RT simulators with multiple capabilities have been developed at different institutions. Still, no single existing testbed can offer full scalability while simultaneously meeting high fidelity requirements for resiliency experimentation. Co-simulating federated testbed assets can provide a scalable experimentation platform that can be leveraged for verification and validation. In this article, an architecture is developed for federated cyber-physical testbed. A local federation with two real-time simulators is developed: real-time digital simulator (RTDS) and OPAL-RT have been interfaced using VILLAS framework for end-to-end testing. Also, a real-time linear predictor is developed and integrated here to address the communication latency impact on geographically allocated federated RT simulation. Finally, resiliency analysis tools are formulated and utilized for T&D systems. Finally as an illustrative use case, the resiliency of a T&D test system is simulated, and the results are analyzed. A 179-bus Western Electricity Coordinating Council (WECC) transmission system is developed using OPAL-RT/ HYPERSIM, and a modified IEEE 13 node feeder system is modeled in RTDS/RSCAD and interfaced for resiliency analysis.

42 ENGINEERING↗
DOE OSTI2020

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗