Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “domination number”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Elevated ionizing photon production efficiency in faint high-equivalent-width Lyman-α emitters

ABSTRACT While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z ≈ 4–5 continuum-faint Lyman-α emitters from deep VLT/MUSE spectroscopy and directly measure their H α emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than −16 (≈ 0.01 L⋆), implying a median Lyman-α equivalent width of 259 Å. By combining the H α measurement with the UV magnitude, we determine the ionizing photon production efficiency, ξion, a first for such faint galaxies. The measurement of log10 (ξion [Hz erg−1]) = 26.28 ($^{+0.28}_{-0.40}$) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower luminosity, Lyman-α-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4 × 10−4 and 0.008, with light-weighted ages less than 3 Myr.

Maseda, Michael V.↗
DOE OSTI2020

Multi-stage preconditioners for thermal–compositional–reactive flow in porous media

We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.

42 ENGINEERING↗
DOE OSTI2025

On resistive interchange, double tearing, and resonant and non-resonant infernal modes in spherical tokamak negative central shear discharges

We examine the linear and nonlinear stability of a sequence of reversed shear NSTX equilibria with the same toroidal current density and pressure profile but with different toroidal field strengths. All equilibria have two q = 2 surfaces and have q > 1 everywhere. For this sequence, which all have β P ∼ 0.5, as the minimum value of q decreases below about 1.4, the localized resistive interchange criterion is strongly violated, and the unstable mode with toroidal mode number n=1 changes its dominant poloidal mode number m from a (m,n)=(2,1) double tearing mode to a non-resonant (1,1) infernal mode. This (1,1) mode nonlinearly flattens both the current density and the pressure near the magnetic axis. The higher-n unstable modes are all resistive infernal modes localized around the surface where q has a minimum and the magnetic shear vanishes. They generally saturate nonlinearly at low amplitude but can cause magnetic surface breakup near the axis.

Spitzer resistivity↗
DOE OSTI2024

Hydrodynamic Mechanism for Clumping along the Equatorial Rings of SN1987A and Other Stars

An explanation for the origin and number of clumps along the equatorial ring of Supernova 1987A has eluded decades of research. Our linear analysis and hydrodynamic simulations of the expanding ring prior to the supernova reveal that it is subject to the Crow instability between vortex cores. Here, the dominant wave number is remarkably consistent with the number of clumps, suggesting that the Crow instability stimulates clump formation. Although the present analysis focuses on linear fluid flow, future nonlinear analysis and the incorporation of additional stellar physics may further elucidate the remnant structure and the evolution of the progenitor and other stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2024

Turbulence spreading effects on the ELM size and SOL width

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width (λ q ) broadening in small ELM regimes. Here, this study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (E r ) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal E r well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width (λ q ) is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width (λ q ) broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux Γ ε at the separatrix can be enhanced by increasing either pedestal E r flow shear or local SOL pressure gradient. The λ q is broadened as the fluctuation energy intensity flux Γ ε at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Causality, intermittence, and crossphase evolution during confinement transitions in the TJ-II stellarator

In this work, we study spontaneous electron to ion root transitions in TJ-II using Langmuir probes. By scanning the probe position on a shot to shot basis, we reconstruct a spatiotemporal map of the evolution of important turbulent quantities in the plasma edge region. We pay particular attention to the evolution of the cross phase between transport-relevant variables, showing the spatiotemporal evolution of this quantity for the first time, revealing the outward propagation of the changes associated with the transition. Additionally, we also compute the intermittence parameter, which allows us to conclude that the turbulence, although its amplitude increases, condenses in a reduced number of dominant modes and becomes less bursty. The causal relationship between variables is studied using the transfer entropy, clarifying the interactions between the main variables and offering a rather complete picture of the complex evolution of the plasma across the confinement transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Plasma self-driven current in tokamaks with magnetic islands

Magnetic island perturbations may cause a reduction in plasma self-driven current that is needed for tokamak operation. A novel effect on tokamak self-driven current revealed by global gyrokinetic simulations is due to magnetic-island-induced 3D electric potential structures, which have the same dominant mode numbers as that of the magnetic island, whereas centered at both the inner and outer edge of the island. The non-resonant potential islands are shown to drive a current through an efficient nonlinear parallel acceleration of electrons. In large aspect ratio (large-A) tokamak devices, this new effect can result in a significant global reduction of the electron bootstrap current when the island size is sufficiently large, in addition to the local current loss across the island region due to the pressure profile flattening. It is shown that there exists a critical magnetic island width for large-A tokamaks beyond which the electron bootstrap current loss is global and increases rapidly with the island size. As such, this process may introduce a size limit for tolerable magnetic islands in large-A tokamak devices in the context of steady state operation. On the other hand, the current loss caused by magnetic islands in low-A tokamaks such as spherical tokamak (ST) NSTX/U is minor. The reduction of the axisymmetric current by magnetic islands scales with the square of island width. However, the loss of the current is mainly local to the island region, and the pace of current loss as the island size increases is substantially slower compared to large-A tokamaks. In particular, the bootstrap current reduction in STs is even smaller in the reactor-relevant high-β p regime where neoclassical tearing modes are more likely to develop.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Systematic analysis of magnetic equilibrium reconstruction with eddy currents on LTX- β

