Search NASA⌕ Search

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.

At least 811 records · Page 45

A Global Repository for Planet-Sized Experiments and Observations

Greenland’s ice loss doubled between 1996 and 2005. Nine of the ten warmest years on record have occurred since 2000. Carbon dioxide levels in the air are at their highest in 650,000 years. Are these phenomena related, and if so, what can they tell us about the Earth’s climate in the coming years and decades? The answer will require sophisticated computer models—and a way of dealing with the data they produce. Past climate patterns can be discerned by examining variables such as temperature and precipitation records, lake sediments, tree rings, and ice cores. Projecting future regional and global climate patterns, however, requires computer models that can simulate the complex interplay between physical, chemical, and biological processes occurring in the atmosphere, ocean, and on land and sea ice. These models are extremely computationally demanding. Over the past several decades, an ongoing revolution in high-capacity, high-performance computing has accelerated climate modeling. In 1998, a 1-year simulation of global climate using a 300- kilometer grid could be run in a day. In 2002, it took 4 or 5 hours. In 2015, this same model can be run on a high-end desktop in minutes.

54 ENVIRONMENTAL SCIENCES↗

The Role and Lifetime of Dissociative Heterogeneous Processes in Improving Simulated Ozone on Mars

Abstract Ozone simulated in Mars Global Climate Models (MGCMs) is used to assess the underlying chemistry occurring in the atmosphere. Currently, ozone total column abundance (TCA) is under‐predicted in MGCMs by up to 120%, implying missing or inaccurate chemistry in models. Heterogeneous reactions of hydroxyl radicals (HO X ) have been offered as an explanation for some of this bias, because they cause ozone to increase at locations where it's currently under‐predicted. We use four simulations to compare modeled ozone TCA with observations from the UVIS spectrometer aboard the ExoMars Trace Gas Orbiter to improve the representation of heterogeneous processes and their impact on ozone. We use a gas‐phase only run, a dissociative scheme, an adsorbed HO X retention scheme, and a hybrid scheme that combines the dissociative mechanism with the retention of HO X on water ice. We find retention of HO X is dependent on water ice sublimation, and ozone abundance increases when water ice persists for longer periods (1–20 sols). Over time, the loss of HO X causes a depletion in H 2 O 2 concentration (HO X reservoir), and thus allows ozone concentration to increase. When adsorbed HO X are desorbed and dissociate into other by‐products, HO X are not immediately available to destroy ozone. This results in larger ozone concentrations than if desorbed HO X are released directly back into their gaseous states. When using the hybrid scheme, ozone TCA is increased up to 50% where the ozone deficit is greatest, demonstrating the best agreement with observations, and implying that HO X radicals are both retained when adsorbed and dissociate.

Brown, M. A. J.↗

Material analysis of metal hydrides for bulk hydrogen storage

The United States Department of Energy (DOE) Hydrogen Storage Engineering Center of Excellence (HSECoE), 2009–2014, was established to evaluate the viability of material-based hydrogen storage for light duty vehicles. The HSECoE concluded that there were no metal hydrides that simultaneously met the DOE targets for gravimetric and volumetric capacity, and rate of onboard refueling. However, it was noted by the DOE that potential applications existed for stationary hydrogen storage systems. An application of particular interest is using a metal-hydride to store a hydrogen supply for fuel cell backup power to a 5 MWe data center. A generalized numerical model that coupled mass and energy conservation, chemical kinetics and thermodynamics was developed to simulate the discharge process for a given metal-hydride. The model was applied to an interstitial and a complex metal-hydride. In conclusion, it was found that the interstitial metal-hydride, which is discharged by a pressure swing is preferrable in a practical sense to the complex metal-hydride, which requires a temperature swing discharge process.

08 HYDROGEN↗

Changes in soil water content and lateral flow exert large effects on soil thermal dynamics across Alaskan landscapes

Both lateral surface and subsurface water flow affect soil moisture dynamics, yet most land surface models only solve subsurface water movement vertically. Here, we use a 3D ecosystem model that considers both land surface and subsurface hydrologic processes to simulate soil moisture, which is then used to drive a 1-D vertical soil thermal model to simulate the soil moisture effects on soil thermal dynamics in central Alaska. Our coupled model improves soil temperature (ST) estimates by 43.5% in comparison with observational data. Soil moisture has little effect on ST during the wet season (-1.5%) and a substantial influence during the dry season (60%). Spatially, water lateral flow has significant impacts on both soil moisture and ST, causing model estimates for thawed areas in the transition season to increase by ~10% in the study area. Our results highlight the importance of considering dynamical soil moisture, as well as lateral flow effects, on soil thermal dynamics in permafrost regions.

