Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer model”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Progressive Damage Analysis of Laminated Composite (PDALC)-A Computational Model Implemented in the NASA COMET Finite Element Code

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete list of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occur during the load history. Residual strength predictions made with this information compared favorably with experimental measurements.

Lo, David C.↗
NASA NTRS1998

Progressive Damage Analysis of Laminated Composite (PDALC) (A Computational Model Implemented in the NASA COMET Finite Element Code)

A method for analysis of progressive failure in the Computational Structural Mechanics Testbed is presented in this report. The relationship employed in this analysis describes the matrix crack damage and fiber fracture via kinematics-based volume-averaged damage variables. Damage accumulation during monotonic and cyclic loads is predicted by damage evolution laws for tensile load conditions. The implementation of this damage model required the development of two testbed processors. While this report concentrates on the theory and usage of these processors, a complete listing of all testbed processors and inputs that are required for this analysis are included. Sample calculations for laminates subjected to monotonic and cyclic loads were performed to illustrate the damage accumulation, stress redistribution, and changes to the global response that occurs during the loading history. Residual strength predictions made with this information compared favorably with experimental measurements.

Coats, Timothy W.↗
NASA NTRS1994

Computer Modeling Of Atomization

Improved mathematical models based on fundamental principles of conservation of mass, energy, and momentum developed for use in computer simulation of atomization of jets of liquid fuel in rocket engines. Models also used to study atomization in terrestrial applications; prove especially useful in designing improved industrial sprays - humidifier water sprays, chemical process sprays, and sprays of molten metal. Because present improved mathematical models based on first principles, they are minimally dependent on empirical correlations and better able to represent hot-flow conditions that prevail in rocket engines and are too severe to be accessible for detailed experimentation.

Giridharan, M.↗
DOE OSTI2023

Discovering Active Subspaces for High-Dimensional Computer Models

Dimension reduction techniques have long been an important topic in statistics, and active subspaces (AS) have received much attention this past decade in the computer experiments literature. The most common approach towards estimating the AS is to use Monte Carlo with numerical gradient evaluation. While sensible in some settings, this approach has obvious drawbacks. Recent research has demonstrated that active subspace calculations can be obtained in closed form, conditional on a Gaussian process (GP) surrogate, which can be limiting in high-dimensional settings for computational reasons. In this paper, we produce the relevant calculations for a more general case when the model of interest is a linear combination of tensor products. These general equations can be applied to the GP, recovering previous results as a special case, or applied to the models constructed by other regression techniques including multivariate adaptive regression splines (MARS). Furthermore, using a MARS surrogate has many advantages including improved scaling, better estimation of active subspaces in high dimensions and the ability to handle a large number of prior distributions in closed form. In one real-world example, we obtain the active subspace of a radiation-transport code with 240 inputs and 9,372 model runs in under half an hour.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1994

Representing spatial information in a computational model for network management

While currently available relational database management systems (RDBMS) allow inclusion of spatial information in a data model, they lack tools for presenting this information in an easily comprehensible form. Computer-aided design (CAD) software packages provide adequate functions to produce drawings, but still require manual placement of symbols and features. This project has demonstrated a bridge between the data model of an RDBMS and the graphic display of a CAD system. It is shown that the CAD system can be used to control the selection of data with spatial components from the database and then quickly plot that data on a map display. It is shown that the CAD system can be used to extract data from a drawing and then control the insertion of that data into the database. These demonstrations were successful in a test environment that incorporated many features of known working environments, suggesting that the techniques developed could be adapted for practical use.

Blaisdell, James H.↗
DOE OSTI2020

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1978

Computer model for refinery operations with emphasis on jet fuel production. Volume 2: Data and technical bases

The FORTRAN computing program predicts the flow streams and material, energy, and economic balances of a typical petroleum refinery, with particular emphasis on production of aviation turbine fuel of varying end point and hydrogen content specifications. The program has provision for shale oil and coal oil in addition to petroleum crudes. A case study feature permits dependent cases to be run for parametric or optimization studies by input of only the variables which are changed from the base case. The report has sufficient detail for the information of most readers.

Dunbar, D. N.↗
NASA NTRS1984

