SEARCH · Search NASA

Results for “Dynamical evolution and stability”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Different interfacial reactivity of lithium metal chloride electrolytes with high voltage cathodes determines solid-state battery performance

A deep understanding of the interaction of the surface of cathode materials with solid electrolytes is crucial to design advanced solid-state batteries (SSBs). This is especially true for the new class of lithium metal chloride (Li-M-Cl) solid electrolytes which are receiving rapidly growing attention due to their very high oxidative stability (>4 V) in combination with good ionic conductivity that can enable long cell cycle life. While Li-M-Cl electrolytes typically contain resource-limited metals (M) such as indium or rare earths, work has focused on substituting M with more abundant elements such as zirconium. Via operando resistance measurements using intermittent current interruption we explore the dynamic evolution of the interphase at the surface of Ni-rich NCM85 or NCM111 cathode particles inside a working SSB with three different Li-(M 1 ,M 2 )-Cl catholytes (Li 3 InCl 6 , Li 2 Sc 1/3 In 1/3 Cl 4 and Li 5/2 Y 1/2 Zr 1/2 Cl 6 ) to reveal the impact of the cationic metal substitution on the interfacial chemistry. We show that the metal plays a critical role in determining high voltage stability, contrary to prior assumptions. Using a combination of cyclic voltammetry and ultraviolet photoelectron spectroscopy measurements of the electronic band structure to assess oxidative stability; coupled with DFT calculations and ToF-SIMS to evaluate products formed at the interface at different upper cutoff potentials and degrees of delithiation, we are able to differentiate between electrochemical and chemical degradation. We find that Li 2 Sc 1/3 In 1/3 Cl 4 yields the highest (and Li 3 InCl 6 the lowest) stability against electrochemical oxidation, while Li 5/2 Y 1/2 Zr 1/2 Cl 6 undergoes a detrimental chemical reaction with oxygen released from Ni-rich NCM85 at high potentials, resulting in fast capacity fading. Altogether, our work establishes a platform for the metrics and an approach that can be utilized to efficiently evaluate the stability of new halide SEs in SSB cells.

25 ENERGY STORAGE↗

DOE OSTI2025

Uncovering the True Active Sites in Ni–N–C Catalysts for CO 2 Electroreduction

Understanding and designing active sites in single-atom catalysts (SACs) requires going beyond static models to capture their dynamic evolution under realistic electrochemical conditions. Here, in this work, we develop an integrated theoretical framework that accounts for operational conditions, by combining grand canonical density functional theory (GC-DFT) with machine-learning-accelerated sampling, to uncover structure–activity–stability relationships in Ni–N–C SACs for the CO 2 reduction reaction (CO 2 RR). A library of NiN x C 4–x (x = 0–4) motifs─representing coordination defects likely formed during high-temperature synthesis─was systematically evaluated. Under working conditions, these sites were found to undergo hydrogenation, and NiN 3 C 1_ H 1 was identified as the most probable active site. At reducing potentials, hydrogen adsorbs spontaneously at C–Ni bridge sites rather than Ni top sites, while subsurface hydrogen facilitates bent CO 2 adsorption crucial for activation. High CO 2 RR selectivity toward CO arises from site separation: Ni centers drive CO2RR, while the hydrogen evolution reaction (HER) occurs at the C–Ni bridge or N sites and from thermodynamic suppression of HER at moderate hydrogen coverage. At more negative potentials, a shift in the CO 2 RR rate-determining process (RDP) and Ni out-of-surface displacement induced by coadsorption of H and H 2 O jointly reduce activity and selectivity. Thus, both the high CO2RR selectivity of Ni–N–C catalysts and its reversal with more negative potentials can be rationalized by accounting for hydrogenated surfaces. This highlights the necessity of modeling realistic; in situ conditions. This framework provides generalizable insights into the dynamic behavior of active sites in SACs, offering guidance for the rational design of active and robust catalysts for a wide range of electrochemical reactions.

25 ENERGY STORAGE↗

DOE OSTI2024

Stabilized Oxygen Vacancy Chemistry toward High-Performance Layered Oxide Cathodes for Sodium-Ion Batteries

