Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational modeling and simulations”

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 883 records · Page 49

An equation for the evolution of solar and stellar flare loops

An ordinary differential equation describing the evolution of a coronal loop subjected to a spatially uniform but time-varying heating rate is discussed. It is assumed that the duration of heating is long compared to the sound transit time through the loop, which is assumed to have uniform cross section area. The form of the equation changes as the loop evolves through three states: 'strong evaporation', 'scaling law behavior', and 'strong condensation'. Solutions to the equation may be used to compute the time dependence of the average coronal temperature and emission measure for an assumed temporal variation of the flare heating rate. The results computed from the model agree reasonably well with recent published numerical simulations and may be obtained with far less computational effort. The model is then used to study the May 21, 1980, solar flare observed by SMM and the giant April 12, 1985, flare observed on the star AD Leo.

Fisher, George H.↗

Evaluation of thermochemical models for particle and continuum simulations of hypersonic flow

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered that are dominated in turn by vibrational relaxation, weak dissociation, strong dissociation and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Boyd, Iain D.↗

Simulation of Local Tie Accuracy on VLBI Antennas

We introduce a new mathematical model to compute the centering parameters of a VLBI antenna. These include the coordinates of the reference point, axis offset, orientation, and non-perpendicularity of the axes. Using the model we simulated how precisely parameters can be computed in different cases. Based on the simulation we can give some recommendations and practices to control the accuracy and reliability of the local ties at the VLBI sites.

Kallio, Ulla↗

Three-Dimensional Computational Model for Flow in an Over-Expanded Nozzle With Porous Surfaces

A three-Dimensional computational model is used to simulate flow in a non-axisymmetric, convergent-divergent nozzle incorporating porous cavities for shock-boundary layer interaction control. The nozzle has an expansion ratio (exit area/throat area) of 1.797 and a design nozzle pressure ratio of 8.78. Flow fields for the baseline nozzle (no porosity) and for the nozzle with porous surfaces of 10% openness are computed for Nozzle Pressure Ratio (NPR) varying from 1.29 to 9.54. The three dimensional computational results indicate that baseline (no porosity) nozzle performance is dominated by unstable, shock-induced, boundary-layer separation at over-expanded conditions. For NPR less than or equal to 1.8, the separation is three dimensional, somewhat unsteady, and confined to a bubble (with partial reattachment over the nozzle flap). For NPR greater than or equal to 2.0, separation is steady and fully detached, and becomes more two dimensional as NPR increased. Numerical simulation of porous configurations indicates that a porous patch is capable of controlling off design separation in the nozzle by either alleviating separation or by encouraging stable separation of the exhaust flow. In the present paper, computational simulation results, wall centerline pressure, mach contours, and thrust efficiency ratio are presented, discussed and compared with experimental data. Results indicate that comparisons are in good agreement with experimental data. The three-dimensional simulation improves the comparisons for over-expanded flow conditions as compared with two-dimensional assumptions.

Abdol-Hamid, K. S.↗

A Summary of Structural Test and Simulation Results for Hawker 4000 Composite Fuselage Hardware

NASA Langley Research Center (LaRC) was able to obtain a partially fabricated Hawker Beechcraft Model 4000 (Hawker 4000) aft aircraft section to use for study. The Hawker 4000 is a medium sized business jet (14 Code of Federal Regulations Part 25 – Transport Category Aircraft) with seating for 8 to 14 people and fabricated out of carbon composite. The Hawker 4000 was designed in the early 2000s to compete with Gulfstream and Bombardier business jet units. The Hawker Beechcraft business unit declared bankruptcy in 2012, and abruptly halted production of the 4000 series aircraft. As a result, only 73 full production aircraft were built. This all-composite aft fuselage section was an ideal candidate to utilize for investigations of composite structures undergoing dynamic loading in order to satisfy NASA Revolutionary Vertical Lift Technology (RVLT) project objectives. Thus, a team at NASA LaRC embarked on a test and simulation program in late 2021 in order to generate test data of composite structures undergoing dynamic loading conditions. As a part of these investigations, methods including reverse engineering the geometry, conducting as-fabricated materials testing and using non-destructive evaluation (NDE) techniques such as ultrasonic through transmission were utilized in order to develop digital twin finite element models for use in computational simulations.

dynamic composite testing↗

User's instructions for the cardiovascular Walters model

