Search NASA⌕ Search

SEARCH · Search NASA

Results for “Taylor approximation”

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

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

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

An energy-stable time-integrator for phase-field models

Here, we introduce a provably energy-stable time-integration method for general classes of phase-field models with polynomial potentials. We demonstrate how Taylor series expansions of the nonlinear terms present in the partial differential equations of these models can lead to expressions that guarantee energy-stability implicitly, which are second-order accurate in time. The spatial discretization relies on a mixed finite element formulation and isogeometric analysis. We also propose an adaptive time-stepping discretization that relies on a first-order backward approximation to give an error-estimator. This error estimator is accurate, robust, and does not require the computation of extra solutions to estimate the error. This methodology can be applied to any second-order accurate time-integration scheme. We present numerical examples in two and three spatial dimensions, which confirm the stability and robustness of the method. The implementation of the numerical schemes is done in PetIGA, a high-performance isogeometric analysis framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Dynamic field theory and equations of motion in cosmology

We discuss a field-theoretical approach based on general-relativistic variational principle to derive the covariant field equations and hydrodynamic equations of motion of baryonic matter governed by cosmological perturbations of dark matter and dark energy. The action depends on the gravitational and matter Lagrangian. The gravitational Lagrangian depends on the metric tensor and its first and second derivatives. The matter Lagrangian includes dark matter, dark energy and the ordinary baryonic matter which plays the role of a bare perturbation. The total Lagrangian is expanded in an asymptotic Taylor series around the background cosmological manifold defined as a solution of Einstein’s equations in the form of the Friedmann–Lemaître–Robertson–Walker (FLRW) metric tensor. The small parameter of the decomposition is the magnitude of the metric tensor perturbation. Each term of the series expansion is gauge-invariant and all of them together form a basis for the successive post-Friedmannian approximations around the background metric. The approximation scheme is covariant and the asymptotic nature of the Lagrangian decomposition does not require the post-Friedmannian perturbations to be small though computationally it works the most effectively when the perturbed metric is close enough to the background FLRW metric. The temporal evolution of the background metric is governed by dark matter and dark energy and we associate the large scale inhomogeneities in these two components as those generated by the primordial cosmological perturbations with an effective matter density contrast δρ/ρ≤1. The small scale inhomogeneities are generated by the condensations of baryonic matter considered as the bare perturbations of the background manifold that admits δρ/ρ≫1. Mathematically, the large scale perturbations are given by the homogeneous solution of the linearized field equations while the small scale perturbations are described by a particular solution of these equations with the bare stress–energy tensor of the baryonic matter. We explicitly work out the covariant field equations of the successive post-Friedmannian approximations of Einstein’s equations in cosmology and derive equations of motion of large and small scale inhomogeneities of dark matter and dark energy. We apply these equations to derive the post-Friedmannian equations of motion of baryonic matter comprising stars, galaxies and their clusters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dynamically Downscaled (WRF) 1km, Hourly Meteorological Conditions 1987-2020. East/Taylor Watersheds

This dataset contains meteorological output from the Weather Research and Forecasting (WRF) version 3.8.1. This dataset has been created to 1) investigate hydrometeorological processes impacting the East River and water-delivery to the Critical Zone, and 2) provide meteorological forcing data for distributed Earth-science modeling applications in the East River watershed. Variables have a 1 kilometer spatial resolution and hourly temporal resolution and encompass a rectangular region encompassing the East and Taylor River watersheds, Colorado, near the town of Crested Butte. WRF was forced using Climate Forecast System Reanalysis (CFSR) lateral boundary conditions. Each .zip file contains one "water year" of data (October 1 -- September 30; i.e. water year 2017 starts October 1, 2016 and ends September 30, 2017). Each zip folder contains 12 netcdf (.nc) files containing one month of hourly data each and are approximately 250mb. Model timestamps are in UTC time.The files contain the following data variables:EAST_MASK: binary mask of the watershed regionTAYLOR_MASK: binary mask of the watershed regionGLW downwelling longwave radiation (w/m2)HR_PRCP: Hourly Precipitation Rate (mm/hr). Includes all hydrometeors (solid+liquid). HFX NoahMP LSM total grid-cell modelled sensible heat flux (w/m2) [positive towards atmosphere]LH NoahMP LSM total grid-cell modelled latent heat flux (w/m2) [positive towards atmosphere; can be converted to ET]PSFC Surface Barometric Pressure (hPa)Q2 Two-meter specific humidity (kg/kg)SWDOWN Downwelling shortwave solar radiation (w/m2)SWNORM Terrain-normal downwelling shortwave radiation (w/m2)T2 Two-meter air temperature (deg K)U10 10-m U-component of wind velocity (m/s)V10 10m V-component of wind velocity (m/s)XLAT Latitude of grid-center point XLONG Longitude of grid-center point XTIME Model timestamp, **in UTC**

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Using Hessian update formulae to construct modified Shepard interpolated potential energy surfaces: Application to vibrating surface atoms

