Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Unity”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Measurement of K 0 S and K *0 in p+p, d+Au, and Cu+Cu collisions at sqrt S NN = 200 GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider has performed a systematic study of K 0 S and K *0 meson production at midrapidity in p+p, d+Au, and Cu+Cu collisions at sqrt S NN = 200 GeV. The K 0 S and K *0 mesons are reconstructed via their K 0 S and π 0 (→γγ)π 0 (→γγ) and K *0 → K ± π ± decay modes, respectively. The measured transverse-momentum spectra are used to determine the nuclear modification factor of K 0 S and K *0 mesons in d+Au and Cu+Cu collisions at different centralities. In the d+Au collisions, the nuclear modification factor of K 0 S and K *0 mesons is almost constant as a function of transverse momentum and is consistent with unity showing that cold-nuclear-matter effects do not play a significant role in the measured kinematic range. In Cu+Cu collisions, within the uncertainties no nuclear modification is registered in peripheral collisions. In central collisions, both mesons show suppression relative to the expectations from the p+p yield scaled by the number of binary nucleon-nucleon collisions in the Cu+Cu system. In the p T range 2–5 GeV/c, the strange mesons ( K 0 S , K *0 ) similarly to the Φ meson with hidden strangeness, show an intermediate suppression between the more suppressed light quark mesons (π 0 ) and the nonsuppressed baryons (p, p-bar). At higher transverse momentum, p T > 5 GeV/c, production of all particles is similarly suppressed by a factor of ≈2. (auth)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Measurement of the EMC effect in the deuteron

We have determined the structure function ratio $R^d_{\rm EMC}=F_2^d/(F_2^n+F_2^p)$ from recently published $F_2^n/F_2^d$ data taken by the BONuS experiment using CLAS at Jefferson Lab. This ratio deviates from unity, with a slope $dR_{\rm EMC}^{d}/dx= -0.10\pm 0.05$ in the range of Bjorken $x$ from 0.35 to 0.7, for invariant mass $W>1.4$ GeV and $Q^2>1$ GeV$^2$. The observed EMC effect for these kinematics is consistent with conventional nuclear physics models that include off-shell corrections, as well as with empirical analyses that find the EMC effect proportional to the probability of short-range nucleon-nucleon correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Flavor dependence of the pion and kaon form factors and parton distribution functions

The separate quark flavor contributions to the pion and kaon valence quark distribution functions are studied, along with the corresponding electromagnetic form factors in the space-like region. The calculations are made using the solution of the Bethe-Salpeter equation for the model of Nambu and Jona-Lasinio with proper-time regularization. Both the pion and kaon form factors and the valence quark distribution functions reproduce many features of the available empirical data. The larger mass of the strange quark naturally explains the empirical fact that the ratio u(K) + (x)/u(pi) + (x) drops below unity at large x, with a value of approximately M$2\atop{u}$/Ms$2\atop{s}$ as x → 1. With regard to the elastic form factors we report a large flavor dependence, with the u-quark contribution to the kaon form factor being an order of magnitude smaller than that of the s-quark at large Q 2 , which may be a sensitive measure of confinement effects in QCD. Surprisingly though, the total K + and π + form factors differ by only 10%. Lastly, in general we find that flavor breaking effects are typically around 20%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Charge symmetry breaking effects in pion and kaon structure

Charge symmetry breaking (CSB) effects associated with the u and d quark mass difference are investigated in the quark distribution functions and spacelike electromagnetic form factors of the pion and kaon. We use a confining version of the Nambu-Jona-Lasimo model, where CSB effects at the infrared scale associated with the model are driven by the dressed u and d quark mass ratio, which because of dynamical chiral symmetry breaking is much closer to unity than the associated current quark mass ratio. The pion and kaon are given as bound states of a dressed quark and a dressed antiquark governed by the Bethe-Salpeter equation, and exhibit the properties of Goldstone bosons, with a pion mass difference given by m(Pi)+(2) - m(Pi)0(2) (m(u) - m(d))(2) as demanded by dynamical chiral symmetry breaking. We find significant CSB effects for realistic current quark mass ratios (m(u)/m(d) similar to 0.5) in the quark flavor-sector electromagnetic form factors of both the pion and kaon. For example, the difference between the u and d quark contributions to the Pi(+) electromagnetic form factors is about 8% at a momentum transfer of Q(2) = 10 GeV2, while the analogous effect for the light quark sector form factors in the K+ and K-0 is about twice as large. For the parton distribution functions we find CSB effects which are considerably smaller than those found in the electromagnetic form factors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Distortions in the surface of last scattering

