Engineering Papers⌕ Search

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.

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Characterizing falling particle curtain receivers at commercially relevant scales: Research Performance Progress Report (RPPR-1)

Sandia will construct a cold flow receiver test platform in order to characterize falling particle curtain receivers at commercially relevant scales. While Sandia has extensive experience in R&D of falling particle curtain receivers, most have been at pilot scale and smaller - on the order of 1 MW th with characteristic dimensions of nominally 1-2 m to adequately collect solar energy from the heliostat field at the National Solar Thermal Test Facility (NSTTF). However, scaling up receivers to commercially relevant scales (25 MW th and above) will require a thorough understanding of particle curtain dynamics at larger scales, especially longer drop heights, for design certainty. The goal of this project will be to construct a cold falling particle curtain test rig capable of simulating particle characteristics that are expected in a commercial scale CSP plant, namely the drop height, curtain thickness, and particle mass flow rate (normalized by length of curtain). This will enable data collection on curtain opacity and spread, both of which are correlated to receiver efficiency and reliable construction, for commercially relevant scales. It will also permit validation of numerical models that will enable detailed receiver characterization and design past currently validated scales.

14 SOLAR ENERGY↗
DOE OSTI2024

Transfer Tube Material Properties Obtained from Test Data

A series of tension and compression tests were performed on typical transfer tubes made of 304L corrosion-resistant steel. The results were compiled and converted to true stress and true plastic strain for comparison to numerical models. The parameters of a Johnson-Cook constitutive model and fracture criterion were calibrated to mimic the test data. Comparisons were made between the converted data and the empirical formulas in the Johnson-Cook model. Simulations were performed to evaluate the implementation of the formula in finite element analysis, and the parameters were adjusted to improve the correlation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

University of Hawai‘i, Shallow Geothermal Resources: Energy Technology Innovation Partnership Project (Final Report)

Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that, over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

SAM User's Guide

The System Analysis Module (SAM) is a modern system analysis tool being developed at Argonne National Laboratory for advanced non-LWR safety analysis. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for Sodium-cooled Fast Reactors (SFR), Lead-cooled Fast Reactors (LFR), and Molten Salt Reactors (MSR) or Fluoride-cooled High-temperature Reactors (FHR). SAM takes advantage of advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability. It utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage the modern advanced software environments and numerical methods. This document provides a user’s guide, which will help users understand the input description and core capabilities of the SAM code. A brief overview of the code is presented, as well as how to obtain and run it. The input syntax for various parts of the code is provided. Additionally, a number of example problems, starting with simple unit component problems to problems with increasing complexity, are provided. Because the code is still under active development, this SAM User’s Guide will evolve with periodic updates.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Q-Fibre Quarterly Performance Progress Report

The project titled Quantitative Strain Monitoring for Well and Seal Integrity through Wellbore Fiber Optic Cable s, or Q-Fibre in short, is a multi-national collaborative project under the Clean Energy Transitional Partnership (CETP) program. The lead organization of the project is the Norwegian Geotechnical Institute (NGI), and the LLNL team joins as a beneficiary member. The LLNL team receives the funding from the U.S. Department of Energy with the amount of $\$$1 million across three years, and the budget is primarily allocated to Work Package 3 (WP3) of the project, which focuses on Numerical Modeling of Fiber Optic Cable Data . Figure 1 presents all work packages involved in the project and corresponding lead organization for each work package. The LLNL team leads WP3.

42 ENGINEERING↗
DOE OSTI2025

Characterizing falling particle curtain receivers at commercially relevant scales: Research Performance Progress Report (RPPR-1)

Sandia will construct a cold flow receiver test platform in order to characterize falling particle curtain receivers at commercially relevant scales. While Sandia has extensive experience in R&D of falling particle curtain receivers, most have been at pilot scale and smaller - on the order of 1 MWth with characteristic dimensions of nominally 1-2 m to adequately collect solar energy from the heliostat field at the National Solar Thermal Test Facility (NSTTF). However, scaling up receivers to commercially relevant scales (25 MWth and above) will require a thorough understanding of particle curtain dynamics at larger scales, especially longer drop heights, for design certainty. The goal of this project will be to construct a cold falling particle curtain test rig capable of simulating particle characteristics that are expected in a commercial scale CSP plant, namely the drop height, curtain thickness, and particle mass flow rate (normalized by length of curtain). This will enable data collection on curtain opacity and spread, both of which are correlated to receiver efficiency and reliable construction, for commercially relevant scales. It will also permit validation of numerical models that will enable detailed receiver characterization and design past currently validated scales.

