Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A method for mapping submicron-scale crystallographic order/disorder applied to human tooth enamel

Tooth enamel, the outermost layer of human teeth, is a complex, hierarchically structured biocomposite. The details of this structure are important in multiple human health contexts, from understanding the progression of dental caries (tooth decay) to understanding the process of amelogenesis and related developmental defects. Enamel is composed primarily of long, nanoscale crystallites of hydroxyapatite that are bundled by the thousands to form micron-scale rods. Studies with transmission electron microscopy show the relationships between small groups of crystallites and X-ray diffraction characterize averages over many rods, but the direct measurement of variations in local crystallographic structure across and between enamel rods has been missing. Here, we describe a synchrotron X-ray-based experimental approach and a novel analysis method developed to address this gap in knowledge. A ~500-nm-wide beam of monochromatic X-rays in conjunction with a sample section only 1 μm in thickness enables 2D diffraction patterns to be collected from small well-separated volumes within the enamel microstructure but still probes enough crystallites (~300 per pattern) to extract population-level statistics on crystallographic features like lattice parameter, crystallite size, and orientation distributions. Furthermore, the development of a quantitative metric to characterize relative order and disorder based on the azimuthal autocorrelation of diffracted intensity enables these crystallographic measurements to be correlated with their location within the enamel microstructure (e.g., between rod and interrod regions). These methods represent a step forward in the characterization of human enamel and will elucidate the variation of the crystallographic structure across and between enamel rods for the first time.

crystallographic order/disorder↗
DOE OSTI2020

Assessment of the Corrosion Risks to a Double Shell Tank During the Selective Removal of a Cesium-rich Supernate Layer - 20227

The Direct Feed Low Activity Waste (DF-LAW) activity will provide the initial feed to the Waste Treatment and Immobilization Plant at the Hanford Tank Farm facility. Feed preparation for DF-LAW will be carried out in multiple double shell waste tanks (DST) utilizing the Tank Side Cesium Removal (TSCR) process. Tank AP-106 will be an integral part of the TSCR process and hence key to the DFLAW activity. Prior to being utilized for this purpose, the Cs-137 concentration levels in the AP-106 supernate must be reduced significantly. The AP-106 re-purposing process manages normal waste tank addition procedures such that waste acceptance criteria objectives for Cs-137 removal are achieved. As this process was being developed, several concerns regarding tank integrity during and after repurposing of Tank AP-106, and the associated tanks (i.e., AP-102, AN-101, and AP-108), were addressed. The potential corrosion degradation mechanisms identified based on process knowledge included: general and pitting corrosion, caustic stress corrosion cracking, concentration cell corrosion, and microbiologically induced corrosion (MIC). A systematic approach was undertaken to assess the risk of corrosion during the AP-106 re-purposing process. The assessment considered the initial and present condition of the tank, the service history of the tank, the environmental conditions created in each tank during each stage of the process, process timing, and corrosion testing and modeling that approximated tank conditions. Based on this information, recommendations were provided to minimize the exposure of the tank to corrosive conditions. These recommendations included time limits for each stage of the process, staging of corrosion inhibitor solutions should process delays occur, and waste sampling with potential chemistry adjustments at defined steps in the process. Additionally, recommendations for post-process supernate chemistry adjustments and tank inspections were put forth. The greatest risk for this process occurs during the manipulation of large quantities of Columbia River water (i.e., greater than 450,000 gallons) within AP-106, the associated tanks and the supporting infrastructure (i.e., transfer line piping). This risk becomes more pronounced if there is a process upset and the tank walls remain exposed to a corrosive condition longer than intended. As long as the river water is blended into the supernate in a timely fashion, and the waste is adjusted appropriately with corrosion inhibitors, corrosion may be minimized and the final condition of the supernate in AP-106 and the associated tanks is expected to be in compliance with corrosion chemistry requirements. Contingency plans will be made to minimize the impact of a process upset. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2007

Implementing an Applied Science Program

