Search NASASearch

SEARCH · Search NASA

Results for “sound velocity anomalies”

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.

160 records · Page 8

Characterization of Flashback and Flame-Holding in a Jet-in-Crossflow Mixing Configuration With Methane–Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence diagnostics techniques. The JICF burner has an optically accessible premixing tube allowing the optical diagnostics. The testing was conducted for varied premixer velocities (V) and equivalence ratio (ϕ) for 90–100% (H 2 , by mole) H 2 /CH 4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashback events were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower momentum flux ratio which bends the jet sharply, closer toward the injection plane. The mean OH-PLIF images characterized the flame-holding behavior, where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port. A large eddy simulation (LES) with detailed chemistry combustion modeling approach was implemented along with an OH* submechanism, and it showed qualitative agreement with the integrated line-of-sight chemiluminescence results as well as the planar OH-PLIF measurements.

chemiluminescence diagnostic

Eddy Covariance Theory: A Review

Eddy covariance (EC), the gold standard for measuring ecosystem scale gas and heat exchanges, has transformed our understanding of the breathing of the biosphere, and thus global change biology. Despite numerous methodological improvements and insights gained from the technique, the community faces persistent challenges that have been present since the first EC measurements. Here, we review the theoretical developments underpinning EC. We present theoretical developments in four important areas that have relevance to EC measurements of the net ecosystem exchanges (NEE) of gases and heat from a single tower: (i) measuring the total vertical flux density, (ii) flux attenuation, (iii) coordinate rotations, and (iv) energy balance closure. Persistent problems with EC measurements, such as the inability to close the energy budget, led us to identify two priorities for revisiting the theory underlying: (i) sensible heat flux calculations, and (ii) constraining the mean vertical wind velocity. We present a framework for improved calculation of sensible heat flux derived from first principles of fluid mechanics and thermodynamics that considers coupled heat and mass transfer so that conservation of both is obeyed. These refinements are motivated by the need for unbiased measurements of energy and mass transfer between the land surface and atmosphere for ecosystem research and to validate satellite observations and land surface models.

ecosystem fluxes

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES

Plasma flow generation and particle acceleration from expanding magnetic bubbles

Impulsive plasma dynamics in the laboratory are often driven by rising electric currents, yet their quantitative plasma response has not been well established. By means of fully kinetic particle-in-cell simulations and laser-driven capacitor-coil experiments, we show that a rising current expels plasma, forming an expanding magnetic bubble and accelerating particles. The expansion front velocity scales with the Alfvén speed determined by the magnetic field at its inner edge and the plasma density at its outer edge. This mechanism establishes impulsive current drive as a fundamental way that generates plasma flows and accelerates particles in laboratory plasmas, with potential relevance to astrophysics.

