Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random projection”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Kinetics of Cyclic Oxidation and Cracking and Finite Element Analysis of MA956 and Sapphire/MA956 Composite System

Sapphire fiber-reinforced MA956 composites hold promise for significant weight savings and increased high-temperature structural capability, as compared to unreinforced MA956. As part of an overall assessment of the high-temperature characteristics of this material system, cyclic oxidation behavior was studied at 1093 C and 1204 C. Initially, both sets of coupons exhibited parabolic oxidation kinetics. Later, monolithic MA956 exhibited spallation and a linear weight loss, whereas the composite showed a linear weight gain without spallation. Weight loss of the monolithic MA956 resulted from the linking of a multiplicity of randomly oriented and closely spaced surface cracks that facilitated ready spallation. By contrast, cracking of the composite's oxide layer was nonintersecting and aligned nominally parallel with the orientation of the subsurface reinforcing fibers. Oxidative lifetime of monolithic MA956 was projected from the observed oxidation kinetics. Linear elastic, finite element continuum, and micromechanics analyses were performed on coupons of the monolithic and composite materials. Results of the analyses qualitatively agreed well with the observed oxide cracking and spallation behavior of both the MA956 and the Sapphire/MA956 composite coupons.

Lee, Kang N.↗
NASA NTRS2005

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

Previous work at NASA Langley Research Center (LaRC) involved fabrication and testing of composite beams with embedded, pre-strained shape memory alloy (SMA) ribbons within the beam structures. That study also provided comparison of experimental results with numerical predictions from a research code making use of a new thermoelastic model for shape memory alloy hybrid composite (SMAHC) structures. The previous work showed qualitative validation of the numerical model. However, deficiencies in the experimental-numerical correlation were noted and hypotheses for the discrepancies were given for further investigation. The goal of this work is to refine the experimental measurement and numerical modeling approaches in order to better understand the discrepancies, improve the correlation between prediction and measurement, and provide rigorous quantitative validation of the numerical analysis/design tool. The experimental investigation is refined by a more thorough test procedure and incorporation of higher fidelity measurements such as infrared thermography and projection moire interferometry. The numerical results are produced by a recently commercialized version of the constitutive model as implemented in ABAQUS and are refined by incorporation of additional measured parameters such as geometric imperfection. Thermal buckling, post-buckling, and random responses to thermal and inertial (base acceleration) loads are studied. The results demonstrate the effectiveness of SMAHC structures in controlling static and dynamic responses by adaptive stiffening. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗
NASA NTRS2016

Attentional Modulation of Eye Torsion Responses

Eye movements generally have both reflexive and voluntary aspects, but torsional eye movements are usually thought of as a reflexive response to image rotation around the line of sight (torsional OKN) or to head roll (torsional VOR). In this study we asked whether torsional responses could be modulated by attention in a case where two stimuli rotated independently, and whether attention would influence the latency of responses. The display consisted of rear-projected radial "pinwheel" gratings, with an inner annulus segment extending from the center to 22 degrees eccentricity, and an outer annulus segment extending from 22 degrees out to 45 degrees eccentricity. The two segments rotated around the center in independent random walks, stepping randomly 4 degrees clockwise or counterclockwise at 60 Hz. Subjects were asked to attend to one or the other while keeping fixation steady at the center of the display. To encourage attention on one or the other segment of the display, subjects were asked to move a joystick in synchrony with the back and forth rotations of one part of the image while ignoring the other. Eye torsion was recorded with the scleral search coil technique, sampled at 500 Hz. All four subjects showed roughly 50% stronger torsion responses to the attended compared to unattended segments. Latency varied from 100 to 150 msec across subjects and was unchanged by attention. These findings suggest that attention can influence eye movement responses that are not typically under voluntary control.

attention↗
NASA NTRS2016

Attentional Modulation of Eye Torsion Responses

