Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization methods”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Evaluating methods to reduce duration of near-threshold fatigue crack growth rate measurements for low-alloy steels in hydrogen gas

Measurement of the near-threshold fatigue crack growth rate (da/dN) vs. stress-intensity factor range (∆K) relationship in hydrogen gas is essential for maximizing the calculated design fatigue life of high-pressure hydrogen storage vessels. However, such measurements are rarely performed, since the low cyclic loading frequencies applied in standard practice lead to prohibitively protracted test durations. The objective of this study was to demonstrate two means for reducing test durations when measuring near-threshold da/dN vs. ΔK relationships under decreasing ΔK for low-alloy pressure vessel steels in hydrogen gas: 1) imposing steeper K-gradients relative to the recommended limits in standards such as ASTM E647, and 2) increasing cyclic loading frequency relative to typical values applied during fatigue crack growth testing of low-alloy steels in hydrogen gas. Recognizing that steeper K-gradients could amplify loading-history effects, test methods employing this approach were designed to mitigate such effects by either maintaining constant K max or gradually increasing the K-gradient as the threshold was approached. Although the varying K-gradient method was vulnerable to loading-history effects in the form of plasticity-induced crack closure, particularly at lower stress ratio (R) and higher starting K max values, these effects could be compensated by applying the adjusted compliance ratio (ACR) method. Here, it was demonstrated that steeper K-gradients in concert with increased cyclic loading frequency reduced the duration of near-threshold fatigue crack growth tests in hydrogen gas by more than 99% relative to standard practices.

Fatigue threshold↗
DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Reliability, stability during long-term storage, and intra-individual variation of circulating levels of osteopontin, osteoprotegerin, vascular endothelial growth factor-A, and interleukin-17A

Background: Studies in many populations have reported associations between circulating cytokine levels and various physiological or pathological conditions. However, the reliability of cytokine measurements in population studies, which measure cytokines in multiple assays over a prolonged period, has not been adequately examined; nor has stability during sample storage or intra-individual variation been assessed. Methods: We assessed (1) analytical reliability in short- and long-term repeated measurements; (2) stability and analytical reliability during long-term sample storage, and (3) variability within individuals over seasons, of four cytokines—osteopontin (OPN), osteoprotegerin (OPG), vascular endothelial growth factor-A (VEGF-A), and interleukin-17A (IL-17A). Measurements in plasma or serum samples were made with commercial kits according to standard procedures. Estimation was performed by fitting a random or mixed effects linear model on the log scale. Results: In repeated assays over a short period, OPN, OPG, and VEGF-A had acceptable reliability, with intra- and inter-assay coefficients of variation (CV) less than 0.11. Reliability of IL-17A was poor, with inter- and intra-assay CV 0.85 and 0.43, respectively. During long-term storage, OPG significantly decayed (-33% per year; 95% confidence interval [-54, -3.7]), but not OPN or VEGF-A (-0.3% or -6.3% per year, respectively). Intra- and inter-assay CV over a long period were comparable to that in a short period except for a slight increase in inter-assay CV of VEGF-A. Within-individual variation was small for OPN and VEGF-A, with intra-class correlations (ICC) 0.68 and 0.83, respectively, but large for OPG (ICC 0.11). Conclusions: We conclude that OPN and VEGF-A can be reliably measured in a large population, that IL-17A is suitable only for small experiments, and that OPG should be assessed with caution due to degradation during storage and intra-individual variation. The overall results of our study illustrate the need for validation under relevant conditions when measuring circulating cytokines in population studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Subsurface mechanical damage of fused silica glass during grinding by various sub-aperture tools with and without ultrasonics