14 SOLAR ENERGY↗
DOE OSTI2025

Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS) - Deliverable Task 7: Risk Assessment

The “Optimizing Alabama’s CO 2 Storage in Shelby County: Project OASIS” CarbonSAFE Phase II Project seeks to build on regional data sets that demonstrate that the subsurface within Shelby County, Alabama has the potential to store commercial volumes of CO 2 safely, permanently, and economically. The primary target CO 2 Storage Complex is the deep Ketona Dolomite located within a 140 square mile area of the Valley and Ridge Region of Alabama. This deep saline reservoir is beneath a confining system encompassing at least 6,500 ft of shales and other low permeability sediments. Project OASIS drilled a deep stratigraphic test well to confirm the geological properties of the confining system and storage reservoir(s) within the Storage Complex. The geological data was incorporated into numerical models to establish the areal extent of the CO 2 plume and help design the storage site and its monitoring system. The project is managed by the Southern States Energy Board (SSEB), an interstate compact organization consisting of governors and state legislative leaders from sixteen southern states, Puerto Rico, and the U.S. Virgin Islands, as well as an appointee by the President of the United States. The organizational compact provides it with access to state government organizations and legislatures. The Board also maintains an Associate Members program comprised of energy resource companies, utilities, trade groups, academic R&D science and technology experts and energy consultants. Further, the SSEB staff is experienced in managing and coordinating complex energy and environmental programs, from research programs to full-scale design and demonstrations of new and innovative technologies. Project OASIS is a public-private partnership of six entities with multiple principal investigators (PIs). SSEB’s Lead PI and Co-PI are responsible for all aspects of project performance in accordance with the DOE-NETL Cooperative Agreement. SSEB has issued subgrants to Advanced Resources International, Inc., Alabama A&M University, Auburn University, Crescent Resource Innovation, and Oklahoma State University. Advanced Resources International, Inc., issued subgrants to Baker Hughes and Loudon Technical Services for field services

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Modular Hydronic Room Conditioning System

This project presents the development and validation of high-fidelity numerical models for fin-tube heat exchangers to enable accurate performance prediction and informed design optimization. The modeling framework integrates detailed geometric specifications, thermophysical property data, and system-level constraints to simulate the heat exchanger’s behavior under a range of operating conditions. The model is calibrated using real-world product specifications and validated against experimental data collected from controlled cooling and heating tests. In cooling mode, the model captures the overall trends in capacity and outlet air temperature but tends to underpredict latent effects, especially at lower air flow rates. In heating mode, the simulation consistently overestimates both the thermal capacity and outlet air temperature, indicating the need for refinement in air-side heat transfer assumptions. Despite these deviations, the model provides a solid foundation for optimization, allowing key design variables to be tuned within physical and performance-based constraints. This research advances the ability to simulate, validate, and optimize fin tube heat exchanger designs with greater confidence and efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation of a High-Confinement Tidal Energy Array (Abstract)

When the projected area of a current turbine array occupies a substantial fraction of a tidal or river channel cross-sectional area, the confinement associated with the channel boundaries substantially increases turbine efficiency. This effect can be exploited to substantially reduce the levelized cost of energy for current turbine arrays. However, the presence of such an array will affect inflow velocity, water depth, and circulation in the surrounding area. This has implications for array performance and environmental impacts. This project combines existing numerical modeling capability from PNNL with an understanding of confined array performance from UW. The work will 1) Simulate tidal currents and water surface elevation in a regional circulation model that includes the Knik Arm site at appropriate resolution; 2) Incorporate a varying thrust coefficient into this model that reflects the outcomes of tidal hydrodynamic simulations of confined turbine arrays; 3) Use the model to determine the effects of array operation on water levels and inflow conditions as a function of different blockage ratios and control strategies; and 4) Identify the “far field” environmental implications of operating an array of this type. UW will use the outputs from this model to evaluate the implications of “feedback effects” and variable inflow conditions on levelized cost of energy estimates for a confinement exploiting array of cross-flow turbines.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Cyber-Physical Reformer Development at NETL