Computational modeling of jet induced mixing of cryogenic propellants in low-G

The SOLA-ECLIPSE Code is being developed to enable computational prediction of jet induced mixing in cryogenic propellant tanks in a low-gravity environment. Velocity fields, predicted for scale model tanks, are presented which compare favorably with the available experimental data. A full scale liquid hydrogen tank for a typical Orbit Transfer Vehicle is analyzed with the conclusion that coupling an axial mixing jet with a thermodynamic vent system appears to be a viable concept for the control of tank pressure.

Hochstein, J. I.↗
NASA NTRS1984

Computational modeling of jet induced mixing of cryogenic propellants in low-G

The SOLA-ECLIPSE Code is being developed to enable computational prediction of jet induced mixing in cryogenic propellant tanks in a low-gravity environment. Velocity fields, predicted for scale model tanks, are presented which compare favorably with the available experimental data. A full scale liquid hydrogen tank for a typical Orbit Transfer Vehicle is analyzed with the conclusion that coupling an axial mixing jet with a thermodynamic vent system appears tobe a viable concept for the control of tank pressure. Previously announced in STAR as N84-25000

Hochstein, J. I.↗
DOE OSTI2021

A comparison of Gaussian processes and neural networks for computer model emulation and calibration

The Department of Energy relies on complex physics simulations for prediction in domains like cosmology, nuclear theory, and materials science. These simulations are often extremely computationally intensive, with some requiring days or weeks for a single simulation. In order to assure their accuracy, these models are calibrated against observational data in order to estimate inputs and systematic biases. Because of their great computational complexity, this process typically requires the construction of an emulator, a fast approximation to the simulation. In this paper, two emulator approaches are compared: Gaussian process regression and neural networks. Their emulation accuracy and calibration performance on three real problems of Department of Energy interest is considered. On these problems, the Gaussian process emulator tends to be more accurate with narrower, but still well-calibrated uncertainty estimates. The neural network emulator is accurate, but tends to have large uncertainty on its predictions. Finally, as a result, calibration with the Gaussian process emulator produces more constrained posteriors that still perform well in prediction.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Computer model for refinery operations with emphasis on jet fuel production. Volume 3: Detailed systems and programming documentation

The FORTRAN computing program predicts flow streams and material, energy, and economic balances of a typical petroleum refinery, with particular emphasis on production of aviation turbine fuels of varying end point and hydrogen content specifications. The program has a provision for shale oil and coal oil in addition to petroleum crudes. A case study feature permits dependent cases to be run for parametric or optimization studies by input of only the variables which are changed from the base case.

Dunbar, D. N.↗
NASA NTRS1987

A new computational model for the prediction of mass loading phenomena for solar wind interactions with cometary and planetary ionospheres

The modified gasdynamic convected-magnetic-field MHD model developed by Spreiter and Stahara (1980) to simulate the supersonic flow of the solar wind past planetary magnetoionospheres is extended to account for cases (such as Venus and comets) in which significant numbers of neutral atmospheric atoms become ionized in the surrounding flow and add to its momentum, energy, and mass. The mathematical model and the solution procedures for the nose and tail regions are explained; typical computational grids are shown; and numerical results for a comet and for the Venus ionosheath and bow shock are presented graphically. It is found that the bow shock weakens and moves further upstream of the obstacle as mass loading is increased and the flow upstream of the bow shock becomes more compressed.

Stahara, Stephen S.↗
DOE OSTI2021

Accelerating computational modeling and design of high-entropy alloys

High-entropy alloys, with N elements and compositions {$c_{ν = 1,N}$} in competing crystal structures, have large design spaces for unique chemical and mechanical properties. In this work, to enable computational design, we use a metaheuristic hybrid Cuckoo search (CS) to construct alloy configurational models on the fly that have targeted atomic site and pair probabilities on arbitrary crystal lattices, given by supercell random approximates (SCRAPs) with S sites. Our Hybrid CS permits efficient global solutions for large, discrete combinatorial optimization that scale linearly in a number of parallel processors, and linearly in sites S for SCRAPs. For example, a four-element, 128-site SCRAP is found in seconds—a more than 13,000-fold reduction over current strategies. Our method thus enables computational alloy design that is currently impractical. We qualify the models and showcase application to real alloys with targeted atomic short-range order. Being problem-agnostic, our Hybrid CS offers potential applications in diverse fields.

