Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “codes”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

TISCC: A Surface Code Compiler and Resource Estimator for Trapped-Ion Processors

We introduce the Trapped-Ion Surface Code Compiler (TISCC), a software tool that generates circuits for a universal set of surface code patch operations in terms of a native trapped-ion gate set. To accomplish this, TISCC manages an internal representation of a trapped-ion system where a repeating pattern of trapping zones and junctions is arranged in an arbitrarily large rectangular grid. Surface code operations are compiled by instantiating surface code patches on the grid and using methods to generate transversal operations over data qubits, rounds of error correction over stabilizer plaquettes, and/or lattice surgery operations between neighboring patches. Beyond the implementation of a basic surface code instruction set, TISCC contains corner movement functionality and a patch translation that is implemented using ion movement alone. Except in the latter case, all TISCC functionality is extensible to alternative grid-like hardware architectures. TISCC output has been verified using the Oak Ridge Quasi-Clifford Simulator (ORQCS).

LeBlond, Tyler↗
DOE OSTI2014

Code Development in Coupled PARCS/RELAP5 for Supercritical Water Reactor

The new capability is added to the existing coupled code package PARCS/RELAP5, in order to analyze SCWR design under supercritical pressure with the separated water coolant and moderator channels. This expansion is carried out on both codes. In PARCS, modification is focused on extending the water property tables to supercritical pressure, modifying the variable mapping input file and related code module for processing thermal-hydraulic information from separated coolant/moderator channels, and modifying neutronics feedback module to deal with the separated coolant/moderator channels. In RELAP5, modification is focused on incorporating more accurate water properties near SCWR operation/transient pressure and temperature in the code. Confirming tests of the modifications is presented and the major analyzing results from the extended codes package are summarized.

Hu, Po↗
DOE OSTI2015

CombAlign: a code for generating a one-to-many sequence alignment from a set of pairwise structure-based sequence alignments

In order to better define regions of similarity among related protein structures, it is useful to identify the residue-residue correspondences among proteins. Few codes exist for constructing a one-to-many multiple sequence alignment derived from a set of structure or sequence alignments, and a need was evident for creating such a tool for combining pairwise structure alignments that would allow for insertion of gaps in the reference structure. This report describes a new Python code, CombAlign, which takes as input a set of pairwise sequence alignments (which may be structure based) and generates a one-to-many, gapped, multiple structure- or sequence based sequence alignment (MSSA). The use and utility of CombAlign was demonstrated by generating gapped MSSAs using sets of pairwise structure-based sequence alignments between structure models of the matrix protein (VP40) and pre-small/secreted glycoprotein (sGP) of Reston Ebolavirus and the corresponding proteins of several other filoviruses. The gapped MSSAs revealed structure-based residue-residue correspondences, which enabled identification of structurally similar versus differing regions in the Reston proteins compared to each of the other corresponding proteins. CombAlign is a new Python code that generates a one-to-many, gapped, multiple structure- or sequence-based sequence alignment (MSSA) given a set of pairwise sequence alignments (which may be structure based). CombAlign has utility in assisting the user in distinguishing structurally conserved versus divergent regions on a reference protein structure relative to other closely related proteins. CombAlign was developed in Python 2.6, and the source code is available for download from the GitHub code repository.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

The Modeling of Advanced BWR Fuel Designs with the NRC Fuel Depletion Codes PARCS/PATHS

The PATHS (PARCS Advanced Thermal Hydraulic Solver) code was developed at the University of Michigan in support of U.S. Nuclear Regulatory Commission research to solve the steady-state, two-phase, thermal-hydraulic equations for a boiling water reactor (BWR) and to provide thermal-hydraulic feedback for BWR depletion calculations with the neutronics code PARCS (Purdue Advanced Reactor Core Simulator). The simplified solution methodology, including a three-equation drift flux formulation and an optimized iteration scheme, yields very fast run times in comparison to conventional thermal-hydraulic systems codes used in the industry, while still retaining sufficient accuracy for applications such as BWR depletion calculations. Lastly, the capability to model advanced BWR fuel designs with part-length fuel rods and heterogeneous axial channel flow geometry has been implemented in PATHS, and the code has been validated against previously benchmarked advanced core simulators as well as BWR plant and experimental data. We describe the modifications to the codes and the results of the validation in this paper.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

