Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling and simulation”

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,693 records · Page 94

Overview of modelling and simulation of metal powder bed fusion process at Lawrence Livermore National Laboratory

The metal laser powder bed fusion additive manufacturing process uses high power lasers to build parts layer upon layer by melting fine metal powders. Qualification of parts produced using this technology is broadly recognised as a significant challenge. Physics based process models have been identified as being foundational to qualification of additively manufactured metal parts. In the present article, a multiscale modelling strategy is described that will serve as the foundation upon which process control and part qualification can be built. This includes a model at the scale of the powder that simulates single track/single multilayer builds and provides powder bed and melt pool thermal data. A second model computationally builds a complete part and predicts manufactured properties (residual stress, dimensional accuracy) in three dimensions. Finally, modelling is tied to experiment through data mining.

36 MATERIALS SCIENCE↗

Modeling and Simulating Diffused Aeration for Hydrogen Removal from Expansion Tanks of Parabolic Trough Solar Thermal Power Plants

High-temperature decomposition of organic heat transfer fluids (HTFs) can cause hydrogen build-up in the receiver annulus of parabolic trough power plants. This build-up increases the receiver thermal losses and results in a decline in power output. Prior work has shown that removal of hydrogen from the expansion tank can be an effective mitigation strategy. This paper presents a modeling methodology and simulation results for passive and diffused aeration removal of hydrogen from the expansion tank. Focus is on estimating the mass transfer of hydrogen from the HTF into the headspace gas across the surface and bubble liquid/gas interfaces. Simulations for the operation conditions and expansion tank geometry of the Nevada Solar One power plant (located near Las Vegas, Nevada, USA) show that passive removal (surface mass transfer only) may be sufficient to obtain hydrogen removal rates of 1.7 x 10-4 mol/s (minimal removal rate necessary). A diffused aeration system (surface plus bubble mass transfer) with gas injection volume rates up to 100 L/s may be necessary to obtain a hydrogen removal rate of 2.0 x 10-4 mol/s (optimal removal rate).

14 SOLAR ENERGY↗

On the modelling and simulation of the dynamics and control of large flexible orbiting systems

This paper attempts to review the steps involved in the development of mathematical models that can be used to simulate the in-orbit dynamics of large flexible systems. The use of graph theoretic techniques can often be used to reduce the computational effort involved for calculating the eigenvalues of large ordered systems. Computer generated graphical techniques may provide additional insight into the understanding of elastic modal shape functions of complex systems. Finally the numerical techniques commonly used to develop shape and attitude control laws will be briefly reviewed.

Bainum, P. M.↗

Modeling and Simulation of Ternary Pumped Storage Hydropower for Power System Studies

As the deployment of wind and solar energy increases in the United States, energy storage (ES) will play an important role in future electric power grids. To help manage the variability from high penetration levels of renewable generation, ES can provide promising power/energy demand to coordinate with renewables generation, as a virtual power plant (VPP). There is an industry need for the capability in power system studies to model ternary pumped storage hydropower (T-PSH) that offers increased system benefits. This paper presents a comprehensive vendor-neutral dynamic model of T-PSH in GE's commercial software positive sequence load flow (PSLF). A new governor model is developed with detailed gate valve modeling and a shared-penstock function. Specifically, this model is designed to simulate the seamless transition among three operation modes: generation mode, pumping mode, and hydraulic short-circuit mode. The developed T-PSH model has been tested and validated on a 10-bus test system. A comparison study of T-PSH to conventional pumped storage hydropower has also been conducted on the Western Interconnection system.

13 HYDRO ENERGY↗

Simulation and modeling of homogeneous compressible turbulence under isotropic mean compression

Compressible homogeneous turbulence subjected to isotropic mean compression is simulated numerically for high and low turbulent Mach numbers, at various compression rates. It is found that at low Mach numbers, the effects of viscosity variations on the development of the turbulent kinetic energy can be significant. A model consistent with the invariance of the Navier-Stokes equations to spherical compression, that takes into account variable viscosity effects is proposed. At high Mach numbers, the distribution of energy between the acoustic and solenoid fields is found to depend upon initial conditions. The simulation results are also used to evaluate two-equation compressible turbulence models.

Coleman, G. N.↗

Development of Macroscale Models of UO 2 Fuel Sintering and Densification using Multiscale Modeling and Simulation

