Sensitivity technologies for large scale simulation.

Sensitivity analysis is critically important to numerous analysis algorithms, including large scale optimization, uncertainty quantification,reduced order modeling, and error estimation. Our research focused on developing tools, algorithms and standard interfaces to facilitate the implementation of sensitivity type analysis into existing code and equally important, the work was focused on ways to increase the visibility of sensitivity analysis. We attempt to accomplish the first objective through the development of hybrid automatic differentiation tools, standard linear algebra interfaces for numerical algorithms, time domain decomposition algorithms and two level Newton methods. We attempt to accomplish the second goal by presenting the results of several case studies in which direct sensitivities and adjoint methods have been effectively applied, in addition to an investigation of h-p adaptivity using adjoint based a posteriori error estimation. A mathematical overview is provided of direct sensitivities and adjoint methods for both steady state and transient simulations. Two case studies are presented to demonstrate the utility of these methods. A direct sensitivity method is implemented to solve a source inversion problem for steady state internal flows subject to convection diffusion. Real time performance is achieved using novel decomposition into offline and online calculations. Adjoint methods are used to reconstruct initial conditions of a contamination event in an external flow. We demonstrate an adjoint based transient solution. In addition, we investigated time domain decomposition algorithms in an attempt to improve the efficiency of transient simulations. Because derivative calculations are at the root of sensitivity calculations, we have developed hybrid automatic differentiation methods and implemented this approach for shape optimization for gas dynamics using the Euler equations. The hybrid automatic differentiation method was applied to a first order approximation of the Euler equations and used as a preconditioner. In comparison to other methods, the AD preconditioner showed better convergence behavior. Our ultimate target is to perform shape optimization and hp adaptivity using adjoint formulations in the Premo compressible fluid flow simulator. A mathematical formulation for mixed-level simulation algorithms has been developed where different physics interact at potentially different spatial resolutions in a single domain. To minimize the implementation effort, explicit solution methods can be considered, however, implicit methods are preferred if computational efficiency is of high priority. We present the use of a partial elimination nonlinear solver technique to solve these mixed level problems and show how these formulation are closely coupled to intrusive optimization approaches and sensitivity analyses. Production codes are typically not designed for sensitivity analysis or large scale optimization. The implementation of our optimization libraries into multiple production simulation codes in which each code has their own linear algebra interface becomes an intractable problem. In an attempt to streamline this task, we have developed a standard interface between the numerical algorithm (such as optimization) and the underlying linear algebra. These interfaces (TSFCore and TSFCoreNonlin) have been adopted by the Trilinos framework and the goal is to promote the use of these interfaces especially with new developments. Finally, an adjoint based a posteriori error estimator has been developed for discontinuous Galerkin discretization of Poisson's equation. The goal is to investigate other ways to leverage the adjoint calculations and we show how the convergence of the forward problem can be improved by adapting the grid using adjoint-based error estimates. Error estimation is usually conducted with continuous adjoints but if discrete adjoints are available it may be possible to reuse the discrete version for error estimation. We investigate the advantages and disadvantages of continuous and discrete adjoints through a simple example

Lawson criterion for ignition exceeded in an inertial fusion experiment

For more than half a century, researchers around the world have been engaged in attempts to achieve fusion ignition as a proof of principle of various fusion concepts. Following the Lawson criterion, an ignited plasma is one where the fusion heating power is high enough to overcome all the physical processes that cool the fusion plasma, creating a positive thermodynamic feedback loop with rapidly increasing temperature. In inertially confined fusion, ignition is a state where the fusion plasma can begin "burn propagation" into surrounding cold fuel, enabling the possibility of high energy gain. While "scientific breakeven" (i.e., unity target gain) has not yet been achieved (here target gain is 0.72, 1.37 MJ of fusion for 1.92 MJ of laser energy), this Letter reports the first controlled fusion experiment, using laser indirect drive, on the National Ignition Facility to produce capsule gain (here 5.8) and reach ignition by nine different formulations of the Lawson criterion

Perspectives on Disconnects Between Scientific Information and Management Decisions on Post-fire Recovery in Western US

Environmental regulations frequently mandate the use of "best available" science, but ensuring that it is used in decisions around the use and protection of natural resources is often challenging. In the Western US, this relationship between science and management is at the forefront of post-fire land management decisions. Recent fires, post-fire threats (e.g. flooding, erosion), and the role of fire in ecosystem health combine to make post-fire management highly visible and often controversial. This paper uses post-fire management to present a framework for understanding why disconnects between science and management decisions may occur. We argue that attributes of agencies, such as their political or financial incentives, can limit how effectively science is incorporated into decision-making. At the other end of the spectrum, the lack of synthesis or limited data in science can result in disconnects between science-based analysis of post-fire effects and agency policy and decisions. Disconnects also occur because of the interaction between the attributes of agencies and the attributes of science, such as their different spatial and temporal scales of interest. After offering examples of these disconnects in post-fire treatment, the paper concludes with recommendations to reduce disconnects by improving monitoring, increasing synthesis of scientific findings, and directing social-science research toward identifying and deepening understanding of these disconnects

Tracking the density evolution in counter-propagating shock waves using imaging X-ray scattering

We present results from time-resolved X-ray imaging and inelastic scattering on collective excitations. These data are then employed to infer the mass density evolution within laser-driven shock waves. In our experiments, thin carbon foils are first strongly compressed and then driven into a dense state by counter-propagating shock waves. The different measurements agree that the graphite sample is about twofold compressed when the shock waves collide, and a sharp increase in forward scattering indicates disassembly of the sample 1 ns thereafter. We can benchmark hydrodynamics simulations of colliding shock waves by the X-ray scattering methods employed

Tracking the density evolution in counter-propagating shock waves using imaging X-ray scattering

We present results from time-resolved X-ray imaging and inelastic scattering on collective excitations. These data are then employed to infer the mass density evolution within laser-driven shock waves. In our experiments, thin carbon foils are first strongly compressed and then driven into a dense state by counter-propagating shock waves. The different measurements agree that the graphite sample is about twofold compressed when the shock waves collide, and a sharp increase in forward scattering indicates disassembly of the sample 1 ns thereafter. We can benchmark hydrodynamics simulations of colliding shock waves by the X-ray scattering methods employed