Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 235 records · Page 13

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Rapid Identification of X-ray Diffraction Patterns Based on Very Limited Data by Interpretable Convolutional Neural Networks

Large volumes of data from material characterizations call for rapid and automatic data analysis to accelerate materials discovery. Herein, we report a convolutional neural network (CNN) that was trained based on theoretical data and very limited experimental data for fast identification of experimental X-ray diffraction (XRD) patterns of metal–organic frameworks (MOFs). To augment the data for training the model, noise was extracted from experimental data and shuffled; then it was merged with the main peaks that were extracted from theoretical spectra to synthesize new spectra. For the first time, one-to-one material identification was achieved. Theoretical MOFs patterns (1012) were augmented to a whole data set of 72 864 samples. It was then randomly shuffled and split into training (58 292 samples) and validation (14 572 samples) data sets at a ratio of 4:1. For the task of discriminating, the optimized model showed the highest identification accuracy of 96.7% for the top 5 ranking on a test data set of 30 hold-out samples. Neighborhood component analysis (NCA) on the experimental XRD samples shows that the samples from the same material are clustered in groups in the NCA map. Analysis on the class activation maps of the last CNN layer further discloses the mechanism by which the CNN model successfully identifies individual MOFs from the XRD patterns. Furthermore, this CNN model trained by the data augmentation technique would not only open numerous potential applications for identifying XRD patterns for different materials, but also pave avenues to autonomously analyze data by other characterization tools such as FTIR, Raman, and NMR spectroscopies.

36 MATERIALS SCIENCE↗

Watching a signaling protein function: What has been learned over four decades of time-resolved studies of photoactive yellow protein

Photoactive yellow protein (PYP) is a signaling protein whose internal p-coumaric acid chromophore undergoes reversible, light-induced trans-to-cis isomerization, which triggers a sequence of structural changes that ultimately lead to a signaling state. Since its discovery nearly 40 years ago, PYP has attracted much interest and has become one of the most extensively studied proteins found in nature. The method of time-resolved crystallography, pioneered by Keith Moffat, has successfully characterized intermediates in the PYP photocycle at near atomic resolution over 12 decades of time down to the sub-picosecond time scale, allowing one to stitch together a movie and literally watch a protein as it functions. But how close to reality is this movie? To address this question, results from numerous complementary time-resolved techniques including x-ray crystallography, x-ray scattering, and spectroscopy are discussed. Emerging from spectroscopic studies is a general consensus that three time constants are required to model the excited state relaxation, with a highly strained ground-state cis intermediate formed in less than 2.4 ps. Persistent strain drives the sequence of structural transitions that ultimately produce the signaling state. Crystal packing forces produce a restoring force that slows somewhat the rates of interconversion between the intermediates. Moreover, the solvent composition surrounding PYP can influence the number and structures of intermediates as well as the rates at which they interconvert. When chloride is present, the PYP photocycle in a crystal closely tracks that in solution, which suggests the epic movie of the PYP photocycle is indeed based in reality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diagnostics for multiple frequency heating and investigation of underlying processes

The development of new facilities routinely challenges ion source designers to build and operate sources that can achieve ever higher beam intensities and energies. As shown in this work, electron cyclotron resonance ion sources have proven to be extremely capable in meeting these challenges through the production of intense beams of medium and high-charge state ions. As performance boundaries are pushed, source stability becomes an issue as does the technology required to meet the challenge. Multiple frequency heating, the simultaneous use of two or more plasma heating frequencies, is a powerful tool in meeting the simultaneous need of intensity and stability. Relatively straightforward to utilize, the technique has been employed at numerous facilities to increase beam current and achievable charge state while also stabilizing the plasma. Its application has expanded the operational boundaries of existing and next generation sources, demonstrating that these devices have not yet achieved their full operational potential. To better understand the underlying physics, the diagnostics used to probe the source operational boundaries and the plasma properties have become increasingly sophisticated. In concert with detailed modeling, they are beginning to provide insight into the heating mechanism and, with that, the prospect of future advances.

47 OTHER INSTRUMENTATION↗

The cosmological analysis of the SDSS/BOSS data from the Effective Field Theory of Large-Scale Structure

