SEARCH · Engineering Papers

Results for “numerical modelling”

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 199 records · Page 11

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Lattice Boltzmann simulation of the dissolution of slag in alkaline solution using real-shape particles

Highlights: • A dissolution numerical model was proposed to capture the real dissolution kinetics of slag in alkaline solution. • The log forward dissolution rate of Si was described as a function of NBO/T and solution pH. • A threshold solid volume fraction of 0.688 was found for a voxel in 3D, 63.8% larger than that for a pixel in 2D. • The proposed dissolution numerical model provides a reliable alternative to study the dissolution kinetics of slag. A dissolution numerical model was proposed in this study to capture the real dissolution kinetics of slag in alkaline solution. It consists of three modules, i.e. (i) simulation of the initial particle parking structure of slag in alkaline solution using real-shape particles of slag, (ii) simulation of the chemical reactions between slag and solution based on the transition state theory, and (iii) simulation of the physical transport of aqueous ions using the lattice Boltzmann method. This dissolution numerical model was verified using experimental results, showing reasonable accuracy. After verification, the dissolution numerical model was implemented to study the influences of temperature and particle shape using a proper recipe of slag in alkaline solution. This recipe was designed to avoid solid phase precipitation or gel formation via thermodynamic analysis. The simulation results showed faster dissolution kinetics of slag when using higher temperatures and more irregular particle shapes.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Techno-economic feasibility of borehole thermal energy storage system connected to geothermal heat pumps for seasonal heating load of two buildings in Fairbanks, Alaska

Borehole thermal energy storage (BTES) is an effective solution for managing imbalanced heating and cooling loads in cold regions. This study evaluated the long-term feasibility of a BTES system in the Fairbanks area, Alaska, through building energy modeling, resource characterization, and numerical modeling. The system was designed to store waste heat from a nearby coal power plant during summer and provide thermal energy during winter to geothermal heat pumps supplying heating loads in two buildings. Heating load profiles were modeled for the buildings using EnergyPlus, and the results indicated the annual heating load was 5.6 times greater than the cooling load. 40 borehole heat exchangers were pre-designed approximately 100 m away from the two buildings in terms of land availability and regulatorily optimized depth. The 20-year performance of the designed BTES system under two operational scenarios—one with a 5-year preheating period and one without—was numerically modeled using subsurface temperature and properties characterized through the literature review and thermal response tests. Both scenarios demonstrated that the BTES has the capacity to fully cover the heating loads in the two buildings throughout the 20-year lifetime. Production temperatures at central wells were 33 % higher on average than at outer wells in both scenarios. The 5-year preheating period increased subsurface and extraction temperatures, and correspondingly annual average and total thermal energy production was higher for 8 years than in the scenario without the preheating period. These results highlight the long-term reliability and sustainability of the BTES system in meeting heating demands over its lifetime, with the preheating period offering potential performance improvements. Implementing the BTES system in cold regions with high heating demand, such as Fairbanks, Alaska, could provide a long-term, sustainable energy solution for managing imbalanced heating and cooling loads.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Scalable foundation models for numerical simulations on HPC platforms

In recent years, foundation models (FMs) have begun to reshape numerical simulations on high-performance computing (HPC) platforms. These large, pre-trained AI models enable rapid predictions across a broad range of physical domains, including Earth system modeling, fluid dynamics, materials science, as well as complex multi-modal simulations in aerospace engineering and fusion research. By training on diverse datasets, FMs learn intricate relationships and underlying physical behavior while also enabling the quantification of uncertainty in their predictions. This capability allows simulations that once required days of numerical calculation to be completed in minutes (FM inference), supporting real-time design optimization, uncertainty-aware decision making, and more comprehensive exploration of complex scenarios.

AI↗

DOE OSTI2020

Using numerical methods to screen magnetocaloric materials in an active magnetic regenerative cycle

A 2-dimensional numerical model of a room temperature operating Active Magnetic Regenerator (AMR) that comprises of a regenerator, hot and cold heat exchangers, heat transfer fluid is developed. The regenerator is made of a magnetocaloric material (MCM) which heats up upon applying a magnetic field, H, and cools down when the field is removed; thus, making it the most essential part of an AMR. The model takes experimentally measured ΔT ad (H,T) and the C p (H,T) data as input and provides quantitative performance metrics of the magnetic cooling system, such as ΔT span and the cooling load, as output. With this model, it is possible to assess a wide range of MCMs in an AMR. 4 different MCMs were investigated using this model in terms of their ΔT span and the cooling loads—LaFe 10.96 Co 0.97 Si 1.07 , MnFeP (1-x) As (x) and AlFe 2 B 2 and Gd. During the screening of MCMs, all the important operating conditions of the device were fixed such as the fluid flow rate, ambient temperature, cycle duration, magnetic field strength. Our results indicate that Gd exhibits the maximum ΔT span with respectable cooling load, while AlFe 2 B 2 generates the lowest ΔT span . Even though LaFeCoSi compound did not perform as well as Gd, it could be the MCM choice for the realization of magnetic refrigeration on a global scale due to its much lower cost.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A two-fluid solar-wind model with intermittent Alfvénic turbulence