Anionic redox has emerged as a transformative paradigm for high-energy layered transition-metal (TM) oxide cathodes, but it is usually accompanied by the formation of anionic redox-mediated oxygen vacancies (OVs) due to irreversible oxygen release. Additionally, external factor-induced OVs (defined as intrinsic OVs) also play a pivotal role in the physicochemical properties of layered TM oxides. However, an in-depth understanding of the interplay between intrinsic and anionic redox-mediated OVs and the corresponding regulation mechanism of the dynamic evolution of OVs is still missing. Herein, we disclose the strong interrelationship between these OVs and demonstrate that the presence of intrinsic OVs in the TMO2 layers could induce weak integrity of the TM-O frameworks and unlock additional diffusion paths to trigger the generation and migration of anionic redox-mediated OVs. Accordingly, an OV stabilization strategy is proposed by deliberately introducing high-valence Nb5+, which could serve as an important building block in anchoring the oxygen sublattice and preventing the formation of a percolating OV migration network, thereby suppressing the formation/diffusion of anionic redox-mediated OVs. Consequently, superb structural integrity and improved electrochemical performance with reversible anionic redox chemistry are achieved. This work advances our understanding of the role of OVs for developing high-performance energy storage systems utilizing anionic redox.

anionic redox↗

DOE OSTI2026

Interfacial Void Formation and Self-Healing in Oxide Scales on Al-containing High-Entropy Alloy

The exceptional high-temperature oxidation resistance of Al-containing high-entropy alloys (HEAs) is often attributed to the formation of a protective α-Al 2 O 3 scale. However, the dynamic, atomic-scale mechanisms governing the stability of this scale—including interfacial void formation and the often-postulated but rarely visualized “self-healing” capacity—remain poorly understood. Herein, we reveal the complex evolution of the triple-layer oxide scale on an Al 10 CoCrFeNi HEA through combined electron microscopy and diffraction study. We show that interfacial voids are an inherent consequence of the scaling process, originating from two distinct mechanisms: the Kirkendall effect at the interface between the γ-Al 2 O 3 /α-Al 2 O 3 and alloy driven by cationic diffusion imbalance and volumetric contraction due to phase transformations at the spinel/Cr 2 O 3 interface. Crucially, we provide microstructural evidence consistent with an intrinsic self-healing response. This process is driven by coupled inward diffusion of oxygen and outward diffusion of metal cations, leading to the in-situ formation of transient θ-Al 2 O 3 and spinel phases that partially fill and seal the voids. Here, these results provide atomic-scale insights into the phase evolution, defect formation, and self-repair of oxide scales in HEAs—highlighting pathways to enhance their oxidation resistance in extreme environments.

High-entropy alloy↗

DOE OSTI2024

Characterization and Rationalization of Microstructural Evolution in GRCop-84 Processed by Laser-Powder Bed Fusion (L-PBF)

In this study, prismatic geometries of GRCop-84 [Cu-8Cr-4Nb (at. pct)] were built with laser-powder bed fusion (L-PBF) process. The samples were sectioned parallel or perpendicular to the build direction and characterized in the as-built and after post-processing with a hot-isostatically pressing (HIP) treatment. The microstructure and phase evolutions were evaluated with optical microscopy, scanning electron microscopy (SEM), electron backscattered diffraction (EBSD), and high-temperature X-ray diffraction (HTXRD) up to 1223 K. The samples in the as-built conditions exhibited predominantly columnar epitaxial and misoriented Cu-FCC grains. The microstructure evolutions are discussed based on locations within the overall build geometry, the dynamics of small melt pool shape and sectioning effects. The above grain structure did not change significantly during post-process HIP treatment. The stability of this FCC grain structure is attributed to the formation of primary stable Cr 2 Nb (Laves phase) during L-PBF, even before the emergence of FCC grains from liquid. The stability of Cr 2 Nb in both as-built and HIPed samples were evaluated using high-temperature X-ray diffraction measurements and compared with that of gas-atomized powder. The significance of these results is discussed with reference to aerospace applications.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Refining Planetary Boundary Layer Height Retrievals From Micropulse‐Lidar at Multiple ARM Sites Around the World

