Search NASA⌕ Search

SEARCH · Search NASA

Results for “outer approximation”

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

Fundamental Regolith Properties, Handling and Water Capture (FLEET) Update

The FLEET project, begun in FY21, was designed to research several concepts for components and subsystems needed for extracting water from planetary ice deposits. This is one of two approaches for providing water off-Earth for direct use and also for feedstock to manufacture spacecraft propellants. The other approach is to extract lattice-bound water and hydroxyl from the minerals that make up the solid grains of rock and regolith. The ultimate intent is to integrate the findings with other efforts, such as, in the creation of complete water-from-ice production systems for various down-stream uses. The project consists of several elements, now in various stages of completion. They are presented here in the approximate order they would generally be applied to a water-extraction system: Digging icy regolith, transporting it, capturing water vapor, and sealing chambers from the outer environment.

regolith excavation↗

GHRS Observations of Cool, Low-Gravity Stars: The Outer Atmosphere and Wind of the Nearby K Supergiant Lambda Velorum - 5

UV spectra of lambda Velorum taken with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope are used to probe the structure of the outer atmospheric layers and wind and to estimate the mass-loss rate from this K5 lb-II supergiant. VLA radio observations at lambda = 3.6 cm are used to obtain an independent check on the wind velocity and mass-loss rate inferred from the UV observations, Parameters of the chromospheric structure are estimated from measurements of UV line widths, positions, and fluxes and from the UV continuum flux distribution. The ratios of optically thin C II] emission lines indicate a mean chromospheric electron density of log N(sub e) approximately equal 8.9 +/- 0.2 /cc. The profiles of these lines indicate a chromospheric turbulence (v(sub 0) approximately equal 25-36 km/s), which greatly exceeds that seen in either the photosphere or wind. The centroids of optically thin emission lines of Fe II and of the emission wings of self-reversed Fe II lines indicate that they are formed in plasma approximately at rest with respect to the photosphere of the star. This suggests that the acceleration of the wind occurs above the chromospheric regions in which these emission line photons are created. The UV continuum detected by the GHRS clearly traces the mean flux-formation temperature as it increases with height in the chromosphere from a well-defined temperature minimum of 3200 K up to about 4600 K. Emission seen in lines of C III] and Si III] provides evidence of material at higher than chromospheric temperatures in the outer atmosphere of this noncoronal star. The photon-scattering wind produces self-reversals in the strong chromospheric emission lines, which allow us to probe the velocity field of the wind. The velocities to which these self-absorptions extend increase with intrinsic line strength, and thus height in the wind, and therefore directly map the wind acceleration. The width and shape of these self-absorptions reflect a wind turbulence of approximately equal 9-21 km/s. We further characterize the wind by comparing the observations with synthetic profiles generated with the Lamers et al. Sobolev with Exact Integration (SEI) radiative transfer code, assuming simple models of the outer atmospheric structure. These comparisons indicate that the wind in 1994 can be described by a model with a wind acceleration parameter beta approximately 0.9, a terminal velocity of 29-33 km/s, and a mass-loss rate approximately 3 x 10(exp -9) solar M/yr. Modeling of the 3.6 cm radio flux observed in 1997 suggests a more slowly accelerating wind (higher beta) and/or a higher mass-loss rate than inferred from the UV line profiles. These differences may be due to temporal variations in the wind or from limitations in one or both of the models. The discrepancy is currently under investigation.

Carpenter, Kenneth G.↗

Vapor condensation on liquid surface due to laminar jet-induced mixing: The effects of system parameters

The effects of system parameters on the interface condensation rate in a laminar jet induced mixing tank are numerically studied. The physical system consists of a partially filled cylindrical tank with a slightly subcooled jet discharged from the center of the tank bottom toward the liquid-vapor interface which is at a saturation temperature corresponding to the constant tank pressure. Liquid is also withdrawn from the outer part of the tank bottom to maintain the constant liquid level. The jet velocity is selected to be low enough such that the free surface is approximately flat. The effect of vapor superheat is assumed to be negligible. Therefore, the interface condensation rate can be determined from the resulting temperature field in the liquid region alone. The nondimensional form of the steady state conservation equations are solved by a finite difference method for various system parameters including liquid height to tank diameter ratio, tank to jet diameter ratio, liquid inflow to outflow area ratio, and a heat leak parameter which characterizes the uniform wall heat flux. Detailed analyses based on the numerical solutions are performed and simplified equations are suggested for the prediction of condensation rate.

Lin, Chin-Shun↗