Zhang, Yang [Princeton University, NJ (United Stat

Λ polarization from vortex rings as the medium response for jet thermalization

We performed a systematic study on the formation of vorticity rings as the process for jet thermalization in the medium created in high-energy nuclear collisions. In this work, we expanded our previous analysis to a more realistic framework by considering noncentral events and fluctuations in the initial condition. We simulate the formation and evolution of the flow vortex structure in a relativistic viscous hydrodynamic model and study the sensitivity of the proposed “ring observable” (ℛ$^{𝑡}_{Λ}$) that can be measured experimentally through the polarization of Λ hyperons. We show that this observable is robust with respect to fluctuating initial conditions to capture the jet-induced vortex flow signal and further study its dependence on different model parameters, such as the jet's velocity, position, the fluid's shear viscosity, and the collision centrality. The proposed observable is associated with the formation of vorticity in a quark-gluon plasma, showing that the measurement of particle polarization can be a powerful tool to probe different properties of jet-medium interactions and to understand better the polarization induced by the transverse and longitudinal expansions of the medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models

The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.

Automation

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]

Lubricating Properties of Ceramic-Bonded Calcium Fluoride Coatings on Nickel-Base Alloys From 75° to 1900° F

The endurance life and the friction coefficient of ceramic-bonded calcium fluoride (CaF2) coatings on nickel-base alloys were determined at temperatures from 75° to 1900° F. The specimen configuration consisted of a hemispherical rider (3/16-in. rad.) sliding against the flat surface of a rotating disk. Increasing the ambient temperature (up to 1500° F) or the sliding velocity generally reduced the friction coefficient and improved coating life. Base-metal selection was critical above 1500° F. For instance, cast Inconel sliding against coated Inconel X was lubricated effectively to 1500° F, but at 1600° F severe blistering of the coatings occurred. However, good lubrication and adherence were obtained for René 41 sliding against coated René 41 at temperatures up to 1900° F; no blisters developed, coating wear life was fairly good, and the rider wear rate was significantly lower than for the unlubricated metals . Friction coefficients were 0.12 at 1500° F, 0.15 at 1700° F, and 0.17 at 1800° and 1900° F. Because of its ready availability, Inconel X appears to be the preferred substrate alloy for applications in which the temperature does not exceed 1500° F. René 41 would have to be used in applications involving higher temperatures. Improved coating life was derived by either preoxidizing the substrate metals prior to the coating application or by applying a very thin (<0.0002 in.) burnished and sintered overlay to the surface of the coating. Preoxidation did not affect the friction coefficient. The overlay generally resulted i n a higher friction coefficient than that obtained without the overlay. The combination of both modifications resulted in longer coating life and in friction coefficients intermediate between those obtained with either modification alone.

Nickel alloy

Host galaxies of ultra-strong Mg ii absorbers at z ∼ 0.7

ABSTRACT We report spectroscopic identification of the host galaxies of 18 ultra-strong Mg ii systems (USMg ii) at 0.6 ≤ z ≤ 0.8. We created the largest sample by merging these with 20 host galaxies from our previous survey within 0.4 ≤ z ≤ 0.6. Using this sample, we confirm that the measured impact parameters ($\rm 6.3\leqslant D[kpc] \leqslant 120$ with a median of 19 kpc) are much larger than expected, and the USMg ii host galaxies do not follow the canonical $\rm {\it W}_{2796}-{\it D}$ anticorrelation. We show that the presence and significance of this anticorrelation may depend on the sample selection. The $\rm {\it W}_{2796}-{\it D}$ anticorrelation seen for the general Mg ii absorbers show a mild evolution at low $\rm W_{2796}$ end over the redshift range 0.4 ≤ z ≤ 1.5 with an increase of the impact parameters. Compared to the host galaxies of normal Mg ii absorbers, USMg ii host galaxies are brighter and more massive for a given impact parameter. While the USMg ii systems preferentially pick star-forming galaxies, they exhibit slightly lower ongoing star-forming rates compared to main sequence galaxies with the same stellar mass, suggesting a transition from star-forming to quiescent states. For a limiting magnitude of mr < 23.6, at least 29 per cent of the USMg ii host galaxies are isolated, and the width of the Mg ii absorption in these cases may originate from gas flows (infall/outflow) in isolated haloes of massive star forming but not starbursting galaxies. We associate more than one galaxy with the absorber in $\ge 21~{{\ \rm per\ cent}}$ cases, where interactions may cause wide velocity spread.

Astronomy & Astrophysics

Self-generated electrokinetic flows from active-charged boundary patterns

We develop a hydrodynamic description of self-generated electrolyte flow in capillaries whose bounding walls feature nonuniform distributions of charge nonuniform active ionic fluxes. The hydrodynamic velocity arising in such a system has components that are forbidden by symmetry in the absence of charge and fluxes. However, when these two boundary mechanisms are simultaneously present, they can lead to a symmetry broken state where steady flows with both unidirectional and circulatory components emerge. We show that these flow states arise when modulated boundary patterns of charge and fluxes are offset by a flux-charge phase difference, which is associated with the separation between sites of their peak densities on the wall. Mismatch in diffusivity of cationic and anionic species can modify the flow states and becomes an enhancing factor when fluxes of both ion species are being produced together at the same site. We demonstrate that this mechanism can be realized with a microfluidic generator that is powered by enzyme-coated patches that catalyze reactants in the solution to produce fluxes of ions. The local ionic elevation or depletion, which disrupts a nonuniform double layer, promotes self-induced gradients yielding persistent body forces to generate bulk fluid motion. Our work quantifies a boundary-driven mechanism behind self-sustained electrolyte flow in confined environments that exists without any external bulk-imposed fields or gradients. It provides a theoretical framework for understanding the combined effect of active and charged boundaries that are relevant in biological or soft matter systems, and can be utilized in electrofluidic and iontronic applications.

active matter

Resistance of Silicon Nitride Turbine Components to Erosion and Hot Corrosion/Oxidation Attack

Silicon nitride turbine components are under intensive development by AlliedSignal to enable a new generation of higher power density auxiliary power systems. In order to be viable in the intended applications, silicon nitride turbine airfoils must be designed for survival in aggressive oxidizing combustion gas environments. Erosive and corrosive damage to ceramic airfoils from ingested sand and sea salt must be avoided. Recent engine test experience demonstrated that NT154 silicon nitride turbine vanes have exceptional resistance to sand erosion, relative to superalloys used in production engines. Similarly, NT154 silicon nitride has excellent resistance to oxidation in the temperature range of interest - up to 1400 C. Hot corrosion attack of superalloy gas turbine components is well documented. While hot corrosion from ingested sea salt will attack silicon nitride substantially less than the superalloys being replaced in initial engine applications, this degradation has the potential to limit component lives in advanced engine applications. Hot corrosion adversely affects the strength of silicon nitride in the 850 to 1300 C range. Since unacceptable reductions in strength must be rapidly identified and avoided, AlliedSignal and the NASA Lewis Research Center have pioneered the development of an environmental life prediction model for silicon nitride turbine components. Strength retention in flexure specimens following 1 to 3300 hour exposures to high temperature oxidation and hot corrosion has been measured and used to calibrate the life prediction model. Predicted component life is dependent upon engine design (stress, temperature, pressure, fuel/air ratio, gas velocity, and inlet air filtration), mission usage (fuel sulfur content, location (salt in air), and times at duty cycle power points), and material parameters. Preliminary analyses indicate that the hot corrosion resistance of NT154 silicon nitride is adequate for AlliedSignal's initial engine applications. Protective coatings and/or inlet air filtration may be required to achieve required ceramic component lives in more aggressive environments.

Nitride

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)

