Search NASA⌕ Search

SEARCH · Search NASA

Results for “fast time simulation”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Penalized ensemble Kalman filters for high dimensional non-linear systems

The ensemble Kalman filter (EnKF) is a data assimilation technique that uses an ensemble of models, updated with data, to track the time evolution of a usually non-linear system. It does so by using an empirical approximation to the well-known Kalman filter. However, its performance can suffer when the ensemble size is smaller than the state space, as is often necessary for computationally burdensome models. This scenario means that the empirical estimate of the state covariance is not full rank and possibly quite noisy. To solve this problem in this high dimensional regime, we propose a computationally fast and easy to implement algorithm called the penalized ensemble Kalman filter (PEnKF). Under certain conditions, it can be theoretically proven that the PEnKF will be accurate (the estimation error will converge to zero) despite having fewer ensemble members than state dimensions. Further, as contrasted to localization methods, the proposed approach learns the covariance structure associated with the dynamical system. These theoretical results are supported with simulations of several non-linear and high dimensional systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Mimicking Natural Photosynthesis: Designing Ultrafast Photosensitized Electron Transfer into Multiheme Cytochrome Protein Nanowires

Efficient nanomaterials for artificial photosynthesis require fast and robust unidirectional electron transfer (ET) from photosensitizers through charge-separation and accumulation units to redox-active catalytic sites. We explored the ultrafast time-scale limits of photo-induced charge transfer between a Ru(II)tris(bipyridine) derivative photosensitizer and PpcA, a 3-heme c-type cytochrome serving as a nanoscale biological wire. Four covalent attachment sites (K28C, K29C, K52C, and G53C) were engineered in PpcA enabling site-specific covalent labeling with expected donor-acceptor (DA) distances of 4–8 Å. X-ray scattering results demonstrated that mutations and chemical labeling did not disrupt the structure of the proteins. Time-resolved spectroscopy revealed three orders of magnitude difference in charge transfer rates for the systems with otherwise similar DA distances and the same number of covalent bonds separating donors and acceptors. All-atom molecular dynamics simulations provided additional insight into the structure-function requirements for ultrafast charge transfer and the requirement of van der Waals contact between aromatic atoms of photosensitizers and hemes in order to observe sub-nanosecond ET. This work demonstrates opportunities to utilize multi-heme c-cytochromes as frameworks for designing ultrafast light-driven ET into charge-accumulating biohybrid model systems, and ultimately for mimicking the photosynthetic paradigm of efficiently coupling ultrafast, light-driven electron transfer chemistry to multi-step catalysis within small, experimentally versatile photosynthetic biohybrid assemblies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Double echo symmetry-based REDOR and RESPDOR pulse sequences for proton detected measurements of heteronuclear dipolar coupling constants

1 H{X} symmetry-based rotational echo double resonance pulse sequences (S-REDOR) and symmetry-based rotational echo saturation pulse double resonance (S-RESPDOR) solid-state NMR experiments have found widespread application for 1 H detected measurements of difference NMR spectra, dipolar coupling constants, and internuclear distances under conditions of fast magic angle spinning (MAS). In these experiments the supercycled $\text{R4}^2_1 (\text{SR4}^2_1)$ symmetry-based recoupling pulse sequence is typically applied to the 1 H spins to reintroduce heteronuclear dipolar couplings. However, the timing of $\text{SR4}^2_1$ and other symmetry-based pulse sequences must be precisely synchronized with the rotation of the sample, otherwise, the evolution of 1 H CSA and other interactions will not be properly refocused. For this reason, significant distortions are often observed in experimental dipolar dephasing difference curves obtained with S-REDOR or S-RESPDOR pulse sequences. Here we introduce a family of double echo (DE) S-REDOR/S-RESPDOR pulse sequences that function in an analogous manner to the recently introduced $t_1$-noise eliminated (TONE) family of dipolar heteronuclear multiple quantum coherence (D-HMQC) pulse sequences. Through numerical simulations and experiments the DE S-REDOR/S-RESPDOR sequences are shown to provide dephasing difference curves similar to those obtained with S-REDOR/S-RESPDOR. However, the DE sequences are more robust to the deviations of the MAS frequency from the ideal value that occurs during typical solid-state NMR experiments. In this work, the DE sequences are shown to provide more reliable 1 H detected dipolar dephasing difference curves for nuclei such as 15 N (with isotopic labelling), 183 W and 35 Cl. The double echo sequences are therefore recommended to be used in place of conventional S-REDOR/S-RESPDOR sequences for measurement of weak dipolar coupling constants and long-range distances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Insight into the Chemistry of PETN Under Shock Compression Through Ultrafast Broadband Mid-Infrared Absorption Spectroscopy

