SEARCH · Search NASA
Results for “process simulation 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.
Collaborative Research: Process-resolving Decomposition of the Global Temperature Response to Modes of Low Frequency Variability in a Changing Climate
El Niño-Southern Oscillation (ENSO) and Annular Modes (AMs) represent respectively the most important modes of low frequency variability in the tropical and extratropical circulations. The future projection of the ENSO and AM variability, however, remains highly uncertain with the state-of-the-art coupled general circulation models. A comprehensive understanding of the factors responsible for the inter-model discrepancies in projecting future changes in the ENSO and AM variability, in terms of multiple feedback processes involved, has yet to be achieved. The proposed research aims to identify sources of such uncertainty and establish a set of process-resolving quantitative evaluations of the existing predictions of the future ENSO and AM variability. The proposed process-resolving evaluations are based on a feedback analysis method formulated in Lu and Cai (2009), which is capable of partitioning 3D temperature anomalies/perturbations into components linked to 1) radiation-related thermodynamic processes such as cloud and water vapor feedbacks, 2) local dynamical processes including convection and turbulent/diffusive energy transfer and 3) non-local dynamical processes such as the horizontal energy transport in the oceans and atmosphere. Taking advantage of the high-resolution, multi-model ensemble products from the Coupled Model Intercomparison Project Phase 5 (CMIP5) soon to be available at the Lawrence Livermore National Lab, we will conduct a process-resolving decomposition of the global three-dimensional (3D) temperature (including SST) response to the ENSO and AM variability in the preindustrial, historical and future climate simulated by these models. Specific research tasks include 1) identifying the model-observation discrepancies in the global temperature response to ENSO and AM variability and attributing such discrepancies to specific feedback processes, 2) delineating the influence of anthropogenic radiative forcing on the key feedback processes operating on ENSO and AM variability and quantifying their relative contributions to the changes in the temperature anomalies associated with different phases of ENSO and AMs, and 3) investigating the linkages between model feedback processes that lead to inter-model differences in time-mean temperature projection and model feedback processes that cause inter-model differences in the simulated ENSO and AM temperature response. Through a thorough model-observation and inter-model comparison of the multiple energetic processes associated with ENSO and AM variability, the proposed research serves to identify key uncertainties in model representation of ENSO and AM variability, and investigate how the model uncertainty in predicting time-mean response is related to the uncertainty in predicting response of the low-frequency modes. The proposal is thus a direct response to the first topical area of the solicitation: Interaction of Climate Change and Low Frequency Modes of Natural Climate Variability. It ultimately supports the accomplishment of the BER climate science activity Long Term Measure (LTM): "Deliver improved scientific data and models about the potential response of the Earth's climate and terrestrial biosphere to increased greenhouse gas levels for policy makers to determine safe levels of greenhouse gases in the atmosphere."
Prediction of grain structure after thermomechanical processing of U-10Mo alloy using sensitivity analysis and machine learning surrogate model
Abstract Hot rolling and annealing are critical intermediate steps for controlling microstructures and thickness variations when fabricating uranium alloyed with 10% molybdenum (U-10Mo), which is highly relevant to worldwide nuclear non-proliferation efforts. This work proposes a machine-learning surrogate model combined with sensitivity analysis to identify and predict U-10Mo microstructure development during thermomechanical processing. Over 200 simulations were collected using physics-based microstructure models covering a wide range of thermomechanical processing routes and initial alloy grain features. Based on the sensitivity analysis, we determined that an increase in rolling reduction percentage at each processing pass has the strongest effect in reducing the grain size. Multi-pass rolling and annealing can significantly improve recrystallization regardless of the reduction percentage. With a volume fraction below 2%, uranium carbide particles were found to have marginal effects on the average grain size and distribution. The proposed stratified stacking ensemble surrogate predicts the U-10Mo grain size with a mean square error four times smaller than a standard single deep neural network. At the same time, with a significant speedup (1000×) compared to the physics-based model, the machine learning surrogate shows good potential for U-10Mo fabrication process optimization.
Modeling Tool Development and Validation for Solar Industry Process Heat Using Particle Thermal Energy Storage
U.S. industry sectors used 26.2 quadrillion Btu and accounted for 33% of total energy consumption in 2021 according to the Energy Information Agency. Industrial process heat accounts for 70% of industrial energy use with application temperatures ranging from 60 degrees -1100 degrees C. Industry processes, heavily relying on fossil fuels of cheap coal or natural gas, differ widely in operating conditions and load requirements which makes them difficult to standardize and imposes great challenges in decarbonization. Industry processes require reliable energy supply and vary widely in temperature ranges. Storing energy from renewable sources is necessary to improve reliability and to mitigate renewable intermittency when replacing carbon fuel-based heat supplies to achieve energy savings and reduce emissions. To this end, we have developed a particle-based thermal energy storage (TES) technology using low-cost and highly stable silica sand as a storage medium. The economic and performance-based analysis is key for renewable energy sources to reliably supply industry process heat and ultimately displace fossil fuels for decarbonization. The diversified industrial processes need case-by-case analysis and design. Therefore, an adaptive modeling tool is key for renewable power with energy storage to meet industry demands. Thus, a modeling tool to simulate a solar industry process heat system using the particle TES has been developed using the object-oriented equation-based language Modelica and the commercial platform of Modelon Impact. The Modelica-based software tool provides a general simulation environment for the design of reliable solar energy sources integrated with TES for various industrial process applications at different temperatures for economic competence with fossil fuels such as coal and natural gases. It uses both customized and standard component modeling modules in Modelon libraries for the flexibility to be adapted to a specific energy demand application. The particle TES system establishes a uniform energy supply platform with an efficient heat exchanger and particle thermal energy reservoir integrated with renewable powers. The particle TES system can provide a wide temperature range and can have a large storage temperature difference that increases storage energy density; therefore, it can be an adaptable energy storage system integrated with renewable power to supply 24/7 heat for industry decarbonization.
A combined kinetic Monte Carlo and phase field approach to model thermally activated dislocation motion
In this work, a kinetic Monte Carlo (KMC) scheme is integrated into a phase field dislocation dynamics (PFDD) model to account for thermal activation. While the implementation is general and could be applied to many different thermally activated processes, we apply the model to simulate a dislocation transmitting through an interface. Two cases are analyzed: a dislocation transmitting through a strong interface, and a dislocation transmitting through a more realistic grain boundary in tungsten. In conclusion, results and trends for both cases are discussed in detail, and ultimately show that the combined KMC-PFDD algorithm is effective in capturing thermally activated dislocation processes.
Data-Driven Mathematical Modeling of the Effect of Particle Size Distribution on the Transitory Reaction Kinetics of Hot Metal Desulfurization
The aim of this work was to develop a prediction model for hot metal desulfurization. More specifically, the study aimed at finding a set of explanatory variables that are mandatory in prediction of the kinetics of the lime-based transitory desulfurization reaction and evolution of the sulfur content in the hot metal. The prediction models were built through multivariable analysis of process data and phenomena-based simulations. The model parameters for the suggested model types are identified by solving multivariable least-squares cost functions with suitable solution strategies. One conclusion we arrived at was that in order to accurately predict the rate of desulfurization, it is necessary to know the particle size distribution of the desulfurization reagent. It was also observed that a genetic algorithm can be successfully applied in numerical parameter identification of the proposed model type. It was found that even a very simplistic parameterized expression for the first-order rate constant provides more accurate prediction for the end content of sulfur compared to more complex models, if the data set applied for the modeling contains the adequate information.
Adiabatic and radiative cooling are both important causes of aerosol activation in simulated fog events in Europe
Aerosol–fog interactions affect the visibility in, and life cycle of, fog and are difficult to represent in weather and climate models. Here we explore processes that impact the simulation of fog droplet number concentrations (N d ) at sub-kilometer scale horizontal grid resolutions in the UK Met Office Unified Model. We modify the parameterization of aerosol activation to include droplet activation by radiative cooling in addition to adiabatic cooling and determine the relative importance of the two cooling mechanisms. We further test the sensitivity of simulated N d to: (a) interception of droplets by trees and buildings, (b) overestimation of updrafts in temperature inversions (which leads to artificially high N d values), and (c) potential mechanisms for droplet deactivation due to downward fluctuations in supersaturation. We evaluate our model against observations from the ParisFog and LANFEX field campaigns, building on evaluation described in the companion paper. Including radiative cooling in the activation mechanism improves how accurately we represent the liquid water path and the vertical structure of the fog in our LANFEX case study. However, with radiative cooling, the N d are overestimated for most of the ParisFog cases and for the LANFEX case. The time-averaged overestimate exceeds a factor of three (the normalized mean bias factor exceeds 2.0) in 4 out of 11 ParisFog cases. Our sensitivity studies demonstrate how these overestimates can be mitigated. Assuming the overestimate affects both radiative and adiabatic cooling, we find that although radiative cooling is more often the dominant source, both cooling sources can sometimes dominate activation.
One-dimensional biomass fast pyrolysis model with reaction kinetics integrated in an Aspen Plus Biorefinery Process Model
A biomass fast pyrolysis reactor model with detailed reaction kinetics and one-dimensional fluid dynamics was implemented in an equation-oriented modeling environment (Aspen Custom Modeler). Portions of this work were detailed in previous publications; further modifications have been made here to improve stability and reduce execution time of the model to make it compatible for use in large process flowsheets. The detailed reactor model was integrated into a larger process simulation in Aspen Plus and was stable for different feedstocks over a range of reactor temperatures. Sample results are presented that indicate general agreement with experimental results, but with higher gas losses caused by stripping of the bio-oil by the fluidizing gas in the simulated absorber/condenser. Lastly, this integrated modeling approach can be extended to other well-defined, predictive reactor models for fast pyrolysis, catalytic fast pyrolysis, as well as other processes.
Penetration through slots in cylindrical cavities with cavity modes overlapping with the first slot resonance
In this work, we analyze the coupling into a slotted cylindrical cavity operating at fundamental cavity modal frequencies overlapping with the slot’s first resonance frequency through an unmatched formulation that accounts for the slot’s absorption and radiation processes. The model is validated through full-wave simulations and experimental data. We then couple the unmatched formulation to a perturbation theory model to investigate an absorber within the cavity to reduce the interior field strength, also validated with full-wave simulations and experiments. These models are pivotal to understanding the physical processes involved in the electromagnetic penetration through slots, and may constitute design tools to mitigate electromagnetic interference effects within cavities.
Melt-Pool Temperature and Size Measurement During Direct Laser Sintering
Additive manufacturing has demonstrated the ability to fabricate complex geometries and components not possible with conventional casting and machining. In many cases, industry has demonstrated the ability to fabricate complex geometries with improved efficiency and performance. However, qualification and certification of processes is challenging, leaving companies to focus on certification of material though design allowable based approaches. This significantly reduces the business case for additive manufacturing. Therefore, real time monitoring of the melt pool can be used to detect the development of flaws, such as porosity or un-sintered powder and aid in the certification process. Characteristics of the melt pool in the Direct Laser Sintering (DLS) process is also of great interest to modelers who are developing simulation models needed to improve and perfect the DLS process. Such models could provide a means to rapidly develop the optimum processing parameters for new alloy powders and optimize processing parameters for specific part geometries. Stratonics’ ThermaViz system will be integrated with the Renishaw DLS system in order to demonstrate its ability to measure melt pool size, shape and temperature. These results will be compared with data from an existing IR camera to determine the best approach for the determination of these critical parameters.
Investigation of Isobaric Mixing as a Mechanism for Boundary‐Layer Cloud Formation
This study investigates the potential role of isobaric mixing in the formation of marine boundary layer clouds. Cloud formation theory emphasizes uplift and adiabatic cooling, but recent observations support the existence of small clouds forming at various altitudes, even below the lifting condensation level. Isobaric mixing of air with different thermodynamic properties can generate localized supersaturation. A Gaussian mixing model is employed to simulate this process, considering the correlation between temperature and water vapor. Cloud droplet size distributions from aircraft measurements show a persistent and prominent mode of small droplets at 9 m, and the size of this mode compares favorably with predictions from the model. The results suggests that isobaric mixing plausibly contributes to the formation of clouds, particularly those observed at multiple altitudes with narrow droplet size distributions. This finding highlights the importance of considering isobaric mixing processes in understanding and modeling cloud formation.
Sensitivity of a Subject-specific Ankle Sprain Simulation to Extrinsic Versus Intrinsic Biomechanical Factors
Ankle sprains are the most common musculoskeletal injury in sport and military activity, despite existing prophylactic strategies. The purpose of this report was to develop a probabilistic simulation of lateral ankle sprains during single-limb drop landing, towards accelerating innovation in ankle sprain prevention. A deterministic, subject-specific musculoskeletal model was extended with automation and probabilistic distributions on sprain-related biomechanical factors. Probabilistic simulations were generated using traditional Monte Carlo techniques and the advanced mean value method, a more computationally-efficient approach. Predicted distributions of peak ankle joint rotations, velocities, and moments borne by supporting passive structures agreed favorably with the deterministic model and with reports of real sprain biomechanics. Parameter sensitivities identified that predictions were most strongly influenced by drop height, subtalar joint posture at contact, invertor/evertor co-activation, and passive ankle stiffness. The advanced mean value method predicted confidence bounds comparable to a 1000-trial Monte Carlo simulation, and required only 14 model evaluations and 4-min processing time. The extended probabilistic simulation may be useful to virtually test new prophylactic strategies for ankle sprains, and is made available for open-source use (https://simtk.org/projects/sprain-sim).
Global patterns and climate drivers of water-use efficiency in terrestrial ecosystems deduced from satellite-based datasets and carbon cycle models
Our aim is to investigate how ecosystem water-use efficiency (WUE) varies spatially under different climate conditions, and how spatial variations in WUE differ from those of transpiration-based water-use efficiency (WUE t ) and transpiration-based inherent water-use efficiency (IWUE t ). LocationGlobal terrestrial ecosystems. We investigated spatial patterns of WUE using two datasets of gross primary productivity (GPP) and evapotranspiration (ET) and four biosphere model estimates of GPP and ET. Spatial relationships between WUE and climate variables were further explored through regression analyses. Global WUE estimated by two satellite-based datasets is 1.9 ± 0.1 and 1.8 ± 0.6g C m -2 mm -1 lower than the simulations from four process-based models (2.0 ± 0.3g C m -2 mm -1 ) but comparable within the uncertainty of both approaches. In both satellite-based datasets and process models, precipitation is more strongly associated with spatial gradients of WUE for temperate and tropical regions, but temperature dominates north of 50 degrees N. WUE also increases with increasing solar radiation at high latitudes. The values of WUE from datasets and process-based models are systematically higher in wet regions (with higher GPP) than in dry regions. WUE t shows a lower precipitation sensitivity than WUE, which is contrary to leaf- and plant-level observations. IWUE t , the product of WUE t and water vapour deficit, is found to be rather conservative with spatially increasing precipitation, in agreement with leaf- and plant-level measurements. In conclusion, WUE, WUE t and IWUE t produce different spatial relationships with climate variables. In dry ecosystems, water losses from evaporation from bare soil, uncorrelated with productivity, tend to make WUE lower than in wetter regions. Yet canopy conductance is intrinsically efficient in those ecosystems and maintains a higher IWUEt. This suggests that the responses of each component flux of evapotranspiration should be analysed separately when investigating regional gradients in WUE, its temporal variability and its trends.
Modeling Heat Loss through Piston and Effects of Thermal Boundary Coatings in Diesel Engine Simulations using Conjugate Heat Transfer models
Heat loss through wall boundaries play a dominant role in the overall performance and efficiency of internal combustion engines. Typical engine simulations use constant temperature wall boundary conditions. These boundary conditions cannot be estimated accurately from experiments due to the complexities involved with engine combustion. As a result they introduce a large uncertainty in engine simulations and serve as a tuning parameter. Modeling the process of heat transfer through the solid walls in an unsteady engine computational fluid dynamics (CFD) simulation can lead to the development of higher fidelity engine calculations. These models can be used to study the impact of heat loss on engine efficiency and explore new design methodologies that can reduce heat losses. In this work, a single cylinder diesel engine is modeled along with the solid piston coupled to the fluid domain. Conjugate heat transfer (CHT) modeling techniques were implemented to model heat losses for a full cycle of a Navistar diesel engine. This CFD model is then validated against experimental data available from thermocouples embedded inside the piston surface. The overall predictions from the model match closely with the experimental observations. The validated model is further used to explore the benefits of thermal barrier coatings (TBC) on piston bowls. The effect of TBC coatings were modeled as a thermal resistance in the heat transfer models. Full cycle 3D engine simulations provide quantitative insights into heat loss and thus calculate the efficiency gain by the use of TBC coatings. The work establishes a validated modeling framework for CHT modeling in reciprocating engine simulations.
Assessment of Buffer-IPyC Thermomechanical Debonding Behavior Using New Experiment Strength Data in BISON
TRIstructural ISOtropic (TRISO) single sub-millimeter-diameter nuclear fuel particle consists of a spherical fuel kernel surrounded by four coating layers. During reactor operation, the TRISO coating layers are subjected to irradiation-induced dimensional changes and the associated thermomechanical behavior of each layer. One of the observed behaviors is gap formation between the buffer and IPyC layer due to the porous buffer’s irradiation-induced shrinkage exceeding that of the IPyC layer. To predict TRISO thermomechanical performance, multi-physics models have been built that are continually updated and refined. The BISON code is a finite element-based nuclear fuel performance code that may be used for TRISO particle simulations. One of the recent additions to the BISON code is the ability to model the process of layer debonding. This paper will focus on the simulation results of the improved BISON debonding model that will utilize updated strengths measured from irradiated TRISO fuel particles. The new experimental strength data from micromechanical tests of irradiated TRISO fuel samples were exercised in the BISON simulations and compared to baseline strength data to assess their applicability in the models. This also includes updated buffer-IPyC bond strengths to simulate layer delamination. Based on current experimental observations it is noted that the buffer-IPyC separation occurs not exactly at the junction of these two layers, but more on the side of the buffer layer. This improved modeling approach using experimental strength data to characterize buffer-IPyC debonding and its potential subsequent cracking will be presented in the paper along with comparisons to available experimental observations.
Improved methodologies for Earth system modelling of atmospheric soluble iron and observation comparisons using the Mechanism of Intermediate complexity for Modelling Iron (MIMI v1.0)
Abstract. Herein, we present a description of the Mechanism of Intermediate complexity for Modelling Iron (MIMI v1.0). This iron processing module was developed for use within Earth system models and has been updated within a modal aerosol framework from the original implementation in a bulk aerosol model. MIMI simulates the emission and atmospheric processing of two main sources of iron in aerosol prior to deposition: mineral dust and combustion processes. Atmospheric dissolution of insoluble to soluble iron is parameterized by an acidic interstitial aerosol reaction and a separate in-cloud aerosol reaction scheme based on observations of enhanced aerosol iron solubility in the presence of oxalate. Updates include a more comprehensive treatment of combustion iron emissions, improvements to the iron dissolution scheme, and an improved physical dust mobilization scheme. An extensive dataset consisting predominantly of cruise-based observations was compiled to compare to the model. The annual mean modelled concentration of surface-level total iron compared well with observations but less so in the soluble fraction (iron solubility) for which observations are much more variable in space and time. Comparing model and observational data is sensitive to the definition of the average as well as the temporal and spatial range over which it is calculated. Through statistical analysis and examples, we show that a median or log-normal distribution is preferred when comparing with soluble iron observations. The iron solubility calculated at each model time step versus that calculated based on a ratio of the monthly mean values, which is routinely presented in aerosol studies and used in ocean biogeochemistry models, is on average globally one-third (34 %) higher. We redefined ocean deposition regions based on dominant iron emission sources and found that the daily variability in soluble iron simulated by MIMI was larger than that of previous model simulations. MIMI simulated a general increase in soluble iron deposition to Southern Hemisphere oceans by a factor of 2 to 4 compared with the previous version, which has implications for our understanding of the ocean biogeochemistry of these predominantly iron-limited ocean regions.
CFD Analysis of RLUOB Zone 1 HEPA Filter Plenum and Testing Manifolds
The purpose of this computational fluid dynamics (CFD) analysis is to ensure that Camfil Farr’s upstream and downstream injection and sampling manifolds can meet or exceed the requirements outlined in the ASME AG-1 1997 a(2000) Code on Nuclear Air and Gas Treatment for the testing of HEPA and adsorbent filters. This paper will present a numerical simulation of airflow in the Radiological Laboratory Utility Office Building (RLUOB) zone 1 HEPA filter plenum and testing manifolds using the commercial CFD software ANSYS FLUENT 2020R1. The CFD analysis focuses on the investigation of the air flow distribution and air-aerosol mixing uniformity. The evaluation was done for all steps of the modeling process: grid generation, physics setup, simulation, and post-processing. The mass flow rate in each section of the zone 1 injection and sampling manifolds is also reported.