Light water reactor fuel pellets are fabricated using sintering to final densities of 95% or greater. During reactor operation, the porosity remaining in the fuel after fabrication decreases further due to irradiation-assisted densification. While empirical models have been developed to describe this densification process, a mechanistic model is needed as part of the ongoing work by the NEAMS program to develop a more predictive fuel performance code. In this work we will develop a phase field model of sintering of UO 2 in the MARMOT code, and validate it by comparing to published sintering data. We will then add the capability to capture irradiation effects into the model, and use it to develop a mechanistic model of densification that will go into the BISON code and add another essential piece to the microstructure-based materials models. The final step will be to add the effects of applied fields, to model field-assisted sintering of UO 2 . The results of the phase field model will be validated by comparing to data from field-assisted sintering. Tasks over three years: 1. Develop a sintering model for UO 2 in MARMOT 2. Expand model to account for irradiation effects 3. Develop a mechanistic macroscale model of densification for BISON

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Robust Dynamical Contribution to Precipitation Extremes in Idealized Warming Simulations across Model Resolutions

The impact of the circulation shift under climate warming on the distribution of precipitation extremes and the associated sensitivity to model resolution are investigated using the aquaplanet Community Atmosphere Model CAM3. The response of the probability density function of the precipitation to a uniform SST warming can be interpreted as superimposition of a dynamically induced poleward shift and a thermodynamically induced upward shift toward higher intensities, which give rise to manyfold increase in the frequency of the most extreme categories of the precipitation events at the poleward side of the midlatitude storm track. Meanwhile, the thermodynamic contribution to the intensification of the precipitation extremes is substantially less than expected from the Clausius-Clapeyron relation, implicative of significant change in the vertical structure of the precipitation processes. While coarser resolutions underestimate the dynamical contribution to the increase of precipitation extremes, a modest increase of the equator-to-pole SST warming gradient can have a significant opposite effect.

dynamical↗

Modeling assessment of storm surge in the Salish Sea

Here we present a modeling assessment of storm surge in the Salish Sea—a large and complex estuarine system in the Pacific Northwest that is exposed to extratropical cyclones. A method was developed to systematically identify storm surge events in the Salish Sea based on water level records at four representative tidal gauges in the Salish Sea. The maximum storm surge distribution in the Salish Sea was calculated based on model simulations of 34 major historical storm events from 1980 to 2016. The model simulations were validated by field observations of water level and velocity using a suite of model performance metrics. Model results suggest that the maximum storm surge in the Salish Sea is generally in the range of 0.8 m–1.03 m and is heterogeneous in space. The model results also show a strong north-south gradient with larger surge magnitudes occurring in the north that is primarily caused by the southerly wind in winter months. The sensitivity analysis of wind forcing suggested that storm surge in the Salish Sea is dominated by the remote surge from the Pacific Ocean, and the local wind contributes up to approximately 20% of the surge variability. This paper systematically elucidates storm surge characteristics and potential risk in the Salish Sea.

54 ENVIRONMENTAL SCIENCES↗

Ocean cavity regime shift reversed West Antarctic grounding line retreat in the late Holocene

Recent geologic and modeled evidence suggests that the grounding line of the Siple Coast of the West Antarctic Ice Sheet (WAIS) retreated hundreds of kilometers beyond its present position in the middle to late Holocene and readvanced within the past 1.7 ka. This grounding line reversal has been attributed to both changing rates of isostatic rebound and regional climate change. Here, we test these two hypotheses using a proxy-informed ensemble of ice sheet model simulations with varying ocean thermal forcing, global glacioisostatic adjustment (GIA) model simulations, and coupled ice sheet-GIA simulations that consider the interactions between these processes. Our results indicate that a warm to cold ocean cavity regime shift is the most likely cause of this grounding line reversal, but that GIA influences the rate of ice sheet response to oceanic changes. This implies that the grounding line here is sensitive to future changes in sub-ice shelf ocean circulation.

58 GEOSCIENCES↗

Hourly temperature and surface water level at the Ikpikpuk delta from 5/2012 to 11/2012

This data package includes the soil temperature and surface water level observations in 2012 at the Ikpikpuk River delta in northern Alaska. The data file consists of the observation table with both the .xlsx and .csv file formats. The files can be accessed through several software tools, like Microsoft Excel, and text editors. This is raw data from the near-surface ground temperature sensors at 15 cm (active layer) and 100 cm (permafrost) soil depth and a surface water sensor at the Ikpikpuk delta site in Alaska, US. The Ikpikpuk River delta is located on the Arctic Coastal Plain of northern Alaska. The Ikpikpuk River drains a land area of 15,330 km2, emptying into Smith Bay, and the delta is fan-shaped with an aerial extent of 106 km2 (Fuchs et al., 2018). This dataset supports the publication entitled "The thermal response of permafrost to coastal floodplain flooding" submitted to Environmental Research Letters (under review), where the observed water level is used as the upper boundary condition for the model simulations and the observed subsurface temperature is compared with the model simulations to answer the question about how coastal flooding affects permafrost thaw.

