SEARCH · Engineering Papers

Results for “Normalization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗

DOE OSTI2020

MHD stability and disruptions in the SPARC tokamak

SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C 2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B 2 , C 2 , O 2 , and F 2 in the local spin density approximation (LSDA), the Perdew–Burke–Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C 2 . Here, this work shows that spin symmetry breaking in singlet C 2 , which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

From Femtoseconds to Gigaseconds: The SolDeg Platform for the Performance Degradation Analysis of Silicon Heterojunction Solar Cells

Heterojunction Si solar cells exhibit notable performance degradation. Here, we modeled this degradation by electronic defects getting generated by thermal activation across energy barriers over time. To analyze the physics of this degradation, we developed the SolDeg platform to simulate the dynamics of electronic defect generation. First, femtosecond molecular dynamics simulations were performed to create a-Si/c-Si stacks, using the machine learning-based Gaussian approximation potential. Second, we created shocked clusters by a cluster blaster method. Third, the shocked clusters were analyzed to identify which of them supported electronic defects. Fourth, the distribution of energy barriers that control the generation of these electronic defects was determined. Fifth, an accelerated Monte Carlo method was developed to simulate the thermally activated time-dependent defect generation across the barriers. Our main conclusions are as follows. (1) The degradation of a-Si/c-Si heterojunction solar cells via defect generation is controlled by a broad distribution of energy barriers. (2) We developed the SolDeg platform to track the microscopic dynamics of defect generation across this wide barrier distribution and determined the time-dependent defect density N(t) from femtoseconds to gigaseconds, over 24 orders of magnitude in time. (3) We have shown that a stretched exponential analytical form can successfully describe the defect generation N(t) over at least 10 orders of magnitude in time. (4) We found that in relative terms, Voc degrades at a rate of 0.2%/year over the first year, slowing with advancing time. (5) We developed the time correspondence curve to calibrate and validate the accelerated testing of solar cells. We found a compellingly simple scaling relationship between accelerated and normal times t normal ∝ t accel T(accel)/T(normal) . (6) We also carried out experimental studies of defect generation in a-Si:H/c-Si stacks. We found a relatively high degradation rate at early times that slowed considerably at longer time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Interfacial Reactivity of Silicon Electrodes: Impact of Electrolyte Solvent and Presence of Conductive Carbon

Silicon (Si) is a promising high-capacity material for lithium-ion batteries; however, its limited reversibility hinders commercial adoption. Approaches such as particle and crystallite size reduction, introduction of conductive carbon, and use of different electrolyte solvents have been explored to overcome these electrochemical limitations. Herein, operando isothermal microcalorimetry (IMC) is used to probe the influence of silicon particle size, electrode composition, and electrolyte additives fluoroethylene carbonate and vinylene carbonate on the heat flow during silicon lithiation. In this work, the IMC data are complemented by X-ray photoelectron and Raman spectroscopies to elucidate differences in solid electrolyte interphase (SEI) composition. Nanosized (~50 nm, n-Si) and micrometer-sized (~4 μm, μ-Si) silicon electrodes are formulated with and without amorphous carbon and electrochemically lithiated in ethylene carbonate (EC), fluoroethylene carbonate (FEC), or vinylene carbonate (VC) based electrolytes. Notably, n-Si electrodes generate 53–61% more normalized heat relative to their μ-Si counterparts, consistent with increased surface area and electrode/electrolyte reactivity. Introduction of amorphous carbon significantly alters the heat flow profile where multiple exothermic peaks and increased normalized heat dissipation are observed for all electrolyte types. Notably, the VC-containing electrolyte demonstrates the greatest normalized heat dissipation of the electrode compositions tested showing as much as a 50% increase compared to the EC or FEC counterparts. The results are relevant to the understanding of silicon negative electrode function in the presence of electrolyte additives and provide insight relative to silicon containing cell reactivity and safety.

25 ENERGY STORAGE↗

DOE OSTI2022

Lower‐Hybrid Wave Structures and Interactions With Electrons Observed in Magnetotail Reconnection Diffusion Regions

