Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling statistics”

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 613 records · Page 34

What’s in My Toolkit? A Review of Technologies for Assessing Changes in Habitats Caused by Marine Energy Development

Marine energy devices are installed in highly dynamic environments and have the potential to affect the benthic and pelagic habitats around them. Regulatory bodies often require baseline characterization and/or post-installation monitoring to determine whether changes in these habitats are being observed. However, a great diversity of technologies is available for surveying and sampling marine habitats, and selecting the most suitable instrument to identify and measure changes in habitats at marine energy sites can become a daunting task. We conducted a thorough review of journal articles, survey reports, and grey literature to extract information about the technologies used, the data collection and processing methods, and the performance and effectiveness of these instruments. We examined documents related to marine energy development, offshore wind farms, oil and gas offshore sites, and other marine industries around the world over the last 20 years. A total of 120 different technologies were identified across six main habitat categories: seafloor, sediment, infauna, epifauna, pelagic, and biofouling. The technologies were organized into 12 broad technology classes: acoustic, corer, dredge, grab, hook and line, net and trawl, plate, remote sensing, scrape samples, trap, visual, and others. Visual was the most common and the most diverse technology class, with applications across all six habitat categories. Technologies and sampling methods that are designed for working efficiently in energetic environments have greater success at marine energy sites. In addition, sampling designs and statistical analyses should be carefully thought through to identify differences in faunal assemblages and spatiotemporal changes in habitats.

16 TIDAL AND WAVE POWER↗

Influence of geometry on columnar to equiaxed transition during electron beam powder bed fusion of IN718

Correlation between spot-melt scan parameters (linear spot-density aka areal energy density), build geometry, and solidification microstructure evolution (columnar vs equiaxed) in a powder bed fusion technology is investigated. It is shown that to maintain the equiaxed solidification microstructure evolution in electron beam powder bed additive manufacturing (AM), the areal energy density per layer needs to be scaled with respect to the 2D cross-sectional area of the layer being melted. Samples with two different cross-sectional areas (40 × 40 mm and 20 × 20 mm) have been fabricated with varying areal energy densities. For a given square cross-section (20 × 20 mm), increasing the areal energy density (4.8 MJ/sq.m to 14.7 MJ/sq.m) transitioned the solidification microstructure from columnar to equiaxed. The observed microstructure data (Electron Back Scattered Diffraction - EBSD) is quantified by calculating the principal component (PC) score using a spatial statistics methodology. The sample with equiaxed grains is found to have a low PC score while the sample with columnar grain had a high PC score. A semi-analytical model is used to simulate the heat transfer and the local solidification conditions as a function of processing parameters (linear spot-density). The result from the heat transfer model is correlated with previously quantified microstructure data. Space-Time analysis of the melt pattern is done and correlated with the observed microstructure. In addition, from the findings, appropriate parameters have been used to additively manufacture a turbine blade with site-specific or hybrid solidification microstructure (traditional fabrication possible via a patented method of localized cold working and heat treatment).

36 MATERIALS SCIENCE↗

Prediction of Probabilistic Detonation Threshold via Millimeter-Scale Microstructure-Explicit and Void-Explicit Simulations

We present an approach and relevant models for predicting the probabilistic shoc-to-detonation transition (SDT) behavior and Pop plot (PP) of heterogeneous energetic materials (HEM) via mesoscopic microstructure-explicit (ME) and void explicit (VE) simulations at the millimeter (mm) sample size scale. Although the framework here is general, the particular material considered in this paper is pressed Octahydro-1,3,5,7-tetranitro-1,2,3,5-tetrazocine (HMX). To systematically delineate the effects of material heterogeneities, four material cases are considered. These cases are homogeneous material, material with granular microstructure but no voids, homogeneous material with voids, and material with both granular microstructure and voids. Statistically equivalent microstructure sample sets (SEMSS) are generated and used. Eulerian hydrocode simulations explicitly resolve the material heterogeneities, voids, and the coupled mechanical-thermal-chemical processes. In particular, it is found that both microstructure and voids strongly influence the SDT behavior and PP. The effects of different combinations of microstructure heterogeneity and voids on the SDT process and PP are quantified and rank-ordered. The overall framework uses the Mie–Grüneisen equation of state and a history variable reactive burn model (HVRB). A novel probabilistic representation for quantifying the PP is developed, allowing the calculation of (1) the probability of observing SDT at a given combination of shock pressure and run distance, (2) the run-distance to detonation under a given combination of shock pressure and prescribed probability, and (3) the shock pressure required for achieving SDT at a given run distance with a prescribed probability. The results are in agreement with general trends in experimental data in the literature.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A KINE-CHEMICAL INVESTIGATION OF THE AB DOR MOVING GROUP 'STREAM'

