Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 469 records · Page 26

Initial Demonstration of New Griffin Technologies for Simulating the Running-In Phase of Pebble Bed Reactors

Griffin is a reactor multiphysics modeling application based on MOOSE (Multiphysics Object-Oriented Simulation Environment) and specifically targeting transient modeling of advanced reactors. Griffin has been used recently to model pebble-bed reactors for the Nuclear Regulatory Commission (NRC) Office of Nuclear Regulatory Research and the Advanced Reactor Technology program. This modeling work has focused thus far on the direct calculation of equilibrium cores. This report documents an initial demonstration of a new running-in simulation capability. The new running-in capability is verified using the existing direct equilibrium core calculation capability. A simplified pebble-bed reactor model is then used to demonstrate the running-in simulation capability. This demonstration shows that Griffin is able to simulate years of operation during the running-in phase efficiently with each depletion step taking only several seconds. Two new technologies are also presented in this report which have been developed in Griffin that will be essential for improved accuracy both of the direct equilibrium core computation and the new running-in simulation capability. The first technology is an online cross section generation capability specifically targeted for pebble-bed reactors. This will improve the accuracy of the depletion calculation as the cross sections are generated at the exact core status. This also avoids the difficult step of pre-generating a separate standalone multigroup cross section set. Secondly, a newly implemented discretization for discontinuous finite element method (DFEM) SN transport in cylindrical (RZ) coordinates, which can be solved efficiently using the existing SN sweep solver, is discussed and some results are shown demonstrating the usefulness of the additional accuracy transport provides over a diffusion approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Comprehensive eVTOL Performance Evaluation Framework in Urban Air Mobility

In this paper, we developed an open-source simulation framework for the evaluation of electric vertical takeoff and landing vehicles (eVTOLs) in the context of Unmanned Traffic Management (UTM) and under the concept of Urban Air Mobility (UAM). Unlike most existing studies, the proposed framework combines the utilization of UTM and eVTOLs to develop a realistic UAM testing platform. For this purpose, we first develop an UTM simulator to simulate the real-world UAM environment. Then, instead of using a simplified eVOTL model, a high-fidelity eVTOL design tool, namely SUAVE, is employed and an dilation sub-module is introduced to bridge the gap between the UTM simulator and SUAVE eVTOL performance evaluation tool to elaborate the complete mission profile. Based on the developed simulation framework, experiments are conducted and the results are presented to analyze the performance of eVTOLs in the UAM environment.

Mrinmoy Sarkar↗

Eulerian simulation of complex suspensions and biolocomotion in three dimensions

Significance Fluid–structure interactions are ubiquitous in many natural and man-made environments. They are difficult to study analytically, and therefore accurate and flexible computational methods are an indispensable tool in the field. Typically, fluids are simulated with a fixed background computational mesh, whereas a solid is simulated with a mesh that moves with it, making it challenging to couple the two. Here we develop a three-dimensional computational method where both fluid and solid can be represented on a fixed computational grid, which simplifies the coupling between the two phases considerably. Our method can simulate scenarios that are difficult to do with other methods, such as complex suspensions containing soft particles that are both heavier and lighter than the fluid.

3D fluid-structure interaction↗

Analytic insights into nonlocal energy transport: Steady state Fokker Planck theory in arbitrary Z plasmas

The generation of energetic electrons in laser fusion in an important issue. The electrons may either arise from a laser plasma instability, or from the uncoupled high temperature tail of a Maxwellian distribution. To study these in a laser fusion context, it is important to find a method accurate enough to be useful, and simple enough to be incorporated into a radiation hydrodynamics numerical simulation, the main workhorse for studying the laser fusion target. That is why analytic insights become important, they allow one to simplify the Fokker Planck theory so that a solution of it can be incorporated into a radiation hydrodynamic simulation. This work develops and analyzes a steady state Fokker Planck theory for plasmas of arbitrary Z. It develops a method of solving the simplified Fokker Planck method with a technique called sparse eigenfunction analysis. As a result, this method appears to work reasonably well when compared to the experimental results from the Rochester/NIF on plastic spherical targets with and without a silicon layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coronal shock waves