The Effective Field Theory of Large-Scale Structure is a formalism that allows us to predict the clustering of Cosmological Large-Scale Structure in the mildly non-linear regime in an accurate and reliable way. In this paper, after validating our technique against several sets of numerical simulations, we perform the analysis for the cosmological parameters of the DR12 BOSS data. We assume Λ CDM, a fixed value of the baryon/dark-matter ratio, Ω b /Ω c , and of the tilt of the primordial power spectrum, n s , and no significant input from numerical simulations. By using the one-loop power spectrum multipoles, we measure the primordial amplitude of the power spectrum, A s , the abundance of matter, Ω m , and the Hubble parameter, H 0 , {to about 13% , 3.2% and 3.2% respectively, obtaining ln (10 10 A s )=2.72± 0.13 , 0Ω m =0.309± 0.01 , H 0 =68.5± 2.2 km/(s Mpc) at 68% confidence level. If we then add a CMB prior on the sound horizon, the error bar on H 0 is reduced to 1.6% .} These results are a substantial qualitative and quantitative improvement with respect to former analyses, and suggest that the EFTofLSS is a powerful instrument to extract cosmological information from Large-Scale Structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Metal-insulator transition and magnetism of SU(3) fermions in the square lattice

We study the SU(3) symmetric Fermi-Hubbard model (FHM) in the square lattice at 1/3-filling using numerically exact determinant quantum Monte Carlo and numerical linked-cluster expansion techniques. We present the different regimes of the model in the T–U plane, which are characterized by local and short-range correlations, and capture signatures of the metal-insulator transition and magnetic crossovers. These signatures are detected as the temperature scales characterizing the rise of the compressibility, and an interaction-dependent change in the sign of the diagonal spin-spin correlation function. The analysis of the compressibility estimates the location of the metal-insulator quantum critical point at U c /t ~ 6, and provides a temperature scale for observing Mott physics at finite T. Furthermore, from the analysis of the spin-spin correlation function we observe that for U/t ≳ 6 and T ~ J = 4⁢t 2 /U there is a development of a short-range two-sublattice (2SL) antiferromagnetic structure, as well as an emerging three-sublattice (3SL) antiferromagnetic structure as the temperature is lowered below T/J ≲ 0.57. This crossover from 2SL to 3SL magnetic ordering agrees with Heisenberg limit predictions, and has observable effects on the density of on-site pairs. Finally, we describe how the features of the regimes in the T–U plane can be explored with alkaline-earth-like atoms in optical lattices with currently achieved experimental techniques and temperatures. Furthermore, the results discussed in this paper provide a starting point for the exploration of the SU(3) FHM upon doping.

74 ATOMIC AND MOLECULAR PHYSICS↗

Particle-Hole Asymmetric Ferromagnetism and Spin Textures in the Triangular Hubbard-Hofstadter Model

In a lattice model subject to a perpendicular magnetic field, when the lattice constant is comparable to the magnetic length, one enters the “Hofstadter regime,” where continuum Landau levels become fractal magnetic Bloch bands. Strong mixing between bands alters the nature of the resulting quantum phases compared to the continuum limit; lattice potential, magnetic field, and Coulomb interaction must be treated on equal footing. Using determinant quantum Monte Carlo and density matrix renormalization group techniques, we study this regime numerically in the context of the Hubbard-Hofstadter model on a triangular lattice. In the field-filling phase diagram, we find a broad wedge-shaped region of ferromagnetic ground states for filling factor ν ≤ 1 , bounded below by filling factor ν = 1 and bounded above by half filling the lowest Hofstadter subband. We observe signatures of SU(2) quantum Hall ferromagnetism at filling factors ν = 1 and ν = 3 . The phases near ν = 1 are particle-hole asymmetric, and we observe a rapid decrease in ground-state spin polarization consistent with the formation of skyrmions only on the electron doped side. At large fields, above the ferromagnetic wedge, we observe a low-spin metallic region with spin correlations peaked at small momenta. We argue that the phenomenology of this region likely results from exchange interaction mixing fractal Hofstadter subbands. The phase diagram derived beyond the continuum limit points to a rich landscape to explore interaction effects in magnetic Bloch bands. Published by the American Physical Society 2024

Ding, Jixun K.↗

Electromagnetic Transient Simulation of Large-Scale Inverter-Based Resources With High-Granularity

