Search NASA⌕ Search

SEARCH · Search NASA

Results for “time evolution”

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 397 records · Page 22

An Investigation of Surge in a High-Speed Centrifugal Compressor Using Digital PIV

Compressor stall is a catastrophic breakdown of the flow in a compressor, which can lead to a loss of engine power, large pressure transients in the inlet/nacelle and engine flameout. The implementation of active or passive strategies for controlling rotating stall and surge can significantly extend the stable operating range of a compressor without substantially sacrificing performance. It is crucial to identify the dynamic changes occurring in the flow field prior to rotating stall and surge in order to successfully control these events. Generally, pressure transducer measurements are made to capture the transient response of a compressor prior to rotating stall. In this investigation, Digital Particle Imaging Velocimetry (DPIV) is used in conjunction with dynamic pressure transducers to simultaneously capture transient velocity and pressure measurements in the non-stationary flow field during compressor surge. DPIV is an instantaneous, planar measurement technique which is ideally suited for studying transient flow phenomena in high speed turbomachinery and has been used previously to successfully map the stable operating point flow field in the diffuser of a high speed centrifugal compressor. Through the acquisition of both DPIV images and transient pressure data, the time evolution of the unsteady flow during surge is revealed.

Wernet, Mark P.↗

The Delta x B = 0 Constraint Versus Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. It has previously been shown by Powell et al. (1995) that an 'almost' equivalent MHD system in non-conservative form can be derived. This non-conservative system has a better conditioned eigensystem. Aside from Powell et al., the MHD equations can be derived from basic principles in either conservative or non-conservative form. The Delta x B = 0 constraint of the MHD equations is only an initial condition constraint, it is very different from the incompressible Navier-Stokes equations in which the divergence condition is needed to close the system (i.e., to have the same number of equations and the same number of unknown). In the MHD formulations, if Delta x B = 0 initially, all one needs is to construct appropriate numerical schemes that preserve this constraint at later time evolutions. In other words, one does not need the Delta x B condition to close the MHD system. We formulate our new scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the divergence of the magnetic condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗

Microphysics, Radiation and Surface Processes in the Goddard Cumulus Ensemble (GCE) Model

One of the most promising methods to test the representation of cloud processes used in climate models is to use observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and realistic representations of cloud microphysical processes, and they can reasonably well resolve the time evolution, structure, and life cycles of clouds and cloud systems (size about 2-200 km). The CRMs also allow explicit interaction between out-going longwave (cooling) and in-coming solar (heating) radiation with clouds. Observations can provide the initial conditions and validation for CRM results. The Goddard Cumulus Ensemble (GCE) Model, a CRM, has been developed and improved at NASA/Goddard Space Flight Center over the past two decades. The GCE model has been used to understand the following: 1) water and energy cycles and their roles in the tropical climate system; 2) the vertical redistribution of ozone and trace constituents by individual clouds and well organized convective systems over various spatial scales; 3) the relationship between the vertical distribution of latent heating (phase change of water) and the large-scale (pre-storm) environment; 4) the validity of assumptions used in the representation of cloud processes in climate and global circulation models; and 5) the representation of cloud microphysical processes and their interaction with radiative forcing over tropical and midlatitude regions. Four-dimensional cloud and latent heating fields simulated from the GCE model have been provided to the TRMM Science Data and Information System (TSDIS) to develop and improve algorithms for retrieving rainfall and latent heating rates for TRMM and the NASA Earth Observing System (EOS). More than 90 referred papers using the GCE model have been published in the last two decades. Also, more than 10 national and international universities are currently using the GCE model for research and teaching. In this talk, five specific major GCE improvements: (1) ice microphysics, (2) longwave and shortwave radiative transfer processes, (3) land surface processes, (4) ocean surface fluxes and (5) ocean mixed layer processes are presented. The performance of these new GCE improvements will be examined. Observations are used for model validation.

Tao, Wei-Kuo↗

Particle Acceleration in Relativistic Jets due to Weibel Instability

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The jitter radiation (Medvedev 2000) from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

Particle Acceleration and Emission in Relativistic Jets

Shock wave acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The 'jitter' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

Particle Acceleration and Radiation associated with Magnetic Field Generation from Relativistic Collisionless Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.↗

Goddard Cumulus Ensemble (GCE) Model: Application for Understanding Precipitation Processes