54 ENVIRONMENTAL SCIENCES↗

Design, tuning, and blackbox optimization of laser systems

Chirped pulse amplification (CPA) and subsequent nonlinear optical (NLO) systems constitute the backbone of myriad advancements in semiconductor manufacturing, communications, biology, defense, and beyond. Accurately and efficiently modeling CPA+NLO-based laser systems is challenging because of the complex coupled processes and diverse simulation frameworks. Our modular start-to-end model unlocks the potential for exciting new optimization and inverse design approaches reliant on data-driven machine learning methods, providing a means to create tailored CPA+NLO systems unattainable with current models. To demonstrate this new, to our knowledge, technical capability, we present a study on the LCLS-II photo-injector laser, representative of a high-power and spectro-temporally non-trivial CPA+NLO system.

47 OTHER INSTRUMENTATION↗

The need for carbon-emissions-driven climate projections in CMIP7

Abstract. Previous phases of the Coupled Model Intercomparison Project (CMIP) have primarily focused on simulations driven by atmospheric concentrations of greenhouse gases (GHGs), for both idealized model experiments and climate projections of different emissions scenarios. We argue that although this approach was practical to allow parallel development of Earth system model simulations and detailed socioeconomic futures, carbon cycle uncertainty as represented by diverse, process-resolving Earth system models (ESMs) is not manifested in the scenario outcomes, thus omitting a dominant source of uncertainty in meeting the Paris Agreement. Mitigation policy is defined in terms of human activity (including emissions), with strategies varying in their timing of net-zero emissions, the balance of mitigation effort between short-lived and long-lived climate forcers, their reliance on land use strategy, and the extent and timing of carbon removals. To explore the response to these drivers, ESMs need to explicitly represent complete cycles of major GHGs, including natural processes and anthropogenic influences. Carbon removal and sequestration strategies, which rely on proposed human management of natural systems, are currently calculated in integrated assessment models (IAMs) during scenario development with only the net carbon emissions passed to the ESM. However, proper accounting of the coupled system impacts of and feedback on such interventions requires explicit process representation in ESMs to build self-consistent physical representations of their potential effectiveness and risks under climate change. We propose that CMIP7 efforts prioritize simulations driven by CO2 emissions from fossil fuel use and projected deployment of carbon dioxide removal technologies, as well as land use and management, using the process resolution allowed by state-of-the-art ESMs to resolve carbon–climate feedbacks. Post-CMIP7 ambitions should aim to incorporate modeling of non-CO2 GHGs (in particular, sources and sinks of methane and nitrous oxide) and process-based representation of carbon removal options. These developments will allow three primary benefits: (1) resources to be allocated to policy-relevant climate projections and better real-time information related to the detectability and verification of emissions reductions and their relationship to expected near-term climate impacts, (2) scenario modeling of the range of possible future climate states including Earth system processes and feedbacks that are increasingly well-represented in ESMs, and (3) optimal utilization of the strengths of ESMs in the wider context of climate modeling infrastructure (which includes simple climate models, machine learning approaches and kilometer-scale climate models).

54 ENVIRONMENTAL SCIENCES↗

A generalized machine learning workflow to visualize mechanical discontinuity

Accurate detection and mapping of mechanical discontinuity in materials has widespread industrial and research applications. Herein, we developed a generalized machine-learning framework for visualizing single mechanical discontinuity embedded in material of any composition, velocity, density, porosity, and size with limited data. The proposed visualization of discontinuity requires accurate estimations of the length, location, and orientation of the embedded discontinuity by processing multipoint wave-transmission measurements. k-Wave simulator is used to create a large dataset of elastic waveforms recorded during multi-point wave-transmission measurements through materials containing single mechanical discontinuity. k-Wave simulator considers the wave attenuation, dispersion, and mode conversion in wave motion. Discrete wavelet transform (DWT) and statistical feature extraction are essential for data preprocessing prior to the data-driven model development. DWT also minimizes the effect of noise. Using hyper-parameter tuning and cross validation, gradient boosting regression can visualize the mechanical discontinuity with an accuracy of 0.85, in terms of coefficient of determination. A double-layered neural network-based regression has better performance with an accuracy of 0.95. Use of convolutional neural network converts the predictive task from a waveform processing to an image processing problem. Convolutional neural network achieved a generalization performance of 0.91. The proposed generalized workflow requires robust simulation of wave propagation, signal processing, feature engineering, and model evaluation. Sensors closest to the source and those located opposite the source are the most significant for the desired visualization. Notably, the sensors closest to the source capture the non-linear associations, whereas the sensor on the border opposite to the source capture the linear associations between the measured waveforms and the properties of the mechanical discontinuity.

