Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Process Modeling”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

BOREAS RSS-8 BIOME-BGC Model Simulations at Tower Flux Sites in 1994

BIOME-BGC is a general ecosystem process model designed to simulate biogeochemical and hydrologic processes across multiple scales (Running and Hunt, 1993). In this investigation, BIOME-BGC was used to estimate daily water and carbon budgets for the BOREAS tower flux sites for 1994. Carbon variables estimated by the model include gross primary production (i.e., net photosynthesis), maintenance and heterotrophic respiration, net primary production, and net ecosystem carbon exchange. Hydrologic variables estimated by the model include snowcover, evaporation, transpiration, evapotranspiration, soil moisture, and outflow. The information provided by the investigation includes input initialization and model output files for various sites in tabular ASCII format.

Hall, Forrest G.↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2024

Techno-economic and life cycle analysis of synthetic natural gas production from low-carbon H 2 and point-source or atmospheric CO 2 in the United States

Synthetic natural gas (SNG) is of great interest in reducing fossil energy consumption while maintaining compatibility with existing NG infrastructure and end-use applications equipment. SNG can be produced using clean H 2 generated from renewable or nuclear energy and CO 2 captured from stationary sources or the atmosphere. In this study, we develop an engineering process model of SNG production using Aspen Plus® and production scales reported by the industry. We examine the levelized cost and life cycle greenhouse gas (GHG) emissions of SNG production under various CO 2 supply scenarios. Considering the higher cost of H 2 transportation compared with CO 2 transportation, we assume that CO 2 feedstock is transported via pipeline to the H 2 production location, which is collocated with the SNG plant. We also evaluate the cost of CO 2 captured from the atmosphere, assuming the direct air capture process can occur near the SNG facility. Depending on the CO 2 supply chain, the levelized cost of SNG is estimated to be in the range of $\$45–76$ per million British thermal units (MMBtu) on a higher heating value (HHV) basis. The SNG production cost may be reduced to $\$27–57$/MMBtu-HHV by applying a tax credit available in the United States for low-carbon H 2 production (45V). With a lower electricity price of 3¢/kWh for water electrolysis and accounting for a 45V tax credit, the SNG cost reaches parity with the cost of fossil NG. Depending on the CO 2 supply chain, SNG can reduce life cycle GHG emissions by 52–88% compared with fossil NG.

03 NATURAL GAS↗
DOE OSTI2019

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗
DOE OSTI2022

Physics-Informed and Data-Driven Prediction of Residual Stress in Three-Dimensional Machining

Efficient and reliable prediction of machining-induced residual stress (RS) is a key requirement for truly integrated computational materials engineering (ICME). Currently available process modeling approaches, including empirical, analytical, and numerical methodologies lack predictive power and require substantial calibration and validation data. Moreover, most model-based approaches consider only two-dimensional (2D) (i.e., orthogonal), cutting processes. Meanwhile, industrial processes such as milling, turning, and drilling are inherently three-dimensional (3D). The present work attempts to bridge the gap between 2D and 3D through careful consideration of the process physics, including geometric, kinematic, and size-effect constraints to realize robust prediction of how RS develops in 3D machining. Using a novel in-situ experimental technique and digital image correlation (DIC) to determine equivalent Hertzian contact widths, contact pressures, and friction coefficients, the proposed methodology leverages a discretized conversion algorithm that includes multi-pass shakedown effects. This paper presents a semi-analytical model to predict machining-induced RS in 3D turning operations, which are used representatively for 3D processes more generally. Rather than follow a ‘brute force’ 3D FEM approach or conduct countless experiments to train a purely data-driven machine learning algorithm, the proposed approach builds on previous 2D modeling work. Through careful consideration of the process physics, including complex geometry/kinematic considerations of 3D turning, the authors demonstrated an experimentally calibrated approach, as well as validation based on published RS data. Model predictions and previously published measurement data of RS depth profiles for turning of Inconel 718 were compared for a range of process parameters. Correlation between the proposed 3D model and validation data was found to be within the margin of experimental error for most conditions. The proposed model appears to capture the overall behavior of 3D RS depth profiles with acceptable accuracy, particularly the key metrics of near-surface stress, peak stress magnitude and location, as well as overall stress profile depth. This report presents a physics-informed, data-driven approach for efficient calibration of a 2D model for machining-induced RS through DIC analysis of in-situ characterized subsurface displacement fields.

42 ENGINEERING↗
NASA NTRS2021