The AB Dor Moving Group consists of a 'nucleus' of {approx}10 stars at d {approx_equal} 20 pc, along with dozens of purported 'stream' members distributed across the sky. We perform a chemical and kinematic analysis of a subsample of AB Dor stream stars to test whether they constitute a physical stellar group. We use the NEMO Galactic kinematic code to investigate the orbits of the stream members, and perform a chemical abundance analysis using high resolution spectra taken with the Magellan Clay 6.5 m telescope. Using a {chi}{sup 2} test with the measured abundances for 10 different elements, we find that only half of the purported AB Dor stream members could possibly constitute a statistically chemically homogeneous sample. Some stream members with three-dimensional velocities were hundreds of parsecs from the AB Dor nucleus {approx}10{sup 8} yr ago, and hence were unlikely to share a common origin. We conclude that the published lists of AB Dor moving group stream members are unlikely to represent the dispersed remnant of a single star formation episode. A subsample of the stream stars appears to be both statistically chemically homogeneous and in the vicinity of the AB Dor nucleus at birth. Their mean metallicity is [Fe/H] = 0.02 {+-} 0.02 dex, which we consider representative for the AB Dor group. Finally, we report a strong lower limit on the age of the AB Dor nucleus of >110 Myr based on the pre-main sequence contraction times for K-type members which have reached the main sequence.

79 ASTRONOMY AND ASTROPHYSICS↗

Combining multispectral and ancillary data in remote sensing using interval-valued probabilities

This paper presents a method for classifying multisource data in remote sensing and geographic information systems using interval-valued probabilities. In this method, each data source is considered as an information source which provides a body of statistical evidence. In order to integrate information obtained from multiple data sources, the method adopts Dempster's rule for combining multiple bodies of evidence. Preliminary experiments have been undertaken to illustrate the use of the method in a supervised ground-cover classification on multispectral data combined with digital elevation data. They demonstrate the ability of the method in capturing information provided by inexact and incomplete evidence when there are not enough training samples to estimate statistical parameters.

Kim, H.↗

MAGNITUDE GAP STATISTICS AND THE CONDITIONAL LUMINOSITY FUNCTION

In a recent preprint, Hearin et al. (H12) suggest that the halo mass-richness calibration of clusters can be improved by using the difference in the magnitude of the brightest and the second brightest galaxy (magnitude gap) as an additional observable. They claim that their results are at odds with the results from Paranjape and Sheth (PS12) who show that the magnitude distribution of the brightest and second brightest galaxies can be explained based on order statistics of luminosities randomly sampled from the total galaxy luminosity function. We find that a conditional luminosity function (CLF) for galaxies which varies with halo mass, in a manner which is consistent with existing observations, naturally leads to a magnitude gap distribution which changes as a function of halo mass at fixed richness, in qualitative agreement with H12. We show that, in general, the luminosity distribution of the brightest and the second brightest galaxy depends upon whether the luminosities of galaxies are drawn from the CLF or the global luminosity function. However, we also show that the difference between the two cases is small enough to evade detection in the small sample investigated by PS12. This shows that the luminosity distribution is not the appropriate statistic to distinguish between the two cases, given the small sample size. We argue in favor of the CLF (and therefore H12) based upon its consistency with other independent observations, such as the kinematics of satellite galaxies, the abundance and clustering of galaxies, and the galaxy-galaxy lensing signal from the Sloan Digital Sky Survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Evolutions and Calibrations of Long Gamma-Ray-burst Luminosity Correlations Revisited