LLM Benchmarking with LLaMA2: Evaluating Code Development Performance Across Multiple Programming Languages

The rapid evolution of large language models (LLMs) has opened new possibilities for automating various tasks in software development. This paper evaluates the capabilities of the LLaMA 2-70B model in automating these tasks for scientific applications written in commonly used programming languages. Using representative test problems, we assess the model's capacity to generate code, documentation, and unit tests, as well as its ability to translate existing code between commonly used programming languages. Our comprehensive analysis evaluates the compilation, runtime behavior, and correctness of the generated and translated code. Additionally, we assess the quality of automatically generated code, documentation, and unit tests. Here, our results indicate that while LLaMA 2-70B frequently generates syntactically correct and functional code for simpler numerical tasks, it encounters substantial difficulties with more complex, parallelized, or distributed computations, requiring considerable manual corrections. We identify key limitations and suggest areas for future improvements to better leverage AI-driven automation in scientific computing workflows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2011

ParaDiS-FEM dislocation dynamics simulation code primer

The ParaDiS code is developed to study bulk systems with periodic boundary conditions. When we try to perform discrete dislocation dynamics simulations for finite systems such as thin films or cylinders, the ParaDiS code must be extended. First, dislocations need to be contained inside the finite simulation box; Second, dislocations inside the finite box experience image stresses due to the free surfaces. We have developed in-house FEM subroutines to couple with the ParaDiS code to deal with free surface related issues in the dislocation dynamics simulations. This primer explains how the coupled code was developed, the main changes from the ParaDiS code, and the functions of the new FEM subroutines.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Code Verification of the HIGRAD Computational Fluid Dynamics Solver

The purpose of this report is to outline code and solution verification activities applied to HIGRAD, a Computational Fluid Dynamics (CFD) solver of the compressible Navier-Stokes equations developed at the Los Alamos National Laboratory, and used to simulate various phenomena such as the propagation of wildfires and atmospheric hydrodynamics. Code verification efforts, as described in this report, are an important first step to establish the credibility of numerical simulations. They provide evidence that the mathematical formulation is properly implemented without significant mistakes that would adversely impact the application of interest. Highly accurate analytical solutions are derived for four code verification test problems that exercise different aspects of the code. These test problems are referred to as: (i) the quiet start, (ii) the passive advection, (iii) the passive diffusion, and (iv) the piston-like problem. These problems are simulated using HIGRAD with different levels of mesh discretization and the numerical solutions are compared to their analytical counterparts. In addition, the rates of convergence are estimated to verify the numerical performance of the solver. The first three test problems produce numerical approximations as expected. The fourth test problem (piston-like) indicates the extent to which the code is able to simulate a 'mild' discontinuity, which is a condition that would typically be better handled by a Lagrangian formulation. The current investigation concludes that the numerical implementation of the solver performs as expected. The quality of solutions is sufficient to provide credible simulations of fluid flows around wind turbines. The main caveat associated to these findings is the low coverage provided by these four problems, and somewhat limited verification activities. A more comprehensive evaluation of HIGRAD may be beneficial for future studies.

17 WIND ENERGY↗
DOE OSTI2012

Sensitivity Analysis of FEAST-Metal Fuel Performance Code: Initial Results