In one of the leading theories for the origin of the solar wind, photospheric motions launch Alfvén waves (AWs) that propagate along open magnetic-field lines through the solar atmosphere and into the solar wind. The radial variation in the Alfvén speed causes some of the AWs to reflect, and counter-propagating AWs subsequently interact to produce Alfveńic turbulence, in which AW energy cascades from long wavelengths to short wavelengths and dissipates, heating the plasma. In this paper we develop a one-dimensional two-fluid solar-wind model that includes Alfvénic turbulence, proton temperature anisotropy and a novel method for apportioning the turbulent heating rate between parallel proton heating, perpendicular proton heating and electron heating. We employ a turbulence model that accounts for recent observations from NASA’s Parker Solar Probe, which find that AW fluctuations in the near-Sun solar wind are intermittent and less anisotropic than in previous models of anisotropic magnetohydrodynamic turbulence. Our solar-wind model reproduces a wide range of remote observations of the corona and in-situ measurements of the solar wind, and our turbulent heating model consists of analytic equations that could be usefully incorporated into other solar-wind models and numerical models of more distant astrophysical plasmas.

astrophysical plasmas↗

DOE OSTI2021

Optimization and supervised machine learning methods for fitting numerical physics models without derivatives

Here, we address the calibration of a computationally expensive nuclear physics model for which derivative information with respect to the fit parameters is not readily available. Of particular interest is the performance of optimization-based training algorithms when dozens, rather than millions or more, of training data are available and when the expense of the model places limitations on the number of concurrent model evaluations that can be performed. As a case study, we consider the Fayans energy density functional model, which has characteristics similar to many model fitting and calibration problems in nuclear physics. We analyze hyperparameter tuning considerations and variability associated with stochastic optimization algorithms and illustrate considerations for tuning in different computational settings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗

DOE OSTI2024

Recent Progress on Surface Water Quality Models Utilizing Machine Learning Techniques

Surface waterbodies are heavily exposed to pollutants caused by natural disasters and human activities. Empowering sensor technologies in water quality monitoring, sufficient measurements have become available to develop machine learning (ML) models. Numerous ML models have quickly been adopted to predict water quality indicators in various surface waterbodies. This paper reviews 78 recent articles from 2022 to October 2024, categorizing water quality models utilizing ML into three groups: Point-to-Point (P2P), which estimates the current target value based on other measurements at the same time point; Sequence-to-Point (S2P), which utilizes previous time series data to predict the target value at one time point ahead; and Sequence-to-Sequence (S2S), which uses previous time series data to forecast sequential target values in the future. The ML models used in each group are classified and compared according to water quality indicators, data availability, and model performance. Widely used strategies for improving performance, including feature engineering, hyperparameter tuning, and transfer learning, are recognized and described to enhance model effectiveness. The interpretability limitations of ML applications are discussed. This review provides a perspective on emerging ML for surface water quality models.

machine learning (ML)↗

DOE OSTI2022

Graphite Oxidation Modeling

Oxidation modeling importance ASME consideration Effect on graphite properties Model Overview Numerical Model Oxidation Diffusion Reaction Kinetics Model Parameterization Example problems Editing the model Simple problem example Conclusions Coupling to stress analysis Oxidation result analysis

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Field Laboratory for Emerging Stacked Unconventional Plays in Central Appalachia

