Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic model”

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 361 records · Page 20

Towards Co-evolution of Membrane Transport and Metabolism

Protocellular boundaries were inextricably connected to the metabolism they encapsulated: to be inheritable, early metabolism must have led to an increased rate of growth and division of vesicles and, similarly, transport through vesicle boundaries must have supported the evolution of metabolism. Even though explaining how this coupling emerged and evolved in the absence of the complex machinery of modern cells is one of the key issues in studies on the origin of life, little is known about the biochemical and biophysical processes that might have been involved. This gap in our knowledge is a major impediment in efforts to construct scenarios for the origin of life and laboratory models of protocells. A combination of experimental and computational studies carried out by us and our collaborators is aimed at helping to close this gap. Properties of membranes might have contributed to the selection of RNA as an early biopolymer. A kinetic mechanism was proposed (Sacerdote & Szostak, 2005) in which ribose was supplied more quickly than other aldopentoses to primordial cells for preferential incorporation of ribonucleotides into nucleic acids. This proposal is based on a finding that ribose permeates membranes an order of magnitude faster than its diastereomers, arabinose and xylose. Our computer simulations, which yield permeation rates in excellent agreement with experiment, and kinetic modeling explain this phenomenon in terms of inter- and intramolecular interactions involving exocyclic hydroxyl groups attached to carbon atoms of the pyranose ring (Wei and Pohorille, 2009). They also constrain scenarios for the formation of the earliest nucleic acids (Wei and Pohorille, 2013). In one scenario, sugars permeate protocellular walls and subsequently are used to synthesize nucleic acids inside protocells. As long as this process proceeds at the rate faster than 6x10(exp -3)/s, ribose derivatives will be available for synthesis easier than their diastereomers. If nucleosides or their activated derivatives are synthesized outside protocells and subsequently transported across protocellular membranes the kinetic mechanism does not apply because all diastereomers, which have their sugars in the furanose rather than pyranose form, permeate the membrane at approximately the same rate. Properties of membranes might have been also coupled to metabolism involving peptides. Recently, Adamala and Szostak (2013) have shown that a dipeptide inside fatty-acid vesicles catalyzes the formation of another dipeptide that binds to vesicle walls and, by doing so, promotes their growth at the expense of other vesicles. This coupling of metabolism, permeability of vesicles and their growth is the first demonstration of evolutionary advantage imparted by small, membrane-bound peptides. Building on this work we have calculated the rate at which different blocked amino acids are delivered to a protocell for synthesis of dipeptides. We have further shown that the dipeptides are located at the water-membrane interface rather than in the center of the bilayer. On these basis it is anticipated that other dipeptides containing aromatic, but not necessarily hydrophobic amino acids (e.g. tyrosine) could have the same catalytic effects. Insight from these studies allows for estimating the rate of vesicle growth and the rates of dipeptide synthesis required to keep the system in balance. These results, in combination with our earlier studies, lead to a general scenario for evolution from membrane-bound dipeptides to ion channels in the origin of life.

Metabolism↗

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗

CFD Analysis of Emissions for a Candidate N+3 Combustor

An effort was undertaken to analyze the performance of a model Lean-Direct Injection (LDI) combustor designed to meet emissions and performance goals for NASA's N+3 program. Computational predictions of Emissions Index (EINOx) and combustor exit temperature were obtained for operation at typical power conditions expected of a small-core, high pressure-ratio (greater than 50), high T3 inlet temperature (greater than 950K) N+3 combustor. Reacting-flow computations were performed with the National Combustion Code (NCC) for a model N+3 LDI combustor, which consisted of a nine-element LDI flame-tube derived from a previous generation (N+2) thirteen-element LDI design. A consistent approach to mesh-optimization, spraymodeling and kinetics-modeling was used, in order to leverage the lessons learned from previous N+2 flame-tube analysis with the NCC. The NCC predictions for the current, non-optimized N+3 combustor operating indicated a 74% increase in NOx emissions as compared to that of the emissions-optimized, parent N+2 LDI combustor.

nitrogen oxides↗

Analysis of Nitrogen Shockwave Radiation at Velocity Around 7km/s

In the last decade, the Electric Arc Shock Tube (EAST) facility has been widely used to provide experimental data aiming to validate radiative heating models used for atmospheric entry, mostly at velocities above 9km/s in air. In more recent years, another experimental campaign aimed to look at velocity below 9 km/s and its analysis resulted in a significant improvement of the air kinetic model using new rates found in the literature but discrepancies were still observed, especially for the N2 second positive system. Despite these previous studies, the overall uncertainties related to surface heating for lunar return remain around 50% To reduce these uncertainties, the last experimental campaign at EAST, named test 62, was conducted in pure nitrogen at velocity ranging from 6 to 11 km/s at a freestream pressure of 0.2 Torr to provide more fundamental data on nitrogen. A previous study analyzed the data obtained at a velocity of 10.3 km/s, where atomic radiation is dominant, and provided spatially resolved density measurements of different states of nitrogen. This study will focus on experimental data obtained at shock velocity of 6.88 and 7.15 km/s where molecular nitrogen is expected to be the dominant source of radiation.