42 ENGINEERING↗

CALIBRATED ULTRA FAST IMAGE SIMULATIONS FOR THE DARK ENERGY SURVEY

Image simulations are becoming increasingly important in understanding the measurement process of the shapes of galaxies for weak lensing and the associated systematic effects. For this purpose we present the first implementation of the Monte Carlo Control Loops (MCCL), a coherent framework for studying systematic effects in weak lensing. It allows us to model and calibrate the shear measurement process using image simulations from the Ultra Fast Image Generator (UFig) and the image analysis software SExtractor. We apply this framework to a subset of the data taken during the Science Verification period (SV) of the Dark Energy Survey (DES). We calibrate the UFig simulations to be statistically consistent with one of the SV images, which covers ∼0.5 square degrees. We then perform tolerance analyses by perturbing six simulation parameters and study their impact on the shear measurement at the one-point level. This allows us to determine the relative importance of different parameters. For spatially constant systematic errors and point-spread function, the calibration of the simulation reaches the weak lensing precision needed for the DES SV survey area. Furthermore, we find a sensitivity of the shear measurement to the intrinsic ellipticity distribution, and an interplay between the magnitude-size and the pixel value diagnostics in constraining the noise model. This work is the first application of the MCCL framework to data and shows how it can be used to methodically study the impact of systematics on the cosmic shear measurement.

79 ASTRONOMY AND ASTROPHYSICS↗

A Modularized Urban Scale Building Energy Modeling Framework Designed with An Open Mind

In recent years, physics-based building energy modeling (BEM) has started being used to evaluate the performance of buildings in the context of connected communities and on an urban scale to study their aggregated energy use, interactions, and impacts on the energy supply infrastructure and environment. The development of urban-scale BEM solutions needs extensive effort. Existing attempts tend to focus on different aspects of BEM on an urban scale, such as collecting as-built building data from different information sources, integrating geometry modeling with geographic information systems (GISs), representing operational and occupancy profiles, automating workflow, processing and visualizing the results, and conducting large-scale simulations. Urban-scale BEM development would benefit from multi-disciplinary research areas and from an open platform to adopt advancements on data sources and tools. For these purposes, this research proposes a modularized bottom-up model creation and simulation framework that is built on the state-of-the-art BEM tools and can accommodate different building stock data. This framework uses a standardized schema to describe building design and operational characteristics, and it can be instantiated from different building survey datasets with heterogeneous structures. The paper demonstrates how thousands of surveyed buildings from the 2012 U.S. Energy Information Administration’s Commercial Buildings Energy Consumption Survey (CBECS) were one-to-one converted to EnergyPlus models through the schema and the model generation process, then simulated with distributed computing, and their results are summarized.

Lei, Xuechen↗

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation: Quarterly Research Performance Progress Report, QR4 (Q4FY24)

We developed an integrated Computational Fluid Dynamics (CFD) model to simulate the multi-physics coupled cooling process of mixed gas by cold water within a Direct Contact Cooler (DCC) equipped with a rotating packing bed (RPB). The model captures the interactions between fluid dynamics, heat transfer, mass transport, and phase transitions, while accounting for key operational variables such as RPB rotational speed and the mass flow rates of both liquid and gas. The CFD model has been validated using experimental data, specifically by comparing predicted outflow gas and liquid temperatures to measured results. Our findings demonstrate the significant effects of RPB rotational speed and mass flow rates on cooling performance, providing valuable insights for optimizing DCC efficiency in industrial applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Integrating a Water Tracer Model Into WRF‐Hydro for Characterizing the Effect of Lateral Flow in Hydrologic Simulations