The goal of the Field Laboratory for Emerging Stacked Unconventional Plays (ESUP) in Central Appalachia project was to investigate and characterize the resource potential for multi-play production of emerging unconventional reservoirs in Central Appalachia. Project activities included drilling, logging and coring of a vertical characterization well drilled to basement to approximately 15,000 feet in depth. The data from drilling and well logs as well recovered core samples were used by the research team to support characterization of the geology and potential pay zones within the Nora Gas Field of southwestern Virginia and the greater Central Appalachian Basin. This project was led by Dr. Nino Ripepi of the Virginia Center for Coal and Energy Research (VCCER) in close collaboration with EnerVest Operating LLC (EnerVest). VCCER is housed in the Department of Mining and Minerals Engineering at Virginia Tech and was created by an Act of the Virginia General Assembly on March 30, 1977, as an interdisciplinary study, research, information, and resource facility for the Commonwealth of Virginia. EnerVest is a top tier, low-cost oil and natural gas company with a long history of traditional operating relationships with institutional investors as well as public and private companies. This report summarizes two major research activities: 1) Log and Core Analysis and 2) Numerical Modeling. Core analysis comprises the characterization efforts of the ESUP Field Laboratory and includes a discussion of the well logs run in the deep characterization wells, an initial analysis of those, the number, size and location of cores recovered, core testing and analysis done to date as well as plans for future analysis. Numerical modeling is comprised of seven major modeling efforts, including (1) Multiphysics shale transport modeling work, (2) the developed in-house compositional simulator studies, (3) reservoir modeling studies using GEM software, (4) fracture modeling work using EFRAC3D software and (5) numerical modeling efforts with ABAQUS software, (6) Numerical modeling efforts using FLAC3D and PFC3D, and (7) Automatic Machine Learning (AutoML) Studies.

02 PETROLEUM↗

DOE OSTI2022

A deformation-dependent coupled Lagrangian/semi-Lagrangian meshfree hydromechanical formulation for landslide modeling

Abstract The numerical modelling of natural disasters such as landslides presents several challenges for conventional mesh-based methods such as the finite element method (FEM) due to the presence of numerically challenging phenomena such as severe material deformation and fragmentation. In contrast, meshfree methods such as the reproducing kernel particle method (RKPM) possess unique features conducive to modelling extreme events such as the absence of a structured mesh and the ease of adaptive refinement, among others. While the semi-Lagrangian reproducing kernel (SL-RK) shape functions of RKPM defined in the current configuration have proven to be effective in extreme event modelling, the computational cost for the re-evaluation of the shape functions at every time step is costly. In this work, a deformation-dependent coupling of the Lagrangian reproducing kernel (L-RK) and SL-RK approximations is proposed for the solution of a hydro-mechanical formulation for effective simulations of landslides. The ramp function is constructed based on an equivalent plastic strain as a deformation-dependent transition from L-RK shape functions to SL-RK ones as the deformation progresses. The particular focus of the paper will be on modelling seepage-induced landslides with a mixed $$u$$ u – $$p$$ p formulation to couple the solid and fluid phases. Examples are presented to examine the effectiveness of this coupled Lagrangian/semi-Lagrangian reproducing kernel (L–SL RK) formulation and to highlight its performance in landslide modelling.

Baek, Jonghyuk↗

DOE OSTI2023

Review of TEAMER Awards for WEC-Sim Support: Preprint

Testing Expertise and Access for Marine Energy Research (TEAMER) is a U.S. Department of Energy Water Power Technologies Office sponsored program, overseen by the Pacific Ocean Energy Trust, which aims to advance the state of marine energy technologies. The program connects technology developers with experts at U.S. facilities, including numerical modeling and analysis facilities, like WEC-Sim. The WEC-Sim facility is supported by the WEC-Sim development team at Sandia National Laboratories and the National Renewable Energy Laboratory. WEC-Sim (Wave Energy Converter SIMulator) is an open-source software for simulating wave energy converters. WEC-Sim can model the multi-body dynamics of devices comprised of bodies, joints, power take-off systems, and mooring systems. Since TEAMER's first round of support in 2020, there have been eighteen TEAMER awards focused on numerical model development in WEC-Sim. TEAMER awards for WEC-Sim support have modeled a wide range of wave energy converter archetypes, including point absorbers, attenuators, oscillating water columns, and many other novel architectures. A wide variety of studies have been conducted, leading to important insights for TEAMER partners and software improvements for WECSim. This article highlights several successful WEC-Sim TEAMER awards. The awards described herein include TEAMER recipients Ocean Motion Technologies, AquaHarmonics, iProTech, East Carolina University, Virginia Tech, Maiden Wave Energy, and the University of Massachusetts Dartmouth. The awards of these seven partners contain a wide range of investigations and cover the creation of baseline hydrodynamic models, PTO modeling, geometry optimization in both boundary element methods and WECSim, and model tuning and validation.

industry support↗

DOE OSTI2024

A Panoramic View of Temperature and Field Distributions of the Structured Catalyst Under Microwave Irradiation Using Experimental and Modeling Approaches