In this study, thin films of pentaerythritol tetranitrate (PETN) were shock compressed using the laser driven shock apparatus at Los Alamos National Laboratory (LANL). Two spectroscopic probes were available to this apparatus: visible white light transient absorption spectroscopy (VIS) from 400 to 700 nm and mid-infrared transient absorption spectroscopy (MIR) from 1150 to 3800 cm –1 . Important PETN vibrational modes are the symmetric and antisymmetric NO 2 stretches at 1280 and 1650 cm –1 , respectively, as well as CH stretches at ~2900 cm –1 . Shock strength was varied from approximately 3 to 55 GPa to span from the chemically unreactive regime to the regime in which fast chemical reaction took place on the 250 ps time scale of the measurements. VIS and MIR results suggest irreversible chemistry was induced in PETN at pressures above 30 GPa. At lower shock pressures, the spectroscopy showed minimal changes attributable to pressure induced effects. Under the higher-pressure reactive conditions, the frequency region at the antisymmetric NO 2 stretch mode had a significantly increased absorption while the region around the symmetric NO 2 stretch did not. No observable increased absorption occurred in the higher frequency regions where CH-, NH-, and OH- bond absorptions would be observed. A broad absorption appeared on the shoulder at the red-edge of the CO 2 vibrational band around 2200 cm –1 . In addition to the experiments, reactive molecular dynamics were carried out under equivalent shock conditions to correlate the evolution of the infrared spectrum to molecular processes. The simulations show results consistent to experiments up to 30 GPa but suggest that NO and NO 2 related features provided the strongest contributions to the shocked infrared changes. Proposed mechanisms for shocked PETN chemistry are analyzed as consistent or inconsistent with the data presented here. Our experimental data suggests C≡O or N 2 O bond formation, nitrite formation, and absence of significant hydroxyl or amine concentrations in the initial chemistry steps in PETN shocked above 30 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. Here, we report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Macro Analysis to Estimate Electric Vehicles Fast-Charging Infrastructure Requirements in Small Urban Areas

Electric vehicles (EVs) are known to reduce emissions and fossil fuel dependency. However, the limited range, long charging time, and inadequate charging infrastructure have hampered the adoption of EVs. The current EV charging infrastructure planning studies and tools require detailed information, extensive resources, and skills that can be a significant barrier to urban areas for finding the required charging infrastructure to support a targeted EV market share. This study generates regression models to estimate the number of direct current fast charging stations and the chargers to support the EV charging demand for urban areas. These models provide macro-level estimates of the required infrastructure investment in urban areas, which can be easily implemented by policy-makers and city planners. This study incorporates data obtained from applying a disaggregate optimization-based charger placement model, developed recently by the same authors, for multiple case studies to generate the required data to calibrate the macro-level models, in the state of Michigan. This simulated data set includes the number of charging stations and chargers for each market share, technology advancement scenario, and the transportation network topology. The results show that the number of charging stations reduces with battery size and charging power and increases with EV market share and the road network lane length. The number of chargers reduces with charging power, whereas it increases with battery size, EV market share, and vehicle miles traveled in the system. The model developed here can be applied to any state having urban characteristics and weather conditions similar to Michigan.

Engineering↗
DOE OSTI2020

Study of the Mu2e sensitivity to the $\mu^-\rightarrow e^+$ conversion process