Circadian Rhythms and Homeostatic Mechanisms for Sleep Regulation

This chapter examines circadian rhythms and homeostatic mechanisms for sleep regulation. It reviews the current evidence describing the two-process model of sleep regulation and how to assess disruption to either of these sleep drives. This chapter also reviews the role of the photic and non-photic resetting of the circadian rhythm and describes how some aspects of modern society can cause sleep and circadian disruption. Furthermore, this chapter describes how misalignment between the circadian rhythm and sleep homeostat, such as occurs during jet lag and shiftwork, can lead to sleep disruption. The short and long-term consequences of circadian misalignment are also reviewed.

Sleep↗
DOE OSTI2022

Spatial and temporal prediction of radiation dose rates near Fukushima Daiichi Nuclear Power Plant

In this paper, we have developed a methodology to estimate the spatiotemporal distribution of radiation air dose rates around the Fukushima Daiichi Nuclear Power Plant (FDNPP). In our exploratory data analysis, we found that (1) the temporal evolution of dose rates is composed of a log-linear decay trend and fluctuations of air dose rates that are spatially correlated among adjacent monitoring posts; and (2) the slope of the log-linear environmental decay trend can be represented as a function of the apparent initial dose rates, coordinate position, land-use type, and soil type. From these observations, we first estimated the log-linear decay trend at each location based on these predictors, using the random forest method. We then developed a modified Kalman filter coupled with a Gaussian process model to estimate the dose-rate time series at a given location and time. We applied this method to the Fukushima evacuation zone (as of March 2017), which included 17 monitoring post locations (with monitoring datasets collected between 2014 and 2018) and generated a time series of dose-rate maps. Our results show that this approach allows us to produce accurate spatial and temporal predictions of radiation dose-rate maps using limited spatiotemporal measurements.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2012

Development of the GPM Observatory Thermal Vacuum Test Model

A software-based thermal modeling process was documented for generating the thermal panel settings necessary to simulate worst-case on-orbit flight environments in an observatory-level thermal vacuum test setup. The method for creating such a thermal model involved four major steps: (1) determining the major thermal zones for test as indicated by the major dissipating components on the spacecraft, then mapping the major heat flows between these components; (2) finding the flight equivalent sink temperatures for these test thermal zones; (3) determining the thermal test ground support equipment (GSE) design and initial thermal panel settings based on the equivalent sink temperatures; and (4) adjusting the panel settings in the test model to match heat flows and temperatures with the flight model. The observatory test thermal model developed from this process allows quick predictions of the performance of the thermal vacuum test design. In this work, the method described above was applied to the Global Precipitation Measurement (GPM) core observatory spacecraft, a joint project between NASA and the Japanese Aerospace Exploration Agency (JAXA) which is currently being integrated at NASA Goddard Space Flight Center for launch in Early 2014. From preliminary results, the thermal test model generated from this process shows that the heat flows and temperatures match fairly well with the flight thermal model, indicating that the test model can simulate fairly accurately the conditions on-orbit. However, further analysis is needed to determine the best test configuration possible to validate the GPM thermal design before the start of environmental testing later this year. Also, while this analysis method has been applied solely to GPM, it should be emphasized that the same process can be applied to any mission to develop an effective test setup and panel settings which accurately simulate on-orbit thermal environments.

Yang, Kan↗
DOE OSTI2023

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗
DOE OSTI2025

Technical, Economic, Energetic, and Environmental Evaluation of Pretreatment Strategies for Scaling Control in Brackish Water Desalination Brine Treatment

Effective pretreatment is essential for achieving long-term stable operation and high water recovery during the desalination of alternative waters. This study developed a process modeling approach for technical, economic, energetic, and environmental assessments of pretreatment technologies to identify the impacts of each technology treating brackish water desalination brine with high scaling propensity. The model simulations evaluated individual pretreatment technologies, including chemical softening (CS), chemical coagulation (CC), electrocoagulation (EC), and ion exchange (IX). In addition, combinations of these pretreatment technologies aiming at the effective reduction of key scaling constituents such as hardness and silica were investigated. The three evaluation parameters in this assessment consist of levelized cost of water (LCOW, $/m 3 ), specific energy consumption and cumulative energy demand (SEC|CED, kWh/m 3 ), and carbon dioxide emissions (CO 2 , kg CO 2-eq /m 3 ). The case study evaluated in this work was the desalination brine from the Kay Bailey Hutchison Desalination Plant (KBHDP) with a total dissolved solids (TDS) concentration of 11,000 mg/L and rich in hardness and silica. The evaluation of individual pretreatment units from the highest to lowest LCOW, SEC|CED, and CO 2 emissions in the KBHDP brine was IX > CS > EC > CC, CS > IX > EC > CC, and CC > CS > EC > IX, respectively. In the case of pretreatment combinations for the KBHDP, the EC + IX treatment combination was shown to be the best in terms of the LCOW and CO 2 emissions. The modeling and evaluation of these pretreatment units provide valuable guidance on the selection of cost-effective, energy-efficient, and environmentally sustainable pretreatment technologies tailored to desalination brine applications for minimal- or zero-liquid discharge.