Luminosity correlations of long gamma-ray bursts (GRBs) are extensively proposed as an effective complementarity to trace the Hubble diagram of the universe at high redshifts, which is of great importance to explore properties of dark energy. Recently, several empirical luminosity correlations have been statistically proposed from GRB observations. However, to treat a GRB as the distance indicator, there are two key issues: the redshift evolution of luminosity correlations and their calibrations. In this paper, we choose the Amati relation, the correlation between the peak spectra energy and the equivalent isotropic energy of GRBs (E{sub p}--E{sub iso}) as an example and find that the current GRB data set implies that there could be a evolution of the luminosity correlation with respect to the redshift. Therefore, we propose an extended Amati relation with two extra redshift-dependent terms to correct the redshift evolution of the GRB relation. Second, we carefully check the reliability of the calibration method using the low-redshift GRB data. Importantly, we find that the low-redshift calibration method does not take whole correlations between Ω{sub m} and coefficients into account. Neglecting this correlation information can break the degeneracies and obtain the biased constraint on Ω{sub m}, which is very sensitive to values of parameters for the calibration. A small shift in the parameters of the “calibrated” relation could significantly change the final constraint on Ω{sub m} in the low-redshift calibration method. Finally, we simulate several GRB samples with different statistical errors and find that, in order to correctly recover the fiducial value of Ω{sub m}, using the low-redshift calibration method, we need a large number of GRB samples with high precisions.

79 ASTRONOMY AND ASTROPHYSICS↗

OzDES reverberation mapping program: Lag recovery reliability for 6-yr C iv analysis

ABSTRACT We present the statistical methods that have been developed to analyse the OzDES reverberation mapping sample. To perform this statistical analysis we have created a suite of customizable simulations that mimic the characteristics of each source in the OzDES sample. These characteristics include: the variability in the photometric and spectroscopic light curves, the measurement uncertainties, and the observational cadence. By simulating the sources in the OzDES sample that contain the C iv emission line, we developed a set of criteria that rank the reliability of a recovered time-lag depending on the agreement between different recovery methods, the magnitude of the uncertainties, and the rate at which false positives were found in the simulations. These criteria were applied to simulated light curves and these results used to estimate the quality of the resulting Radius–Luminosity relation. We grade the results using three quality levels (gold, silver, and bronze). The input slope of the R–L relation was recovered within 1σ for each of the three quality samples, with the gold standard having the lowest dispersion with a recovered a R–L relation slope of 0.454 ± 0.016 with an input slope of 0.47. Future work will apply these methods to the entire OzDES sample of 771 AGN.

79 ASTRONOMY AND ASTROPHYSICS↗

A CATALOG OF ROTATION AND ACTIVITY IN EARLY-M STARS

We present a catalog of rotation and chromospheric activity in a sample of 334 M dwarfs of spectral types M0-M4.5 populating the parameter space around the boundary to full convection. We obtain high-resolution optical spectra for 206 targets and determine projected rotational velocity, vsin i, and H{alpha} emission. The data are combined with measurements of vsin i in field stars of the same spectral type from the literature. Our sample adds 157 new rotation measurements to the existing literature and almost doubles the sample of available vsin i. The final sample provides a statistically meaningful picture of rotation and activity at the transition to full convection in the solar neighborhood. We confirm a steep rise in the fraction of active stars at the transition to full convection known from earlier work. In addition, we see a clear rise in rotational velocity in the same stars. In very few stars, no chromospheric activity but a detection of rotational broadening is reported. We argue that all of them are probably spurious detections; we conclude that in our sample all significantly rotating stars are active, and all active stars are significantly rotating. The rotation-activity relation is valid in partially and in fully convective stars. Thus, we do not observe any evidence for a transition from a rotationally dominated dynamo in partially convective stars to a rotation-independent turbulent dynamo in fully convective stars; turbulent dynamos in fully convective stars of spectral types around M4 are still driven by rotation. Finally, we compare projected rotational velocities of 33 stars to rotational periods derived from photometry in the literature and determine inclinations for a few of them.

79 ASTRONOMY AND ASTROPHYSICS↗

Equibiaxial Flexure Strength of a Superfine-Grained Nuclear Graphite