Abstract Most current land models approximate terrestrial hydrological processes as one‐dimensional vertical flow, neglecting lateral water movement from ridges to valleys. Such lateral flow is fundamental at catchment scales and becomes crucial for finer‐scale land models. To test the effect of incorporating lateral flow toward three‐dimensional representations of hydrological processes in the next generation land models, we integrate a water tracer model into the WRF‐Hydro framework to track water movement from precipitation to discharge and evapotranspiration. This hydrologic‐tracer integrated system allows us to identify the key mechanisms by which lateral flow affects the flow paths and transit times in WRF‐Hydro. By comparing modeling experiments with and without lateral routing in two contrasting catchments, we determine the impacts of lateral flow on the transit times of precipitation event‐water. Results show that with limited hydrologic connectivity, lateral flow extends the transit times by reducing (increasing) event‐water drainage loss (accumulation) in ridges (valleys) and allowing reinfiltration of infiltration‐excess flow, which is missing in most land models. On the contrary with high hydrologic connectivity, lateral flow can effectively accelerate the water release to streams and reduce the transit time. However, the transit times are substantially underestimated by the model compared with isotope‐derived estimates, indicating model limitations in representing flow paths and transit times. This study provides some insights on the fundamental differences in terrestrial hydrology simulated by land models with and without lateral flow representation.

54 ENVIRONMENTAL SCIENCES↗

MASKED AREAS IN SHEAR PEAK STATISTICS: A FORWARD MODELING APPROACH

The statistics of shear peaks have been shown to provide valuable cosmological information beyond the power spectrum, and will be an important constraint of models of cosmology in forthcoming astronomical surveys. Surveys include masked areas due to bright stars, bad pixels etc., which must be accounted for in producing constraints on cosmology from shear maps. We advocate a forward-modeling approach, where the impacts of masking and other survey artifacts are accounted for in the theoretical prediction of cosmological parameters, rather than correcting survey data to remove them. We use masks based on the Deep Lens Survey, and explore the impact of up to 37% of the survey area being masked on LSST and DES-scale surveys. By reconstructing maps of aperture mass the masking effect is smoothed out, resulting in up to 14% smaller statistical uncertainties compared to simply reducing the survey area by the masked area. We show that, even in the presence of large survey masks, the bias in cosmological parameter estimation produced in the forward-modeling process is ≈1%, dominated by bias caused by limited simulation volume. We also explore how this potential bias scales with survey area and evaluate how much small survey areas are impacted by the differences in cosmological structure in the data and simulated volumes, due to cosmic variance.

79 ASTRONOMY AND ASTROPHYSICS↗

1D model of tokamak scrape-off layer and divertor plasmas: Dimensionless parameters and the results of numerical simulations

Solutions of a 1D model of the scrape-off layer and a divertor plasma in a tokamak are analyzed on the basis of dimensionless parameters governing the physics of both plasma and neutral transport and atomic physics processes. Numerical simulations based on fluid plasma and Monte Carlo neutral modeling confirm earlier results obtained with simplified analytic considerations and numerical simulations performed with fluid neutral description and demonstrate bifurcation of plasma parameters. Further, these simulations also confirm that the major physical reasons for the rollover of the plasma flux to the target, which is the key feature of the phenomenon of divertor plasma detachment, are the impurity radiation loss and plasma recombination. It is suggested that simulations of the edge plasma based on dimensionless parameters provide a powerful tool for checking the consistency of the physics incorporated in existing complex 2D codes used for edge plasma modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Scattering Outcomes of Kepler Circumbinary Planets: Planet Mass Ratio

Recent studies reveal that the free eccentricities of Kepler-34b and Kepler-413b are much larger than their forced eccentricities, implying that scattering events may take place in their formation. The observed orbital configuration of Kepler-34b cannot be well reproduced in disk-driven migration models, whereas a two-planet scattering scenario can play a significant role of shaping the planetary configuration. These studies indicate that circumbinary planets discovered by Kepler may have experienced scattering process. In this work, we extensively investigate the scattering outcomes of circumbinary planets focusing on the effects of planet mass ratio . We find that the planetary mass ratio and the the initial relative locations of planets act as two important parameters that affect the eccentricity distribution of the surviving planets. As an application of our model, we discuss the observed orbital configurations of Kepler-34b and Kepler-413b. We first adopt the results from the disk-driven models as the initial conditions, then simulate the scattering process that occurs in the late evolution stage of circumbinary planets. We show that the present orbital configurations of Kepler-34b and Kepler-413b can be well reproduced when considering a two unequal-mass planet ejection model. Our work further suggests that some of the currently discovered circumbinary single-planet systems may be survivors of original multiple-planet systems. The disk-driven migration and scattering events occurring in the late stage both play an irreplaceable role in sculpting the final systems.

79 ASTRONOMY AND ASTROPHYSICS↗

