Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation modelling”

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 1,225 records · Page 68

NASA ETM Modeling and Simulation Upper E Traffic Management Meeting December 2020

This is a slide set to be used as part of a meeting series with members of a working group aimed at the development of a concept that addresses needs and gaps in the management of high altitude airspace operations. This concept leverages elements developed through the UAS Traffic Management project with respect to a cooperative, service-based approach that provides services and capabilities in areas (e.g., Upper Class E airspace) that currently receive no or limited service from Air Traffic Control. This concept is meant to provide a safe, fair, and scalable approach to management of Upper Class E operations that reduces the burden on ATC while providing the flexibility and access desired by current and new users of the airspace. This set of slides includes an overview of the past accomplishments, current project structure, and a simulation roadmap for planned work ahead.

NASA ETM modeling and simulation↗

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗

Modeling and Simulation of Air-Source CO2 Heat Pump Water Heater

Carbon dioxide (CO2) has been widely used as working fluid for the vapor-compression refrigeration systems in large marine device. Due to the potential energy efficiency and the favorable environmental properties of CO2 as a working fluid, CO2 heat pump water heater (HPWH) systems are regarded a promising technology for centralized domestic hot water (DHW) heating in residential and commercial buildings. However, there is still at the early stage of appropriately optimizing and improving the energy performance of CO2 HPWH. This requires CO2 HPWH simulation tools capable of capturing the accurate impact of the emerging compressor, throttle device, and heat exchanger technology on CO2 heat transfer and energy efficiency. In this study, high efficiency components (compressors, pumps, fans, heat exchangers) were identified and applied to the state-of-art CO2 HPWH designs and analyzed their performance by using numerical simulation. This was done by simulating the performance of CO2 HPWH using ACMODEL design model combined with the component models developed at Oak Ridge National Laboratory (ORNL) for orifice tube, map-based compressor, and tube-in-tube gas cooler. ACMODEL is an equipment design model for CO2-based air conditioners and heat pumps developed by Purdue University to account for the details of each component. The simulated CO2 HPWH performance was then compared with the heat pump water heater using conventional refrigerants.

Gao, Zhiming↗

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3, documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation is to describe the application of the hydraulic property fields and recharge parameters documented in ECF-HANFORD-20-0027, Null Space Monte Carlo Evaluation of the Plateau to River Model, to the CA flow and fate and transport simulation results to quantify the uncertainty in the simulated results due to input parameter selection. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, models lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. The result of NSCM analysis is a set of F&T simulations that provide an estimate of the range of possible outcomes that are used to quantify the uncertainty in simulated concentrations produced using the base case simulations. The simulated concentrations from all simulations will support calculation of the uncertainty of the total dose calculated in a separate calculation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computer modeling and simulation of a 20kHz ac distribution system for Space Station

A computer model of a 20 kHz, ac distribution testbed for Space Station is presented. The system consists of six resonant inverters, a one-hundred-meter transmission line, and three load receivers: a dc receiver, a bidirectional receiver, and an ac receiver. A model library is generated characterizing all system components. The system's local and global behaviors are investigated using the EASY5 dynamic analysis program.

Tsai, Fu-Sheng↗

Automated Modeling and Simulation Using the Bond Graph Method for the Aerospace Industry

Bond graph modeling was originally developed in the late 1950s by the late Prof. Henry M. Paynter of M.I.T. Prof. Paynter acted well before his time as the main advantage of his creation, other than the modeling insight that it provides and the ability of effectively dealing with Mechatronics, came into fruition only with the recent advent of modern computer technology and the tools derived as a result of it, including symbolic manipulation, MATLAB, and SIMULINK and the Computer Aided Modeling Program (CAMPG). Thus, only recently have these tools been available allowing one to fully utilize the advantages that the bond graph method has to offer. The purpose of this paper is to help fill the knowledge void concerning its use of bond graphs in the aerospace industry. The paper first presents simple examples to serve as a tutorial on bond graphs for those not familiar with the technique. The reader is given the basic understanding needed to appreciate the applications that follow. After that, several aerospace applications are developed such as modeling of an arresting system for aircraft carrier landings, suspension models used for landing gears and multibody dynamics. The paper presents also an update on NASA's progress in modeling the International Space Station (ISS) using bond graph techniques, and an advanced actuation system utilizing shape memory alloys. The later covers the Mechatronics advantages of the bond graph method, applications that simultaneously involves mechanical, hydraulic, thermal, and electrical subsystem modeling.

Granda, Jose J.↗

Power system model parameter conditioning tool

A power system model parameter conditioning tool including a server control processor in communication with phasor measurement unit monitored data records of multiple disturbance events, a model calibration unit providing event screening, power system model simulation, and simultaneous tuning of model parameters. The model calibration performing a simulation using default model parameters, the processor comparing the simulation results to the monitored data. If the prediction is within threshold, then terminating conditioning; else performing parameter identifiability analysis to determine differing effects of various model parameters on power system model accuracy, selecting a parameter set causing a degradation in power system model prediction, and updating the default model parameters corresponding to members of the parameter set with values selected to reduce the degradation. A method and a non-transitory computer readable medium are also disclosed.