Integrated energy systems are considered one of promising technologies to provide efficient, reliable and resilient power generation. The U.S. Department of Energy, National Energy Technology Laboratory is building an automated reformer-solid oxide fuel cell-gas turbine integrated energy system using a cyber-physical systems (CPS) approach, which exploits the advantage of both numerical model and physical system, as well as gapping the inaccessible technologies. Both the fuel reformer and the fuel cell stack are designed to be CPS components, of which the hardware representations are physically integrated with the gas turbine. A compact design is used with the consideration of future commercialization by placing all components inside of a pressure chamber, which is pressurized by the compressor discharge. (Presented virtually at the MILLENNIUM CLEAN and SUSTAINABLE POWER Workshop 2025)

cyber-physical↗
DOE OSTI2025

System Analysis of an Internal Combustion Engine (ICE) - Solid Oxide Fuel Cell (SOFC) Hybrid Cycle

The variability of renewable energy sources poses challenges for reliable grid operation. Conventional thermal power sources, though reliable, often lack operational flexibility. Hybrid energy systems that integrate Solid Oxide Fuel cells (SOFC) with Internal Combustion Engines (ICE) offer a promising solution by achieving relatively higher efficiency and grid-following capability. This study investigates performance of a 100-kW pressurized SOFC-ICE hybrid cycle. In this configuration, unutilized SOFC fuel is used to drive the engine, with a turbocharger providing air supply and an external reformer generating syngas. System components were numerically modeled using MATLAB/SIMULINK for the SOFC and reformer, and EBSILON® for the ICE and balance of plant. Parametric studies varied fuel utilization (70-90%), reformer temperature (600 – 1000K), anode off-gas recirculation (0 – 70%), and current density (0.2 – 0.55 A/cm2). Results show that the SOFC and ICE operate as thermally independent topping and bottoming cycles, achieving peak efficiency of 62% under optimized conditions.

hybrid↗
DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Evaluation of 1/100-Scale Mooring Systems for Wave Energy Converters

This study evaluates the feasibility, accuracy, and limitations of using 1/100th scale physical mooring systems to represent full-scale mooring behavior for wave energy converters (WECs) during small-scale tank testing. The work focuses on an RM3-style point absorber deployed in conditions representative of the PacWave South test site and examines whether small-scale physical testing can reliably inform numerical modeling, design decisions, and future prototype development. Overall, this study concludes that small-scale (1/100th) physical mooring models can provide valuable qualitative insights, including relative comparisons between mooring types and trends in device behavior, but cannot reliably replicate full-scale mooring loads or dynamic response without significant scaling distortion. Physical tests at this scale are most appropriate for motion characterization and model validation within known limitations, not for deriving absolute mooring loads or final engineering design values. Key lessons learned emphasize the importance of improved instrumentation strategies, iterative wave tank tuning, more precisely manufactured scaled mooring components, and enhanced anchoring systems to reduce experimental uncertainty.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Final report: Insights from ARM observations into aerosol processing and transport by extratropical cyclones and aerosol effects on cyclone clouds

This project advanced understanding of aerosol–cloud interactions in extratropical cyclones using observations from the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) observatory and ACE-ENA field campaigns, along with satellite datasets and numerical modeling. We have quantified the susceptibility of frontal clouds to aerosols, showing that this susceptibility is similar to that in non-frontal clouds despite substantial (factor 2) differences in droplet concentrations between these cloud types, and also significant differences in cloud albedo. A key finding was that Aitken-mode (here, 60-100nm diameter) aerosols play a disproportionately important role in droplet activation in frontal clouds, and concentrations of these particles are frequently strongly biased in climate models. Activation of aerosols depends strongly on updraft speeds, and via collaboration with Argonne National Laboratory we contributed to quantifying these better with the ARM radar wind profiler at the ENA site. Our studies in the North Atlantic have analogs in the Southern Ocean, and to draw these out we examined data from the MARCUS and SOCRATES field campaigns. We identified new particle formation within extratropical cyclones as a potentially important source of cloud condensation nuclei in these pristine marine environments. By combining ARM observations with perturbed parameter ensemble simulations, we demonstrated that surface observations can effectively constrain aerosol–cloud adjustments and reduce uncertainty in simulated cloud liquid water path responses by approximately 15%, yielding stronger constraints on historical aerosol cooling effects. The work further established precipitation processes as a dominant control on aerosol–cloud interactions and Earth system predictability.

