Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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 667 records · Page 37

Convex Relaxation of Grid-Connected Energy Storage System Models With Complementarity Constraints in DC OPF

Including complementarity constraints in energy storage system (ESS) models in optimization problems ensure an optimal solution will not produce a physically unrealizable control strategy where there is simultaneous charging and discharging. However, the current approaches to impose complementarity constraints require the use of non-convex optimization methods. Here, we propose a convex relaxation for a common ESS model that has terms for both charging and discharging based on a penalty reformulation for use in a model predictive control (MPC) based optimal power flow (DC OPF) problem. In this approach, the complementarity constraints are omitted and a penalty term is added to the optimization objective function. For the DC OPF problem, we provide analysis for the conditions under which the convex relaxation of the complementarity constraint ensures that a solution with simultaneous ESS charging and discharging operation is suboptimal. Simulation results demonstrating ESS behavior with and without the penalty reformulation are provided for an MPC-based DC OPF problem on multiple IEEE test systems.

25 ENERGY STORAGE↗

Surface Method Raytracer

Matlab based experiments that implement multiple ray tracing methods for deterministic radiation transport, including the Surface Method as described in the PhD Thesis, "CONSERVATIVE FIRSTCOLLISION SOURCE TREATMENT FOR RAY EFFECT MITIGATION IN DISCRETE-ORDINATE RADIATION TRANSPORT SOLUTIONS" by Alex Christensen.

Christensen, AlexB↗

Model Inputs, Outputs, and Scripts associated with: “Spatial microbial respiration variations in the hyporheic zones within the Columbia River Basin”

This data package is associated with the publication “Spatial microbial respiration variations in the hyporheic zones within the Columbia River Basin” published in the Journal of Geophysical Research: Biogeosciences (Son et al. 2022) available at doi: 10.1029/2021JG006654. This data package includes the key model inputs/outputs of the river corridor model for the Columbia River Basin (CRB) and the model source codes, which were used in the manuscript. The model is a carbon-nitrogen-coupled river corridor model (RCM), and the model is used to quantify hyporheic zone (HZ) aerobic and anaerobic respiration at the NHDPLUS stream reach scales. The RCM used in this study combines empirical substrate models derived from observations and three microbially driven reactions to compute respiration of the HZ for each National Hydrography Dataset (NHD) reach within the CRB. The reactions in HZs of each NHD reach include anaerobic respiration and two-step anaerobic respiration via denitrification. Our HZ respiration estimates are limited to the lotic (or flowing) stream/river systems, and do not account for the respiration process in water column. Note that the RCM only simulates the HZ’s contribution to the dissolved carbon dioxide (CO2) concentrations in the streams, and the CO2 emissions to the atmosphere are not modelled. The model computes at hourly timesteps because of the fast reaction rates. The key input data of the model are exchange flux, residence time, and stream solute (dissolved organic carbon (DOC), dissolved oxygen (DO), and nitrate concentrations). These inputs are constant over time and represent long-term averaged values.This modeling framework successfully quantified HZ respiration components over multiple scales. It revealed key mechanisms driving the spatial variation of HZ aerobic and anaerobic respiration in reaches with varying hydrologic and substrate conditions. Thus, this modeling study offers a testing hypothesis in different river system (e.g., climate and biomes) for the HZ respiration processes, and can be used as a sampling design tool for large-scale HZ experimental studies.This dataset contains five folders: (1) model_inputs, (2) model_outputs, (3) Rscripts, (4) figures, and (5) model_codes. It also contains a readme, file level metadata (FLMD), and data dictionary (dd). Please see the FLMD for a list of all the files contained in this data package and descriptions for each. The model_inputs folder contains the model inputs used to drive the model simulations. The model_outputs folder contains key model output files from the river corridor model. The Rscripts folder contains the Rscripts for pre- and post- processing model results. The figures folder contains the raw figures associated with the manuscript. The model_codes folder includes key model source codes/input files. All files are .jpg, .jpeg, .out, .e, .od, .dat, .sub, .F90, .0, .R, .sbx, .cpg, .sbn, .shx, .shp, .dbf, .prj, .tfw, .tif, .xml, .pdf, or .csv.

54 ENVIRONMENTAL SCIENCES↗

