Engineering PapersSearch

SEARCH · Engineering Papers

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.

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

At least 73 records · Page 4

Modeling betatron radiation using particle-in-cell codes for plasma wakefield accelerator diagnostics

The analysis of plasma wakefield acceleration experimental measurements, particularly in the characterization of photons emitted through the betatron radiation mechanism, requires the development of accurate numerical models. These computational models are crucial for supporting modern instrumentation designed to measure the single-shot, double-differential angular-energy radiation spectra resulting from interactions between beams and plasmas. Motivated by the needs of such applications, this paper presents detailed numerical models of betatron radiation generated in beam-plasma acceleration experiments. These models are based on the integration of the Liénard-Wiechert (LW) potentials, applied to computed particle trajectories. The particle trajectories are generated using three distinct methods: first, by tracking particles through idealized fields in the blowout regime of PWFA; second, by obtaining trajectories using the fast quasistatic particle-in-cell (PIC) code quickpic; and third, obtaining trajectories from the fully self-consistent PIC code osiris. To ensure the accuracy and reliability of these models, the paper includes various benchmark tests using analytical expressions, as well as employing the PIC code epoch, which takes an alternative approach by using a Monte Carlo quantum electrodynamics (QED)-based radiation model. Additionally, the paper presents simulations of the expected experimental betatron radiation spectra, taking into account parameters relevant to PWFA and plasma photocathode experiments at the SLAC FACET-II facility.