Augustin Tibere-Inglesse↗

Introducing the VRT gas turbine combustor

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual function use of the incoming air. The feasibility of the concept was demonstrated by water analogue tests and 3-D computer modeling. The computer model predicted a 50 percent reduction in pattern factor when compared to a state of the art conventional combustor. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. The results of the chemical kinetics model require further investigation, as the NO(x) predictions did not correlate with the available experimental and analytical data base.

Melconian, Jerry O.↗

Introducing the VRT gas turbine combustor

An innovative annular combustor configuration is being developed for aircraft and other gas turbine engines. This design has the potential of permitting higher turbine inlet temperatures by reducing the pattern factor and providing a major reduction in NO(x) emission. The design concept is based on a Variable Residence Time (VRT) technique which allows large fuel particles adequate time to completely burn in the circumferentially mixed primary zone. High durability of the combustor is achieved by dual function use of the incoming air. The feasibility of the concept was demonstrated by water analogue tests and 3-D computer modeling. The computer model predicted a 50 percent reduction in pattern factor when compared to a state of the art conventional combustor. The VRT combustor uses only half the number of fuel nozzles of the conventional configuration. The results of the chemical kinetics model require further investigation, as the NO(x) predictions did not correlate with the available experimental and analytical data base.

Melconian, Jerry O.↗

Explicit least squares system parameter identification for exact differential input/output models

The equation error for a class of systems modeled by input/output differential operator equations has the potential to be integrated exactly, given the input/output data on a finite time interval, thereby opening up the possibility of using an explicit least squares estimation technique for system parameter identification. The paper delineates the class of models for which this is possible and shows how the explicit least squares cost function can be obtained in a way that obviates dealing with unknown initial and boundary conditions. The approach is illustrated by two examples: a second order chemical kinetics model and a third order system of Lorenz equations.

Pearson, A. E.↗

Acceleration mechanisms of energetic ion debris in laser-driven tin plasma EUV sources

Laser-driven tin plasmas are driving new-generation nanolithography as sources of extreme ultraviolet (EUV) radiation centered at 13.5 nm. A major challenge facing industrial EUV source development is predicting energetic ion debris produced during the plasma expansion that may damage the sensitive EUV channeling multilayer optics. Gaining a detailed understanding of the plasma dynamics and ion acceleration mechanisms in these sources could provide critical insights for designing debris mitigation strategies in future high-power EUV sources. Here, we develop a fully kinetic model of tin-EUV sources using one-dimensional particle-in-cell simulations to study ion debris acceleration, which will be valuable for cross-validation of radiation-hydrodynamic simulations. An inverse-bremsstrahlung heating operator is used to model the interaction of a tin target with an Nd:YAG laser, and thermal conduction is included through a Monte Carlo Coulomb collision operator. While the large-scale evolution is in reasonable agreement with analogous hydrodynamic simulations, the significant timescale for collisional equilibration between electrons and ions allows for the development of prominent two-temperature features. A collimated flow of energetic ions is produced with a spectrum that is significantly enhanced at high energies compared to fluid simulations. The dominant acceleration mechanism is found to be a large-scale electric field supported mainly by the electron pressure gradient, which is enhanced in the kinetic simulations due to the increased electron temperature. We discuss the implications of these results for future modeling of tin-EUV sources and the development of debris mitigation schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Induced nucleation of carbon dust in red giant stars

This study quantitatively tests the proposed model of induced nucleation of carbonaceous grains in carbon-rich red giant stars. Induced nucleation is the process of grain growth initiated by the presence of reactive surfaces provided by seed particles. The numerical study was performed using a deailed chemical kinetic model of carbon deposition, grain coagulation, and homogeneous nucleation of polycyclic aromatic hydrocarbons (PAHs). The model uses a method of moments to keep track of developing grain population in the forming dust shell. We test the efficiency of grain formation for large ranges of dust shell parameters typical for carbon stars. Our model is capable of producing a range of optically thick and thin dust shells in carbon stars. Results are in accord with (IRAS) spectral classes of carbon stars. The resulting composite grains produced are consistent with those recently found in ancient meteorites. This model also provides a realistic explanation for high abundances of (PAHs) in the interstellar medium and some planetary nebulae.

Cadwell, Brian J.↗