NMDQi Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST)

The Thermal-hydraulic Optimization, Analysis, and Scoping Tool (TOAST) is developed to support the thermal-hydraulic aspects for the design of optimal irradiation vehicles for the Nuclear Materials Discovery and Qualification initiative (NMDQi). TOAST is currently a MATLAB-based tool that utilizes a simplified steady-state analytical model where a radial thermal resistance circuit is utilized to account for conduction through up to 2 capsules, a gas gap, a thermal bond, a sample cladding, and a sample, as well as the convection from a water coolant outside the capsules. The geometry is discretized axially along the user-defined height of a basket with a user-defined number of points/nodes, essentially turning TOAST into a semi-2D heat transfer analysis utility for cylindrical irradiation vehicles with practically any configuration, solving for temperatures radially at multiple axial locations. TOAST utilizes a Graphical User Interface (GUI) in-which a user can define the geometrical layout, the materials utilized for each component, the coolant characteristics, and the required solution. The user can define the geometry by inputting the diameters, thicknesses, and heights for each component, whereas the materials are defined via the user-provided constant or variable thermal conductivities. The user can select the coolant's inlet temperature, pressure, and the pressure drop across the height of the problem (which is used to calculate the velocity of the flow). Finally, the user can choose to solve for a sample heat generation rate limit by inputting temperature limits for the sample and capsules or can choose to purely solve for the axial temperature distributions of each component by inputting a pre-selected heat generation rate. Either way, TOAST provides the user with axial temperature distributions of all the components including the coolant, and the heat generation rate limit as well as the maximum outlet coolant temperature. The user can also choose to do one of 6 sensitivity analyses in TOAST. The sensitivity analyses yield plots of the sensitivity of the heat generation rate limit, the maximum coolant temperature, and the pressure limit for the annular components due to perturbations in 1-2 unknown variables based on the selected sensitivity analysis. Benchmarks between computations in TOAST and equivalent 2D axis symmetric ABAQUS models are presented, showing that TOAST results are within less than 3% of ABAQUS results in most cases, and a maximum of 8% difference in some cases. The benchmarks also revealed that this uncertainty is tied to the selection of an appropriate Nusselt number correlation and appropriate thermal conductivities, which the user can do from the GUI. Regardless, TOAST is demonstrated as a computationally efficient, highly accessible, and accurate utility for optimization and scoping s of studies different irradiation vehicles.

42 ENGINEERING↗

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗

Whole-Voltage-Range Solid-Solution Reaction in Layered Oxide Cathode of Sodium-Ion Batteries

Layered manganese-based oxides (LMOs) are promising cathode materials for sodium-ion batteries (SIBs) due to their versatile structures. However, the Jahn–Teller effect of Mn 3+ induces severe distortion of MnO 6 octahedra, and the resultant low symmetry is responsible for the gliding of MnO 2 layers and then inferior multiple-phase transitions upon Na + extraction/insertion. Here, hexagonal P2-Na 0.643 Li 0.078 Mn 0.827 Ti 0.095 O 2 is synthesized through the incorporation of Li and Ti into the distorted orthorhombic P'2-Na 0.67 MnO 2 to function as a phase-transition-free oxide cathode. Further, it is revealed that Li in both the transition-metal and Na layers enhances the covalency of Mn–O bonds and allows degeneracy of Mn 3d e g orbitals to favor the formation of hexagonal phase, and the high strength of Ti–O bonds reduces the electrostatic interaction between Na and O for suppressed Na + /vacancy rearrangements. These collectively lead to a whole-voltage-range solid-solution reaction between 1.8 and 4.3 V with a small volume variation of 1.49%. This rewards its excellent cycling stability (capacity retention of 90% after 500 cycles) and rate capability (89 mAh g -1 at 2000 mA g -1 ).

25 ENERGY STORAGE↗

NeRVI: Compressive neural representation of visualization images for communicating volume visualization results