The subsurface mechanical damage (SSD) depth after grinding fused silica glass with a comprehensive set of sub-aperture fixed abrasive grinding tools [cup, wheel, belt, pad, and rotary face mill (with and without ultrasonics)] and process parameters has been statistically measured using the taper wedge technique and evaluated. Consistent with a previously reported grinding model [J. Non-Cryst. Solids 352, 5601–5617 (2006) Crossref , Materials Science and Technology of Optical Fabrication (Wiley & Sons, 2018)], based on the sliding indentation fracture by sliding particles or asperities where the normal load per particle determines the depth of the fracture (and ultimately the overall SSD depth), the dominant factor controlling SSD depth was found to be the abrasive size regardless of the tool type and process conditions. Compared to full aperture grinding methods, the overall SSD depth was higher using the sub-aperture tools, likely due to the higher effective pressure and higher load per particle distribution. Here, in addition to abrasive size, a significant reduction in SSD depth was achieved by: (1) reducing the load distribution on the abrasive particles via increase in contact area and/or decrease in mechanical loading; (2) using a more compliant host tool medium; and (3) in what we believe is a more novel way, using ultrasonics. Combining low abrasive size, larger contact area, and a compliant host, the 6 µm diamond in a resin matrix (Trizact) on a foam pad led to very low SSD depth (~ 4.6 µm), relatively fast grinding rate (186 mm 3 /h), and little or no figure degradation. This grinding tool/process is an attractive choice for final grind, resulting in significantly reduced polish out (i.e., “grey out”) times. With the rotary face mill tool, the use of ultrasonics consistently led to a SSD depth reduction (ranging from 17%–34%). A new fracture mechanics-based model, to the best of our knowledge, where the relevant normal load is parallel to the feed direction, has been developed to explain how ultrasonics leads to lower SSD depth. The key factors, supported by finite element stress analysis and load measurements, are (1) the initiation of fractures at higher z heights during the tool’s ultrasonic vertical oscillations, thus propagating less deep into workpiece; (2) reduction in load (and therefore reduction in fracture propagation distance) due to smaller tool-workpiece feed direction contact area (again caused by higher heights relative to depth of cut); (3) upward movement of the tool during oscillation leads to fracturing toward the surface instead into the depth; and finally (4) at tool’s lowest point of oscillation cycle, there may not be enough time for the fracture to propagate to its full length.

Optics and optical instruments↗
DOE OSTI2025

Development of a TSR-based method for understanding structural relationships of cofactors and local environments in photosystem I

All chemical forms of energy and oxygen on Earth are generated via photosynthesis where light energy is converted into redox energy by two photosystems (PS I and PS II). There is an increasing number of PS I 3D structures deposited in the Protein Data Bank (PDB). The Triangular Spatial Relationship (TSR)-based algorithm converts 3D structures into integers (TSR keys). A comprehensive study was conducted, by taking advantage of the PS I 3D structures and the TSR-based algorithm, to answer three questions: (i) Are electron cofactors including P700, A -1 and A 0 , which are chemically identical chlorophylls, structurally different? (ii) There are two electron transfer chains (A and B branches) in PS I. Are the cofactors on both branches structurally different? (iii) Are the amino acids in cofactor binding sites structurally different from those not in cofactor binding sites? The key contributions and important findings include: (i) a novel TSR-based method for representing 3D structures of pigments as well as for quantifying pigment structures was developed; (ii) the results revealed that the redox cofactor, P700, are structurally conserved and different from other redox factors. Similar situations were also observed for both A -1 and A 0 ; (iii) the results demonstrated structural differences between A and B branches for the redox cofactors P700, A -1 , A 0 and A 1 as well as their cofactor binding sites; (iv) the tryptophan residues close to A 0 and A 1 are structurally conserved; (v) The TSR-based method outperforms the Root Mean Square Deviation (RMSD) and the Ultrafast Shape Recognition (USR) methods. The structural analyses of redox cofactors and their binding sites provide a foundation for understanding the unique chemical and physical properties of each redox cofactor in PS I, which are essential for modulating the rate and direction of energy and electron transfers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗
DOE OSTI2022

The Strange Case of Ground-Coupled Airwaves on Seismoacoustic Stations at Local to Near-Regional Scales

Here we investigate the application of ground-coupled airwaves observed by seismoacoustic stations at local to near-regional scales to detect signals of interest and determine back-azimuth information. Ground-coupled airwaves are created from incident pressure waves traveling through the atmosphere that couple to the earth and transmit as a seismic wave with retrograde elliptical motion. Previous studies at sub-local scales (<10 km from a source of interest) found the back-azimuth to the source could be accurately determined from seismoacoustic signals recorded by acoustic and 3-component seismic sensors spatially separated on the order of 10 to 150 m. The potential back-azimuth directions are estimated from the coherent signals between the acoustic and vertical seismic data, via a propagation-induced phase shift of the seismoacoustic signal. A unique solution is then informed by the particle motion of the 3-component seismic station, which was previously found to be less accurate than the seismoacoustic-sensor method. We investigate the applicability of this technique to greater source-receiver distances, from 50-100 km and up to 400 km, which contains pressure waves with tropospheric and stratospheric ray paths, respectively. Specifically, we analyze seismoacoustic sources with ground truth from rocket motor fuel elimination events at the Utah Test and Training Range (UTTR) as well as a 2020 rocket launch in Southern California. From these sources we observe evidence that while coherent signals can be seen from both sources on multiple seismoacoustic station pairs, the determined ground-coupled airwave back-azimuths are more complicated than results at more local scales. Our findings suggest more complex factors including incidence angle, coupling location, subsurface material, and atmospheric propagation effects need to be fully investigated before the ground-coupled airwave back-azimuth determination method can be applied or assessed at these further distances.

