Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time integrator”

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,459 records · Page 81

WVU’s contribution within the ZAPP collaboration on the Z Facility at Sandia National Laboratories: Final Technical Report

In Year 1, photoionized silicon spectroscopy experiments were designed and executed, with charge-state distribution and ionic column densities comparable to the circumsource environments of accretion-powered black holes. Charge-state distribution was measured with absorption spectroscopy, using the Z-pinch itself as backlight, and the photoionizing radiation flux was determined with x-ray power measurements, accounting appropriately for geometrical dilution. The emergent emission spectrum was measured with a high-sensitivity time-integrated spectrometer. The emission has been measured from plasmas with varying size to elucidate the radiation transport effects. Data and knowledge of the controlled experimental conditions were provided to collaborators who used different atomic kinetics spectral synthesis codes to predict the emission spectrum. In Years 2 and 3, multi-element (in layers) foil was used for diagnostic development. A framing mask was implemented into the configuration to aperture the radiation incident on the foil. The foil distance from the z-pinch was adjusted to achieve designed electron temperature and ultimately mounted directly on the Return-Current Cannister. A much higher ionization state was achieved, which is desirable for diagnostic development. In Year 4, significant data analysis took place. Astrophysics related projects were introduced "Seeking a full theory of Stark broadening at HED in the presence of crossed E and B fields for estimating surface gravity from line profiles in Magnetic White Dwarf (MWD) stars" and "Plasma collective effects in atomic kinetics calculations of lineshape, line position, and line amplitude" for inclusion into the ZAPP collaboration, but were declined by ZAPP leadership. During the grant period, WVU contributed to three proposals for a ZAPP Center of Excellence that were prepared and/or submitted in part or in full. (1) NSF Physics Frontier Center solicitation 16-561, PI / Director: A. Wootton, UTexas, "Center for experimental astrophysic"s (declined), (2) NSF Science and Technology Center solicitation 2017, PI/Director: D. Winget, UTexas, "Formation and Evolution of Planets, Stars and Galaxies and Extreme-Parameter Laboratory Astrophysics" (awaiting restoration of the solicitation after unexpected postponement), and (3) NNSA Stewardship Science Grant Program solicitation 81-112, PI/Director: D. Winget, UTexas, "Center for Astrophysical Plasma Properties" (submitted in 2017, awarded in 2018). Based on the work performed during the grant period, a new NNSA grant application to support a WVU project in ZAPP, "Spectroscopic methods to obtain electron density and electron temperature applied to soft x-ray absorption spectra from radiatively heated Z-pinch plasmas", was submitted in 2017 and funded in 2018.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ventilation and Indoor Air Quality in New California Homes with Gas Appliances and Mechanical Ventilation

Substantial energy is used to condition the air that enters California homes through leaks in the building envelope and ductwork - typically about a third of all heating and cooling. Reducing this through air sealing is essential to California achieving zero energy homes. However, this outdoor air also dilutes pollutants emitted inside homes and contributes to a healthy indoor environment and acceptable indoor air quality (IAQ). To address this IAQ issue, California’s Title 24 Building Standards have required mechanical ventilation in new homes since 2008. This report presents a comprehensive study of the impacts of these requirements in recently constructed homes with natural gas appliances. The study included a survey about satisfaction and activities that impact IAQ; a field study of homes built to 2008 or later; and simulations assessing how various ventilation rates would impact chronic exposures to an indoor emitted pollutant as air tightness improves in California. The report focuses on the field study; the webbased survey and simulation elements are described in appendices. The field study characterized 70 homes built between 2011 and 2017. Each home was monitored over roughly one week with the dwelling unit mechanical ventilation system operating and windows closed. Pollutant measurements included time-resolved fine particulate matter (PM 2.5 ) indoors and outdoors, and formaldehyde, nitrogen dioxide (NO 2 ), and carbon dioxide (CO 2 ) indoors. Time-integrated measurements were made for formaldehyde, NO 2 and nitrogen oxides (NO X ) indoors and outdoors at all homes. Activity monitoring devices were installed on the cooktop, range hood and other exhaust fans, and the heating and cooling system. The field study found that most homes met most ventilation requirements and the dwelling unit ventilation fans on average moved 50% more airflow than the minimum specified in Title 24. Air pollutant concentrations were similar or lower than those reported in a study of recent construction California new homes conducted in 2007-08. Notably, the median formaldehyde level was 38% lower than in the prior study. Measured concentrations were below health guidelines for most pollutants, indicating that IAQ is acceptable in new California homes when dwelling unit mechanical ventilation is used. However, the dwelling unit mechanical ventilation fans were only operating in one quarter of the homes when first visited and the control switches in many homes did not have informative labels as required by the standards. Corrective action needs to be taken to improve labeling and controls for ventilation systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Non-Hydrostatic Variable Resolution Atmospheric Model in ACME (Final Technical Report)