Abstract We investigate waves close to the lower‐hybrid frequency in 12 magnetotail reconnection electron diffusion region (EDR) events with guide field levels of near‐zero to 30%. In about half of the events, the wave vector has a small component along the current sheet normal, consistent with known lower‐hybrid drift wave properties, but the perpendicular magnetic field fluctuations can be comparable or greater than the parallel component, a feature unique to the waves inside and adjacent to EDRs. Another new wave property is that the wave vector has a significant component along the current sheet normal in some events and completely along the normal for one event. In 1/4 of the events, the term has a significant contribution to the wave electric field, possibly a feature of lower‐hybrid waves more likely to exist in the diffusion region than further away from the X‐line. Electron temperature variations are correlated with the wave potential, due to wave electric field acceleration and crossings at the corrugated separatrix region with different amounts of mixing between reconnection inflowing and outflowing populations. The latter also leads to the anti‐correlation between parallel and perpendicular temperature components. Using four‐spacecraft measurements, the magnetic field line twisting is demonstrated by the correlated fluctuations in and . The lower‐hybrid wave in the EDR of weak guide field reconnection may be generated near separatrices and penetrate to the mid‐plane or locally generated, and the latter possibility is beyond the prediction of previous reconnection simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Amplified Extreme Floods and Shifting Flood Mechanisms in the Delaware River Basin in Future Climates

Abstract Historical records in the Delaware River Basin reveal complex and spatially diverse flood generating mechanisms influenced by the region's mountains‐to‐plains gradients. This study focuses on predicting future flood hazards and understanding the underlying drivers of changes across the region. Using a process‐based hydrological model, we analyzed the hydrometeorological condition of each historical and future flood event. For each event, at the subbasin scale, we identified the dominant flood generating mechanism, including snowmelt, rain‐on‐snow, short‐duration rain, and long‐duration rain. The rain‐induced floods are further categorized based on the soil's Antecedent Moisture Condition (AMC) before the event, whether dry, normal, or wet. Our historical analysis suggests that rain‐on‐snow is the primary flood mechanism of the Upper Basin. Although most frequent, the magnitude of rain‐on‐snow floods is often less severe than short rain floods. In contrast, historical floods in the Lower Basin are primarily caused by short rain under normal AMC. Given the uncertainties in climate projections, we used an ensemble of future climate scenarios for flood projections. Despite variations in regional climate projections, coherent perspectives emerge: the region will shift toward a warmer, wetter climate, with a projected intensification of extreme floods. The Upper Basin is projected to experience a marked decrease in rain‐on‐snow floods, but a substantial increase in short rain floods with wet AMC. The largest increase in flood magnitude will be driven by short rains with wet AMC in the Upper Basin and by short rains with normal AMC in the Lower Basin.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Resolution of a Few Problems in the Application of Quasilinear Theory to Calculating Diffusion Coefficients in Heliophysics

A theory of quasilinear diffusion for obliquely propagating electromagnetic waves was developed in the 1960's and applied in the 1970's to model scattering of relativistic electrons by a prescribed distribution of waves. In the latter work, a transformation of variables, from wavevector space to the temporal frequency and tangent of the wave normal angle, was used so that simple Gaussian functions of frequency and tangent of the wave normal angle could be multiplied together to define the distribution of wave power, although arbitrary distributions of power in these two variables is also permitted. Finally, in 2005, previous work was consolidated and has been widely used in heliophysics studies that require computation of quasilinear diffusion coefficients. Here it is shown that this transformation is inppropriate when the precise wave vector distribution is known. The correct transformation is derived and used to produce diffusion coefficients that can differ by orders of magnitude from those computed using the inappropriate transformation. The differences are largest when the distribution of wave power extends to wave normal angles near the resonance cone. When the ratio of the plasma frequency to the gyrofrequency is large, only low energies (keV) are affected, but as the ratio decreases higher energies (MeV) also show differences. It is also shown that the derivation from the 1960's uses a notation that results in the diffusion coefficients depending on the distribution of wave power with respect to the wave azimuthal angle whereas there should be no such dependence.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

The Turbulent Pressure Spectrum Within the Roughness Sublayer of a Subarctic Forest Canopy

