Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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

Structure-Activity Relationships of Hydrothermally Aged Titania-Supported Vanadium-Tungsten Oxide Catalysts for SCR of NOx Emissions with NH3

Supported V2O5-WO3/TiO2 materials employed as selective catalytic reduction (SCR) catalysts for NOx emission control from power plants. Fresh SCR catalysts usually receive exposure to harsh treatments in industry to accelerate catalyst activation (calcination in air at 650oC) and catalyst aging (hydrothermal aging at 650oC) in a way that represents various points in the catalyst/product lifetime. The present study investigates the catalyst structural and chemical changes occurring during such harsh treatments. Three series of supported V2O5-WO3/TiO2 catalysts were prepared by incipient-wetness impregnation of aqueous ammonium metavanadate and metatungstate precursors. The catalysts were subsequently dried and calcined at 550oC in O2, 650oC in O2, and hydrothermal conditions (10% O2, 8% H2O, 7% CO2 and 75% N2) at 650oC. The resulting catalysts were physically characterized by numerous techniques (in situ Raman, in situ IR, in situ High Field-High Spinning solid-state 51V MAS NMR, in situ EPR, XRD, BET surface area and ICP) and chemically probed with adsorbed ammonia, SCR-TPSR, and the SCR reaction. The surface WOx sites on the TiO2 support behave as a textural promoter that stabilizes the TiO2 (anatase) phase from sintering and transforming to the undesirable crystalline TiO2 (rutile) phase that can lead to formation of a Ti1-xVxO2 (rutile) solid solution with reduced V4+ cations (~7-15%). The surface VOx sites are mostly oligomerized as surface V+5Ox sites (~50-85% oligomers) and the extent of oligomerization tends to increase with surface WOx coverage and calcination temperature. A major difference between the calcined and hydrothermally treated catalysts was the low concentration of surface NH3* species on Lewis acid sites for the hydrothermally treated catalysts, yet the SCR activity was almost comparable for both catalysts. This finding suggests that surface NH4+*, primarily associated with the surface VOx sites, are able to efficiently perform the SCR reaction. Given that multiple catalyst parameters were simultaneously varying during these treatments, it was difficult to correlate the SCR activity with any single catalyst parameter. A correlation, however, was found between the SCR TOF/activity and the sum of the surface NH3* and NH4+* species, which is dominated by the surface NH4+* species.

SCR, NO, NH3, catalyst, V2O5, WO3, TiO2, promotion↗
DOE OSTI2023

Offshore-onshore record of Last Glacial Maximum–to–present grounding line retreat at Pine Island Glacier, Antarctica

Pine Island Glacier, West Antarctica, is the largest Antarctic contributor to global sea-level rise and is vulnerable to rapid retreat, yet our knowledge of its deglacial history since the Last Glacial Maximum is based largely on marine sediments that record a retreat history ending in the early Holocene. Using a suite of 10 Be exposure ages from onshore glacial deposits directly adjacent to Pine Island Glacier, we show that this major glacier thinned rapidly in the early to mid-Holocene. Our results indicate that Pine Island Glacier was at least 690 m thicker than present prior to ca. 8 ka. We infer that the rapid thinning detected at the site farthest downstream records the arrival and stabilization of the retreating grounding line at that site by 8–6 ka. By combining our exposure ages and the marine record, we extend knowledge of Pine Island Glacier retreat both spatially and temporally: to 50 km from the modern grounding line and to the mid-Holocene, providing a data set that is important for future numerical ice-sheet model validation.

58 GEOSCIENCES↗
DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Alternating and Gaussian Fermionic Isometric Tensor Network States