The project goals are divided into five main tasks to be accomplished over three years: Task 1: Implement a nonhydrostatic atmospheric dynamical core in ACME (via HOMME) Task 2: Improve load balancing in variable resolution simulations by extending Zoltan. Task 3: Optimize implicit-explicit time integrators for overcoming the vertical acoustic wave CFL limitation Task 4: Validate the NH dynamical core with DCMIP and AMIP tests on uniform and variable resolution grids Task 5: Apply the NH dynamical core to study and parameterize turbulent mountain stress and orographic gravity wave drag with ultra-high variable resolution grids This report only concerns contributions from the UC Davis team, which includes Co-PI Paul Ullrich, GSR/Postdoc Jorge Guerra, GSR Alan Rhoades, GSR Shiheng Duan, and GSR Jennifer Piaseczny. From the above list, the UC Davis team is responsible for contributions to Tasks 1, 3, 4 and 5.

54 ENVIRONMENTAL SCIENCES↗

High Energy Micron Scale Pixel Hybrid Detector

The objective of this technology development project has been to fabricate and characterize a new X-ray area image sensor innovation designed for high efficiency detection and imaging of hard X-rays (>20 KeV) with micron scale pixel resolution. A key technology component consists of a monolithic hybrid X-ray detector built by layering an amorphous Selenium X-ray photoconductor film directly deposited on a small pitch (7.8μm) CMOS active pixel readout array. A specific goal of this work is to provide a new detector tool to researchers working with synchrotron light sources, to extend the depth of exploration into nanoscale structures, and do so cost effectively. Early generation X-ray imaging detector tools fail to provide adequate resolving power and efficiency at high energies and our innovation provides a solution. A majority of proposed phase I tasks were completed. Those tasks included, in collaboration with KA Imaging, successful fabrication of specified a-Se layers on a 1Mpixel CMOS readout array to create the direct hybrid X-ray imaging array, packaging of the imaging array, assembly of electronics into a shielded detector box, loading user interface software for camera control and image capture, and performing initial (no X-ray) detector evaluation at Farrier Microengineering facilities. After initial tests, extensive testing using the X-ray beam line 1-BM source at Argonne National Laboratory Advanced Photon Source(ANL-APS) was performed. Imaging performance studies were conducted at multiple X-ray beam energies. Studies included X-ray response vs. X-ray energy, vs. applied detector field strength (HV), response linearity vs. integration time, MTF vs. beam energy, image lag measurements, and micron scale resolution pattern target imaging. This high energy micron scale X-ray detection technology is a first of a kind to be evaluated at the ANL-APS and is supported by Dr. Antonino Miceli, head of the detector physics group. Accomplishments include the first successful demonstration of a 1Mpix amorphous Selenium on CMOS readout direct hybrid X-ray imaging detector at a DOE sponsored National Laboratory. Modeled performance predictions were confirmed, proving micron scale resolution as well as an order of magnitude increase in conversion efficiency over conventional scintillator detector performance. A discovery of evidence of propagation-based phase contrast edge enhancement was unplanned and will be exploited in future experiments. This phase 1 technology demonstration verifies the efficacy of proposed high energy micron scale pixel X-ray hybrid a-Se imaging array products. The new detector products will cost effectively replace early generation detector technologies that currently fail to address increasing demand for tools to study nano structures and physical processes at high X-ray energies above 20keV. For example, current systems for transmission X-ray microscopy have low efficiency and available large pixel (~55µm) direct detectors do not offer sufficient resolution for Bragg fringe imaging at high energies. This detector innovation will solve these inadequacies. The new detector technology will compete in a $50M scientific X-ray detector market, growing at 6%, and will benefit Microtomography(µCT), X-ray diffractometry (XRD), Coherent Diffraction Imaging (CDI) applications. One application specified by Dr. Miceli at ANL is the 3D reconstruction of compact crystals at high energies requires small pixel high efficiency X-ray detectors. Farrier Microengineering, working with colleagues at KA Imaging and the University of Waterloo, successfully achieved the critical objectives of the Phase I proposal. We have proven the technical feasibility of this innovation by building and successfully testing a 1Mpixel a-Se/CMOS hybrid direct X-ray imaging array. Our commercialization roadmap of detector products will be implemented during a phase II project. We will demonstrate a migration of the technology to larger detector array sizes, improved temporal performance with detector layer engineering, and improved environmental and operational stability. Manufacturability and cost effectiveness are key objectives and product designs will leverage current state of the art CMOS wafer fabrication quality and volume manufacturing economics.

