Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “models and simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modeling of RMS Current in CSI Filter Capacitor and Minimum Conduction Loss Operation of CSI-Fed PMSM Drives for Traction Applications

This paper presents a new minimum conduction loss (MCL) torque control algorithm for current source inverter- (CSI-) fed permanent-magnet synchronous machine drives that improves the overall machine drive efficiency by minimizing the combined conduction losses in the inverter and machine at each operating point. First, analytical models of conduction losses in CSI-fed motor drives are presented, and a closed-form expression for optimal d-axis stator current to achieve MCL operation is derived. An expression for the rms current in both wye- and delta-connected CSI output capacitors is derived. A detailed simulation model that emulates the operation of a 100 kW SiC CSI-fed integrated machine drive (IMD) has been developed based on experimental results and finite element analysis. This simulation model is used to evaluate the proposed MCL control algorithm applied to a CSI-fed IMD system for a battery-electric vehicle traction drive, and the predicted total drive system loss with MCL control is compared to predicted losses with maximum-torque-per-ampere (MTPA) and minimum dc-link current (MDCC) control. Results show that MCL control can achieve loss reductions compared to the other two control algorithms over a wide range of operating conditions, with significant loss reductions >20% in the medium-speed regime.

42 ENGINEERING↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2026

A Framework for Parametric and Predictive Uncertainty Quantification in the E3SM Land Model: Assessing Site and Observable Generalizability

Quantifying parametric uncertainty using observations from individual sites provides a critical foundation for Earth system modeling, serving as a necessary first step before scaling up to regional or global applications. This study introduces a novel computational framework designed to enhance model predictability by reducing parametric uncertainty and assessing site and observable generalizability using various observational constraints. The framework integrates five components: Model Simulation, Statistical Emulation, Global Sensitivity Analysis (GSA), Model Calibration, and Model Prediction. Using the E3SM land model, we simulated site-level land-atmosphere carbon and energy fluxes from 2003 to 2007 across five evergreen needleleaf FLUXNET sites, perturbing 26 vegetation-related model parameters. Gaussian process emulators were employed to expedite GSA and model calibration. Four critical parameters that strongly influence selected land-atmosphere fluxes were identified by GSA. Bayesian approaches were used to infer parameter probability distributions leveraging synthetic data and FLUXNET observations. The results reveal that posterior parameter distributions vary significantly across different sites and observables within the same plant functional type. Probabilistic predictions indicate that parameters calibrated at one site can enhance predictive accuracy at other sites, although site heterogeneity may sometimes outweigh parametric uncertainty. Additionally, the probabilistic predictions demonstrate that calibration for one variable can also improve predictability for other variables, thereby maximizing predictive capabilities with limited observations. This framework provides a powerful approach for reducing parametric uncertainty in Earth system models and deepening our understanding of carbon dynamics and energy cycles. Its adaptability makes it a valuable tool for broader applications in Earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Model and remote-sensing-guided experimental design and hypothesis generation for monitoring snow-soil–plant interactions

In this study, we develop a machine-learning (ML)-enabled strategy for selecting hillslope-scale ecohydrological monitoring sites within snow-dominated mountainous watersheds, with a particular focus on snow-soil–plant interactions. Data layers rely on spatial data layers from both remote sensing and hydrological model simulations. Specifically, a Landsat-based foresummer drought sensitivity index is used to define the dependency of the annual peak plant productivity on the Palmer drought severity index in the early growing season. Hydrological simulations provide the spatiotemporal dynamics of near-surface soil moisture and snow depth. In this framework, a regression analysis identifies the key hydrological variables relevant to the spatial heterogeneity of drought sensitivity. We then apply unsupervised clustering to these key variables, using the Gaussian mixture model, to group hillslopes into several zones that have divergent relationships regarding soil moisture, snow dynamics, and drought sensitivity. Using the datasets collected in the East River Watershed (Crested Butte, Colorado, United States), results show that drought sensitivity is significantly correlated with model-derived soil moisture and snow-free timing over space and time. The relationship is, however, non-linear, such that the correlation decreases above a threshold elevation and in a heavy snow year due to large snowpacks, lateral flow, and soil storage limitations. Clustering is then able to define the zones that have high or low sensitivity to drought, as well as the mid-elevation regions where sensitivity is associated with the topographic aspect and net potential radiation. In addition, the algorithm identifies the most representative hillslopes with road/trail access within each zone for installing monitoring sites. Our method also aims to significantly increase the use of ML and model-simulation results to guide critical zone and watershed monitoring activities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Design, modeling and simulation of nuclear-powered integrated energy systems with cascaded heating applications