Yadav, M. [University of California, Los Angeles,

Utah FORGE: LBNL Reports on VEMP Electromagnetic Data Collection and Processing - 2024

This archive contains reports related to Vertical Electromagnetic Profiling (VEMP) tool data collection and processing at Utah FORGE in 2024. The first report describes LBNL's effort to collect electromagnetic geophysical data with the tool in well 78-32B and a downhole electrode in well 16A. The second report describes the final data acquisition and processing of the VEMP electromagnetic data collected at the Utah FORGE site in May of 2024. Also included are a noise analysis as well as a comparison of the data to numerical models. This was originally presented as a paper at the 2025 Stanford Geothermal Workshop.

15 GEOTHERMAL ENERGY

Analysis on Evaluations of Monterey Bay Aquarium Research Institute’s Wave Energy Converter’s Field Data Using WEC-Sim and Gazebo: A Simulation Tool Comparison

Although many studies have validated wave energy converter (WEC) numerical models against scaled prototype experimental data, there remains a notable lack of validation using data from full-scale deployed WECs. This paper compares two numerical models of Monterey Bay Aquarium Research Institute’s Wave Energy Converter (MBARI-WEC), a two-body point absorber with an electro-hydraulic power take-off system (PTO). The models are implemented in WEC-Sim/Simscape and Gazebo Simulator. A statistical analysis of the models was performed, and field results were obtained to compare the models’ accuracy in predicting the RMS piston velocity, RMS motor speed, and mean electric power compared to field data for 56 observations across varying sea states. The Gazebo model demonstrated a closer agreement across all three parameters for a majority of the observations. When compared to the field data, the Gazebo and WEC-Sim models exhibited average mean electric power overestimations of 13% and 22%, respectively.

16 TIDAL AND WAVE POWER

Data-driven global ocean modeling for seasonal to decadal prediction

Accurate modeling of ocean dynamics is crucial for enhancing our understanding of complex ocean circulation processes, predicting climate variability, and tackling challenges posed by climate change. Although great efforts have been made to improve traditional numerical models, predicting global ocean variability over multiyear scales remains challenging. Here, we propose ORCA-DL (Oceanic Reliable foreCAst via Deep Learning), a data-driven three-dimensional ocean model for seasonal to decadal prediction of global ocean dynamics. ORCA-DL accurately simulates the three-dimensional structure of global ocean dynamics with high physical consistency and outperforms state-of-the-art numerical models in capturing extreme events, including El Niño–Southern Oscillation and upper ocean heat waves. Moreover, ORCA-DL stably emulates ocean dynamics at decadal timescales, demonstrating its potential even for skillful decadal predictions and climate projections. Our results demonstrate the high potential of data-driven models for providing efficient and accurate global ocean modeling and prediction.

Science & Technology - Other Topics

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS

Modeling the influence of the solid electrolyte interphase on the sand’s time and dendrite formation on lithium metal electrodes

Lithium metal is a sought after battery material for its high energy density due to the low electrochemical potential and density. However, lithium metal is also highly reactive, which results in a strong propensity for dendrite formation. The Sand’s time has previously been used to predict the time of dendrite initiation on metals that do not form a solid-electrolyte interphase (SEI), but it has been shown that the Sand’s time is not accurate for lithium electrodes when using transport parameters associated with the electrolyte. Thus, we built a numerical model to simulate lithium ion transport through a growing SEI to predict the Sand’s time. The numerical model is shown to be more accurate than previous analytical solutions, especially for low current densities. We then analyze the sensitivity of the Sand’s time to different SEI properties and the chemical potential gradients present in the SEI, driving lithium transport. The results showed that high lithium concentration has a greater impact at high current density, while fast diffusivity is more important at low current density. Lastly, we modeled the influence of surface roughness on the plating evolution and chemical potential gradients when an SEI is present in comparison to the electrolyte. As a result, we demonstrate that the SEI plays a critical role in lithium electrode stability, and that improved characterization techniques are needed to better understand transport through the SEI and increase lithium metal utilization in energy storage devices.

Chemistry

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham

Demonstration of gold nanorod systems for enhanced total efficiency: Experimental and numerical analysis

This study investigates the photothermal performance of gold nanorods engineered to exhibit longitudinal plasmon resonances at 695 nm, 780 nm, and 970 nm. The work combines synthesis, structural characterization, extinction measurements, numerical modeling, and controlled temperature experiments to quantify how nanorod geometry, resonance tuning, concentration, and chamber shape jointly influence heat generation. Transmission electron microscopy confirms that increasing nanorod aspect ratio systematically shifts the longitudinal plasmon peak toward the near-infrared region. Extinction measurements show strong agreement with theoretical predictions, validating the numerical model across two independent datasets. Three chamber geometries were tested under laser excitation at 640 nm, 808 nm, and 980 nm: an ascending stepped base, a flat base, and a descending stepped base. Without nanorods, the ascending geometry produced the highest efficiency due to enhanced natural convection. After introducing gold nanorods, all geometries exhibited substantial thermal enhancement, with total efficiencies exceeding 20%. The strongest improvement was obtained for nanorods resonant at 780 nm with a mass concentration of 4.6 mg/mL implemented on the descending stepped-base geometry. This performance resulted from the combined effect of spectral overlapping with the 808 nm laser, the highest nanorod concentration, and localized heat accumulation that intensified buoyancy-driven flow. The findings demonstrate that total efficiency is governed by a synergistic interplay between optical resonance, nanoparticle concentration, and macroscopic chamber design, revealing the system-level coupling between nanoscale plasmonic absorption and macroscale heat-transfer phenomena. The results provide a validated framework for tuning nanoscale plasmonic absorbers and optimizing thermal systems for applications requiring efficient light-to-heat conversion.

77 NANOSCIENCE AND NANOTECHNOLOGY

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Solar Decarbonization of Paraffin Dehydrogenation Through Particle Heat Carriers (Final Technical Report)

This project focuses on solutions to decarbonize high-temperature catalytic processes using solar thermal heat. The primary project goal is to show the validity of a moving packed bed reactor for propane dehydrogenation using catalyst particles as the heat carrier for the reaction, which can be heated by concentrated solar energy in a particle receiver. This concept, if further developed, may provide a cost-effective pathway for converting lower value gases to important chemical precursors for industrial materials using only renewable energy. The project was divided into six tasks. In Task 1, DFT calculations were performed to understand the role of Pt and Sn in the catalytic dehydrogenation reaction. In Task 2 chemical kinetics measurements were made for several catalyst formulations at high temperatures. In Task 3, the solar absorptance of catalyst particles was compared to the absorptance of commonly used materials in particle receivers. In Task 4, numerical models were developed which could predict performance of the complete system and predict specific temperatures in the system. In Task 5, a prototype system was designed, fabricated, and tested to show the validity of the concept. Task 6 concerned project management activities. Experiments with the prototype showed repeatable thermal performance at temperatures targeted for the reaction. A limited set of tests were done with active catalyst and propane dehydrogenation, showing conversion of propane to propylene with a range of conversions and selectivities. The results are promising, and the prototype designed was reliable during testing, and the team expects that further development of the prototype would yield improved results. A numerical model framework based on coupled fluid and particle mechanics was developed with high computational efficiency using GPU calculations. The model may prove highly useful for evaluating other high-temperature particle systems. However, it was determined that simpler porous media models were good fits for the needs of the current moving packed bed concept. Data showing strong solar absorption of the particles validates the plan of using existing solar particle receivers with only a change in the particle type. Catalyst investigation showed that Pt 1 Sn 3 is the most viable candidate for developing PtSn catalysts for high temperature propane dehydrogenation, considering the balance of activity, selectivity, and deactivation. This project completed an initial study of various factors needed to incorporate a moving bed catalytic reactor for propane dehydrogenation into a concentrated solar thermal particle system. Future developments may allow this technology to be scaled up and help to use solar thermal energy to decarbonize not only the propane dehydrogenation reaction, but other gas-solid catalytic reactions at similar temperatures.

14 SOLAR ENERGY

TEAMER: Crossflow Turbine Fairing Geometry Optimization - Report and CFD Modeling Files

The dataset includes computational fluid dynamics (CFD) models and simulation files for crossflow turbines as well as a detailed project report. The report documents the project undertaken by the Ocean Renewable Power Company (ORPC) to design and optimize a modular fairing for the Modular RivGen Marine Hydrokinetic (MHK) turbine, which enhances the efficient deployment and operation of turbine arrays. The project focused on optimizing the hydrodynamic performance of the fairing using CFD, with an emphasis on two key geometric parameters: the fairing's cross-sectional shape and the spacing between the rotor and the fairing. The analysis aimed to maximize net power output while also assessing discretized loading to evaluate ultimate and fatigue loads on the turbine components. The numerical modeling was conducted using both the commercial CFD software Star-CCM+ and the open-source code openFOAM, with the latter utilizing the actuator line library, turbinesFOAM. This dual-code approach was intended to increase confidence in the results and demonstrate the viability of using open-source tools for high-fidelity marine energy modeling. This dataset includes all necessary files for actuator line simulations in openFOAM, as well as 2D blade-resolved CFD results, along with Python and Java scripts for setting up and post-processing simulations.

16 TIDAL AND WAVE POWER

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES

Air Force Institute of Technology Dust Cloud Model Vertical Dispersion

The Air Force Institute of Technology (AFIT) has developed and refined a nuclear dust cloud model for modeling nuclear fallout, and the dust mass and radiation encountered by aircraft flying through a nuclear dust cloud. Recently, a question was raised requesting an explanation of why the vertical mass and activity distributions had an anomalous shape for yields above 1 Mt. This report explains that the anomalous shape occurs because the dust cloud model is a discrete numeric model, and not a continuous model.

54 ENVIRONMENTAL SCIENCES

DECOVALEX-2023: Task A Final Report

Task A, also known as HGFrac, examines the fracturing processes that may occur in the Callovo-Oxfordian claystone (COx) in the context of the high-level (HLW) and intermediate-level long-lived (ILW-LL) radioactive waste repository in France. Understanding these processes and improving numerical models to reproduce them will aid in the design, optimization, and safety of the repository. Heat and gas fracturing are studied in two independent subtasks following a stepwise approach: laboratory tests/benchmark exercises, in-situ experiment, and finally, an application case. The in-situ heater experiment aimed to thermally induce a hydraulic fracture; temperature and pore pressure were monitored to detect any evidence of fracturing. The in-situ gas injection experiment aimed to study the effect of the stress orientation and gas injection kinetics on the gas fracturing process. The occurrence of fracturing was monitored by gas pressure measured in the injection interval. In both experiments, the excavation-induced fracture network around the heater/injection boreholes played an important role in the reduction of the compressive stress state, leading to both a tensile and shear failure response of the COx. During the first two years of the project, the research teams working on each task developed and/or proposed numerical approaches for reproducing the occurrence of fracturing in the in-situ experiments. The failure criteria were defined by reproducing the measurements from laboratory extension tests for the heat fracturing subtask. In the gas fracturing subtask, their approaches were used for simulating several benchmark exercises, and an inter-comparison between models was carried out. In both tasks, the developed approaches were compared with a simplified approach considering poro-elasticity for the mechanical behaviour of the COx. Most of the developed approaches are based on a continuous medium that takes into account variations in hydraulic properties due to mechanical degradation, such as plastic deformation or damage. Other approaches implicitly modelled weak planes or embedded discontinuities to reproduce fracture propagation. The potential for fracture initiation was also studied through of a discrete approach. In the second half of the project, the research teams mainly focused on interpretative modelling of two in-situ experiments and a blind prediction exercise to test their respective approaches. The models developed by the research teams were also applied at the repository scale to evaluate fracture initiation in a case study under vi unfavourable conditions, particularly in terms of spacing between High-Level Waste cells. The results showed that the poro-elasticity approach could be an efficient tool for understanding the main processes occurring in the COx. One example is the explicit representation of the excavation-induced fracture network around the boreholes, which yielded acceptable results compared to the measurement data. However, advanced approaches were needed to evaluate the potential increase of the excavation-induced fracture network extend and better understand fracture initiation. The stress analyses carried out by the teams revealed that hydraulic boundary conditions had a strong impact on fracture initiation in the heater experiment. Furthermore, in most cases, the results required higher pore pressure increments to reach fracturing than those measured in the experiment. This implies that the measurements may have been biased by the packer’s capacity to fully isolate the piezometric chambers, leading to lower pressures. On the contrary, there was no agreement on the fracturing mode, as some reported either shear or tensile fracturing, while others reported a combination of the two modes. In the case of gas fracturing, the research teams were limited to the comparison of a single point, which complicated their task. Nonetheless, the numerical results were able to reproduce the measurements and capture processes such as longitudinal gas flow through the excavation-induced fracture network, as suggested by some evidence in the observation piezometric chambers. The numerical models also agreed with the measurements in the sense of higher probability of developing along the injection borehole than radially towards the sound rock. The approaches developed by the research teams showed that they are capable of analysing and reproducing fracture initiation in the COx. However, areas of future work should focus on the fracture propagation and fracture aperture, which were out of the scope of this task. To this end, additional data must be gathered for the parameter characterisation and validation of the numerical models. Nonetheless, various approaches showed promising results as they were able to reproduce fracture development under certain conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Modeling and Experimental Demonstration of Flux Spreading in Light Trapping Planar-Cavity Solar Enclosed Particle Receivers

This study experimentally validates and numerically models the flux-spreading effect in a light-trapping planar-cavity solar receiver) for particle-based concentrating solar power systems. The receiver's shallow cavity with vertical planar walls redistributes concentrated solar flux, reducing peak intensity and achieving uniform heat flux. On-sun tests at National Renewable Energy Laboratory's High-Flux Solar Furnace under flux up to 1500 kW/m2 measured cavity wall temperatures, which were compared with Monte Carlo ray-tracing (SolTrace) and computational numerical simulations. Three angular absorptance models were evaluated: constant absorptance, a Pyromark-based directional model, and a Fresnel-based Cr2O3 model. The Fresnel-derived model showed the best agreement with experiments, achieving high correlation (PC > 0.85), structural similarity (SSIM > 0.98), and signal-to-noise ratios (PSNR > 40 dB), with temperature prediction errors of 1-11%. Results confirm that flux spreading mitigates local overheating and validate the integrated modeling approach, supporting the solar receiver scalability for high-efficiency, high-temperature concentrating solar power applications.