14 SOLAR ENERGY↗
NASA NTRS2013

The Importance of Earth as a Meteor Detector

In order to understand the complex meteoroid environment that threatens spacecraft in near Earth space, ground based observational techniques provide us with an effective option for evaluating the important parameters of flux, speed, mass, density and directionality. While no observational method (in ]situ or ground based detectors) provides bias free measurements of direct quantities, there are several benefits to using the Earth as a detector for calibrating meteoroid models that extend to the inner Solar System and beyond Earth fs orbit. The Earth constantly interacts with meteoritic particles from comets and asteroids as they evolve from their parent objects through dynamical processes. Modeling these dynamical processes with a physics ]based approach can help reveal the causes of observational features but developing those models involves constraining sequences of poorly defined parameters. These parameters must be tuned to a particular detector that offers direct or indirect measurements of mass, density, flux and speed - the quantities used in evaluating penetration risk. The following topics will be discussed: a) Consistency between ground based and in.situ measurements of flux in Earth orbit. b) Consistency between radar and electro.optical speed determinations and agreement of these with dynamical studies. c) The fact that there are six sporadic sources which cannot be determined by integrated line of sight emissivity: Observations using IRAS or COBE cannot give you all important directionality needed for risk evaluations. d) The large collecting area of the Earth s atmosphere compared to in.situ detectors gives better number statistics in the threat regime. The NASA Meteoroid Environment Office utilizes ground based radar and electro.optical measurements to calibrate the Meteoroid Engineering Model (MEM). As MEM is an engineering design model our main emphasis is on threat size particle, speed, flux and directionality. However, MEM has acknowledged limitations particularly with regard to source strength as a function of mass and work is proceeding with improved distributions that will rectify many of these shortcomings.

Koehler, Heather↗
NASA NTRS2008

Cloud-System Resolving Models: Status and Prospects

Cloud-system resolving models (CRM), which are based on the nonhydrostatic equations of motion and typically have a grid-spacing of about a kilometer, originated as cloud-process models in the 1970s. This paper reviews the status and prospects of CRMs across a wide range of issues, such as microphysics and precipitation; interaction between clouds and radiation; and the effects of boundary-layer and surface-processes on cloud systems. Since CRMs resolve organized convection, tropical waves and the large-scale circulation, there is the prospect for several advances in both basic knowledge of scale-interaction requisite to parameterizing mesoscale processes in climate models. In superparameterization, CRMs represent convection, explicitly replacing many of the assumptions necessary in contemporary parameterization. Global CRMs have been run on an experimental basis, giving prospect to a new generation of climate weather prediction in a decade, and climate models due course. CRMs play a major role in the retrieval of surface-rain and latent heating from satellite measurements. Finally, enormous wide dynamic ranges of CRM simulations present new challenges for model validation against observations.

Tao, Wei-Kuo↗
DOE OSTI2020

The EGS Collab Project: Learnings from Experiment 1

The primary objective of the EGS Collab Project sponsored by DOE is to increase the understanding needed to efficiently implement enhanced geothermal systems (EGS). One goal of the EGS Collab project is to create a collaborative research environment in which to study stimulation of crystalline rock at the 10 meter scale. Key to this effort is the collection of high quality data to allow comparison to numerical coupled process models in an effort to build confidence in the codes and modeling techniques used. In response to this, the EGS Collab team has created an underground test bed at the Sanford Underground Research Facility (SURF) in Lead SD at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently designing a second test bed aimed at investigating shear stimulation (Experiment 2). At the Experiment 1 location, we have characterized our host rock using laboratory testing and numerous field-based geophysical and geological techniques, and created a well-instrumented test bed to allow us to carefully monitor stimulation events and flow tests. In addition to the installed geophysical sensors, we have used tracer tests, differences in the ambient microbial communities at flow collection locations, and cold water injection to inform us about dynamic flow pathways. In Experiment 1, we have hydraulically stimulated the host rock in a number times at several locations in one well, creating new fractures that connect to existing fractures between the injection and production boreholes. We have performed long-term ambient and chilled water injection tests as an analog to EGS, and have monitored system changes resulting from these water injections through geophysical monitoring, flow and pressure measurements, tracer tests, and microbiology. Here, we summarize the tests performed, issues identified including poroelastic and thermoelastic effects, Joule-Thomson effects, restarting effects, indications of flow channeling, and the primary learnings from Experiment 1.