This memo documents the completion of the LANL milestone, M3FT-12LA0202041, describing methodologies and initial results using FEAST-Metal. The FEAST-Metal code calculations for this work are being conducted at LANL in support of on-going activities related to sensitivity analysis of fuel performance codes. The objective is to identify important macroscopic parameters of interest to modeling and simulation of metallic fuel performance. This report summarizes our preliminary results for the sensitivity analysis using 6 calibration datasets for metallic fuel developed at ANL for EBR-II experiments. Sensitivity ranking methodology was deployed to narrow down the selected parameters for the current study. There are approximately 84 calibration parameters in the FEAST-Metal code, of which 32 were ultimately used in Phase II of this study. Preliminary results of this sensitivity analysis led to the following ranking of FEAST models for future calibration and improvements: fuel conductivity, fission gas transport/release, fuel creep, and precipitation kinetics. More validation data is needed to validate calibrated parameter distributions for future uncertainty quantification studies with FEAST-Metal. Results of this study also served to point out some code deficiencies and possible errors, and these are being investigated in order to determine root causes and to improve upon the existing code models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

ASME Code Case Development Strategy for Ferritic/Martensitic HT-9 Steel

The US Department of Energy (DOE), Office of Nuclear Energy (NE), is supporting the development of an ASME Code Case for adoption of 12Cr-1Mo-VW ferritic/martensitic (F/M) steel, commonly known as HT-9, primarily for use in design of liquid-metal fast reactors (LMFR) and components. In 2011, Los Alamos National Laboratory (LANL) nuclear engineering staff began assisting in the development of a small modular reactor (SMR) design concept, previously known as the Hyperion Module, now called the Gen4 Module, being developed by Gen4 Energy, Inc. located in Denver, Colorado. LANL engineering staff had initially proposed either T91 (9Cr-1Mo-VNb) or HT-9 for the reactor vessel and components, as well as fuel clad and ducting. Although T91 is currently an approved/authorized material for elevated temperature design in the ASME Code, the preference is utilization of HT-9 due to its extensive irradiation database and similar creep resistance as T91. In addition, HT-9 has relatively good corrosion properties for fast reactor environment, such as liquid lead bismuth eutectic coolant. Although the ASME material Code Case, for adoption of HT-9 as an approved elevated temperature material for LMFR service, is the ultimate goal of this project, several key deliverables have been accomplished in FY2012 to provide assurance to DOE/NE of the viability and feasibility of pursuing this goal. An important key deliverable was the recently completed, Gap Analysis of Material Properties Data for Ferritic/Martensitic HT-9 (1). This report provided a cursory research, accumulation, and documentation of specific material parameters; physical, mechanical, and environmental, which are necessary for ensuring a successful ASME Code Case. Time-independent tensile and ductility data and time-dependent creep and creep-rupture behavior are some of the material properties required for a successful ASME Code case. Unfortunately, results of the Gap Analysis report show that many of the material properties have limited data availability in open-source documentation, or are contained in Applied Technology documents, or there is insufficient data to support the ASME requirement from 3-separate heats of material.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Implementation of a coupled treatment of the one-turn mapping in the ray-tracing code Zgoubi