This is a presentation covers NETL's recent research outcome in numerical modeling of microwave-assisted catalytic reaction. In this presentation, the intrinsic microwave-material interaction was demonstrated in this study through numerical modeling and experimental measurement. Different material combinations and reactor setups were investigated and compared. By combining numerical modeling with advanced measurement techniques, the microwave-material interaction and the relationship between electromagnetic field and heating can be better understood, which further benefits the development of microwave-assisted reactors and their application in process intensification.

Bai, Xinwei↗

DOE OSTI2022

Yield Estimation of the August 2020 Beirut Explosion by Using Physics‐Based Propagation Simulations of Regional Infrasound

Abstract The August 2020 Beirut explosion is the largest single‐fired ammonium nitrate explosion documented in history. The massive explosion excited loud infrasound in the atmosphere, and clear waveforms were recorded by a regional infrasound array at an epicentral distance of 100 km, allowing for accurate measurements of explosion energy. We estimate the explosion size based on the infrasound waveform inversion. Unlike conventional inversions using empirical models, we perform full 3‐D finite‐difference simulations to obtain a physics‐based propagation model for the inversion. Accurate numerical modeling of infrasound is challenging as the propagation is substantially affected by the turbulent atmosphere. Instead of a single deterministic prediction, we provide a range of waveform predictions by running multiple simulations with stochastic weather forecast models, which allows for comprehensive uncertainty analysis of numerical modeling and estimated yields. Finally, we expand the yield estimation technique for seismoacoustic analysis and demonstrate the substantial advantage of the joint approach.

58 GEOSCIENCES↗

DOE OSTI2022

CondiDiag1.0: a flexible online diagnostic tool for conditional sampling and budget analysis in the E3SM atmosphere model (EAM)

Abstract. Numerical models used in weather and climate prediction take into account a comprehensive set of atmospheric processes (i.e., phenomena) such as the resolved and unresolved fluid dynamics, radiative transfer, cloud and aerosol life cycles, and mass or energy exchanges with the Earth's surface. In order to identify model deficiencies and improve predictive skills, it is important to obtain process-level understanding of the interactions between different processes. Conditional sampling and budget analysis are powerful tools for process-oriented model evaluation, but they often require tedious ad hoc coding and large amounts of instantaneous model output, resulting in inefficient use of human and computing resources. This paper presents an online diagnostic tool that addresses this challenge by monitoring model variables in a generic manner as they evolve within the time integration cycle. The tool is convenient to use. It allows users to select sampling conditions and specify monitored variables at run time. Both the evolving values of the model variables and their increments caused by different atmospheric processes can be monitored and archived. Online calculation of vertical integrals is also supported. Multiple sampling conditions can be monitored in a single simulation in combination with unconditional sampling. The paper explains in detail the design and implementation of the tool in the Energy Exascale Earth System Model (E3SM) version 1. The usage is demonstrated through three examples: a global budget analysis of dust aerosol mass concentration, a composite analysis of sea salt emission and its dependency on surface wind speed, and a conditionally sampled relative humidity budget. The tool is expected to be easily portable to closely related atmospheric models that use the same or similar data structures and time integration methods.

58 GEOSCIENCES↗

DOE OSTI2022

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Development, Monitoring, and Control of Fracture Thermal Energy Storage (FTES) in Crystalline Rock Formations (DEMO-FTES) (CRADA Final Report)

The DEMO-FTES project sought to demonstrate the thermal efficiency of fracture thermal energy storage (FTES) through numerical simulations, laboratory and meso-scale field tests. A detailed dimensional and scaling analysis was performed to identify key parameters and how they can be most effectively scaled to the laboratory and decameter scale. Numerical models were developed and used for three purposes: 1. Before field experiments, numerical modelling can be used to estimate fracture properties based on previous data from the EGS Collab experiment and then predict thermal hydrological behaviors of the fracture system with hot water injection/withdrawal, therefore, help to design the experiments (e.g., to decide the duration of the cycles based on the flow rate the pump can provide, and the estimated fracture properties); 2. After the field experiment, to estimate the system properties during the experiment (as the size and shape of a fracture could change over time), and understand why system performance is different than what has been predicted, i.e., to help understand the meso-scale test; and, 3. To model the lab experiments and estimate fracture properties and storage efficiency. Ultimately, the experiment and numerical models could shed light on the processes and uncertainty happening during fracture activation and help understand the scaling between lab and field test, and finally, the design and optimization of potential fracture thermal energy storage systems.

25 ENERGY STORAGE↗

DOE OSTI2021

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