Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional thermal state”

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 271 records · Page 15

Infinite temperature at zero energy

We construct a family of static, geometrically local Hamiltonians that inherit eigenstate properties of periodically-driven (Floquet) systems. Our construction is a variation of the Feynman-Kitaev clock -- a well-known mapping between quantum circuits and local Hamiltonians -- where the clock register is given periodic boundary conditions. Assuming the eigenstate thermalization hypothesis (ETH) holds for the input circuit, our construction yields Hamiltonians whose eigenstates have properties characteristic of infinite temperature, like volume-law entanglement entropy, across the whole spectrum -- including the ground state. We then construct a family of exactly solvable Floquet quantum circuits whose eigenstates are shown to obey the ETH at infinite temperature. Combining the two constructions yields a new family of local Hamiltonians with provably volume-law-entangled ground states, and the first such construction where the volume law holds for all contiguous subsystems.

FOS: Physical sciences↗

Further Development of the Tamped Richtmyer-Meshkov Instability Method and Application to Molybdenum Dynamic Strength Calibration and Tabulation

The high pressure and high strain rate dynamic strength of Mo is experimentally and computationally investigated in the 3–20 GPa stress, 50–600 C temperature, and 10 5 –10 6 /s strain rate regimes using a modified tamped Richtmyer-Meshkov instability (RMI) method. Modifications to the tamped RMI method include a method to determine loading states during strain, a new strength calibration function based on interface shape, and a robust uncertainty quantification method. These modifications improve fidelity of the tamped RMI method, allowing evaluation of the compensating effects of pressure hardening, strain rate hardening, strain hardening, and thermal softening. The new calibration function based on interface shape is not limited to sinusoidal corrugations and could be applied to additional interface shapes. Plate impact experiments are performed at Argonne National Laboratory’s Advanced Photon Source’s Dynamic Compression Sector operated by Washington State University (DCS), driving a planar shock front through a corrugated Mo-D 2 O or Mo-C 8 F 18 interface, forcing the corrugation to significantly deform. The extent of interfacial deformation, RMI growth, is experimentally observed using X-ray phase contrast imaging at the DCS. RMI jet lengths and jet shapes are extracted from the experimental radiographs, then used to calibrate numerical simulations performed with the Sandia National Laboratories (SNL) hydrocode CTH. Mo yield strength, Y, as a function of shock pressure, P, strain rate, $\dot{\varepsilon }$, accumulated strain, ϵ, relative volumetric compression, RD , and temperature, T , is determined for each impact experiment and presented. The calibrated Mo yield strength values range 1.2–1.8 GPa, with strength generally decreasing as the impact stress increases. This trend is likely caused by thermal softening or strain localization. The tabular yield strength versus loading condition data presented in this paper can be used to fit complex strength models.