36 MATERIALS SCIENCE↗

Interpenetrating Plasma Simulations

This project investigated multifluid models for simulating interpenetrating plasmas. These models solve the Euler equations for each plasma stream along with an appropriate model for electron dynamics. The project developed a novel code, Euclid, that simulates these models using high-order conservative spatial discretization methods and explicit and semi-implicit time integration methods. Euclid was used to solve several multifluid test cases representative of plasma dynamics in high-energy density physics experiments. Kinetic simulations of interpenetrating flows using the LOKI code revealed novel physical effects of the ion-ion streaming instability. A new multi-flow collision operator was developed for the LOKI code. A simple model of the ion-ion streaming instability was also developed. Several presentations were given on the project results; three journal articles have been published, one has been submitted, and two are planned.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Foams in Hohlraums

Foam materials are starting to find application in laser-heated hohlraums used to drive inertial confinement fusion (ICF) implosions. The disparate scales of the solid foam elements (nm) compared to the typical problem scale (mm) make it expensive to resolve the individual foam elements in hydrodynamic calculations. In this project we developed two new models for treating foam materials – a 3D hydrodynamic model that treats the foam as bits of pre-expanded solid material separated by voids, and a subgrid model that captures the essential physics of laser-heated foams without requiring resolution of the individual solid elements. We also explored the microphysics of the transition from interpenetrating kinetic behavior to fully-collisional fluid-like behavior during the early time heating and expansion of foam elements through a combination of single-fluid and multi-fluid hydrodynamics calculations and particle-in-cell calculations. The model was benchmarked against these calculations and against experimental data on laser-heated chemical and additive-manufactured (AM) foams. Since there was no existing data on laser-heated AM foams, we performed the first experiments on such foams at the Jupiter Laser Facility. Samples of three different types of printed AM foams were heated using a single 0.5 micron laser. The foam densities ranged from 10 to 100 mg/cc (supercritical for 0.5 micron light). We measured the backscattered light (power and spectrum), the transmitted light, and the time-integrated temperature of the foam samples. This benchmarked subgrid foam model is a new capability that should allow routine calculation of ICF designs using foam elements, potentially opening up a new design space.

36 MATERIALS SCIENCE↗

Coarse-Grained Reactive Molecular Dynamics Simulations of Heterogeneities in Shocked Energetic Materials (LDRD Final Report)

This report summarizes the result of LDRD project 16-0161, titled "Coarse-Grained Re- active Molecular Dynamics Simulations of Heterogeneities in Shocked Energetic Materials." The purpose of the project was to develop a coarse-grained reactive molecular dynamics capability in LAMMPS enabling simulations of initiation in energetic materials comparable in accuracy to what is currently possible using large-scale reactive molecular dynamics, but with greatly reduced computational cost. The starting point for this work was the reactive dissipative particle dynamics (DPD) approach, which has been implemented as the new USER-DPD package in LAMMPS by researchers at Army Research Laboratory. Using modified versions of the examples provided with the new package, we examined the computational efficiency of the method, as well as its ability to model energy release in energetic materials. We observed that the Shardlow splitting method provides a great speed and accuracy advantage over conventional velocity Verlet time integration. We observed that the generic model of an energetic material provided with the USER-DPD package exhibited incomplete reaction under constant volume conditions. This was caused by quenching of the internal temperature of the molecules due to a rapid build-up of repulsive interactions between product gas components under constant volume conditions. Under constant pressure conditions, complete reaction was observed, as volume expansion prevented the buildup of strong repulsive interactions. Finally, a more realistic model calibrated to reproduce the equation of state of the RDX molecular crystal was examined. This model exhibited much less quenching of the internal temperature under constant volume conditions and reacted very rapidly under constant pressure conditions.

