Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ROUGHNESS”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological implications from two decades of spectroscopic surveys at the Apache Point Observatory

We present the cosmological implications from final measurements of clustering using galaxies, quasars, and Lyα forests from the completed Sloan Digital Sky Survey (SDSS) lineage of experiments in large-scale structure. These experiments, composed of data from SDSS, SDSS-II, BOSS, and eBOSS, offer independent measurements of baryon acoustic oscillation (BAO) measurements of angular-diameter distances and Hubble distances relative to the sound horizon, r d , from eight different samples and six measurements of the growth rate parameter, fσ 8 , from redshift-space distortions (RSD). This composite sample is the most constraining of its kind and allows us to perform a comprehensive assessment of the cosmological model after two decades of dedicated spectroscopic observation. We show that the BAO data alone are able to rule out dark-energy-free models at more than eight standard deviations in an extension to the flat, ΛCDM model that allows for curvature. When combined with Planck Cosmic Microwave Background (CMB) measurements of temperature and polarization, under the same model, the BAO data provide nearly an order of magnitude improvement on curvature constraints relative to primary CMB constraints alone. Independent of distance measurements, the SDSS RSD data complement weak lensing measurements from the Dark Energy Survey (DES) in demonstrating a preference for a flat ΛCDM cosmological model when combined with Planck measurements. The combined BAO and RSD measurements indicate σ 8 =0.85±0.03, implying a growth rate that is consistent with predictions from Planck temperature and polarization data and with General Relativity. When combining the results of SDSS BAO and RSD, Planck, Pantheon Type Ia supernovae (SNe Ia), and DES weak lensing and clustering measurements, all multiple-parameter extensions remain consistent with a ΛCDM model. Regardless of cosmological model, the precision on each of the three parameters, Ω Λ , H 0 , and σ 8 , remains at roughly 1%, showing changes of less than 0.6% in the central values between models. In a model that allows for free curvature and a time-evolving equation of state for dark energy, the combined samples produce a constraint Ω K =-0.0022±0.0022. The dark energy constraints lead to w 0 =-0.909±0.081 and w a =$-0.49_{-0.30}^{+0.35}$, corresponding to an equation of state of w p =-1.018±0.032 at a pivot redshift z p =0.29 and a Dark Energy Task Force Figure of Merit of 94. The inverse distance ladder measurement under this model yields H 0 =68.18±0.79 km s -1 Mpc -1 , remaining in tension with several direct determination methods; the BAO data allow Hubble constant estimates that are robust against the assumption of the cosmological model. In addition, the BAO data allow estimates of H 0 that are independent of the CMB data, with similar central values and precision under a ΛCDM model. Our most constraining combination of data gives the upper limit on the sum of neutrino masses at Σm ν <0.115 eV (95% confidence). Finally, we consider the improvements in cosmology constraints over the last decade by comparing our results to a sample representative of the period 2000-2010. We compute the relative gain across the five dimensions spanned by w, Ω k , Σm ν , H 0 , and σ 8 and find that the SDSS BAO and RSD data reduce the total posterior volume by a factor of 40 relative to the previous generation. Finally, adding again the Planck, DES, and Pantheon SN Ia samples leads to an overall contraction in the five-dimensional posterior volume of 3 orders of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Does inhomogeneous big bang nucleosynthesis produce an inhomogeneous element distribution today?

Inhomogeneous big bang nucleosynthesis (BBN) produces a spatially inhomogeneous distribution of element abundances at $T \sim 10^9$ K , but subsequent element diffusion will tend to erase these inhomogeneities. In this work, we calculate the cosmological comoving diffusion length for the BBN elements. This diffusion length is limited by atomic scattering and is therefore dominated by diffusion when the atoms are neutral, between the redshifts of recombination and reionization. We find that the comoving diffusion length today is $d_{\text{com}} ≈ 70$ pc for all of the elements of interest except 7 Li , for which $d_{\text{com}}$ is an order of magnitude smaller because 7 Li remains ionized throughout the relevant epoch. This comoving diffusion length corresponds to a substellar baryonic mass scale and is roughly equal to the horizon scale at BBN. These results lend support to the possibility that inhomogeneities on scales larger than the horizon at BBN could lead to a spatially inhomogeneous distribution of elements today, while purely subhorizon fluctuations at BBN can result only in a homogeneous element distribution at present.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Nuclear two point correlation functions on a quantum computer