The spin is an intrinsic property of elementary particles and is considered as a quantum angular momentum. It is defined by a direction which is either up or down and an absolute value which is proportional to the elementary quantum of action. It has no analogue in classical mechanics but, as the angular momentum in classic mechanics, the spin of composite systems such as the proton is a combination of its spin components. Although knowing the proton components and their spin, their contribution to the proton spin is unknown and the physicists are still trying to unravel the mystery. A recent theoretical work suggests that the contribution carried by the gluons may be large but it has to be confirmed experimentally. The main accelerator facility which is looking for answers about the origin of the proton spin is located at the Brookhaven National Laboratory (Upton, New York, United States of America). This laboratory has the biggest accelerator complex dedicated to the colliding of polarized proton beam. The protons are first supplied by the Optically Pumped Polarized Ion Source (OPPIS), followed by a 200-MeV Linac and are then transferred to the Booster synchrotron. The protons go into the Alternating Gradient Synchrotron (AGS) which accelerate and prepair them to be injected into the Relativistic Heavy Ion Collider (RHIC). Finally the proton beam is stored into the RHIC ring for several hours and supplies collisions in two detectors of particle experiments placed. The data collected are used for studying the decay products of polarized proton collisions. Although the current RHIC spin experiments are already providing the first information about the proton's spin origin, the complex RHIC is still under improvement. Indeed it is working at minimizing undesirable effects which damage the beam polarization along the acceleration. The polarized proton beam is prone to depolarization due to spin resonances. They could be partially avoided by improving the control of some beam parameters. For instance coupling between the two transverse motions occurs in AGS and modifies the beam parameters predicted by the simulations. The locations of all the coupling sources in the AGS are not really known which makes their removal or compensation difficult to elaborate. This requires improving the simulations of beam dynamics made in the AGS. They need to be performed with, on one hand, a very accurate model and, on the other hand, a code capable of analyzing a coupled motion. Two codes are mainly used: MADX and a ray-tracing code called ZGOUBI 4 which is developing by Dr. François MÉOT. It is scheduled to install ZGOUBI onto the computers of the RHIC complex main control room in order to facilitate and extend its use. Nevertheless ZGOUBI is not completely able to deal with coupled one-turn mapping produced by the ray-tracing so it is necessary to find a way of treating the coupled systems and to implement it within the existing code, that will allow producing coupled optical functions and transporting them using the ray-tracing method. The purpose of this report is to present the results of my research work performed during this 5-month internship. It mainly consisted in searching and synthesizing the appropriate formalism in order to fully treat the coupled which outcomes from ZGOUBI. It required to nd a way of implementing it into the existing ZGOUBI code then to test it. Finally after validation of these add-ons some simulations have been handled by ZGOUBI for the AGS.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

TART2012 An Overview of A Coupled Neutron Photon 3-D, Combinatorial Geometry Time Dependent Monte Carlo Transport Code

This document is designed to provide an overview of TART2012, a coupled neutron photon, 3 Dimensional, combinatorial geometry, time dependent Monte Carlo radiation transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART2012 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

Verification and validation of the PLTEMP/ANL code for thermal hydraulic analysis of experimental and test reactors

The document compiles in a single volume several verification and validation works done for the PLTEMP/ANL code during the years of its development and improvement. Some works that are available in the open literature are simply referenced at the outset, and are not included in the document. PLTEMP has been used in conversion safety analysis reports of several US and foreign research reactors that have been licensed and converted. A list of such reactors is given. Each chapter of the document deals with the verification or validation of a specific model. The model verification is usually done by comparing the code with hand calculation, Microsoft spreadsheet calculation, or Mathematica calculation. The model validation is done by comparing the code with experimental data or a more validated code like the RELAP5 code.

Kalimullah, M.↗
DOE OSTI2013

Using Third-Party Inspectors in Building Energy Codes Enforcement in India