The work implied in the NASA Applied Science Program requires a delicate balancing act for the those doing it. At the implementation level there are multiple tensions intrinsic to the program. For example each application of an existing product to a decision support process requires deep knowledge about the data and deep knowledge about the decision making process. It is highly probable no one person has this range of knowledge. Otherwise the decision making process would already be using the data. Therefore, a team is required. But building a team usually requires time, especially across agencies. Yet the program mandates efforts of relatively short duration. Further, those who know the data are scientists, which makes them essential to the program. But scientists are evaluated on their publication record. Anything which diverts a scientist from the research for his next publication is an anathema to him and potential death to their career. Trying to get another agency to use NASA data does not strike most scientists as material inherently suitable for publication. Also, NASA wishes to rapidly implement often substantial changes to another agency's process. For many reasons, such as budget and program constraints, speed is important. But the owner of a decision making process is tightly constrained, usually by law, regulation, organization and custom. Changes when made are slow, cautious, even hesitant, and always done according a process specific to the situation. To manage this work MSFC must balance these and other tensions. Some things we have relatively little control over, such as budget. These we try to handle by structural techniques. For example by insisting all of our people work on multiple projects simultaneously we inherently have diversification of funding for all of our people. In many cases we explicitly use some elements of tension to be productive. For example the need for the scientists to constantly publish is motivation to keep tasks short and to the immediate point. As another example, the DSS agency and it personnel are treated as integral to the effort at every step. Thus, their ways of doing things, their problems and their assets become part of the solution. MSFC also manages the Applied Science work' within a strategic framework.' First, the scientists are necessarily at the core of all of this work as well as all of the work within the larger' organization. We therefore strive to keep a roughly 50-50 balance between the work done on tasks funded directly by the research side of NASA's Science Mission Directorate and funding from the applied side. Done at both the organization level and the level of the individual, this keeps the scientist both happy and productive over the long term. We also try diligently to remove as much burden from the scientist as practical by employing people such as the authors and others to do tasks not requiring; scientific knowledge. We also have designed our effort to take full use of external partners. We actively seek and support, including fund, people from multiple organizations to join us as committed collaborators. In this we use today's money and today's problem to help us diversify and strengthen for tomorrow. MSFC also considers the Applied Science work holistically. Each element is viewed as a step in a larger process. At a management level we can chose to emphasize or encourage certain areas which service long term goals. Thus, if we think work in a particular area should be developed, we can start with the smaller, less costly elements and grow.

Rickman, Doug↗
NASA NTRS2012

Absolute Position of Targets Measured Through a Chamber Window Using Lidar Metrology Systems

Lidar is a useful tool for taking metrology measurements without the need for physical contact with the parts under test. Lidar instruments are aimed at a target using azimuth and elevation stages, then focus a beam of coherent, frequency modulated laser energy onto the target, such as the surface of a mechanical structure. Energy from the reflected beam is mixed with an optical reference signal that travels in a fiber path internal to the instrument, and the range to the target is calculated based on the difference in the frequency of the returned and reference signals. In cases when the parts are in extreme environments, additional steps need to be taken to separate the operator and lidar from that environment. A model has been developed that accurately reduces the lidar data to an absolute position and accounts for the three media in the testbed air, fused silica, and vacuum but the approach can be adapted for any environment or material. The accuracy of laser metrology measurements depends upon knowing the parameters of the media through which the measurement beam travels. Under normal conditions, this means knowledge of the temperature, pressure, and humidity of the air in the measurement volume. In the past, chamber windows have been used to separate the measuring device from the extreme environment within the chamber and still permit optical measurement, but, so far, only relative changes have been diagnosed. The ability to make accurate measurements through a window presents a challenge as there are a number of factors to consider. In the case of the lidar, the window will increase the time-of-flight of the laser beam causing a ranging error, and refract the direction of the beam causing angular positioning errors. In addition, differences in pressure, temperature, and humidity on each side of the window will cause slight atmospheric index changes and induce deformation and a refractive index gradient within the window. Also, since the window is a dispersive media, the effect of both phase and group indices have to be considered. Taking all these factors into account, a method was developed to measure targets through multiple regions of different materials and produce results that are absolute measurements of target position in three-dimensional space, rather than simply relative position. The environment in which the lidar measurements are taken must be broken down into separate regions of interest and each region solved for separately. In this case, there were three regions of interest: air, fused silica, and vacuum. The angular position of the target inside the chamber is solved using only phase index and phase velocity, while the ranging effects due to travel from air to glass to vacuum/air are solved with group index and group velocity. When all parameters are solved simultaneously, an absolute knowledge of the position of each target within an environmental chamber can be derived. Novel features of this innovation include measuring absolute position of targets through multiple dispersive and non-dispersive media, deconstruction of lidar raw data from a commercial off-the-shelf unit into reworkable parameters, and use of group velocities to reduce range data. Measurement of structures within a vacuum chamber or other harsh environment, such as a furnace, may now be measured as easily as if they were in an ambient laboratory. This analysis permits transformation of the raw data into absolute spatial units (e.g., mm). This technique has also been extended to laser tracker, theodolite, and cathetometer measurements through refractive media.