Gordon, Hamish [Department of Chemical Engineering↗
DOE OSTI2026

Measurements of TRACER pre-convective conditions and mesoscale circulations using small unmanned aircraft systems (sUAS)

Improved comprehension of the physical processes governing convective cloud formation and lifecycle are of critical importance for understanding and predicting future climate states. The influence of these clouds on the planetary energy budget, including on precipitation, is significant. Things are particularly complex in coastal regimes, where gradients in aerosol particle properties, localized circulations such as sea breezes, and large population centers are found. To date, numerical models struggle to accurately represent these critical clouds and are therefore challenged to provide a realistic view on the planetary energy budget. Through the proposed research, we deployed two uncrewed aircraft systems (UAS) equipped with a variety of instruments alongside sensors deployed by the US Department of Energy Atmospheric Radiation Measurement (ARM) program for the TRACER (Tracking Aerosol Convection Interactions Experiment) field campaign. The two small UAS platforms consisted of a CU RAAVEN fixed-wing airplane and an OU CopterSonde system, with the copter collecting frequent vertical profiles of thermodynamic and kinematic variables such as temperature, pressure, wind and humidity. At the same time, the fixed-wing captured horizontal gradients of these quantities and aerosol size distribution. These systems were deployed south of the Houston metro area, in an area that is impacted by the Gulf of Mexico sea breeze on a daily basis. These observations offer enhanced and complementary perspectives to those provided by the DOE ARM Mobile Facility (AMF) which is was deployed in southeast Houston, and an ancillary site in a more rural location west of the urban Houston area. Quality-controlled versions of the UAS data were collected and posted on the DOE ARM data archive after the conclusion of the campaign where they are accessible by the research community and general public. The UAS perspective offers revolutionary insight into key spatial and temporal effects that have not been evaluated previously.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Fast and Thermal Neutron Spectrum Dosimetry Measurements in the Advanced Test Reactor large-B and small-I Positions Following the Sixth Core Internals Change-out [Slides]

The Advanced Test Reactor (ATR) has a wide variety of irradiation positions that have had experiments that were developed by users from around the world. Most experiment irradiations rely on thoroughly benchmarked numerical models. However, some irradiation positions in ATR are not as well- characterized and are complicated by spectral perturbations from control cylinder orientation. The “small I” and “large-B” irradiation positions are located nearby control cylinders and suffer from these flux perturbations from control cylinder orientations that change during an irradiation cycle to maintain the desired core power distribution. Models predicted that these types of position would exhibit both spectral shifts and amplitude changes in neutron flux, but few measurements have been conducted to benchmark these predictions. Recently, requalification testing was performed to confirm the operational readiness of the ATR following the completion of the Core Internals Change-out (CIC). These tests provided a unique opportunity to validate the analytical methods used to simulate the ATR because nearly all components in the reactor were in a clean as-built state, significantly reducing modelling uncertainties. One subset of the testing included characterization of the fast and thermal neutron flux in the “small-I” and “large-B” positions using silver, cobalt, and nickel neutron dosimetry. In contrast to typical irradiation cycles, the control cylinders were held in position during the post-CIC nuclear testing. The specific activity of these dosimeter wires was measured following two low-power tests with different control cylinder positions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2024

Supporting Greater Interactivity in the IPython Visualization Ecosystem

Interactive visualizations are invaluable tools for building intuition and supporting rapid exploration of datasets and models. Numerous libraries in Python support interactivity, and workflows that combine Jupyter and IPyWidgets in particular make it straightforward to build data analysis tools on the fly. However, the field is missing the ability to arbitrarily overlay widgets and plots on top of others to support more flexible details-on-demand techniques. This work discusses some limitations of the base IPyWidgets library, explains the benefits of IPyVuetify and how it addresses these limitations, and finally presents a new open-source solution that builds on IPyVuetify to provide easily integrated widget overlays in Jupyter.

Martindale, Nathan↗