Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time-resolved dynamic anti-counterfeiting”

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.

983 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

An Efficient Approach to Collect Inlet Characterization Data

This report considers a novel process to acquire performance data for supersonic mixed compression inlets that is faster, and therefore less costly, than the previous one. In the previous process, measurements were recorded in discrete increments. The backpressuring actuator was advanced, there was a pause to allow for airflow fluctuations to diminish, and then the measurements were recorded. While this process was methodical, reliable and accurate, it consumed significant wind-on time. In general, high-speed wind tunnels are expensive to operate considering maintenance, energy and human resource requirements. Furthermore, many facilities, such as blowdown wind tunnels, are limited in the duration that they can maintain the desired test condition. Consequently, a quick method to collect the measurements is desired. A favored approach has the actuator moving in a continuous fashion with measurements recorded as the actuator continuously progresses through the increment points. This approach, identified as the dynamic inlet characteristic data acquisition procedure, uses in situ dynamic high-speed pressure sensors to measure pressure signals. While faster, continuous movement adds dynamic effects to the data. To better understand these dynamic effects and other potential influences, a study was undertaken to compare wind tunnel measurements taken at discrete actuator position points with data gathered during continuous actuator movement. Specially, plots of inlet pressure recovery versus mass capture ratio, or total pressure characteristic curves, were examined. The data were measured during experiments in the NASA Glenn Research Center Abe Silverstein 10- by 10-Foot Supersonic Wind Tunnel (SWT) with an inlet propulsion test article. In this report, the wind tunnel experiment and the study objectives are described. The processes to reduce the data and for analysis are explained. A discussion of the results along with features noted in the data is given. Finally, conclusions from this study that may be used to guide future wind tunnel experiment planning are presented.

Inlet↗
NASA NTRSDate not supplied

An Efficient Approach to Collect Inlet Characterization Data

This report considers a novel process to acquire performance data for supersonic mixed compression inlets that is faster, and therefore less costly, than the previous one. In the previous process, measurements were recorded in discrete increments. The backpressuring actuator was advanced, there was a pause to allow for airflow fluctuations to diminish, and then the measurements were recorded. While this process was methodical, reliable and accurate, it consumed significant wind-on time. In general, high-speed wind tunnels are expensive to operate considering maintenance, energy and human resource requirements. Furthermore, many facilities, such as blowdown wind tunnels, are limited in the duration that they can maintain the desired test condition. Consequently, a quick method to collect the measurements is desired. A favored approach has the actuator moving in a continuous fashion with measurements recorded as the actuator continuously progresses through the increment points. This approach, identified as the dynamic inlet characteristic data acquisition procedure, uses in situ dynamic high-speed pressure sensors to measure pressure signals. While faster, continuous movement adds dynamic effects to the data. To better understand these dynamic effects and other potential influences, a study was undertaken to compare wind tunnel measurements taken at discrete actuator position points with data gathered during continuous actuator movement. Specially, plots of inlet pressure recovery versus mass capture ratio, or total pressure characteristic curves, were examined. The data were measured during experiments in the NASA Glenn Research Center Abe Silverstein 10- by 10-Foot Supersonic Wind Tunnel (SWT) with an inlet propulsion test article. In this report, the wind tunnel experiment and the study objectives are described. The processes to reduce the data and for analysis are explained. A discussion of the results along with features noted in the data is given. Finally, conclusions from this study that may be used to guide future wind tunnel experiment planning are presented.

Data Acquisition↗
NASA NTRS2026

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗
DOE OSTI2026

The Future Polarized Target Program at Jefferson Lab

Polarized targets have played a crucial role in Jefferson Lab's exploration of nuclear structure over the past four decades. The three original experimental halls have seen 19 separate installations of polarized solid or gas targets for use in the particle physics scattering experiments, and this trend will continue in the next decade. Five polarized target systems are in preparation for use at JLab in the coming years. Hall B will see the use of two solid polarized targets, one longitudinally polarized to the beam, the other transversely, as well as a novel 3He gas polarized target. In Hall C, new experiments will augment the tensor polarization in dynamically polarized solids. Plans are under development to bring a polarized solid target to Jefferson Lab's photon beam hall, Hall D, for the first time. To support these efforts, the JLab polarized target group is building a test laboratory to develop dynamic nuclear polarization techniques, as well as an apparatus to irradiate target material using electrons from JLab's injector test facility. In this talk, we will explore the development progress and plans for each of these efforts.