14 SOLAR ENERGY

Modeling the formation of Sedan Crater using the FLAG and HOSS codes

Numerical modeling of explosion crater formation requires accounting for complex physical processes. Numerical validation of explosion cratering is an important step in modeling and requires experimental data for comparison. Models using discrete elements and continuum models have both benefits and drawbacks to their approaches. In this work, we consider both an arbitrary Lagrangian–Eulerian (ALE) hydrocode and a finite discrete element method (FDEM) approach to modeling the formation of the Sedan crater, the largest human-made crater in the United States. The Sedan crater formed from an underground nuclear detonation in the Nevada desert as part of Project Plowshare. Our models show that the continuum approach of the hydrocode matched well compared to early test time prior to the mound rupture and subsequent fireball venting, when most of the alluvium exhibited fluid behavior. Our FDEM approach matched the final crater dimensions well, after material had settled back into the crater, when material strength and solid mechanics play key roles. Our work shows how leveraging the benefits of multiple numerical approaches can lead to better understanding of complex physical problems, especially problems with limited experimental data. By using a continuum approach to early-time hydrodynamics and an FDEM approach to later-time solid mechanics, we can better understand the different physical regimes of explosion crater formation.

36 MATERIALS SCIENCE

Characterizing Spatial and Temporal Variability of California’s Coastal Upwelling Using a Satellite-Derived CUTI Index