Isometric tensor networks in two dimensions enable efficient and accurate study of quantum many-body states, yet the effect of the isometric restriction on the represented quantum states is not fully understood. We address this question in two main contributions. First, we introduce an improved variant of isometric tensor network states (isoTNS) in two dimensions, where the isometric arrows on the columns of the network alternate between pointing upward and downward; hence the name alternating isometric tensor network states. Second, we introduce a numerical tool—the isometric Gaussian fermionic TNS (isoGfTNS)—that incorporates isometric constraints into the framework of Gaussian fermionic tensor network states. We demonstrate in numerous ways that alternating isoTNSs represent many-body ground states of two-dimensional quantum systems significantly better than the original isoTNSs. First, we show that the entanglement in an isoTNS is mediated along the isometric arrows and that alternating isoTNSs mediate entanglement more efficiently than conventional isoTNSs. Second, alternating isoTNSs correspond to a deeper, and thus more representative, sequential-circuit construction of depth 𝒪⁢(𝐿𝑥 ⋅𝐿𝑦) compared to the original isoTNSs of depth 𝒪⁢(𝐿𝑥 +𝐿𝑦). Third, using the Gaussian framework and gradient-based energy minimization, we provide numerical evidence of better bond-dimension scaling and variational energy of alternating isoGfTNSs for ground states of various free-fermionic models, including the Fermi surface, the band insulator, and the 𝑝𝑥 +𝑖⁢𝑝𝑦 mean-field superconductor. Finally, benchmarking on the transverse-field Ising model, we demonstrate that an alternating isoTNS provides substantially improved performance and stability relative to the original isoTNS for the ground-state search algorithm in interacting systems.