58 GEOSCIENCES↗
DOE OSTI2024

SISGR: The regulation of carbon fixation in plant and green algae: Rubisco activase and the origin of heat inactivation of CO 2 assimilation

Rubisco activase (Rca) is a critical AAA+ ATPase protein complex that remodels and promotes the Rubisco enzyme, a key player in photosynthetic performance and carbon fixation. The assembly and function of the Rca protein complex are regulated by a range of factors, including subunit concentration, nucleotide-binding states, thermal conditions, metal-ion coordination, and post-translational modifications, such as phosphorylation. Despite its importance in photosynthesis, the detailed molecular mechanisms underlying the regulation of plant Rca and how it activates Rubisco remain elusive. This project aims to bridge this knowledge gap by integrating sophisticated enzymology tools with single-molecule methods and high-resolution electron microscopy to elucidate the structure and function of plant Rca. Through these multiple approaches, we have systematically investigated how the activity of plant Rca is impacted by various factors, such as phosphorylation and metal-ion coordination. The Rca complex assembly/disassembly dynamics were captured using anti-Brownian electrokinetic (ABEL) trap-based measurements, providing unprecedented insight into its structural flexibility and diverse assembly states. Furthermore, the structural analysis of Rca through electron crystallography and single-particle cryogenic electron microscopy (cryo-EM) reveals novel assembly states of the spinach Rca, providing insight into the mechanistic action for Rubisco remodeling. By combining cutting-edge tools and approaches, this work uncovers critical aspects of Rca’s regulation and assembly, paving the way for a deeper understanding of its role in photosynthetic efficiency and the potential for enhancing carbon fixation in crops.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Measurement of the cross section of $$ {e}^{+}{e}^{-}\to {\Xi}^{-}{\overline{\Xi}}^{+} $$ at center-of-mass energies between 3.510 and 4.843 GeV

Abstract Usinge + e − collision data corresponding to a total integrated luminosity of 12.9 fb −1 collected with the BESIII detector at the BEPCII collider, the exclusive Born cross sections and the effective form factors of the reaction$$ {e}^{+}{e}^{-}\to {\Xi}^{-}{\overline{\Xi}}^{+} $$ e + e − → Ξ − Ξ ¯ + are measured via the single baryon-tag method at 23 center-of-mass energies between 3.510 and 4.843 GeV. Evidence for the decay$$ \psi (3770)\to {\Xi}^{-}{\overline{\Xi}}^{+} $$ ψ 3770 → Ξ − Ξ ¯ + is observed with a significance of 4.5σby analyzing the measured cross sections together with earlier BESIII results. For the other charmonium(-like) statesψ(4040),ψ(4160),Y(4230),Y(4360),ψ(4415), andY(4660), no significant signal of their decay to$$ {\Xi}^{-}{\overline{\Xi}}^{+} $$ Ξ − Ξ ¯ + is found. For these states, upper limits of the products of the branching fraction and the electronic partial width at the 90% confidence level are provided.

Physics↗
DOE OSTI2020

Hazard analysis for identifying common cause failures of digital safety systems using a redundancy-guided systems-theoretic approach

Replacing the existing aging analog instrumentation and control (I&C) systems with modern safety control and protection, digital technology offers one of the foremost means of performance improvements and cost reductions for the existing nuclear power plants (NPPs). However, the qualification of digital I&C systems remains a challenge, especially considering the issue of software common-cause failures (CCFs), which are difficult to address. With the application and upgrades of advanced digital I&C systems, software CCFs have become a potential threat to plant safety because most redundant designs use similar digital platforms or software in the operating and application systems. With complex designs of multilayer redundancy to meet the single-failure criterion, digital I&C safety systems (e.g., engineered safety-features actuation system [ESFAS]) are of a particular concern in the U.S. Nuclear Regulatory Commission (NRC) licensing procedures. Here, this paper applies a modularized approach to conduct redundancy-guided systems-theoretic hazard analysis for an advanced digital ESFAS with multilevel redundancy designs. Systematic methods and risk-informed tools are incorporated to address both hardware and software CCFs, which provide guidance to eliminate the causal factors of potential single points of failure in the design of digital safety systems in advanced plant designs.

42 ENGINEERING↗
DOE OSTI2020