Vapor condensation on liquid surface due to laminar jet-induced mixing - The effects of system parameters

The effects of system parameters on the interface condensation rate in a laminar jet induced mixing tank are numerically studied. The physical system consists of a partially filled cylindrical tank with a slightly subcooled jet discharged from the center of the tank bottom toward the liquid-vapor interface which is at a saturation temperature corresponding to the constant tank pressure. Liquid is also withdrawn from the outer part of the tank bottom to maintain the constant liquid level. The jet velocity is selected to be low enough such that the free surface is approximately flat. The effect of vapor superheat is assumed to be negligible. Therefore, the interface condensation rate can be determined from the resulting temperature field in the liquid region alone. The nondimensional form of the steady state conservation equations are solved by a finite difference method for various system parameters including liquid height to tank diameter ratio, tank to jet diameter ratio, liquid inflow to outflow area ratio, and a heat leak parameter which characterizes the uniform wall heat flux. Detailed analyses based on the numerical solutions are performed and simplified equations are suggested for the prediction of condensation rate.

Lin, Chin-Shun↗

Bubble Eliminator Based on Centrifugal Flow

The fluid bubble eliminator (FBE) is a device that removes gas bubbles from a flowing liquid. The FBE contains no moving parts and does not require any power input beyond that needed to pump the liquid. In the FBE, the buoyant force for separating the gas from the liquid is provided by a radial pressure gradient associated with a centrifugal flow of the liquid and any entrained bubbles. A device based on a similar principle is described in Centrifugal Adsorption Cartridge System (MSC- 22863), which appears on page 48 of this issue. The FBE was originally intended for use in filtering bubbles out of a liquid flowing relatively slowly in a bioreactor system in microgravity. Versions that operate in normal Earth gravitation at greater flow speeds may also be feasible. The FBE (see figure) is constructed as a cartridge that includes two concentric cylinders with flanges at the ends. The outer cylinder is an impermeable housing; the inner cylinder comprises a gas-permeable, liquid-impermeable membrane covering a perforated inner tube. Multiple spiral disks that collectively constitute a spiral ramp are mounted in the space between the inner and outer cylinders. The liquid enters the FBE through an end flange, flows in the annular space between the cylinders, and leaves through the opposite end flange. The spiral disks channel the liquid into a spiral flow, the circumferential component of which gives rise to the desired centrifugal effect. The resulting radial pressure gradient forces the bubbles radially inward; that is, toward the inner cylinder. At the inner cylinder, the gas-permeable, liquid-impermeable membrane allows the bubbles to enter the perforated inner tube while keeping the liquid in the space between the inner and outer cylinders. The gas thus collected can be vented via an endflange connection to the inner tube. The centripetal acceleration (and thus the radial pressure gradient) is approximately proportional to the square of the flow speed and approximately inversely proportional to an effective radius of the annular space. For a given FBE geometry, one could increase the maximum rate at which gas could be removed by increasing the rate of flow to obtain more centripetal acceleration. In experiments and calculations oriented toward the original microgravitational application, centripetal accelerations between 0.001 and 0.012 g [where g normal Earth gravitation (.9.8 m/s2)] were considered. For operation in normal Earth gravitation, it would likely be necessary to choose the FBE geometry and the rate of flow to obtain centripetal acceleration comparable to or greater than g.

Gonda, Steve R.↗

Quasi-Sun-Pointing of Spacecraft Using Radiation Pressure

A report proposes a method of utilizing solar-radiation pressure to keep the axis of rotation of a small spin-stabilized spacecraft pointed approximately (typically, within an angle of 10 deg to 20 deg) toward the Sun. Axisymmetry is not required. Simple tilted planar vanes would be attached to the outer surface of the body, so that the resulting spacecraft would vaguely resemble a rotary fan, windmill, or propeller. The vanes would be painted black for absorption of Solar radiation. A theoretical analysis based on principles of geometric optics and mechanics has shown that torques produced by Solar-radiation pressure would cause the axis of rotation to precess toward Sun-pointing. The required vane size would be a function of the angular momentum of the spacecraft and the maximum acceptable angular deviation from Sun-pointing. The analysis also shows that the torques produced by the vanes would slowly despin the spacecraft -- an effect that could be counteracted by adding specularly reflecting "spin-up" vanes.

Spilker, Thomas↗

Thermal Inspection of a Composite Fuselage Section Using a Fixed Eigenvector Principal Component Analysis Method