Magnetic equilibrium reconstruction in the Lithium Tokamak Experiment-Beta (LTX-β) is complicated by strong eddy currents and toroidal asymmetries arising from its segmented, close-fitting conducting shell. In the earlier experiment LTX, these three-dimensional (3D) conductor effects significantly distorted diagnostic signals and challenged conventional axisymmetric reconstruction methods. In this work, we demonstrate that accurate plasma equilibria can be recovered on LTX-β by incorporating a small number of dominant eddy current modes derived from realistic conductor models. Using the open-source $\tt{TokaMaker}$ Grad–Shafranov solver, we reconstruct equilibria across a systematically selected set of LTX-β discharges and validate them against the legacy $\tt{PSI-Tri}$ hybrid 2D-3D code. Our results show that fully 2D $\tt{TokaMaker}$ reconstructions achieve significantly improved agreement with flux loop and Mirnov probe measurements, reducing total chi-squared fitting errors, especially during startup. Finally, we develop new hybrid 2D-3D $\tt{TokaMaker}$ reconstructions by integrating the $\tt{ThinCurr}$ 3D eddy current model, which we optimize for significant further reduction of chi-squared in most scenarios. These findings underscore the importance of realistic wall-current modeling in short-pulse tokamaks, and establish a physics-based reconstruction framework that is extensible to devices with complex passive structures and 3D wall interactions.

Grad–Shafranov solver↗
DOE OSTI2022

Identification of Distribution of Elastic Constants of SiC/SiC Braided Tubes

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiCf/SiCm) composite tubes due to their superior mechanical properties. There are many parameters in the manufacturing of SiCf/SiCm composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. A sensitivity analysis will be carried out to determine the dominant parameters among the abovementioned ones. We intend to determine the variability in elastic constants of SiCf/SiCm composite tubes considering the variability in dominant parameters. A chosen number of combinations of dominant variables will be sampled through Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in circumferential and longitudinal direction, Poisson's ratio and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used a surrogate to FE analysis. once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The variability can be calculated based on MCS-generated data which will be presented at the conference.

42 ENGINEERING↗
DOE OSTI2024

A New Track Trigger for Characterization of the Antiproton-Induced Background in the Mu2e Experiment

The Mu2e experiment at Fermilab will enable the search for the neutrinoless muon to electron conversion in the field of an Al nucleus, a charged lepton flavor violating process. If observed, there would be a clear indication of physics beyond the Standard Model. Mu2e aims to reach a single event sensitivity of $3 /times 10^{-17}$, improving from the previous limit by 4 orders of magnitude. This improvement relies on the development of trigger selection systems, designed to discard data from background-induced events by placing kinematic, topological cuts on a particle’s reconstructed track. One of the largest sources of background Mu2e faces is proton-antiproton annihilation. These annihilations produce a 2 GeV shower of particles, among which there could be an electron that mimics the conversion electron signal, with an expected number of 0.010 ± 0.010. The large uncertainty on this number is dominated by the systematic uncertainty associated with the theoretical production model. To better characterize this background, we have developed an antiproton trigger selection by taking advantage of the track multiplicity and topology of these events. We discuss the steps taken in this development and the first performance study of this trigger, evaluating the signal efficiency and background rate. This trigger is essential to enable a data-driven analysis targeting the reduction of the systematic uncertainty of the antiproton-induced background in the Mu2e experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Hidden sector monopole dark matter with matter domination

The thermal freeze-out mechanism for relic dark matter heavier than O (10 – 100 TeV) requires cross-sections that violate perturbative unitarity. Yet the existence of dark matter heavier than these scales is certainly plausible from a particle physics perspective, pointing to the need for a non-thermal cosmological history for such theories. Topological dark matter is a well-motivated scenario of this kind. Here the hidden-sector dark matter can be produced in abundance through the Kibble-Zurek mechanism describing the non-equilibrium dynamics of defects produced in a second order phase transition. We revisit the original topological dark matter scenario, focusing on hidden-sector magnetic monopoles, and consider more general cosmological histories. We find that a monopole mass of order (1–10 5 ) PeV is generic for the thermal histories considered here, if monopoles are to entirely reproduce the current abundance of dark matter. In particular, in a scenario involving an early era of matter domination, the monopole number density is always less than or equal to that in a pure radiation dominated equivalent provided a certain condition on critical exponents is satisfied. This results in a larger monopole mass needed to account for a fixed relic abundance in such cosmologies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