74 ATOMIC AND MOLECULAR PHYSICS↗

TEAM project review, year 1

This note summarizes our research activities within the TEAM project between October 2019 and November 2020, funded by the ASCR Advanced Research in Quantum Computing program. Three research efforts have begun during the review period. The first effort, resulting in the JuQBox software, implemented a novel algorithm for optimal control in closed quantum systems based on symplectic time integration and exact calculations of the gradient of the discretized objective function. In our second effort, resulting in a code called Quandary, we consider the optimal control problem for open quantum systems. Here the evolution of the quantum system is modeled by the Lindblad master equation. In this code, we are considering both the realization of quantum gates as well as the optimal reset problem in which a general mixed quantum state is driven towards its ground state. In our third research effort we have begun the development of computational software for what we call Lindblad learning, i.e. the process of characterizing the equations governing our quantum computing hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling the Vertical Response of an Aeroseismometer Balloon System

Aeroseismometery is a novel, cutting edge capability that involves balloon based systems for detecting and geolocating sources of infrasound. The incident infrasound from a range of sources such as volcanos, earthquakes, explosions, supersonic aircraft impinges upon the balloon system causing it to respond dynamically. The dynamic response is post-processed to locate the infrasound source. This report documents the derivation of an analytical model that predicts the balloon dynamics. Governing equations for the system are derived as well as a transfer function relating the infrasound signal to the net force on the balloon components. Experimental measurements of the infrasound signals are convolved with the transfer function and the governing equations numerically time integrated to obtain predictions of the displacement, velocity and acceleration of the balloon system. The predictions are compared to the experimental measurements with good agreement observed. The derivation focuses only on the vertical dynamics of the balloon system. Future work will develop governing equations for the swinging response of the balloon to the incident infrasound.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

(U) Marker Lagrangian Simulations using Geometry from CAD Assembly Files [Slides]

The Eulerian code PAGOSA now contains an option for the forward time integration of any user selected material in the Lagrangian frame of reference. This option is accomplished using a Multifield version of Brackbill’s FLIP plus Sulsky’s MPM. FLIP+MPM is a Marker Lagrangian scheme that is not to be confused with well known Grid Lagrangian schemes. This Marker Lagrangian method allows materials to undergo completely arbitrary deformation, up to and including fracture and separation. Using this option makes PAGOSA into a mixed frame simulation tool such that fluids can be integrated in the Eulerian frame and solids integrated in the Lagrangian frame. Boundaries between the fields are said to be ’immersed’ in the computational domain, and their nature is formulated using Multifield Theory. Each Marker in a field of FLIP+MPM material tracks the thermodynamic state of a small (but finite) piece of mass. Initiation of the Markers using the geometry and material information from CAD Assembly files is the subject of the talk. The ease and speed with which the dynamics of the assembly can be studied is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

ALP Experiment Sensitivity Exploration (Preliminary Findings)

This report presents an exploration of sensitivity estimates for multiple prominent proposed techniques to search for axion-like particles (ALPs). The sensitivities are plotted over particle mass ($m_a$) and coupling ($g_{aγγ}$) parameter space. Explorations include variation over detector type, integration time, background cosmology, and size of the apparatus.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High temperature inelastic constitutive models for the ASME Section III, Division 5 Class A materials

This report describes the exploration of a universal high temperature inelastic constitutive model for use with the ASME Boiler & Pressure Vessel Code Section III, Division 5, Class A design rules. The idea, developed based on feedback from reactor vendors, is to simplify the current bespoke material models for Grade 91, 316H, and Alloy 617 and the new model for Alloy 800H into a simple, single model form. The report describes a new parallel time integration technique implemented in the pyoptmat package which allowed us to explore a wide variety of model forms, searching for a suitable common model. The report then describes preliminary work on models for monontonic deformation and then the development of a set of models suitable for capturing high temperature cyclic deformation. These models are suitable for use with the ASME design rules, except potentially for a narrow, material-specific band of temperatures at the edge of the creep regime. An incremental improvement the current models could overcome this discrepancy and provide a new, simpler set of constitutive model for all four materials.