Modified Shepard interpolation based on second order Taylor series expansions has proven to be a flexible tool for constructing potential energy surfaces in a range of situations. Extending this to gas–surface dynamics where surface atoms are allowed to move represents a substantial increase in the dimensionality of the problem, reflected in a dramatic increase in the computational cost of the required Hessian (matrix of second derivatives) evaluations. This work demonstrates that using approximate Hessians derived from well known Hessian update formulae and a single accurate Hessian can provide an effective way to avoid this expensive accurate Hessian determination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Gradient Flow and Scale Setting on MILC HISQ Ensembles

We report on a scale determination with gradient-flow techniques on the N f = 2 + 1 + 1 HISQ ensembles generated by the MILC collaboration. The ensembles include four lattice spacings, ranging from approximately 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales p √t 0 /a and w 0 /a and their tree-level improvements,√t 0;imp and w 0;imp , are computed on each ensemble using Symanzik ow and the cloverleaf definition of the energy density E. Using a combination of continuum chiral perturbation theory and a Taylor-series ansatz for the lattice-spacing and strong-coupling dependence, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. We also determine the scales p t 0 = 0:1416( +8 -5 ) fm and w 0 = 0:1717( +12 -11 ) fm, where the errors are sums, in quadrature, of statistical and all systematic errors. The precision of w 0 and √t 0 is comparable to or more precise than the best previous estimates, respectively. We also find the continuum mass-dependence of w 0 that will be useful for estimating the scales of other ensembles. Furthermore, we estimate the integrated autocorrelation length of . For long flow times, the autocorrelation length of appears to be comparable to or smaller than that of the topological charge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1973

Comparison of Te and Ti from Ogo 6 and from various incoherent scatter radars.

Langmuir probe and retarding potential analyzer (RPA) data on the electron and ion temperatures Te and Ti obtained from Ogo 6 are compared with Te and Ti values obtained from the incoherent scatter network. The satellite to radar temperature ratio TeS/TeR is 1.15 on the average for these comparisons. This discrepancy is larger than the uncertainties usually placed on the probe and radar Te values. The ion temperature ratio TiS/TiR approximately 1.0, independent of the particular radar examined. This comparison serves as an intercalibration of the incoherent scatter network.

Mcclure, J. P.↗
NASA NTRS2011

Scale Dependence of Magnetic Helicity in the Solar Wind

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first pioneered by Matthaeus et al. by which, under the assumption of spatial homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hypothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the symmetric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change of magnetic helicity at wavenumber k approximately equal to 2AU(sup -1) (or frequency nu approximately equal to 2 microHz) at distances below 2.8AU and at k approximately equal to 30AU(sup -1) (or nu approximately equal to 25 microHz) at larger distances. At small scales the magnetic helicity is positive at northern heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is argued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux integrated separately over one hemisphere amounts to about 10(sup 45) Mx(sup 2) cycle(sup -1) at large scales and to a three times lower value at smaller scales.

Brandenburg, Axel↗
NASA NTRS1993

The granulite suite: Impact melts and metamorphic breccias of the early lunar crust