The strength of advanced graphite is reported in accordance with ASTM D7846-16 using a two-parameter Weibull distribution for uniaxial strength testing. The rule of thumb for Weibull distributions is to use a minimum of 30 strength measurements to have a high level of confidence in the Weibull characteristic strength and the Weibull modulus. These large sample sets for statistical confidence are easily obtained for as-manufactured graphite, but determination of the Weibull two-factor parameters on graphite that has been exposed to neutron irradiation is nearly impossible due to irradiation testing restraints on specimen size and space limitations in test reactors for the accommodation of specimens. These restrictions have resulted in irradiation programs that use sub-sized specimens to measure uniaxial strength change, but the specimen geometries still hinder the number of replicate specimens for each irradiation condition. In contrast, this work used an advanced ceramic standard, ASTM C1499-05, as the foundation for equibiaxial strength testing of a superfine-grained nuclear graphite to investigate the irradiation-induced strength change and any changes to the two-parameter Weibull distribution. The goal was to investigate the effect of specimen size on the two-parameter Weibull distribution parameters in support of the small specimens and geometry desired for irradiation. A study of the effect of a reduced sample population was also undertaken to investigate the reliability of using the two-parameter Weibull distribution to study changes to the materials properties caused by neutron irradiation damage. The results of this work suggest that these equibiaxial specimens would be ideal for statistically significant studies of the effects of irradiation on the Weibull characteristic strength (but not the modulus) and for providing a method for surveillance specimen campaigns for future operating commercial reactors.

Campbell, Anne↗

MULTIPLICITY AMONG F-TYPE STARS. II

In continuation of our previous study we present an updated census of new companions and model atmosphere analyses for some 50 southern dwarfs, mostly in the mass range 0.90⩽M⩽1.10 M{sub ⊙}. For the common-proper-motion companions μ Vir B, HR 2225 B, HD 67199 B, and HD 114853 B, we confirm their physical association from their radial velocities. We report the discovery of the F-type visual binary α For as a field blue straggler and confirm (ζ Ret, HR 5864) or identify (HD 67199, HR 4013, HR 8843) another five mass transfer systems or candidates. For the F stars τ{sup 1} Eri and 111 Tau, we present 10σ and 7σ cases for astrometric binaries by virtue of the very accurate van Leeuwen Hipparcos parallaxes. Following the work of Shaya and Olling, we suggest the F-type star ι Vir to be a wide (0.37 pc) hierarchical quadruple system. We confirm the visual binary NLTT 23781/2 as a common-proper-motion object to the very wide (0.54 pc) F star 40 Leo, but discard the G star HD 128987 as an ultra-wide (1.01 pc) physical companion to the α Lib quadruple system on account of a diverse metallicity. The improved statistics of our sample establishes the previously discovered positive correlation of stellar multiplicities with primary mass. For the F star multiplicity census in the mass range 1.10⩽M⩽1.70 M{sub ⊙}, we find that at least a quarter consists of triple or higher level systems and at least two out of three F stars are non-single.

79 ASTRONOMY AND ASTROPHYSICS↗

Aeroelastic Uncertainty Quantification Studies Using the S4T Wind Tunnel Model

This paper originates from the joint efforts of an aeroelastic study team in the Applied Vehicle Technology Panel from NATO Science and Technology Organization, with the Task Group number AVT-191, titled "Application of Sensitivity Analysis and Uncertainty Quantification to Military Vehicle Design." We present aeroelastic uncertainty quantification studies using the SemiSpan Supersonic Transport wind tunnel model at the NASA Langley Research Center. The aeroelastic study team decided treat both structural and aerodynamic input parameters as uncertain and represent them as samples drawn from statistical distributions, propagating them through aeroelastic analysis frameworks. Uncertainty quantification processes require many function evaluations to asses the impact of variations in numerous parameters on the vehicle characteristics, rapidly increasing the computational time requirement relative to that required to assess a system deterministically. The increased computational time is particularly prohibitive if high-fidelity analyses are employed. As a remedy, the Istanbul Technical University team employed an Euler solver in an aeroelastic analysis framework, and implemented reduced order modeling with Polynomial Chaos Expansion and Proper Orthogonal Decomposition to perform the uncertainty propagation. The NASA team chose to reduce the prohibitive computational time by employing linear solution processes. The NASA team also focused on determining input sample distributions.