One of the most promising methods to test the representation of cloud processes used in climate models is to use observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and realistic representations of cloud microphysical processes, and they can reasonably well resolve the time evolution, structure, and life cycles of clouds and cloud systems (size about 2-200 km). The CRMs also allow explicit interaction between out-going longwave (cooling) and incoming solar (heating) radiation with clouds. Observations can provide the initial conditions and validation for CRM results. The Goddard Cumulus Ensemble (GCE) Model, a cloud-resolving model, has been developed and improved at NASA/Goddard Space Flight Center over the past two decades. Dr. Joanne Simpson played a central role in GCE modeling developments and applications. She was the lead author or co-author on more than forty GCE modeling papers. In this paper, a brief discussion and review of the application of the GCE model to (1) cloud interactions and mergers, (2) convective and stratiform interaction, (3) mechanisms of cloud-radiation interaction, (4) latent heating profiles and TRMM, and (5) responses of cloud systems to large-scale processes are provided. Comparisons between the GCE model's results, other cloud-resolving model results and observations are also examined.

Tao, Wei-Kuo↗

External and Internal Reconnection in Two Filament-Carrying Magnetic-Cavity Solar Eruptions

We observe two near-limb solar filament eruptions, one of 2000 February 26 and the other of 2002 January 4. For both we use 195 Angstroms Fe XII images from the Extreme-Ultraviolet (EUV) Imaging Telescope (EIT) and magnetograms from the Michelson Doppler Imager (MDI), both on the Solar and Heliospheric Observatory (SOHO) satellite. For the earlier event we also use soft X-ray telescope (SXT), hard X-ray telescope (HXT), and Bragg Crystal Spectrometer (BCS) data from the Yohkoh satellite, and hard X-ray data from the BATSE experiment on the Compton Gamma Ray Observatory (CGRO). Both events occur in quadrupolar magnetic regions, and both have coronal features that we infer belong to the same magnetic-cavity structures as the filaments. In both cases the cavity and filament first rise slowly at approximately 10 kilometers per second prior to eruption, and then accelerate to approximately 100 kilometers per second during the eruption, although the slow-rise movement for the higher-altitude cavity elements is clearer in the later event. We estimate both filaments and both cavities to contain masses of approximately 10(exp 14-15) g and approximately 10(exp 15-16) g respectively. We consider whether two specific magnetic-reconnection-based models for eruption onset, the tether cutting and the breakout models, are consistent with our observations. In the earlier event soft X-rays from SXT show an intensity increase during the 12-minute interval over whch fast eruption begins, which is consistent with tether-cutting-model predictions. Substantial hard X-rays, however, do not occur until after fast eruption is underway, and so this is a constraint the tether-cutting model must satisfy. During the same 12-minute interval over which fast eruption begins, there are brightenings and topological changes in the corona indicative of high-altitude reconnection early in the eruption, and this is consistent with breakout predictions. In both eruptions, however, the onset of the filament fast-rise phase occurs while overlying cavity-related coronal loops are still evolving from a "closed" state to an "open" state, and our observations of the time evolution of these loops compared to the eruption state in both events are constraints the breakout model must meet. Therefore our findings are consistent with runaway-tether-cutting-type reconnection and fast breakout-type reconnection both occurring early in the fast phase of the February eruption and with both types of reconnection being important in unleashing the explosion, but we are not able to say which, if either, type of reconnection actually triggered the fast phase. We have also found specific constraints that either model, or any other model, must satisfy if correct.

Sterling, Alphonse C.↗

Particle Acceleration in Relativistic Jets Due to Weibel Instability

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a three-dimensional relativistic electromagnetic particle code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The "jitter" radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

On Nonlinear Combustion Instability in Liquid Propellant Rocket Motors

All liquid propellant rocket instability calculations in current use have limited value in the predictive sense and serve mainly as a correlating framework for the available data sets. The well-known n-t model first introduced by Crocco and Cheng in 1956 is still used as the primary analytical tool of this type. A multitude of attempts to establish practical analytical methods have achieved only limited success. These methods usually produce only stability boundary maps that are of little use in making critical design decisions in new motor development programs. Recent progress in understanding the mechanisms of combustion instability in solid propellant rockets"' provides a firm foundation for a new approach to prediction, diagnosis, and correction of the closely related problems in liquid motor instability. For predictive tools to be useful in the motor design process, they must have the capability to accurately determine: 1) time evolution of the pressure oscillations and limit amplitude, 2) critical triggering pulse amplitude, and 3) unsteady heat transfer rates at injector surfaces and chamber walls. The method described in this paper relates these critical motor characteristics directly to system design parameters. Inclusion of mechanisms such as wave steepening, vorticity production and transport, and unsteady detonation wave phenomena greatly enhance the representation of key features of motor chamber oscillatory behavior. The basic theoretical model is described and preliminary computations are compared to experimental data. A plan to develop the new predictive method into a comprehensive analysis tool is also described.

Sims, J. D.↗

