SEARCH · Search NASA

Results for “numerical 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.

SourceResearch connectionFrom yearThrough year

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

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Progressive collapse analysis of stiffened panel structure with consideration of actual welding distortion and residual stresses

To date, inclusion of welding distortion and residual stresses in determining the strength of marine structures is based on empirical formulations due to welding complexity. In this study, a more feasible numerical model was developed from the perspective of inherent strain theory. The accuracy of this simplified numerical model was confirmed by transient thermo-elastic–plastic finite element analysis using a novel dynamic mesh refining and iterative substructure method which is 50 times faster than conventional finite element method. Welding distortion by the two numerical approaches was validated with experimental measurement. The collapse analyses indicated that the inherent strain model predicted a compressive strength and load-displacement curve with high accuracy, as validated by the transient thermo-mechanical model. The global stiffness and the ultimate strength of the welded structure would be overestimated by about 22% if no welding effects were considered. The proposed computation procedure for determining the ultimate strength of welded structures is accurate, simple, and practical and enables improvement of ultimate strength through optimization of welding joint design to reduce the welding residual stress and distortion.

Huang, Hui↗

DOE OSTI2022

Can reanalysis products outperform mesoscale numerical weather prediction models in modeling the wind resource in simple terrain?

Mesoscale numerical weather prediction (NWP) models are generally considered more accurate than reanalysis products in characterizing the wind resource at heights of interest for wind energy, given their finer spatial resolution and more comprehensive physics. However, advancements in the latest ERA-5 reanalysis product motivate an assessment on whether ERA-5 can model wind speeds as well as a state-of-the-art NWP model – the Weather Research and Forecasting (WRF) Model. We consider this research question for both simple terrain and offshore applications. Specifically, we compare wind profiles from ERA-5 and the preliminary WRF runs of the Wind Integration National Dataset (WIND) Toolkit Long-term Ensemble Dataset (WTK-LED) to those observed by lidars at a site in Oklahoma, United States, and in a United States Atlantic offshore wind energy area. We find that ERA-5 shows a significant negative bias (~-1ms-1) at both locations, with a larger bias at the land-based site. WTK-LED-predicted wind speed profiles show a limited negative bias (~-0.5ms-1) offshore and a slight positive bias (~+0.5ms-1) at the land-based site. On the other hand, we find that ERA-5 outperforms WTK-LED in terms of the centered root-mean-square error (cRMSE) and correlation coefficient, for both the land-based and offshore cases, in all atmospheric stability conditions. We find that WTK-LED's higher cRMSE is caused by its tendency to overpredict the amplitude of the wind speed diurnal cycle. At the land-based site, this is partially caused by wind plant wake effects not being accurately captured by WTK-LED.

17 WIND ENERGY↗

DOE OSTI2023

Thermohydraulic experiments on a supercritical carbon dioxide–air microtube heat exchanger

Heat exchangers are critical components in supercritical CO 2 Brayton cycles and provide necessary heat input, recovery, and dissipation. Tubular heat exchangers with unconventionally small tube sizes (tube diameters less than 5 mm) are promising components for supercritical CO 2 cycles and potentially provide excellent structural stability with wide scope of application. This paper provides essential design and fabrication guidelines for a compact shell-and-tube heat exchanger with microtubes (with an inner diameter equal to 1.75 mm). A heat exchanger test rig is used to evaluate the thermohydraulic performance of this heat exchanger with supercritical CO 2 and air as working fluids. Thermohydraulic data are reported for more than forty sets of experiments with varying Reynolds numbers for shell and tube flows. Critical performance metrics are calculated from the data and compared with predictions from a previously described numerical model. The average deviations between the experimental and model results fall within 10% for all critical metrics. This excellent agreement validates the numerical model for supercritical CO 2 heat exchanger optimization and scale-up.

42 ENGINEERING↗

DOE OSTI2020

Analytic Thermal Model of an Optical Fiber Based Gamma Thermometer and its Application in a University Research Reactor

This paper describes and validates, by comparison with numerical modeling results, an analytical model of thermal transport in an optical fiber based gamma thermometer (OFBGT) that is appropriate for use in university research reactors. The maximum temperature difference between the thermal mass and the outer sheath ( Δ T ) for the OFBGT design that we have considered is approximately 50 ∘ C , for the OFBGT in the Central Irradiation Facility of the Ohio State University Research Reactor (OSURR) with the reactor operating at full power (450 kW). The maximum value of Δ T that is predicted by the analytic model for the OFBGT design is smaller by approximately 1.1 °C than the value of Δ T which is predicted by the numerical model, for the same OFBGT design, but including all the details of the design. We have used the analytical model of thermal transport in an OFBGT to determine a normalized Modulation Transfer Function M T F ′ ( k ) ) for the OFBGT. We conclude that MTF ′ ( k OSU ) > 0.99 , where k OSU is the spatial frequency for the axial dependence of the reactor power distribution in the

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Riverine Plastic Pollution: Field Sampling Protocol and Implementation in U.S. Rivers