High-Fidelity Analysis of ULF Wave Mode Structure Following Interplanetary Shock Compression of the Dayside Magnetopause Using MMS Multi-Point Observations

During the 9 March 2018 event with two consecutive interplanetary shocks compressing the dayside magnetosphere, the azimuthal mode structure and frequency spectrum of ultra low frequency magnetic pulsations are resolved using a cross-spectral analysis based on high-fidelity multi-probe Magnetospheric Multiscale Mission (MMS) magnetometer data. The results based on the MMS 4 and MMS 3 pair of measurements show that shock arrival leads to low mode (|m| ≤ 3) magnetic fluctuations in the Pc4-5 regimes, and smaller spatial scale fluctuations implied by the dominant high mode numbers are observed after both shock signatures hit and passed the magnetosphere. Detailed evolution of the mode structure is also shown for the first shock to reveal the development of high mode structure from a bump-on-tail distribution at m≈20 to a dominant peak at m≈50 in about 10 min. In addition, an interesting change of sign in m from negative to positive is observed as MMS crosses ~11 MLT pre-noon, which is consistent with the picture of wave generation by dayside magnetopause compression and then anti-sunward propagation. Finally, for both shocks, the contribution of higher frequency waves (Pc-4 range compared with Pc-5) to the total wave power is found to be negligible before and after the shock impact, but it becomes more significant during the shock impact.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Percolation and nanosecond fluctuators in V 2 O 3 films within the metal–insulator transition

Vanadium sesquioxide (V 2 O 3 ) exhibits a metal–insulator transition (MIT) at 160 K between a low temperature, monoclinic, antiferromagnetic Mott insulator and a high temperature, rhombohedral, paramagnetic, metallic phase. In thin films, a percolative transition takes place over a finite temperature range of phase coexistence. We study the fluctuating dynamics of this percolative MIT by measuring voltage noise spectra at both low frequencies (up to 100 kHz) and radio frequencies (between 10 MHz and 1 GHz). Noise intensity quadratic in bias is observed in the MIT region, as expected for resistive fluctuations probed nonperturbatively by the current. The low frequency noise resembles flicker-type 1/f β noise, often taking on the form of Lorentzian noise dominated by a small number of fluctuators as the volume fraction of the insulating phase dominates. Radio frequency noise intensity also quadratic in the bias current allows the identification of resistance fluctuations with lifetimes below 1 ns, approaching timescales seen in non-equilibrium pump–probe studies of the transition. We find quantitative consistency with a model for fluctuations in the percolative fraction. The thermodynamics of the MIT suggests that dominant fluctuations are ones that alter small volumes affecting the connectivity of domain boundaries. This noise serves as a sensitive and nonperturbative probe for the dynamics of switching phenomena in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cloud Morphology Evolution in Arctic Cold‐Air Outbreak: Two Cases During COMBLE Period

Abstract Cloud feedbacks play an important role in Arctic warming. Cloud morphology, for example, cloud size and spatial distributions, is among key factors that directly impact their radiative effects. In this work, we use two cases observed during the Cold‐air Outbreak (CAO) in the Marine Boundary Layer Experiment (COMBLE) to study the evolution of cloud size distributions as an air mass is advected from the Arctic over a comparatively warm ocean and cloud mesoscale organization changes from rolls to cells. Cloud objects are identified from Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance images through an object segmentation procedure and roll breakup is identified by homogeneities in cloud water path (CWP). Roll breakup is found to be accompanied by a local minimum in wind shear and local maxima in cloud size and marine cold air outbreak index. The mean cloud horizontal aspect ratio has weak fetch dependency and is around 2 in roll, transition, and cell regimes. Regardless of distance from the ice edge, smaller clouds (<10 km 2 ) dominate the population number but not cloud cover. Cloud size distributions show bimodality in transition and cell regimes. For clouds with comparable sizes, mean nearest neighbor distances normalized by equivalent cloud radius converge to a single value for all regimes and for all but the smallest clouds, suggesting that clouds of comparable sizes in CAOs are separated by distance proportional to their sizes. The presented statistical results pave the way to evaluating model simulated cloud organizations during CAO events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