Nuclear-renewable integrated energy systems (IES) consist of a variety of energy generation and conversion technologies and can be used to meet heterogeneous end uses (e.g., electricity, heat, and cooling demands). In addition to supply-demand balance, end-use heat demands usually require heat supply of certain temperature ranges. The effective and efficient utilization of heat produced within an IES is, therefore, a critical challenge. Here, this paper examines design options of an IES that includes heating processes of multiple temperature grades. We investigate a cascaded design configuration, where the remaining residual heat after high-grade heating processes [e.g., hydrogen production through high-temperature steam electrolysis (HTSE)] is recovered to meet the low-grade heating needs [e.g., district heating (DH)]. Additionally, a thermal energy storage system is integrated into the DH system to address the imbalance between heat supply and demand. This paper primarily focuses on the design and modeling of the proposed system and evaluates its operation with a 24-h transient process simulation using a DH demand profile with hourly resolution. The results indicate that the residual heat from the HTSE exhaust is insufficient for the DH demand, and additional topping heat directly from the reactor process steam is needed. Furthermore, the inclusion of thermal energy storage within the DH system provides the necessary balance between thermal generation and demand, thereby ensuring a consistent rated temperature of the DH supply water. This approach helps minimize the control actions needed on the reactor side.

08 HYDROGEN↗
DOE OSTI2023

Evaluation of enzymatic depolymerization of PET, PTT, and PBT polyesters

Millions of tons of waste polyester plastics, including polyethylene terephthalate (PET), polytrimethylene terephthalate (PTT), and polybutylene terephthalate (PBT), end up in the environment as soil and water contaminants. Recent advances in enzymatic polyester degradation have motivated researchers toward biorecycling of plastic wastes. Leaf-branch compost cutinase (LCC) enzymes have been proven to be effective in the biodegradation of PET. This study focuses on enzymatic depolymerization of waste PET, PTT, and PBT materials by using an ICCG variant of LCC (LCC ICCG ) produced from Escherichia coli BL21(DE3). The degradation efficiency of the polyesters was determined by the monomer terephthalic acid (TPA) released from the depolymerization reaction. It was found that the most efficient depolymerization was achieved for PET, followed by PTT and PBT. A kinetic model based on Langmuir adsorption and the Michaelis-Menten equation was developed to describe the enzymatic depolymerization of PET, PTT, and PBT with various enzyme and substrate loadings. The model simulation results revealed that the LCCICCG enzyme loading should be linearly increased as the work capacity of the polyester substrate increases. A specific enzyme loading of 0.91 mg/g PET is suggested to achieve 90% depolymerization of PET within three days. Furthermore, the experimental data and model simulation results can be used to help further engineer the enzyme and process to achieve a complete biodegradation of polyester wastes at large scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Elimination LArTPC Simulation Uncertainty

Liquid Argon Time Projection Chambers (LArTPC) are crucial for measuring muons and neutrinos by capturing the paths of fast-moving particles through argon gas. However, these detectors face challenges such as electron-ion recombination, diffusion, and attenuation, which introduce uncertainties in simulation models. This study, conducted by Ka ren Mkrtchyan at FERMILAB, aims to reduce these uncertainties by adjusting the amplitude and width of signals detected by the TPC wires. Initial findings indicate that the current modification algorithm requires further refinement to better align simulations with observed data. Ongoing work focuses on correcting computational bugs and enhancing the simulation model for improved accuracy and statistical confidence.