India is experiencing fast income growth and urbanization, and this leads to unprecedented increases in demand for building energy services and resulting energy consumption. In response to rapid growth in building energy use, the Government of India issued the Energy Conservation Building Code (ECBC) in 2007, which is consistent with and based on the 2001 Energy Conservation Act. ECBC implementation has been voluntary since its enactment and a few states have started to make progress towards mandatory implementation. Rajasthan is the first state in India to adopt ECBC as a mandatory code. The State adopted ECBC with minor additions on March 28, 2011 through a stakeholder process; it became mandatory in Rajasthan on September 28, 2011. Tamil Nadu, Gujarat, and Andhra Pradesh have started to draft an implementation roadmap and build capacity for its implementation. The Bureau of Energy Efficiency (BEE) plans to encourage more states to adopt ECBC in the near future, including Haryana, Uttar Pradesh, Karnataka, Maharashtra, West Bengal, and Delhi. Since its inception, India has applied the code on a voluntary basis, but the Government of India is developing a strategy to mandate compliance. Implementing ECBC requires coordination between the Ministry of Power and the Ministry of Urban Development at the national level as well as interdepartmental coordination at the state level. One challenge is that the Urban Local Bodies (ULBs), the enforcement entities of building by-laws, lack capacity to implement ECBC effectively. For example, ULBs in some states might find the building permitting procedures to be too complex; in other cases, lack of awareness and technical knowledge on ECBC slows down the amendment of local building by-laws as well as ECBC implementation. The intent of this white paper is to share with Indian decision-makers code enforcement approaches: through code officials, third-party inspectors, or a hybrid approach. Given the limited capacity and human resources available in the state and local governments, involving third-party inspectors could rapidly expand the capacity for plan reviews and broad implementation. However, the procedures of involving third-parties need to be carefully designed in order to guarantee a fair process. For example, there should be multiple checks and certification requirements for third-party inspectors, and the government should have the final approval when third-party inspectors are used in a project. This paper discusses different approaches of involving third-parties in ECBC enforcement; the Indian states may choose the approaches that work best in their given circumstances.

Yu, Sha↗
DOE OSTI2012

Integral and Separate Effects Tests for Thermal Hydraulics Code Validation for Liquid-Salt Cooled Nuclear Reactors

The objective of the 3-year project was to collect integral effects test (IET) data to validate the RELAP5-3D code and other thermal hydraulics codes for use in predicting the transient thermal hydraulics response of liquid salt cooled reactor systems, including integral transient response for forced and natural circulation operation. The reference system for the project is a modular, 900-MWth Pebble Bed Advanced High Temperature Reactor (PB-AHTR), a specific type of Fluoride salt-cooled High temperature Reactor (FHR). Two experimental facilities were developed for thermal-hydraulic integral effects tests (IETs) and separate effects tests (SETs). The facilities use simulant fluids for the liquid fluoride salts, with very little distortion to the heat transfer and fluid dynamics behavior. The CIET Test Bay facility was designed, built, and operated. IET data for steady state and transient natural circulation was collected. SET data for convective heat transfer in pebble beds and straight channel geometries was collected. The facility continues to be operational and will be used for future experiments, and for component development. The CIET 2 facility is larger in scope, and its construction and operation has a longer timeline than the duration of this grant. The design for the CIET 2 facility has drawn heavily on the experience and data collected on the CIET Test Bay, and it was completed in parallel with operation of the CIET Test Bay. CIET 2 will demonstrate start-up and shut-down transients and control logic, in addition to LOFC and LOHS transients, and buoyant shut down rod operation during transients. Design of the CIET 2 Facility is complete, and engineering drawings have been submitted to an external vendor for outsourced quality controlled construction. CIET 2 construction and operation continue under another NEUP grant. IET data from both CIET facilities is to be used for validation of system codes used for FHR modeling, such as RELAP5-3D. A set of numerical models were developed in parallel to the experimental work. RELAP5-3D models were developed for the salt-cooled PB-AHTR, and for the simulat fluid CIET natural circulation experimental loop. These models are to be validated by the data collected from CIET. COMSOL finite element models were used to predict the temperature and fluid flow distribution in the annular pebble bed core; they were instrumental for design of SETs, and they can be used for code-to-code comparisons with RELAP5-3D. A number of other small SETs, and numerical models were constructed, as needed, in support of this work. The experiments were designed, constructed and performed to meet CAES quality assurance requirements for test planning, implementation, and documentation; equipment calibration and documentation, procurement document control; training and personnel qualification; analysis/modeling software verification and validation; data acquisition/collection and analysis; and peer review.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

User's Manual for the Hydrate v1.5 Option of TOUGH+ v1.5: A Code for the Simulation of System Behavior in Hydrate-Bearing Geologic Media