The turbulent static pressure spectrum E pp (k x ) as a function of longitudinal wavenumber k x in the roughness sublayer of forested canopies is of interest to a plethora of problems such as pressure transport in the turbulent kinetic energy budget, pressure pumping from snow or forest floor, and coupling between flow within and above canopies. Long term static pressure measurements above a sub-arctic forested canopy for near-neutral conditions during the winter and spring were collected and analyzed for three snow cover conditions: trees and ground covered with snow, trees are snow free but the ground is covered with snow, and snow free cover. In all three cases, it is shown that E pp (k x ) obeys the attached eddy hypothesis at low wavenumbers (1/δ < k x < 1/z) —with E pp (k x ) ∝ u$^4_*$k$^{-1}_x$ and Kolmogorov scaling in the inertial subrange at higher wavenumbers—with E pp (k x ) ∝ ϵ 4/3 k x -7/3 , where u * is the friction velocity at the canopy top, is the mean turbulent kinetic energy dissipation rate, z is the distance from the snow top, and is the boundary layer depth. The implications of these two scaling laws to the normalized root-mean squared pressure C p = o p /u$^2_*$ and its newly proposed logarithmic scaling with normalized wall-normal distance z/δ are discussed for snow covered and snow free vegetation conditions. The work here also shows that k$^{-1}_x$ the in the E pp (k x ) appears more extensive and robust than its longitudinal velocity counterpart.