A composite fuselage aircraft forward section was inspected with flash thermography. The fuselage section is 24 feet long and approximately 8 feet in diameter. The structure is primarily configured with a composite sandwich structure of carbon fiber face sheets with a Nomex(Trademark) honeycomb core. The outer surface area was inspected. The thermal data consisted of 477 data sets totaling in size of over 227 Gigabytes. Principal component analysis (PCA) was used to process the data sets for substructure and defect detection. A fixed eigenvector approach using a global covariance matrix was used and compared to a varying eigenvector approach. The fixed eigenvector approach was demonstrated to be a practical analysis method for the detection and interpretation of various defects such as paint thickness variation, possible water intrusion damage, and delamination damage. In addition, inspection considerations are discussed including coordinate system layout, manipulation of the fuselage section, and the manual scanning technique used for full coverage.

Zalameda, Joseph N.↗

The radial variation of corotating energetic particle streams in the inner and outer solar system

The radial gradient of long-lived, corotating energetic particle streams was measured using observations of .9-2.2 MeV protons from Helios 1 and 2, IMP 7, Pioneer 10 and Pioneer 11. A positive gradient of approximately 350% per AU is found between .3 AU and 1 AU. Between 1 AU and some 3-5 AU, the gradient is variable with an average value of 100% per AU which is consistent with earlier statistical results. A comparison between measurements at 9 AU and approximately 4 AU shows a negative gradient which is variable from -40 to -100% per AU. Possible solar latitudinal effects on these gradient studies are also discussed. Using solar wind and magnetic field data from Helios 1 between 1 AU and .3 AU, the relation between corotating energetic particle events in the inner solar system and the interplanetary medium is examined. It is found that the energetic particles are contained inside the high speed solar wind stream in a region adjacent to the interaction region between low speed and high speed streams.

Vanhollebeke, M. A. I.↗

Distribution function of continuously created newborn and pickup ions in outer cometary exospheres

The time evolution of the distribution function of newborn ions in the solar wind is investigated using a quasi-linear-type diffusion equation. The initial distribution is taken to be a ring beam, which is approximated by delta function in pitch angle and velocity, and it is assumed that the ions are created at a constant rate with a similar distribution. A long-time asymptotic form of ion distribution is obtained, which is a mixture of newborn ions and ions generated throughout the entire process. It is shown that the time asymptotic distribution function exists even in the presence of a continuous ionization process. The stability of the long-time asymptotic distribution was examined for the case of parallel propagation, and the results show that the distribution function can be unstable to low-frequency hydromagnetic waves. The results of the analysis were found to agree with recent satellite observations.

Gaffey, J. D., Jr.↗

The nature of dark-etching rims in meteoritic taenite.

Taenite fields when etched develop a cloudy brown rim with approximate compositional limits of 25 and 40% Ni. In iron meteorites this cloudy zone is only a few microns wide, with a sharp, high-Ni edge about 1 micron from the kamacite interface and a diffuse edge several microns from the central plessite. It is always present in irons unless the meteorite has been cosmically or terrestrially reheated. X-ray and electron diffraction of grains scratched from exceptionally larger areas of cloudy taenite in the mesosiderite Estherville show that this etching zone contains a fine exsolution of kamacite. Detailed electron-probe investigations of taenite in Estherville show that there is a step in the M-shaped Ni profile at the sharp, high-Ni edge of the cloudy region, the Ni dropping suddenly from approximately 45 to 42%. It is proposed that exsolution in the cloudy region effectively froze in the Ni profile at that temperature. On subsequent cooling only the clear outer taenite continued to equilibrate with the kamacite matrix, producing the kink in the M profile. Cloudy taenite is therefore a variety of plessite differing from the usual varieties in that it forms at lower temperatures in areas much richer in Ni, and the morphology is not crystallographically oriented. Its absence can provide a sensitive indication of reheating.

Scott, E. R. D.↗

Nonadiabatic and three-dimensional effects in compressible turbulent boundary layers

A defect stream function formulation for nonadiabatic flow with small crossflow is developed. The first-integral property of this formulation provides for two removal of the streamline curvature term in the governing equation so that the form of the reduced equation for small crossflow is the same as that for two dimensional flow. The combined law of the wall and wake is used in place of the no-slip boundary condition. The tangential velocity equation for law-of-the-wall flow is shown to be the same for three-dimensions as for two when the Boussinesq approximation applies, and a closed form solution for the crossflow angle in the inner region is obtained. Analytic solutions for nonadiabatic, compressible, equilibrium flow with a Clauser outer-region eddy-viscosity model are obtained, and excellent agreement with experimental skin friction and velocity profile data for nonadiabatic, compressible flat-plate flow is achieved. An analytic solution for a linear inner-region eddy-viscosity model is also obtained; the wake function part of this solution is found to be inconsistent with the empirically established law of the wake.