Visualization of Subsurface Defects in Composites using a Focal Plane Array Infrared Camera

A technique for enhanced defect visualization in composites via transient thermography is presented in this paper. The effort targets automated defect map construction for multiple defects located in the observed area. Experimental data were collected on composite panels of different thickness with square inclusions and flat bottom holes of different depth and orientation. The time evolution of the thermal response and spatial thermal profiles are analyzed. The pattern generated by carbon fibers and the vignetting effect of the focal plane array camera make defect visualization difficult. An improvement of the defect visibility is made by the pulse phase technique and the spatial background treatment. The relationship between a size of a defect and its reconstructed image is analyzed as well. The image processing technique for noise reduction is discussed.

Plotnikov, Yuri A.↗

Adsorption and solvation of ethanol at the water liquid-vapor interface: a molecular dynamics study

The free energy profiles of methanol and ethanol at the water liquid-vapor interface at 310K were calculated using molecular dynamics computer simulations. Both alcohols exhibit a pronounced free energy minimum at the interface and, therefore, have positive adsorption at this interface. The surface excess was computed from the Gibbs adsorption isotherm and was found to be in good agreement with experimental results. Neither compound exhibits a free energy barrier between the bulk and the surface adsorbed state. Scattering calculations of ethanol molecules from a gas phase thermal distribution indicate that the mass accommodation coefficient is 0.98, and the molecules become thermalized within 10 ps of striking the interface. It was determined that the formation of the solvation structure around the ethanol molecule at the interface is not the rate-determining step in its uptake into water droplets. The motion of an ethanol molecule in a water lamella was followed for 30 ns. The time evolution of the probability distribution of finding an ethanol molecule that was initially located at the interface is very well described by the diffusion equation on the free energy surface.

NASA Discipline Exobiology↗

The 2.5-12 micrometers spectrum of comet Halley from the IKS-VEGA experiment

The infrared instrument IKS flown on board the VEGA space probes was designed for the detection of emission bands of parent molecules, and for a measurement of the size and temperature of the thermal emitting nuclear region. The instrument had three channels with cooled detectors: an "imaging channel" designed to modulate the signal of the nucleus and two spectroscopic channels operating at 2.5-5 and 6-12 micrometers, respectively, equipped with circular variable filters of resolving power approximately 50. This paper presents and discusses the results from the spectral channels. On VEGA 1, usable spectra were obtained at distances D from the comet nucleus ranging from 250,000 to 40,000 km corresponding to fields of view 4000 and 700 km in diameter, respectively. The important internal background signal caused by the instrument itself, which could not be cooled, had to be eliminated. Since no sky chopping was performed, we obtain difference spectra between the current spectrum and a reference spectrum with little or no cometary signal taken at the beginning of the observing sequence (D approximately 200,000 km). Final discrimination between cometary signal and instrumental background is achieved using their different time evolution, since the instrumental background is proportional to the slow temperature drift of the instrument, and the cometary signal due to parent molecules or dust grains is expected to vary in first order as D-1. The 2.5-5 micrometers IKS spectra definitely show strong narrow signals at 2.7 and 4.25 micrometers, attributed to the nu 3 vibrational bands of H2O and CO2, respectively, and a broader signal in the region 3.2-3.5 micrometers, which may be attributed to CH-bearing molecules. All these signals present the expected D-1 intensity variation. Weaker emission features at 3.6 and 4.7 micrometers could correspond to the nu 1 and nu 5 bands of H2CO and the (1 - 0) band of CO, respectively. Molecular production rates are derived from the observed emissions, assuming that they are due to resonance fluorescence excited by the Sun's infrared radiation. For the strong bands of H2O and CO2, the rovibrational lines are optically thick, and radiative transfer is taken into account. We derive production rates, at the moment of the VEGA 1 flyby, of approximately 10(30) sec-1 for H2O, approximately 2.7 x 10(28) sec-1 for CO2, approximately 5 x 10(28) sec-1 for CO, and 4 x 10(28) sec-1 for H2CO, if attributions to CO and H2CO are correct. The production rate of carbon atoms in CH-bearing molecules is approximately 9 x 10(29) sec-1 assuming fluorescence of molecules in the gas phase, but could be much less if the 3.2-3.5 micrometers emission is attributed to C-H stretch in polycyclic aromatic hydrocarbons or small organic grains. In addition, marginal features are present at 4.85 and 4.45 micrometers, tentatively attributed to OCS and molecules with the CN group, respectively. Broad absorption at 2.8-3.0 micrometers, as well as a narrow emission at 3.15 micrometers, which follow well the D-1 intensity variation, might be due to water ice. Emission at 2.8 micrometers is also possibly present, and might be due to OH created in vibrationally excited states after water photodissociation. The 6-12 micrometers spectrum does not show any molecular emission, nor emission in the 7.5-micrometers region. The spectrum is dominated by silicate emission showing a double structure with maxima at 9.0 and 11.2 micrometers, which suggests the presence of olivine.

