Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 757 records · Page 42

Coupling a Computational Fluid Dynamics Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop® (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

Heat transfer↗

Coupling a Computational Fluid Dynamics (CFD) Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

heat transfer↗

Need A Boost? A Comparison of Traditional Commuting Models with the XGBoost Model for Predicting Commuting Flows

Commuting models estimate the number of commuting trips from home to work locations in a given area. Since their infancy, they have been increasingly used in a variety of fields to reduce traffic and pollution, drive infrastructure choices, and solve a variety of other problems. Traditional commuting models, such as gravity and radiation models, typically have a strict structural form and limited number of input variables, which may limit their ability to predict commuting flows as well as machine learning models that might better capture the complex dynamics of the commuting process. To determine whether machine learning models might add value to the field of commuter flow prediction, we compare and discuss the performance of two standard traditional models with the XGBoost machine learning algorithm for predicting home to work commuter flows from a well-known United States commuting dataset. We find that the XGBoost model outperforms the traditional models on three commonly used metrics, indicating that machine learning models may add value to the field of commuter flow prediction.

Morton, April↗

Comparison of Surrogate Models with Physical Models for Dynamic Probabilistic Risk Analysis Using the RAVEN Code

Dynamic probabilistic risk analysis (DPRA) has emerged in the nuclear industry as an effective tool to identify points of weakness in nuclear reactors and create a comprehensive approach to quantifying risk. While DPRA allows for detailed analysis, it usually requires thousands of runs due to the necessity of exploring low probability, high consequence events, which can be computationally prohibitively difficult to assess. Risk Analysis Virtual Environment (RAVEN) is a software package being developed under the Nuclear Energy Advanced Modeling and Simulation program at the Idaho National Laboratory. It is able to manage complex control logic to drive a RELAP-7 simulation and to handle sampling strategies including Monte Carlo, Latin Hypercube, Grid and Adaptive samplers. By performing a thorough investigation of the possibility space of an accident scenario with RELAP-7, RAVEN is able to identify consequences of accident scenarios that have not been encountered. The most computational resource intensive step of a RAVEN analysis is running the RELAP-7 code modeling the physical system behavior. Using surrogate models to represent the outcomes of more complex simulations rather than the full model can greatly reduce the computation effort required. In situations where the uncertainty distributions of a model are unknown, or when there are multiple sets of uncertainty distributions, DPRA generally requires rerunning the physical model for each set of uncertainty distributions used. By taking an unbiased sample of the physical model and creating a surrogate model from the results, the surrogate model can be sampled under a different set of uncertainty distributions to provide equivalent results. This paper presents an example of this approach. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Drivers of Change in Managed Water Resources: Modeling the Impacts of Climate and Socioeconomic Changes Using the US Midwest as a Case Study

A global integrated assessment model including a water-demand model driven by socio-economics, is coupled in a one-way fashion with a land surface hydrology – routing – water resources management model. The integrated modeling framework is applied to the U.S. Upper Midwest (Missouri, Upper Mississippi, and Ohio) to advance understanding of the regional impacts of climate and socio-economic changes on integrated water resources. Implications for future flow regulation, water supply, and supply deficit are investigated using climate change projections with the B1 and A2 emission scenarios, which affect both natural flow and water demand. Changes in water demand are driven by socio-economic factors, energy and food demands, global markets and prices. The framework identifies the multiple spatial scales of interactions between the drivers of changes (natural flow and water demand) and the managed water resources (regulated flow, supply and supply deficit). The contribution of the different drivers of change are quantified regionally, and also evaluated locally, using covariances. The integrated framework shows that water supply deficit is more predictable over the Missouri than the other regions in the Midwest. The predictability of the supply deficit mostly comes from long term changes in water demand although changes in runoff has a greater contribution, comparable to the contribution of changes in demand, over shorter time periods. The integrated framework also shows that spatially, water demand drives local supply deficit. Using elasticity, the sensitivity of supply deficit to drivers of change is established. The supply deficit is found to be more sensitive to changes in runoff than to changes in demand regionally. It contrasts with the covariance analysis that shows that water demand is the dominant driver of supply deficit over the analysed periods. The elasticity indicates the level of mitigation needed to control the demand in order to reduce the vulnerability of the integrated system in future periods. The elasticity analyses also emphasize the need to address uncertainty with respect to changes in natural flow in integrated assessment.

climate change, water resources management, socio-↗

Local composition activity coefficient model for mixed-gas adsorption equilibria

Taking into consideration the adsorbate–adsorbent interactions, a novel activity coefficient model is derived from the non-random two-liquid theory for mixed-gas adsorption equilibria. In contrast with the conventional activity coefficient models developed for bulk liquids, the new model correctly predicts negative deviations from ideality for adsorbed phase mixtures including azeotropic behavior exhibited by selected gas-adsorbent systems. Here, the requiring a single binary interaction parameter per adsorbate–adsorbate pair, the model successfully correlates wide varieties of binary adsorption isotherm data and it should be a powerful engineering thermodynamic tool in correlating and predicting mixed-gas adsorption equilibria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring the use of dynamic linear panel data models for evaluating energy/economy/environment models — an application for the transportation sector

This paper evaluates the RoSE transportation sector scenarios of the GCAM and REMIND models for the U.S. region by comparing these models’ projected future final energy use levels (1) to historical trends and (2) across each other. To estimate and compare the intrinsic model elasticities, we explore the use of dynamic linear panel data models. On the basis of 26 scenarios (panels) from each of the two energy system models, our analysis suggests that most intrinsic model elasticities with respect to final energy use are roughly comparable to each other and to those observed historically. The key difference is the income elasticity of final energy use, which is significantly smaller for the two energy models than that derived from the historical development. This difference and other minor inconsistencies are interpreted through key assumptions underlying both energy system models.

54 ENVIRONMENTAL SCIENCES↗

Capturing the variable reactivity of goethites in surface complexation modeling by correlating model parameters with specific surface area

The development of predictive tools capable of simulating adsorption behavior on goethite (α-FeOOH) surfaces has evolved over the past several decades. Recent research has focused on capturing the variation of surface reactivity and adsorption behavior for different preparations of synthetic goethites. Here in this research, an approach to estimate modeling parameters based on the specific surface area (SSA) of goethites was established and combined with a Charge Distribution – Multi-site Complexation (CD-MUSIC) surface complexation model. It was assumed that differences in crystal face contributions (CFC), inner-Helmholtz capacitances (C 1 ), and protonation constants of surface oxygens (pK a ’s) among different preparations of goethite contribute to the variation of surface reactivity. Thus, correlations of these modeling parameters with the SSA of a range of goethites were identified by comprehensively combining microscopic, macroscopic, and modeling data reported in literature along with adsorption data collected from macroscopic experiments conducted as part of this study. A CD-MUSIC model shown to be capable of accurately predicting adsorption for several metal ions and oxyanion selenite on a 62.9 m 2 /g SSA goethite served as the basis for the development of the adsorption modeling approach described in this paper. The predictive capability of the newly developed modeling approach was verified by comparing predictions of protons, SeO 3 2- , and Cd 2+ adsorption with experimental data from goethites with SSAs ranging between 21.3 and 105 m 2 /g; in total, 40 adsorption data sets from 13 different goethites were used for model verification. The results from this study indicate that the variation of surface reactivity of goethites can be successfully simulated by changing CFC, C 1 , and pK 2,singly (i.e. second protonation constant of singly-coordinated oxygens on the goethite surface) based on the SSA, which implies that self-consistent adsorption parameters can be applied to predict adsorption behavior over a range of goethite preparations.

58 GEOSCIENCES↗

A review of predictive nonlinear theories for multiscale modeling of heterogeneous materials

Since the beginning of the industrial age, material performance and design have been in the midst of innovation of many disruptive technologies. Today's electronics, space, medical, transportation, and other industries are enriched by development, design and deployment of composite, heterogeneous and multifunctional materials. As a result, materials innovation is now considerably outpaced by other aspects from component design to product cycle. In this article, we review predictive nonlinear theories for multiscale modeling of heterogeneous materials. Deeper attention is given to multiscale modeling in space and to computational homogenization in addressing challenging materials science questions. Moreover, we discuss a state-of-the-art platform in predictive image-based, multiscale modeling with co-designed simulations and experiments that executes on the world's largest supercomputers. Such a modeling framework consists of experimental tools, computational methods, and digital data strategies. Once fully completed, this collaborative and interdisciplinary framework can be the basis of Virtual Materials Testing standards and aids in the development of new material formulations. Moreover, it will decrease the time to market of innovative products.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying the Influence of Size, Shape, and Density of Microplastics on Their Transport Modes: A Modeling Approach

Microplastics (MPs) pose significant risks to marine ecosystems and human health, necessitating accurate predictions of their distributions in aquatic environments for effective risk mitigation. However, understanding MP transport dynamics is challenging because of the inadequate representation of MP characteristics such as size, shape, and density in numerical models. Further, the accuracy of the MP vertical profiles in existing models has not been thoroughly validated. Thus, we developed an MP transport model within the Finite Volume Community Ocean Model framework (FVCOM-MP) by integrating MP characteristics. We validated FVCOM-MP against experimental and analytical data, focusing on various MP transport modes and transitions. FVCOM-MP successfully replicates MP profiles in different transport modes, including the bedload, surface-load, suspended-load, and mixed-load modes. Additionally, we introduce phase diagrams for classifying MP transport modes based on particle characteristics, enhancing our understanding of MP dynamics in aquatic systems. The transport modes for a number of real-world MP particles, including fishing line, plastic bag/bottle fragments, synthetic fibers, tire wear particles, polyvinyl chloride and expanded polystyrene pellets, were analyzed with these phase diagrams.

Microplastic transport, Settling velocity, Rising ↗

A robust Bayesian approach to modeling epistemic uncertainty in common-cause failure models

In a standard Bayesian approach to the alpha-factor model for common-cause failure, a precise Dirichlet prior distribution models epistemic uncertainty in the alpha-factors. This Dirichlet prior is then updated with observed data to obtain a posterior distribution, which forms the basis for further inferences. In this paper, we adapt the imprecise Dirichlet model of Walley to represent epistemic uncertainty in the alpha-factors. In this approach, epistemic uncertainty is expressed more cautiously via lower and upper expectations for each alpha-factor, along with a learning parameter which determines how quickly the model learns from observed data. For this application, we focus on elicitation of the learning parameter, and find that values in the range of 1 to 10 seem reasonable. The approach is compared with Kelly and Atwood's minimally informative Dirichlet prior for the alpha-factor model, which incorporated precise mean values for the alpha-factors, but which was otherwise quite diffuse. Next, we explore the use of a set of Gamma priors to model epistemic uncertainty in the marginal failure rate, expressed via a lower and upper expectation for this rate, again along with a learning parameter. As zero counts are generally less of an issue here, we find that the choice of this learning parameter is less crucial. Finally, we demonstrate how both epistemic uncertainty models can be combined to arrive at lower and upper expectations for all common-cause failure rates. Thereby, we effectively provide a full sensitivity analysis of common-cause failure rates, properly reflecting epistemic uncertainty of the analyst on all levels of the common-cause failure model.

99 GENERAL AND MISCELLANEOUS↗

Scalar dark matter in the B−L model

The U(1){sub B−L} extension of the Standard Model requires the existence of right-handed neutrinos and naturally realizes the seesaw mechanism of neutrino mass generation. We study the possibility of explaining the dark matter in this model with an additional scalar field, ϕ{sub DM}, that is a singlet of the Standard Model but charged under U(1){sub B−L}. An advantage of this scenario is that the stability of ϕ{sub DM} can be guaranteed by appropriately choosing its B−L charge, without the need of an extra ad hoc discrete symmetry. We investigate in detail the dark matter phenomenology of this model. We show that the observed dark matter density can be obtained via gauge or scalar interactions, and that semi-annihilations could play an important role in the latter case. The regions consistent with the dark matter density are determined in each instance and the prospects for detection in future experiments are analyzed. If dark matter annihilations are controlled by the B−L gauge interaction, the mass of the dark matter particle should lie below 5 TeV and its direct detection cross section can be easily probed by XENON1T; if instead they are controlled by scalar interactions, the dark matter mass can be much larger and the detection prospects are less certain. Finally, we show that this scenario can be readily extended to accommodate multiple dark matter particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cogenerating and pre-annihilating dark matter by a new gauge interaction in a unified model

Grand unified theories based on large groups (with rank ≥6) are a natural context for dark matter models. They contain Standard-Model-singlet fermions that could be dark matter candidates, and can contain new non-abelian interactions whose sphalerons convert baryons, leptons, and dark matter into each other, “cogenerating' a dark matter asymmetry comparable to the baryon asymmetry. In this paper it is shown that the same non-abelian interactions can “pre-annihilate' the symmetric component of heavy dark matter particles χ, which then decay late into light stable dark matter particles ζ that inherit their asymmetry. We derive cosmological constraints on the parameters of such models. The mass of χ must be <3000 TeV and their decays must happen when 2×10{sup −7}<T{sub dec}/m{sub χ}<10{sup −4}. It is shown that such decays can come from d=5 operators with coefficients of order 1/M{sub GUT} or 1/M{sub Pℓ}. We present a simple realization of our model based on the group SU(7).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

PFLOTRAN modeling data and scripts associated with “Refining the Hydrogeologic Framework of a Large River Corridor Model Using Waterborne Transient Electromagnetics”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the publication “Refining the Hydrogeologic Framework of a Large River Corridor Model Using Waterborne Transient Electromagnetics” submitted to Water Resources Research (Terry et al. 2025). The data package contains the groundwater modeling dataset from PFLOTRAN software. It includes the python script for mesh generation, boundary condition setting, PFLOTRAN input deck formation and postprocessing. It couples groundwater flow and species transport for Hanford Reach river corridor and pipelines the model generation and processing. This model can be used to easily generate the model and analysis for Hanford site. It can also be adjusted to other hydrologic area with ease. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package consists of 6 folders: (1) “data” contains all necessary data as input and intermediate data for processing; (2) “mesh” contains all mesh related files to generate mesh in Hanford Reach river corridor; (3) “model_run” contains the generated script for PFLOTRAN modeling; (4) “notebooks” contains all the Python script to generate the model; (5) “output” contains all the output from the computation; (6) “postprocessing” contains the Python script to generate scientific figure for manuscript. All files are .csv (comma-separated values), .h5 (HDF5 format), .in (input files), .ipynb (Jupyter notebooks), .p (Python pickle), .png (images), .PNG (images), .py (Python scripts), .pyc (Python bytecode), .r (R scripts), .sh (shell scripts), .txt (text files), .vtu (3D mesh/visualization format), .xz (compressed archive), or .zip (compressed archive).

54 ENVIRONMENTAL SCIENCES↗

Data-model files associated with the manuscript "Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California"

This package contains the data, simulation setups, notebooks and figures used in “Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California” (Xu et al., 2025). In this study, we selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise on coastal hydrology using the process-based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site-specific data. We designed a novel modeling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modeled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected sea level rise over the next few decades. The low-lying Yampah Marsh, the reference site, is likely to be inundated due to future sea level rise when highest tides arrive; while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.Several files can be found from this data package.1. README.md: This file describes the title, journal, co-authors, abstract, repository structure and model version.2. Simulation_Setups.zip: The file contains the model configuration files (XML format) for ATS. 3. Notebooks.zip: The file contains the Jupyter notebooks for generating the pre- and post-restoration meshes and the meshes of future scenarios. 4. Figures.zip: The file contains the figures used in the manuscript.5. Data.zip: The file contains the data used to drive the model simulations, including watershed and wetlands boundaries, mesh files and references to additional datasets (e.g., meteorological forcing, tidal dataset, DEMs, land cover, soil properties). Also, it contains water level observations at the restored wetland.

54 ENVIRONMENTAL SCIENCES↗

Fractional Modeling in Action: A Survey of Nonlocal Models for Subsurface Transport, Turbulent Flows, and Anomalous Materials

Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully-resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas – subsurface transport, turbulence, and anomalous materials – in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.

97 MATHEMATICS AND COMPUTING↗

On approaches to analyze the sensitivity of simulated hydrologic fluxes to model parameters in the community land model

Here, effective sensitivity analysis approaches are needed to identify important parameters or factors and their uncertainties in complex Earth system models composed of multi-phase multi-component phenomena and multiple biogeophysical-biogeochemical processes. In this study, the impacts of 10 hydrologic parameters in the Community Land Model on simulations of runoff and latent heat flux are evaluated using data from a watershed. Different metrics, including residual statistics, the Nash-Sutcliffe coefficient, and log mean square error, are used as alternative measures of the deviations between the simulated and field observed values. Four sensitivity analysis (SA) approaches, including analysis of variance based on the generalized linear model, generalized cross validation based on the multivariate adaptive regression splines model, standardized regression coefficients based on a linear regression model, and analysis of variance based on support vector machine, are investigated. Results suggest that these approaches show consistent measurement of the impacts of major hydrologic parameters on response variables, but with differences in the relative contributions, particularly for the secondary parameters. The convergence behaviors of the SA with respect to the number of sampling points are also examined with different combinations of input parameter sets and output response variables and their alternative metrics. This study helps identify the optimal SA approach, provides guidance for the calibration of the Community Land Model parameters to improve the model simulations of land surface fluxes, and approximates the magnitudes to be adjusted in the parameter values during parametric model optimization.

58 GEOSCIENCES↗

CMIP5 model selection for ISMIP6 ice sheet model forcing: Greenland and Antarctica

Abstract. The ice sheet model intercomparison project for CMIP6 (ISMIP6) effort brings together the ice sheet and climate modeling communities to gain understanding of the ice sheet contribution to sea level rise. ISMIP6 conducts stand-alone ice sheet experiments that use space- and time-varying forcing derived from atmosphere–ocean coupled global climate models (AOGCMs) to reflect plausible trajectories for climate projections. The goal of this study is to recommend a subset of CMIP5 AOGCMs (three core and three targeted) to produce forcing for ISMIP6 stand-alone ice sheet simulations, based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century. The selection is performed separately for Greenland and Antarctica. Model evaluation over the historical period focuses on variables used to generate ice sheet forcing. For stage (i), we combine metrics of atmosphere and surface ocean state (annual- and seasonal-mean variables over large spatial domains) with metrics of time-mean subsurface ocean temperature biases averaged over sectors of the continental shelf. For stage (ii), we maximize the diversity of climate projections among the best-performing models. Model selection is also constrained by technical limitations, such as availability of required data from RCP2.6 and RCP8.5 projections. The selected top three CMIP5 climate models are CCSM4, MIROC-ESM-CHEM, and NorESM1-M for Antarctica and HadGEM2-ES, MIROC5, and NorESM1-M for Greenland. This model selection was designed specifically for ISMIP6 but can be adapted for other applications.

54 ENVIRONMENTAL SCIENCES↗