Wu, Yantao [Chinese Academy of Sciences, Beijing (↗
DOE OSTI2024

Energetic Particles - microturbulence interaction thrust (EP SIWG white paper)

This white paper outlines a strategic approach to Energetic Particle (EP) research, aimed at applying the results to the future burning devices such as ITER. Our goal is to deepen the understanding of the interplay between plasma background microturbulence and diverse species of energetic ions in burning plasma devices. This strategy is poised to advance our knowledge in this crucial area of fusion science, paving the way for significant breakthroughs in plasma physics and fusion energy research. We consider the confinement of energetic ions such as auxiliary heating beams, minority ions from ICRH, and fusion product alpha particles. Recent studies have underscored the potential of EP/microturbulence interaction to enhance fusion plasma performance by stabilizing microturbulence. However, it also poses risks, such as exacerbating Alfvénic eigenmode (AE) instabilities, which could compromise the sustainability of plasma discharge by causing EP losses. The interplay between EPs and background microturbulence is evident in the phenomenon of effective pitch angle scattering, a crucial aspect of the quasilinear (QL) theory. This theory is integral to developing numerically-efficient yet comprehensive and self-consistent approaches, recently employed to investigate the relaxation of energetic particle populations in the holistic modeling of fusion-grade plasmas. The primary objectives of this white paper thrust encompass a dual focus. First is that we need to delve into the mechanisms through which microturbulence engenders effective pitch angle scattering. This investigation will entail the examination of the presence of known or self-consistently predicted spectra of modes accountable for microturbulence. In these explorations, electrostatic microturbulence serves as a logical initial stride towards achieving the trust's objectives. This endeavor is anticipated to yield formulations expressing the parametric dependencies of the effective pitch angle scattering frequency on variables such as thermal plasma electron and/or ion temperatures, as well as their respective thermal conductivities. In our second objective, we aim to delve into the intricate formation of zonal flow (ZF) structures amidst the complex interplay of microturbulence and Alfvenic eigenmodes (AEs). This endeavor poses greater challenges as we seek to unravel the macroscopic manifestations influenced by microturbulence, termed as zonal structures (ZS), stemming from microturbulence-induced ZF. Drawing from this understanding, we anticipate employing a QL approach to yield comprehensive insights into the distribution function of energetic particles (EP) within phase space. This method entails resolving the dominant multidimensional phase space diffusion processes while effectively averaging over the rapid ballistic responses. Nevertheless, substantial strides remain imperative to realize a comprehensive whole-device modeling framework. This entails meticulous verification and validation exercises against experimental observations, as well as rigorous benchmarking against theoretical frameworks and numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2023

Photochemical Decomposition of Y‐Series Non‐Fullerene Acceptors Is Responsible for Degradation of High‐Efficiency Organic Solar Cells

Abstract Organic photovoltaic cells that employ Y‐series non‐fullerene acceptors (NFAs) have recently achieved impressive power‐conversion efficiencies (>18%). To fulfill their commercial promise, it is important to quantify their operational lifetimes and understand their degradation mechanisms. In this work, the spectral‐dependent photostability of films and solar cells comprising several Y‐series acceptors and the donor polymer PM6 is investigated systematically. By applying longpass filters during aging, it is shown that UV/near‐UV photons are responsible for the photochemical decomposition of Y‐series acceptors; this degradation is the primary driver of early solar cell performance losses. Using mass spectrometry, the vinylene linkage between the core and electron‐accepting moieties of Y‐series acceptors is identified as the weak point susceptible to cleavage under UV‐illumination. Employing a series of device characterization, along with numerical simulations, the efficiency losses in organic photovoltaic cells are attributed to the formation of traps, which reduces charge extraction efficiency and facilitates non‐radiative recombination as the Y‐series acceptors degrade. This study provides new insights for molecular degradation of organic photovoltaic absorber materials and highlights the importance of future molecular design and strategies for improved solar cell stability.

14 SOLAR ENERGY↗
DOE OSTI2021

Evaluation of pressure reconstruction techniques for Model Order Reduction in incompressible convective heat transfer

This paper compares pressure reconstruction strategies in Model Order Reduction for incompressible flows with convective heat transfer. The Navier-Stokes equation are reduced along with the passive scalar transport equation for the temperature using the POD-Galerkin technique. Six different pressure reconstruction methods are evaluated, two of which are novel to the best of the authors’ knowledge. Accurate pressure reconstruction is key to avoid error buildup in when solving for the conservation of linear momentum at the reduced level. The six approaches are compared using Direct Numerical Simulations of convective heat exchange processes in a 3D Backward Facing Step with a heated cylinder. Additionally, when comparing time-averaged metrics, we observe that the reconstruction methods that approximate the reduced pressure field using techniques borrowed from full order models (mechanical analogy, pressure Poisson, and velocity supremizers) yield higher errors than the methods that seek to stabilize the reduced systems (reduced residual stabilization, artificial divergence, and Uzawa operator).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Absolute calibration of the Lyman-α measurement apparatus at DIII-D

The LLAMA (Lyman-Alpha Measurement Apparatus) diagnostic was recently installed on the DIII-D tokamak [Rosenthal et al., Rev. Sci. Instrum. (submitted) (2020)]. LLAMA is a pinhole camera system with a narrow band Bragg mirror, a bandpass interference filter, and an absolute extreme ultraviolet photodiode detector array, which measures the Ly-α brightness in the toroidal direction on the inboard, high field side (HFS) and outboard, low field side (LFS). This contribution presents a setup and a procedure for an absolute calibration near the Ly-α line at 121.6 nm. The LLAMA in-vacuum components are designed as a compact, transferable setup that can be mounted in an ex situ vacuum enclosure that is equipped with an absolutely calibrated Ly-α source. The spectral purity and stability of the Ly-α source are characterized using a vacuum ultraviolet spectrometer, while the Ly-α source brightness is measured by a NIST-calibrated photodiode. The non-uniform nature of the Ly-α source emission was overcome by performing a calibration procedure that scans the Ly-α source position and employs a numerical optimization to determine the emission pattern. Nominal and measured calibration factors are determined and compared, showing agreement within their uncertainties. A first conversion of the measured signal obtained from DIII-D indicates that the Ly-α brightness on the HFS and LFS is on the order of 1020 Ph sr-1 m-2 s-1. The established calibration setup and procedure will be regularly used to re-calibrate the LLAMA during DIII-D vents to monitor possible degradation of optical components and detectors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Three-Magnon Bound State in the Quasi-One-Dimensional Antiferromagnet α-NaMnO 2

Here we report on the formation of a three-magnon bound state in the quasi-one-dimensional antiferromagnet α-NaMnO 2 , where the single-ion, uniaxial anisotropy inherent to the Mn 3+ ions in this material provides a binding mechanism capable of stabilizing higher order magnon bound states. While such states have long remained elusive in studies of antiferromagnetic chains, neutron scattering data presented here demonstrate that higher order n > 2 composite magnons exist, and, speci fically, that a weak three-magnon bound state is detected below the antiferromagnetic ordering transition of NaMnO 2 . We corroborate our findings with exact numerical simulations of a one-dimensional Heisenberg chain with easy-axis anisotropy using matrix-product state techniques, nding a good quantitative agreement with the experiment. Furthermore, these results establish NaMnO 2 as a unique platform for exploring the dynamics of composite magnon states inherent to a classical antiferromagnetic spin chain with Ising-like single ion anisotropy.

1-dimensional spin chains↗
DOE OSTI2025

Evaluation of data driven low-rank matrix factorization for accelerated solutions of the Vlasov equation

Low-rank methods have shown success in accelerating simulations of a collisionless plasma described by the Vlasov equation, but still rely on computationally costly linear algebra every time step. We propose a data-driven factorization method using artificial neural networks, specifically with convolutional layer architecture, that trains on existing simulation data. At inference time, the model outputs a low-rank decomposition of the distribution field of the charged particles, and we demonstrate that this step is faster than the standard linear algebra technique. Numerical experiments show that the method achieves comparable reconstruction accuracy for interpolation tasks, generalizing to unseen test data in a manner beyond just memorizing training data; patterns in factorization also inherently followed the same numerical trend as those within algebraic methods (e.g., truncated singular-value decomposition). However, when training on the first 70% of a time-series data and testing on the remaining 30%, the method fails to meaningfully extrapolate. Despite this limiting result, the technique may have benefits for simulations in a statistical steady-state or otherwise showing temporal stability. These results suggest that while the model offers a computationally efficient alternative for datasets with temporal stability, its current formulation is best suited for interpolation rather than for predicting future states in time-evolving systems. This study thus lays the groundwork for further refinement of neural network-based approaches to low-rank matrix factorization in high-dimensional plasma simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modeling Vapor Transport Deposition of Metal-halide Perovskite Thin Films for Photovoltaic and Optoelectronic Devices

Over the past decade, metal halide perovskites (MHPs) have emerged as a promising materials platform for high-efficiency solar cells and low-cost optoelectronics. However, there are challenges that frustrate the large-scale manufacturing of MHP devices, including difficulty in controlling film composition, interface formation and their device instability under ambient conditions. Vapor processing offers an attractive path to manufacturability, while also opening the door to new opportunities in device design that could favorably impact ultrahigh efficiency tandem solar cells or overall stability. Emphasis here is on the use of an alternate processing methodology, vapor transport deposition (VTD), to deposit the MHP layer. In experimental work, we found that our deposition system had significant run-to-run variations in film thickness and composition. To understand and resolve these issues, we use COMSOL Multiphysics software to model the precursor deposition rate numerically, trying to guide the broad range of parameters in the system. The result can be matched with experimental data and provides insight into system hydrodynamics and molar transport effects. Finally, we show that chamber pressure is a key factor to scale up the VTD technique.

Hsu, Wan-Ju↗
DOE OSTI2021

Dependence of alpha-particle-driven Alfven eigenmode linear stability on device magnetic field strength and consequences for next-generation tokamaks

Recently-proposed tokamak concepts use magnetic fields up to 12 T, far higher than in conven- tional devices, to reduce size and cost. Theoretical and computational study of trends in plasma behavior with increasing field strength is critical to such proposed devices. This paper considers trends in Alfven eigenmode (AE) stability. Energetic particles, including alphas from D-T fusion, can destabilize AEs, possibly causing loss of alpha heat and damage to the device. AEs are sensitive to device magnetic field via the field dependence of resonances, alpha particle beta, and alpha orbit width. We describe the origin and effect of these dependences analytically and by using recently- developed numerical techniques (Rodrigues et al. 2015 Nucl. Fusion 55 083003). The work suggests high-field machines where fusion-born alphas are sub-Alfvenic or nearly sub-Alfvenic may partially cut off AE resonances, reducing growth rates of AEs and the energy of alphas interacting with them. High-field burning plasma regimes have non-negligible alpha particle beta and AE drive, but faster slowing down time, provided by high electron density, and higher field strength reduces this drive relative to low-field machines with similar power densities. The toroidal mode number of the most unstable modes will tend to be higher in high magnetic field devices. The work suggests that high magnetic field devices have unique, and potentially advantageous, AE instability properties at both low and high densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling Vapor Transport Deposition of Metal-halide Perovskite Thin Films for Photovoltaic and Optoelectronic Devices

Over the past decade, metal halide perovskites (MHPs) have emerged as a promising materials platform for high-efficiency solar cells and low-cost optoelectronics. However, there are challenges that frustrate the large-scale manufacturing of MHP devices, including difficulty in controlling film composition, interface formation and their device instability under ambient conditions. Vapor processing offers an attractive path to manufacturability, while also opening the door to new opportunities in device design that could favorably impact ultrahigh efficiency tandem solar cells or overall stability. Emphasis here is on the use of an alternate processing methodology, vapor transport deposition (VTD), to deposit the MHP layer. In experimental work, we found that our deposition system had significant run-to-run variations in film thickness and composition. To understand and resolve these issues, we use COMSOL Multiphysics software to model the precursor deposition rate numerically, trying to guide the broad range of parameters in the system. The result can be matched with experimental data and provides insight into system hydrodynamics and molar transport effects. Finally, we show that chamber pressure is a key factor to scale up the VTD technique.

Hsu, Wan-Ju↗
DOE OSTI2023

A multifidelity deep operator network approach to closure for multiscale systems

Projection-based reduced order models (PROMs) have shown promise in representing the behavior of multiscale systems using a small set of generalized (or latent) variables. Despite their success, PROMs can be susceptible to inaccuracies, even instabilities, due to the improper accounting of the interaction between the resolved and unresolved scales of the multiscale system (known as the closure problem). In the current work, we interpret closure as a multifidelity problem and use a multifidelity deep operator network (DeepONet) framework to address it. In addition, to enhance the stability and/or accuracy of the multifidelity-based closure, we employ the recently developed ``in-the-loop" training approach from the literature on coupling physics and machine learning models. The resulting approach is tested on shock advection for the one-dimensional viscous Burgers equation and vortex merging for the two-dimensional Navier-Stokes equations. Furthermore, the numerical experiments show significant improvement of the predictive ability of the closure-corrected PROM over the un-corrected one both in the interpolative and the extrapolative regimes.

42 ENGINEERING↗
DOE OSTI2021

Behavior of solitary waves of coupled nonlinear Schrödinger equations subjected to complex external periodic potentials with odd-$\mathcal{PT}$ symmetry

In this work, we discuss the response of both moving and trapped solitary wave solutions of a two-component nonlinear Schrödinger system in 1 + 1 dimensions to an odd-$\mathcal{PT}$ external periodic complex potential. The dynamical behavior of perturbed solitary waves is explored by conducting numerical simulations of the nonlinear system and using a collective coordinate variational approximation. We present case examples corresponding to choices of parameter values and initial conditions involved therein. The results of the collective coordinate approximation are compared against numerical simulations where we observe qualitatively good agreement between the two. Unlike the case for a single-component solitary wave in a complex periodic $\mathcal{PT}$-symmetric potential, the collective coordinate equations do not have a small oscillation regime, and initially the height of the two components changes in opposite directions often causing instability. We find that the dynamic stability criteria we have used in the one-component case are a good indicator for the onset of dynamic instabilities in the present setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