HYDRATE v1.5 is a numerical code that for the simulation of the behavior of hydrate-bearing geologic systems, and represents the third update of the code since its first release [Moridis et al., 2008]. It is an option of TOUGH+ v1.5 [Moridis, 2014], a successor to the TOUGH2 [Pruess et al., 1999] family of codes for multi-component, multiphase fluid and heat flow developed at the Lawrence Berkeley National Laboratory. HYDRATE v1.5 needs the TOUGH+ v1.5 core code in order to compile and execute. It is written in standard FORTRAN 95/2003, and can be run on any computational platform (workstation, PC, Macintosh) for which such compilers are available. By solving the coupled equations of mass and heat balance, the fully operational TOUGH+HYDRATE code can model the non-isothermal gas release, phase behavior and flow of fluids and heat under conditions typical of common natural CH4-hydrate deposits (i.e., in the permafrost and in deep ocean sediments) in complex geological media at any scale (from laboratory to reservoir) at which Darcy’s law is valid. TOUGH+HYDRATE v1.5 includes both an equilibrium and a kinetic model of hydrate formation and dissociation. The model accounts for heat and up to four mass components, i.e., water, CH4, hydrate, and water-soluble inhibitors such as salts or alcohols. These are partitioned among four possible phases (gas phase, liquid phase, ice phase and hydrate phase). Hydrate dissociation or formation, phase changes and the corresponding thermal effects are fully described, as are the effects of inhibitors. The model can describe all possible hydrate dissociation mechanisms, i.e., depressurization, thermal stimulation, salting-out effects and inhibitor-induced effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

The Cost of Enforcing Building Energy Codes: Phase 2

The purpose of this study is to present key findings regarding costs associated with enforcing building energy code compliance–primarily focusing on costs borne by local government. Building codes, if complied with, have the ability to save a significant amount of energy. However, energy code compliance rates have been significantly lower than 100%. Renewed interest in building energy codes has focused efforts on increasing compliance, particularly as a result of the 2009 American Recovery and Reinvestment Act (ARRA) requirement that in order for states to receive additional energy grants, they must have “a plan for the jurisdiction achieving compliance with the building energy code…in at least 90 percent of new and renovated residential and commercial building space” by 2017 (Public Law 111-5, Section 410(2)(C)). One study by the Institute for Market Transformation (IMT) estimated the costs associated with reaching 90% compliance to be $810 million, or $610 million in additional funding over existing expenditures, a non-trivial value. [Majersik & Stellberg 2010] In this context, Lawrence Berkeley National Laboratory (LBNL) conducted a study to better pinpoint the costs of enforcement through a two-phase process.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2015

The Initial Atmospheric Transport (IAT) Code: Description and Validation

The Initial Atmospheric Transport (IAT) computer code was developed at Sandia National Laboratories as part of their nuclear launch accident consequences analysis suite of computer codes. The purpose of IAT is to predict the initial puff/plume rise resulting from either a solid rocket propellant or liquid rocket fuel fire. The code generates initial conditions for subsequent atmospheric transport calculations. The Initial Atmospheric Transfer (IAT) code has been compared to two data sets which are appropriate to the design space of space launch accident analyses. The primary model uncertainties are the entrainment coefficients for the extended Taylor model. The Titan 34D accident (1986) was used to calibrate these entrainment settings for a prototypic liquid propellant accident while the recent Johns Hopkins University Applied Physics Laboratory (JHU/APL, or simply APL) large propellant block tests (2012) were used to calibrate the entrainment settings for prototypic solid propellant accidents. North American Meteorology (NAM )formatted weather data profiles are used by IAT to determine the local buoyancy force balance. The IAT comparisons for the APL solid propellant tests illustrate the sensitivity of the plume elevation to the weather profiles; that is, the weather profile is a dominant factor in determining the plume elevation. The IAT code performed remarkably well and is considered validated for neutral weather conditions.

97 MATHEMATICS AND COMPUTING↗