Kubalak, David↗
DOE OSTI2024

Tracking the reaction networks of acetaldehyde oxide and glyoxal oxide Criegee intermediates in the ozone-assisted oxidation reaction of crotonaldehyde

The reaction of unsaturated compounds with ozone (O 3 ) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH 3 CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O 2 , at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O 3 product, C 4 H 6 O 4 , is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H 2 O 2 ), methyl hydroperoxide (CH 3 OOH), and ethyl hydroperoxide (C 2 H 5 OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

Bird Migration Under Climate Change - A Mechanistic Approach Using Remote Sensing

The broad-scale reductions and shifts that may be expected under climate change in the availability and quality of stopover habitat for long-distance migrants is an area of increasing concern for conservation biologists. Researchers generally have taken two broad approaches to the modeling of migration behaviour to understand the impact of these changes on migratory bird populations. These include models based on causal processes and their response to environmental stimulation, "mechanistic models", or models that primarily are based on observed animal distribution patterns and the correlation of these patterns with environmental variables, i.e. "data driven" models. Investigators have applied the latter technique to forecast changes in migration patterns with changes in the environment, for example, as might be expected under climate change, by forecasting how the underlying environmental data layers upon which the relationships are built will change over time. The learned geostatstical correlations are then applied to the modified data layers.. However, this is problematic. Even if the projections of how the underlying data layers will change are correct, it is not evident that the statistical relationships will remain the same, i.e. that the animal organism may not adapt its' behaviour to the changing conditions. Mechanistic models that explicitly take into account the physical, biological, and behaviour responses of an organism as well as the underlying changes in the landscape offer an alternative to address these shortcomings. The availability of satellite remote sensing observations at multiple spatial and temporal scales, coupled with advances in climate modeling and information technologies enable the application of the mechanistic models to predict how continental bird migration patterns may change in response to environmental change. In earlier work, we simulated the impact of effects of wetland loss and inter-annual variability on the fitness of migratory shorebirds in the central fly ways of North America. We demonstrated the phenotypic plasticity of a migratory population of Pectoral sandpipers consisting of an ensemble of 10,000 individual birds in response to changes in stopover locations using an individual based migration model driven by remotely sensed land surface data, climate data and biological field data. With the advent of new computing capabilities enabled hy recent GPU-GP computing paradigms and commodity hardware, it now is possible to simulate both larger ensemble populations and to incorporate more realistic mechanistic factors into migration models. Here, we take our first steps use these tools to study the impact of long-term drought variability on shorebird survival.

Smith, James A.↗
DOE OSTI2021

Constitutive Model Development for Aging Polymer Encapsulants (ASC P&EM FY2021 L2 Milestone 7836)

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Numerical Modeling of the Chilldown of Cryogenic Transfer Lines Using a Sinda/GFSSP Integrated Solver

An important first step in cryogenic propellant loading is the chilldown of transfer lines. During the chilldown of the transfer line, the flow is two-phase and unsteady, with solid to fluid heat transfer and therefore a coupled thermo-fluid analysis is necessary to model the system. This paper describes a numerical model of pipe chilldown that utilizes the Sinda/GFSSP Conjugate Integrator (SGCI). SGCI is a new analysis tool developed at NASA's Marshall Space Flight Center (MSFC). SGCI facilitates the solution of thermofluid problems in interconnected solid-fluid systems. The solid component of the system is modeled in MSC Patran and translated into an MSC Sinda thermal network model. The fluid component is modeled in GFSSP, the Generalized Fluid System Simulation Program. GFSSP is a general network flow solver developed at NASA/MSFC. GFSSP uses a finite-volume approach to model fluid systems that can include phase change, multiple species, fluid transients, and heat transfer to simple solid networks. SGCI combines the GFSSP Fortran code with the Sinda input file and compiles the integrated model. Sinda solves for the temperatures of the solid network, while GFSSP simultaneously solves the fluid network for pressure, temperature, and flow rate. The two networks are coupled by convection heat transfer from the solid wall to the cryogenic fluid. The model presented here is based on a series of experiments conducted in 1966 by the National Bureau of Standards (NBS). A vacuum-jacketed, 200 ft copper transfer line was chilled by liquid nitrogen and liquid hydrogen. The predictions of transient temperature profiles and chilldown time of the integrated Sinda/GFSSP model will be compared to the experimental measurements.

LeClair, Andre↗
NASA NTRS2010

Real-Time Feature Tracking Using Homography

This software finds feature point correspondences in sequences of images. It is designed for feature matching in aerial imagery. Feature matching is a fundamental step in a number of important image processing operations: calibrating the cameras in a camera array, stabilizing images in aerial movies, geo-registration of images, and generating high-fidelity surface maps from aerial movies. The method uses a Shi-Tomasi corner detector and normalized cross-correlation. This process is likely to result in the production of some mismatches. The feature set is cleaned up using the assumption that there is a large planar patch visible in both images. At high altitude, this assumption is often reasonable. A mathematical transformation, called an homography, is developed that allows us to predict the position in image 2 of any point on the plane in image 1. Any feature pair that is inconsistent with the homography is thrown out. The output of the process is a set of feature pairs, and the homography. The algorithms in this innovation are well known, but the new implementation improves the process in several ways. It runs in real-time at 2 Hz on 64-megapixel imagery. The new Shi-Tomasi corner detector tries to produce the requested number of features by automatically adjusting the minimum distance between found features. The homography-finding code now uses an implementation of the RANSAC algorithm that adjusts the number of iterations automatically to achieve a pre-set probability of missing a set of inliers. The new interface allows the caller to pass in a set of predetermined points in one of the images. This allows the ability to track the same set of points through multiple frames.

Clouse, Daniel S.↗
DOE OSTI2020

Leveraging 3D CAD and Virtual Reality in Design of the Calcine Disposition Project - 20512

Fluor Idaho's Calcine Disposition Project (CDP) is leveraging the power of three-dimensional (3D) visualization, light detection and ranging (Lidar), and virtual reality to improve designs and worker safety. With the use of computer-aided design (CAD) modeling software, the CDP has created 3D models of the facility and systems to aid in the design, development, and operation of its retrieval and transfer system. This by itself is nothing new and has been done in industry for years. The CDP project however is taking this a few steps farther with the use of Lidar and virtual reality software. The CRP is tasked with removing a radioactive granular material from stainless steel bins located in concrete vaults that were constructed in the late 1950's and early 1960's. While construction drawings are available, it is not certain that the drawings are as-built or how accurate they truly are. Part of the project requires precise placement of equipment on to bins around multiple pipe, electrical lines, and bin stiffeners. Relying solely on the existing drawings and models created from these drawings is highly risky and prone to failure. The first step to determine the precise location and orientation of all obstructions in the vault the project will use a Lidar system to scan the vault structure and precisely locate all vault and bin components. Once this information is obtained, the data will be integrated into the CAD model. The model will then be verified, and a true as-built model developed. Precise component location, size, and orientation will be used for final design and placement of retrieval components. The second step will be to take the verified 3D model and scan data and import it into virtual reality software creating a virtual world. While some work on the project requires all work to be completed remotely due to high radiation fields, some areas have low enough radiation levels that personnel can enter and perform work. With the use of a virtual reality headset, operators and technicians will be able to enter our virtual world and become familiar with the surroundings and perform work prior to ever entering the radiation environment. This training is expected to pay dividends by improving worker efficiency, reducing errors, and improving confidence that the work can be performed as expected. It will play a key aspect in as low as reasonably achievable (ALARA) principles. An additional aspect in using the Lidar scans and virtual environment, is that it will give engineers and managers an opportunity to revisit the facility, especially radiation areas, at any time to obtain design information and measurements that may be costly or impossible to obtain once the facility transitions into radiation operations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

HT-SIP: a semi-automated stable isotope probing pipeline identifies cross-kingdom interactions in the hyphosphere of arbuscular mycorrhizal fungi

Abstract Background Linking the identity of wild microbes with their ecophysiological traits and environmental functions is a key ambition for microbial ecologists. Of many techniques that strive for this goal, Stable-isotope probing—SIP—remains among the most comprehensive for studying whole microbial communities in situ. In DNA-SIP, actively growing microorganisms that take up an isotopically heavy substrate build heavier DNA, which can be partitioned by density into multiple fractions and sequenced. However, SIP is relatively low throughput and requires significant hands-on labor. We designed and tested a semi-automated, high-throughput SIP (HT-SIP) pipeline to support well-replicated, temporally resolved amplicon and metagenomics experiments. We applied this pipeline to a soil microhabitat with significant ecological importance—the hyphosphere zone surrounding arbuscular mycorrhizal fungal (AMF) hyphae. AMF form symbiotic relationships with most plant species and play key roles in terrestrial nutrient and carbon cycling. Results Our HT-SIP pipeline for fractionation, cleanup, and nucleic acid quantification of density gradients requires one-sixth of the hands-on labor compared to manual SIP and allows 16 samples to be processed simultaneously. Automated density fractionation increased the reproducibility of SIP gradients compared to manual fractionation, and we show adding a non-ionic detergent to the gradient buffer improved SIP DNA recovery. We applied HT-SIP to 13 C-AMF hyphosphere DNA from a 13 CO 2 plant labeling study and created metagenome-assembled genomes (MAGs) using high-resolution SIP metagenomics (14 metagenomes per gradient). SIP confirmed the AMF Rhizophagus intraradices and associated MAGs were highly enriched (10–33 atom% 13 C), even though the soils’ overall enrichment was low (1.8 atom% 13 C). We assembled 212 13 C-hyphosphere MAGs; the hyphosphere taxa that assimilated the most AMF-derived 13 C were from the phyla Myxococcota, Fibrobacterota, Verrucomicrobiota, and the ammonia-oxidizing archaeon genus Nitrososphaera . Conclusions Our semi-automated HT-SIP approach decreases operator time and improves reproducibility by targeting the most labor-intensive steps of SIP—fraction collection and cleanup. We illustrate this approach in a unique and understudied soil microhabitat—generating MAGs of actively growing microbes living in the AMF hyphosphere (without plant roots). The MAGs’ phylogenetic composition and gene content suggest predation, decomposition, and ammonia oxidation may be key processes in hyphosphere nutrient cycling.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2016

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗
DOE OSTI2021

Modeling Framework for Bulk Electric Grid Impacts from HEMP E1 and E3 Effects (Tasks 3.1 Final Report)

This report presents a framework to evaluate the impact of a high-altitude electromagnetic pulse (HEMP) event on a bulk electric power grid. This report limits itself to modeling the impact of EMP E1 and E3 components. The co-simulation of E1 and E3 is presented in detail, and the focus of the paper is on the framework rather than actual results. This approach is highly conservative as E1 and E3 are not maximized with the same event characteristics and may only slightly overlap. The actual results shown in this report are based on a synthetic grid with synthetic data and a limited exemplary EMP model. The framework presented can be leveraged and used to analyze the impact of other threat scenarios, both manmade and natural disasters. This report d escribes a Monte-Carlo based methodology to probabilistically quantify the transient response of the power grid to a HEMP event. The approach uses multiple fundamental steps to characterize the system response to HEMP events, focused on the E1 and E3 components of the event. 1) Obtain component failure data related to HEMP events testing of components and creating component failure models. Use the component failure model to create component failure conditional probability density function (PDF) that is a function of the HEMP induced terminal voltage. 2) Model HEMP scenarios and calculate the E1 coupled voltage profiles seen by all system components. Model the same HEMP scenarios and calculate the transformer reactive power consumption profiles due to E3. 3) Sample each component failure PDF to determine which grid components will fail, due to the E1 voltage spike, for each scenario. 4) Perform dynamic simulations that incorporate the predicted component failures from E1 and reactive power consumption at each transformer affected by E3. These simulations allow for secondary transients to affect the relays/protection remaining in service which can lead to cascading outages. 5) Identify the locations and amount of load lost for each scenario through grid dynamic simulation. This can be an indication of the immediate grid impacts from a HEMP event. In addition, perform more detailed analysis to determine critical nodes and system trends. 6) To help realize the longer-term impacts, a security constrained alternating current optimal power flow (ACOPF) is run to maximize critical load served. This report describes a modeling framework to assess the systemic grid impacts due to a HEMP event. This stochastic simulation framework generates a large amount of data for each Monte Carlo replication, including HEMP location and characteristics, relay and component failures, E3 GIC profiles, cascading dynamics including voltage and frequency over time, and final system state. This data can then be analyzed to identify trends, e.g., unique system behavior modes or critical components whose failure is more likely to cause serious systemic effects. The proposed analysis process is demonstrated on a representative system. In order to draw realistic conclusions of the impact of a HEMP event on the grid, a significant amount of work remains with respect to modeling the impact on various grid components.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Striped Data Analysis Framework