Mkrtchyan, Ka'ren↗
DOE OSTI2024

Coal-derived conductive pavement for winter de-icing: prototype, modeling, and simulation

Existing pavement de-icing methods result in high installation and maintenance costs, traffic delays, excessive weariness and corrosion, and negative environmental and safety impacts. To overcome these challenges, this paper presents an innovative pathway to designing and constructing smart self-heating pavements using a low-cost coal-char bearing asphalt material. The conductive asphalt incorporates coal char, i.e., a key byproduct of the coal pyrolysis process, into the Stone Mastic Asphalt (SMA) mixture. This asphalt containing coal-derived solid carbon exhibits highly tailorable electrical conductivity, satisfactory mechanical and thermophysical properties, and superior cost-efficiency as compared to other conductive pavement materials. The de-icing performance was also demonstrated by laboratory experiments on a bench-scale prototype. Furthermore, an efficient thermal network model was developed and validated by experiments to investigate the transient thermal behavior and energy performance of the Ohmic heating pavement system. Furthermore, the whole-year energy simulation case studies were conducted on a bridge pavement with an annual energy use of 24.6–1444.8 kWh/m 2 , showcasing its potential in field applications across cool humid, cold humid, and subarctic/arctic climate zones.

42 ENGINEERING↗
DOE OSTI2020

A generalized approximate control variate framework for multifidelity uncertainty quantification

We describe and analyze a variance reduction approach for Monte Carlo (MC) sampling that accelerates the estimation of statistics of computationally expensive simulation models using an ensemble of models with lower cost. These lower cost models — which are typically lower fidelity with unknown statistics — are used to reduce the variance in statistical estimators relative to a MC estimator with equivalent cost. We derive the conditions under which our proposed approximate control variate framework recovers existing multifidelity variance reduction schemes as special cases. We demonstrate that existing recursive/nested strategies are suboptimal because they use the additional low-fidelity models only to efficiently estimate the unknown mean of the first low-fidelity model. As a result, they cannot achieve variance reduction beyond that of a control variate estimator that uses a single low-fidelity model with known mean. However, there often exists about an order-of-magnitude gap between the maximum achievable variance reduction using all low-fidelity models and that achieved by a single low-fidelity model with known mean. We show that our proposed approach can exploit this gap to achieve greater variance reduction by using non-recursive sampling schemes. The proposed strategy reduces the total cost of accurately estimating statistics, especially in cases where only low-fidelity simulation models are accessible for additional evaluations. Several analytic examples and an example with a hyperbolic PDE describing elastic wave propagation in heterogeneous media are used to illustrate the main features of the methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

PyFLEXTRKR: a flexible feature tracking Python software for convective cloud analysis

Abstract. This paper describes the new open-source framework PyFLEXTRKR (Python FLEXible object TRacKeR), a flexible atmospheric feature tracking software package with specific capabilities to track convective clouds from a variety of observations and model simulations. This software can track any atmospheric 2D objects and handle merging and splitting explicitly. The package has a collection of multi-object identification algorithms, scalable parallelization options, and has been optimized for large datasets including global high-resolution data. We demonstrate applications of PyFLEXTRKR on tracking individual deep convective cells and mesoscale convective systems from observations and model simulations ranging from large-eddy resolving (∼100s m) to mesoscale (∼10s km) resolutions. Visualization, post-processing, and statistical analysis tools are included in the package. New Lagrangian analyses of convective clouds produced by PyFLEXTRKR applicable to a wide range of datasets and scales facilitate advanced model evaluation and development efforts as well as scientific discovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

FEM Modeling and Simulation of 2-D High Specific Heat Nb 3 Sn Wires