The model is a combined, steady-state cardiovascular and thermal model. It was originally developed for interactive use, but was converted to batch mode simulation for the Sigma 3 computer. The model has the purpose to compute steady-state circulatory and thermal variables in response to exercise work loads and environmental factors. During a computer simulation run, several selected variables are printed at each time step. End conditions are also printed at the completion of the run.

Croston, R. C.↗

Assessment of Fast Reactor Hot Channel Factor Calculation Capability in Griffin and NekRS

The DOE-NE Advanced Modeling and Simulation (NEAMS) program develops, demonstrates, and deploys predictive computational modeling tools for nuclear reactor design and analysis. Under the Fast Reactor Applications activity within the Application Drivers Technical Area of NEAMS, the ability to perform high fidelity “hot channel factor” (HCF) coupled multiphysics simulations has been identified as a valuable simulation capability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Model of Mixing Layer With Multicomponent Evaporating Drops

A mathematical model of a three-dimensional mixing layer laden with evaporating fuel drops composed of many chemical species has been derived. The study is motivated by the fact that typical real petroleum fuels contain hundreds of chemical species. Previously, for the sake of computational efficiency, spray studies were performed using either models based on a single representative species or models based on surrogate fuels of at most 15 species. The present multicomponent model makes it possible to perform more realistic simulations by accounting for hundreds of chemical species in a computationally efficient manner. The model is used to perform Direct Numerical Simulations in continuing studies directed toward understanding the behavior of liquid petroleum fuel sprays. The model includes governing equations formulated in an Eulerian and a Lagrangian reference frame for the gas and the drops, respectively. This representation is consistent with the expected volumetrically small loading of the drops in gas (of the order of 10 3), although the mass loading can be substantial because of the high ratio (of the order of 103) between the densities of liquid and gas. The drops are treated as point sources of mass, momentum, and energy; this representation is consistent with the drop size being smaller than the Kolmogorov scale. Unsteady drag, added-mass effects, Basset history forces, and collisions between the drops are neglected, and the gas is assumed calorically perfect. The model incorporates the concept of continuous thermodynamics, according to which the chemical composition of a fuel is described probabilistically, by use of a distribution function. Distribution functions generally depend on many parameters. However, for mixtures of homologous species, the distribution can be approximated with acceptable accuracy as a sole function of the molecular weight. The mixing layer is initially laden with drops in its lower stream, and the drops are colder than the gas. Drop evaporation leads to a change in the gas-phase composition, which, like the composition of the drops, is described in a probabilistic manner

Bellan, Josette↗

Inferring the temperature profile of the radiative shock in the COAX experiment with shock radiography, Dante, and spectral temperature diagnostics

Predicting and modeling the behavior of experiments with radiation waves propagating through low-density foams require a detailed quantification of the numerous uncertainties present. In regimes where a prominent radiative shock is produced, key dynamical features include the shock position, temperature, and curvature and the spatial distribution and temperature of the corresponding supersonic radiation wave. The COAX experimental platform is designed to constrain numerical models of such a radiative shock propagating through a low-density foam by employing radiography for spatial and shock information, Dante for characterizing the x-ray flux from the indirectly driven target, and a novel spectral diagnostic designed to probe the temperature profile of the wave. In this work, we model COAX with parameterized 2D simulations and a Hohlraum-laser modeling package to study uncertainties in diagnosing the experiment. The inferred temperature profile of the COAX radiation transport experiments has been shown to differ from simulations more than expected from drive uncertainties that have been constrained by simultaneous soft x-ray flux and radiography measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The MOOSE fluid properties module

The Fluid Properties module within the Multiphysics Object-Oriented Simulation Environment (MOOSE) is used to compute fluid properties for numerous applications, ranging from nuclear reactor thermal hydraulics to geothermal energy. Those applications drove the development of the module to enable numerous different fluid equations of states, property lookups with primitive and conserved flow variable to cater to pressure and density-driven solvers, and an object-oriented design facilitating expansion and maintenance. Each fluid property is implemented in its own class but inherits capabilities such as automatic differentiation, automated out-of-bounds handling or variable conversion capabilities. Here, this paper presents the module, its design, its user and developer interface, its content in terms of fluids and properties, and several of its applications showing its major role in the MOOSE simulation ecosystem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High Performance Parallel Computational Nanotechnology