36 MATERIALS SCIENCE↗

Ventilation and Indoor Air Quality in Recently Constructed U.S. Homes: Measured Data From Select Southeastern States

Little data is available that quantifies IAQ or associates contaminant levels with air exchange in U.S. homes. To address this gap, the U.S. Department of Energy (DOE) Building America Program conducted a study to characterize IAQ in U.S. homes constructed since 2013, along with presence, functionality, and occupant use of control measures. Specific objectives of this study effort included: Measuring time-integrated concentrations and temporal profiles of established contaminants of concern; monitoring the use of ventilation equipment; and tracking activities that impact air pollutant emissions and removal processes in typical homes in various climate zones; Characterizing the prevalence, type, and installed performance of mechanical ventilation equipment in new homes; exploring regional variations in system designs and performance; Investigating associations between contaminant levels and the presence of control measures including whole-house mechanical ventilation (WHMV). The study engaged research teams to collect data at a regional level. This report introduces the study protocol and presents and discusses high-level results obtained from Florida Solar Energy Center (FSEC) data collection efforts in the southeastern U.S. states of Florida, Georgia, and South Carolina.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

PETSc/TAO Users Manual V.3.21

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations (PDEs) and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for implementing large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and optimizers that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started ). The library is organized hierarchically, enabling users to employ the abstraction level most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users.

97 MATHEMATICS AND COMPUTING↗

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range

This assessment uses screening level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program Mixing-Cell Model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS).

99 GENERAL AND MISCELLANEOUS↗

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range

This assessment uses screening level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface soil impacts were calculated with the computer program Mixing-Cell Model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS).

99 GENERAL AND MISCELLANEOUS↗

Development of High-Fidelity and Efficient Modeling Capabilities for Enabling Co-Optimization of Fuels and Multi-Mode Engines

The objective of this closely coordinated research project is to develop improved physical sub-models and innovative numerical algorithms to enable accurate and efficient predictions of multi-mode combustion to support EERE’s Co-Optima program. To this end, a research pro-gram is proposed that seeks • to develop accurate submodels for predicting multi-mode combustion regimes, wall-heat transfer, non-equilibrium plasma ignition, and combustion-mode transition, • to develop dynamic load-rebalancing algorithms and efficient semi-implicit time-integration schemes for enabling exascale computing, • to validate the computational submodels against experiments and direct numerical simulation (DNS) data in conjunction with sensitivity analysis, and • to support the direct transition of the resulting sub-models into other engine-simulation tools through the development of self-contained modules. All sub-models will be integrated into the open-source exascale engine-simulation platform Nek5000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Theory and Implementation of the Spectacular Nonlinear Viscoelastic Constitutive Model

This report is a comprehensive guide to the nonlinear viscoelastic Spectacular model, which is an isotropic, thermo-rheologically simple constitutive model for glass-forming materials, such as amorphous polymers. Spectacular is intermediate in complexity to the previous PEC and SPEC models (Potential Energy Clock and Simplified Potential Energy Clock models, respectively). The model form consists of two parts: a Helmholtz free energy functional and a nonlinear material clock that controls the rate of viscoelastic relaxation. The Helmholtz free energy is derived from a series expansion about a reference state. Expressions for the stress and entropy functionals are derived from the Helmholtz free energy following the Rational Mechanics approach. The material clock depends on a simplified expression for the potential energy, which itself is a functional of the temperature and strain histories. This report describes the thermo-mechanical theory of Spectacular, the numerical methods for time-integrating the model, model verification for its implementation in LAMÉ, a user guide for its implementation in LAMÉ, and ideas for future work. A number of appendices provide supplementary mathematical details and a description of the procedure used to derive the simplified potential energy from the full expression for the potential energy. The goal of this report is create a convenient point-of-entry for engineers who wish to learn more about Spectacular, but also to serve as a reference manual for advanced users of the model.

36 MATERIALS SCIENCE↗