A columnar data representation is known to be an efficient way for data storage, specifically in cases when the analysis is often done based only on a small fragment of the available data structures. A data representation like Apache Parquet is a step forward from a columnar representation, which splits data horizontally to allow for easy parallelization of data analysis. Based on the general idea of columnar data storage, working on the [LDRD Project], we have developed a striped data representation, which, we believe, is better suited to the needs of High Energy Physics data analysis. A traditional columnar approach allows for efficient data analysis of complex structures. While keeping all the benefits of columnar data representations, the striped mechanism goes further by enabling easy parallelization of computations without requiring special hardware. We will present an implementation and some performance characteristics of such a data representation mechanism using a distributed no-SQL database or a local file system, unified under the same API and data representation model. The representation is efficient and at the same time simple so that it allows for a common data model and APIs for wide range of underlying storage mechanisms such as distributed no-SQL databases and local file systems. Striped storage adopts Numpy arrays as its basic data representation format, which makes it easy and efficient to use in Python applications. The Striped Data Server is a web service, which allows to hide the server implementation details from the end user, easily exposes data to WAN users, and allows to utilize well known and developed data caching solutions to further increase data access efficiency. We are considering the Striped Data Server as the core of an enterprise scale data analysis platform for High Energy Physics and similar areas of data processing. We have been testing this architecture with a 2TB dataset from a CMS dark matter search and plan to expand it to multiple 100 TB or even PB scale. We will present the striped format, Striped Data Server architecture and performance test results.