The Mu2e experiment at Fermilab (Batavia, IL) will search for the coherent neutrinoless conversion of a muon into an electron in the electric field of an aluminum atom. The conversion process results in a monochromatic electron with energy slightly below the muon rest mass ($E_e= 104.97$ MeV). The experiment goal is to improve the previous measurement by four orders of magnitude and reach a single event sensitivity of $3 \times 10^{-17} $ on the conversion rate with respect to the muon capture rate. \\ The experiment will deliver a very intense pulsed negative muon beam to an aluminum target to achieve a total number of $10^{18}$ stopped muons in three years of running. Production and transport of the muons are accomplished with a sophisticated magnetic system comprised of a Production, a Transport and a Detector solenoid, with the last one hosting a Straw-tube Tracker and a Crystal Calorimeter. The entire detector region is surrounded by a Cosmic Ray Veto system. Mu2e is und er construction at the Muon Campus of Fermilab. \\ The INFN italian group is responsible for the design and construction of the crystal calorimeter, which plays a crucial role in the Mu2e measurement, providing particle identification capabilities fundamental for the rejection of the different background sources that can mimic a conversion electron. The calorimeter information allows also to improve the tracking performances it provides a fast standalone trigger for the experiment.\\ After a long R\&D phase, the final drawings see a calorimeter comprised of 1348 pure CsI crystals arranged in two annular disks. Each crystal is readout by two custom UV-extended Silicon Photomultipliers (SiPMs). Performances of the selected components were verified with an electron beam test, carried out in 2017 at the Beam Test Facility (BTF) at the National Laboratories of Frascati, on a large size prototype of 51 crystals and 102 SiPMs, confirming excellent energy ($\sigma_{E}/E = \mathcal{O }(10\%)$) and timing ($\sigma_{T} <500$ ps) resolutions in a! greement with the Mu2e specifications.\\ In the first chapters of this thesis, I describe my involvement in the calorimeter activities, with particular emphasis on the data analysis of the prototype energy response and resolution, that I was personally leading. Moreover, a description of the tuning of the simulation and the careful data-MC comparison is also reported.\\ An auxiliary measurement that can be performed with Mu2e is the search for a Lepton Number Violation (LNV) process where the negatively charged muon converts, without neutrinos, into a positron, $\mu ^{-} +N(A,Z) \rightarrow $ e$^+ +N(A,Z-2)$. At the moment of writing, Neutrinoless Double Beta Decay (0$\nu \beta \beta$) has set the most stringent limit on the LNV processes (T$^{0\nu \beta \beta}_{1/2} > 1.1 \times 10^{26} $ yr at 90$\%$ CL) with the GERDA experiment. In 0$\nu \beta \beta$ a pair of free electrons is created in the transformation from a nucleus $(A, Z)$ into its daughter $(A, Z + 2)$, namely: $(A,Z)\rightarrow (A, Z+2) =2 e^-$. Similar transitions could proceed through the emission of a pair of positrons, double electron capture (EC), or EC and single emission of a positron, with the nucleus changing from $(A, Z)$ to $(A, Z - 2)$. All these variations are equally interesting when discussing of new Physics, since they all manifest a non-conservation of the number of leptons.\\ Moreover, recent studies show that, because of flavor effect, off-diagonal processes, such as the $\mu^- -e^+$ conversion, can have their rate enhanced with respect to processes occurring in the $e-e$ sector.\\ Mu2e will search for $\mu^- + Al(27,13) \rightarrow e^+ +Na(27,11)$, only considering ground state transition, producing a conversion positron with energy $E_{e^+}= 92.32 $ MeV.\\ This thesis is focused on this search. In particular, I developed an algorithm to account for the radiative corrections to be applied to the conversion positrons that resulted in a few percent reduction in s ignal acceptance.\\ One of the most prominent background to ! this sear! ch is the Radiative Muon Capture (RMC) $ \mu+ N(Z,A)\rightarrow \gamma + \nu_{\mu} + N(Z-1,A)$. The produced photon can convert, virtually or in the detector materials (in the following referred as internal or external conversion respectively) in an electron-positron pair, producing positrons with energies approaching the expected conversion energy. Experimental results from previous experiments indicate the RMC endpoint ($\kmax$) to be of $\sim$ 90 MeV. A precise modeling of this background is relevant due to the proximity of the photon spectrum endpoint to the positron conversion energy. The calorimeter group is obviously interested to reconstruct the RMC distribution having also the opportunity to improve the current knowledge of the process. \\ A toy MC has been developed to evaluate the number of $e^+$ from external RMC background in the signal window. Knowing this value and ev...

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

MOOSE Framework Meshing Enhancements to Support Reactor Analysis

MOOSE-based physics codes require an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The use of difficult-touse, external licensed software is often required to generate high quality meshes for reactor geometries. High-fidelity geometry modeling also requires elaborate tracking of groups of elements for material property assignment and output reporting which can be considerably complex for the user to identify and maintain. Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, several meshingrelated enhancements have been developed for the MOOSE framework to address user challenges in creating finite element meshes for advanced reactor geometries. MOOSE mesh generators have been developed to mesh hexagonal geometries (pins, ducted assemblies, and cores) commonly found in liquid-metal cooled fast reactor concepts. The mesh generator used for hexagonal pin cells is generic for regular polygons and therefore may also be used for Cartesian pin cells. Hexagonal pin cells can be stitched into ducted assemblies, and assemblies can be stitched together into a core. The user may specify region ids, region names, and other preferences on the mesh. This control is useful for later material mapping in the MOOSE-based physics codes input. A capability was also developed for meshing rotating control drums including determination of material volume fractions in each mesh element as a function of time. Control drum meshes may be stitched to other hexagonal assemblies to create a core configuration. Additional mesh generators were developed that wrap around the hexagonal meshing capabilities and utilize “extra element integer” ID values on each element. In regular Cartesian or hexagonal assemblies or cores, the bookkeeping of element groups for both material assignment and output reporting can now be automated through assignment of pin, assembly, core, axial and depletion id values stored as extra element integers. The extra element tags on the mesh greatly speed the reactor analyst’s efforts to map materials to meshes, track depletion zones, and parse output such as axial pin power distributions. At the highest level, pin, assembly, and core mesh generators (with this reactor terminology) have also been developed to easily generate regular Cartesian and hexagonal cores, including axial extrusion. These reactor geometry builders call upon the previously mentioned capabilities to produce analysis-ready 3D meshes including material assignments. Open source mesh triangulation capabilities were also investigated for integration into the MOOSE framework to address the need for meshing the core periphery region which extends from the irregular outer assembly border to a cylindrical boundary. Options are limited due to licensing constraints, and the recommendation is pursue building a native MOOSE Delaunay triangulator routine with full functionality. Finally, a series of verification problems were performed with NEAMS physics tools. All developed capabilities will be available in the new open-source “Reactor” module of the MOOSE framework, which is accessible to any MOOSE-based NEAMS physics tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A comprehensive and fair comparison of two neural operators (with practical extensions) based on $\mathrm{FAIR}$ data