Baone, Chaitanya Ashok↗

Environmental Thermodynamic Effects on Future Changes of Mesoscale Convective Systems and Isolated Deep Convection Over the US: A Case Study

Mesoscale convective systems (MCSs) and isolated deep convection (IDC) significantly influence local weather conditions and the hydrological cycle across the United States. Their future changes under warming is critical to the society. However, the exact impacts of thermodynamic environmental changes on them are difficult to resolve using climate model simulations which integrate both dynamic and thermodynamic factors. This study explores a theoretical modeling approach to isolate environmental thermodynamic effects on their future changes through a case study. For a 50-day period during the 2020 summer, a pair of convection-permitting model simulations indicate contrasting changes in MCSs and IDC between inland and coastal regions under pseudo global warming. Driven by the thermodynamic environments of these simulations, a single-column parcel model indicates a decrease in the frequency of IDC occurrences, along with increases in duration and precipitation amount under warming, attributable to rising most unstable convective available potential energy (MUCAPE), convective inhibition (MUCIN), and precipitable water (PW). A multi-column parcel model reveals contrasting changes in the frequency and mean area of MCSs between inland and coastal regions, underscoring the increase in mean MUCIN over inland regions and increases in mean MUCAPE and PW over coastal regions. The increase in mean MCS area over inland regions is linked to the interplay between accelerated gust fronts and enhanced subsidence strength, which is further traced to the unchanged mean MUCIN but increased mean MUCAPE in large-scale environments. In conclusion, these case study results suggest a major role of thermodynamic environmental changes in controlling characteristics of MCSs and IDC under warming.

Mesoscale convective system↗

High-Fidelity Roadway Modeling and Simulation

Roads are an essential feature in our daily lives. With the advances in computing technologies, 2D and 3D road models are employed in many applications, such as computer games and virtual environments. Traditional road models were generated by professional artists manually using modeling software tools such as Maya and 3ds Max. This approach requires both highly specialized and sophisticated skills and massive manual labor. Automatic road generation based on procedural modeling can create road models using specially designed computer algorithms or procedures, reducing the tedious manual editing needed for road modeling dramatically. But most existing procedural modeling methods for road generation put emphasis on the visual effects of the generated roads, not the geometrical and architectural fidelity. This limitation seriously restricts the applicability of the generated road models. To address this problem, this paper proposes a high-fidelity roadway generation method that takes into account road design principles practiced by civil engineering professionals, and as a result, the generated roads can support not only general applications such as games and simulations in which roads are used as 3D assets, but also demanding civil engineering applications, which requires accurate geometrical models of roads. The inputs to the proposed method include road specifications, civil engineering road design rules, terrain information, and surrounding environment. Then the proposed method generates in real time 3D roads that have both high visual and geometrical fidelities. This paper discusses in details the procedures that convert 2D roads specified in shape files into 3D roads and civil engineering road design principles. The proposed method can be used in many applications that have stringent requirements on high precision 3D models, such as driving simulations and road design prototyping. Preliminary results demonstrate the effectiveness of the proposed method.

Wang, Jie↗

Intercomparison of Large-Eddy Simulations of Arctic Mixed-Phase Clouds: Importance of Ice Size Distribution Assumptions

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP), in agreement with earlier studies. In contrast to previous intercomparison studies, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSDs) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case. Sensitivity tests indicate LWP and IWP are much closer to the bin model simulations when a modified shape factor which is similar to that predicted by bin model simulation is used in bulk scheme. These results demonstrate the importance of representation of ice PSD in determining the partitioning of liquid and ice and the longevity of mixed-phase clouds.

clouds↗

Modelling and simulation of Space Station Freedom berthing dynamics and control

A large-angle, flexible, multibody, dynamic modeling capability has been developed to help validate numerical simulations of the dynamic motion and control forces which occur during berthing of Space Station Freedom to the Shuttle Orbiter in the early assembly flights. This paper outlines the dynamics and control of the station, the attached Shuttle Remote Manipulator System, and the orbiter. The simulation tool developed for the analysis is described and the results of two simulations are presented. The first is a simulated maneuver from a gravity-gradient attitude to a torque equilibrium attitude using the station reaction control jets. The second simulation is the berthing of the station to the orbiter with the station control moment gyros actively maintaining an estimated torque equilibrium attitude. The influence of the elastic dynamic behavior of the station and of the Remote Manipulator System on the attitude control of the station/orbiter system during each maneuver was investigated. The flexibility of the station and the arm were found to have only a minor influence on the attitude control of the system during the maneuvers.

Cooper, Paul A.↗

MBMS1.0: An Open-Source Code for Modeling and Simulation of Membrane-Based Dehumidification and Energy Recovery