Barnwell, Richard W.↗

On the theory of the turbulent boundary layer

As a rule, a division of the turbulent boundary layer is admissible: a division into a part near the wall, where the flow is governed only by the wall effects, and into an outer part, where the wall roughness and the viscosity of the flow medium affects only the wall shearing stress occurring as boundary condition but does not exert any other influence on the flow. Both parts may be investigated to a large extent independently. Under certain presuppositions there result for the outer part "similar" solutions. The theoretical considerations give a cue how to set up, by appropriate experiments and their evaluation, generally valid connections which are required for the approximate calculation of the turbulent boundary layer according to the momentum and energy theorem.

TURBULENCE↗

Determination of nitrogen to carbon abundance ratios from transition layer emission lines

We have finished studying the nitrogen to carbon abundance ratios for stars with different effective temperatures T(sub eff) and luminosities using transition layer emission lines and using spectra available in the IUE archives. The N/C abundance ratio determinations using transition layer emission lines are as accurate as the photospheric abundance determinations as found by comparison of results obtained by both methods for the same stars. Our measurements confirm photospheric abundance determinations in regions of the HR diagram where they can be obtained. Our studies have extended the temperature range to higher temperatures. They have shown the exact positions in the HR diagram where the mixing due to the outer convection zones reaches deep enough to bring nuclear processed material to the surface. This occurs at effective temperatures which are higher by delta log T(sub eff) approximately 0.04 or roughly 400 K than expected theoretically. Since the depth of the convection zone increases rapidly with decreasing T(sub eff) this may indicate considerable overshoot beyond the lower boundary of the convection zone. Our N/C abundance ratio determinations from transition layer emission lines have confirmed that the actual enrichment observed for some cool giants is larger than expected theoretically, again indicating a larger degree of mixing in several stars either from below or from above. For the supergiants it probably indicates overshoot above the convective core in the progenitor main sequence stars. For the more massive giants this may also be the case, though we did not find a correlation between delta log N/C and the absolute magnitudes, but these are rather uncertain. As byproducts of these studies we also found anomalies in Si/C and N/C abundance ratios for F giants which can be understood as the relict of surface abundance changes for their main sequence progenitors due to diffusion. This anomaly disappears for G giants, for which the depths of the convection zones are apparently deep enough to wipe out these element separations (Bohm-Vitense 1992).

Boehm-Vitense, Erika↗

Radiation-Driven Warping: Nonisothermal Disks - 2