Neural operators can learn nonlinear mappings between function spaces and offer a new simulation paradigm for real-time prediction of complex dynamics for realistic diverse applications as well as for system identification in science and engineering. Herein, we investigate the performance of two neural operators, which have shown promising results so far, and we develop new practical extensions that will make them more accurate and robust and importantly more suitable for industrial-complexity applications. The first neural operator, DeepONet, was published in 2019 (Lu et al., 2019), and its original architecture was based on the universal approximation theorem of Chen & Chen (1995). The second one, named Fourier Neural Operator or FNO, was published in 2020, and it is based on parameterizing the integral kernel in the Fourier space. DeepONet is represented by a summation of products of neural networks (NNs), corresponding to the branch NN for the input function and the trunk NN for the output function; both NNs are general architectures, e.g., the branch NN can be replaced with a CNN or a ResNet. According to Kovachki et al. (2021), FNO in its continuous form can be viewed conceptually as a DeepONet with a specific architecture of the branch NN and a trunk NN represented by a trigonometric basis. In order to compare FNO with DeepONet computationally for realistic setups, we develop several extensions of FNO that can deal with complex geometric domains as well as mappings where the input and output function spaces are of different dimensions. We also develop an extended DeepONet with special features that provide inductive bias and accelerate training, and we present a faster implementation of DeepONet with cost comparable to the computational cost of FNO, which is based on the Fast Fourier Transform. Here we consider 16 different benchmarks to demonstrate the relative performance of the two neural operators, including instability wave analysis in hypersonic boundary layers, prediction of the vorticity field of a flapping airfoil, porous media simulations in complex-geometry domains, etc. We follow the guiding principles of FAIR (Findability, Accessibility, Interoperability, and Reusability) for scientific data management and stewardship. The performance of DeepONet and FNO is comparable for relatively simple settings, but for complex geometries the performance of FNO deteriorates greatly. We also compare theoretically the two neural operators and obtain similar error estimates for DeepONet and FNO under the same regularity assumptions.

42 ENGINEERING↗
DOE OSTI2020

Novel Fault Location Method for Power Systems Based on Attention Mechanism and Double Structure GRU Neural Network

Fault location is one of the most essential techniques to maintain the stable operation of power systems. A fast and accurate fault location allows operators to restore power grids faster and avoid economic losses. Conventional methods rely on expert knowledge to extract the necessary features (e.g. DWT, DFT). For large systems, more coupling effects of transmission lines require more complex feature engineering, and incomplete features can easily introduce large errors. To overcome this, a deep learning approach without manual feature extraction is introduced to the fault location model under big data application. Towards this end, in the proposed method, the attention mechanism, the Bi-GRU and a dual structure network are applied to analyze the current data from different perspectives. Complete information for the fault features is extracted for the fault location. Based on a time series model and benefit from the ability to internally acquire the information architecture of faulty line, the established model is adaptive to the power grids with very complex topologies. Simulation results indicate that the proposed double-structure model reduces the maximum error and is less affected by noise. In comparison with different structures and different models, the proposed method shows better performance in IEEE 39-bus system.

42 ENGINEERING↗
DOE OSTI2024

Forced Oscillation Source Location in Power Systems Using MVMD-assisted DEF in TF Plane

Fast and accurate forced oscillation source location (FOSL) is critical for mitigating forced oscillations in power systems. Here, to overcome the shortcomings of traditional dissipating energy flow (DEF) in decomposing the forced oscillation components from the measurement responses, a multivariate variational mode decomposition (MVMD)-assisted FOSL is proposed in this paper to locate the forced oscillation sources from measurements in time-frequency (TF) plane. The multi-channel measurement matrix of each generator is formed, then the multi-channel modes are decomposed from the formed multi-channel measurement matrix using the MVMD. Further, the intrinsic mode functions (IMFs) associated with the forced oscillation mode are distinguished from the decomposed modes according to the relative energy weights. With the distinguished IMFs, the MVMD-assisted DEF calculation method is developed to locate the forced oscillation source. Finally, the performance of the proposed method is evaluated using the simulation data of the WECC 179-bus test system and IEEE-NASPI Oscillation Source Location Contest as well as the actual PMU data of ISO New England. The results validated the effectiveness and accuracy of the proposed method in the FOSL.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation

In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Lyapunov stability of smart inverters using linearized distflow approximation

Fast-acting smart inverters that utilize preset operating conditions to determine real and reactive power injection/consumption can create voltage instabilities (over-voltage, voltage oscillations and more) in an electrical distribution network if set-points are not properly configured. In this work, linear distribution power flow equations and droop-based Volt–Var and Volt–Watt control curves are used to analytically derive a stability criterion using Lyapunov analysis that includes the network operating condition. The methodology is generally applicable for control curves that can be represented as Lipschitz functions. The derived Lipschitz constants account for smart inverter hardware limitations for reactive power generation. A local policy is derived from the stability criterion that allows inverters to adapt their control curves by monitoring only local voltage, thus avoiding centralized control or information sharing with other inverters. The criterion is independent of the internal time-delays of smart inverters. Simulation results for inverters with and without the proposed stabilization technique demonstrate how smart inverters can mitigate voltage oscillations locally and mitigate real and reactive power flow disturbances at the substation under multiple scenarios. The study concludes with illustrations of how the control policy can dampen oscillations caused by solar intermittency and cyberattacks.