Riverine plastic pollution is increasingly being recognized as a serious environmental concern. The flux of plastic from U.S. rivers into the oceans is estimated to be on the order of hundreds of tons per year. Plastic pollution has been detected in all major U.S. rivers but the observations did not measure metrics that could be used for mass flux calculations. The Department of Energy is addressing the problem of riverine plastic pollution with the Waterborne Plastics Resource Assessment and Debris Characterization (WaterPACT) project. Phase one of the WaterPACT project includes the development of a sampling protocol, field sampling events that implement the protocol, and numerical modeling studies. The field sampling will quantify and characterize the plastic pollution and the numerical modeling studies will analyze the movement of the plastic pollution by rivers. This report details the WaterPACT sampling protocol and its implementation in four major U.S. rivers. The WaterPACT sampling protocol was implemented in the Mississippi, Columbia, Delaware, and Los Angeles rivers with at least three sampling events in each river. The rivers represent four types ranging from the large Mississippi with its vast agricultural watershed to the small Los Angeles with its highly engineered urban waterway. The three sampling events covered a variety of discharge conditions from low drought flow to extreme storm events. The data obtained by the WaterPACT sampling of the four rivers will be combined to provide an invaluable dataset to researchers studying the source and chemical composition of plastic pollution and modelers estimating the flux of plastic being transported to and released into the oceans.

13 HYDRO ENERGY↗

DOE OSTI2023

Non-linear Wire Rope Isolator Model to Enhance Transportation Simulation of Fragile Equipment

Abstract Introduction Wire rope isolators are often used as means of vibration isolation for the transportation of fragile machine equipment. When deployed as suspension systems in a transportation scenario, a simple spring–damper model is often used to predict the dynamic loads that act on the structure being transported. However, this simple model has proved to be too simplistic to be used in the development of an accurate numerical model that can predict the vibration levels experienced by the inner components of fragile equipment during its transportation. This paper describes the experimental tests conducted on a wire rope isolator used for the transportation of the prototype SSR1 cryomodule. Materials and methods A hysteretic Bouc–Wen model has been used to analytically describe the force–deformation relationship of the wire rope isolators. The developed model of the isolator has been implemented in a larger model to simulate the actual transportation of the prototype SSR1 cryomodule, a section of the new PIP-II linear accelerator under construction at Fermilab. A series of multibody dynamic simulations with rigid and flexible components was used to numerically determine the acceleration of some critical components. Results and conclusions An actual experimental transportation was simulated using two numerical models: the developed Bouc–Wen model and a conventional spring–damper model. It is shown how the Bouc–Wen formulation of the isolator characteristics drastically improves the correspondence between experimental and simulated results if compared to a spring–damper model, especially in the range of 0–30 Hz which is the most critical for transportation problems.

42 ENGINEERING↗

DOE OSTI2021

Leveraging machine learning to improve understanding and predictability of weather/climate extremes and the resilience of human systems