At a recent press conference, NASA Administrator Dan Goldin encouraged NASA Ames Research Center to take a lead role in promoting research and development of advanced, high-performance computer technology, including nanotechnology. Manufacturers of leading-edge microprocessors currently perform large-scale simulations in the design and verification of semiconductor devices and microprocessors. Recently, the need for this intensive simulation and modeling analysis has greatly increased, due in part to the ever-increasing complexity of these devices, as well as the lessons of experiences such as the Pentium fiasco. Simulation, modeling, testing, and validation will be even more important for designing molecular computers because of the complex specification of millions of atoms, thousands of assembly steps, as well as the simulation and modeling needed to ensure reliable, robust and efficient fabrication of the molecular devices. The software for this capacity does not exist today, but it can be extrapolated from the software currently used in molecular modeling for other applications: semi-empirical methods, ab initio methods, self-consistent field methods, Hartree-Fock methods, molecular mechanics; and simulation methods for diamondoid structures. In as much as it seems clear that the application of such methods in nanotechnology will require powerful, highly powerful systems, this talk will discuss techniques and issues for performing these types of computations on parallel systems. We will describe system design issues (memory, I/O, mass storage, operating system requirements, special user interface issues, interconnects, bandwidths, and programming languages) involved in parallel methods for scalable classical, semiclassical, quantum, molecular mechanics, and continuum models; molecular nanotechnology computer-aided designs (NanoCAD) techniques; visualization using virtual reality techniques of structural models and assembly sequences; software required to control mini robotic manipulators for positional control; scalable numerical algorithms for reliability, verifications and testability. There appears no fundamental obstacle to simulating molecular compilers and molecular computers on high performance parallel computers, just as the Boeing 777 was simulated on a computer before manufacturing it.

Saini, Subhash↗

Numerical simulation of a relativistic magnetron using a fluid electron model

An approach to numerically modeling relativistic magnetrons, in which the electrons are represented with a relativistic fluid, is described. A principal effect in the operation of a magnetron is space-charge-limited (SCL) emission of electrons from the cathode. We have developed an approximate SCL emission boundary condition for the fluid electron model. This boundary condition prescribes the flux of electrons as a function of the normal component of the electric field on the boundary. Here we show the results of a benchmarking activity that applies the fluid SCL boundary condition to the one-dimensional Child–Langmuir diode problem and a canonical two-dimensional diode problem. Simulation results for a two-dimensional A6 magnetron are then presented. Computed bunching of the electron cloud occurs and coincides with significant microwave power generation. Numerical convergence of the solution is considered. Sharp gradients in the solution quantities at the diocotron resonance, spanning an interval of three to four grid cells in the most well-resolved case, are present and likely affect convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Baycal

Bayesian Model Calibration (BayCal) toolkit is a software plugin for Risk Analysis Virtual Environment (RAVEN) framework, arming at inversely quantifying the uncertainties associated with simulation model parameters based on available experiment data. BayCal seeks statistical inference of the uncertain input parameters that are consistent with the available measurement data or observed data. The unique feature of BayCal is the capability to be linked with RAVEN to build corresponding calibration workflows for complex multi-physics simulations. In addition to be able to use the machine learning capability of RAVEN to significantly reduce the computational cost of expensive simulation models, another distinctive feature of BayCal is the capability to deal with high-dimensional correlated model outputs, such as time series observations at multiple locations, via principal component analysis (PCA) technique.

Wang, Congjian↗

Reduced-order modeling of near-field THMC coupled processes for nuclear waste repositories in shale

Performance assessment (PA) of geologic radioactive waste repositories requires three-dimensional simulation of highly nonlinear, thermo-hydro-mechanical-chemical (THMC), multiphase flow and transport processes across many kilometers and over tens to hundreds of thousands of years. Integrating the effects of a near-field geomechanical process (i.e. buffer swelling) into coupled THC simulations through reduced-order modeling, rather than through fully coupled geomechanics, can reduce the dimensionality of the problem and improve computational efficiency. In this study, PFLOTRAN simulations model a single waste package in a shale host rock repository, where re-saturation of a bentonite buffer causes the buffer to swell and exert stress on a highly fractured disturbed rock zone (DRZ). Three types of stress-dependent permeability functions (exponential, modified cubic, and Two-part Hooke’s law models) are implemented to describe mechanical characteristics of the system. Overall, our modeling study suggests that compressing fractures reduces DRZ permeability, which could influence the rate of radionuclide transport and exchange with corrosive species in host rock groundwater that could accelerate waste package degradation. Less permeable shale host rock delays buffer swelling, consequently retarding DRZ permeability reduction as well as chemical transport within the barrier system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ordering Unstructured Meshes for Sparse Matrix Computations on Leading Parallel Systems