Improving air leakage prediction of buildings using the fan pressurization method with the Weighted Line of Organic Correlation

In many countries, the fan pressurization method is the most frequently chosen approach for measuring the air leakage of houses. The measurements are usually performed at pressures that far exceed pressures to which buildings are exposed to under normal conditions. A fit of these tests to the power-law formulation allows an extrapolation to data points outside the measured pressure range. With the Ordinary Least Square (OLS) fitting method, the pressure exponent and flow coefficient can be determined. However, the measurement results are highly sensitive to uncertainties induced by external factors like changing wind conditions during the tests, which is neglected by OLS. This may lead to errors in the prediction of flows at lower pressures. The Weighted Line of Organic Correlation (WLOC) is an alternative approach and takes measurement uncertainty into account. In this paper, a statistical analysis of an extensive data set of pressurization measurements has been performed. Both regression techniques have been compared for almost 7500 fan pressurization measurements of six houses in 109 different house leak configurations. The variability in predicting pressure exponent and flow coefficient for both WLOC and OLS regression was analyzed using probability density functions. It was found that the Weighted Line of Organic Correlation significantly decreases the uncertainty in predicting pressure exponent, flow coefficient, and other low-pressure air leakage metrics compared to the Ordinary Least Square fitting. Here, the authors highly recommend the implementation of WLOC in current measurement standards and test equipment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Microbial community diversity changes during voltage reversal repair in a 12-unit microbial fuel cell

Microbial fuel cell stacks (MFC-Stack) are often confronted with voltage reversals, likely due to an interplay between microbial community dynamics and insufficient electric circuit balancing. Herein, we provide new insight into voltage reversals by examining the microbiomes of twelve MFC units of a 12-liter Pilot-MFC-Stack during repair. Different biofilm repair methods (self-healing, electrostimulation, and re-acclimatization upon cross-inoculation) were used to evaluate the microbial community response. In addition, MFC-Stack simulation was performed based on Kirchhoff’s Second Law to predict values for source potentials and post-evaluate internal resistances. Analysis of the 16S rRNA amplicon sequencing data suggests that the biofilm repair methods could slowly heal damaged biofilms. Notably, severely voltage reversed MFC units had low electrogen relative abundances (18%) and positive anode potentials, while strong bioanodes and contained more than 50% electrogens and had negative anode potentials. Between-community analyses (beta diversity ordination and multinomial regression) of the voltage reversed MFC units revealed differences among biofilms in contrast to healthy/strong MFC units. Permutational multivariate analysis of variance (PERMANOVA) confirmed that reversed biofilms were, indeed, significantly (p < 0.05) different from stronger ones. Overall, these analyses demonstrated the utility of combining electrotechnical and microbial community analyses, especially beta diversity ordination and multinomial regression, to understand problematic MFC units and the potential success of a biofilm repair method. Finally, thicker biofilms were usually healthier and stronger, although thickness was no guarantee for proper structure and power function as all factors were interdependent. There was an evolutionary trend that strong anodes became stronger/healthier and others weaker. This spontaneous trend has to be considered to avoid irreversible voltage reversals and to repair electrogenic biofilms in an MFC-Stack.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗
DOE OSTI2024

Evaluation and development of flow condensation correlations using the data from low GWP refrigerants in an axial micro-fin aluminum tube

To mitigate global warming, the world is transitioning to refrigerants with low global warming potential (GWP). Supporting this shift requires a model that can accurately predict the heat transfer and pressure drop of new refrigerants, crucial for designing efficient heat exchangers. Existing models, however, are largely based on currently deployed refrigerants and primarily developed for unexpanded micro-fin tubes with spiral angles of 6° to 30°. Their applicability to new refrigerants, especially in expanded micro-fin tubes, is uncertain. This study assesses the performance of four well-known condensation models for six emerging refrigerants—R-32, R-454B, R-454C, R-455A, R-1234yf, and R-1234ze(E)—against experimental data. Initially, the Han and Lee (2005) model shows the best prediction accuracy with a mean absolute deviation (MAD) of 22.1 %. To enhance the accuracy of heat transfer models for new refrigerants and geometries with large temperature glides, two approaches are proposed. Here, the first approach applies a simple correction factor, reducing the MAD of the Cavallini et al. (2009) model from 68.2 % to 15.4 %. The second approach uses the variable metric method for minimization, fitting new constants to the data. This optimization results in the Kedzierski and Goncalves (1997) model achieving the highest accuracy, with a MAD of 13.1 %. For pressure drop models, the Cavallini et al. (1997) model is the most accurate with a MAD of 6.4 %, followed by the Haraguchi et al. (1993) model with a MAD of 9.4 %. Due to its simplicity, the Haraguchi et al. (1993) model is a practical option for predicting frictional pressure drop.