Recent work by Pringle and by Maloney, Begelman, & Pringle has shown that geometrically thin, optically thick, accretion disks are unstable to warping driven by radiation torque from the central source. This work was confined to isothermal (i.e., surface density Sigma varies as R(sup -3/2) disks. In this paper we generalize the study of radiation-driven warping to include general power-law surface density distributions, Sigma varies as R(sup -delta).We consider the range from Delta = 3/2 (the isothermal case) to Delta = -3/2, which corresponds to a radiation-pressure-supported disk; this spans the range of surface density distributions likely to be found in real astrophysical disks. In all cases there are an infinite number of zero-crossing solutions (i.e., solutions that cross the equator), which are the physically relevant modes if the outer boundary of the disk is required to lie in a specified plane. However, unlike the isothermal disk, which is the degenerate case, the frequency eigenvalues for Delta does not equal 3/2 are all distinct. In all cases the location of the zero moves outward from the steady state (pure precession) value with increasing growth rate; thus, there is a critical minimum size for unstable disks. Modes with zeros at smaller radii are damped. The critical radius and the steady state precession rate depend only weakly on Delta. An additional analytic solution has been found for Delta = 1. The case Delta = 1 divides the solutions into two qualitatively different regimes. For Delta greater than or equal to 1, the fastest growing modes have maximum warp amplitude, close to the disk outer edge, and the ratio of Beta(sub max) to the warp amplitude at the disk inner edge, Beta(sub o), is much greater than 1. For Delta less than 1, Beta(sub max/Beta(sub o) approximately equals 1, and the warp maximum steadily approaches the origin as Delta decreases. This implies that nonlinear effects must be important if the warp extends to the disk inner edge for Delta greater than or equal to 1, but for Delta less than 1 nonlinearity will be important only if the warp amplitude is large at the origin. Because of this qualitative difference in the shapes of the warps, the effects of shadowing of the central source by the warp will also be very different in the two regimes of Delta. This has important implications for radiation-driven warping in X-ray binaries, for which the value of Delta characterizing the disk is likely to be less than unity. In real accretion disks the outer boundary condition is likely to be different from the zero-crossing condition that we have assumed. In accretion disks around massive black holes in active galactic nuclei, the disk will probably become optically thin before the outer disk boundary is reached, whereas in X-ray binaries there will be an outer disk region (outside the circularization radius) in which the inflow velocity is zero but angular momentum is still transported. We show that in both these cases the solutions are similar to the zero-crossing eigenfunctions.

Maloney, Philip R.↗

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

Characterizing the seismic response of a molten salt nuclear reactor

Abstract One Generation IV nuclear reactor, which uses a fluoride salt as a coolant, graphite reflector blocks as a moderator, and circulating buoyant TRISO pebbles as fuel is at an advanced stage of development. To characterize the seismic behavior of components of this reactor, validate numerical models for analysis, and develop recommendations for design, a set of earthquake‐simulator experiments on a scaled model of the reactor vessel and its internals was executed on a six‐degree‐of‐freedom earthquake simulator. The model was seismically isolated at its base using two types of spherical sliding bearings. The scaled model involved representations of the prototype reactor vessel, core barrel, reflector blocks, coolant, and spherical fuel pebbles. The material and geometric properties of different test components were selected based on a dynamic similitude scaling analysis and an approximate length scale of 0.4. Four sets of three‐component earthquake motions were used as inputs for testing. Instrumentation on the test specimen recorded the dynamic responses of the outer vessel, core barrel, and reflector‐block assembly, the hydrodynamic responses (sloshing and hydrodynamic pressure) of the liquid coolant, pebble consolidation under earthquake shaking, and the behavior of the isolation systems. This paper describes the design of the experiments and presents key results from the tests. The dynamic responses of the outer vessel, core barrel, and the reflector blocks revealed that the components responded as a unit for the intense shaking used in the experiments. The sloshing response of the fluid in a thin annulus near the perimeter of the vessel was heavily damped. The change in the packing fraction of the pebble bed under repeated, intense 3D earthquake shaking was less than 3%. Seismically isolating the vessel substantially reduced demands on its internal components.

Engineering↗

Direct stellarator coil optimization for nested magnetic surfaces with precise quasi-symmetry

We present a robust optimization algorithm for the design of electromagnetic coils that generate vacuum magnetic fields with nested flux surfaces and precise quasi-symmetry. The method is based on a bilevel optimization problem, where the outer coil optimization is constrained by a set of inner least squares optimization problems whose solutions describe magnetic surfaces. The outer optimization objective targets coils that generate a field with nested magnetic surfaces and good quasi-symmetry. The inner optimization problems identify magnetic surfaces when they exist, and approximate surfaces in the presence of magnetic islands or chaos. We show that this formulation can be used to heal islands and chaos, thus producing coils that result in magnetic fields with precise quasi-symmetry. We show that the method can be initialized with coils from the traditional two-stage coil design process, as well as coils from a near-axis expansion optimization. We present a numerical example where island chains are healed and quasi-symmetry is optimized up to surfaces with aspect ratio 6. Another numerical example illustrates that the aspect ratio of nested flux surfaces with optimized quasi-symmetry can be decreased from 6 to approximately 4. In conclusion, the last example shows that our approach is robust and a cold-start using coils from a near-axis expansion optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nucleosynthesis of an 11.8 M ⊙ Supernova with 3D Simulation of the Inner Ejecta: Overall Yields and Implications for Short-lived Radionuclides in the Early Solar System

Based on a 3D supernova simulation of an 11.8 M ⊙ progenitor model with an initial solar composition, we study the nucleosynthesis using tracers covering the innermost 0.1 M ⊙ of the ejecta. These ejecta are mostly proton-rich and contribute significant amounts of 45 Sc and 64 Zn. The production of heavier isotopes is sensitive to the electron fraction and hence the neutrino emission from the proto-neutron star. The yields of these isotopes are rather uncertain due to the approximate neutrino transport used in the simulation. In order to obtain the total yields for the whole supernova, we combine the results from the tracers with those for the outer layers from a suitable 1D model. Using the yields of short-lived radionuclides (SLRs), we explore the possibility that an 11.8 M ⊙ supernova might have triggered the formation of the solar system and provided some of the SLRs measured in meteorites. In particular, we discuss two new scenarios that can account for at least the data on 41 Ca, 53 Mn, and 60 Fe without exceeding those on the other SLRs.

Astronomy & Astrophysics↗