Eye movements generally have both reflexive and voluntary aspects, but torsional eye movements are usually thought of as a reflexive response to image rotation around the line of sight (torsional OKN) or to head roll (torsional VOR). In this study we asked whether torsional responses could be modulated by attention in a case where two stimuli rotated independently, and whether attention would influence the latency of responses. The display consisted of rear-projected radial pinwheel gratings, with an inner annulus segment extending from the center to 22 degrees eccentricity, and an outer annulus segment extending from 22 degrees out to 45 degrees eccentricity. The two segments rotated around the center in independent random walks, stepping randomly 4 degrees clockwise or counterclockwise at 60 Hz. Subjects were asked to attend to one or the other while keeping fixation steady at the center of the display. To encourage attention on one or the other segment of the display, subjects were asked to move a joystick in synchrony with the back and forth rotations of one part of the image while ignoring the other. Eye torsion was recorded with the scleral search coil technique, sampled at 500 Hz. All four subjects showed roughly 50 stronger torsion responses to the attended compared to unattended segments. Latency varied from 100 to 150 msec across subjects and was unchanged by attention. These findings suggest that attention can influence eye movement responses that are not typically under voluntary control.

ocular torsion↗
NASA NTRS2022

Environmental Testing of the OVEN System for Lunar Water Extraction and Prospecting

Introduction: The presence of water ice in permanently shadowed regions on the lunar surface [1] may enable a sustained human presence on the Moon with minimal need for consumables. However, in order to develop a long term utilization plan that includes the usage of in-situ water we must first understand the abundance, stratigraphy and distribution of this re-source. Multiple space agencies currently have plans for lunar water prospecting missions. The Optimized Volatile Extraction Node (OVEN) was designed for water prospecting missions that require samples to be weighed, sealed, and heated as the means of determining water concentration. This method of water quantification necessitates a fair amount of automation, so a rigorous environmental test program was performed in order to build confidence in the performance of the OVEN design. The work presented here describes the OVEN environmental test program as well as ongoing efforts to improve on the design. Vibration: The OVEN participated in two rounds of random vibration tests. The first test was a stand-alone test performed at the Energy Systems Test Area of the Johnson Space Center. The second test was an integrated test with the mobile platform developed for the Resource Prospector project. The OVEN survived both tests without damage, but the tests did provide valuable lessons learned with regards to specific operations. Thermal Vacuum: The OVEN was successfully demonstrated at a temperature range of -50 to 75 C in a thermal vacuum chamber. The need to heat motor gearboxes at lower temperatures was predetermined so this test program was completed by implementing a method of gearbox heating that used the existing circuitry within the motors. Dust: A custom dynamometer was built in order to determine the torque required to move the various mechanisms within the OVEN at a range of temperatures. The OVEN system was coated with lunar dust simulant in order to determine mechanism torques under a worst-case operating condition. Sublimation: Sublimation losses within the OVEN were quantified through a series of test configurations, including an integrated test in a thermal vacuum chamber at Glenn Research Center [2]. Current Work: The OVEN subsystem is currently not a component of any existing prospecting missions, but work continues that will take the lessons learned from previous environmental tests and improve on the design in order to be considered for future prospecting opportunities. References: [1] Colaprete, A., Schultz, P., Heldmann, J., Wooden, D., Shirley, M., Ennico, K., ... & Sollitt, L. (2010). Detection of water in the LCROSS ejecta plume. science, 330(6003), 463-468. [2] Kleinhenz, J., Smith, J., Roush, T., Colaprete, A., Zacny, K., Paulsen, G., ... & Paz, A. (2018). Volatiles Loss from water bearing regolith simulant at Lunar Environments. In Earth and Space 2018: Engineering for Extreme Environments (pp. 454-466). Reston, VA: American Society of Civil Engineers.

A Paz↗
NASA NTRS2024

Jefferson County Ecological Conservation: Quantifying the Effects of Hydrologic Restoration in the Camas National Wildlife Refuge and Mud Lake Wildlife Management Area

Wetlands in the Camas National Wildlife Refuge and Mud Lake Wildlife Management Area, a rare landscape feature in the Intermountain West, are crucial for migratory birds along the Pacific Flyway, yet these wetlands have experienced a noticeable decline in extent and inundation over the last 40+ years. Partners at the U.S. Fish and Wildlife Service and Idaho Department of Fish and Game have ongoing restoration efforts that are yet to be quantified. To assist these partners, we used NASA Earth observations to quantify change in wetland extent and determine if the restoration projects met the intended impact. The QA_PIXEL band from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) was used to identify water. Landsat 8 OLI imagery was also used for classifying landcover in ArcGIS Pro for 2016 and 2020 and to forecast landcover in 2060 using TerrSet’s Land Change Modeler. We also compared lidar data from 2011 and 2019 to detect changes in vegetation height, ground elevation, and surface water levels, and Sentinel-2 Multispectral Instrument (MSI) imagery to detect changes in vegetation health. Classification results indicated overall accuracies over 84% and Kappa indices above 0.80, well surpassing random classification performance. From 2016 to 2020 our classifications showed an 825 acre increase in wetland extent, and our forecasted 2060 model showed a 2795 acre decrease. Combining our partners’ knowledge of the area with our analyses, we created a remote sensing workflow to enhance future monitoring and decision-making. This study demonstrated the feasibility of NASA Earth observations in quantifying wetlands, informing restoration projects, and supporting avian conservation.