Voorhees, T. J. [Sandia National Lab. (SNL-CA), Li↗

Complete Initial Studies in Microstructural Studies of Nuclear Graphite

Nuclear graphite is a critical material for the performance and safety of high-temperature reactors, making its study essential for the successful deployment of next-generation reactors in the United States. During the past year, Oak Ridge National Laboratory (ORNL) has contributed extensively to advancing the understanding of nuclear graphite’s properties, behavior, and applications; particular focus has been given to irradiation effects, oxidation mechanisms, thermal and mechanical performance, and long-term stability under reactor conditions. This report provides an update on the ongoing publications being prepared at ORNL. These studies are collectively intended to build a comprehensive knowledge base on nuclear graphite. The studies focus on key aspects such as microstructural characterization, oxidation behavior, advanced analytical techniques, molten salt infiltration, technique development, and billet-to-billet variability. By summarizing the progress of these publications, this report outlines ORNL’s contributions to the field, identifies knowledge gaps, and proposes future research directions. The findings presented in these studies will help enhance the understanding of nuclear graphite’s behavior, support the qualification of nuclear graphite materials, further the development of new techniques and guidelines for using microscopy techniques, and improve reactor safety and performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT↗

Structural and Chemical Evolution of Highly Fluorinated Li‐Rich Disordered Rocksalt Oxyfluorides as a Function of Temperature

Abstract Li‐rich disordered rocksalt (DRS) oxyfluorides have emerged as promising high‐energy cathode materials for lithium‐ion batteries. While a high level of fluorination in DRS materials offers performance advantages, it can only be achieved via mechanochemical synthesis, which poses challenges of reproducibility and scalability. The definition of relationships between fluorination and structural stability is required to devise alternative methods that overcome these challenges. In this study, the thermal evolution of three highly fluorinated phases, Li 2 TMO 2 F (TM = Mn, Co, and Ni), is investigated in an inert atmosphere. Diffraction and spectroscopic techniques are utilized to examine their electronic and chemical states up until conditions of decomposition. The analysis reveals that the materials phase‐separate above 400 °C, at most. It is also observed that heat‐treated DRS materials exhibit intricate changes in the local coordination of the metals, including their spin, and ordering compared to the pristine states. The changes upon annealing are accompanied by a modulation of the voltage profile, including reduced hysteresis, when used as electrodes. These results provide an in‐depth understanding of the fundamental crystal chemistry of DRS oxyfluorides in view of their promising role as the next generation of Li‐ion cathodes.

25 ENERGY STORAGE↗

Wind and structural loads data measured on parabolic trough solar collectors at an operational power plant

Abstract Wind loading is a primary contributor to structural design costs of concentrating solar-thermal power collectors, such as heliostats and parabolic troughs. These structures must resist the mechanical forces generated by turbulent wind, while the reflector surfaces must maintain optimal optical performance. Studying wind-driven loads at a full-scale, operational concentrating solar-thermal power plant provides insights into the wind impact on the solar collector field beyond the capabilities of wind tunnel tests or state-of-the-art simulations. We conducted comprehensive field measurements of the atmospheric turbulent wind conditions and the resulting structural wind loads on parabolic troughs at the Nevada Solar One plant over a two-year period. The measurement setup included meteorological masts and structural load sensors on four trough rows. Additionally, a lidar scanned the horizontal plane above the trough field. In this study, we describe the high-resolution dataset characterizing the complex flow field and resulting structural loads. This first-of-its-kind dataset will enhance the understanding of wind loading on collector structures and will help in designing the next-generation solar collectors and photovoltaic trackers.

14 SOLAR ENERGY↗

Pressurized Water Reactor Control Rod Ejection Analysis Using PARCS, RELAP5-3D, and BISON for High Burnup Fuel

This paper describes safety assessment of high burnup fuel (HBF) during loss of coolant accidents (LOCA), reactivity-initiated accidents (RIA), and anticipated operational occurrences (AOOs) for Southern Nuclear’s Vogtle Electric Generating plant (VEGP). A goal of this project is to create a multiphysics modeling approach in aspects of thermal hydraulics, neutronics, and fuel performance, as well as to use it in combination with experimental studies (Figure 1). In support of the project, this paper seeks to address fuel performance of a HBF pin during steady state conditions and during a control rod ejection (CRE) scenario using a loose coupling of PARCS, RELAP5-3D, and BISON.. The primary goal of this work was to analyze high burnup fuel behavior during a CRE scenario. The results indicate fuel failure is not expected to occur under the conditions presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Electron Inversion and Tunneling at Silicon Thermal Oxide Interfaces for Solar-Driven Molecular Catalysis to Syngas

Semiconductor photoelectrodes are regularly coupled to solid-state heterogeneous catalysts to perform solar-driven reduction of CO 2 . Less frequently, molecular catalysts are employed to better control the reactivity toward desired products, yet the development of robust semiconductor/molecule interfaces has proven challenging. Here, we demonstrate that a 2–3 nm thermal oxide layer on Si exhibits stability in aqueous solution, high photovoltage, and a photocurrent density of ∼10 mA/cm 2 for the solar-driven photoelectrochemical reduction of a homogeneous molecular catalyst, producing syngas with an ∼2:1 H 2 to CO ratio. Because of a low defect density, the oxide interface forms an electron inversion layer with metal-like electron density at cathodic potentials. This inversion layer facilitates electron transfer to redox-active molecules via tunneling even if the molecule’s reduction potential is beyond the semiconductor’s conduction band edge. Using an electrolyte solution composed of a homogeneous cobalt bis(terpyridine) catalyst in a water/organic solvent mixture, stable photoelectrochemistry was observed under 1-sun illumination, exhibiting an ∼30% Faradaic efficiency for CO that was similar to a glassy carbon electrode under comparable conditions. Furthermore, the results demonstrate that an ultrathin thermal oxide interface is a robust platform for development of aqueous-stable, molecule-driven photoelectrocatalysis.

Catalysts↗

Mechanistic transition of stainless steel 304L corrosion in MgCl 2 phase change material under thermal cycling

Magnesium chloride (MgCl 2 )-based molten salt phase change materials (PCMs) are promising candidates for thermal energy storage (TES) owing to their advantageous thermophysical properties, however, their high corrosivity in the molten state limits material compatibility. In this study, the corrosion behavior of 304?L stainless steel (SS304L) exposed to MgCl2 PCM was investigated as a function of thermal cycling (up to 9000 cycles) relevant to PCM-based TES systems. Microstructural changes and elemental redistribution at the alloy-salt interface under thermal cycling conditions were characterized using synchrotron-based X-ray fluorescence (XRF) and scanning electron microscopy coupled with energy dispersive X-ray spectroscopy (SEM-EDS). Results show that corrosion is governed by selective Cr-dissolution and elemental diffusion near the interface and grain boundaries. A clear transition in the dominant corrosion mechanism is observed with increasing thermal cycles. Surface corrosion is the dominant mechanism up to 1000 cycles, whereas intergranular corrosion becomes the primary degradation mode beyond 2000 cycles. These mechanistic changes are interpreted in terms of evolving interfacial chemistry, microstructure and the diffusion-driven redistribution of alloying elements.

36 MATERIALS SCIENCE↗

Temperature measurement of Quark-Gluon plasma at different stages

In a Quark-Gluon Plasma (QGP), the fundamental building blocks of matter, quarks and gluons, are under extreme conditions of temperature and density. A QGP could exist in the early stages of the Universe, and in various objects and events in the cosmos. The thermodynamic and hydrodynamic properties of the QGP are described by Quantum Chromodynamics (QCD) and can be studied in heavy-ion collisions. Despite being a key thermodynamic parameter, the QGP temperature is still poorly known. Thermal lepton pairs (e + e − and μ + μ − ) are ideal penetrating probes of the true temperature of the emitting source, since their invariant-mass spectra suffer neither from strong final-state interactions nor from blue-shift effects due to rapid expansion. Here we measure the QGP temperature using thermal e+e− production at the Relativistic Heavy Ion Collider (RHIC). The average temperature from the low-mass region (in-medium ρ 0 vector-meson dominant) is (2.01 ± 0.23) × 10 12 K, consistent with the chemical freeze-out temperature from statistical models and the phase transition temperature from Lattice QCD. The average temperature from the intermediate mass region (above the ρ0 mass, QGP dominant) is significantly higher at (3.25 ± 0.60) × 10 12 K. This work provides essential experimental thermodynamic measurements to map out the QCD phase diagram and understand the properties of matter under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Model predictive control of heating, ventilation, and air conditioning (HVAC) systems: A state-of-the-art review

Due to the fast advancement of communication and information technology, intelligent buildings have garnered great interest. These buildings can forecast weather, ambient temperature, and sun irradiation and can modify heating, ventilation, and air conditioning (HVAC) operations appropriately, based on current and previous data. This change is intended to reduce HVAC system energy usage while maintaining an appropriate degree of thermal comfort and indoor air quality. Since its inception, model predictive control (MPC) has been one of the prospective solutions for HVAC management systems to reduce both costs and energy usage. Additionally, MPC is becoming increasingly practical as the processing capacity of building automation systems increases and a large quantity of monitored building data becomes available. MPC also provides the potential to improve the energy efficiency of HVAC systems via its capacity to consider limitations, to predict disruptions, and to factor in multiple competing goals such as interior thermal comfort and building energy consumption. Although substantial research has been conducted on MPC in building HVAC systems, there is a shortage of critical reviews and a lack of a comprehensive framework that formulates and defines the applications. Here, this article provides a comprehensive state-of-the-art overview of MPC in HVAC systems. Detailed discussions of modeling approaches and optimization algorithms are included. Numerous design aspects such as prediction horizon, occupancy behavior, building type, and cost function, that impact MPC performance are discussed in detail. The technical characteristics, advantages, and disadvantages of various types of modeling software are discussed. The primary objective of this work is to highlight critical design characteristics for the MPC control scheme and to give improved suggestions for future research. Moreover, numerous prospective scenarios have been suggested that might provide future research direction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Enhancement of operational Flexibility of Power Plants Using IN740

Headers are crucial components within diverse industries, especially in the energy sector, as they enable the efficient transfer of fluids. The selection of materials for headers is determined by their specific applications; for example, Grade 91 and Grade 92 steels are commonly utilized. In our research, we recommend the adoption of INCONEL alloy 740 due to its exceptional robustness and heat-resistant characteristics. One critical parameter in header systems is the heat transfer coefficient, which directly influences the efficiency of heat exchange processes and can thereby impact the structural integrity of the headers. This coefficient is closely linked with factors such as the Nusselt number, which is influenced by fluid flow characteristics and thermal properties. The flow direction within headers, whether unidirectional or multidirectional—significantly affects the overall dynamics of heat transfer. This research focuses on exploring the behavior of headers, a specific pipeline system component, using ANSYS simulation software. The study aims to predict heat transfer and mechanical behavior within headers under various conditions through steady-state and transient simulations (parts 1 and 2 of the report). Key parameters such as heat transfer coefficient, velocity, and temperature are examined with the goal of optimizing header design. Part 3 of the report addresses the critical yet underexplored relationship between pressure drop and heat transfer coefficients in the transient flow regime within headers, specifically focusing on smooth horizontal circular tubes. Limited experimental work has been conducted in this area, prompting the need for comprehensive analysis. By leveraging machine learning techniques, this research aims to establish a correlation between pressure drop and heat transfer across various flow conditions, including laminar, transient, quasi-turbulent, and turbulent regimes. The data utilized for this analysis were meticulously gathered from existing literature, capturing simultaneous measurements of pressure drops and heat transfer. Part 4: In the context of power plants, header pipes are essential components that significantly influence system performance by facilitating the collection and distribution of steam. This report highlights the critical role of header pipes in enhancing reliability, efficiency, and overall power plant performance. A key aspect of this investigation is shape optimization, which aims to maximize performance while minimizing material usage. By focusing on shape optimization, this research contributes to improved efficiency and a reduced environmental footprint for power plant installations. The methodology developed in this study emphasizes optimizing header shapes to decrease reliance on expensive alloy materials and lower maintenance costs. Furthermore, a case study was conducted using a header from an operational power plant to validate the proposed optimization techniques. Finally, part 5 addresses the challenges associated with flexible operations in boiler systems, the study outlines several strategies for enhancing the durability and reliability of steam headers. Key approaches include material selection, which involves utilizing advanced materials with superior high-temperature properties and enhanced fatigue resistance to extend the lifespan of steam headers. Design modifications are also recommended, focusing on implementing changes that mitigate thermal stress and cyclic loading to reduce the likelihood of failures under varying operating conditions. Additionally, the study emphasizes the importance of regular inspection and monitoring by establishing rigorous protocols to detect early signs of damage, allowing for timely maintenance and minimizing the risks of catastrophic failures. Furthermore, operational guidelines are developed to minimize the frequency and severity of thermal transients, ensuring stable and efficient operations. Collectively, these strategies provide a robust framework for overcoming the unique challenges posed by flexible operation in boiler systems, ultimately contributing to improved reliability and performance in power plant operations. Through this comprehensive approach, the study not only enhances the understanding of pressure drop and heat transfer relationships but also promotes advancements in design and operation that will benefit the entire power generation industry.

01 COAL, LIGNITE, AND PEAT↗

Microstructural stability of irradiated yttrium hydride under thermal cycles

Hydrides that retain their hydrogen in service are critical for advanced microreactor and space nuclear systems. Hydrogen retention is affected by evolving microstructure under temperature and irradiation extremes. This study investigated microstructural features in neutron-irradiated yttrium hydride and the stability of these features under thermal cycles. Differential scanning calorimetry was used to determine hydrogen desorption and temperatures at which significant phase changes are occurring, and transmission electron microscopy was used for microstructural analysis after each thermal cycle. Cyclic heating led to crystallization and epitaxial growth of surface oxidation and dehydriding of the matrix. Other features such as the bulk matrix crystal structure, irradiation-induced cavities, and matrix precipitation remained constant after thermal cycling. Results allow us to design better hydride alloys through microstructure tailoring and potential irradiation conditioning treatments.

36 MATERIALS SCIENCE↗

Cold Weather Impacts on Electric School Bus Performance in Aurora, Colorado

This brief highlight details the key takeaways from a project that utilized NLR's Fleet Research, Energy Data, and Insights (FleetREDI) data analysis pipeline related to electric school bus (ESB) operation. ESBs using battery energy as their primary heating source have a higher energy consumption rate in cold weather, which fleet managers can account for when planning ESB purchases and making dispatching and charging decisions. Researchers found that electric school buses operate 2-5 times more efficiently than conventional buses, on average. Cold weather can double electric school bus energy demands, but strategies such as thermal pre-conditioning significantly reduce this effect. Understanding these impacts can help fleets plan charging, dispatching, and purchase decisions.

33 ADVANCED PROPULSION SYSTEMS↗

Comparative analysis of thermal management systems in electric vehicles at extreme weather conditions: Case study on Nissan Leaf 2019 Plus, Chevrolet Bolt 2020 and Tesla Model 3 2020

With the surge in electric vehicle (EV) adoption and the need for extended driving ranges, optimizing energy efficiency, particularly through thermal management, is critical, especially in extreme weather. Managing the substantial energy needed for cabin climate control and battery temperature regulation can increase energy demands by over 50 %, severely limiting range. This study conducts a comparative analysis of thermal management systems (TMS) in three popular EV vehicles, 2020 Chevrolet Bolt, 2019 Nissan Leaf Plus, and 2020 Tesla Model 3, evaluating their distinct TMS configurations and performance under varied weather conditions. Using both numerical simulations and experimental data collected on a controlled test bench at Argonne National Laboratory, we assess how TMS architecture and operational modes influence energy consumption and range. A comprehensive TMS model was developed, integrating cabin and battery thermal sub-models in the Autonomie software platform, to simulate temperature fluctuations and range impacts. Cabin climate was modeled using a mono-zonal approach, while battery cell temperature distribution was estimated through a 2D nodal structure. Each vehicle's distinct TMS setup was evaluated: the Chevrolet Bolt and Tesla Model 3 use a dual evaporator vapor compression cycle with a PTC heater for the cabin and a coolant loop for battery thermal management; the Nissan Leaf Plus employs a heat pump with a PTC heater for the cabin and air-cooling for the battery. Tests conducted at ambient temperatures of 35°C, 22°C, -7°C, and -18°C reveal significant differences in energy use and range reduction across both configurations and conditions. At 35°C, the Tesla Model 3, Chevrolet Bolt, and Nissan Leaf Plus have a range reduction of 8%, 9%, and 13%, respectively, due to air conditioning. In winter, heating technology is paramount; at -7°C, the Nissan Leaf's heat pump configuration achieves a lower range reduction (19.3%) compared to the Tesla and Chevrolet Bolt PTC heaters, which reduce range by 28.3% and 31%, respectively. Further, this study provides valuable insights for automotive engineers, EV technology researchers, and thermal management system designers aiming to enhance electric vehicle performance by understanding how different weather conditions and TMS architectures impact energy consumption and driving range.

33 ADVANCED PROPULSION SYSTEMS↗

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Phase Transitions and Thermal Equation of State of Fe‐9wt.%Si Applied to the Moon and Mercury

Abstract Accurate knowledge of the phase transitions and thermoelastic properties of candidate iron alloys, such as Fe‐Si alloys, is essential for understanding the nature and dynamics of planetary cores. The phase diagrams of some Fe‐Si alloys between 1 atm and 16 GPa have been back‐extrapolated from higher pressures, but the resulting phase diagram of Fe 83.6 Si 16.4 (9 wt.% Si) is inconsistent with temperature‐induced changes in its electrical resistivity between 6 and 8 GPa. This study reports in situ synchrotron X‐ray diffraction (XRD) measurements on pre‐melted and powder Fe 83.6 Si 16.4 samples from ambient conditions to 60 GPa and 900 K using an externally heated diamond‐anvil cell. Upon compression at 300 K, the bcc phase persisted up to ∼38 GPa. The hcp phase appeared near 8 GPa in the pre‐melted sample, and near 17 GPa in the powder sample. The appearance of the hcp phase in the pre‐melted sample reconciles the reported changes in electrical resistivity of a similar sample, thus resolving the low‐pressure region of the phase diagram. The resulting high‐temperature Birch‐Murnaghan equation of state (EoS) and thermal EoS based on the Mie‐Gruneisen‐Debye model of the bcc and hcp structures are consistent with, and complement the literature data at higher pressures. The calculated densities based on the thermal EoS of Fe‐9wt.%Si indicate that both bcc and hcp phases agree with the reported core density estimates for the Moon and Mercury.

Berrada, Meryem↗