Nikbay, Melike↗

Strong effect of liquid Fe–S on elastic wave velocity of olivine aggregate: Implication for the low velocity anomaly at the base of the lunar mantle

The origin of the seismic low velocity zone (LVZ) observed at the base of the lunar mantle is important for understanding the nature of the lunar mantle. Liquid Fe–S is considered as a possible candidate of the cause of the LVZ. However, effect of liquid Fe–S on elastic wave velocity of mantle rock has been considered to be minor, according to theoretical estimation based on the high dihedral angle of liquid Fe–S in olivine aggregate observed in previous microstructural analysis of polished sample cross sections. Here, in this work, we carry out direct measurements of elastic wave velocities of liquid Fe–S-bearing olivine aggregate and report strong reduction of the elastic wave velocities in the presence of liquid Fe–S. The effects of liquid Fe–S in olivine aggregates are much more pronounced than the theoretical estimation. Three-dimensional X-ray microtomography analysis reveals that the geometry of liquid Fe–S strongly depends on the size, with larger liquid Fe–S blobs having low aspect ratios, which may predominantly affect elastic wave velocity. Since smaller spherical Fe–S blobs are much more abundant than the larger elongated blobs, conventional dihedral angle analyses based on scanning electron microscope images tend to statistically over-sample smaller blobs and under-sample larger blobs. Our results suggest that the LVZ at the base of the lunar mantle can be explained by the presence of 4.3–5.0 vol.% liquid Fe–S. The LVZ may be a hidden reservoir of highly siderophile elements (HSEs), which may explain the low abundance of HSEs in the lunar mantle compared to the Earth.

58 GEOSCIENCES↗

Sampling for Confidence: Guidance for Evaluating the Quality of Canister Sampling Plans

There is a need for an inspection frequency that demonstrates reasonable assurance that the confinement function of dual-purpose-canisters (DPCs) is maintained as these DPCs are returned to storage service following transportation. The DPCs contain spent fuel from shutdown reactor sites. The currently-proposed 100% DPC inspection frequency is overly conservative and unnecessary to demonstrate reasonable assurance that the confinement function of all DPCs is maintained for return to storage service. Testing all DPCs does not appropriately reflect the very low risk of leakage (i.e., low probability based on design and low consequences based on analysis). Currently, there are 1066 DPCs forecast to be loaded at reactor sites which either are shutdown or plan to cease operations by 2025. These DPCs may be transported by rail in sets of 3 to 5 DPCs per shipment. Assuming each rail consist originates at a single reactor site, an alternative proposal to test one canister/consist/site would result in 356 to 214 canisters, respectively, being tested given the population of 1,066 canisters This Savannah River National Laboratory (SRNL) report outlines statistically-based, random sampling plans that provide guidance on the number of samples of these canisters that are needed to establish bounds on the quality of the outcome of such sampling. The quality of each proposed sampling plan is expressed in the form of a statement that the number of undetected, unacceptable canisters in the population is less than some number, call it D, at a high level of confidence. In addition, evaluations using a similar approach for random sampling at the proposed levels of 214 and 356 canisters are provided.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Sampling for Confidence: Guidance for Evaluating the Quality of Canister Sampling Plans