In the past few years, new high specific heat Nb 3 Sn wires have gained much focus at FNAL. Indeed, they have proved to be more stable against thermal perturbations with respect to standard wires. Nevertheless, a trade off exists between their thermal efficiency and production feasibility. In this report I describe the thermal and structural models that I have developed by exploiting ANSYS Mechanical APDL ® , which I got acquainted with at the beginning of my training. The aim has been to optimize the location of high specific heat elements in order to obtain an optimal thermal stability, while minimizing the risk of wire breakage during drawing. FEM results are compared with experimental ones made beforehand. Other minimum quench energy (MQE) thermal models have been developed for expected new experimental results. New data may improve the understanding of the embedded physical uncertainties in the model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Variable Pressure Multi-Pressure Rail System Design for Agricultural Applications

This paper presents a solution for reducing energy loss in the hydraulic control system of agricultural tractors and their implements. The solution is referred to as a multi-pressure rail (MPR) and provides power to the hydraulic functions following a pressure control logic, as opposed to the traditional flow control logic typical of hydraulic systems used in off-road vehicles. The proposed hydraulic control system allows for elimination of redundant flow control valves in the state-of-the-art system, which cause excessive throttling losses leading to poor overall energy efficiency. Related work on MPR technology targets construction vehicles, where the MPR solution can allow energy recovery during overrunning loads and better engine management. This paper alternatively addresses the case of agricultural applications where functions mostly operate under resistive load conditions with slow dynamics, which offers an opportunity to target throttle losses. For this purpose, the paper introduces a variable pressure control strategy to handle the instantaneous pressure at each rail. To develop both the controller and the hydraulic system architecture, a stationary test rig is conceived and used to validate a numerical simulation model of the MPR system and its control strategy. Particular focus is given to the dynamic behavior of the system during the switches of a function between different pressure rails, which needs to ensure reduced oscillations of the flow provided to each hydraulic function. Once validated, the simulation model is used to predict the energy savings of the MPR solution in an actual application: a 435 hp hydraulic tractor powering a 16-row planter, for which operating features during typical drive cycles were available to the authors. The results show up to 59% total power reduction at the pump shaft, corresponding to 89.8% system efficiency gain.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling and simulation of oxygen transport in high burnup LWR fuel

We have developed a formulation for oxygen transport in uranium dioxide nuclear fuel that accounts for the effects of irradiation. The overall simulation combines the evolving isotopic composition, thermochemistry, and oxygen transport in irradiated fuel. The driving forces for oxygen transport are computed from local thermodynamic equilibrium calculations and account for the effects of temperature gradients and composition, including fission products. The proposed method provides a mechanism for including complex thermodynamic models of nuclear fuel in modeling of mass redistribution, and alleviates difficulties associated with the common thermodiffusion formulation. Here, the transport model has been implemented within the nuclear fuel performance code BISON utilizing the thermochemistry code Thermochimica, with burnup calculations provided by the ORIGEN isotopic transmutation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Evaluating Pretreatment Strategies with Modeling for Reducing Scaling Potential of Reverse Osmosis Concentrate: Insights from Ion Exchange and Activated Alumina

Reverse osmosis concentrate (ROC) treatment is critical for enhancing water recovery and minimizing concentrate volume for disposal, especially in regions facing water scarcity. This study investigates the application of ion exchange (IX) resins and activated alumina (AA) as pretreatment strategies to mitigate scaling in ROC due to high concentrations of total dissolved solids, hardness (Ca 2+ and Mg 2+ ), and silica. Through a series of Langmuir isotherms, continuous column experiments, and model simulation, two types of strong acid cation IX resins and three types of strong base anion IX resins alongside three types of AA were evaluated. Results indicate that AA exhibits superior performance in silica removal, achieving up to a 65% reduction and maintaining performance for up to 800 bed volume without reaching saturation. Model simulation of a secondary reverse osmosis treating ROC after the IX and AA pretreatment indicated an additional water recovery of ~70% using antiscalants. This study demonstrates the potential for achieving higher water recovery while also identifying opportunities for pretreatment improvement. Challenges such as the limited IX capacity treating ROC, which requires frequent regeneration and increases operational costs, along with the restricted regeneration capacity of AA, underscore the importance of innovation. These findings emphasize the critical need for developing advanced materials and optimized strategies to further enhance the efficiency of ROC treatment processes.

activated alumina↗