Gutsche, Oliver↗
NASA NTRS2014

Pan-Continental Droughts in North America over the Last Millennium

Regional droughts are common in North America, but pan-continental droughts extending across multiple regions, including the 2012 event, are rare relative to single-region events. Here, the tree-ring-derived North American Drought Atlas is used to investigate drought variability in four regions over the last millennium, focusing on pan-continental droughts. During the Medieval Climate Anomaly (MCA), the central plains (CP), Southwest (SW), and Southeast (SE) regions experienced drier conditions and increased occurrence of droughts and the Northwest (NW) experienced several extended pluvials. Enhanced MCA aridity in the SW and CP manifested as multidecadal megadroughts. Notably, megadroughts in these regions differed in their timing and persistence, suggesting that they represent regional events influenced by local dynamics rather than a unified, continental-scale phenomena. There is no trend in pan-continental drought occurrence, defined as synchronous droughts in three or more regions. SW, CP, and SE (SW+CP+SE) droughts are the most common, occurring in 12 percent of all years and peaking in prevalence during the twelfth and thirteenth centuries; patterns involving three other regions occur in about 8 percent of years. Positive values of the Southern Oscillation index (La Nina conditions) are linked to SW, CP, and SE (SW+CP+SE) droughts and SW, CP, and NW (SW+CP+NW) droughts, whereas CP, NW, and SE (CP+NW+SE) droughts are associated with positive values of the Pacific decadal oscillation and Atlantic multidecadal oscillation. While relatively rare, pan-continental droughts are present in the paleo record and are linked to defined modes of climate variability, implying the potential for seasonal predictability. Assuming stable drought teleconnections, these events will remain an important feature of future North American hydroclimate, possibly increasing in their severity in step with other expected hydroclimate responses to increased greenhouse gas forcing.

anomalies↗
NASA NTRS2016

Manufacturing Challenges and Benefits when Scaling the HIAD Stacked-Torus Aeroshell to a 15m Class System

Over a decade of work has been conducted in the development of NASA's Hypersonic Inflatable Aerodynamic Decelerator (HIAD) deployable aeroshell technology. This effort has included multiple ground test campaigns and flight tests culminating in the HIAD project's second generation (Gen-2) aeroshell system. The HIAD project team has developed, fabricated, and tested stacked-torus inflatable structures (IS) with flexible thermal protection systems (F-TPS) ranging in diameters from 3-6m, with cone angles of 60 and 70 deg. To meet NASA and commercial near term objectives, the HIAD team must scale the current technology up to 12-15m in diameter. The HIAD project's experience in scaling the technology has reached a critical juncture. Growing from a 6m to a 15m class system will introduce many new structural and logistical challenges to an already complicated manufacturing process. Although the general architecture and key aspects of the HIAD design scale well to larger vehicles, details of the technology will need to be reevaluated and possibly redesigned for use in a 15m-class HIAD system. These include: layout and size of the structural webbing that transfers load throughout the IS, inflatable gas barrier design, torus diameter and braid construction, internal pressure and inflation line routing, adhesives used for coating and bonding, and F-TPS gore design and seam fabrication. The logistics of fabricating and testing the IS and the F-TPS also become more challenging with increased scale. Compared to the 6m aeroshell (the largest HIAD built to date), a 12m aeroshell has four times the cross-sectional area, and a 15m one has over six times the area. This means that fabrication and test procedures will need to be reexamined to account for the sheer size and weight of the aeroshell components. This will affect a variety of steps in the manufacturing process, such as: stacking the tori during assembly, stitching the structural webbing, initial inflation of tori, and stitching of F-TPS gores. Additionally, new approaches and hardware will be required for handling and ground testing of both individual tori and the fully assembled HIADs. There are also noteworthy benefits of scaling up the HIAD aeroshell to 15m-class system. Two complications in working with handmade textiles structures are the non-linearity of the materials and the role of human accuracy during fabrication. Larger, more capable, HIAD structures should see much larger operational loads, potentially bringing the structural response of the materials out of the non-linear regime and into the preferred linear response range. Also, making the reasonable assumption that the magnitude of fabrication accuracy remains constant as the structures grow, the relative effect of fabrication errors should decrease as a percentage of the textile component size. Combined, these two effects improve the predictive capability and the uniformity of the structural response for a 12-15m class HIAD. In this paper, the challenges and associated mitigation plans related to scaling up the HIAD stacked-torus aeroshell to a 15m class system will be discussed. In addition, the benefits of enlarging the structure will be further explored.

Decellerators↗
NASA NTRS2016

Manufacturing Challenges and Benefits when Scaling the HIAD Stacked-Torus Aeroshell to a 15m-Class System

Over a decade of work has been conducted in the development of NASAs Hypersonic Inflatable Aerodynamic Decelerator (HIAD) deployable aeroshell technology. This effort has included multiple ground test campaigns and flight tests culminating in the HIAD projects second generation (Gen-2) aeroshell system. The HIAD project team has developed, fabricated, and tested stacked-torus inflatable structures (IS) with flexible thermal protection systems (F-TPS) ranging in diameters from 3-6m, with cone angles of 60 and 70 deg. To meet NASA and commercial near term objectives, the HIAD team must scale the current technology up to 12-15m in diameter. Therefore, the HIAD projects experience in scaling the technology has reached a critical juncture. Growing from a 6m to a 15m-class system will introduce many new structural and logistical challenges to an already complicated manufacturing process.Although the general architecture and key aspects of the HIAD design scale well to larger vehicles, details of the technology will need to be reevaluated and possibly redesigned for use in a 15m-class HIAD system. These include: layout and size of the structural webbing that transfers load throughout the IS, inflatable gas barrier design, torus diameter and braid construction, internal pressure and inflation line routing, adhesives used for coating and bonding, and F-TPS gore design and seam fabrication. The logistics of fabricating and testing the IS and the F-TPS also become more challenging with increased scale. Compared to the 6m aeroshell (the largest HIAD built to date), a 12m aeroshell has four times the cross-sectional area, and a 15m one has over six times the area. This means that fabrication and test procedures will need to be reexamined to ac-count for the sheer size and weight of the aeroshell components. This will affect a variety of steps in the manufacturing process, such as: stacking the tori during assembly, stitching the structural webbing, initial inflation of tori, and stitching of F-TPS gores. Additionally, new approaches and hardware will be required for handling and ground testing of both individual tori and the fully assembled HIADs.There are also noteworthy benefits of scaling up the HIAD aeroshell to a 15m-class system. Two complications in working with handmade textile structures are the non-linearity of the material components and the role of human accuracy during fabrication. Larger, more capable, HIAD structures should see much larger operational loads, potentially bringing the structural response of the material components out of the non-linear regime and into the preferred linear response range. Also, making the reasonable assumption that the magnitude of fabrication accuracy remains constant as the structures grow, the relative effect of fabrication errors should decrease as a percentage of the textile component size. Combined, these two effects improve the predictive capability and the uniformity of the structural response for a 12-15m HIAD.In this presentation, a handful of the challenges and associated mitigation plans will be discussed, as well as an update on current 12m aeroshell manufacturing and testing that is addressing these challenges

Infatables↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