Cassidy Bromka↗
NASA NTRSDate not supplied

Characterizing Spatiotemporal Uncertainty in Interpolated Meteorological Data

Interpolated meteorological data invariably contain errors. These errors have structure in time and space, particularly autocorrelation, which can cause the effects of errors to compound when model outputs are aggregated temporally or spatially. One way to account for this uncertainty is with a probabilistic model from which samples can be drawn that are coherent with respect to underlying spatial and temporal covariance structure. This work describes a probabilistic method for spatial interpolation of point-wise meteorological time series. Observational data from weather stations are generally sparse in space and dense in time (but sometimes missing). The method works by projecting time series onto orthogonal basis vectors and spatially interpolating each resulting component independently. Under suitable assumptions, and data transformations to better satisfy those assumptions, Gaussian process regression provides a complete description of the joint predictive distribution over a Gaussian random field. Spatiotemporally coherent realizations are generated as the sum of conditional (spatial) simulations of each orthogonal (temporal) component. Data-derived and generic orthogonal bases are considered. In addition to spatial interpolation, imputation of missing observational data is examined. The method is applied using near-surface air temperature over the Western United States and validated by comparing theoretical versus actual coverage of predictive distributions and analyzing the degree to which spatial and temporal covariance structure is reproduced. Computational considerations, relating to conditional simulation of random fields, are also addressed.

Conor T Doherty↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗
NASA NTRS1981

Development of alternative data analysis techniques for improving the accuracy and specificity of natural resource inventories made with remote sensing data

The applicability of digital LANDSAT data in updating the Minnesota Lake Land Management Information System and in assessing the trophic status of lakes was investigated. Data from various combinations of training, classification, and geometric correction techniques were compared to a photointerpreted reference data set consisting of ground-verified samples and six randomly selected photographs covering approximately one quarter of the study area. For a reconnaissance inventory, the training via polygons selected from aerial photographs with a canonical analysis minimum distance classifier is the most accurate and efficient analysis technique. The applications test project resulted in disappointing thematic classification results per se, but did contribute to the development of mechanisms for implementing a LANDSAT data processing capability in the State Planning Agency on a permanent basis.

Lillesand, T. M.↗
NASA NTRS2011

Application of the Semi-Empirical Force-Limiting Approach for the CoNNeCT SCAN Testbed

The semi-empirical force-limited vibration method was developed and implemented for payload testing to limit the structural impedance mismatch (high force) that occurs during shaker vibration testing. The method has since been extended for use in analytical models. The Space Communications and Navigation Testbed (SCAN Testbed), known at NASA Glenn Research Center (GRC) as, the Communications, Navigation, and Networking re-Configurable Testbed (CoNNeCT) project utilized force-limited testing and analysis following the semi-empirical approach. This presentation presents the steps in performing a force-limited analysis and then compares the results to test data recovered during the CoNNeCT force-limited random vibration qualification test that took place at NASA Glenn Research Center (GRC) in the Structural Dynamics Laboratory (SDL) December 19, 2010 - January 7, 2011. A compilation of lessons learned and considerations for future force-limited tests is also included.

Staab, Lucas↗
NASA NTRS2012

Application of the Semi-Empirical Force-Limiting Approach for the CoNNeCT SCAN Testbed

The semi-empirical force-limiting vibration method was developed and implemented for payload testing to limit the structural impedance mismatch (high force) that occurs during shaker vibration testing. The method has since been extended for use in analytical models. The Space Communications and Navigation Testbed (SCAN Testbed), known at NASA as, the Communications, Navigation, and Networking re-Configurable Testbed (CoNNeCT), project utilized force-limiting testing and analysis following the semi-empirical approach. This paper presents the steps in performing a force-limiting analysis and then compares the results to test data recovered during the CoNNeCT force-limiting random vibration qualification test that took place at NASA Glenn Research Center (GRC) in the Structural Dynamics Laboratory (SDL) December 19, 2010 to January 7, 2011. A compilation of lessons learned and considerations for future force-limiting tests is also included.

Staab, Lucas D.↗
DOE OSTI2021

An Experimental Analysis on the CNS Simulator Comparing Human Performance between Operators and Students

This study, as a previous study of the SHEEP project, aims to identify differences in human performance data through experiments with students and operators using CNS (Compact Nuclear Simulator). To achieve this goal, this experiment seeks to analyze the correlation of measured human performance between students and operators through a randomized factorial experiment design. Five human performances are measured. In addition, several scenarios and related procedures have been developed for the CNS experiments.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2020

South Carolina Water Resources Project Summary - Implementing the Unvegetated-Vegetated Ratio to Assess Salt Marsh Vulnerability in South Carolina Using Airborne and Space-Based Remote Sensing Imagery

Among the most productive ecosystems on earth, salt marshes provide crucial ecosystem services including water filtration, shoreline protection, storm surge buffering, and flood mitigation. Marshes are largely dependent on their sediment budget which can significantly vary across a region. Upstream land use change near Charleston, South Carolina, along with rising sea levels, are expected to alter sediment budgets and threaten marsh stability and long-term health. The unvegetated-vegetated ratio (UVVR) is a scalable and efficient method to assess vulnerability. This NASA DEVELOP project collaborated with the South Carolina Department of Natural Resources, the South Carolina Department of Health and Environmental Control, and the United States Geological Survey Woods Hole Coastal and Marine Science Center. Marsh vulnerability was analyzed using UVVR derived from Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper (ETM+) in conjunction with National Agriculture Imagery Program (NAIP) high-resolution aerial imagery. A Landsat random forest regression showed a low correlation (r2 = 0.247) between Landsat 7 ETM+ bands and NAIP aggregated UVVR suggesting the need for a more complex model and higher resolution sensors. Google Earth Engine scripting provided a novel approach to UVVR methodology that will allow decision makers to input new marsh areas and easily calculate UVVR without external data downloading.

DEVELOP Project Summary↗
DOE OSTI2021

Review: Magnesium Sheet Alloy Development for Room Temperature Forming

Sheet metal forming operations in the automotive industry, including stamping, flanging, bending, hemming and trimming, are dominantly done at room temperature (RT). Unfortunately, the poor RT formability of magnesium due to its hexagonal close packed (HCP) structure and generally strong texture has limited the use of these processes in high-volume automotive production. However, the formability of magnesium can be improved via fine grain structure and random texture to enable some RT forming operations. Here, this paper presents the latest magnesium alloy development and evaluation by the United States Automotive Materials Partnership (USAMP) in collaboration with its university partners. A new sheet alloy developed in a recent USAMP project, ZAXME11100 (USAMP Alloy 2 Plus), offers excellent ductility (31% tensile elongation) and RT formability (7.8 mm Erichsen Index) in solution-treated condition (T4), and a high yield strength (270 MPa) upon post-forming aging treatment (T6), promising RT forming for automotive applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluating the Interplay between Trajectory Segmentation and Mode Inference Error

Travel behavior changes are essential to transportation decarbonization. Travel diaries, consisting of sequences of trips between places, are typically used to instrument human travel behavior. However, these diaries are only as accurate as the underlying methods used to construct them. Travel diary algorithms have been a popular research topic since the advent of Global Positioning System tracking surveys. These algorithms have typically been validated using prompted recall of presegmented trips, thus disregarding the continuity of mode inference. Phone operating systems have adopted battery-conserving techniques, but the resulting data collection errors have not been studied extensively. We introduce a framework to evaluate the accuracy of trip length computations and mode inference by analyzing continuous mode-segmented trajectories for groups of trips. We then use the framework to identify the input data quality and the impact of postprocessing. Our primary inputs to this evaluation are MobilityNet, a public dataset containing information from three artificial timelines covering 15 different travel modes, and sample open-source travel diary creation algorithms from the OpenPATH project. Our framework concretely shows that the variance of the distance error drops from (0.217, 0.0848) to (o.011, 0.0407) (Android, iOS) after postprocessing. Similarly, the weighted F-scores for mode inference increase from (0.25, 0.29) to (0.60, 0.74) (iOS, Android) between random forest and geographic information system-based models. We hope that this standardized method will be adapted to evaluate other, potentially proprietary, travel diary algorithms. Finally, the results can be used to understand and improve the state of the art in the travel diary creation field.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Accelerating Neutrino Event Generation in MARLEY Using CUDA-Based RNG and GPU Parallelization

MARLEY is a simulation tool that helps scientists study how low-energy neutrinos interact with matter. To work properly, MARLEY uses random numbers thousands of times in each simulation. These random numbers are important for modeling things like how neutrinos collide with atoms and what particles they produce. Right now, MARLEY runs on a regular computer processor (CPU) and uses a built-in random number generator called the Mersenne Twister. This setup works, but it can be slow, especially when trying to simulate many events. This research focuses on making MARLEY run faster by moving the random number generation and some of the repetitive calculations from the CPU to a graphics processing unit (GPU), which can handle many tasks at the same time. We use CUDA (a tool for programming NVIDIA GPUs) and cuRAND (a GPU-based random number library) to test faster alternatives to the current random number system. We compare different GPU-based generators, like curand_mtgp32, xorwow, and philox, to see which ones are the quickest and still give reliable results. Early tests show that using the GPU can make MARLEY simulations much faster. This project not only helps improve current simulation performance but also moves closer to a full simulation chain where all stages can run on modern GPU hardware.

Dunkley, Kimieka [Florida A-M]↗
NASA NTRS1986

The effect of a row of Helmholtz resonators on the turbulence in a boundary layer

Detailed acoustic measurements were made of the resonator's response as the free stream speed was varied from 30 mph to 155 mph. The average sound pressure level (SPL) and peak frequency response clearly show the existence of strong tuning between the boundary layer and the resonator for the Helmholtz (f sub 0) mode as well as for the first standing wave (f sub 1) mode. The narrow speed range for tuning and the gap between the strong tuning for the Helmholtz mode and the first standing wave mode is evident. The frequency at which the peak SPL response occurred at each speed is given. The peak response for the Helmholtz mode occurred at a free stream velocity of 26 m/s (at which RE sub theta = 6,560) with a resonant frequency of 570 Hz and a sound pressure level of 141 dB. The peak response for the first standing wave occurs at the maximum wind tunnel speed of 70 m/s (at which Re sub theta = 14,900) Hz with a resonant frequency of 1,890 Hz and a sound pressure level of 154 dB. At resonance the microphones do not maintain a constant SPL, but rather the latter fluctuates in a seemingly random fashion. The values presented here are time averaged rms values. Of interest was the occurrence of phase locking between adjacent resonators witha phase lag of 180 deg. This report covers the period from 5 May to 30 September on the subject project.

Panton, R. L.↗
DOE OSTI2023

2D k -th nearest neighbour statistics: a highly informative probe of galaxy clustering

ABSTRACT Beyond standard summary statistics are necessary to summarize the rich information on non-linear scales in the era of precision galaxy clustering measurements. For the first time, we introduce the 2D k-th nearest neighbour (kNN) statistics as a summary statistic for discrete galaxy fields. This is a direct generalization of the standard 1D kNN by disentangling the projected galaxy distribution from the redshift-space distortion signature along the line-of-sight. We further introduce two different flavours of 2D kNNs that trace different aspects of the galaxy field: the standard flavour which tabulates the distances between galaxies and random query points, and a ‘DD’ flavour that tabulates the distances between galaxies and galaxies. We showcase the 2D kNNs’ strong constraining power both through theoretical arguments and by testing on realistic galaxy mocks. Theoretically, we show that 2D kNNs are computationally efficient and directly generate other statistics such as the popular two-point correlation function (2PCF), voids probability function, and counts-in-cell statistics. In a more practical test, we apply the 2D kNN statistics to simulated galaxy mocks that fold in a large range of observational realism and recover parameters of the underlying extended halo occupation distribution (HOD) model that includes velocity bias and galaxy assembly bias. We find unbiased and significantly tighter constraints on all aspects of the HOD model with the 2D kNNs, both compared to the standard 1D kNN, and the classical redshift-space 2PCF.

79 ASTRONOMY AND ASTROPHYSICS↗