Aslan, Toprak [Finnish Meteorological Inst., Helsi↗

DOE OSTI2020

Deep neural network improves the estimation of polygenic risk scores for breast cancer

Polygenic risk scores (PRS) estimate the genetic risk of an individual for a complex disease based on many genetic variants across the whole genome. Here, we compared a series of computational models for estimation of breast cancer PRS. A deep neural network (DNN) was found to outperform alternative machine learning techniques and established statistical algorithms, including BLUP, BayesA, and LDpred. In the test cohort with 50% prevalence, the Area Under the receiver operating characteristic Curve (AUC) were 67.4% for DNN, 64.2% for BLUP, 64.5% for BayesA, and 62.4% for LDpred. BLUP, BayesA, and LPpred all generated PRS that followed a normal distribution in the case population. However, the PRS generated by DNN in the case population followed a bimodal distribution composed of two normal distributions with distinctly different means. This suggests that DNN was able to separate the case population into a high-genetic-risk case subpopulation with an average PRS significantly higher than the control population and a normal-genetic-risk case subpopulation with an average PRS similar to the control population. This allowed DNN to achieve 18.8% recall at 90% precision in the test cohort with 50% prevalence, which can be extrapolated to 65.4% recall at 20% precision in a general population with 12% prevalence. Interpretation of the DNN model identified salient variants that were assigned insignificant p values by association studies, but were important for DNN prediction. These variants may be associated with the phenotype through nonlinear relationships.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

The breakdown of both strange metal and superconducting states at a pressure-induced quantum critical point in iron-pnictide superconductors

Here we report the first observation of the concurrent breakdown of the strange metal (SM) normal state and superconductivity at a pressure-induced quantum critical point in Ca 10 (Pt 4 As 8 )((Fe 0.97 Pt 0.03 ) 2 As 2 ) 5 superconductor. We find that, upon suppressing the superconducting state, the power exponent (α) changes from 1 to 2, and the slope of the temperature-linear resistivity per FeAs layer (A $\Box$ ) gradually diminishes. At a critical pressure, A $\Box$ and superconducting transition temperature (T c ) go to zero concurrently, where a quantum phase transition from a superconducting state with a SM normal state to a non-superconducting Fermi liquid state occurs. Scaling analysis reveals that the change of A $\Box$ with T c obeys the relation of T c ~ (A $\Box$ ) 0.5 , similar to what is seen in other chemically doped unconventional superconductors. These results suggest that there is a simple but powerful organizational principle of connecting the SM normal state with the high-T c superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Base excision repair and double strand break repair cooperate to modulate the formation of unrepaired double strand breaks in mouse brain

Abstract We lack the fundamental information needed to understand how DNA damage in the brain is generated and how it is controlled over a lifetime in the absence of replication check points. To address these questions, here, we integrate cell-type and region-specific features of DNA repair activity in the normal brain. The brain has the same repair proteins as other tissues, but normal, canonical repair activity is unequal and is characterized by high base excision repair (BER) and low double strand break repair (DSBR). The natural imbalance creates conditions where single strand breaks (SSBs) can convert to double strand breaks (DSBs) and reversibly switch between states in response to oxidation both in vivo and in vitro. Our data suggest that, in a normal background of repair, SSBs and DSBs are in an equilibrium which is pushed or pulled by metabolic state. Interconversion of SSB to DSBs provides a physiological check point, which would allow the formation of unrepaired DSBs for productive functions, but would also restrict them from exceeding tolerable limits.

Science & Technology - Other Topics↗

DOE OSTI2024

Transport phase diagram and anomalous metallicity in superconducting infinite-layer nickelates

Despite obvious similarities in their electronic and crystallographic structures, it remains unclear whether the interactions that shape the normal and superconducting (SC) state properties of high-T c cuprates and infinite-layer nickelates (ILNs) have the same origin. This question has been brought into sharper focus with recent studies on ILNs of improved crystallinity that reveal a SC dome of comparable extent and similar transport properties above T c as the hole-doped cuprates. The evolution of these properties in the magnetic-field-induced normal state, however, has yet to be determined. Here, we examine the magnetotransport properties of new-generation Nd 1−x Sr x NiO 2 films in the T → 0 limit across the phase diagram in fields up to 54 T. This extensive study reveals that the limiting low-T form of the normal-state resistivity in ILNs exhibits non-Fermi-liquid behaviour over an extended doping range inside the SC dome, rather than at a singular quantum critical point. While there are clear differences in the charge dynamics of ILNs and cuprates, most notably in the magnetoresistance, our findings reveal that both systems exhibit anomalous metallicity characteristic of a quantum critical phase.

Hsu, Yu-Te↗

DOE OSTI2025

Joint neutrino oscillation analysis from the T2K and NOvA experiments

The landmark discovery that neutrinos have mass and can change type (or flavour) as they propagate—a process called neutrino oscillation—has opened up a rich array of theoretical and experimental questions being actively pursued today. Neutrino oscillation remains the most powerful experimental tool for addressing many of these questions, including whether neutrinos violate charge-parity (CP) symmetry, which has possible connections to the unexplained preponderance of matter over antimatter in the Universe. Oscillation measurements also probe the mass-squared differences between the different neutrino mass states (Δ$m^2$), whether there are two light states and a heavier one (normal ordering) or vice versa (inverted ordering), and the structure of neutrino mass and flavour mixing. Here we carry out the first joint analysis of datasets from NOvA and T2K, the two currently operating long-baseline neutrino oscillation experiments (hundreds of kilometres of neutrino travel distance), taking advantage of our complementary experimental designs and setting new constraints on several neutrino sector parameters. This analysis provides new precision on the $Δm_{32}^2$ mass difference, finding $-2.43_{−0.03}^{+0.04}$ x $10^{-13}$ $\textrm{eV}^2$ in the normal ordering and $-2.43^{+0.03}_{-0.04}$ x $10^{-13}$ $\textrm{eV}^2$ in the inverted ordering, as well as a 3σ interval on δ CP of [−1.38π, 0.30π] in the normal ordering and [−0.92π, −0.04π] in the inverted ordering. The data show no strong preference for either mass ordering, but notably, if inverted ordering were assumed true within the three-flavour mixing model, then our results would provide evidence of CP symmetry violation in the lepton sector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

A new platform for ultra-high dose rate radiobiological research using the BELLA PW laser proton beamline

Abstract Radiotherapy is the current standard of care for more than 50% of all cancer patients. Improvements in radiotherapy (RT) technology have increased tumor targeting and normal tissue sparing. Radiations at ultra-high dose rates required for FLASH-RT effects have sparked interest in potentially providing additional differential therapeutic benefits. We present a new experimental platform that is the first one to deliver petawatt laser-driven proton pulses of 2 MeV energy at 0.2 Hz repetition rate by means of a compact, tunable active plasma lens beamline to biological samples. Cell monolayers grown over a 10 mm diameter field were exposed to clinically relevant proton doses ranging from 7 to 35 Gy at ultra-high instantaneous dose rates of 10 7 Gy/s. Dose-dependent cell survival measurements of human normal and tumor cells exposed to LD protons showed significantly higher cell survival of normal-cells compared to tumor-cells for total doses of 7 Gy and higher, which was not observed to the same extent for X-ray reference irradiations at clinical dose rates. These findings provide preliminary evidence that compact LD proton sources enable a new and promising platform for investigating the physical, chemical and biological mechanisms underlying the FLASH effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Transactive energy systems for distributed blackstart and service recovery

Abstract Current transactive controls use marginal benefits and marginal costs to achieve an economic market efficiency during normal grid operations. However, the transactive mechanisms designed for normal economic operations cannot be applied directly for the contingencies because the grid operations during contingencies are often dictated by technical needs rather than purely economic criteria. For instance, one of the key technical requirements for the blackstart is to have at least one blackstart capable resource cleared which cannot be ensured by the transactive mechanism designed for normal economic operations because they work primarily based on the marginal benefit and marginal cost of the participants. This article presents one of the first attempts to develop a transactive mechanism to be used during grid contingencies. A distributed blackstart and service recovery is used as an example contingency to evaluate the performance of the proposed transactive mechanism. The performance of the proposed transactive mechanism is demonstrated for various use cases using a modified IEEE‐123 node test system. The simulation results demonstrated the proof of concept of applying a transactive mechanism to enable distributed blackstart and service recovery by engaging the mix of blackstart capable and non‐capable distributed energy resources.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

ZTF SN Ia DR2 follow-up: Exploring the origin of the Type Ia supernova host galaxy step through Si II velocities

The relation between Type Ia supernovae (SNe Ia) and the stellar masses of their host galaxy is well documented. In particular, Hubble residuals display a distinct luminosity shift based on host mass. This is known as the mass step. This effect is widely used as an additional correction factor in the standardisation of SN Ia luminosities. We investigate the Hubble residuals and the mass step of normal SNe Ia in the context of Si IIλ6355 velocities based on 277 normal SNe Ia that are near their peak in the second data release (DR2) of the Zwicky Transient Facility (ZTF). We divided the sample into high-velocity (HV) and normal-velocity (NV) SNe Ia, separated at 12,000 km s −1 . This produced a sample of 70 HV and 207 NV objects. We then explored potential environment- and/or progenitor-related effects by investigating the Si IIλ6355 velocities with parameters such as the light-curve stretch x 1 , the colour c, and the host galaxy properties. Although we only find a marginal difference between the Hubble residuals of HV and NV SNe Ia, the NV mass step is 0.149 ± 0.024 mag (6.3σ). The HV mass step is smaller, 0.046 ± 0.041 mag (1.1σ), and is consistent with zero. The difference between the NV and HV mass steps is modest, at ∼2.2σ. Moreover, the clearest subtype difference appears for SNe in central regions (d DLR < 1), where NV SNe Ia show a large mass step, whereas HV SNe Ia are consistent with no step, yielding a difference of 3.1–3.6σ between NV and HV SNe Ia. We observe a host-colour step for both subtypes. NV SNe Ia show a step of 0.142 ± 0.024 mag (5.9σ), while HV SNe Ia show a step of 0.158 ± 0.042 mag (3.8σ), where the HV SNe Ia step appears to be larger, but the significance is lower because the sample size is smaller. Overall, the NV and HV colour steps are statistically consistent. HV SNe Ia also show modest (∼2.5–3σ) steps in certain subsets, such as those in outer regions (d DLR > 1), whereas NV SNe display stronger environmental trends. Our results indicate that NV SNe Ia appear to be more environmentally sensitive, particularly in central likely metal-rich and older regions, while HV SNe Ia show weaker and subset-dependent trends. This suggests that applying a universal mass-step correction might introduce biases, and that incorporating refined classifications and/or environment-dependent factors, such as the location within the host, might improve future cosmological analyses beyond the standard x 1 and c cuts.

supernovae: general↗

DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is γ τ A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of γ τ A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