The formation of magnetohydrodynamic (MHD) shocks near the leading edge of Coronal Mass Ejections (CME's) is considered using results from both shock theory and numerical simulations. By using just the shock theory results and a simplified model corona, the following is shown: a slow shock front should be concave upward (away from the solar surface); a configuration containing intermediate shocks should also be concave upward near the center of the CME and concave downward at some distance from the CME centerline; and a fast shock configuration should be concave downward. These results are verified with numerical simulations of the MHD equations for the propagation of shocks through a quiescent coronal streamer. By examining the characteristic wave speeds in the ambient streamer, it is estimated that the various shock configurations should form for CME speeds within the following speed increments: slow shocks, 200 to 300 km/sec; intermediate shocks, 300 to 900 km/sec; and fast shocks, greater than 900 km/sec.

Steinolfson, Richard S.↗

Improving Predictive Capability in REHEDS Simulations with Fast, Accurate, and Consistent Non-Equilibrium Material Properties

Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.

36 MATERIALS SCIENCE↗

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The problem of neuron modeling

Neuron modeling - simulation of neuronal behavior, model based on Hodgkin and Huxley equations, and simplified model using properties of pulse frequency modulated systems

BEHAVIOR↗

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercia2 CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and for these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were m to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire lee wing without the wheel well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was normal to the spar. It included all of the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was Found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas↗

Understanding Columbia's Reentry Problems

Soon after the Space Shuttle Columbia accident occurred last year, a group of CFD analysts from NASA centers and private industry was organized to help determine the cause of the accident. This group was under the direction of the Applied Aeroscience and CFD Branch of the Aeroscience and Flight Mechanics Division at the Johnson Space Center. For external flow simulations, noncommercial CFD codes that specialize in hypersonic or high Mach number flows were used. These tools were used to determine heating rates, pressures, and temperatures for a large number of vehicle damage scenarios. Lockheed Martin Space Operations was called upon to provide CFD support in the area of internal flows within the shuttle wing cavity, and or these simulations, FLUENT 6.1 was chosen. Two large-scale, simplified models were run to understand the flow patterns once a breach of the internal wing cavity was initiated. The results were primarily used to visualize flow patterns within the wing cavity. The first CFD model included the entire left wing without the wheel. well cavity. The purpose of the first model, which did not include the reinforced carbon-carbon (RCC) cavity along the wing leading edge, was to visualize the flow field within the wing cavity immediately after the leading edge spar breach, This model assumed that the flow coming into the wing cavity was nom1 to the spar. It included all o f the primary vents that allow for flow between the main cavities of the wing. A six-inch diameter hole was modeled in the spar at the approximate location where the spar breach was judged to have occurred, which was between RCC panels 8 and 9. The results of the modeling showed that at this location, the high temperature, high velocity gas stream entering the wing cavity impinged on the outboard wheel well cavity. Instrumentation in the Shuttle wheel well cavity registered abnormal temperatures during reentry, so the FLUENT results helped support the conclusion of the accident investigation team that the spar breach was in the RCC panel 8-9 area, and that the initial spar breach was likely entering the wing cavity normal to the spar. This model also showed that the flow entering the wing cavity tended to swirl within the cavity just outboard of the wheel well, and did not initially penetrate further into the rear cavities of the wing. The second CFD model was a 2-D simulation of the left wing cavity and the RCC cavity. It was used to visualize the flow through the RCC breach, through the wing spar breach, and into the wing cavity directly outboard of the wheel well. The purpose of this model was to verify whether or not it was possible for the flow to come into the wing cavity normal to the leading edge spar. The results from this 2-D model showed that the internal structure behind the RCC panels probably deflected the flow entering the RCC cavity so that it impinged normal to the spar. As in the 3-D model, this deflected flow stream was found to impinge on the wheel well outer wall. This model again supported the conclusions regarding the location of the spar breach and how the flow behaved inside the wing cavity.

Paul, Thomas H.↗

How Should Life Support Be Modeled and Simulated?

Why do most space life support research groups build and investigate large models for systems simulation? The need for them seems accepted, but are we asking the right questions and solving the real problems? The modeling results leave many questions unanswered. How then should space life support be modeled and simulated? Life support system research and development uses modeling and simulation to study dynamic behavior as part of systems engineering and analysis. It is used to size material flows and buffers and plan contingent operations. A DoD sponsored study used the systems engineering approach to define a set of best practices for modeling and simulation. These best practices describe a systems engineering process of developing and validating requirements, defining and analyzing the model concept, and designing and testing the model. Other general principles for modeling and simulation are presented. Some specific additional advice includes performing a static analysis before developing a dynamic simulation, applying the mass and energy conservation laws, modeling on the appropriate system level, using simplified subsystem representations, designing the model to solve a specific problem, and testing the model on several different problems. Modeling and simulation is necessary in life support design but many problems are outside its scope.

Jones, Harry W.↗

Hurricane‐Like Vortices in Conditionally Unstable Moist Convection

Abstract This study investigates the emergence of hurricane‐like vortices in idealized simulations of rotating moist convection. A Boussinesq atmosphere with simplified thermodynamics for phase transitions is forced by prescribing the temperature and humidity at the upper and lower boundaries. The governing equations are solved numerically using a variable‐density incompressible Navier‐Stokes solver with adaptive mesh refinement to explore the behavior of moist convection under a broad range of conditions. In the absence of rotation, convection aggregates into active patches separated by large unsaturated regions. Rotation modulates this statistical equilibrium state so that the self‐aggregated convection organizes hurricane‐like vortices. The warm and saturated air converges to the center of the vortices, and the latent heat released through the upwelling, forms the warm core structure. These hurricane‐like vortices share characteristics similar to tropical cyclones in the earth's atmosphere. The hurricane‐like vortices occur under conditionally unstable conditions where the potential energy given at the boundaries is large enough, corresponding to a moderate rate of rotation. This regime shares many similar characteristics to the tropical atmosphere indicating that the formation of intense meso‐scale vortices is a general characteristic of rotating moist convection. The model used here does not include any interactions with radiation, wind‐evaporation feedback, or cloud microphysics, indicating that, while these processes may be relevant for tropical cyclogenesis in the Earth atmosphere, they are not its primary cause. Instead, our results confirm that the formation and maintenance of hurricane‐like vortices involve a combination of atmospheric dynamics under the presence of rotation and of phase transitions.

54 ENVIRONMENTAL SCIENCES↗

Multi-fidelity modeling to predict the rheological properties of fiber suspensions

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications like biofuel production. The 3D numerical simulations of the suspension of fibers are often computationally expensive and time-consuming. Machine learning methods such as neural networks can simplify the prediction of rheological behavior; however, they require a relatively large training data set. Multi -fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified physical equations, can lead to optimized predictions. Here, we focus on a neural network with two levels of fidelity, i.e., high and low fidelity networks. To produce high-fidelity data, we perform direct numerical simulations to model the fibers as one-dimensional inextensible slender bodies that obey the Euler- Bernoulli beam equation. The Navier-Stokes equations govern the suspended fluid, and an immersed boundary method is used to couple the fluid and solid motion. The low-fidelity data is produced by using constitutive equations. Noticeable improvements have been observed in the accuracy of predicting the rheological behavior when a multi-fidelity network is used compared to the single-fidelity network.

Boodaghidizaji, Miad↗

Detached Eddy Simulation of the UH-60 Rotor Wake Using Adaptive Mesh Refinement

Time-dependent Navier-Stokes flow simulations have been carried out for a UH-60 rotor with simplified hub in forward flight and hover flight conditions. Flexible rotor blades and flight trim conditions are modeled and established by loosely coupling the OVERFLOW Computational Fluid Dynamics (CFD) code with the CAMRAD II helicopter comprehensive code. High order spatial differences, Adaptive Mesh Refinement (AMR), and Detached Eddy Simulation (DES) are used to obtain highly resolved vortex wakes, where the largest turbulent structures are captured. Special attention is directed towards ensuring the dual time accuracy is within the asymptotic range, and verifying the loose coupling convergence process using AMR. The AMR/DES simulation produced vortical worms for forward flight and hover conditions, similar to previous results obtained for the TRAM rotor in hover. AMR proved to be an efficient means to capture a rotor wake without a priori knowledge of the wake shape.

Chaderjian, Neal M.↗

Cognitive simulation models for inertial confinement fusion: Combining simulation and experimental data

The design space for inertial confinement fusion (ICF) experiments is vast, and experiments are extremely expensive. Researchers rely heavily on computer simulations to explore the design space in search of high-performing implosions. However, ICF multiphysics codes must make simplifying assumptions, and thus deviate from experimental measurements for complex implosions. For more effective design and investigation, simulations require input from past experimental data to better predict future performance. In this work, we describe a cognitive simulation method for combining simulation and experimental data into a common, predictive model. This method leverages a machine learning technique called “transfer learning,” the process of taking a model trained to solve one task, and partially retraining it on a sparse dataset to solve a different, but related task. In the context of ICF design, neural network models are trained on large simulation databases and partially retrained on experimental data, producing models that are far more accurate than simulations alone. Here, we demonstrate improved model performance for a range of ICF experiments at the National Ignition Facility and predict the outcome of recent experiments with less than 10% error for several key observables. We discuss how the methods might be used to carry out a data-driven experimental campaign to optimize performance, illustrating the key product—models that become increasingly accurate as data are acquired.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TERRA Battery Thermal Control Anomaly - Simulation and Corrective Actions

The TERRA spacecraft was launched in December 1999 from Vandenberg Air Force Base, becoming the flagship of NASA's Earth Observing System program to gather data on how the planet's processes create climate. Originally planned as a 5 year mission, it still provides valuable science data after nearly 10 years on orbit. On October 13th, 2009 at 16:23z following a routine inclination maneuver, TERRA experienced a battery cell failure and a simultaneous failure of several battery heater control circuits used to maintain cell temperatures and gradients within the battery. With several cells nearing the minimum survival temperature, preventing the electrolyte from freezing was the first priority. After several reset attempts and power cycling of the control electronics failed to reestablish control authority on the primary side of the controller, it was switched to the redundant side, but anomalous performance again prevented full heater control of the battery cells. As the investigation into the cause of the anomaly and corrective action continued, a battery thermal model was developed to be used in determining the control ability remaining and to simulate and assess corrective actions. Although no thermal model or detailed reference data of the battery was available, sufficient information was found to allow a simplified model to be constructed, correlated against pre-anomaly telemetry, and used to simulate the thermal behavior at several points after the anomaly. It was then used to simulate subsequent corrective actions to assess their impact on cell temperatures. This paper describes the rapid development of this thermal model, including correlation to flight data before and after the anomaly., along with a comparative assessment of the analysis results used to interpret the telemetry to determine the extent of damage to the thermal control hardware, with near-term corrective actions and long-term operations plan to overcome the anomaly.

Grob, Eric W.↗

Simulation model of a twin-tail, high performance airplane

The mathematical model and associated computer program to simulate a twin-tailed high performance fighter airplane (McDonnell Douglas F/A-18) are described. The simulation program is written in the Advanced Continuous Simulation Language. The simulation math model includes the nonlinear six degree-of-freedom rigid-body equations, an engine model, sensors, and first order actuators with rate and position limiting. A simplified form of the F/A-18 digital control laws (version 8.3.3) are implemented. The simulated control law includes only inner loop augmentation in the up and away flight mode. The aerodynamic forces and moments are calculated from a wind-tunnel-derived database using table look-ups with linear interpolation. The aerodynamic database has an angle-of-attack range of -10 to +90 and a sideslip range of -20 to +20 degrees. The effects of elastic deformation are incorporated in a quasi-static-elastic manner. Elastic degrees of freedom are not actively simulated. In the engine model, the throttle-commanded steady-state thrust level and the dynamic response characteristics of the engine are based on airflow rate as determined from a table look-up. Afterburner dynamics are switched in at a threshold based on the engine airflow and commanded thrust.

Buttrill, Carey S.↗