A unified framework for polycrystal plasticity with grain boundary evolution

Plastic deformation in polycrystals is governed by the interplay between intra-granular slip and grain boundary-mediated plasticity. However, while the role played by bulk dislocations is relatively well-understood, the contribution of grain boundaries (GBs) has only recently begun to be studied. GB plasticity is known to play a key role along with bulk plasticity under a wide range of conditions, such as dynamic recovery, superplasticity, severe plastic deformation, etc., and developing models capable of simultaneously capturing GB and bulk plasticity has become a topic of high relevance. In this paper we develop a thermodynamically-consistent polycrystal plasticity model capable of simulating a variety of grain boundary-mediated plastic processes in conjunction with bulk dislocation slip. The model starts from the description of a single crystal and creates lattice strain-free polycrystalline configurations by using a specially-designed multiplicative decomposition developed by the authors. This leads to the introduction of a particular class of geometrically necessary dislocations (GND) that define fundamental GB features such as misorientation and inclination. The evolution of the system is based on an energy functional that uses a non-standard function of the GND tensor to account for the grain boundary energy, as well as for the standard elastic energy. Our implementation builds on smooth descriptions of GBs inspired on diffuse-interface models of grain evolution for numerical convenience. Here we demonstrate the generality and potential of the methodology by simulating a wide variety of phenomena such as shear-induced GB sliding, coupled GB motion, curvature-induced grain rotation and shrinkage, and polygonization via dislocation sub-grain formation.

36 MATERIALS SCIENCE↗

Time-of-Flight Estimation for Nondestructive Evaluation

Ultrasonic testing (UT) for nondestructive evaluation (NDE) is a critical entity necessary to resolve both the quality and precision questions of complex parts evolving from the classical approaches to the innovative additive manufacturing (AM) process. This modality provides the essential quantitative information for acceptance and potential flaw detectionof apart under investigation. A primary ingredient in UT besides the required precision robotic hardware for theacquisition of high-qualitymeasurement data is the underlying signal processing. It is here that much of the system performance capability resides. In this report,we discuss the basic steps in UT signal processing along with current capabilities that must be achieved in order to satisfy the criticaldemands of the various LLNL programs.In itsbasic foundations, signal processing is essentially the “extraction of critical information (signals) from noisy, uncertain data.” Clearly, in a perfect world the best signal processing is none---just make a reliable, noise-free, uncertainty-free, measurement. Unfortunately, even the best of systems is still straddled within the confines and limitations of physical instrumentation. With this in mind, we discuss the development of signal processing techniques toextract the desired UT information from noisy measurement data. We start with a discussion of a simple homogeneous representation of a “part” under investigation and its insonification by an ultrasonicexcitation from a physics-based perspective. Once developed, we briefly discussthe various choices of excitation signals and their tradeoffs.Next,we discuss a simulation approach employing simple models from the signal processing perspective and then move on to the development of the basic signal processing approach to ultrasonic signal processing all based on estimating the pulse arrival estimation. Starting with simple peak detection techniques, progressing to the processing workhorse so-called “matched-filter” and finally to the more sophisticated “model-based matched filter.”The underlying pre-and post-processing of noisy measurement data along with the application of these techniques are subsequently demonstrated on experimental ultrasonic data.Finally, we present some results for processing “weld” data when the ultrasonic excitation signal is not directly available.

42 ENGINEERING↗

Powder spreading, densification, and part deformation in binder jetting additive manufacturing

Binder jetting additive manufacturing (AM) can print complex structures in economical and scalable manner. Binder jetting AM comprises of deposition and weak binding of particles, known as green part, at room temperature and subsequent binder removal and sintering densification at high temperatures. However, during the densification (i.e., sintering), the part significantly deforms due to volume shrinkage. The deformation during sintering is difficult to predict, which prevents the widespread application of this technology. In this research, powder spreading simulation using discrete element method (DEM) was performed first to capture local variations in powder bed configuration. Second, finite element method (FEM) with a phenomenological continuum constitutive model was used to predict part shrinkage during the sintering process. DEM simulation showed variations in packing density, particle segregation, formation of uneven powder bed surface, and shift in particle size distribution (PSD). The sintering simulation modeled part deformation with a reasonable accuracy of 3% for solid-state sintering and intermediate liquid phase sintering. A demonstration case with non-uniform initial packing density showed that inhomogeneous green part density and PSD should be accounted for prediction of part deformation in binder jetting AM.

36 MATERIALS SCIENCE↗

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