The granulite suite consists of two major types of rocks. One is coarse-grained and poikilitic with many euhedral crystals of olivine and plagioclase. These characteristics indicate crystallization from a melt; the poikilitic granulites are impact melt breccias. The other group is finer-grained and granoblastic, with numerous triple junctions; the granoblastic granulites are metamorphic rocks. Compositional groups identified by Lindstrom and Lindstrom contain both textural types. Two pyroxene thermometry indicates that both groups equilibrated at 1000 to 1150 C. Calculations suggest that the granoblastic group, which has an average grain size of about 80 microns, was annealed for less than 6 x 10 exp 4 y at 1000 C, and for less than 2500 y at 1150 C. Similar equilibration temperatures suggest that both groups were physically associated after impact events produced the poikilitic melts. Granulitic impactites hold important information about the pre-Nectarian bombardment history of the Moon, and the composition and thermal evolution of the early lunar crust. Granulitic impactites are widely considered to be an important rock type in the lunar crust, but how they formed is poorly understood. Metal compositions and elevated concentrations of meteoritic siderophile elements suggest that most lunar granulites are impact breccias. Their occurrence as clasts in approximately 3.9 Ga breccias, and Ar-(40-39) ages greater than or = 4.2 Ga for some granulites show that they represent a component of the lunar crust which formed prior to the Nectarian cataclysm. Petrographic characteristics of lunar granulites indicate at least two endmember textural variants which apparently formed in fundamentally different ways. One type has granoblastic textures consisting of equant, polygonal to rounded grains, and abundant triple junctions with small dispersions around 120 degrees indicating a close approach to textural equilibrium. As suggested by many authors, granoblastic granulites probably formed by subsolidus annealing and recrystallization of fragmental or glassy protoliths. Examples of this type include 15418, 78155, and 79215. The other textural type consists of poikilitic to poikiloblastic rocks with euhedral to subhedral plagioclase and olivine enclosed by interstitial pyroxene. In some cases, the texture resembles that of an orthocumulate. Examples of this type include 60035, 67955, and 77017. Rounding of grain edges is common in poikilitic granulites, but the regular crystal shapes and widely dispersed dihedral angles show they are far from textural equilibrium. The textures of poikilitic granulites are more consistent with the formation of these rocks by crystallization from a melt than by subsolidus metamorphism. A few samples have been recognized with textural characteristics transitional between those of the granoblastic and poikiloblastic endmembers (e.g., 72559, 78527). Pyroxene compositions taken from the literature and determined for this study by electron microprobe were used to calculate equilibration temperatures. The Kretz Ca transfer (solvus) thermometer and the Lindsley and Anderson graphical method both give similar temperatures, which range from approximately 1000 to 1150 C. There is no apparent temperature difference between granoblastic and poikilitic varieties, but there is a hint in these data that the more ferroan varieties equilibrated to lower temperatures. Additional studies are in progress to test this possibility.

Cushing, J. A.↗
DOE OSTI2014

Semiconductor neutron detectors using depleted uranium oxide

This paper reports on recent attempts to develop and test a new type of solid-state neutron detector fabricated from uranium compounds. It has been known for many years that uranium oxide (UO2), triuranium octoxide (U3O8) and other uranium compounds exhibit semiconducting characteristics with a broad range of electrical properties. We seek to exploit these characteristics to make a direct-conversion semiconductor neutron detector. In such a device a neutron interacts with a uranium nucleus, inducing ssion. The ssion products deposit energy-producing, detectable electron-hole pairs. The high energy released in the ssion reaction indicates that noise discrimination in such a device has the potential to be excellent. Schottky devices were fabricated using a chemical deposition coating technique to deposit UO2 layers a few microns thick on a sapphire substrate. Schottky devices have also been made using a single crystal from UO2 samples approximately 500 microns thick. Neutron sensitivity simulations have been performed using GEANT4. Neutron sensitivity for the Schottky devices was tested experimentally using a 252Cf source.

neutron detection, uranium oxide, semiconductor↗
NASA NTRS1972

Survival of bacterial isolates exposed to simulated Jovian trapped radiation belt electrons and solar wind protons