Membrane-based dehumidification is currently being considered as a promising solution for the building application due to its low cost and very limited energy consumption. Developing a simple and efficient open-source code simulation tool is important for boosting the optimization and evaluation of such device in HVAC community. This paper reports a first-order physics based model which accounts for the fundamental heat and mass transfer of humid-air vapor at feed side to flow stream at permeate side. The current model comprises two membrane mass transfer submodels (i.e. microstructure model and performance map model); and it adopts a segment-by-segment methodology for discretizing heat and mass transfer governing equations. The model is capable of simulating both dehumidifiers and energy recovery ventilators with parallel-flow cross-flow, and counter-flow configurations. The model was validated with the measurements at appropriate device. The practices in dehumidification and energy recovery exchangers are also discussed. The model and open-source codes are expected to become a solid fundament for developing a more comprehensive and accurate membrane-based dehumidification in the future.

Gao, Zhiming↗

Studsvik R2 materials test reactor Ad Hoc depletion strategy for the derivation of the fuel isotopic composition of the MPCMIV benchmark

The Ad Hoc Depletion Strategy elaborated by the NINE company, developed in support of the organization of the MPCMIV (Multi-physics Pellet Cladding Mechanical Interaction Validation) benchmark input and output specifications, is presented. This work aims at illustrating the strategy itself and then showing the results obtained with its application over the Studsvik R2 Testing Reactor, which is analyzed in the benchmark. The objective of the application of the strategy is to compute the fuel elements isotopic compositions at the beginning of some core loadings of interest for the benchmark. To this objective, it is necessary to implement first the simulation model of the three single assembly types and perform the infinite lattice depletions, then, to build the full core model and to perform the simulation of the core cycle. All the models and simulations were carried out with the use of the Monte Carlo particle transport code Serpent 2. Finally, the simulations results are assessed against Studsvik isotopic compositions of the fuel elements discharged from the R2 Reactor at the end of the core loading. Several assumptions were necessary during all the steps of the strategy, to overcome the lack of information regarding the core management. For this reason, the solution found at the end of the current analysis may not be completely optimized and further improvements regarding the model assumptions will be tested in a future work. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Modeling and Simulation of Power-Electronic Inverters in Analog Electronic Circuit Simulators

This paper demonstrates how equivalent-circuit representations of grid-following power-electronic inverters can be realized within a SPICE-based development environment using common circuit components and VerilogA code. This facilitates computationally lean simulations of inverter networks leveraging the strengths of SPICE in large-scale simulations. We validate the approach with time-domain simulations for a modified version of the IEEE standard 118-bus system modeled in Virtuoso (a SPICE-based solver). Simulation results are compared-focusing on accuracy and computation speed-with results from a commonly used power-electronics simulation package. We note a significant decrease in simulation time with comparable signal resolution when simulating the network in SPICE.

Billmeyer, Ryan↗

Experiments with the Mesoscale Atmospheric Simulation System (MASS) using the synthetic relative humidity

This study is intended to examine the impact of the synthetic relative humidity on the model simulation of mesoscale convective storm environment. The synthetic relative humidity is derived from the National Weather Services surface observations, and non-conventional sources including aircraft, radar, and satellite observations. The latter sources provide the mesoscale data of very high spatial and temporal resolution. The synthetic humidity data is used to complement the National Weather Services rawinsonde observations. It is believed that a realistic representation of initial moisture field in a mesoscale model is critical for the model simulation of thunderstorm development, and the formation of non-convective clouds as well as their effects on the surface energy budget. The impact will be investigated based on a real-data case study using the mesoscale atmospheric simulation system developed by Mesoscale Environmental Simulations Operations, Inc. The mesoscale atmospheric simulation system consists of objective analysis and initialization codes, and the coarse-mesh and fine-mesh dynamic prediction models. Both models are a three dimensional, primitive equation model containing the essential moist physics for simulating and forecasting mesoscale convective processes in the atmosphere. The modeling system is currently implemented at the Applied Meteorology Unit, Kennedy Space Center. Two procedures involving the synthetic relative humidity to define the model initial moisture fields are considered. It is proposed to perform several short-range (approximately 6 hours) comparative coarse-mesh simulation experiments with and without the synthetic data. They are aimed at revealing the model sensitivities should allow us both to refine the specification of the observational requirements, and to develop more accurate and efficient objective analysis schemes. The goal is to advance the MASS (Mesoscal Atmospheric Simulation System) modeling expertise so that the model output can provide reliable guidance for thunderstorm forecasting.

Chang, Chia-Bo↗

Multibody Modeling and Simulation for the Mars Phoenix Lander Entry, Descent and Landing

A multi-body flight simulation for the Phoenix Mars Lander has been developed that includes high fidelity six degree-of-freedom rigid-body models for the parachute and lander system. The simulation provides attitude and rate history predictions of all bodies throughout the flight, as well as loads on each of the connecting lines. In so doing, a realistic behavior of the descending parachute/lander system dynamics can be simulated that allows assessment of the Phoenix descent performance and identification of potential sensitivities for landing. This simulation provides a complete end-to-end capability of modeling the entire entry, descent, and landing sequence for the mission. Time histories of the parachute and lander aerodynamic angles are presented. The response of the lander system to various wind models and wind shears is shown to be acceptable. Monte Carlo simulation results are also presented.

Queen, Eric M.↗