Search NASA⌕ Search

SEARCH · Search NASA

Results for “boiling two-phase flow”

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.

627 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Entropy-defect synergy for dual luminescence mechanism in spinel: Time-resolved anti-counterfeiting and fingerprint visualization

Multimodal luminescent materials, while promising for anti-counterfeiting, often lack dynamic time-dependent responses and controllable spatial distribution, limiting their encryption capabilities in the spatiotemporal dimension. Here, this work presents a coordinated control strategy based on entropy and defect engineering, and uses a backpropagation (BP) neural network for material screening to successfully prepare spinel Mg 0.8 (Fe 0.04 Co 0.04 Ni 0.04 Cu 0.04 Zn 0.04 )Cr 2 O 4 (MgA 5 CO) phosphors with time-dependent dynamic luminescence behavior. This phosphor simultaneously activated the d-d transition luminescence (∼618 nm) derived from Co 2+ /Cr 3+ and the defect luminescence (∼398 nm) related to zinc vacancies (V Zn ) in a single-phase solid solution. The phosphor exhibits a time-dependent color evolution from pink to purple under fixed-wavelength excitation, due to the different excited-state dynamics and decay lifetimes associated with the d-d transition and defect luminescence. Structural characterization and spectral analysis confirmed the existence of V Zn and its significant role in defect luminescence process. The fluorescent and dynamic luminescent properties of entropy-based spinel oxide enable its use in advanced anti-counterfeiting applications like fingerprint recognition and color-changing dedicated anti-counterfeiting mark, showing promise in high-end and time-dynamic anti-counterfeiting fields. This research not only developed a new type of fluorescent dynamic anti-counterfeiting material, but also provided a new idea for constructing advanced optical functional materials with multiple luminescence mechanisms.

Defect project↗
DOE OSTI2024

Theoretical design and performance of three-dimensional, pillared FeS 2 cathodes

Three-dimensional (3D) electrode design can provide improved capacities and rate capabilities over conventional two-dimensional electrodes by enhancing electrical and ionic transport. Here, expanding upon our previous modeling efforts for conversion chemistry lithium-ion batteries, we develop a pseudo-four-dimensional (P4D) approach that is subsequently used to investigate the design of a pillared FeS 2 electrode. The model considers transport in three dimensions with an additional “fourth” dimension corresponding to the solid-state lithium transport within the active material particles. Additionally, we allow for expansion of the active material during the conversion reaction to understand how internal stresses impact the electrochemical performance of the cell. By optimizing the model with respect to areal capacity, we are able to predict areal capacities up to 16.8 mAh/cm 2 for an areal current density of 1.78 mA/cm 2 and a 103% improvement for the three-dimensional electrodes over planar electrodes of equal volume. Despite the promising results, our simulations suggest that 3D design may be difficult for conversion cathode materials due to the large internal stresses that arise during conversion. Nevertheless, the model is robust and adaptable to other materials that may be more suitable for 3D electrodes due to a lesser change in volume during discharge.

Conversion cathode materials↗
DOE OSTI2022

U.S. Pacific Coast Workshop Report on Preconstruction Research Recommendations (U.S. Offshore Wind Synthesis of Environmental Effects Research (SEER) Project)

