Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dilatancy”

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 253 records · Page 14

Adsorption space for microporous polymers with diverse adsorbate species

Abstract The enormous number of combinations of adsorbing molecules and porous materials that exist is known as adsorption space. The adsorption space for microporous polymers has not yet been systematically explored, especially when compared with efforts for crystalline adsorbents. We report molecular simulation data for the adsorptive and structural properties of polymers of intrinsic microporosity with a diverse set of adsorbate species with 345 distinct adsorption isotherms and over 240,000 fresh and swollen structures. These structures and isotherms were obtained using a sorption-relaxation technique that accounts for the critical role of flexibility of the polymeric adsorbents. This enables us to introduce a set of correlations that can estimate adsorbent swelling and fractional free volume dilation as a function of adsorbate uptake based on readily characterized properties. The separation selectivity of the 276 distinct binary molecular pairs in our data is reported and high-performing adsorbent systems are identified.

36 MATERIALS SCIENCE↗

Deep learning of structural morphology imaged by scanning X-ray diffraction microscopy

Scanning X-ray nanodiffraction microscopy is a powerful technique for spatially resolving nanoscale structural morphologies by diffraction contrast. One of the critical challenges in experimental nanodiffraction data analysis is posed by the convergence angle of nanoscale focusing optics which creates simultaneous dependency of the far-field scattering data on three independent components of the local strain tensor-corresponding to dilation and two potential rigid body rotations of the unit cell. All three components are in principle resolvable through a spatially mapped sample tilt series; however, traditional data analysis is computationally expensive and prone to artifacts. In this study, we implement NanobeamNN, a convolutional neural network specifically tailored to the analysis of scanning probe X-ray microscopy data. NanobeamNN learns lattice strain and rotation angles from simulated diffraction of a focused X-ray nanobeam by an epitaxial thin film and can directly make reasonable predictions on experimental data without the need for additional fine-tuning. We demonstrate that this approach represents a significant advancement in computational speed over conventional methods, as well as a potential improvement in accuracy over the current standard.