Ocean upwelling off California is a wind-driven coastal process with complex interactions between spatial wind patterns and sea surface temperatures. Large-scale spatial and temporal upwelling characterizations are traditionally based on numerical models, while quantifying upwelling with direct observations is not straightforward. To complement model-based upwelling characterizations, here we use satellite estimates of winds, sea surface temperatures, and ocean currents along the California coast to calculate the coastal upwelling transport index (CUTI), originally calculated from a numerical ocean model. The CUTI framework accounts for both Ekman and geostrophic offshore-directed ocean transport, the latter of which is often neglected in satellite-based upwelling assessments. The approach yields the CUTI over a 30-yr period at a daily temporal and 0.25° spatial resolution to analyze upwelling patterns, with focus on potential Morro Bay and Humboldt wind energy areas. Morro Bay (CUTI median of 0.63 m 2 s −1 ) exhibits stronger, more consistent upwelling throughout the year, whereas Humboldt (CUTI median of 0.45 m 2 s −1 ) exhibits a pronounced seasonal cycle with strongest upwelling in the summer and occasional downwelling in winter. The cluster analysis of 1993–2022 satellite maps of CUTI identifies four characteristic spatial upwelling patterns along the coast, corresponding to distinct spatial distributions of sea surface temperature anomalies, Ekman transport, and geostrophic ocean water transport. This dataset provides an observation-based foundation for characterizing natural coastal upwelling off the California coast, showing good agreement with the original model CUTI while being based on satellite data and offering increased spatial resolution.

16 TIDAL AND WAVE POWER