The surface of last scattering of the photons in the cosmic microwave background is not a spherical shell. Apart from its finite width, each photon experiences a different gravitational potential along its journey to us, leading to different travel times in different directions. Photons passing through an underdense region started farther away from us than those traveling through an overdense region. Thus, the surface of last scattering is corrugated, a deformed spherical shell. Here, we present an estimator quadratic in the temperature and polarization fields that could provide a map of the time delays as a function of position on the sky. The signal to noise of this map could exceed unity for the dipole, thereby providing a rare insight into the universe on the largest observable scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Relic cosmological vector fields and inflationary gravitational waves

Here we show that relic vector fields can significantly impact a spectrum of primordial gravitational waves in the postinflationary era. We consider a triplet of U(1) fields in a homogeneous, isotropic configuration. The interaction between the gravitational waves and the vector fields, from the end of reheating to the present day, yields novel spectral features. The amplitude, tilt, shape, and net chirality of the gravitational wave spectrum are shown to depend on the abundance of the electric- and magneticlike vector fields. Our results show that even a modest abundance can have strong implications for efforts to detect the imprint of gravitational waves on the cosmic microwave background polarization. We find that a vector field comprising less than 2% of the energy density during the radiation-dominated era can have a greater than order unity effect on the predicted inflationary gravitational wave spectrum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Gravitational waves from fully general relativistic oscillon preheating

As long-lived quasisolitons from the fragmentation of a scalar condensate, oscillons may dominate the preheating era after inflation. During this period, stochastic gravitational waves can also be generated. Here we quantify the gravitational-wave production in this period with simulations accounting for full general relativity to capture all possible nonperturbative effects. We compute the gravitational-wave spectra across a range of choices of the oscillon preheating models and compare our results to a conventional perturbative approach on a Friedmann-Lemaître-Robertson-Walker (FLRW) background. We clarify the gauge ambiguities in computing induced gravitational waves from scenarios where dense nonperturbative objects such as oscillons are being formed. In particular, we find that the synchronous gauge tends to contain large artificial enhancements in the gravitational-wave spectrum due to gauge modes if gravity plays an important role in the formation of the oscillons, while other gauge choices, such as the radiation gauge or a suitably chosen “1+log” gauge, can efficiently reduce the contributions of gauge modes. The full general-relativistic simulations indicate that gravitational-wave spectra obtained from the perturbative approach on the FLRW background are fairly accurate, except when oscillon formation induces strong gravitational effects, for which case there can be an order unity enhancement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Measurements of the branching fractions B ( B ¯ 0 → D * + π - ) and B ( B ¯ 0 → D * + K - ) and tests of QCD factorization