Petermann Glacier on the Brink: Progress, Challenges and Insights

Petermann Glacier, the largest marine-terminating glacier in northern Greenland based on catchment area and ice discharge, plays a key role in regulating ice discharge from the Greenland Ice Sheet into the Arctic Ocean. With an upstream catchment connected to the ice sheet interior via a deep subglacial canyon, its future stability has major implications for sea level rise. In this review, we synthesize recent advances in understanding Petermann’s dynamics across three critical interfaces: the ice–ocean, ice–atmosphere, and ice–bed boundaries. At the surface, observations show that reanalysis products underestimate air temperatures and melt, while regional climate models diverge significantly in their estimates of surface mass balance, underscoring the need for improved in situ data and models. At the ocean boundary, enhanced basal melting driven by both subglacial runoff and Atlantic water intrusions is identified as the dominant driver of recent mass loss of Petermann Glacier, with continued warming posing a serious threat to the stability of the floating tongue. At the bed, new geophysical synthesis reveals complex geology, likely spatial variability in geothermal heat flux, and the influence of the megacanyon on seasonal hydrology and velocity fluctuations. Petermann’s mass balance has been negative in recent decades without corresponding flow acceleration. However, the glacier has undergone significant calving events, and the current rifting that began in September 2025 highlights its vulnerability. This upcoming calving event underscores the timeliness of this review, as it will put Petermann Glacier’s terminus at its most retreated position since records began in 1923. The anticipated retreat of the ice tongue also reduces buttressing and brings the terminus closer to the grounding zone, and modeling studies suggest that calving within 12 km of the grounding zone could potentially trigger dynamic retreat, accelerating ice discharge, and a doubling of flow speeds. We conclude that improved observations, sustained monitoring of oceanographic and atmospheric properties, and high-resolution modeling are critical to constraining projections of Petermann Glacier’s future and its role in the stability of the Greenland Ice Sheet.