The calculation of dynamic response functions is expected to be an early application benefiting from rapidly developing quantum hardware resources. The ability to calculate real-time quantities of strongly correlated quantum systems is one of the most exciting applications that can easily reach beyond the capabilities of traditional classical hardware. Response functions of fermionic systems at moderate momenta and energies corresponding roughly to the Fermi energy of the system are a potential early application because the relevant operators are nearly local, and the energies can be resolved in moderately short real time, reducing the spatial resolution and gate depth required. This is particularly the case in quasielastic electron and neutrino scattering from nuclei, a topic of great interest in the nuclear and particle physics communities and directly related to experiments designed to probe neutrino properties. In this paper we use current quantum hardware and error mitigation protocols to calculate response functions for a highly simplified nuclear model through calculations of a 2-point real time correlation function for a modified Fermi-Hubbard model in two dimensions with three distinguishable nucleons on four lattice sites.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Atacama Cosmology Telescope: Constraints on prerecombination early dark energy

The early dark energy (EDE) scenario aims to increase the value of the Hubble constant (H0) inferred from cosmic microwave background (CMB) data over that found in the standard cosmological model (Λ CDM ), via the introduction of a new form of energy density in the early Universe. The EDE component briefly accelerates cosmic expansion just prior to recombination, which reduces the physical size of the sound horizon imprinted in the CMB. Previous work has found that nonzero EDE is not preferred by Planck CMB power spectrum data alone, which yield a 95% confidence level (C.L.) upper limit f EDE < 0.087 on the maximal fractional contribution of the EDE field to the cosmic energy budget. In this paper, we fit the EDE model to CMB data from the Atacama Cosmology Telescope (ACT) data release 4. We find that a combination of ACT, large-scale Planck TT (similar to WMAP), Planck CMB lensing, and BAO data prefers the existence of EDE at >99.7 % C .L . : f EDE = 0.091$^{+0.020}_{-0.036}$, with H 0 = 70.9$^{+1.0}_{-2.0}$ km/s/Mpc (both 68% C.L.). From a model-selection standpoint, we find that EDE is favored over Λ CDM by these data at roughly 3 σ significance. In contrast, a joint analysis of the full Planck and ACT data yields no evidence for EDE, as previously found for Planck alone. We show that the preference for EDE in ACT alone is driven by its TE and EE power spectrum data. The tight constraint on EDE from Planck alone is driven by its high-ℓ TT power spectrum data. Understanding whether these differing constraints are physical in nature, due to systematics, or simply a rare statistical fluctuation is of high priority. Here, the best-fit EDE models to ACT and Planck exhibit coherent differences across a wide range of multipoles in TE and EE, indicating that a powerful test of this scenario is anticipated with near-future data from ACT and other ground-based experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Neutron detection and application with a novel 3D-projection scintillator tracker in the future long-baseline neutrino oscillation experiments

Neutrino oscillation experiments require a precise measurement of the neutrino energy. However, the kinematic detection of the final-state neutron in the neutrino interaction is missing in current neutrino oscillation experiments. The missing neutron kinematic detection results in a smaller detected neutrino energy than the true neutrino energy. A novel 3D-projection scintillator tracker, which consists of roughly ten million active cubes covered with an optical reflector, is capable of measuring the neutron kinetic energy and direction on an event-by-event basis using the time-of-flight technique thanks to the fast timing, fine granularity, and high light yield. The $\overline{v}$ μ interactions tend to produce neutrons in the final state. By measuring the neutron kinetic energy, the $\overline{v}$ μ energy can be reconstructed better, allowing a tighter incoming neutrino flux constraint. This article shows the detector's ability to reconstruct neutron kinetic energy and the $\overline{v}$ μ flux constraint achieved by selecting the charged-current interactions without mesons or protons in the final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Standard model and the lattice

The S U ( 3 ) ⊗ S U ( 2 ) ⊗ U ( 1 ) standard model maps smoothly onto a conventional lattice gauge formulation, including the parity violation of the weak interactions. The formulation makes use of the pseudoreality of the weak group and requires the inclusion a full generation of both leptons and quarks. As in continuum discussions, chiral eigenstates of the Dirac operator generate known anomalies, although with rough gauge configurations these are no longer exact zero modes of the Dirac operator. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Atacama Cosmology Telescope: High-resolution component-separated maps across one third of the sky

Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-y distortion due to the thermal Sunyaev-Zel’dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multiwavelength observations to spectrally isolate one component. Here, in this work, we present a new arc-minute-resolution Compton-y map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 deg 2 ). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) data release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the Planck satellite at frequencies between 30 and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing Planck component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Toward a neutrino-limited dark matter search with crystalline xenon

Experiments searching for weakly interacting massive particle dark matter are now detecting background events from solar neutrino-electron scattering. However, the dominant radioactive background in state-of-the-art experiments such as LZ and XENONnT is beta decays from radon contamination. In spite of careful detector material screening, radon progenitor atoms are ubiquitous and long-lived, and radon is extremely soluble in liquid xenon. We propose a change of phase and demonstrate that crystalline xenon offers more than a factor × 500 exclusion against radon ingress, compared with the liquid state. This level of radon exclusion would allow crystallized versions of existing experiments to probe spin-independent cross sections near 10 − 47 cm 2 in roughly 11 years, as opposed to the 35 years required otherwise. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Constraints on metastable superheavy dark matter coupled to sterile neutrinos with the Pierre Auger Observatory