Enhanced Geothermal Systems, EGS Collab, stimulati↗
DOE OSTI2025

Trusted Simulation: Considering Model Quality in the Context of User Trust

A high‐quality simulation model should help its users to easily and appropriately calibrate their trust in the model. Traditional evaluation metrics such as validation and robustness are necessary but insufficient for this task. Trust calibration depends on factors like the model's transparency, applicability to intended use, usability, reputation, and consideration of potential bias. This article proposes a framework for designing and evaluating system dynamics models by considering factors that contribute to the proper calibration of user trust. This framework takes inspiration from trusted artificial intelligence, broadening our traditional concept of model quality and explicitly focusing on what users need to consider a model trustworthy and to understand the model's relevance to its intended purpose. The trusted simulation framework can improve our integration of model quality activities throughout the modeling process, leading to more impactful and better‐targeted model design, development, and evaluation.

Naugle, Asmeret Bier [Sandia National Laboratories↗
DOE OSTI2021

Regional trends and drivers of the global methane budget

Abstract The ongoing development of the Global Carbon Project (GCP) global methane (CH 4 ) budget shows a continuation of increasing CH 4 emissions and CH 4 accumulation in the atmosphere during 2000–2017. Here, we decompose the global budget into 19 regions (18 land and 1 oceanic) and five key source sectors to spatially attribute the observed global trends. A comparison of top‐down (TD) (atmospheric and transport model‐based) and bottom‐up (BU) (inventory‐ and process model‐based) CH 4 emission estimates demonstrates robust temporal trends with CH 4 emissions increasing in 16 of the 19 regions. Five regions—China, Southeast Asia, USA, South Asia, and Brazil—account for >40% of the global total emissions (their anthropogenic and natural sources together totaling >270 Tg CH 4 yr −1 in 2008–2017). Two of these regions, China and South Asia, emit predominantly anthropogenic emissions (>75%) and together emit more than 25% of global anthropogenic emissions. China and the Middle East show the largest increases in total emission rates over the 2000 to 2017 period with regional emissions increasing by >20%. In contrast, Europe and Korea and Japan show a steady decline in CH 4 emission rates, with total emissions decreasing by ~10% between 2000 and 2017. Coal mining, waste (predominantly solid waste disposal) and livestock (especially enteric fermentation) are dominant drivers of observed emissions increases while declines appear driven by a combination of waste and fossil emission reductions. As such, together these sectors present the greatest risks of further increasing the atmospheric CH 4 burden and the greatest opportunities for greenhouse gas abatement.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

HERON

HERON is a modeling toolset and plugin for RAVEN to accelerate stochastic technoeconomic assessment of the economic viability of various grid-energy system configurations, especially with application to electrical grids and integrated energy systems (IES). The frameworks enabling this toolset are INL's RAVEN and its TEAL plugin. The toolset includes automated RAVEN workflow generation tools for performing stochastic analysis and optimization driven by technoeconomic factors. The workflows generated employ RAVEN for parametric or optimizing study of grid component characteristics, and again for stochastic sampling of regulated and deregulated markets, load, and weather histories. The toolset further provides methods for economically-optimal dispatch of system components given the stochastic samples of market, load, and weather histories. The toolset allows assembling of existing and new models through the RAVEN interface to the Modelica process modeling language.

Talbot, PaulW↗
NASA NTRS1988

IDSE Version 1 User's Manual

The integrated development support environment (IDSE) is a suite of integrated software tools that provide intelligent support for information modelling. These tools assist in function, information, and process modeling. Additional tools exist to assist in gathering and analyzing information to be modeled. This is a user's guide to application of the IDSE. Sections covering the requirements and design of each of the tools are presented. There are currently three integrated computer aided manufacturing definition (IDEF) modeling methodologies: IDEF0, IDEF1, and IDEF2. Also, four appendices exist to describe hardware and software requirements, installation procedures, and basic hardware usage.

Mayer, Richard↗