We present NeRVI, a new deep-learning approach that compresses a large collection of visualization images generated from time-varying data for communicating volume visualization results. Based on an image-based implicit neural representation, our approach represents tens of thousands of high-resolution rendering images parametrized by different parameters via a hybrid model of multilayer perceptrons and convolutional neural networks. Here, our model predicts images and corresponding masks, and the masks are utilized for loss computation and network training to capture fine structural details and small components. In conjunction with model quantization and weight encoding, NeRVI yields highly compact compressive neural representations while preserving the image fidelity well. We demonstrate the effectiveness of NeRVI with isosurface rendering and direct volume rendering images generated from multiple data sets and compare NeRVI with other state-of-the-art deep learning-based (InSituNet, SIREN, NeRF, and NeRV) methods. Quantitative and qualitative results show that NeRVI provides an alternative solution that augments domain scientists' ability to manage, represent, and communicate scientific visualization output.

97 MATHEMATICS AND COMPUTING↗

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗

A geospatial model for the analysis of time-dependent land subsidence induced by reservoir depletion

Land subsidence due to fluid depletion is an outcome of physical processes operating across a wide range of time and length scales. Although geomechanical models are crucial to simulate reservoir compaction and predict its long-term fate, their use across large regions often bears prohibitive computational costs. To overcome this obstacle, this paper proposes a simplified modelling framework consisting of (i) a near-field numerical solver simulating the coupling between fluid flow and rock deformation with reference to a simplified one-dimensional geometry and (ii) a far-field geospatial algorithm mapping ground settlements across a region through the superposition of poroelastic computations at multiple wells. The model computes the delay between depletion history and reservoir compaction in proximity of a producing well by assuming basal depletion of a fluid-saturated deformable disk, while the Geertsma solution of nucleus of strain is used to extrapolate the impact of such time-varying reservoir compaction around the well. This approach has been used to back-analyze the spatio-temporal progression of subsidence at the Groningen gas field. The results are compared against measurements collected over 50 year of production at 25 benchmark locations scattered over an area of 900 km 2 . It is shown that coupled simulations based on average values of rock compressibility and permeability lead to nonlinear trends of subsidence evolution in good agreement with field measurements, while uncoupled analyses overpredict settlements by more than 70%. Lastly, synthetic forecasts based on different rates of depletion were provided. The results suggest that slower depletion rates lead to lower subsidence at a given time, and that residual subsidence may continue to develop for several decades after interruption of production activities.

58 GEOSCIENCES↗

Encapsulating Semiconductor Quantum Dots in Supramolecular Cages Enables Ultrafast Guest–Host Electron and Vibrational Energy Transfer

In the field of supramolecular chemistry, host- guest systems have been extensively explored to encapsulate a wide range of substrates, owing to emerging functionalities in nano-confined space that cannot be achieved in dilute solutions. However, host-guest chemistry is still limited to encapsulation of small guests. Herein, we construct a water-soluble metallo-supramolecular hexagonal prism with a large hydrophobic cavity by anchoring multiple polyethylene glycol chains onto the building blocks. Then, assembled prisms are able to encapsulate quantum dots (QDs) with diameters of less than 5.0 nm. Furthermore, we find that the supramolecular cage around each QD strongly modifies the photophysics of the QD by universally increasing the rates of QD relaxation processes via ultrafast electron and vibrational energy transfer. In conclusion, taken together, these efforts expand the scope of substrates in host-guest systems and provide a new approach to tune the optical properties of QDs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Short-lived metal-centered excited state initiates iron-methionine photodissociation in ferrous cytochrome c

Abstract The dynamics of photodissociation and recombination in heme proteins represent an archetypical photochemical reaction widely used to understand the interplay between chemical dynamics and reaction environment. We report a study of the photodissociation mechanism for the Fe(II)-S bond between the heme iron and methionine sulfur of ferrous cytochrome c . This bond dissociation is an essential step in the conversion of cytochrome c from an electron transfer protein to a peroxidase enzyme. We use ultrafast X-ray solution scattering to follow the dynamics of Fe(II)-S bond dissociation and 1 s 3 p (Kβ) X-ray emission spectroscopy to follow the dynamics of the iron charge and spin multiplicity during bond dissociation. From these measurements, we conclude that the formation of a triplet metal-centered excited state with anti-bonding Fe(II)-S interactions triggers the bond dissociation and precedes the formation of the metastable Fe high-spin quintet state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigating the electrical transport properties and electronic structure of Zr 2 ⁢CuSb 3