Dark matter particles could be superheavy, provided their lifetime is much longer than the age of the Universe. Using the sensitivity of the Pierre Auger Observatory to ultrahigh energy neutrinos and photons, we constrain a specific extension of the Standard Model of particle physics that meets the lifetime requirement for a superheavy particle by coupling it to a sector of ultralight sterile neutrinos. Our results show that, for a typical dark coupling constant of 0.1, the mixing angle θ m between active and sterile neutrinos must satisfy, roughly, θ m ≲ 1.5 × 10 − 6 ( M X / 10 9 GeV ) − 2 for a mass M X of the dark-matter particle between 10 8 GeV and 10 11 GeV . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Cosmology From CMB Lensing and Delensed EE Power Spectra Using 2019-2020 SPT-3G Polarization Data

From CMB polarization data alone we reconstruct the CMB lensing power spectrum, comparable in overall constraining power to previous temperature-based reconstructions, and an unlensed E -mode power spectrum, with clear detections of the third through tenth acoustic peaks. The observations, taken in 2019 and 2020 with the South Pole Telescope (SPT) and the SPT-3G camera, cover 1500 deg 2 at 95, 150, and 220 GHz with arcminute resolution and roughly 4.9 µ K-arcmin coadded noise in polarization. The power spectrum estimates, together with systematic parameter estimates and a joint covariance matrix, follow from a Bayesian analysis using the Marginal Unbiased Score Expansion (MUSE) method. The E -mode spectrum at ℓ > 2000 and lensing spectrum at L > 350 are the most precise to date. Assuming the ΛCDM model, and using only these SPT data and priors on τ and absolute calibration from Planck, we find H 0 = 66.81 ± 0.81 km/s/Mpc, comparable in precision to the Planck determination and in 5.4 σ tension with the most precise H 0 inference derived via the distance ladder. We also find S 8 ≡ σ 8 (Ω m /0.3) 0.5 = 0.850 ± 0.017, providing further independent evidence of a slight tension with low-redshift structure probes. The ΛCDM model provides a good simultaneous fit to the combined Planck, ACT, and SPT data, and thus passes a powerful test. Combining these CMB datasets with BAO observations, we explore extensions to the ΛCDM model. We find that the effective number of neutrino species, spatial curvature, and primordial helium fraction are consistent with standard model values, and that the 95% confidence upper limit on the neutrino mass sum is 0.075 eV, close to the minimum sum expected from observations of solar and atmospheric neutrino oscillations. The SPT data are consistent with the somewhat weak (< 3 σ ) preference for excess lensing power seen in Planck and ACT data relative to predictions of the ΛCDM model given the combined Planck, ACT, and BAO data sets. Finally, we also detect at greater than 3 σ the influence of non-linear evolution in the CMB lensing power spectrum and discuss it in the context of the S 8 tension. Forthcoming SPT-3G analyses will feature deeper and wider observations in temperature and polarization, providing even tighter constraints and more powerful tests of the ΛCDM model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Hadronic vacuum polarization for the muon 𝑔 − 2 from lattice QCD: Complete short and intermediate windows

We present complete results for the hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment 𝑎 𝜇 in the short- and intermediate-distance window regions, which account for roughly 10% and 35% of the total HVP contribution to 𝑎 𝜇 , respectively. In particular, we perform lattice-QCD calculations for the isospin-symmetric connected and disconnected contributions, as well as corrections due to strong-isospin breaking. For the short-distance window observables, we investigate the so-called log-enhancement effects as well as the significant oscillations associated with staggered quarks in this region. For the dominant, isospin-symmetric light-quark-connected contribution, we obtain 𝑎$^{𝑙⁢𝑙,SD}_{𝜇⁡}$(conn) = 48.139⁢(11) stat ⁢(91) syst ⁢[92] total × 10 −10 and 𝑎$^{𝑙⁢𝑙,W}_{𝜇}$⁡(conn) = 206.90⁢(14) stat ⁢(61) syst ⁢[63] total × 10 −10 . We use Bayesian model averaging to fully estimate the covariance matrix between the individual contributions. Our determinations of the complete window contributions are 𝑎$^{SD}_{𝜇}$ = 69.05⁢(1) stat⁢ (21) syst ⁢[21] total × 10 −10 and 𝑎$^{W}_{𝜇}$ = 236.45⁢(17) stat⁢ (83) syst ⁢[85] total × 10 −10 . This work is part of our ongoing effort to compute all contributions to HVP with an overall uncertainty at the few-permille level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Measurements of the temperature and $E$-mode polarization of the cosmic microwave background from the full 500-square-degree SPTpol dataset