There is a need for an inspection frequency that demonstrates reasonable assurance that the confinement function of dual-purpose-canisters (DPCs) is maintained as these DPCs are returned to storage service following transportation. The DPCs contain spent fuel from shutdown reactor sites. The currently-proposed 100% DPC inspection frequency is overly conservative and unnecessary to demonstrate reasonable assurance that the confinement function of all DPCs is maintained for return to storage service. Testing all DPCs does not appropriately reflect the very low risk of leakage (i.e., low probability based on design and low consequences based on analysis). Currently, there are 1066 DPCs forecast to be loaded at reactor sites which either are shutdown or plan to cease operations by 2025. These DPCs may be transported by rail in sets of 3 to 5 DPCs per shipment. Assuming each rail consist originates at a single reactor site, an alternative proposal to test one canister/consist/site would result in 356 to 214 canisters, respectively, being tested given the population of 1,066 canisters This Savannah River National Laboratory (SRNL) report outlines statistically-based, random sampling plans that provide guidance on the number of samples of these canisters that are needed to establish bounds on the quality of the outcome of such sampling. The quality of each proposed sampling plan is expressed in the form of a statement that the number of undetected, unacceptable canisters in the population is less than some number, call it D, at a high level of confidence. In addition, evaluations using a similar approach for random sampling at the proposed levels of 214 and 356 canisters are provided. This is a technical report that does not take into account the contractual limitations under the Standard Contract for Disposal of Spent Nuclear Fuel and/or High-Level Radioactive Waste (Standard Contract) (10 CFR Part 961). Under the provisions of the Standard Contract, DOE does not consider spent nuclear fuel in multi-assembly canisters to be an acceptable waste form, absent a mutually agreed upon contract amendment. To the extent discussions or recommendations in this report conflict with the provisions of the Standard Contract, the Standard Contract provisions prevail.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Optimizing human activity patterns using global sensitivity analysis

Implementing realistic activity patterns for a population is crucial for modeling, for example, disease spread, supply and demand, and disaster response. Using the dynamic activity simulation engine, DASim, we generate schedules for a population that capture regular (e.g., working, eating, and sleeping) and irregular activities (e.g., shopping or going to the doctor). We use the sample entropy (SampEn) statistic to quantify a schedule’s regularity for a population. We show how to tune an activity’s regularity by adjusting SampEn, thereby making it possible to realistically design activities when creating a schedule. The tuning process sets up a computationally intractable high-dimensional optimization problem. To reduce the computational demand, we use Bayesian Gaussian process regression to compute global sensitivity indices and identify the parameters that have the greatest effect on the variance of SampEn. Here we use the harmony search (HS) global optimization algorithm to locate global optima. Our results show that HS combined with global sensitivity analysis can efficiently tune the SampEn statistic with few search iterations. We demonstrate how global sensitivity analysis can guide statistical emulation and global optimization algorithms to efficiently tune activities and generate realistic activity patterns. Finally, though our tuning methods are applied to dynamic activity schedule generation, they are general and represent a significant step in the direction of automated tuning and optimization of high-dimensional computer simulations.

97 MATHEMATICS AND COMPUTING↗

Addendum to Sampling and Analysis Plan (SAP) for Assessment of LANL-Derived Residual Radionuclides in Soils within Tract A-16-d for Land Conveyance and Transfer for Sewage Treatment Facility Area

This report summarizes the sampling design used, associated statistical assumptions, as well as general guidelines for conducting post-sampling data analysis. Sampling plan components presented here include how many sampling locations to choose and where within the sampling area to collect those samples. The type of medium to sample (i.e., soil, groundwater, etc.) and how to analyze the samples (in-situ, fixed laboratory, etc.) are addressed in other sections of the sampling plan.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The moving group targets of the seeds high-contrast imaging survey of exoplanets and disks: Results and observations from the first three years

We present results from the first three years of observations of moving group (MG) targets in the Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS) high-contrast imaging survey of exoplanets and disks using the Subaru telescope. We achieve typical contrasts of ∼10{sup 5} at 1'' and ∼10{sup 6} beyond 2'' around 63 proposed members of nearby kinematic MGs. We review each of the kinematic associations to which our targets belong, concluding that five, β Pictoris (∼20 Myr), AB Doradus (∼100 Myr), Columba (∼30 Myr), Tucana-Horogium (∼30 Myr), and TW Hydrae (∼10 Myr), are sufficiently well-defined to constrain the ages of individual targets. Somewhat less than half of our targets are high-probability members of one of these MGs. For all of our targets, we combine proposed MG membership with other age indicators where available, including Ca II HK emission, X-ray activity, and rotation period, to produce a posterior probability distribution of age. SEEDS observations discovered a substellar companion to one of our targets, κ And, a late B star. We do not detect any other substellar companions, but do find seven new close binary systems, of which one still needs to be confirmed. A detailed analysis of the statistics of this sample, and of the companion mass constraints given our age probability distributions and exoplanet cooling models, will be presented in a forthcoming paper.

79 ASTRONOMY AND ASTROPHYSICS↗