In May 2022, the U.S. Offshore Wind Synthesis of Environmental Effects Research (SEER) project team hosted a stakeholder workshop focused on preconstruction (baseline) research needs for potential floating offshore wind (OSW) energy development on the U.S. Pacific Coast, including California, Oregon, and Washington. Prior to the workshop, the SEER team developed a set of initial synthesized research recommendations that were identified based on a review of relevant, publicly available resources and with advisory group input. The workshop covered three marine life breakout groups on subsequent days to discuss research recommendations related to 1) marine mammals and sea turtles, 2) fish and invertebrates, and 3) birds and bats. As part of the workshop, over a hundred participants from the public and private sectors provided feedback on various aspects of the initial research recommendations, including associated data and knowledge gaps, benefits/limitations of available methods and technologies, and technological advancements or infrastructure needed to address the recommendation. Approximately 1,000 total comments were received on the workshop MURAL boards and were synthesized in this report. Based on workshop feedback, SEER developed a final database of over 500 specific research recommendations based on more than 40 resources. In Fall 2022, the full database and a tool with updated synthesized research recommendations were disseminated on Tethys (https://tethys.pnnl.gov) to assist with informing future funding opportunities and research programming. There is a continued need to improve awareness of the potential environmental effects, monitoring technologies, and management strategies for floating OSW energy development on the U.S. Pacific Coast. Coordination of these activities will require the sustained involvement of multiple stakeholders from across sectors. Beyond the baseline considerations discussed in this workshop, future state-of-the-science activities should be planned to consider research needs across wind energy life cycle phases for all relevant wildlife taxa and associated habitat and ecosystem processes.

17 WIND ENERGY↗
DOE OSTI2026

Novel additive manufacturing for plasma facing materials ‐ creating a research pathway for minority students

This project addresses two critical and intertwined challenges in fusion energy, namely the shortage of a broadly trained scientific workforce and the lack of scalable manufacturing solutions for plasma-facing components (PFCs). Through a collaboration among Florida International University (FIU), Miami Dade College (MDC), and Purdue University, the project established structured, reproducible educational and research pathways that recruit and advance students from institutions historically outside the fusion energy enterprise, building the human capital that this field urgently needs. The project integrates the complementary research strengths of FIU and Purdue to investigate flash sintering as a transformative processing route for tungsten-based PFCs. Unlike conventional sintering approaches, flash sintering offers rapid densification at significantly reduced thermal budgets, making it a compelling candidate for fabricating complex tungsten geometries that must withstand extreme plasma-facing environments. Systematic experimental and modeling efforts will elucidate the fundamental mechanisms governing microstructure evolution, grain boundary chemistry, and thermomechanical response during flash sintering — knowledge that is presently lacking but essential for translating this technology into reliable manufacturing practice. The convergence of workforce development and cutting-edge manufacturing research positions this project to deliver measurable, durable impact: a pipeline of fusion-ready researchers cultivated through expanded institutional partnerships, and a validated materials processing framework that accelerates domestic readiness for next-generation fusion reactor construction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Feasibility of Nitrogen as a Carrier Gas for Inconel Cold Spray in Hydropower Application

Cold spray deposition of Inconel nickel alloys has emerged as a promising strategy for mitigating cavitation erosion in hydropower facilities. Most prior developments employ helium (He) as the carrier gas because it enables high particle velocities and the formation of dense coatings. However, He is nearly two orders of magnitude more expensive than nitrogen (N2), which limits its widespread adoption in the hydropower sector. Although He-based cold spray can be justified for high value repairs, the lower cost and broad availability of N2 make it an attractive alternative. This study assesses the feasibility of N2-based cold spray of Inconel 625 powders for hydropower applications. Cavitation erosion testing shows that He-based coating exhibits cavitation resistance of ~398% relative to the 304L stainless steel (SS304L) substrate. In contrast, N2-based coating shows substantially lower cavitation resistance, reaching only ~69% of the SS304L substrate, even when higher carrier gas temperature and pressure are applied. When the Inconel 625 powder size is reduced from 44 to 22 µm under N2-based cold spray conditions, cavitation resistance increases significantly to ~148% of the SS304L substrate. Additional improvement can be obtained by incorporating fine chromium carbide powders into the feedstock, resulting in cavitation resistance of ~172% of the SS304L substrate. Overall, with optimized feedstock design, N2-based cold spray offers a practical and cost-effective approach for producing cavitation resistant coatings on hydropower components, although He-based cold spray continues to deliver higher cavitation resistance.

Wang, Tianhao↗
DOE OSTI2026

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

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

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials↗
DOE OSTI2025

Robust and tunable oxide nanoscrolls for solar-driven H 2 generation and storage

Hydrogen gas is a promising alternative to fossil fuels due to its high energy output and environmentally safe byproducts. Various morphologies of photocatalytic materials have been explored for high-efficiency H 2 production, for instance, quasi-1D nanoscroll structures that provide a larger surface-to-volume ratio. Recently, we predicted layer-by-layer formation of stable oxide nanoscrolls directly from dichalcogenide precursors, eliminating the need for costly formation of two-dimensional oxides for a roll-up synthesis of nanoscrolls. Here, in this study, we evaluate the suitability of those oxide nanoscroll materials—MoO 3 , WO 3 , PdO 2 , HfO 2 , and GeO 2 —for solar-driven photocatalytic H 2 production and storage. Using excited state theory coupled with Bethe–Salpeter equation simulations, we discern their electronic and optical properties as a function of interlayer scroll spacing and find them to be highly conducive for solar-driven photocatalysis. Additionally, using ab initio molecular dynamics simulations, we show that they are also suitable for H 2 storage as the nanoscrolls exhibit an effective trapping of hydrogen, even in the presence of defects and vacancies in the oxides. This work thus demonstrates the discovery of robust and tunable oxide nanoscrolls as materials for advancing solar-driven hydrogen technologies.

Chemical compounds↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
DOE OSTI2026

Fracture Analysis of Cohesive Zone Models for Modeling Residual Stress Induced Delamination in Composite Structures

A fracture study of coupon-scale composite cylinders with embedded defects was conducted with an objective to assess and validate a modeling approach using two available cohesive material models. The study included experimental and simulation evaluations of initiation of crack growth and progression. Interrupted thermal experiments used acoustic emissions monitoring to identify the onset of crack progression during each cooling interval and ultrasonic scanning provided images of defect growth. Verification, validation, and uncertainty quantification (VVUQ) processes were performed in the assessment of the simulation predicted temperature at which crack propagation begins (quantity of interest). The Sobol sensitivity analysis identified the hoop direction elastic modulus in the carbon fiber reinforced polymer (CFRP) plies as the most influential parameter for simulations using both cohesive models, accounting for at least 70% of the variation in the temperature at crack propagation. The UQ temperature range for the Tvergaard-Hutchinson model was higher (more conservative) than the experimental acoustic measurement indicators of crack progression, while the temperature range for the Thouless-Parmigiani model enveloped the experimental data points for the primary defect size of 0.75 x 1 in. The simulations could not capture the stable crack growth indicated in the experiments. This is likely due to the models’ inability to represent anisotropic fracture toughness attributed to the structure of the orthotropic fiber weave in a woven composite laminate.

42 ENGINEERING↗
DOE OSTI2026

A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing

Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.

Kong, Xiangxu [Nanjing Univ. of Information Scienc↗
DOE OSTI2025

Topological prethermal strong zero modes on superconducting processors

Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.

Science & Technology - Other Topics↗
DOE OSTI2026

Peptoid-Based Nanosheets Exhibiting Broad Antiviral Activity Against Enveloped RNA Viruses

Enveloped RNA viruses, such as Influenza A (H1N1) and Sindbis virus, pose persistent global health threats due to their high mutation rates, efficient transmission, and frequent drug resistance. By mimicking host cell membrane receptors, multivalent virus inhibitors can block viral attachment, making them promising broad-spectrum antiviral agents. However, most of existing antivirals are often limited by strain specificity, short-lived efficacy, and toxicity. Here, we introduce a broad-spectrum antiviral platform based on highly tunable and biocompatible two-dimensional nanomembranes (2DNMs) self-assembled from amphiphilic peptoids, operating via a non-genomic, mutation-insensitive mechanism. By varying peptoid sequence, we design and synthesize over twenty different 2DNMs with various surface charge and high density of viral-attachment ligands (VALs). The self-assembled architecture of these stable 2DNMs provides cooperative noncovalent multivalent binding to virus particles that result in effective inhibition of viral infection. Screening of variants identified three leads that potently suppressed Influenza A (H1N1) and Sindbis virus infection across median tissue culture infectious dose (TCID50), plaque, RT–qPCR, and immunofluorescence assays, while maintaining >90% cell viability. These nanosheets significantly reduced infectious titers, viral RNA replication, and intracellular viral protein expression, indicating inhibition at early stages of viral entry and propagation. The sequence programmability, chemical robustness, and mutation-insensitive antiviral activity distinguish 2DNMs from traditional antivirals and positions them as a versatile materials platform for antiviral coatings, protective barriers, and prophylactic biomedical applications.

Influenza A virus↗
DOE OSTI2026

Hund's coupling governed orbital-selective superconductivity in Ba 1−𝑥 ⁢K 𝑥⁢ Fe 2 ⁢As 2

Understanding how strong electronic correlations shape superconductivity remains a central challenge in quantum materials. In multiorbital systems, correlations driven by Hund's coupling can differentiate the behavior of individual orbitals, producing the so-called Hund's metal state. How such orbital-selectivity also governs superconducting pairing, however, has remained largely unexplored experimentally. Here, in this study, we use high-resolution angle-resolved photoemission spectroscopy to systematically map the superconducting gap structure across the phase diagram of the representative iron-based superconductor Ba 1−x K x Fe 2 As 2 . We find that superconductivity evolves in a strongly orbital-dependent manner: the gap associated with the d xy orbital collapses beyond optimal doping while pairing on the d xz /d yz orbitals persists. This behavior mirrors the orbital-selective correlations observed in the normal state and reveals a direct connection between Hund's metal physics and the superconducting pairing landscape. Our results demonstrate that superconducting gaps themselves can serve as a sensitive probe of orbital-dependent correlations and suggest that Hund's coupling plays a central role in shaping pairing in multiorbital superconductors.

Corbae, Elena [SLAC National Accelerator Laborator↗
DOE OSTI2026

Topographic effects on reflected acoustic waves from the OSIRIS-REx reentry observed from stratospheric balloons

During long-distance sound propagation in planetary atmospheres, acoustic waves may reflect off the air/surface interface one or more times. For low sound frequencies and flat interfaces, the incident and reflected wave tend to be nearly identical. However, this may not be the case when the downgoing acoustic wave encounters topography. Here, we describe a set of direct and reflected acoustic signals recorded on free-flying balloons during the hypersonic entry of the OSIRIS-REx sample return capsule (SRC). In two of the three cases presented here, an impulsive reflected arrival similar in form to the direct sonic boom of the SRC is observed, followed by a diffuse coda. In contrast, one of the floating stations lacked an impulsive reflection entirely, with only the coda present. We use full waveform modeling to show how reflection in the presence of complex topography can explain coda seen in all three examples as well as the lack of impulsive arrival on the third. Our results indicate that the complex signals often observed in long range acoustic propagation could be due, in part, to interactions with topography during transmission.

Lees, Jonathan M. [University of North Carolina, C↗
DOE OSTI2026

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

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

chemiluminescence diagnostic↗
DOE OSTI2024

Fast, Nondestructive and Precise Biomass Measurements Are Possible Using Lidar-Based Convex Hull and Voxelization Algorithms

Light detection and ranging (lidar) scanning tools are available that can make rapid digital estimations of biomass. Voxelization and convex hull are two algorithms used to calculate the volume of the scanned plant canopy, which is correlated with biomass, often the primary trait of interest. Voxelization splits the scans into regular-sized cubes, or voxels, whereas the convex hull algorithm creates a polygon mesh around the outermost points of the point cloud and calculates the volume within that mesh. In this study, digital estimates of biomass were correlated against hand-harvested biomass for field-grown corn, broom corn, and energy sorghum. Voxelization (r = 0.92) and convex hull (r = 0.95) both correlated well with plant dry biomass. Lidar data were also collected in a large breeding trial with nearly 900 genotypes of energy sorghum. In contrast to the manual harvest studies, digital biomass estimations correlated poorly with yield collected from a forage harvester for both voxel count (r = 0.32) and convex hull volume (r = 0.39). However, further analysis showed that the coefficient of variation (CV, a measure of variability) for harvester-based estimates of biomass was greater than the CV of the voxel and convex-hull-based biomass estimates, indicating that poor correlation was due to harvester imprecision, not digital estimations. Overall, results indicate that the lidar-based digital biomass estimates presented here are comparable or more precise than current approaches.

Environmental Sciences & Ecology↗
DOE OSTI2024

Nature of the diffuse emission sources in the H i supershell in the galaxy IC 1613

ABSTRACT We present a study of the nearby low-metallicity dwarf galaxy IC 1613, focusing on the search for massive stars and related feedback processes, as well as for faint supernova remnants (SNR) in late stages of evolution. We obtained the deepest images of IC 1613 in the narrow-band H α, He ii and [S ii] emission lines and new long-slit spectroscopy observations using several facilities (6-m BTA, 2.5m SAI MSU, and 150RTT telescopes), in combination with the multiwavelength archival data from MUSE/VLT, VLA, XMM–Newton, and Swift/XRT. Our deep narrow-band photometry identifies several faint shells in the galaxy, and we further investigate their physical characteristics with the new long-slit spectroscopy observations and the archival multiwavelength data. Based on energy balance calculations and assumptions about their possible nature, we propose that one of the shells is a possible remnant of a supernova explosion. We study five out of eight Wolf–Rayet (WR) star candidates previously published for this galaxy using the He ii emission line mapping, MUSE/VLT archival spectra, and new long-slit spectra. Our analysis discards the considered WR candidates and finds no new ones. We found P Cyg profiles in H α line in two stars, which we classify as Luminous Blue Variable (LBV) star candidates. Overall, the galaxy IC 1613 may have a lower rate of WR star formation than previously suggested.

Astronomy & Astrophysics↗