Focal Area(s): Developing hybrid models for predicting weather/climate extremes at the intersection of the urban environment (focal area #2). Science Challenge: Limited by many sources of uncertainty, numerical models are unable to accurately simulate weather and climate processes and this problem is particularly acute for modeling extreme events, including hurricanes. Moreover, computational limits within these numerical models bound our ability to develop large ensembles to explore this uncertainty and the implications relative to characterization of risk and resilience within the human systems that are impacted by changes in extremes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Progress Report on Model Development for the Transport of Aerosol Through Microchannels

Accuracy of consequence assessments could be improved by accounting for the leak path deposition of aerosol in the source term. Recent experimental measurements from Sandia (2018) provided evidence of partial microchannel plugging due to particulate deposition. This filtration effect is relevant to the source term assessments and it is not unreasonable to anticipate that aerosol deposition in more realistic geometries could even lead to complete plugging of the leak path. This report summarizes current progress (as of end of Q3 FY 2020) on the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation. This report includes (1) an overview of recent additions and improvements on the recently developed numerical model to analyze the various deposition processes in leak paths, to provide quantitative estimates of penetration factors, and to gain an understanding of the variables that affect them and (2) a summary of recent benchmarking results of the model’s validity with recent experimental measurements from Sandia.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (< km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ~20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

58 GEOSCIENCES↗

DOE OSTI2020

Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones

Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one-way coupling between the surface and subsurface domains. Such one-way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam, was developed using the open-source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod-type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three-dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.

58 GEOSCIENCES↗

DOE OSTI2021

Machine Learning Analysis of Hydrologic Exchange Flows and Transit Time Distributions in a Large Regulated River

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. High-resolution numerical models were often used to resolve the spatial and temporal variations of exchange flows, which are computationally expensive. In this study, we adopt Random Forest (RF) and Extreme Gradient Boosting (XGB) approaches for deriving reduced order models of hydrologic exchange flows and associated transit time distributions, with integrated field observations (e.g., bathymetry) and hydrodynamic simulation data (e.g., river velocity, depth). The setup allows an improved understanding of the influences of various physical, spatial, and temporal factors on the hydrologic exchange flows and transit times. The predictors also contain those derived using hybrid clustering, leveraging our previous work on river corridor system hydromorphic classification. The machine learning-based predictive models are developed and validated along the Columbia River Corridor, and the results show that the top parameters are the thickness of the top geological formation layer, the flow regime, river velocity, and river depth; the RF and XGB models can achieve 70% to 80% accuracy and therefore are effective alternatives to the computational demanding numerical models of exchange flows and transit time distributions. Each machine learning model with its favorable configuration and setup have been evaluated. The transferability of the models to other river reaches and larger scales, which mostly depends on data availability, is also discussed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Compaction of crushed salt for safe containment – overview of the KOMPASS project

Abstract. In Germany, rock salt formations are possible host rock candidates for a repository for heat-emitting radioactive waste. The safety concept of a repository in salt bases on a multibarrier system consisting mainly of the geological barrier salt and geotechnical seals ensuring safe containment. Crushed salt will be used for backfilling of cavities and sealing measures in drifts and shafts due to its favourable properties and its easy availability (mined-off material). The creep of the rock salt leads to crushed salt compaction with time. Thereby, the crushed salts' porosity is reduced from the initial porosity of 30 %–40 % to a value comparable to the porosity of undisturbed rock salt (≤1 %). In such low porosity ranges, technical impermeability is assumed. The compaction behaviour of crushed salt is rather complex and involves several coupled THM processes (Kröhn et al., 2017; Hansen et al., 2014). It is influenced by internal properties like humidity and grain size distribution, as well as boundary conditions such as temperature, compaction rate or stress state. However, the current process understanding has some important gaps referring to the material behaviour, experimental database and numerical modelling. It needs to be extended and validated, especially in the low porosity range. The objective of the KOMPASS project was development of methods and strategies for the reduction of deficits in the prediction of crushed salt compaction leading to an improvement of the prognosis quality. Key results are as follows (KOMPASS Phase 1, 2020): selection of an easily available and permanently producible synthetic crushed salt mixture, acting as a reference material for generic investigations; development and proof of different techniques for producing pre-compacted samples for further investigations; establishment of a tool of microstructure investigation methods to demonstrate the comparability of grain structures of pre-compacted samples with in-situ compacted material for future investigations; execution of various laboratory experiments using pre-compacted samples, e.g. long-term creep tests which deliver reliable information about time- and stress-dependent compaction behaviour; development of a complex experimental investigation strategy to derive necessary model parameters considering individual functional dependencies. Its technical feasibility was successfully verified; benchmarking with various existing numerical models using datasets from three different triaxial long-term tests. The result was not entirely satisfactory; however, the number of influencing factors is small and further validation work has to be done. Overall, the KOMPASS project has made significant progress in the approaches to solving the outstanding question, building the basis for further investigations.

Friedenberg, Larissa↗

DOE OSTI2020

Investigation of Failure Modes, Mechanisms and Driving Forces for Electrically Conductive Adhesives as Interconnects in PV Modules

In this work we describe our current progress into investigating the reliability of ECA interconnects for PV modules. Preliminary thermal cycling experiments demonstrated that both cohesive and adhesive failure of the ECA joint are relevant failure modes that result in an increase of the joints electrical resistance. A thorough viscoelastic characterization of the candidate ECA was conducted to enable accurate numerical modeling of the joint’s response to thermal cycling. Combined with critical and subcritical fracture measurements, the numerical model elucidates that an additional fatigue mechanism exists which can result in fracture well below subcritical debonding threshold values. It is also demonstrated that fracture of the ECA is more sensitive to humidity than temperature and subcritical debonding mechanisms do exist which are again more sensitive to humidity.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2022

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Teaching Freight Mode Choice Models New Tricks Using Interpretable Machine Learning Methods

Understanding and forecasting the intricate freight mode choice behavior under various industry, policy, and technology contexts is essential in freight planning and policymaking. Numerous models have been developed in prior studies to provide insights into freight mode selection, the majority of which use discrete choice models such as multinomial logit (MNL) models. However, logit models often rely on linear specifications of independent variables, despite potential nonlinear relationships in the data. Moreover, there often lacks a heuristic and efficient approach to identify such complex relationships to define the logit model specifications. To fill this gap, we developed an MNL model for freight mode choice using the insights from state-of-the- art machine learning (ML) models. ML models can capture the nonlinear nature of the complex decision-making process, and recent advances in 'explainable AI' have greatly improved their interpretability. The interpretable ML methods help enhance the performance of MNL models and advance knowledge of freight mode choice. Specifically, the influential factors and their relationship with individual modes are identified using SHapley Additive exPlanations (SHAP) to improve the MNL's performance. The workflow is demonstrated in a case study of Austin, Texas, and the SHAP results reveal multiple nonlinear relationships predicted by ML models. Incorporating those relationships into MNL model specifications improves the interpretability and accuracy of the MNL model compared to a conventional MNL model. Findings from this study can be used to guide freight planning and inform policymakers and practitioners on how key factors affect freight decision-making.

ADVANCED PROPULSION SYSTEMS,MATHEMATICS AND COMPUT↗

DOE OSTI2023

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗

DOE OSTI2023

Physics-constrained Gaussian process model for prediction of hydrodynamic interactions between wave energy converters in an array

To improve the efficiency of wave farms and achieve maximum power generation, the layout of wave energy converters (WECs) in an array needs to be carefully designed so that the hydrodynamic interactions can be positively exploited. For this, the hydrodynamic characteristics of the WEC array in different layouts need to be calculated. However, such calculations using numerical models usually entail significant computational cost, especially for large arrays of WECs. To address the computational challenge, a physics-constrained Gaussian process (GP) model is proposed to replace the original expensive numerical model and predict the hydrodynamic characteristics of the WECs for any array layout. By exploring the relationship between the WEC array (i.e., the input) and different hydrodynamic characteristics (i.e., the output), here we summarize a set of physical constraints/features, including invariance, symmetry, and additivity. This prior knowledge about the input-output relationship is then directly embedded in the constructed GP model through the design of physics-constrained kernels. In particular, a double-sum invariant kernel is first developed to incorporate the invariance and symmetry features, and then an additive kernel is developed to incorporate the additive feature of the problem. The invariant kernel and the additive kernel are then integrated to construct the physics-constrained GP model. Compared to the standard GP model, the proposed physics-constrained GP models require less training data to achieve the desired accuracy in predicting the hydrodynamic characteristics and are also less vulnerable to the curse of dimensionality (i.e., good scalability for large arrays) due to the use of an additive kernel. The efficiency, accuracy, and scalability of the proposed approach are demonstrated through an application to predict the hydrodynamic characteristics for WEC arrays of different sizes and layouts.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

Spontaneous imbibition in a square tube with corner films: theoretical model and numerical simulation

Spontaneous imbibition in an angular tube with corner films is a fundamental problem in many scientific and engineering processes. In this paper, a modified interacting capillary bundle model is developed to describe the liquid imbibition dynamics in a square tube with corner films. The square tube is decomposed into several interacting sub-capillaries and the local capillary pressure in each sub-capillary is derived based on the specific shape of its meniscus. The conductance of each sub-capillary is calculated using single-phase lattice Boltzmann simulation. The modified interacting capillary bundle model and color-gradient lattice Boltzmann method are used to simulate the liquid imbibition dynamics in the square tube with different fluid properties. The predictions by the modified interacting capillary bundle model match well with the lattice Boltzmann simulation results for different conditions, demonstrating the accuracy and robustness of the interacting capillary bundle model to describe the imbibition dynamics with corner films. In addition, the interacting capillary bundle model is helpful to investigate the mechanisms during spontaneous imbibition and the influences of fluid viscosity, surface tension, wetting phase contact angle and gravity on imbibition dynamics. Finally, a universal scaling law of imbibition dynamics for the main meniscus is developed and the scaling law for arc meniscus is also analyzed.

42 ENGINEERING↗