Luo, Aileen [Cornell Univ., Ithaca, NY (United Sta↗

Structure of iso-density sets in supersonic isothermal turbulence

Context: The gas density structure of the cold molecular phase of the interstellar medium is the main controller of star formation. Aims: A theoretical framework is proposed to describe the structural content of the density field in isothermal supersonic turbulence. Methods: It makes use of correlation and structure functions of the phase indicator field defined for different iso-density values. The relations between these two-point statistics and the geometrical features of iso-density sets such as the volume fraction, the surface density, the curvature, and fractal characteristics are provided. An exact scale-by-scale budget equation is further derived revealing the role of the turbulent cascade and dilation on the structural evolution of the density field. Although applicable to many flow situations, this tool is here first invoked for characterising supersonic isothermal turbulence, using data from the currently best-resolved numerical simulation. Results: We show that iso-density sets are surface fractals rather than mass fractals, with dimensions that markedly differ between dilute, neutral, and dense regions. The surface–size relation is established for different iso-density values. We further find that the turbulent cascade of iso-density sets is directed from large towards smaller scales, in agreement with the classical picture that turbulence acts to concentrate more surface into smaller volumes. Intriguingly, there is no range of scales that complies with a constant transfer rate in the cascade, challenging our fundamental understanding of interstellar turbulence. Finally, we recast the virial theorem in a new formulation drawing an explicit relation between the aforementioned geometrical measures and the dynamics of iso-density sets.

79 ASTRONOMY AND ASTROPHYSICS↗

Dam-Break Flow Measurements of Liquid–Granular Mixtures in a Channelized Reservoir

Dam-break flows of granular materials exhibit a complex non-Newtonian rheology, which is strongly influenced by various factors, such as the characteristics of the granular material, the degree of initial compaction, and saturation level of the granular matrix. Specially adapted nonintrusive techniques are required to investigate these flows in the laboratory. This paper investigated the nonintrusive technique of image analysis. The methodologies used in the measurement of bulk flow parameters, such as flow depth, front velocity, and flow rate, and the extraction of detailed flow characteristics such as velocity profiles, were discussed in the context of a series of granular dam-break flow experiments. These experiments were performed under three types of initial conditions: dry, saturated (water filling the pore space), and submerged (a layer of saturated granular material submerged under water) using four types of granular materials—PET pellets, crushed walnut shells, urea, and sand—with grain sizes ranging from 200 μm to 2.9 mm. High-speed cameras were used to capture the evolution and local properties of the flow, such as velocity fields using particle image velocimetry (PIV). Bulk properties of the flow were obtained by tracking the surface profile of the granular matrix and the water table. The extracted profiles were used to measure the evolution of the volume of the granular matrix and the fluid phase separately upstream and downstream of the gate, which allowed measuring the dilatancy (increase in pore space). The study showed that the image analysis techniques are versatile and can be adapted to measure different flow parameters in a variety of granular flows to better understand and analyze their complex behavior. Finally, the advantages and drawbacks of each image analysis technique were discussed and best-practice recommendations were provided.

42 ENGINEERING↗

Biomimetic composite architecture achieves ultrahigh rate capability and cycling life of sodium ion battery cathodes

Sodium ion batteries are an emerging candidate to replace lithium ion batteries in large-scale electrical energy storage systems due to the abundance and widespread distribution of sodium. Despite the growing interest, the development of high-performance sodium cathode materials remains a challenge. In particular, polyanionic compounds are considered as a strong cathode candidate owing to their better cycling stability, a flatter voltage profile, and stronger thermal stability compared to other cathode materials. Here, we report the rational design of a biomimetic bone-inspired polyanionic Na3V2(PO4)3-reduced graphene oxide composite (BI-NVP) cathode that achieves ultrahigh rate charging and ultralong cycling life in a sodium ion battery. At a charging rate of 1 C, BI-NVP delivers 97% of its theoretical capacity and is able to retain a voltage plateau even at the ultra-high rate of 200 C. It also shows long cycling life with capacity retention of 91% after 10 000 cycles at 50 C. The sodium ion battery cells with a BI-NVP cathode and Na metal anode were able to deliver a maximum specific energy of 350 W h kg−1 and maximum specific power of 154 kW kg−1. In situ and postmortem analyses of cycled BI-NVP (including by Raman and XRD spectra) HRTEM, and STEM-EELS, indicate highly reversible dilation–contraction, negligible electrode pulverization, and a stable NVP-reduced graphene oxide layer interface. The results presented here provide a rational and biomimetic material design for the electrode architecture for ultrahigh power and ultralong cyclability of the sodium ion battery full cells when paired with a sodium metal anode.

25 ENERGY STORAGE↗

Characterization of the hardened single line of sight camera at the National Ignition Facility

The hardened single line of sight camera has been recently characterized in preparation for its deployment on the National Ignition Facility. The latest creation based on the pulse-dilation technology leads to many new features and improvements over the previous-generation cameras to provide better quality measurements of inertial confinement fusion experiments, including during high neutron yield implosions. Here, we present the characterization data that illustrate the main performance features of this instrument, such as extended dynamic range and adjustable internal magnification, leading to improved spatial resolution.

Trosseille, C. (ORCID:0000000289430482)↗

Design of the high-yield time-gated x-ray hot-spot imager for OMEGA

Time-resolved x-ray self-emission imaging of hot spots in inertial confinement fusion experiments along several lines of sight provides critical information on the pressure and the transient morphology of the hot spot on the University of Rochester’s OMEGA Laser System. At least three quasi-orthogonal lines of sight are required to infer the tomographic information of the hot spots of deuterium–tritium cryogenic layered implosions. OMEGA currently has two time-gated x-ray hot-spot imagers: the time-resolved Kirkpatrick–Baez x-ray microscope and the single-line-of-sight, time-resolved x-ray imager (SLOS-TRXI). The time-gated x-ray hot-spot imager (XRHSI) is being developed for use on OMEGA as the third line of sight for the high-yield operation of up to 4 × 10 14 neutrons. XRHSI follows the SLOS-TRXI concept; however, it will have improved spatial and temporal resolutions of 5 μm and 20 ps, respectively. The simultaneous operation of the three instruments will provide 3-D reconstructions of the assembled hot-spot fuel at various times through peak thermonuclear output. The technical approach consists of a pinhole array imager and demagnifying time-dilation drift tube that are coupled to two side-by-side hybrid complementary metal–oxide semiconductor (hCMOS) image sensors. To minimize the background and to harden the diagnostics, an angled drift-tube assembly shifting the hCMOS sensors out of the direct line of sight and neutron shielding will be applied. As a result, the technical design space for the instrument will be discussed and the conceptual design will be presented.

47 OTHER INSTRUMENTATION↗

Elucidating ion capture and transport mechanisms of Preyssler anions in aqueous solutions using biased MACE-accelerated MD simulations

Equilibrium and biased multi-atomic cluster expansion (MACE) accelerated molecular dynamics (MD) simulations in aqueous solutions are performed to investigate the ion capture and transport mechanisms of the {P 5 W 30 } Preyssler anion (PA) as the smallest representative member of the extended polyoxometalate (POM) family with an internal cavity. The unique interatomic interactions present in the internal cavity vs the exterior of PA are carefully investigated using equilibrium MACE MD simulations for two representative Na(H 2 O)@PA and Na@PA complexes in aqueous solutions. Our careful analyses of radial distribution functions and coordination numbers show that the presence of confined water in Na(H 2 O)@PA has profound modulating effects on the nature of the interactions of the encapsulated ion with the oxygens of the PA cavity. Using well-converged MACE-accelerated multiple walker well-tempered metadynamics simulations with nanosecond timescales, two different associative (ion exchange) and dissociative (ion ejection) ion transport mechanisms were carefully investigated for Na + as one of the most abundant and representative ions present in seawater and saline solutions. By comparing systems with and without confined water, it was found that the presence of only one pre-encapsulated confined water in Na(H 2 O)@PA dramatically changes the free energy landscape of ion transport processes. It was also found that the contraction and dilation of the two windows present in PA directly influence the Na + and H 2 O transport. Furthermore, the results from this work are helpful, as they show a viable path toward tuning the ion exchange and transport phenomena in aqueous solutions of POM molecular clusters and frameworks.

25 ENERGY STORAGE↗

Molecular basis of force-pCa relation in MYL2 cardiomyopathy mice: Role of the super-relaxed state of myosin

Here, in this study, we investigated the role of the super-relaxed (SRX) state of myosin in the structure–function relationship of sarcomeres in the hearts of mouse models of cardiomyopathy-bearing mutations in the human ventricular regulatory light chain (RLC, MYL2 gene). Skinned papillary muscles from hypertrophic (HCM–D166V) and dilated (DCM–D94A) cardiomyopathy models were subjected to small-angle X-ray diffraction simultaneously with isometric force measurements to obtain the interfilament lattice spacing and equatorial intensity ratios (I 11 /I 10 ) together with the force-pCa relationship over a full range of [Ca 2+ ] and at a sarcomere length of 2.1 μm. In parallel, we studied the effect of mutations on the ATP-dependent myosin energetic states. Compared with wild-type (WT) and DCM–D94A mice, HCM–D166V significantly increased the Ca 2+ sensitivity of force and left shifted the I 11 /I 10 -pCa relationship, indicating an apparent movement of HCM–D166V cross-bridges closer to actin-containing thin filaments, thereby allowing for their premature Ca 2+ activation. The HCM–D166V model also disrupted the SRX state and promoted an SRX-to-DRX (super-relaxed to disordered relaxed) transition that correlated with an HCM-linked phenotype of hypercontractility. While this dysregulation of SRX ↔ DRX equilibrium was consistent with repositioning of myosin motors closer to the thin filaments and with increased force-pCa dependence for HCM–D166V, the DCM–D94A model favored the energy-conserving SRX state, but the structure/function–pCa data were similar to WT. Our results suggest that the mutation-induced redistribution of myosin energetic states is one of the key mechanisms contributing to the development of complex clinical phenotypes associated with human HCM–D166V and DCM–D94A mutations.

59 BASIC BIOLOGICAL SCIENCES↗

Microstructure and creep strength of simulated intercritical heat-affected zone of grade 91 steel

The intercritical heat affected zone (ICHAZ) has been reported as one of the most Type IV cracking susceptible regions in 9Cr creep-resistant steel weldments. However, creep degradation mechanisms within the ICHAZ itself need further clarifications. In this work, two ICHAZ specimens of Grade 91 steel, low-temperature ICHAZ (LT-ICHAZ) and high-temperature ICHAZ (HT-ICHAZ), were simulated using the Gleeble thermomechanical system by exposing to two peak temperatures (860 °C and 900 °C) between $A_{C1}$ and $A_{C3}$. Dramatically different creep strengths of two simulated ICHAZs were observed and studied. We report dilation curve analysis indicates a high fraction of newly transformed martensite formed in the HT-ICHAZ, which results in a much higher hardness (360 HV0.5) of the HT-ICHAZ than 266 HV0.5 of LT-ICHAZ. Precipitates, especially $M_{23}C_6$ carbides, were not fully dissolved in both ICHAZ specimens. After a typical postweld heat treatment (760 °C-2 hours), the faster recovery of low-carbon martensite and reduced precipitation strengthening due to $M_{23}C_6$ carbides coarsening in the HT-ICHAZ led to a significant hardness reduction. These microstructural degradations in the HT-ICHAZ made its creep lifetime about 36 times shorter than that of the LT-ICHAZ tested at 650 °C. The remaining tempered martensite from base metal in the LT-ICHAZ was the primary contributor to maintain its high creep resistance.

36 MATERIALS SCIENCE↗

Magnetoelastic coupling and spin contributions to entropy and thermal transport in biferroic yttrium orthochromite

Direct engineering of material properties through exploitation of spin, phonon, and charge-coupled degrees of freedom is an active area of development in materials science. However, the relative contribution of the competing orders to controlling the desired behavior is challenging to decipher. In particular, the independent role of phonons, magnons, and electrons, quasiparticle coupling, and relative contributions to the phase transition free energy largely remain unexplored, especially for magnetic phase transitions. Here, we study the lattice and magnetic dynamics of biferroic yttrium orthochromite using Raman, infrared, and inelastic neutron spectroscopy techniques, supporting our experimental results with first-principles lattice dynamics and spin-wave simulations across the antiferromagnetic transition at T N ~ 138 K. Spectroscopy data and simulations together with the heat capacity (C p ) measurements, allow us to quantify individual entropic contributions from phonons (0.01 ± 0.01k B atom -1 ), dilational (0.03 ± 0.01k B atom -1 ), and magnons (0.11 ± 0.01k B atom -1 ) across TN. High-resolution phonon measurements conducted in a magnetic field show that anomalous T-dependence of phonon energies across TN originates from magnetoelastic coupling. Additionally, phonon scattering is primarily governed by the phonon–phonon coupling, with little contribution from magnon–phonon coupling, short-range spin correlations, or magnetostriction effects; a conclusion further supported by our thermal conductivity measurements conducted up to 14 T, and phenomenological modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Synchronizing the consistency relation

We study the N-point function of the density contrast to quadratic order in the squeezed limit during the matter-dominated (MD) and radiation-dominated (RD) eras in synchronous gauge. Since synchronous gauge follows the free-fall frame of observers, the equivalence principle dictates that in the gradient approximation for the long-wavelength mode there is only a single, manifestly time-independent consistency relation for the N-point function. This simple form is dictated by the initial mapping between synchronous and local coordinates, unlike Newtonian gauge and its correspondingly separate dilation and Newtonian consistency relations. Dynamical effects only appear at quadratic order in the squeezed limit and are again characterized by a change in the local background, also known as the separate universe approach. We show that for the 3-point function the compatibility between these squeezed-limit relations and second-order perturbation theory requires both the initial and dynamical contributions to match, as they do in single-field inflation. This clarifies the role of evolution or late-time projection effects in establishing the consistency relation for observable bispectra, which is especially important for radiation acoustic oscillations and for establishing consistency below the matter-radiation equality scale in the MD era. Defining an appropriate angle and time average of these oscillations is also important for making separate universe predictions of spatially varying local observables during the RD era, which can be useful for a wider range of cosmological predictions beyond N-point functions.

79 ASTRONOMY AND ASTROPHYSICS↗

Catalog-level blinding on the bispectrum for DESI-like galaxy surveys

We evaluate the performance of the catalog-level blind analysis technique (blinding) presented in Brieden et al. (2020) in the context of a fixed template power spectrum and bispectrum analysis. This blinding scheme, which is tailored for galaxy redshift surveys similar to the Dark Energy Spectroscopic Instrument (DESI), has two components: the so-called “AP blinding” (concerning the dilation parameters α$_{∥}$, α$_{⊥}$) and “RSD blinding” (redshift space distortions, affecting the growth rate parameter f). Through extensive testing, including checks for the RSD part in cubic boxes, the impact of AP blinding on mocks with realistic survey sky coverage, and the implementation of a full AP+RSD blinding pipeline, our analysis demonstrates the effectiveness of the technique in preserving the integrity of cosmological parameter estimation when the analysis includes the bispectrum statistic. We emphasize the critical role of sophisticated — and difficult to accidentally unblind — blinding methods in precision cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Fiducial-cosmology-dependent systematics for the DESI 2024 BAO analysis

When measuring the Baryon Acoustic Oscillations (BAO) scale from galaxy surveys, one typically assumes a fiducial cosmology when converting redshift measurements into comoving distances and also when defining input parameters for the reconstruction algorithm. A parameterised template for the model to be fitted is also created based on a (possibly different) fiducial cosmology. This model reliance can be considered a form of data compression, and the data is then analysed allowing that the true answer is different from the fiducial cosmology assumed. In this study, we evaluate the impact of the fiducial cosmology assumed in the BAO analysis of the Dark Energy Spectroscopic Instrument (DESI) survey Data Release 1 (DR1) on the final measurements in DESI 2024 III. We utilise a suite of mock galaxy catalogues with survey realism that mirrors the DESI DR1 tracers: the bright galaxy sample (BGS), the luminous red galaxies (LRG), the emission line galaxies (ELG) and the quasars (QSO), spanning a redshift range from 0.1 to 2.1. We compare the four secondary AbacusSummit cosmologies against DESI's fiducial cosmology (Planck 2018). The secondary cosmologies explored include a lower cold dark matter density, a thawing dark energy universe, a higher number of effective species, and a lower amplitude of matter clustering. The mocks are processed through the BAO pipeline by consistently iterating the grid, template, and reconstruction reference cosmologies. We determine a conservative systematic contribution to the error of 0.1% for both the isotropic and anisotropic dilation parameters αiso and αAP. We then directly test the impact of the fiducial cosmology on DESI DR1 data.

79 ASTRONOMY AND ASTROPHYSICS↗

Physics-inspired spatiotemporal-graph AI ensemble for the detection of higher order wave mode signals of spinning binary black hole mergers

We present a new class of AI models for the detection of quasi-circular, spinning, non-precessing binary black hole mergers whose waveforms include the higher order gravitational wave modes ($\ell$, |m|) = {(2,2), (2,1), (3,3), (3,2), (4,4)}, and mode mixing effects in the $\ell$ = 3, |m| = 2 harmonics. These AI models combine hybrid dilated convolution neural networks to accurately model both short- and long-range temporal sequential information of gravitational waves; and graph neural networks to capture spatial correlations among gravitational wave observatories to consistently describe and identify the presence of a signal in a three detector network encompassing the Advanced LIGO and Virgo detectors. We first trained these spatiotemporal-graph AI models using synthetic noise, using 1.2 million modeled waveforms to densely sample this signal manifold, within 1.7 h using 256 NVIDIA A100 GPUs in the Polaris supercomputer at the Argonne Leadership Computing Facility. This distributed training approach exhibited optimal classification performance, and strong scaling up to 512 NVIDIA A100 GPUs. With these AI ensembles we processed data from a three detector network, and found that an ensemble of 4 AI models achieves state-of-the-art performance for signal detection, and reports two misclassifications for every decade of searched data. We distributed AI inference over 128 GPUs in the Polaris supercomputer and 128 nodes in the Theta supercomputer, and completed the processing of a decade of gravitational wave data from a three detector network within 3.5 h. Finally, we fine-tuned these AI ensembles to process the entire month of February 2020, which is part of the O3b LIGO/Virgo observation run, and found 6 gravitational waves, concurrently identified in Advanced LIGO and Advanced Virgo data, and zero false positives. This analysis was completed in one hour using one NVIDIA A100 GPU.

79 ASTRONOMY AND ASTROPHYSICS↗

Observational window effects on multi-object reverberation mapping

ABSTRACT Contemporary reverberation mapping campaigns are employing wide-area photometric data and high-multiplex spectroscopy to efficiently monitor hundreds of active galactic nuclei (AGNs). However, the interaction of the window function(s) imposed by the observation cadence with the reverberation lag and AGN variability time-scales (intrinsic to each source over a range of luminosities) impact our ability to recover these fundamental physical properties. Time dilation effects due to the sample source redshift distribution introduce added complexity. We present comprehensive analysis of the implications of observational cadence, seasonal gaps, and campaign baseline duration (i.e. the survey window function) for reverberation lag recovery. We find that the presence of a significant seasonal gap dominates the efficacy of any given campaign strategy for lag recovery across the parameter space, particularly for those sources with observed-frame lags above 100 d. Using the Australian Dark Energy Survey as a baseline, we consider the implications of this analysis for the 4MOST/Time-Domain Extragalactic Survey campaign providing concurrent follow-up of the Legacy Survey of Space and Time deep-drilling fields, as well as upcoming programmes. We conclude that the success of such surveys will be critically limited by the seasonal visibility of some potential field choices, but show significant improvement from extending the baseline. Optimizing the sample selection to fit the window function will improve survey efficacy.

79 ASTRONOMY AND ASTROPHYSICS↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: cosmological implications from multitracer BAO analysis with galaxies and voids

We construct cosmic void catalogues with the dive void finder upon SDSS BOSS DR12 and eBOSS DR16 galaxy samples with BAO reconstruction, and perform a joint BAO analysis using different types of galaxies and the corresponding voids. The BAO peak is evident for the galaxy–galaxy, galaxy–void, and void–void correlation functions of all data sets, including the ones cross-correlating LRG and ELG samples. Two multitracer BAO fitting schemes are tested, one combining the galaxy and void correlation functions with a weight applied to voids, and the other using a single BAO dilation parameter for all clustering measurements. Both methods produce consistent results with mock catalogues, and on average ~10 per cent improvements of the BAO statistical uncertainties are observed for all samples, compared to the results from galaxies alone. By combining the clustering of galaxies and voids, the uncertainties of BAO measurements from the SDSS data are reduced by 5–15 per cent, yielding 0.9 per cent, 0.8 per cent, 1.1 per cent, 2.3 per cent, and 2.9 per cent constraints on the distance $D_{_{\rm V}}(z)$, at effective redshifts 0.38, 0.51, 0.70, 0.77, and 0.85, respectively. When combined with BAO measurements from SDSS MGS, QSO, and Lyα samples, as well as the BBN results, we obtain $H_0 = 67.58 \pm 0.91\, {\rm km}\, {\rm s}^{-1}\, {\rm Mpc}^{-1}$, Ω m = 0.290 ± 0.015, and $\Omega _\Lambda h^2 = 0.3241 \pm 0.0079$ in the flat-ΛCDM framework, where the 1σ uncertainties are around 6 per cent, 6 per cent, and 17 per cent smaller respectively, compared to constraints from the corresponding anisotropic BAO measurements without voids and LRG–ELG cross-correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

Mixing bispectrum multipoles under geometric distortions

We derive general expressions for how the Alcock–Paczynski distortions affect the power spectrum and the bispectrum of cosmological fields. We compute explicit formulas for the mixing coefficients of bispectrum multipoles in the linear approximation. The leading-order effect for the bispectrum is the uniform dilation of all three wavevectors. The mixing coefficients depend on the shape of the bispectrum triplet. Our results for the bispectrum multipoles are framed in terms of the ‘natural’ basis of the lengths of three wavevectors but can be easily generalized for other bases and reduction schemes. Our validation tests confirm that the linear approximation is extremely accurate for all power spectrum multipoles. The linear approximation is accurate for the bispectrum monopole but results in sub-per cent level inaccuracies for the bispectrum quadrupole and fails for the bispectrum hexadecapole. Our results can be used to simplify the analysis of the bispectrum from galaxy surveys, especially the measurement of the baryon acoustic oscillation peak position. They can be used to replace numeric schemes with exact analytical formulae.

79 ASTRONOMY AND ASTROPHYSICS↗