Maxwell, James [Thomas Jefferson National Accelera↗
DOE OSTI2026

Cave-based chronology of valley incision in the glaciated Northern Carpathians

While glaciers leave a mark on mountain landscapes, the importance of glacial erosion in mountains at a global scale remains insufficiently constrained despite over a century of research aimed at quantifying the rates and magnitudes of such erosion. Measuring rates of landscape change across the Quaternary in repeatedly glaciated mountain landscapes requires the ability to accurately partition the relative role of glacial versus fluvial processes in forming relief; however, the modification of landforms by surface erosion hinders quantification of process rates over the last few million years. This study presents 32 paired 26 Al and 10 Be cosmogenic nuclide burial ages from allochthonous, quartz-bearing sediment deposited within a paleophreatic karst system in the northern Carpathians. These data quantify Pliocene-to-present timing of valley incision and karst system reactivation, as well as erosion rates over time. Post-burial, muon production-correction yields 24 fully resolved ages, six minimum ages, and two unresolvable ages. Burial ages are interpreted to document the tempo of valley incision via phreatic passage formation and reveal that ∼80% of current valley relief (∼700 m) had developed by approximately 1.6 Ma. Paleoerosion rates derived from this dataset suggest little change in the processes generating sediment washed into the karst networks since the Pliocene. These results suggest that river incision dominated relief creation in Carpathian valleys that bear the imprint of glaciation. Our findings suggest the importance of investigating other mid-latitude mountainous regions worldwide to determine if they have similar erosion and topographic histories.

Carpathians↗
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗
NASA NTRS2026

Surprisingly Enriched CO and CH4 in Halley-type Comet 13P/Olbers: Clues to Its Interstellar Heritage

The native ice composition of Halley-type Comet 13P/Olbers (hereafter 13P) was measured post-perihelion, on UT 2024 August 16 using iSHELL, the 1–5 μm cross-dispersed, high-resolution facility spectrograph at the NASA InfraRed Telescope Facility. At that time, the heliocentric distance (Rh) of 13P was 1.37 au and the geocentric distance (Δ) was 1.99 au. Using a slit that delivered spectral resolving power RP λ/Δλ 4.5 × 104, two instrument settings targeted rovibrational emissions of seven parent molecules released into the coma through sublimation of ices in the nucleus. These settings provided broad and contiguous or nearly contiguous spectral coverage, one simultaneously measuring H2O, CO, and OCS, and the other simultaneously measuring CH4, C2H6, CH3OH, and H2CO. Meaningful abundances were obtained for all seven species. Our study revealed substantial enrichments of CO and CH4, the two most volatile species systematically targeted at IR wavelengths in comets. Relative to their respective mean abundance ratios among measured comets from the Oort cloud reservoir, CO in 13P was enriched by a factor of about six and CH4 was enriched by a factor of about three. Together with near-average abundances measured for CH3OH and C2H6, this suggests large endowments of CO and CH4 on interstellar grains, yet it could indicate inhibited surface chemistry on the pre-cometary grains subsequently contained in the nucleus or 13P. Assuming that this represents the primordial composition of 13P, it could reflect very low temperatures in its natal environment (∼10 K or lower).

Michael A Disanti↗
NASA NTRSDate not supplied

Enabling Interoperability in Earth System Digital Twins (ESDT): Integrating Observations, Models, and AI for Actionable Insights Through NASA'S Intelligent Systems Technology Program

NASA’s Intelligent Systems Technology Program (IST) is driving a paradigm shift in Earth science through the development of Earth System Digital Twins (ESDT). These integrated information systems create a dynamic "digital replica" of the Earth by harmonizing continuous, multi-source observations with high-fidelity models and state-of-the-art artificial intelligence (AI) that enable “What now?”, “What next?”, and “What if?” scenario building. These scenarios are reflected in NASA IST’s series of ESDTs, from the Coastal Zone Digital Twin that integrates complex data on the current state of the Chesapeake Bay to the Terrestrial Environmental Rapid-Replication and Assimilation Hydrometeorological (TerraHydro) AI-based ESDT that forecasts water movement across Earth’s surface, to the Agriculture Land Information System (AgLIS) which can be used to assess optimal planting dates and crop yield estimates. By bridging the gap between vast data archives and actionable insights, these projects enable a system-of-systems approach to understanding complex, interacting Earth processes. This poster will highlight recent innovations and future directions from NASA’s ESDT initiatives: Continuous Data Assimilation & Multi-Source Fusion. A core requirement of the ESDT work is the transition from static models to dynamic "living" replicas. This involves creating frameworks for the continual assimilation of near-real-time data from uncoordinated, heterogeneous sources, including satellite observations and airborne assets, and ground-based Internet of Things (IoT) sensors. These systems link design, operational status, and environmental data, ensuring the digital twin accurately reflects the current state of the physical Earth system. High-Fidelity Hybrid Modeling & Computational Acceleration to enable interactive "what-if" explorations, programs are moving beyond traditional, slow physical solvers by developing fast surrogate machine learning models and Deep Generative Models (DGMs). These hybrid approaches use neural networks to emulate complex physics, such as cloud feedback or ocean dynamics, at a fraction of the original computing cost, often leveraging advanced hardware like Graphics Processing Units (GPUs) to achieve the necessary scale. Federated Ecosystems & Interoperable Frameworks rather than building isolated tools, NASA IST is moving toward federated ESDTs and reusable analytic collaborative frameworks. This theme focuses on interoperability standards and common ontologies that allow specialized digital twins to interact and share data. This system-of-systems architecture supports multi-discipline investigations, such as analyzing how upstream watershed changes impact downstream urban flooding or how wildfire emissions affect regional air quality. By leveraging these advancements, ESDTs empower researchers and decision-makers to conduct real-time analysis and run complex hypothetical scenarios, ultimately improving our understanding of Earth’s evolving systems and informing critical real-world applications.

Earth System↗
NASA NTRS2005

Chapter 19A – Radar Cross Section

The material in this Section will provide the novice Flight Test Engineer (FTE) with a brief overview of the radar cross section (RCS) concept. Brief details of radar range descriptions, operations, target calibrations, and data reduction are included. It is noted that these range operations concentrate on the "static" test range rather than the much more sophisticated "dynamic" test range. The dynamic test range is what the FTE will probably work with. Radar reflectivity measurement has developed over the last two decades from a relatively simple endeavor involving the measurement of target RCS amplitude statistics to involving wide band, coherent systems that can measure high resolution images of targets as well as, in many cases, the polarization and phasing properties. This rapid growth in RCS technology has occurred because of the increased use of radar in today's commercial and military systems. In general, the goal for commercial systems is to enhance radar reflectivity, whereas the military goal is to reduce radar reflectivity. Also, classification and identification are important for military purposes because in adverse weather radar may be the only system that may be capable of separating enemy targets from friendly ones. This Section will discuss the fundamentals of RCS. Since the flight test techniques associated with RCS measurements are very similar to those of antenna pattern measurements, both will be presented at the same time in Section 19B Antenna Radiation Pattern Measurements.

Robert W Borek↗
DOE OSTI2025

Polaron catastrophe within quantum acoustics

The quantum acoustic framework has recently emerged as a nonperturbative, coherent approach to electron–lattice interactions, uncovering rich physics often obscured by perturbative methods with incoherent scattering events. Here, we model the strongly coupled dynamics of electrons and acoustic lattice vibrations within this framework, representing lattice vibrations as coherent states and electrons as quantum wave packets, in a manner distinctively different from tight-binding or discrete hopping-based approaches. We derive and numerically implement electron backaction on the lattice, providing both visual and quantitative insights into electron wave packet evolution and the formation of acoustic polarons. We investigate polaron binding energies across varying material parameters and compute key observables—including mean square displacement, kinetic energy, potential energy, and vibrational energy—over time. Our findings reveal the conditions that favor polaron formation, which is enhanced by low temperatures, high deformation potential constants, slow sound velocities, and high effective masses. Additionally, we explore the impact of external electric and magnetic fields, showing that while polaron formation remains robust under moderate fields, it is weakly suppressed at higher field strengths. These results deepen our understanding of polaron dynamics and pave the way for future studies into nontrivial transport behavior in quantum materials.

Science & Technology - Other Topics↗
NASA NTRS2026

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗
DOE OSTI2025

Distinguishing Thermal Fluctuations from Polaron Formation in Halide Perovskites

Recent angle-resolved photoelectron spectroscopy (ARPES) measurements of the hole effective mass in CsPbBr$_3$ revealed an enhancement of $\sim$50 % compared to the bare mass computed from first principles for CsPbBr$_3$ at $T = 0 K$. This large enhancement was interpreted as evidence of polaron formation. Employing accurate finite-temperature first-principles calculations, we show that the calculated hole effective mass of CsPbBr$_3$ at $T = 300 K$ can explain experimental results without invoking polarons. Thermal fluctuations are particularly strong in halide perovskites compared to conventional semiconductors such as Si and GaAs, and cannot be ignored when comparing with experiment. We not only resolve the debate on polaron formation in halide perovskites, but also demonstrate the general importance of including thermal fluctuations in first-principles calculations for strongly anharmonic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Anisotropic multi- Q order in Co x TaS 2

The cobalt-intercalated transition metal dichalcogenide Co x TaS 2 hosts a rich landscape of magnetic phases that depend sensitively on x. While the stoichiometric compound with x = 1/3 exhibits a single magnetic transition, samples with x≤0.325 display two transitions with an anomalous Hall effect (AHE) emerging in the lower temperature phase. Here, we resolve the spin structure in each phase by employing a suite of magneto-optical probes that include the discovery of anomalous magneto-birefringence: a spontaneous time-reversal sensitive rotation of the principal optic axes. A symmetry-based analysis identifies the AHE-active phase as an anisotropic (2+1)Q state, in which magnetic modulation at one wavevector (Q) differs in symmetry from that at the remaining two. The (2+1)Q state naturally exhibits scalar spin chirality as a mechanism for the AHE and expands the classification of multi-Q magnetic phases.

Kruppe, Jonathon [University of California, Berkel↗
NASA NTRS2026

CFD-Assisted Nodal Modeling of Sloshing in a Cryogenic Propellant Tank

During autogenous pressurization, tank sloshing causes a significant increase in pressurant consumption to maintain constant ullage pressure during draining of the tank. This increase in pressurant consumption is caused by a significant increase in condensation at the liquid-vapor interface. Sloshing strongly affects the liquid side heat transfer coefficient and thereby the condensation rate. Traditionally, sloshing is modeled by CFD code using the VOF (Volume of Fluid) method to track the liquid vapor interface during sloshing. CFD calculations require a very fine grid to accurately compute the heat and mass transfer at the interface. Therefore, computations are time consuming and prohibit performing many parametric studies often needed during the design of a new system. This paper describes an alternative approach by developing a CFD-assisted nodal model to predict system parameters more economically with reasonable accuracy. In this approach, a multi-node model of tank pressurization was developed using GFSSP. A multi-node model was needed to account for stratification. The model computes heat and mass transfer at the interface to calculate the condensation rate. The liquid side heat transfer is computed using the parameters of sloshing dynamics such as frequency, wave amplitude, and interface area. The parameters of sloshing dynamics are computed by the CFD code LOCI-Stream. The model predictions were compared with test data for several cases.

Cryogenic Tank Sloshing↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2024

Making Plasticized Polymer Electrolytes Stable Against Sodium Metal for High‐Energy Solid‐State Sodium Batteries

Solid polymer electrolytes based on plastic crystals are promising for solid-state sodium metal (Na 0 ) batteries, yet their practicality has been hindered by the notorious Na 0 -electrolyte interface instability issue, the underlying cause of which remains poorly understood. Here, in this study, by leveraging a model plasticized polymer electrolyte based on conventional succinonitrile plastic crystals, we uncover its failure origin in Na 0 batteries is associated with the formation of a thick and non-uniform solid electrolyte interphase (SEI) and whiskery Na 0 nucleation/growth. Furthermore, we design a new additive-embedded plasticized polymer electrolyte to manipulate the Na 0 deposition and SEI formulation. For the first time, we demonstrate that introducing fluoroethylene carbonate (FEC) additive into the succinonitrile-plasticized polymer electrolyte can effectively protect Na 0 against interfacial corrosion by facilitating the growth of dome-like Na 0 with thin, amorphous, and fluorine-rich SEIs, thus enabling significantly improved performances of Na//Na symmetric cells (1,800 h at 0.5 mA cm −2 ) and Na//Na 3 V 2 (PO 4 ) 3 full cells (93.0 % capacity retention after 1,200 cycles at 1 C rate in coin cells and 93.1 % capacity retention after 250 cycles at C/3 in pouch cells at room temperature). Our work provides valuable insights into the interfacial failure of plasticized polymer electrolytes and offers a promising solution to resolving the interfacial instability issue.

25 ENERGY STORAGE↗
DOE OSTI2026

Resonant scattering at the center of the galaxy cluster PKS 0745-191 with XRISM

We report evidence of the resonant scattering effect at the center of the galaxy cluster PKS 0745-191 with XRISM. We analyzed XRISM/Resolve commissioning-phase observations of the distant cluster PKS 0745-191 ($z = 0.103$) with a $54$ ks exposure. The gain drift was corrected using the onboard modulated X-ray source (MXS), and spectra were extracted from all the pixels that were well illuminated by MXS, the core region (four central pixels, $\sim$100 kpc), and the surrounding region. A single-temperature collisional ionization equilibrium (CIE) model fits the full field-of-view spectrum with $kT \approx 6$ keV and a turbulent velocity of ${\approx}120$ km s$^{-1}$. From the core ($r < 50$ kpc) spectrum, we detected a ${\approx}22\%$ suppression of the Fe xxv He$\alpha$ resonance (w) line relative to the CIE prediction. We performed a Monte Carlo simulation to calculate the resonant scattering (RS) effect using radial profiles from Chandra data. The RS-inferred turbulence agrees with that by Resolve line-broadening, demonstrating that RS provides an independent and consistent constraint on ICM turbulence. These results highlight XRISM/Resolve’s potential for turbulence studies in galaxy clusters.

Astronomy and AstroPhysics↗