54 ENVIRONMENTAL SCIENCES↗

Space Station Freedom Data Assessment Study

The SSF Data Assessment Study was initiated to identify payload and operations data requirements to be supported in the Space Station era. To initiate the study payload requirements from the projected SSF user community were obtained utilizing an electronic questionnaire. The results of the questionnaire were incorporated in a personal computer compatible database used for mission scheduling and end-to-end communications analyses. This paper discusses data flow paths and associated latencies, communications bottlenecks, resource needs versus availability, payload scheduling 'warning flags' and payload data loading requirements for each major milestone in the Space Station buildup sequence. This paper also presents the statistical and analytical assessments produced using the data base, an experiment scheduling program, and a Space Station unique end-to-end simulation model. The modeling concepts and simulation methodologies presented in this paper provide a foundation for forecasting communication requirements and identifying modeling tools to be used in the SSF Tactical Operations Planning (TOP) process.

Johnson, Anngienetta R.↗

Combining Satellite Microwave Radiometer and Radar Observations to Estimate Atmospheric Latent Heating Profiles

In this study, satellite passive microwave sensor observations from the TRMM Microwave Imager (TMI) are utilized to make estimates of latent + eddy sensible heating rates (Q1-QR) in regions of precipitation. The TMI heating algorithm (TRAIN) is calibrated, or "trained" using relatively accurate estimates of heating based upon spaceborne Precipitation Radar (PR) observations collocated with the TMI observations over a one-month period. The heating estimation technique is based upon a previously described Bayesian methodology, but with improvements in supporting cloud-resolving model simulations, an adjustment of precipitation echo tops to compensate for model biases, and a separate scaling of convective and stratiform heating components that leads to an approximate balance between estimated vertically-integrated condensation and surface precipitation. Estimates of Q1-QR from TMI compare favorably with the PR training estimates and show only modest sensitivity to the cloud-resolving model simulations of heating used to construct the training data. Moreover, the net condensation in the corresponding annual mean satellite latent heating profile is within a few percent of the annual mean surface precipitation rate over the tropical and subtropical oceans where the algorithm is applied. Comparisons of Q1 produced by combining TMI Q1-QR with independently derived estimates of QR show reasonable agreement with rawinsonde-based analyses of Q1 from two field campaigns, although the satellite estimates exhibit heating profile structure with sharper and more intense heating peaks than the rawinsonde estimates. 2

Grecu, Mircea↗

Clumped-Isotope Constraint on Upper-Tropospheric Cooling During the Last Glacial Maximum

Ice cores and other paleotemperature proxies, together with general circulation models, have provided information on past surface temperatures and the atmosphere's composition in different climates. Little is known, however, about past temperatures at high altitudes, which play a crucial role in Earth's radiative energy budget. Paleoclimate records at high-altitude sites are sparse, and the few that are available show poor agreement with climate model predictions. These disagreements could be due to insufficient spatial coverage, spatiotemporal biases, or model physics; new records that can mitigate or avoid these uncertainties are needed. Here, we constrain the change in upper-tropospheric temperature at the global scale during the Last Glacial Maximum (LGM) using the clumped-isotope composition of molecular oxygen trapped in polar ice cores. Aided by global three-dimensional chemical transport modeling, we exploit the intrinsic temperature sensitivity of the clumped-isotope composition of atmospheric oxygen to infer that the upper troposphere (effective mean altitude 10–11 km) was 6–9°C cooler during the LGM than during the late preindustrial Holocene. A complementary energy balance approach supports a minor or negligible steepening of atmospheric lapse rates during the LGM, which is consistent with a range of climate model simulations. Proxy-model disagreements with other high-altitude records may stem from inaccuracies in regional hydroclimate simulation, possibly related to land-atmosphere feedbacks.

clumped isotope↗

Development of a Discrepancy Checker for the Digital Twin in a Supervisory Control System for a Thermal Energy Delivery System