Using (771.6 ± 10.6) × 10 6 $B\overline{B}$ meson pairs recorded by the Belle experiment at the KEKB e + e - collider, we report the branching fractions $\mathscr{B}$($\overline{B}$ 0 → D* + π - ) = (2.62 ± 0.02 ± 0.09 X 10 -3 and $\mathscr{B}$($\overline{B}$ 0 → D* + K - ) = (2.22 ± 0.06 ± 0.08) × 10 -4 ; the quoted uncertainties are statistical and systematic, respectively. A measurement of the ratio of these branching fractions is also presented, R K/π = $\mathscr{B}$($\overline{B}$ → D* + K - )/ $\mathscr{B}$($\overline{B}$ → D* + π - ) = (8.41 ± 0.24 ± 0.13) × 10 -2 , where systematic uncertainties due to the D* + meson reconstruction cancel out. Furthermore, we report a new QCD factorization test based on the measured ratios for $\overline{B}$ → D* + h - and $\overline{B}$ → D* + ℓ - $v$ decays at squared momentum transfer values equivalent to the mass of the h = π or K hadron. The parameters |a 1 (h)| are measured to be |a 1 (π)| = 0.884 ± 0.004 ± 0.003 ± 0.016 and |a 1 (K)| = 0.913 ± 0.019 ± 0.008 ± 0.013, where the last uncertainties account for all external inputs. These values are approximately 15% lower than those expected from theoretical predictions. Subsequently, flavor SU(3) symmetry is tested by measuring the ratios for pions and kaons, |a 1 (K)| 2 /|a 1 (π) 2 | = 1.066 ± 0.042 ± 0.018 ± 0.023, as well as for different particle species. The ratio is consistent with unity and therefore no evidence for SU(3) symmetry breaking effects is found at the 5% precision level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Current data are consistent with flat spatial hypersurfaces in the Λ CDM cosmological model but favor more lensing than the model predicts

Here, we study the performance of three pairs of tilted, and a pair of untilited, ΛCDM cosmological models, with three of these four pairs allowing for non-flat spatial hypersurfaces, against cosmic microwave background (CMB) temperature and polarization power spectrum data (P18), measurements of the Planck 2018 lensing potential power spectrum (lensing), and a large compilation of non-CMB data (non-CMB). For the eight models, we measure cosmological parameters and study whether or not pairs of the data sets (as well as subsets of them) are mutually consistent in these models. Half of these models allow the lensing consistency parameter A L , which re-scales the gravitational potential power spectrum, to be an additional free parameter to be determined from data, while the other three have A L = 1 which is the theoretically expected value. The pair of untilted non-flat ΛCDM models are incompatible with P18 data. The tilted spatially-flat models assume the usual primordial spatial inhomogeneity power spectrum that is a power law in wave number. The tilted non-flat models assume either the primordial power spectrum used in the Planck group anal yses [Planck P(q)], that has recently been numerically shown to be a good approximation to what is quantum-mechanically generated from a particular choice of closed inflation model initial conditions, or a recently computed power spectrum [new P(q)] that quantum-mechanically follows from a different set of non-flat inflation model initial conditions. In the tilted non-flat models with A L = 1 we find differences between P18 data and non-CMB data cosmological parameter constraints, which are large enough to rule out the Planck P(q) model at 3σ but not the new P(q) model. No significant differences are found when cosmological parameter constraints obtained with two different data sets are compared within the standard tilted flat ΛCDM model. While both P18 data and non-CMB data separately favor a closed geometry, with spatial curvature density parameter Ω k < 0, when P18+non-CMB data are jointly analyzed the evidence in favor of non-flat hypersurfaces subsides. Differences between P18 data and non-CMB data cosmological constraints subside when A L is allowed to vary. From the most restrictive P18+lensing+non-CMB data combination we get almost model-independent constraints on the cosmological parameters and find that the A L > 1 option is preferred over the Ω k < 0 one, with the A L parameter, for all models, being larger than unity by ~ 2.5σ. According to the deviance information criterion, in the P18+lensing+non-CMB analysis, the varying A L option is on the verge of being strongly favored over the A L = 1 one, which could indicate a problem for the standard tilted flat ΛCDM model. These data are consistent with flat spatial hypersurfaces but more and better data could improve the constraints on Ω k and might alter this conclusion. Error bars on some cosmological parameters are significantly reduced when non-CMB data are used jointly with P18+lensing data. For example, in the tilted flat ΛCDM model for P18+lensing+non-CMB data the Hubble constant H 0 = 68.09 ± 0.38 km s -1 Mpc -1 , which is consistent with that from a median statistics analysis of a large compilation of H 0 measurements as well as with a number of local measurements of the cosmological expansion rate. This H 0 error bar is 31% smaller than that from P18+lensing data alone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Illuminating all-hadronic final states with a photon: Exotic decays of the Higgs boson to four bottom quarks in vector boson fusion plus gamma at hadron colliders

We investigate the potential to detect Higgs boson decays to four bottom quarks through a pair of pseudoscalars, a final state that is predicted by many theories beyond the Standard Model. For the first time, the signal sensitivity is evaluated for the final state using the vector boson fusion (VBF) production with and without an associated photon, for the Higgs at m H = 125 GeV , at hadron colliders. The signal significance is 0.5 σ to 6 σ , depending on the pseudoscalar mass m a , when setting the Higgs decay branching ratio to unity, using an integrated luminosity of 150 fb − 1 at s = 13 TeV . This corresponds to an upper limit of 0.3, on the Higgs branching ratio to four bottom quarks, with a nonobservation of the decay. We also consider several variations of selection requirements—input variables for the VBF tagging and the kinematic variables for the photon—that could help guide the design of new triggers for the Run-3 period of the LHC and for the HL-LHC. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Particle production as a function of charged-particle flattenicity in 𝑝⁢𝑝 collisions at $\sqrt{s}$ = 13 TeV

This paper reports the first measurement of the transverse momentum (𝑝 T ) spectra of primary charged pions, kaons, (anti)protons, and unidentified particles as a function of the charged-particle flattenicity in pp collisions at $\sqrt{s}$ =13 TeV. Flattenicity is a novel event shape observable that is measured in the pseudorapidity intervals covered by the V0 detector, 2.8 < 𝜂 < 5.1 and −3.7 < 𝜂 < −1.7. According to QCD-inspired phenomenological models, it shows sensitivity to multiparton interactions and is less affected by biases toward larger 𝑝 T due to local multiplicity fluctuations in the V0 acceptance than multiplicity. The analysis is performed in minimum-bias (MB) as well as in high-multiplicity events up to 𝑝 T = 20 GeV/⁢𝑐. The event selection requires at least one charged particle produced in the pseudorapidity interval |𝜂| < 1. The measured 𝑝 T distributions, average 𝑝 T , kaon-to-pion and proton-to-pion particle ratios, presented in this paper, are compared to model calculations using pythia 8 based on color strings and EPOS LHC. The modification of the 𝑝 T -spectral shapes in low-flattenicity events that have large event activity with respect to those measured in MB events develops a pronounced peak at intermediate 𝑝 T (2 < 𝑝 T < 8 GeV/𝑐), and approaches the vicinity of unity at higher 𝑝 T . The results are qualitatively described by pythia, and they show different behavior than those measured as a function of charged-particle multiplicity based on the V0M estimator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Observation of D → K ¯ 1 ( 1270 ) μ + ν μ and test of lepton flavor universality with D → K ¯ 1 ( 1270 ) ℓ + ν ℓ

By analyzing 7.93 fb − 1 of e + e − collision data collected at the center-of-mass energy of 3.773 GeV with the BESIII detector operated at the BEPCII collider, we report the observation of the semimuonic decays of D + → K ¯ 1 ( 1270 ) 0 μ + ν μ and D 0 → K 1 ( 1270 ) − μ + ν μ with statistical significances of 12.5 σ and 6.0 σ , respectively. Their decay branching fractions are determined to be B [ D + → K ¯ 1 ( 1270 ) 0 μ + ν μ ] = ( 2.36 ± 0.2 0 − 0.27 + 0.18 ± 0.48 ) × 10 − 3 and B [ D 0 → K 1 ( 1270 ) − μ + ν μ ] = ( 0.78 ± 0.1 1 − 0.09 + 0.05 ± 0.15 ) × 10 − 3 , where the first and second uncertainties are statistical and systematic, respectively, and the third originates from the input branching fraction of K ¯ 1 ( 1270 ) 0 → K − π + π 0 or K 1 ( 1270 ) − → K − π + π − . Combining our branching fractions with the previous measurements of B [ D + → K ¯ 1 ( 1270 ) 0 e + ν e ] and B [ D 0 → K 1 ( 1270 ) − e + ν e ] , we determine the branching fraction ratios to be B [ D + → K ¯ 1 ( 1270 ) 0 μ + ν μ ] / B [ D + → K ¯ 1 ( 1270 ) 0 e + ν e ] = 1.03 ± 0.14 − 0.15 + 0.11 and B [ D 0 → K 1 ( 1270 ) − μ + ν μ ] / B [ D 0 → K 1 ( 1270 ) − e + ν e ] = 0.74 ± 0.13 − 0.13 + 0.08 . Using the branching fractions measured in this work and the world-average lifetimes of the D + and D 0 mesons, we determine the semimuonic partial decay width ratio to be Γ [ D + → K ¯ 1 ( 1270 ) 0 μ + ν μ ] / Γ [ D 0 → K 1 ( 1270 ) − μ + ν μ ] = 1.22 ± 0.10 − 0.09 + 0.06 , which is consistent with unity as predicted by isospin conservation. Published by the American Physical Society 2025

Ablikim, M.↗
DOE OSTI2016

Unitarity and the three flavor neutrino mixing matrix

Unitarity is a fundamental property of any theory required to ensure we work in a theoretically consistent framework. In comparison with the quark sector, experimental tests of unitarity for the 3x3 neutrino mixing matrix are considerably weaker. It must be remembered that the vast majority of our information on the neutrino mixing angles originates from v - e and v μ disappearance experiments, with the assumption of unitarity being invoked to constrain the remaining elements. New physics can invalidate this assumption for the 3x3 subset and thus modify our precision measurements. We also perform a reanalysis to see how global knowledge is altered when one refits oscillation results without assuming unitarity, and present 3σ ranges for allowed U PMNS elements consistent with all observed phenomena. We calculate the bounds on the closure of the six neutrino unitarity triangles, with the closure of the v - e and v μ triangle being constrained to be ≤0.03, while the remaining triangles are significantly less constrained to be ≤ 0.1 - 0.2. Similarly for the row and column normalization, we find their deviation from unity is constrained to be ≤ 0.2 - 0.4, for four out of six such normalizations, while for the v μ and v e row normalization the deviations are constrained to be ≤0.07, all at the 3σCL. Additionally, we emphasize that there is significant room for new low energy physics, especially in the v τ sector which very few current experiments constrain directly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Light flavor-singlet scalars and walking signals in N f = 8 QCD on the lattice

Based on the highly improved staggered quark action, we perform in this paper lattice simulations of N f = 8 QCD and confirm our previous observations, both of a flavor-singlet scalar meson (denoted as σ) as light as the pion and of various “walking signals” through the low-lying spectra, with higher statistics, smaller fermion masses m f , and larger volumes. We measure M π , F π , M ρ , M a 0 , M a 1 , M b 1 , M N , M σ , F σ , (φφ) (both directly and through the Gell-Mann-Oakes-Renner relation), and the string tension. The data are consistent with the spontaneously broken phase of the chiral symmetry, in agreement with the previous results: Ratios of the quantities to M π monotonically increase in the smaller m f region towards the chiral limit similarly to N f = 4 QCD, in sharp contrast to N f = 12 QCD where the ratios become flattened. We perform fits to chiral perturbation theory, with the value of F π found in the chiral limit extrapolation reduced dramatically to roughly 2/3 of the previous result, suggesting the theory is much closer to the conformal window. In fact, each quantity obeys the respective hyperscaling relation throughout a more extensive m f region compared with earlier works. The hyperscaling relation holds with roughly a universal value of the anomalous dimension, γ m ≃ 1, with the notable exception of M π with γ m ≃ 0.6 as in the previous results, which reflects the above growing up of the ratios towards the chiral limit. This is a salient feature (walking signal) of N f = 8, unlike either N f = 4, which has no hyperscaling relation at all, or N f = 12 QCD, which exhibits universal hyperscaling. The effective γ m Ξ γ m (m f ) of M π defined for each m f region has a tendency to grow towards unity near the chiral limit, in conformity with the Nambu-Goldstone boson nature, as opposed to the case of N f = 12 QCD where it is almost constant. We further confirm the previous observation of the light σ with mass comparable to the pion in the studied m f region. In a chiral limit extrapolation of the σ mass using the dilaton chiral perturbation theory and also using the simple linear fit, we find the value consistent with the 125 GeV Higgs boson within errors. Finally, our results suggest that the theory could be a good candidate for walking technicolor model, having anomalous dimension γ m ≃ 1 and a light flavor-singlet scalar meson as a technidilaton, which can be identified with the 125 GeV composite Higgs in the N f = 8 one-family model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Primordial black holes and local non-Gaussianity in canonical inflation

Primordial black holes (PBHs) cannot be produced abundantly enough to be the dark matter in canonical single-field inflation under slow roll. This conclusion is robust to local non-Gaussian correlations between long- and short-wavelength curvature modes, which we show have no effect in slow roll on local primordial black hole abundances. For the prototypical model which evades this no go, ultra-slow roll (USR), these squeezed non-Gaussian correlations have at most an order unity effect on the variance of PBH-producing curvature fluctuations for models that would otherwise fail to form sufficient PBHs. Moreover, the transition out of USR, which is necessary for a success fulmodel, suppresses even this small enhancement unless it causes a large increase in the inflaton kinetic energy in a fraction of ane-fold, which we call a large and fast transition. Along the way we apply the in-in formalism, the $δN$ formalism, and gauge transformations to compute non-Gaussianities and illuminate different aspects of the physical origin of these results. Local non-Gaussianity in the squeezed limit does not weaken the Gaussian conclusion that PBHs as dark matter in canonical single-field inflation require a complicated and fine-tuned potential shape with an epoch where slow roll is transiently violated

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.

58 GEOSCIENCES↗
DOE OSTI2022

Mean-field model for the Curie-Weiss temperature dependence of coherence length in metallic liquids

The coherence length of the medium-range order (MRO) in metallic liquids is known to display a Curie-Weiss temperature dependence; its inverse is linearly related to temperature, and when extrapolated from temperatures above the glass transition, the coherence length diverges at a negative temperature with a critical exponent of unity. We propose a mean-field pseudospin model that explains this behavior. Specifically, we model the atoms and their local environment as Ising spins with antiferromagnetic exchange interactions. We further superimpose an exchange interaction between dynamical heterogeneities, or clusters of atoms undergoing cooperative motion. The coherence length in the metallic liquid is thus the correlation length between dynamical heterogeneities. Overall, our results reaffirm the idea that the MRO coherence length is a measure of point-to-set correlations, and that local frustrations in the interatomic interactions are prominent in metallic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Turbulence generation by shock interaction with a highly nonuniform medium

An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, $\sigma_\rho/\bar{\rho}\sim1$) and strong shocks (shock Mach numbers $\bar{M}_s\gtrsim10$). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