Using the full four-year SPTpol 500 deg 2 dataset in both the 95 and 150 GHz frequency bands, we present measurements of the temperature and E-mode polarization of the cosmic microwave background (CMB), as well as the E-mode polarization autopower spectrum (EE) and temperature-E-mode cross-power spectrum (TE) in the angular multipole range 50 < ℓ < 8000. We find the SPTpol dataset to be self-consistent, passing several internal consistency tests based on maps, frequency bands, bandpowers, and cosmological parameters. The full SPTpol dataset is well-fit by the ΛCDM model, for which we find H 0 = 70.48 ± 2.16 km s -1 Mpc -1 and Ω m = 0.271 ± 0.026, when using only the SPTpol data and a Planck-based prior on the optical depth to reionization. The ΛCDM parameter constraints are consistent across the 95 GHz-only, 150 GHz-only, TE-only, and EE-only data splits. Between the ℓ < 1000 and ℓ > 1000 data splits, the ΛCDM parameter constraints are borderline consistent at the ∼2σ level. This consistency improves when including a parameter A L , the degree of lensing of the CMB inferred from the smearing of acoustic peaks. When marginalized over A L , the ΛCDM parameter constraints from SPTpol are consistent with those from Planck. In conclusion, the power spectra presented here are the most sensitive measurements of the lensed CMB damping tail to date for roughly ℓ > 1700 in TE and ℓ > 2000 in EE.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Integral equation theory of thermodynamics, pair structure, and growing static length scale in metastable hard sphere and Weeks-Chandler-Andersen fluids

Here, we employ the Ornstein-Zernike integral equation theory with the Percus-Yevick (PY) and modified-Verlet (MV) closures to study the equilibrium structural and thermodynamic properties of metastable monodisperse hard sphere and continuous repulsion Weeks-Chandler-Andersen (WCA) fluids under density and temperature conditions where the system is strongly overcompressed or supercooled, respectively. The theoretical results are compared to crystal-avoiding simulations of these dense monodisperse model one-component fluids. The equation of state (EOS) and dimensionless compressibility are computed using both the virial and compressibility routes. For hard spheres, the MV-based virial route EOS and dimensionless compressibility are in very good agreement with simulation for all packing fractions, much better than the PY analogs. The corresponding MV-based predictions for the static structure factor are also very good. The amplitude of density fluctuations on the local cage scale and in the long wavelength limit, and three technically different measures of the density correlation length, are studied with both closures. All five properties grow in a roughly exponential manner with density in the metastable regime up to packing fractions of 58% with no sign of saturation. The MV-based results are in good agreement with our crystal-avoiding simulations. Interestingly, the density dependences of long and short wavelength quantities are closely related. The MV-based theory is also quite accurate for the thermodynamics and structure of supercooled monodisperse WCA fluids. Overall our findings are also relevant as critical input to microscopic theories that relate the equilibrium pair correlation function or static structure factor to dynamical constraints, barriers, and activated relaxation in glass-forming liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Applicability of semiclassical theories in the strong-field plasma regime

For many purposes, classical plasma dynamics models can work surprisingly well, even for strong electromagnetic fields, approaching the Schwinger critical fields, and high frequencies, approaching the Compton frequency. However, the applicability of classical models tends to depend rather sensitively on the details of the problem. In the present paper, we study the specific case of plasma oscillations to draw a line between the classical and quantum relativistic regimes. Here, due to the field geometry of study, mechanisms like radiation reaction and Breit-Wheeler pair production, which tend to be important for electromagnetic fields, are rather effectively suppressed. Moreover, we find that the polarization current due to the electron spin is generally negligible for frequencies below the Compton frequency, compared with the free current, whose magnitude is well-approximated by the classical Vlasov theory. However, we show that pair creation due to the Schwinger mechanism can sometimes be important for surprisingly modest field strengths, of the order of 10% of the critical field or even smaller. A rough guideline for when the classical Vlasov theory can be applied is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Paradigm for global gyrokinetic turbulence

Using a novel wave-number-advection algorithm, we describe a method to add global profile curvature terms to the local gyrokinetic equations. This approach enables a high-precision global simulation capability without sacrificing the efficiency and spectral accuracy of local simulations. Preliminary numerical results suggest that there are two distinct types of global correction to the local result: (1) a trivial shift of the local result in radius, (2) a non-trivial decorrelation that is roughly equivalent to the well-known E×B shear effect. This interpretation provides a new, intuitive framework for understanding global gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
DOE OSTI2023

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