Abstract Knowledge of the planetary boundary layer height (PBLH) is crucial for various applications in atmospheric and environmental sciences. Lidar measurements are frequently used to monitor the evolution of the PBLH, providing more frequent observations than traditional radiosonde‐based methods. However, lidar‐derived PBLH estimates have substantial uncertainties, contingent upon the retrieval algorithm used. In addressing this, we applied the Different Thermo‐Dynamic Stabilities (DTDS) algorithm to establish a PBLH data set at five separate Department of Energy's Atmospheric Radiation Measurement sites across the globe. Both the PBLH methodology and the products are subject to rigorous assessments in terms of their uncertainties and constraints, juxtaposing them with other products. The DTDS‐derived product consistently aligns with radiosonde PBLH estimates, with correlation coefficients exceeding 0.77 across all sites. This study delves into a detailed examination of the strengths and limitations of PBLH data sets with respect to both radiosonde‐derived and other lidar‐based estimates of the PBLH by exploring their respective errors and uncertainties. It is found that varying techniques and definitions can lead to diverse PBLH retrievals due to the inherent intricacy and variability of the boundary layer. Our DTDS‐derived PBLH data set outperforms existing products derived from ceilometer data, offering a more precise representation of the PBLH. This extensive data set paves the way for advanced studies and an improved understanding of boundary‐layer dynamics, with valuable applications in weather forecasting, climate modeling, and environmental studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Oxygen vacancy dynamics in monoclinic metallic VO 2 domain structures

It was demonstrated recently that the nano-optical and nanoelectronic properties of VO 2 can be spatially programmed through the local injection of oxygen vacancies by atomic force microscope writing. In this work, we study the dynamic evolution of the patterned domain structures as a function of the oxygen vacancy concentration and the time. A threshold doping level is identified that is critical for both the metal–insulator transition and the defect stabilization. The diffusion of oxygen vacancies in the monoclinic phase is also characterized, which is directly responsible for the short lifetimes of sub-100 nm domain structures. This information is imperative for the development of oxide nanoelectronics through defect manipulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1988

Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode

A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).

Steinolfson, R. S.↗

NASA NTRS1971

Dynamics of plane stellar systems.

The evolution of initially balanced rotating disks of stars is investigated with a computer model for isolated disks of stars. An isolated, initially cold balanced disk is found to be violently unstable. Adding a sufficient amount of velocity dispersion will stabilize all small-scale disturbances. However, the disks are still unstable against slowly growing long wave-length modes and after about two rotations most disks tend to assume a bar-shaped structure. It is found that the final mass distribution over most of the disk can be closely approximated by an exponential variation, irrespective of the initial mass distribution. The gravitational two-stream instability is investigated by means of a modified computer model for infinite doubly periodic stellar systems.

Hohl, F.↗

NASA NTRS1972

Dynamics of plane stellar systems.

Hohl, F.↗

DOE OSTI2021

A Comparison of Solid Electrolyte Interphase Formation and Evolution on Highly Oriented Pyrolytic and Disordered Graphite Negative Electrodes in Lithium–Ion Batteries

The presence and stability of solid electrolyte interphase (SEI) on graphitic electrodes is vital to the performance of lithium-ion batteries (LIBs). However, the formation and evolution of SEI remain the least understood area in LIBs due to its dynamic nature, complexity in chemical composition, heterogeneity in morphology, as well as lack of reliable in situ/operando techniques for accurate characterization. In addition, chemical composition and morphology of SEI are not only affected by the choice of electrolyte, but also by the nature of the electrode surface. While introduction of defects into graphitic electrodes has promoted their electrochemical properties, how such structural defects influence SEI formation and evolution remains an open question. In this work, utilizing nondestructive operando electrochemical atomic force microscopy (EChem-AFM) the dynamic SEI formation and evolution on a pair of representative graphitic materials with and without defects, namely, highly oriented pyrolytic and disordered graphite electrodes, are systematically monitored and compared. Complementary to the characterization of SEI topographical and mechanical changes during electrochemical cycling by EChem-AFM, chemical analysis and theoretical calculations are conducted to provide mechanistic insights underlying SEI formation and evolution. The results provide guidance to engineer functional SEIs through design of carbon materials with defects for LIBs and beyond.

25 ENERGY STORAGE↗

DOE OSTI2024

Genomic factors shaping codon usage across the Saccharomycotina subphylum

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. It has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. We characterized codon usage across 1,154 strains from 1,051 species from the fungal subphylum Saccharomycotina to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns. We found a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is distinct between the 12 orders to such a degree that yeasts can be classified with an accuracy >90% using a machine learning algorithm. We also characterized the degree to which codon usage bias is impacted by translational selection. We found it was influenced by a combination of features, including the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that avoidance of the CGN codons is associated with a decline in arginine tRNA function. Consistent with previous findings, codon usage bias within the Saccharomycotina is shaped by genomic features and GC bias. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗

NASA NTRS2005

Numerical Analysis of Flow Evolution in a Helium Jet Injected into Ambient Air

A computational model to study the stability characteristics of an evolving buoyant helium gas jet in ambient air environment is presented. Numerical formulation incorporates a segregated approach to solve for the transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum using a staggered grid method. The operating parameters correspond to the Reynolds number varying from 30 to 300 to demarcate the flow dynamics in oscillating and non-oscillating regimes. Computed velocity and concentration fields were used to analyze the flow structure in the evolving jet. For Re=300 case, results showed that an instability mode that sets in during the evolution process in Earth gravity is absent in zero gravity, signifying the importance of buoyancy. Though buoyancy initiates the instability, below a certain jet exit velocity, diffusion dominates the entrainment process to make the jet non-oscillatory as observed for the Re=30 case. Initiation of the instability was found to be dependent on the interaction of buoyancy and momentum forces along the jet shear layer.

Satti, Rajani P.↗

DOE OSTI2026

Influence of the drive design on ablation front hydrodynamic instability growth in a capsule implosion at the National Ignition Facility

A critical aspect of inertial confinement fusion experiments lies in the control of ablation front instabilities during the implosion process. The growth of Rayleigh–Taylor and Richtmyer–Meshkov instabilities, seeded by target defects, can significantly degrade the performance of the implosion. Here, this study explores the influence of drive design on ablation front stability using the hydro-growth radiography platform at the National Ignition Facility. This platform allows the measurement of the ablation front hydrodynamic instability growth. Because these measurements are restricted to early convergence, their primary role is to constrain the initial instability growth and benchmark the simulations used to infer shell integrity and performance at peak velocity and ignition. Three ignition designs, Hybrid-E (HyE), SQ-n (“S” for scaling and “Q” for quality), and High temperature, High thickness (HiT) were analyzed. The results demonstrated a high dependence of the growth factor on the choice of design, capsule scale as well as the hohlraum conditions, mainly the radiative temperature and the gold M-band emission. HyE and SQ-n display a similar growth factor range, but their evolution is different due to the design differences in hohlraum conditions. HiT has the lowest growth factor of the studied designs, which is a result of a higher radiative temperature. These results highlight the importance of design choices in controlling instability dynamics. These insights inform future drive design strategies to enhance stability and efficiency in fusion ignition experiments.

Physics↗

DOE OSTI2021

Dynamically Stable Active Sites from Surface Evolution of Perovskite Materials during the Oxygen Evolution Reaction

Perovskite oxides are an important class of oxygen evolution reaction (OER) catalysts in alkaline media, despite the elusive nature of their active sites. In this work, we demonstrate that the origin of the OER activity in a La 1-x SrxCoO 3 model perovskite arises from a thin surface layer of Co hydr(oxy)oxide (CoO x H y ) that interacts with trace-level Fe species present in the electrolyte, creating dynamically stable active sites. Generation of the hydr(oxy)oxide layer is a consequence of a surface evolution process driven by the A-site dissolution and O-vacancy creation. In turn, this imparts a 10-fold improvement in stability against Co dissolution and a 3-fold increase in the activity-stability factor for CoO x H y / LSCO when compared to nanoscale Co-hydr(oxy)oxides clusters. Our results suggest new design rules for active and stable perovskite oxide-based OER materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Delving into the depths of NGC 3783 with XRISM V. Broad-band modelling of ionized outflows

The Seyfert 1 galaxy NGC 3783 hosts a multi-phase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000−2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We performed joint spectral modelling of the XMM-Newton /RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derived the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigated their thermal stability and absorption measure distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compared these results with the 2000−2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log ξ = 1.08 − 3.38 and outflow velocities of 480−1230 km s −1 . The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000−2001, but the earlier data contained more low-ionization, high-velocity components. The total column density of all the outflows in 2024 is 1.5 times larger than in 2000−2001, which means that it has been replenished by fresh material. The dominant unresolved transition array (UTA) absorber (component B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

79 ASTRONOMY AND ASTROPHYSICS↗