Physical and numerical sources of computational inefficiency in integration of chemical kinetic rate equations: Etiology, treatment and prognosis

The design of a very fast, automatic black-box code for homogeneous, gas-phase chemical kinetics problems requires an understanding of the physical and numerical sources of computational inefficiency. Some major sources reviewed in this report are stiffness of the governing ordinary differential equations (ODE's) and its detection, choice of appropriate method (i.e., integration algorithm plus step-size control strategy), nonphysical initial conditions, and too frequent evaluation of thermochemical and kinetic properties. Specific techniques are recommended (and some advised against) for improving or overcoming the identified problem areas. It is argued that, because reactive species increase exponentially with time during induction, and all species exhibit asymptotic, exponential decay with time during equilibration, exponential-fitted integration algorithms are inherently more accurate for kinetics modeling than classical, polynomial-interpolant methods for the same computational work. But current codes using the exponential-fitted method lack the sophisticated stepsize-control logic of existing black-box ODE solver codes, such as EPISODE and LSODE. The ultimate chemical kinetics code does not exist yet, but the general characteristics of such a code are becoming apparent.

Pratt, D. T.↗

Production, Outflow Velocity, and Radial Distribution of H2O and OH in the Coma of Comet C/1995 O1 [Hale-Bopp] from Wide Field Imaging of OH

Observations of OH are a useful proxy of the water production rate (Q(sub H2O)) and outflow velocity (V(sub out)) in comets. From wide field images taken on 03/28/1997 and 04/08/1997 that capture the entire scale length of the OH coma of comet C/1995 O1 (Hale-Bopp), we obtain Q(sub H2O) from the model-independent method of aperture summation. With an adaptive ring summation algorithm, we extract the radial brightness distribution of OH 0-0 band emission out to cometocentric distances of up to 10(exp 6) km, both as azimuthal averages and in quadrants covering different position angles relative to the comet-Sun line. These profiles are fit using both fixed and variable velocity 2-component spherical expansion models to estimate V(sub OH) with increasing distance from the nucleus. The OH coma of Hale-Bopp was more spatially extended than previous comets, and this extension is best matched by a variable acceleration of H2O and OH that acted across the entire coma, but was strongest within 1-2 x 10(exp 4) km from the nucleus. Our models indicate that V(sub OH) at the edge of our detectable field of view (10(exp 6) km) was approx. 2-3 times greater in Hale-Bopp than for a 1P/Halley-class comet at 1 AU, which is consistent with the results of more sophisticated gas-kinetic models, extrapolation from previous observations of OH in comets with Q(sub H2O) greater than 10(exp 29)/s , and direct radio measurements of the outer coma Hale-Bopp OH velocity. The most probable source of this acceleration is thermalization of the excess energy of dissociation of H2O and OH over an extended collisional coma. When the coma is broken down by quadrants in position angle, we find an azimuthal asymmetry in the radial distribution that is characterized by an increase in the spatial extent of OH in the region between the orbit-trailing and anti-sunward directions. Model fits specific to this area and comparison with radio OH measurements suggest greater acceleration here, with V(sub OH) approx. 1.5 times greater at a 10(exp 6) km cometocentric distance than elsewhere in the coma. We discuss several mechanisms that may have acted within the coma to produce the observed effect.

Harris, Walter M.↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction, aircraft performance model,↗

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction↗

A Detailed Reaction Mechanism for Thiosulfate Oxidation by Ozone in Aqueous Environments

The ozone oxidation, or ozonation, of thiosulfate is an important reaction for wastewater processing, where it is used for remediation of mining effluents, and for studying aerosol chemistry, where its fast reaction rate makes it an excellent model reaction. Although thiosulfate ozonation has been studied since the 1950s, challenges remain in developing a realistic reaction mechanism that can satisfactorily account for all observed products with a sequence of elementary reaction steps. Here, we present novel measurements using trapped microdroplets to study the pH-dependent thiosulfate ozonation kinetics. We detect known products and intermediates, including SO32-, SO42-, S3O62-, and S4O62-, establishing agreement with the literature. However, we identify S2O42- as a new reaction intermediate and find that the currently accepted mechanism does not directly explain observed pH effects. Thus, we develop a new mechanism, which incorporates S2O42- as an intermediate and uses elementary steps to explain the pH dependence of thiosulfate ozonation. The proposed mechanism is tested using a kinetic model benchmarked to the experiments presented here, then compared to literature data. We demonstrate good agreement between the proposed thiosulfate ozonation mechanism and experiments, suggesting that the insights in this paper can be leveraged in wastewater treatment and in understanding potential climate impacts.

Deal, Alexandra M↗

Integrating Carbon Capture, Utilization, & Sequestration into Chemical Pulp Mills

The U.S. pulp and paper industry presents a unique and largely untapped opportunity for large- scale carbon dioxide removal (CDR). Unlike most industrial sectors, pulp mills rely heavily on biomass, meaning that much of their carbon emissions originate from atmospheric CO₂ that was recently captured by plants. If this biogenic CO₂ can be captured and permanently stored, pulp mills can be transformed from carbon emitters into net carbon removal facilities. This project was motivated by that opportunity and aimed to develop and evaluate integrated, low-cost strategies for capturing, utilizing, and sequestering CO₂ within existing chemical pulping operations. The scope of this work focused on four complementary innovations designed to integrate seamlessly into kraft pulp mill infrastructure: (1) in situ CO₂ capture within the recovery cycle, (2) oxy-fuel retrofitting of the rotary lime kiln to produce a high-purity CO₂ stream, (3) ex situ CO₂ capture and mineralization using pulp mill residues (dregs, grits, and lime mud), and (4) beneficial reuse of these residues as mineral carbonate fertilizers. The project combined process modeling, laboratory experimentation, life cycle assessment (LCA), and field trials to evaluate the technical feasibility, economic viability, and environmental impact of these approaches. The results demonstrate that pulp mills can serve as effective platforms for carbon removal when equipped with integrated carbon capture systems. Process modeling showed that combining sodium spiking with oxy-fuel calcination significantly enhances CO₂ capture efficiency while reducing costs by up to 31% compared to conventional configurations. Experimental work further revealed that calcination behavior in high-CO₂ environments differs substantially from traditional systems, leading to the development of a new kinetic model that predicts reaction rates under these conditions. This model provides essential design guidance for next-generation decarbonized lime kilns. In parallel, the project demonstrated that alkaline mineral residues generated during pulping operations can be repurposed as a sustainable alternative to agricultural lime. Across a wide range of soils in the southeastern United States, these materials performed equivalently to commercial lime in adjusting soil pH while offering lower greenhouse gas emissions and reduced cost. Field and greenhouse studies confirmed that crop and tree growth responses were comparable, supporting their viability as a drop-in replacement. This co-product pathway provides a practical utilization strategy that offsets costs and improves overall system economics. A major contribution of this project is the first comprehensive life cycle assessment of carbon removal in pulp and paper systems across multiple system boundaries. Results show that retrofitted mills can achieve carbon removal efficiencies ranging from 12% to 92%, depending on how the system is defined. This finding highlights a critical issue in carbon accounting: reported performance is highly sensitive to methodological choices. By explicitly quantifying these differences, this work provides valuable guidance for policymakers, carbon registries, and project developers working to standardize carbon removal metrics. From a commercialization perspective, the technologies investigated in this project are well- aligned with existing industrial infrastructure, minimizing the need for entirely new facilities. 3 DE-EE0009413 Industry engagement throughout the project—including collaboration with pulp and paper companies, equipment manufacturers, and carbon removal developers—has accelerated the transition from research to deployment. Notably, a commercial developer is actively pursuing carbon capture projects at pulp mills in the southeastern United States and has cited this research as a contributing foundation. The emergence of voluntary carbon markets and long-term offtake agreements further strengthens the business case for implementation. The broader public benefits of this work are significant. By enabling large-scale carbon removal using existing industrial systems, this approach offers a near-term pathway to reduce atmospheric CO₂ concentrations while supporting domestic manufacturing and rural economies. The reuse of industrial residues as fertilizers reduces reliance on mined materials, lowers costs for farmers, and decreases environmental impacts associated with conventional lime production. In addition, the project has supported workforce development by training graduate students and researchers in carbon capture technologies, helping to build capacity in a critical area of national interest. In conclusion, this project demonstrates that integrated carbon capture, utilization, and sequestration in pulp mills is both technically feasible and economically promising. By combining process innovation, experimental validation, and systems-level analysis, the work advances the understanding of how biomass-based industries can contribute to climate mitigation. The findings provide a strong foundation for commercial deployment and offer a scalable solution for transforming a major U.S. industry into a source of durable carbon removal.

09 BIOMASS FUELS↗

Modeling of hypersonic reacting flows

The effect of nonequilibrium thermochemical processes on the aerodynamic and heat transfer characteristics of a hypersonic vehicle is illustrated with examples. It is shown that the conventional method of predicting chemical reactions always predict that the flow is closer to equilibrium than it actually is, leading to incorrect predictions of aerodynamic characteristics of a vehicle. A method is presented for predicting chemical processes using a two-temperature kinetic model, which is developed on the basis of combining conservation equations for vibrational energy and for electron electronic energy, which are also derived. The validity of the two-temperature model is proven using two types of experimental data: radiation and shock shapes.

Park, Chul↗

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗

The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure-temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

planets and satellites↗