Axial micro-fin tubes↗
DOE OSTI2026

Multi-Spectroscopic Determination of Exchange Coupling, Zero-Field Splitting, and g-Matrices in Radical-Bridged Dinuclear Fe(III) Complexes

When the energy gap, Δ, between the lowest-lying spin manifolds within a spin-exchange coupled molecule approaches Δ/k B ≈ 300 K, the traditional temperature-dependence (T < 400 K) of the molar magnetic susceptibility is not always a reliable way to obtain a good estimate of intramolecular exchange couplings. We develop a spectroscopic approach capable of accurately parametrizing complex magnetic Hamiltonians by exploiting the separation of the anisotropy and exchange energy scales in strongly coupled magnetic molecules. Specifically, we combine inelastic neutron scattering, high-frequency electron paramagnetic resonance, far-infrared magneto-spectroscopy and magnetometry, and obtain detailed information about the magnetic properties of a series of diiron complexes derived from [[Fe(cth)] 2 (dxbq)] 3+ (H 2 dxbq: 2,5-dihydroxy-1,4-benzoquinone (x = h) or 3,6-dichloro-2,5-dihydroxy-1,4-benzoquinone (x = c), cth: 5,5,7,12,12,14-hexamethyl-1,4,8,11-tetraazacyclotetradecane). Well-isolated S = 9/2 ground states emerge due to strong direct antiferromagnetic exchange between the Fe 3+ centers (S = 5/2) and the radical bridging benzoquinone ligand (S = 1/2). The specific sensitivities and transition selection rules of the applied methods allow us to determine the parameters of the microscopic Hamiltonian including exchange coupling, fourth-order Stevens operators and g-factors. Our methodology is directly portable to other strongly coupled molecular compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Galaxy cluster matter profiles - I. Self-similarity, mass calibration, and observable-mass relation validation employing cluster mass posteriors

We present a study of the weak lensing inferred matter profiles ΔΣ(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel’dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles ΔΣ(R/R200c)/(ρcritR200c) with a lower dispersion than the unscaled ΔΣ(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M200|ζ̂, λ̂, z), which describes the individual cluster masses given their tSZE (ζ̂) and optical (λ̂, z) observables. This method enables simultaneous constraints on richness λ-mass and tSZE detection significance ζ-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPT×DES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPT×DES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the λ-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical ΛCDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.Key words: gravitational lensing: weak / galaxies: clusters: general / large-scale structure of Universe

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Investigation of Gas Turbine Internal Cooling Using Supercritical CO 2 —Effect of Surface Roughness and Channel Aspect Ratio

In this paper, an experimental and numerical investigation of internal cooling channels with rib turbulators is presented with sCO 2 as the working fluid at process conditions (pressure-20.7 MPa and temperature up to 150 °C). The effect of channel aspect ratio up to 2:1 on thermal-hydraulic performance is explored in additively manufactured rectangular channels and square channels, both with and without 60 deg ribs on the top and bottom sides. The Wilson-plot method is employed to experimentally measure channel-averaged Nusselt number over a Reynolds number range up to 370,000. The friction factor is calculated from pressure drop and mass flow rate and additionally, the overall thermal performance factor (TPF) is reported. A companion computational fluid dynamics (CFD) simulation is performed for the rib turbulated cooling configurations reported in the experiments using the Reynolds average Navier–Stokes-based turbulence model. The objective of the numerical study is to gain insight into the local heat transfer augmentation in the ribbed channels as a result of varying the aspect ratio, channel configuration (square versus rectangular), operating conditions (Reynolds number) and the surface roughness, an inherent outcome of the additive manufacturing process. Surface roughness is simulated using sand grain roughness height (K S ) calculated from the experimental data, and a comparison is presented with the corresponding channel configuration with varying surface roughness heights starting from smooth surfaces (K S = 0). Experimental results indicate that the heat transfer augmentation is negligible in the rectangular channels with ribs on the long side compared to the square channel. However, it is enhanced by 60% in comparison to placing ribs on the shorter side. The TPF remains constant at around 1 for the entire range of Reynolds numbers consistent with prior work at the National Energy Technology Laboratory (NETL). Here, the simulation results highlight that increased surface roughness can have a favorable considerable influence on Nusselt number and overall thermal performance enhancement.

42 ENGINEERING↗