36 MATERIALS SCIENCE↗
NASA NTRS1990

A strategy to study regional hydrology and terrestrial ecosystem processes using satellite remote sensing, ground-based data and computer modeling

A strategy is presented for integrating scanning multichannel microwave radiometer data from the Nimbus-7 satellite with meteorological station records and computer simulations of land surface hydrology, terrestrial nutrient cycling, and trace gas emission. Analysis of the observations together with radiative transfer analysis shows that in the tropics the temporal and spatial variations of the polarization difference are determined primarily by the structure and phenology of vegetation and seasonal inundations of major rivers and wetlands. It is concluded that the proposed surface hydrology model, along with climatological records, and, potentially, 37-GHz data for phenology, will provide inputs to a terrestrial ecosystem model that predicts regional net primary production and CO2 gas exchange.

Vorosmarty, C.↗
NASA NTRS1988

Computer modeling of large asteroid impacts into continental and oceanic sites: Atmospheric, cratering, and ejecta dynamics

Numerous impact cratering events have occurred on the Earth during the last several billion years that have seriously affected our planet and its atmosphere. The largest cratering events, which were caused by asteroids and comets with kinetic energies equivalent to tens of millions of megatons of TNT, have distributed substantial quantities of terrestrial and extraterrestrial material over much or all of the Earth. In order to study a large-scale impact event in detail, computer simulations were completed that model the passage of a 10 km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics associated with impact of the asteroid into two different targets, i.e., an oceanic site and a continental site. The calcuations were designed to broadly represent giant impact events that have occurred on the Earth since its formation and specifically represent an impact cratering event proposed to have occurred at the end of Cretaceous time. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock that expanded radially outward. Behind the shock front was a region of highly shock compressed and intensely heated air. Behind the asteroid, rapid expansion of this shocked air created a large region of very low density that also expanded away from the impact area. Calculations of the cratering events in both the continental and oceanic targets were carried to 120 s. Despite geologic differences, impacts in both targets developed comparable dynamic flow fields, and by approx. 29 s similar-sized transient craters approx. 39 km deep and approx. 62 km across had formed. For all practical purposes, the atmosphere was nearly completely removed from the impact area for tens of seconds, i.e., air pressures were less than fractions of a bar out to ranges of over 50 km. Consequently, much of the asteroid and target materials were ejected upward into a near vacuum. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to add substantial solid and vaporized material to the atmosphere, but these conditions were not studied.

Roddy, D. J.↗
NASA NTRS1998

Computational Modeling of Semiconductor Dynamics at Femtosecond Time Scales

The main objective of the Joint-Research Interchange NCC2-5149 was to develop computer codes for accurate simulation of femtosecond pulse propagation in semiconductor lasers and semiconductor amplifiers [I]. The code should take into account all relevant processes such as the interband and intraband carrier relaxation mechanisms and the many-body effects arising from the Coulomb interaction among charge carriers [2]. This objective was fully accomplished. We made use of a previously developed algorithm developed at NASA Ames [3]-[5]. The new algorithm was tested on several problems of practical importance. One such problem was related to the amplification of femtosecond optical pulses in semiconductors. These results were presented in several international conferences over a period of three years. With the help of a postdoctoral fellow, we also investigated the origin of instabilities that can lead to the formation of femtosecond pulses in different kinds of lasers. We analyzed the occurrence of absolute instabilities in lasers that contain a dispersive host material with third-order nonlinearities. Starting from the Maxwell-Bloch equations, we derived general multimode equations to distinguish between convective and absolute instabilities. We find that both self-phase modulation and intensity-dependent absorption can dramatically affect the absolute stability of such lasers. In particular, the self-pulsing threshold (the so-called second laser threshold) can occur at few times the first laser threshold even in good-cavity lasers for which no self-pulsing occurs in the absence of intensity-dependent absorption. These results were presented in an international conference and published in the form of two papers.

Agrawal, Govind P.↗