The power grid is undergoing a significant transformation with the rapid increase in inverter-based resources (IBRs), including large-scale photovoltaic (PV) plants. Ensuring reliable and resilient grid operation in this new paradigm necessitates high-granularity electromagnetic transient (EMT) modeling that accurately captures the behavior of individual inverters and their interactions within IBR plants. Central to this approach is the detailed representation of both the IBR plant’s collector system and the dynamics of individual inverters. To achieve this, a high-granularity EMT model of a large-scale PV plant has been developed using advanced simulation algorithms, including matrix splitting and the Schur complement. These proposed techniques significantly enhance simulation speed, numerical stability, and accuracy while improving the modularity and efficiency of the collector system’s representation. The effectiveness of the proposed methods is validated through simulations of a representative large-scale PV plant consisting of 125 individual PV inverters, 25 IBR unit transformers, and a 52-bus collector system.

Choi, Jongchan [Oak Ridge National Laboratory (ORN↗

STEPS: A Portable Numerical Simulation Toolkit for Electrical Power System Dynamic Studies

Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.

42 ENGINEERING↗

Failure Probability Constrained AC Optimal Power Flow

Despite cascading failures being the central cause of blackouts in power transmission systems, existing operational and planning decisions are made largely by ignoring their underlying cascade potential. This paper posits a reliability-aware AC Optimal Power Flow formulation that seeks to design a dispatch point which has a low operator-specified likelihood of triggering a cascade starting from any single component outage. By exploiting a recently developed analytical model of the probability of component failure, our Failure Probability-constrained ACOPF (FP-ACOPF) utilizes the system's expected first failure time as a smoothly tunable and interpretable signature of cascade risk. Here, we use techniques from bilevel optimization and numerical linear algebra to efficiently formulate and solve the FP-ACOPF using off-the-shelf solvers. Extensive simulations on the IEEE 118-bus case show that, when compared to the unconstrained and N-1 security-constrained ACOPF, our probability-constrained dispatch points can significantly lower the probabilities of long severe cascades and of large demand losses, while incurring only minor increases in total generation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Force Reconstruction Using the Inverse of the Mode-Shape Matrix [Book Chapter]

A new technique for force reconstruction is developed. To estimate the externally applied force, this technique sums the weight-scaled acceleration signals, and is referred to as the Sum of Weighted Accelerations Technique (SWAT). To obtain the scalar weights the inverse of the mode shape matrix is used. Application of this technique is illustrated with both numerical calculations using a mass-spring model and experimental data from a structure impacting a rigid barrier.

055001 -- Nuclear Fuels-- Safeguards, Inspection, ↗

Autonomous Aerial Power Plant Inspection in GPS-denied Environments

Inspection of coal-fired power plants is frequently dangerous, includes difficult places to reach, and can turn expensive due to the downtimes and cost of inspection crew. Robotic systems have shown capabilities to address some of these issues, but most of the current robotic inspection technology in power plants is designed for specific components. Conversely, recent advances in machine vision have empowered aerial platforms for long-range, remotely-controlled, GPS-based inspections of industrial plants. This capability has led to wide spread utilization of aerial robots (commonly termed Drones, UVS or UAS) platforms for inspection in less challenging environments where both collision avoidance, and GPS reception are not significant issues. The challenge in adapting airborne technology for power plant inspection lies in internal structures and the complex network of piping, and distribution systems, which impose significant risks for collision and can hinder the reception and transmission of GPS signals. The current state of the art in aerial inspection technology within the energy sector is controlled via radio control, and utilizes GPS-based navigation, for inspection of large-scale plants such as offshore platforms and wind turbine parks. Nevertheless, close-range and autonomous inspection in the GPS-denied environments of power plants has not yet been achieved, as it requires precise guidance and navigation with real-time situational awareness and obstacle avoidance capabilities. This endeavor introduced the use of rotary wing flying robots, due to their station keeping and vertical take-off capabilities for power plant components inspection. To enable close quarter inspection two methods were used. One method uses the 3D CAD (Three-dimensional Computer-Aided Design) model of the asset to inspect to generate the UAV’s inspection path. To acquire, analyze and process the 3D model, first, the STL file is produced to obtain surface points and vectors normal to the surface. Later, by introducing other variables such as wall offset and a controlled trajectory between each outline and each subsequent layer, the flight path is generated. The proposed framework will generate a path that will pass as close as desired from the surface and navigate in intricate environments. A second method, use advanced manufacturing techniques such as CNC (Computer Numerical Control) and additive manufacturing. Once the inspection flight path is obtained, vision-based navigation systems are employed to have the UAV autonomously tracking the provided trajectory. Finally, Artificial Intelligence-enabled developments are in charge of detecting cracks and corrosion in structural components of power plants. The proposed methods are validated in simulations, laboratory and industrial setups, where it is shown that the developed systems acting together enable close-quarter autonomous aerial inspection and mapping in power plant assets. The system can be further improved by adding more sensors to navigate in different GPS-denied environments, with non-homogeneous lighting conditions, dust and in general situations where vision-based systems may fail.

01 COAL, LIGNITE, AND PEAT↗

SVD Orbit Correction in Zgoubi And Examples: SATURNE, RHIC, EIC RCS

This Tech. Note describes the SVD technique in the Lorentz force numerical solver zgoubi. It has been used as a tool to study the effect of vertical orbit perturbation on polarization transmission in RHIC, and for the study of the compensation of orbit defects resulting from stray fields in the EIC RCS.

43 PARTICLE ACCELERATORS↗

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↗

03.02.02.45: Energy Efficiency Improvement Approaches in Ice Related Processes

The primary objective of this project was to facilitate the dislocation of the interfacial ice layer by employing advanced materials and ultrasonic vibration to reduce ice adhesion strength. This work employed two technical approaches that were thoroughly investigated and previously reported. The effectiveness of these approaches, both individually and in combination, has been quantified, demonstrating notable energy savings. The projected payback period for these enhancements is approximately 2.2 years or less, contingent upon specific energy costs. Furthermore, these advancements hold significant promise for reducing carbon emissions across various equipment scales. This study particularly focused on the ultrasonic deicing technique for diverse structures, utilizing numerical simulations to evaluate performance and potential benefits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Community-Informed Urban Flood Modeling for Impact Mitigation

The intensification of the hydrologic cycle due to climate change poses a threat to aging and under-designed water infrastructure systems which cannot adequately manage intense storm events. Developing a comprehensive plan for managing rain-driven flooding events is challenging due to uncertainties in the magnitude and frequency of future storm events and conflicting stakeholder objectives. In the City of Baltimore, Maryland, stormwater infrastructure is struggling to keep up with rainfall-driven (pluvial) flooding events, which regularly damage housing and disrupt transportation for residents. In this study, a hybrid of community engagement, numerical modeling, and artificial intelligence techniques are employed to explore prospective urban flooding adaptations. Community engagement drives the development of an urban flooding model (EPA Storm Water Management Model) for the Baltimore Harbor watershed. The model integrates complex surface and subsurface stormwater infrastructure data from the City, high-resolution spatial data, insights from local public works experts, and the lived experiences of City residents. This co-developed model simulates adaptations of interest to stakeholders in the city, including green and grey infrastructure and operational management strategies. Stormwater management scenarios focused on inlet cleaning and spatially concentrated green infrastructure are found to be the most effective in reducing flood depths in community priority locations. Together, these adaptations can reduce the duration of intersection inundation by more than twenty minutes, allowing for quicker emergency response and restoration of typical transportation systems. Future work will combine this community engaged flooding model with the Deep Uncertainties Pathways framework to explore tradeoffs between adaptations and develop dynamic adaptations which align with community objectives, enhance climate resilience in Baltimore, and can be adjusted in response to changing future conditions.

Ava, Spangler [Pennsylvania State University]↗

Atomic Layer Deposition for Materials-Based H2 Storage: Opportunities and Limitations

The transportation demands in our growing hydrogen (H2) economy requires robust storage systems. The current commercially implemented technology in fuel cell cars relies on the well-established technology of hydrogen gas compressed to 350-700 bar, depending on the application. The compressed gas tanks in use today are bulky and cost intensive. To address this challenge, material-based storage is one of the long-term alternatives considered and constitutes the focus of this talk. Material-based storage is broadly defined as hydrogen bound to solid materials, with its binding strength varying from physisorption to porous materials, such as zeolites and metal organic frameworks, to chemisorption in (complex) metal hydrides. The ultimate targets set by the U.S. Department of Energy for this technology include a system gravimetric capacity of 6.5 wt% and volumetric capacity of 40 g/L at 100 bar, operating temperatures ranging between -40 C and +40 C and adsorption/desorption timescales of < 5 min. Storage in the form of physisorbed or chemisorbed hydrogen has guided the materials research, which metal- organic framework and (complex) metal hydrides being the most promising materials classes. A variety of these materials have met one or more of the targets, but it has remained elusive for a single material system to meet all these stringent requirements. Nano-encapsulation and low-concentration chemical additives have previously been employed separately to overcome such challenges. Functionalization via atomic layer deposition (ALD), however, offers unique characteristics that make it suitable for both, nano-encapsulation and "doping" with low-concentration additives. This deposition technique has sub-monolayer thickness control, is highly conformal in high-surface area materials and is self-limiting, i.e., once the gas-phase precursor reacts with the available surface sites, the surface reactions stop. In this presentation, we will showcase examples where ALD, more generally vapor-phase functionalization, on (complex) metal hydrides and organic frameworks has improved the material properties for H2 storage. In our first study, Al2O3 was deposited on magnesium borohydride, Mg(BH4)2, at room temperature using trimethylaluminum (TMA) and water. From our findings, encouraging initial results were obtained: the H2 desorption temperature was lowered by 60-120 degrees C, the desorbed gravimetric H2 capacity at temperatures < 250 degrees C was doubled, and the desorption kinetics increased by a factor of ~6 compared to uncoated Mg(BH4)2. However, hydrolysis reactions caused by residual surface -OH groups from the ALD water-pulse degraded the sample substantially. Through this study, the use of TMA was observed to be highly reactive with the Mg(BH4)2 surface species, leading to the strategy of exposing the sample only to TMA. The relative mole fraction of the vapor-phase additive can be precisely controlled with the number of TMA pulses and pulse duration. By tuning these parameters, we show in our second study that 10 pulses of TMA at ambient conditions were able to decrease the H2 desorption temperature by ~100 degrees C while retaining <95 % of its H2 capacity. We applied this approach with other additives such as BBr3, TiCl4 and tetrahydrofuran, and demonstrate that this unique approach opens the door to a new class of molecular additives and catalysts which cannot easily be introduced with conventional mechano-chemical or solvent-based techniques for (complex) metal-hydrides. In the case of metal organic frameworks, ALD is a promising technique to functionalize the pores with metal atoms or functional groups able to tune the gas selectivity and binding energy, for which ~15 kJ/mol has been established as the optimal value for H2 storage in sorbent materials. Vapor-phase techniques such as ALD have numerous benefits over other functionalization tools for H2 storage materials opening innumerable opportunities to the field. To exploit these opportunities, the current limitations on room temperature and water-less ALD processes needs to be overcome which will greatly expand the possibilities of encapsulation and incorporation of additives for organic frameworks and (complex) metal hydrides.

atomic layer deposition↗

Hybridization effect on the x-ray absorption spectra for actinide materials: Application to PuB 4

Studying the local moment and 5 f -electron occupations sheds insight into the electronic behavior in actinide materials. X-ray absorption spectroscopy (XAS) has been a powerful tool to reveal the valence electronic structure when assisted with theoretical calculations. However, the analysis currently taken in the community on the branching ratio of the XAS spectra generally does not account for the hybridization effects between local f orbitals and conduction states. In this paper, we discuss an approach which employs the density functional theory plus Gutzwiller rotationally invariant slave boson method to obtain a local Hamiltonian for the single-impurity Anderson model, and calculates the XAS spectra by the exact diagonalization (ED) method. A customized numerical routine was implemented for the ED XAS part of the calculation. By applying this technique to the recently discovered 5 f -electron topological Kondo insulator Pu B 4 , we determined the signature of 5 f -electronic correlation effects in the theoretical x-ray spectra. Furthermore, we found that the Pu 5 f - 6 d hybridization effect provides an extra channel to mix the j = 5 / 2 and 7 / 2 orbitals in the 5 f valence. As a consequence, the resulting electron occupation number and spin-orbit coupling strength deviate from the intermediate-coupling regime.

36 MATERIALS SCIENCE↗