skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Optimization of Geometric Perturbations on a Rod Moving Through a High Explosive Target

Technical Report ·
DOI:https://doi.org/10.2172/1633518· OSTI ID:1633518
 [1]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

After completing a study to ensure the simulation results were converged, several high resolution 3D Smoothed Particle Hydrodynamic (SPH) simulations of copper rods impacting a high explosive (LX14) target were performed. This was then formulated into an optimization problem: I wanted to find the optimum shape and location of a perturbation on the rod that would maximize its erosion after it left the target. The shape of the perturbation was modeled as a 2D Gaussian bump and parameterized by its location along the rod axis (z0) and amplitude (A). The final mass of the coherent part of the rod as it leaves the target was used as a metric to represent the erosion of the rod, and the optimization was formulated to maximize this metric with respect to the aforementioned design variables. Due to the expensive nature of the high-fidelity 3D SPH simulations, a surrogate model needed to be chosen so that many function calls to the optimizer would be feasible. Thus, a strategic full factorial sampling plan was chosen to build a dataset, which consisted of 24 high-fidelity simulations. Two surrogate models, a third order polynomial regression model and a Gaussian Process Model, were analyzed using a 14%/86% test/train holdout technique. The root mean square and R2 score of the test set was used to determine the best model, and the third order polynomial regression model was chosen as the surrogate model. Finally, the Nelder-Mead Simplex and Basin-hopping optimization algorithms were implemented, and it was found that the two algorithms gave slightly different optimum values. Nelder-Mead gave an optimum point of [z*0;A*] = [9:9;0:4] and Basin-Hopping gave an optimum value of x* = [z*0;A*] = [9:2;0:1].

Research Organization:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA)
DOE Contract Number:
AC52-07NA27344
OSTI ID:
1633518
Report Number(s):
LLNL-TR-811405; 1018187; TRN: US2200872
Country of Publication:
United States
Language:
English

Similar Records

Evaluating Optimization Strategies for Engine Simulations Using Machine Learning Emulators
Journal Article · Thu Jun 20 00:00:00 EDT 2019 · Journal of Engineering for Gas Turbines and Power · OSTI ID:1633518

Gaussian processes with built-in dimensionality reduction: Applications to high-dimensional uncertainty propagation
Journal Article · Thu Sep 15 00:00:00 EDT 2016 · Journal of Computational Physics · OSTI ID:1633518

Co-Optimization of CO2-EOR and Storage Processes in Mature Oil Reservoirs
Journal Article · Tue Aug 02 00:00:00 EDT 2016 · Greenhouse Gases: Science and Technology · OSTI ID:1633518