42 ENGINEERING↗
DOE OSTI2025

Low-Cost, Screen-Printed Silver Metal Complex Inks for Silicon Heterojunction Solar Cells

Screen printing using metal particle pastes, the current photovoltaic industry metallization standard, provides fast and reliable metal grids for silicon solar cells. Recently, metal complex or reactive metal inks are attracting research interest due to their significantly low cost and higher performance compared to traditional nanoparticle silver pastes. In this work, we demonstrate, for the first time, screen-printed high-efficiency silicon heterojunction solar cells metallized by silver metal complex inks on industrial G1-size (158.75 x 158.75 mm2) wafers. We demonstrate screen-printed Ag metal complex ink grid patterns with continuous fingers ~100-120 ..mu..m wide. The printed Ag grid is very thin (~1 ..mu..m), which is an order of magnitude thinner than the current ~20-30 ..mu..m fingers printed with low-temperature nanoparticle-based pastes. Double printing allows silicon heterojunction devices with efficiencies >20%. This is the highest efficiency so far, to our knowledge, of industrial solar cell precursors using this metallization technology. Simulation results suggested that increasing the thickness of the metal film does not significantly improve efficiency due to the dense, highly conductive films. So, a single print of ~1 ..mu..m finger would be enough to produce cells that perform similarly to a ~20 ..mu..m thick nanoparticle paste printed cells. Additionally, solar cells printed on G1 wafers with silver metal complex ink required more than 10 times less silver (~0.03 g) compared to those using silver/copper nanoparticle paste (~0.4 g of Ag). These results indicate that metal complex inks are a very promising replacement for silver nanoparticle pastes for industrial-scale metallization in an age of resource scarcity and high costs of noble metals.

14 SOLAR ENERGY↗
DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗
DOE OSTI2026

Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations

Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Recovered supernova Ia rate from simulated LSST images

Aims.TheVera C. RubinObservatory’s Legacy Survey of Space and Time (LSST) will revolutionize time-domain astronomy by detecting millions of different transients. In particular, it is expected to increase the number of known type Ia supernovae (SN Ia) by a factor of 100 compared to existing samples up to redshift ∼1.2. Such a high number of events will dramatically reduce statistical uncertainties in the analysis of the properties and rates of these objects. However, the impact of all other sources of uncertainty on the measurement of the SN Ia rate must still be evaluated. The comprehension and reduction of such uncertainties will be fundamental both for cosmology and stellar evolution studies, as measuring the SN Ia rate can put constraints on the evolutionary scenarios of different SN Ia progenitors. Methods.We used simulated data from the Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) and LSST Data Preview 0 to measure the SN Ia rate on a 15 deg 2 region of the “wide-fast-deep” area. We selected a sample of SN candidates detected in difference images, associated them to the host galaxy with a specially developed algorithm, and retrieved their photometric redshifts. We then tested different light-curve classification methods, with and without redshift priors (albeit ignoring contamination from other transients, as DC2 contains only SN Ia). We discuss how the distribution in redshift measured for the SN candidates changes according to the selected host galaxy and redshift estimate. Results.We measured the SN Ia rate, analyzing the impact of uncertainties due to photometric redshift, host-galaxy association and classification on the distribution in redshift of the starting sample. We find that we are missing 17% of the SN Ia, on average, with respect to the simulated sample. As 10% of the mismatch is due to the uncertainty on the photometric redshift alone (which also affects classification when used as a prior), we conclude that this parameter is the major source of uncertainty. We discuss possible reduction of the errors in the measurement of the SN Ia rate, including synergies with other surveys, which may help us to use the rate to discriminate different progenitor models.

Astronomy & Astrophysics↗