Defined as a virtual representation of a physical object, process, or service, and used to support real-world decision-making, a digital twin (DT) can be utilized to combine classical and novel frameworks in sensors, state predictions, and multi-input/multi-output systems, and to enable optimal autonomous operations. However, a DT’s usefulness largely depends on its ability to adequately mirror the state of its physical counterpart, and this adequacy should be reflected by the level of uncertainty in the underlying simulation models when estimating and predicting quantities of interest (QOIs). Moreover, simulation models in a DT may involve multiple fidelities of representations—ranging from physics-based models to data-driven ones—but classical uncertainty quantification (UQ) methods struggle to handle numerous uncertainty sources, nor are they designed for real-time applications. This work presents a UQ-based discrepancy checking and diagnosis tool for a DT-based supervisory control system applied to a thermal energy delivery system (TEDS) at Idaho National Laboratory. The discrepancy checker was developed using metadata from an automated DT development process, and these metadata included different combinations of physical model forms and model parameters, training data and hyperparameters for surrogate models, and design parameters for supervisory control systems. Next, correlations between the uncertainty results and the metadata were established and then applied to the DT operations. The discrepancy checker evaluates the discrepancies between model predictions from virtual and sensor measurements and backtraces them to the corresponding major sources of uncertainty. The discrepancy checker showed reasonable performance in detecting discrepancies and diagnosing sources of uncertainty in testing scenarios.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

User's manual for a parameter identification technique

A digital computer program written in FORTRAN is presented that implements the system identification theory for deterministic systems using input-output measurements. The user supplies programs simulating the mathematical model of the physical plant whose parameters are to be identified. The user may choose any one of three options. The first option allows for a complete model simulation for fixed input forcing functions. The second option identifies up to 36 parameters of the model from wind tunnel or flight measurements. The third option performs a sensitivity analysis for up to 36 parameters. The use of each option is illustrated with an example using input-output measurements for a helicopter rotor tested in a wind tunnel.

Kanning, G.↗

A Composite View of Ozone Evolution in the 1995-96 Northern Winter Polar Vortex Developed from Airborne Lidar and Satellite Observations

The processes which contribute to the ozone evolution in the high latitude lower stratosphere are evaluated using a three dimensional model simulation and ozone observations. The model uses winds and temperatures from the Goddard Earth Observing System Data Assimilation System. The simulation results are compared with ozone observations from three platforms: the differential absorption lidar (DIAL) which was flown on the NASA DC-8 as part of the Vortex Ozone Transport Experiment; the Microwave Limb Sounder (MLS) on the Upper Atmosphere Research Satellite; and the Polar Ozone and Aerosol Measurement (POAM II) solar occulation instrument, on board the French Satellite Pour I'Observations de la Terre. Comparisons of the different data sets with the model simulation are shown to provide complementary information and a consistent view of the ozone evolution. The model ozone in December and January is shown to be sensitive to the ozone vertical gradient and the model vertical transport, and only weakly sensitive to the model photochemistry. The most consistent comparison between observed and modeled ozone evolution is found for a simulation where the vertical profiles between 12 and 20 km within the polar vortex closely match December DIAL observations. Diabatic trajectory calculations are used to estimate the uncertainty due to vertical advection quantitatively. The transport uncertainty is significant, and should be accounted for when comparing observations with model ozone. The model ozone evolution during December and January is broadly consistent with the observations when these transport uncertainties are taken into account.

Douglass, Anne R.↗

Modeling and Simulation of Tank Pressure Control using Zero-Boiloff Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗

High-Accuracy Simulations to Model Pyrometallurgical Processes in a Secondary Lead Reverberatory Furnace

The US manufacturing industry produces about 1.3 million tons of refined lead each year using secondary sources consisting mainly of lead batteries. ORNL is partnering with Gopher resource, the second largest lead recycling company in the United States, and GTI, to develop a high-fidelity CFD model of a directly fired, reverberatory-style, secondary lead furnace. These High Performance Computing (HPC) simulations are aimed to use first principles modeling for combustion and melting processes of the secondary lead feed while accounting for complex interphase interactions between the gas, solid charge (lead) material, slag, and metal phases. Through validation against operating plant data, this effort will enable significant improvements in design, operational parameters, and energy efficiency, thus improving productivity and refractory lifetime of secondary lead melting furnaces. Estimated savings/reduction of, at least, 1 trillion BTU, 1 million ton/year of greenhouse gas emissions, and $\$50$ million/year to the US lead industry can be expected. ORNL resources and expertise in high-performance computing and multicomponent, multiphase flows were utilized to realize this goal while advancing the understanding of the smelting and melting processes occurring within the furnace.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