The checkerboard lattice has been proposed to host topological flat bands as a result of destructive interference among its various electronic hopping terms. However, it has proven challenging to realize experimentally due to the difficulty of isolating this structure from any significant out-of-plane bonding while maintaining structural integrity. Here, in this study, single crystals of Zr 2 CuSb 3 , a potential candidate for the checkerboard lattice, were synthesized using the solution (self-flux) method, and their structure was confirmed via x-ray diffraction. Electrical-transport measurements indicate metallic behavior with electron-dominated carriers. Angle-resolved photoemission spectroscopy reveals multiple electron pockets and significant k z broadening due to its large c axis and low dispersion features in k z . Density-functional theory (DFT) calculations further disentangle the contributions from each high-symmetry plane, providing a comprehensive characterization of electronic behavior. The DFT calculations were then used to determine the orbital contributions of the bands and detect the out-of-plane bonding which prevented the flat bands from forming.

Downey, Eoghan [Univ. of Michigan, Ann Arbor, MI (↗

Towards robust laser beam propagation in atmospheric turbulence

High-fidelity optical propagation through the atmosphere is essential for free-space optical technologies, including laser-based remote sensing and optical communication. However, atmospheric turbulence severely distorts beams and compromises system performance. In this work, we employ hypergeometric-Gaussian (HyGG) vortex beams as probes to characterize and mitigate atmospheric turbulence. Using over 250,000 experimental and simulated frames, we show that refining the power spectrum density (PSD) can reduce numerical prediction errors by up to 79.8%. Concurrently, experimental observations supported by numerical simulations demonstrate that HyGG beams exhibit superior turbulence resilience across multiple metrics compared to conventional Gaussian beams, particularly in their ability to withstand over 5 times stronger turbulence while maintaining similar intensity fluctuations. These dual investigations, on both turbulence mitigation and robust beam solutions, converge to form a unified strategy for enhancing free-space optical system performance. Collectively, our findings provide new insights into light–turbulence interactions and highlight the practical utility of vortex beams under atmospheric conditions.

Zhang, Boyu↗

Idaho Falls Power Black Start Field Demonstration (Preliminary Outcomes Paper)

This April 2021 field demonstration builds upon a 2017 field demonstration in which it was determined IFP’s HPPs, on their own, can support islanded black start and operation up to 2.5 MW loading. Modeling and hardware-in-the-loop testing was used in the intervening period to design an energy storage solution, specifically using ultracapacitors, to reduce likelihood of generators tripping during the field demonstration. Overall, the 2021 field demonstration tested three different options for meeting IFP’s requirements: innovating the hydropower controls, synchronizing multiple HPPs on the system, and integrating an ultracapacitor energy storage system. This report documents the testing performed. Follow-on analysis will provide additional insights, for example, including a complete table of comparative results between the tests and scenarios. The analysis will also propose a refined design for an energy storage system to meet IFP’s grid islanded needs.

13 HYDRO ENERGY↗

Gold Nanoparticles AuNP Decorated on Fused Graphene-like Materials for Application in a Hydrogen Generation

The search for a sustainable, alternative fuel source to replace fossil fuels has led to an increased interest in hydrogen fuel. This combustible gas is not only clean-burning but can readily be produced via the hydrolysis of sodium borohydride. The main drawback of this reaction is that the reaction occurs relatively slowly and requires a catalyst to improve efficiency. This study explored a novel composite material made by combining gold nanoparticles and fused graphene-like materials (AuFGLM) as a catalyst for generating hydrogen via sodium borohydride. The novel fused graphene-like material (FGLM) was made with a sustainable dextrose solution and by using a pressure-processing method. Imaging techniques showed that FGLM appears to be an effective support template for nanoparticles. Transmission electron microscopy (TEM), scanning electron microscopy (SEM), energy-dispersive X-ray spectroscopy (EDS), Fourier-transform infrared spectroscopy (FTIR), X-ray diffraction (XRD) and Raman spectroscopy were used to characterize and determine the size, shape, and structure of nanoparticles and composites. The TEM study characterized the fused carbon backbone as it began to take on a rounder shape. The TEM images also revealed that the average diameter of the gold nanoparticle was roughly 23 nm. The FTIR study confirmed O-H, C-C, and C=O as functional groups in the materials. The EDS analysis showed that the composite contained approximately 6.3% gold by weight. The crystal structures of FGLM and AuFGLM were identified via P-XRD analysis. Various reaction conditions were used to test the catalytic ability of AuFGLM, including various solution pHs, temperatures, and doses of NaBH 4 . It was observed that optimal reaction conditions included high temperature, an acidic solution pH, and a higher dose of NaBH 4 . The activation energy of the reaction was determined to be 45.5 kJ mol -1 , and it was found that the catalyst could be used multiple times in a row with an increased volume of hydrogen produced in ensuing trials. The activation energy of this novel catalyst is competitive compared to similar catalysts and its ability to produce hydrogen over multiple uses makes the material an exciting choice for catalyzing the hydrolysis of NaBH 4 for use as a hydrogen fuel source.

08 HYDROGEN↗

Irradiation-Driven Restructuring of UO 2 Thin Films: Amorphization and Crystallization

Combustion synthesis in uranyl nitrate–acetylacetone–2-methoxyethanol solutions was used to deposit thin UO 2 films on aluminum substrates to investigate the irradiation-induced restructuring processes. Thermal analysis revealed that the combustion reactions in these solutions are initiated at ~160 °C. The heat released during the process and the subsequent brief annealing at 400 °C allow the deposition of polycrystalline films with 5–10 nm UO 2 grains. The use of multiple deposition cycles enables tuning of the film thicknesses in the 35–260 nm range. Irradiation with Ar 2+ ions (1.7 MeV energy and a fluence of up to 1 × 10 17 ions/cm 2 ) is utilized to generate a uniform distribution of atomic displacements within the films. X-ray fluorescence (XRF) and alpha-particle emission spectroscopy showed that the films were stable under irradiation and did not undergo sputtering degradation. X-ray photoelectron spectroscopy (XPS) showed that the stoichiometry and uranium ionic concentrations remain stable during irradiation. The high-resolution electron microscopy imaging and electron diffraction analysis demonstrated that at the early stages of irradiation (below 1 × 10 16 ion/cm 2 ) UO 2 films show complete amorphization and beam-induced densification (sintering), resulting in a pore-free disordered film. Prolonged irradiation (5 × 10 16 ion/cm 2 ) is shown to trigger a crystallization process at the surface of the films that moves toward the UO 2 /Al interface, converting the entire amorphous material into a highly crystalline film. This work reports on an entirely different radiation-induced restructuring of the nanoscale UO 2 compared to the coarse-grained counterpart. The preparation of thin UO 2 films deposited on Al substrates fills an area of national need within the stockpile stewardship program of the National Nuclear Security Administration and fundamental research with actinides. Here, the method reported in this work produces pure, robust, and uniform thin-film actinide targets for nuclear science measurements

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radioactive Cleaning Robot

Here at Fermi National Accelerator Laboratory, particle accelerators are constantly running, giving off excess radioactive dust. This dust can be extremely harmful to people and the environment, posing numerous issues. The goal of this project is to conceptualize an automated cleaning solution. With this main goal in mind, the project started with a brainstorming session. The main concept relied on the use of a powder distributing system paired with a vacuum to easily attract and trap the dust. From this main concept, multiple sketches were produced, showing a robot with pieced-together systems. With the initial idea down, each subsystem could be brought further with the inner workings being sketched out. After sketches and proof of concept were polished, the next step was to begin 3D modeling, using NX, to bring life into the sketches. With every day came new ideas, modifying the previous sketches and models as necessary. Soon after the bigger structures were fleshed out, the focus could be shifted onto the less pressing subsystems. Adding a conveyor belt and spinning brush attachment to the powder hopper allowed for more even powder distribution. Along with that, the addition of spring-loaded brushes on the vacuum nozzle ensures constant dust agitation, improving overall functionality. A few other quality of life improvements, such as a removable dust reservoir with a ramp gave a push towards real life functionality. After considering the logistics of this robot functioning in real-life, pressure plates were provided to stop the robot from crashing into its surroundings. Having all of the individual models lead to an overall assembly. Making this robot programmable allowed for this project to meet the design requirements, providing a solution to the previous dilemma.

43 PARTICLE ACCELERATORS↗