Dominik Fahrner

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu

Carbon Tetrachloride Degradation Results for 200-ZP-1 Operable Unit

Carbon tetrachloride (CT) contamination in the 200-ZP-1 Operable Unit (OU) at the Hanford Site originated from large-volume discharges to the subsurface during plutonium production operations between 1955 and 1973. Contamination migrated through more than 70 meters of unsaturated sediment to reach the underlying unconfined aquifer, where it persists as a large and complex groundwater plume. The 200-ZP-1 OU Record of Decision (ROD) requires that groundwater CT concentrations be reduced to 3.4 µg/L within 125 years. Current groundwater modeling projections estimate that the existing pump-and-treat, even when combined with monitored natural attenuation (specifically hydrolysis), will not achieve this target within the designated timeframe. A fundamental contributor to this shortfall is the extremely slow rate of CT hydrolysis under Hanford aquifer conditions, which has been estimated to have a half-life of 630 years. If faster-acting biotic and abiotic degradation processes are operating within the aquifer, their contribution to CT mass reduction could have a meaningful impact. However, site-specific measurements of these processes and their rates have not previously been performed. This report documents the results of a two-phase laboratory investigation designed to characterize and quantify the capacity of site-specific 200-ZP-1 OU sediments and groundwater to support natural attenuation of CT through biotic and abiotic pathways. In this context, degradation capacity is defined as the intrinsic potential of the subsurface matrix to transform CT under optimized, controlled conditions. System capacity is evaluated in two ways: (1) as rate-limited capacity, which establishes the maximum kinetic velocity of CT transformation and is measured using half-lives and first order rate constants; and (2) as mass limited capacity, which defines the total contaminant mass the batch experimental system can degrade before reactants are exhausted, representing the maximum amount of contaminant the microbial community and reactive mineral phases can transform under the experimental conditions.

abiotic degradation

Robustness of pairwise kinematic Sunyaev–Zel’dovich effect to optical-cluster-selection bias

The pairwise kinematic Sunyaev–Zel’dovich (kSZ) effect measures both the pairwise motion between galaxy groups and clusters and the amount of gas within them, providing a tracer for cosmic growth. To interpret the cosmological information in the kSZ measurements, it is crucial to understand the optical-cluster-selection bias on the kSZ observables. Line-of-sight structures that contribute to both the optical observable (e.g., richness) and the cosmological signal can induce a correlation between these two quantities at a fixed cluster mass. The selection bias arising from this correlation is a key systematic effect for cosmological analyses. For cosmological observables such as cluster abundance and weak lensing, controlling this selection bias may help explain the tension between the DES-Y1 results and the Planck constraints. In order to test for a kSZ effect equivalent of such a bias, we adopted an alternative mock richness based on galaxy counts within cylindrical volumes along the line of sight. We applied the cylindrical count method to hydrodynamical simulations across a wide range of galaxy-selection criteria, assigning richness consistent with DES-Y1 to the mock clusters. When comparing optically selected clusters to mass-selected halos, we find no significant bias on pairwise kSZ signals, pairwise velocities, or optical depth within our uncertainty limits of approximately 16, 10, and 8%, respectively.

79 ASTRONOMY AND ASTROPHYSICS

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation