Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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 631 records · Page 35

Photoionization cross sections for atomic chlorine using an open-shell random phase approximation

The use of the Random Phase Approximation with Exchange (RPAE) for calculating partial and total photoionization cross sections and photoelectron angular distributions for open shell atoms is examined for atomic chlorine. Whereas the RPAE corrections in argon (Z=18) are large, it is found that those in chlorine (Z=17) are much smaller due to geometric factors. Hartree-Fock calculations with and without core relaxation are also presented. Sizable deviations from the close coupling results of Conneely are also found.

Starace, A. F.↗

Yield stress of cemented tungsten carbide

Cemented tungsten carbide yields plastically at room temperature in the presence of a large hydrostatic pressure component. By approximate analysis of the state of stress in supported opposed anvils and by measurement of the pressure at which the anvil tips exhibit a permanent deviation from planarity, we have obtained the yield stress of such materials. Our value for the yield stress of a 3% cobalt cemented tungsten carbide is 86 kbar.

Ruoff, A. L.↗

Statistical analysis of the 70 meter antenna surface distortions

Statistical analysis of surface distortions of the 70 meter NASA/JPL antenna, located at Goldstone, was performed. The purpose of this analysis is to verify whether deviations due to gravity loading can be treated as quasi-random variables with normal distribution. Histograms of the RF pathlength error distribution for several antenna elevation positions were generated. The results indicate that the deviations from the ideal antenna surface are not normally distributed. The observed density distribution for all antenna elevation angles is taller and narrower than the normal density, which results in large positive values of kurtosis and a significant amount of skewness. The skewness of the distribution changes from positive to negative as the antenna elevation changes from zenith to horizon.

Kiedron, K.↗

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗

Application of Commercial Non-Dispersive Infrared Spectroscopy Sensors for Sub-Ambient Carbon Dioxide Detection

Monitoring carbon dioxide (CO2) concentration within a spacecraft or spacesuit is critically important to ensuring the safety of the crew. Carbon dioxide uniquely absorbs light at wavelengths of 3.95 micrometers and 4.26 micrometers. As a result, non-dispersive infrared (NDIR) spectroscopy can be employed as a reliable and inexpensive method for the quantification of CO2 within the atmosphere. A multitude of commercial off-the-shelf (COTS) NDIR sensors exist for CO2 quantification. The COTS sensors provide reasonable accuracy as long as the measurements are attained under conditions close to the calibration conditions of the sensor (typically 21.1 C (70.0 F) and 1 atmosphere). However, as pressure deviates from atmospheric to the pressures associated with a spacecraft (8.0{10.2 pounds per square inch absolute (psia)) or spacesuit (4.1{8.0 psia), the error in the measurement grows increasingly large. In addition to pressure and temperature dependencies, the infrared transmissivity through a volume of gas also depends on the composition of the gas. As the composition is not known a priori, accurate sub-ambient detection must rely on iterative sensor compensation techniques. This manuscript describes the development of recursive compensation algorithms for sub-ambient detection of CO2 with COTS NDIR sensors. In addition, the source of the exponential loss in accuracy is developed theoretically. The basis of the loss can be explained through thermal, Doppler, and Lorentz broadening effects that arise as a result of the temperature, pressure, and composition of the gas mixture under analysis. This manuscript provides an approach to employing COTS sensors at sub-ambient conditions and may also lend insight into designing future NDIR sensors for aerospace application.

Swickrath, Michael J.↗

Objective Motion Cueing Criteria Investigation Based on Three Flight Tasks

This paper intends to help establish fidelity criteria to accompany the simulator motion system diagnostic test specified by the International Civil Aviation Organization. Twelve air- line transport pilots flew three tasks in the NASA Vertical Motion Simulator under four different motion conditions. The experiment used three different hexapod motion configurations, each with a different tradeoff between motion filter gain and break frequency, and one large motion configuration that utilized as much of the simulator's motion space as possible. The motion condition significantly affected: 1) pilot motion fidelity ratings, and sink rate and lateral deviation at touchdown for the approach and landing task, 2) pilot motion fidelity ratings, roll deviations, maximum pitch rate, and number of stick shaker activations in the stall task, and 3) heading deviation after an engine failure in the takeoff task. Significant differences in pilot-vehicle performance were used to define initial objective motion cueing criteria boundaries. These initial fidelity boundaries show promise but need refinement.

pilot training↗

NASA NDE Fracture Critical Detectable Flaw Sizes History and Methodology

NASA requires that NDE methods and inspectors demonstrate 90% Probability of Detection (POD) with 95% confidence for critical flaw sizes when inspecting fracture critical metallic components. NASA addresses the known variability of NDE inspector capability in two ways. The first, Special NDE, requires that every inspector demonstrate the required 90/95 POD, which is resource intensive. The second approach is Standard NDE for which conservative flaw sizes for different NDE methods are provided such that it is expected that most properly trained inspectors will exceed the 90/95 POD requirement. As such, individual POD demonstration testing is not required. The origin of NASA Standard NDE dates to the start of the Space Shuttle Program in the early 1970’s. In the first study to quantitatively assess NDE methods and inspectors, the performance of multiple inspectors was evaluated for different NDE methods using a large set of fatigue cracked specimens. A rudimentary POD analysis was performed to estimate the 90/95 POD flaw size for each inspector for each method. Additionally, the average and standard deviation of the 90/95 POD flaw size across the multiple inspectors was calculated to estimate the flaw size for which 95 percent of inspectors would provide the 90/95 POD capability. These estimated 90/95/95 POD flaw sizes evolved into the NASA Standard NDE flaw sizes still in use for structural analysis five decades later. The methodology for performing Standard NDE POD studies was never documented in NASA requirements. Furthermore, POD analysis methods have significantly evolved since this seminal study. Likewise, NDE methods have improved and there has been a push to reassess Standard NDE flaw sizes for existing methods, and to develop Standard NDE flaw sizes for new methods such as digital radiography. In this study, a Standard NDE POD methodology was developed and baselined using the historical data. This reanalysis of the historical data identified several deficiencies in the original test plan as well as an overall lack of conservatism in the estimated 90/95/95 POD flaw sizes. The results of this historical review and the new methodology are being incorporated into an update of NASA NDE POD requirements.

Nondestructive Evaluation↗

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗

Production of Large-Particle-Size Monodisperse Latexes in Microgravity

A latex is a suspension of very tiny (micrometer-size) plastic spheres in water, stabilized by emulsifiers. The growth of billions of these tiny plastic spheres to sizes larger than can be grown on Earth is attempted while keeping all of them exactly the same size and perfectly spherical. Thus far on several of the Monodisperse Latex Reactor (MLR) flights, the latex spheres have been returned to Earth with standard deviations of better than 1.4%. In microgravity the absence of buoyancy effects has allowed growth of the balls up to 30 micrometers in diameter thus far. The MLR has now flown 5 times on the Shuttle. The MLR has now produced the first commercial space product; that is the first commercial material ever manufactured in space and marketed on Earth. Once it is demonstrated that these large-size-monodisperse latexes can be routinely produced in quantity and quality, they can be marketed for many types of scientific applications. They can be used in biomedical research for such things as drug carriers and tracers in the body, human and animal blood flow studies, membrane and pore-sizing in the body, and medical diagnostic tests.

Vanderhoff, J. W.↗

Surface diagnosis of large reflector antennas using microwave holographic metrology - An iterative approach

The radiation pattern of an antenna can be significantly affected by reflector surface errors. A minimization of losses due to surface inaccuracies is, therefore, required. A suitable procedure for the determination of the deviations of the surface from its ideal geometry is based on the microwave holographic technique. In this procedure, the complex pattern of the antenna is measured, and the Fourier transform relationship, existing between the radiation pattern and the function related to the current distribution on the reflector surface, is applied. It is naturally desirable to limit the amount of measured data. The present investigation is concerned with the concept of the Fourier transform relationship between the induced current and the far field. It is demonstrated that the relationship, in general, can be interpreted as a summation of many Fourier transforms. An iterative scheme is employed to continue the far-field pattern outside the measurement window.

Rahmat-Samii, Y.↗

Typical and Unusual Properties of Magnetic Clouds during the WIND Era

A list of 33 magnetic clouds as identified in WIND magnetic field and plasma data has been compiled. The intervals for these events are provided as part of NASA/GSFC, WIND-MFI's Website under the URL http://lepmfi.qsfc.nasa.gov/mfi/mag_cloud publ.html#table The period covered in this study is from early 1995 to November 1998 which primarily occurs in the quiet part of the solar cycle. A force free, cylindrically symmetric, magnetic field model has been applied to the field data in 1-hour averaged form for all of these events (except one small event where 10 min avg's were used) and the resulting fit-parameters examined. Each event was provided a semi-quantitatively determined quality factor (excellent, good or poor). A set of 28 good or better cases, spanning a surprisingly large range of values for its various properties, was used for further analysis. These properties are, for example, durations, attitudes, sizes, asymmetries, axial field strengths, speeds, and relative impact parameters. They will be displayed and analyzed, along with some related derived quantities, with emphasis on typical vs unusual properties and on the magnetic fields magnetic clouds' relationships to the Sun and to upstream interplanetary shocks, where possible. For example, it is remarkable how narrowly distributed the speeds of these clouds are, and the overall average speed (390 techniques km/s) is less than that normally quoted for the average solar wind speed (420 km/s) despite the fact that many of these clouds are d"drivers" of interplanetary shocks. On average, a cloud appears to be a little less symmetric when the spacecraft is able to pass close to the cloud's axis as compared to a farther out passage. The average longitude and latitude (in GSE) of the axes of the clouds are 85 degrees and 8 degrees, respectively, with standard deviations near 40 degrees. Also, the half=yearly averaged axial magnetic flux has approximately tripled. almost monotonically, from about 6 to 17 X 10(exp 29) Mx over the first 3.5 years of consideration, but with a large uncertainty on each of the half-year estimates, because of small sampling. If true,this finding implies an approximate tripling of the events' solar fluxes over this period as it goes into solar maximum.

Lepping, R. P.↗

Investigation of Non-Equilibrium Radiation for Earth Entry

For Earth re-entry at velocities between 8 and 11.5 km/s, the accuracy of NASA's computational uid dynamic and radiative simulations of non-equilibrium shock layer radiation is assessed through comparisons with measurements. These measurements were obtained in the NASA Ames Research Center's Electric Arc Shock Tube (EAST) facility. The experiments were aimed at measuring the spatially and spectrally resolved radiance at relevant entry conditions for both an approximate Earth atmosphere (79% N2 : 21% O2 by mole) as well as a more accurate composition featuring the trace species Ar and CO2 (78.08% N2 : 20.95% O2 : 0.04% CO2 : 0.93% Ar by mole). The experiments were configured to target a wide range of conditions, of which shots from 8 to 11.5 km/s at 0.2 Torr (26.7 Pa) are examined in this paper. The non-equilibrium component was chosen to be the focus of this study as it can account for a significant percentage of the emitted radiation for Earth re-entry, and more importantly, non-equilibrium has traditionally been assigned a large uncertainty for vehicle design. The main goals of this study are to present the shock tube data in the form of a non-equilibrium metric, evaluate the level of agreement between the experiment and simulations, identify key discrepancies and to examine critical aspects of modeling non-equilibrium radiating flows. Radiance pro les integrated over discreet wavelength regions, ranging from the Vacuum Ultra Violet (VUV) through to the Near Infra-Red (NIR), were compared in order to maximize both the spectral coverage and the number of experiments that could be used in the analysis. A previously defined non-equilibrium metric has been used to allow comparisons with several shots and reveal trends in the data. Overall, LAURA/HARA is shown to under-predict EAST by as much as 40% and over-predict by as much as 12% depending on the shock speed. DPLR/NEQAIR is shown to under-predict EAST by as much as 50% and over-predict by as much as 20% depending on the shock speed. The one standard deviation scatter in the EAST results was calculated to be 31%. An estimate for the upper bound of the absolute error in wall-directed heat flux was calculated. Below 9 km/s, where the relative difference is large, the absolute error in radiative heat flux due to non-equilibrium models is estimated to be less then 1 W/sq cm. At the highest shock speed of 11 km/s, the error in non-equilibrium is estimated to be less than 20 W/sq cm.

Shock Tube↗

A Multi-Wavelength Investigation of the Radio-Loud Supernova PTF11qcj and its Circumstellar Environment

We present the discovery, classification, and extensive panchromatic (from radio to X-ray) follow-up observations of PTF11qcj, a supernova (SN) discovered by the Palomar Transient Factory (PTF). Our observations with the Karl G. Jansky Very Large Array show that this event is radio-loud: PTF11qcj reached a radio peak luminosity comparable to that of the famous gamma-ray-burst-associated SN 1998bw (L (sub 5 gigahertz) approximates 10 (sup 29) ergs per second per hertz). PTF11qcj is also detected in X-rays with the Chandra Observatory, and in the infrared band with Spitzer. Our multi-wavelength analysis probes the SN interaction with circumstellar material. The radio observations suggest a progenitor mass-loss rate of approximately 10 (sup -4) the mass of the sun per year times (v (sub omega) divided by 1000 kilometers per second), and a velocity of approximately 0.3-0.5 c (the speed of light in a vacuum) for the fastest moving ejecta (at around 10 days after explosion). However, these estimates are derived assuming the simplest model of SN ejecta interacting with a smooth circumstellar wind, and do not account for possible inhomogeneities in the medium and asphericity of the explosion. The radio data show deviations from such a simple model, as well as a late-time re-brightening. The X-ray flux from PTF11qcj is compatible with the high-frequency extrapolation of the radio synchrotron emission (within the large uncertainties). A light echo from pre-existing dust is in agreement with our infrared data. Our pre-explosion data from the PTF suggest that a precursor eruption of absolute magnitude M (sub r) approximately equal to -13 mag may have occurred approximately 2.5 years prior to the SN explosion. Overall, PTF11qcj fits the expectations from the explosion of a Wolf-Rayet star. Precursor eruptions may be a feature characterizing the final pre-explosion evolution of such stars.

Radio↗

Transformations of Ti-5Al-5V-5Cr-3Mo powder due to reuse in laser powder bed fusion: A surface analytical approach

Reactive Ti-5Al-5 V-5Cr-3Mo (Ti5553) alloy feedstock powder is used in laser powder-bed fusion (LPBF). Due to the large quantity of powder necessary for LPBF, powder is reused numerous times under oxidizing environments. Transformations to the powder’s native surface oxide may impact the LPBF process and lead to deviations from expected behavior of printed parts. In this work, we present a multimodal characterization of new and reused powder by combining X-ray photoelectron spectroscopy (XPS), time-of-flight secondary ion mass spectrometry (ToF-SIMS), and transmission electron microscopy (TEM) of focused ion beam (FIB)-prepared cross sections of individual powder particles to understand how the surface oxide composition changes, grows, and becomes hydroxylated and hydrated due to reuse. We show that the native oxide film of Ti5553 powder is a hydroxylated mixed oxide composed of TiO 2 and Al 2 O 3 that grows from 5.6 +/- 0.7 nm to 8.3 +/- 1.1 nm due to reuse. Al oxide regions of the surface oxide become more hydroxylated during reuse in comparison to Ti oxide regions. This indicates that an outer oxide of pure TiO 2 may enhance recyclability of Ti5553 powder while the presence of Al oxides may decrease recyclability. Further characterization of surface transformations due to recyclability may eliminate unnecessary variation in LPBF.

36 MATERIALS SCIENCE↗

Hydrogen effects on the deformation and slip localization in a single crystal austenitic stainless steel

Hydrogen is known to embrittle austenitic stainless steels, which are widely used in high-pressure hydrogen storage and delivery systems, but the mechanisms that lead to such material degradation are still being elucidated. The current work investigates the deformation behavior of single crystal austenitic stainless steel 316L through combined uniaxial tensile testing, characterization and atomistic simulations. Thermally precharged hydrogen is shown to increase the critical resolved shear stress (CRSS) without previously reported deviations from Schmid’s law. Molecular dynamics simulations further expose the statistical nature of the hydrogen and vacancy contributions to the CRSS in the presence of alloying. Slip distribution quantification over large in-plane distances (> 1 mm), achieved via atomic force microscopy (AFM), highlights the role of hydrogen increasing the degree of slip localization in both single and multiple slip configurations. The most active slip bands accumulate significantly more deformation in hydrogen precharged specimens, with potential implications for damage nucleation. For $\langle$110$\rangle$ tensile loading, slip localization further enhances the activity of secondary slip, increases the density of geometrically necessary dislocations and leads to a distinct lattice rotation behavior compared to hydrogen-free specimens, as evidenced by electron backscatter diffraction (EBSD) maps. Finally, the results of this study provide a more comprehensive picture of the deformation aspect of hydrogen embrittlement in austenitic stainless steels.

36 MATERIALS SCIENCE↗

Observation of electroweak W + W – pair production in association with two jets in proton-proton collisions at √s = 13 TeV

An observation is reported of the electroweak production of a W + W – pair in association with two jets, with both W bosons decaying leptonically. The data sample corresponds to an integrated luminosity of 138 fb –1 of proton-proton collisions at √s = 13 TeV, collected by the CMS detector at the CERN LHC. Events are selected by requiring exactly two opposite-sign leptons (electrons or muons) and two jets with large pseudorapidity separation and high dijet invariant mass. Events are categorized based on the flavor of the final-state leptons. A signal is observed with a significance of 5.6 standard deviations (5.2 expected) with respect to the background-only hypothesis. The measured fiducial cross section is 10.2 ± 2.0 fb and this value is consistent with the standard model prediction of 9.1 ± 0.6 fb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

BSE versus StarTrack: Implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N-body evolution program NBODY7. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. All the above newly implemented schemes are also shared with the standalone versions of SSE and BSE. Furthermore, all these demonstrations are performed with both the updated standalone BSE and the updated NBODY7/BSE. When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing. Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

QUIJOTE scientific results – IV. A northern sky survey in intensity and polarization at 10–20 GHz with the multifrequency instrument

ABSTRACT We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17, and 19 GHz, and covering approximately 29 000 deg2, including most of the northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE multifrequency instrument (MFI), and have angular resolutions of around 1°, and sensitivities in polarization within the range 35–40 µK per 1° beam, being a factor ∼2–4 worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterize the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5 per cent. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A, and 3C274), which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11 GHz and WMAP 23 GHz map peaks at β = −3.09 with a standard deviation of 0.14. The measured BB/EE ratio at scales of ℓ = 80 is 0.26 ± 0.07 for a Galactic cut |b| > 10°. We find a positive TE correlation for 11 GHz at large angular scales (ℓ ≲ 50), while the EB and TB signals are consistent with zero in the multipole range 30 ≲ ℓ ≲ 150. The maps discussed in this paper are publicly available.

Astronomy & Astrophysics↗