The ability of computers to solve hitherto intractable problems and simulate complex processes using mathematical models makes them an indispensable part of modern science and engineering. Computer simulations of large-scale realistic applications usually require solving a set of non-linear partial differential equations (PDES) over a finite region. For example, one thrust area in the DOE Grand Challenge projects is to design future accelerators such as the SpaHation Neutron Source (SNS). Our colleagues at SLAC need to model complex RFQ cavities with large aspect ratios. Unstructured grids are currently used to resolve the small features in a large computational domain; dynamic mesh adaptation will be added in the future for additional efficiency. The PDEs for electromagnetics are discretized by the FEM method, which leads to a generalized eigenvalue problem Kx = AMx, where K and M are the stiffness and mass matrices, and are very sparse. In a typical cavity model, the number of degrees of freedom is about one million. For such large eigenproblems, direct solution techniques quickly reach the memory limits. Instead, the most widely-used methods are Krylov subspace methods, such as Lanczos or Jacobi-Davidson. In all the Krylov-based algorithms, sparse matrix-vector multiplication (SPMV) must be performed repeatedly. Therefore, the efficiency of SPMV usually determines the eigensolver speed. SPMV is also one of the most heavily used kernels in large-scale numerical simulations.

Oliker, Leonid↗

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.↗

Atmospheric modeling related to Thematic Mapper scan geometry

A simulation study was carried out to characterize atmospheric effects in LANDSAT-D Thematic Mapper data. In particular, the objective was to determine if any differences would result from using a linear vs. a conical scanning geometry. Insight also was gained about the overall effect of the atmosphere on Thematic Mapper signals, together with the effects of time of day. An added analysis was made of the geometric potential for direct specular reflections (sun glint). The ERIM multispectral system simulation model was used to compute inband Thematic Mapper radiances, taking into account sensor, atmospheric, and surface characteristics. Separate analyses were carried out for the thermal band and seven bands defined in the reflective spectral region. Reflective-region radiances were computed for 40 deg N, 0 deg, and 40 deg S latitudes; June, Mar., and Dec. days; and 9:30 and 11:00 AM solar times for both linear and conical scan modes. Also, accurate simulations of solar and viewing geometries throughout Thematic Mapper orbits were made. It is shown that the atmosphere plays an important role in determining Thematic Mapper radiances, with atmospheric path radiance being the major component of total radiances for short wavelengths and decreasing in importance as wavelength increases. Path radiance is shown to depend heavily on the direct radiation scattering angle and on haze content. Scan-angle-dependent variations were shown to be substantial, especially for the short-wavelength bands.

Malila, W. A.↗

Graphical User Interface (GUI) Implementation for Agent-Based Microbial Radiobiology Model

Sending human life past the Low Earth Orbit (LEO) to explore the Moon and Mars will be challenging. The Earth’s magnetic field naturally protects life from deep-space particle radiation such as Galactic Cosmic Rays (GCR) and Solar Particle Events (SPE); these will pose health risks to humans in deep space. Research has been done to investigate these effects, like BioSentinel, the first biological CubeSat to fly beyond the LEO, designed to culture yeast in a microfluidic device and record optical measurements of growth and metabolism. However, experiments can only report cell damage as bulk growth curves, while deep-space radiation causes damage that is heterogeneous among individual cells. AMMPER is an open-source, agent-based, computational model coded in Python to simulate the effects of deep-space radiation on individual yeast cells (Saccharomyces cerevisiae) to facilitate interpretation of biological radiation experiments. Version 1.0 of the code ran in a command line interface (CLI), limiting use to those familiar with modularization, object-oriented programming, and computational models. Here we present a graphical user interface (GUI) for AMMPER to increase its accessibility. GUI development included converting input points and UI files, designing an application and logo, and expanding program packages. Additionally, we added optical assistance that corresponded with simulation parameters, which included simulation type, cell type, ROS model, and radiation dosage, as well as customizable display and file exportation features. Following a pilot testing period, its structure was updated further to enhance abilities, adding increased runs, video visualization, data plotting, and an educational/tutorial component. Future work will include creating a bit installer and runtime environment for AMMPER. Ultimately, the creation of the GUI has two main goals: to facilitate the integration of computational models into the work of researchers in microbial radiobiology, and to act as an interactive and visual resource for space biology education.

yeast↗