With missions to Jupiter, the spacecraft will be exposed for extended duration to solar wind radiation and the Jovian trapped radiation belt. This study is designed to determine the effect of these radiation environments on spacecraft bacterial isolates. The information can be used in the probability of contamination analysis for these missions. A bacterial subpopulation from Mariner Mars 1971 spacecraft (nine sporeforming and three nonsporeforming isolates) plus two comparative organisms, Staphylococcus epidermidis ATCC 17917 and a strain of Bacillus subtilis var. niger, were exposed to 2-, 12-, and 25-MeV electrons at different doses with simultaneous exposure to a vacuum of 0.0013 N/sqm at 20 and -20 C. The radioresistance of the subpopulation was dependent on the isolate, dose, and energy of electrons. Temperature affected the radioresistance of only the sporeforming isolates. Survival data indicated that spores were reduced approximately 1 log/1500 J/kg, while nonsporeforming isolates (micrococci) were reduced 1.5 to 2 logs/1500 J/kg with the exception of an apparent radioresistant isolate whose resistance approached that of the spores. The subpopulation was found to be less resistant to lower energy than to higher energy electrons.

Taylor, D. M.↗
DOE OSTI2021

Guide field effects on the distribution of plasmoids in multiple scale reconnection

The effects of a finite guide field on the distribution of plasmoids in high-Lundquist-number current sheets undergoing magnetic reconnection in large plasmas are investigated with statistical models. Merging of plasmoids is taken into account either assuming that guide field flux is conserved resulting in nonforce-free profiles in general, or that magnetic helicity is conserved and Taylor relaxation occurs to convert part of the summed guide field flux into reconnecting field flux toward minimum energy states resulting in force-free profiles. It is found that the plasmoid distribution in terms of reconnecting field flux follows a power law with index 7/4 or 1 depending on whether merger frequencies are independent of or dependent on their relative velocity to the outflow speed, respectively. This result is approximately the same for the force-free and nonforce-free models, with nonforce-free models exhibiting indices of 2 and 1 for the same velocity dependencies. Distributions in terms of guide field flux yield indices of 3/2 for the nonforce-free model regardless of velocity dependence. This is notably distinct from the indices of 11/8 and 1 for the force-free models independent of and dependent on velocity, respectively. At low guide field fluxes, the force-free models exhibit a second power law index of 1/2 due to nonconstant flux growth rates. The velocity-dependent force-free model predicts the production of slightly more rapidly moving large guide field flux plasmoids which are supported by observational evidence of flux ropes with strong core fields. Implications are discussed on particle acceleration via Fermi processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Theory for Glassy Behavior of Supercooled Liquid Mixtures

Here, we introduce a model for glassy dynamics in supercooled liquid mixtures. Given the relaxation behavior of individual supercooled liquids, the model predicts the relaxation times of their mixtures as temperature is decreased. The model is based on dynamical facilitation theory for glassy dynamics, which provides a physical basis for relaxation and vitrification of a supercooled liquid. This is in contrast to empirical linear interpolations such as the Gordon-Taylor equation typically used to predict glass transition temperatures of liquid mixtures. To understand the behavior of supercooled liquid mixtures we consider a multicomponent variant of the kinetically constrained East model in which components have a different energy scale and can also diffuse when locally mobile regions, i.e., excitations, are present. Using a variational approach we confirm an effective single component model with a single effective energy scale that best approximates a mixture. When scaled by this single effective energy, we show that experimental relaxation times of many liquid mixtures all collapse onto the "parabolic law" predicted by dynamical facilitation theory. The model can be used to predict transport properties and glass transition temperatures of mixtures of glassy materials, with implications in atmospheric chemistry, biology, and pharmaceuticals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Evaluating the grain-scale deformation behavior of a single-phase FCC high entropy alloy using synchrotron high energy diffraction microscopy