NASA Discipline Number 52-10↗

Response of a Panel Structure Forced by the Noise from a Nearly Sonic Jet

A model of a high subsonic jet with a nearby array of exible, aircraft-type panels is studied numerically in two dimensions. The jet is excited by a limited duration, spatially localized starter pulse in the potential core. The long time evolution of unsteady disturbances in the jet, the responses of the panels and the ensuing radiation are computed. The results show that the spectral response of both the jet and the panels is concentrated in a relatively narrow frequency band centered at a Strouhal number (based on jet exit velocity) of approximately 0.25 and associated harmonics. The loading on the panels generally increases with downstream distance. Panel radiation is weakest in upstream directions. Interior zones of silence, due to destructive interference of radiation from the panels, are observed.

Fenno, C. C., Jr.↗

The High-Excitation Planetary Nebula NCG 246: Optical and Near-Ultraviolet Observations and Two-dimensional Numerical Models

We have imaged the planetary nebula (PN) NGC 246 in the near-ultraviolet wavelengths [Ne v] 342.6 nm, the Bowen fluorescence line of 0 111 at 344.4 nm, and a nearby line-free region centered on 338.6 nm, as well as H(alpha), [O III] 500.7 nm, and [S II] 673.0 and 671.5 nm. Imaging in the 344.4 nm line is necessary to deconvolve contamination of the [Ne v] images by O III 342.9 nm. The emission from the shell and inner parts of the nebula is detected in [Ne v]. The radial profiles of the [Ne v] brightness decrease with radius from the exciting star, indicating that the bulk of the emission from this ion is due to the hard UV stellar radiation field, with a (probably) marginal contribution from collisional ionization in a shock between the PN shell and the interstellar medium (ISM). In contrast, the radial profiles of the emission in H(alpha), [0 III] 500.7 nm, and [S II] are flatter and peak at the location of the shell. The emission of [S II] probably traces the interaction of the PN with the ambient ISM. We also present two-dimensional numerical simulations for this PN-ISM interaction. The simulations consider the stellar motion with respect to the ambient ISM, with a velocity of 85 km/s , and include the time evolution of the wind parameters and UV radiation field from the progenitor star.

Szentgyorgyi, Andrew↗

Particle acceleration, magnetic field generation, and emission in relativistic pair jets

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Recent simulations show that the Weibel instability created by relativistic pair jets is responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet propagating through an ambient plasma with and without initial magnetic fields. The growth rates of the Weibel instability depends on the distribution of pair jets. The Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. This instability is also responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The jitter radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nishikawa, K.-I.↗

Particle acceleration magnetic field generation, and emission in Relativistic pair jets

Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) are responsible for particle acceleration in relativistic pair jets. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic pair jet propagating through a pair plasma. Simulations show that the Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. Simulation results show that this instability generates and amplifies highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The "jitter' I radiation from deflected electrons can have different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. The growth rate of the Weibel instability and the resulting particle acceleration depend on the magnetic field strength and orientation, and on the initial particle distribution function. In this presentation we explore some of the dependencies of the Weibel instability and resulting particle acceleration on the magnetic field strength and orientation, and the particle distribution function.

Nishikawa, K.-I.↗

Fluid-structural dynamics of ground-based and microgravity caloric tests

Microgravity caloric tests aboard the 1983 SpaceLab1 mission produced nystagmus results with an intensity comparable to those elicited during post- and pre- flight tests, thus contradicting the basic premise of Barany's convection hypothesis for caloric stimulation. In this work, we present a dynamic fluid structural analysis of the caloric stimulation of the lateral semicircular canal based on two simultaneous driving forces for the endolymphatic flow: natural convection driven by the temperature-dependent density variation in the bulk fluid and expansive convection caused by direct volumetric displacement of the endolymph during the thermal irrigation. Direct numerical simulations indicate that on earth, the natural convection mechanism is dominant. But in the microgravity environment of orbiting spacecraft, where buoyancy effects are mitigated, expansive convection becomes the sole mechanism for producing cupular displacement. A series of transient 1 g and microgravity case studies are presented to delineate the differences between the dynamics of the 1 g and microgravity endolymphatic flows. The impact of these different flow dynamics on the endolymph-cupula fluid-structural interactions is also analyzed based on the time evolutions of cupular displacement and velocity and the transcupular pressure differences.

Gravitation↗