We report although the deformation behavior of high-entropy alloys (HEAs) has been extensively studied at the macroscale, many important properties have yet to be explored for these alloys at the microscale, thus hampering accurate prediction of damage and failure processes. Synchrotron high-energy diffraction microscopy (HEDM) and fast-Fourier transform-based crystal plasticity modeling was conducted to investigate the three-dimensional (3D) grain-resolved micromechanical response for approximately 1,900 constituent grains within a single-phase FCC HEA up to 1% applied strain. The evolution of grain-resolved elastic strains, lattice reorientations, and maximum resolved shear stresses (mRSS) were evaluated to quantify elastic, yield, and fully plastic behavior. Overall, the initial critical resolved shear stress (CRSS), determined via in situ HEDM and companion modeling, was found to be > 20% higher than estimated using the classical polycrystalline Taylor factor (M = 3.06). However, a descriptive parameter based on the average grain-resolved Taylor factor ($\overline{M}$) was found to show excellent agreement with plastic yielding behavior observed within HEDM datasets. Noticeable deviations in HEDM lattice reorientations compared to both EVP-FFT simulations and classical predictions for FCC polycrystals were discovered, highlighting the complexity in correlating local lattice reorientations, Taylor, and Schmid factors with plastic response for this material at the grain-scale. Therefore, it is anticipated that the overall trends and parameter identification of 3D grain-resolved properties in this study can serve as an important foundation for continued mesoscale investigation on both well-established and newly developed Cantor-like HEAs.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Accelerated panel methods using the fast multipole method

Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Leathrum, James F., Jr.↗
DOE OSTI2014

An enhanced droplet-based liquid microjunction surface sampling system coupled with HPLC-ESI-MS/MS for spatially resolved analysis

Droplet-based liquid microjunction surface sampling coupled with high-performance liquid chromatography (HPLC)-electrospray ionization (ESI)-tandem mass spectrometry (MS/MS) for spatially resolved analysis provides the possibility of effective analysis of complex matrix samples and can provide a greater degree of chemical information from a single spot sample than is typically possible with a direct analysis of an extract. Described here is the setup and enhanced capabilities of a discrete droplet liquid microjunction surface sampling system employing a commercially available CTC PAL autosampler. The system enhancements include incorporation of a laser distance sensor enabling unattended analysis of samples and sample locations of dramatically disparate height as well as reliably dispensing just 0.5 μL of extraction solvent to make the liquid junction to the surface, wherein the extraction spot size was confined to an area about 0.7 mm in diameter; software modifications improving the spatial resolution of sampling spot selection from 1.0 to 0.1 mm; use of an open bed tray system to accommodate samples as large as whole-body rat thin tissue sections; and custom sample/solvent holders that shorten sampling time to approximately 1 min per sample. Lastly, the merit of these new features was demonstrated by spatially resolved sampling, HPLC separation, and mass spectral detection of pharmaceuticals and metabolites from whole-body rat thin tissue sections and razor blade (“crude”) cut mouse tissue.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Triple and quadruple excitation contributions to the binding in Be clusters: Calibration calculations on Be3

The contribution of connected triple and quadruple excitations to the binding in Be3 is investigated by comparing various coupled-cluster (CC) and truncated configuration interaction (CI) treatments with multireference CI (MRCI) and full CI(FCI) calculations. The CC method with single and double excitations (CCSD) produces results that differ substantially from more elaborate treatments, but most extensions to CCSD that account approximately for connected triple excitations perform very well. In constrast, good agreement with FCI for Be2 can be achieved only with the highest level CC and MRCI methods.

Watts, John D.↗
NASA NTRS1987

Polarization aberrations in optical systems

This paper examines the polarization aberrations in optical systems of isotropic radially symmetric weak polarizers in the paraxial approximation. Optical systems which fit this description are composed of uncoated and coated lenses and mirrors. Polarization aberrations are variations in phase, amplitude, and polarization state of the electromagnetic field across the exit pupil. Some are dependent on the incident polarization state and some are not. Expressions through fourth order for phase, amplitude, linear polarization, and linear retardance aberrations are derived in terms of the chief and marginal ray angles of incidence and Taylor series expansion coefficients of the Fresnel equations for reflection and transmission at interfaces.

Mcguire, James P., Jr.↗