Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural evolution”

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

Comparative In Situ Measurements of Plasma Instabilities in the Equatorial and Auroral Electrojets

This presentation provides a comparison of in situ measurements of plasma instabilities gathered by rocket-borne probes in the equatorial and auroral electrojets. Specifically, using detailed measurements of the DC electric fields, current density, and plasma number density within the unstable daytime equatorial electrojet from Brazil (Guara Campaign) and in the auroral electrojet from Sweden (ERRIS Campaign), we present comparative observations and general conclusions regarding the observed physical properties of Farley-Buneman two-stream waves and large scale, gradient drift waves. The two stream observations reveal coherent-like waves propagating near the E x B direction but at reduced speeds (nearer to the presumed acoustic velocity) with wavelengths of approximately 5-10m in both the equatorial and auroral electrojet, as measured using the spaced-receiver technique. The auroral electrojet data generally shows extensions to shorter wavelengths, in concert with the fact that these waves are driven harder. With respect to gradient-drift driven waves, observations of this instability are much more pronounced in the equatorial electrojet, given the more favorable geometry for growth provided by the vertical gradient and horizontal magnetic field lines. We present new analysis of Guara rocket observations of electric field and plasma density data that reveal considerable structuring in the middle and lower portion of the electrojet (90-105 km) where the ambient plasma density gradient is unstable. Although the electric field amplitudes are largest (approximately 10-15 mV/m) in the zonal direction, considerable structure (approximately 5-10 mV/m) is also observed in the vertical electric field component as well, implying that the dominant large scale waves involve significant vertical interaction and coupling within the narrow altitude range where they are observed. Furthermore, a detailed examination of the phase of the waveforms show that on some, but not all occasions, locally enhanced eastward fields are associated with locally enhanced upwards (polarization) electric fields. The measurements are discussed in terms of theories involving the non-linear evolution and structuring of plasma waves.

Pfaff, Robert F.↗

Model-Based Fatigue Prognosis of Fiber-Reinforced Laminates Exhibiting Concurrent Damage Mechanisms

Prognostics of large composite structures is a topic of increasing interest in the field of structural health monitoring for aerospace, civil, and mechanical systems. Along with recent advancements in real-time structural health data acquisition and processing for damage detection and characterization, model-based stochastic methods for life prediction are showing promising results in the literature. Among various model-based approaches, particle-filtering algorithms are particularly capable in coping with uncertainties associated with the process. These include uncertainties about information on the damage extent and the inherent uncertainties of the damage propagation process. Some efforts have shown successful applications of particle filtering-based frameworks for predicting the matrix crack evolution and structural stiffness degradation caused by repetitive fatigue loads. Effects of other damage modes such as delamination, however, are not incorporated in these works. It is well established that delamination and matrix cracks not only co-exist in most laminate structures during the fatigue degradation process but also affect each other's progression. Furthermore, delamination significantly alters the stress-state in the laminates and accelerates the material degradation leading to catastrophic failure. Therefore, the work presented herein proposes a particle filtering-based framework for predicting a structure's remaining useful life with consideration of multiple co-existing damage-mechanisms. The framework uses an energy-based model from the composite modeling literature. The multiple damage-mode model has been shown to suitably estimate the energy release rate of cross-ply laminates as affected by matrix cracks and delamination modes. The model is also able to estimate the reduction in stiffness of the damaged laminate. This information is then used in the algorithms for life prediction capabilities. First, a brief summary of the energy-based damage model is provided. Then, the paper describes how the model is embedded within the prognostic framework and how the prognostics performance is assessed using observations from run-to-failure experiments

Prognostics↗

The shape of spectral lines: The importance of the far wings

Spectroscopy, the study of the interaction of radiation and matter, provides most of the information we have gleaned about the composition, structure, and evolution of the universe. As is well known, by measuring the frequencies of spectral lines in absorption or emission, one can uniquely infer the presence of atoms or molecules as well as their physical state and environment (e.g., solid or gaseous, neutral or ionized, moving or stationary, etc.). Furthermore, by studying the intensities of these lines, one can determine the abundance (i.e., number of a particular species per unit volume). Although less well known, the shape of the spectral lines, in particular, the structure of the far wings, plays a very important role in many important atmospheric phenomena such as the greenhouse effect or the absorption of harmful ultraviolet radiation. Although first measured more than 50 years ago, the anomalous absorption of radiation by water vapor in the earth's atmosphere was postulated to be due to far wings of allowed lines. However, only within the past few years has a quantitative verification of this hypothesis been possible through the development of an accurate theoretical description of the shape of self-broadened water lines. During the summer, work has been done on improving this theory and in comparing the results to other theories valid near the center of the lines. The relevance of this work to measurements of greenhouse gases, of earth-based measurements of the 3 K cosmic background radiation, of satellite-based measurements of the atmospheres of the earth and other planets, and other similar problems will be discussed briefly.

Tipping, Richard↗

Report of the Workshop on Geologic Applications of Remote Sensing to the Study of Sedimentary Basins

The Workshop on Geologic Applications of Remote Sensing to the Study of Sedimentary Basins, held January 10 to 11, 1985 in Lakewood, Colorado, involved 43 geologists from industry, government, and academia. Disciplines represented ranged from vertebrate paleontology to geophysical modeling of continents. Deliberations focused on geologic problems related to the formation, stratigraphy, structure, and evolution of foreland basins in general, and to the Wind River/Bighorn Basin area of Wyoming in particular. Geological problems in the Wind River/Bighorn basin area that should be studied using state-of-the-art remote sensing methods were identified. These include: (1) establishing the stratigraphic sequence and mapping, correlating, and analyzing lithofacies of basin-filling strata in order to refine the chronology of basin sedimentation, and (2) mapping volcanic units, fracture patterns in basement rocks, and Tertiary-Holocene landforms in searches for surface manifestations of concealed structures in order to refine models of basin tectonics. Conventional geologic, topographic, geophysical, and borehole data should be utilized in these studies. Remote sensing methods developed in the Wind River/Bighorn Basin area should be applied in other basins.

Lang, H. R.↗

Martian Polar Vortices: Comparison of Reanalyses

The structure and evolution of the Martian polar vortices is examined using two recently available reanalysis systems: version 1.0 of the Mars Analysis Correction Data Assimilation (MACDA) and a preliminary version of the Ensemble Mars Atmosphere Reanalysis System (EMARS). There is quantitative agreement between the reanalyses in the lower atmosphere, where Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) data are assimilated, but there are differences at higher altitudes reflecting differences in the free-running general circulation model simulations used in the two reanalyses. The reanalyses show similar potential vorticity (PV) structure of the vortices: There is near-uniform small PV equatorward of the core of the westerly jet, steep meridional PV gradients on the polar side of the jet core, and a maximum of PV located off of the pole. In maps of 30 sol mean PV, there is a near-continuous elliptical ring of high PV with roughly constant shape and longitudinal orientation from fall to spring. However, the shape and orientation of the vortex varies on daily time scales, and there is not a continuous ring of PV but rather a series of smaller scale coherent regions of high PV. The PV structure of the Martian polar vortices is, as has been reported before, very different from that of Earth's stratospheric polar vortices, but there are similarities with Earth's tropospheric vortices which also occur at the edge of the Hadley Cell, and have near-uniform small PV equatorward of the jet, and a large increase of PV poleward of the jet due to increased stratification.

Waugh, D. W.↗

Synoptic Mapping of Chemical Composition, Thermal Structure, and Air Motion from UARS Observations

An operational mapping algorithm was developed to process measurements of individual species observed by different satellite instruments on board UARS. Based on Fast Fourier Synoptic Mapping (FFSM), the algorithm accounts for the precessing orbit of UARS, the error of individual instruments, and gaps associated with instrument duty cycle and the satellite yaw maneuver. It provides synoptic structure and evolution on periods as short as 1 day, derived collectively from all observations of an individual species. The algorithm was applied to synoptically map temperature, thickness, and several chemical constituents observed by the instruments: MLS, CLAES, ISAMS, and HALOE. Each field variable observed by these instruments was mapped twice-daily in continuous global time series several months long. Mapped behavior produced via FFSM was compared against standard archived products generated via Kalman filtering. The standard map products reveal similar behavior, but are limited to features of larger scale. Tracer structure that develops through flow deformation and attending transport is therefore represented in those products more coarsely, eventually being sheared down to scales that are no longer represented properly. The synoptic time series also reveal a diurnal cycle for several of the constituents, one that emerges clearly in their space-time spectra. The zonal and meridional structure of diurnal variability was mapped by filtering the space-time spectrum, an intermediate product of FFSM, to those scales resolved by the UARS sampling. Geographical variations of the diurnal cycle have also been evaluated in its seasonal-mean structure by compositing distributions at individual local times. Both were compared to diurnal variations in chemical models and in dynamical models that account for transport by the diurnal tide. Diurnal variations were found to introduce spurious behavior into the archived products generated via Kalman filtering, behavior that is misrepresented in those synoptic time series. Aliasing by under-sampled diurnal variations introduces a bias into the time-mean structure of archived products. That systematic error is manifested as anomalous wave structure in the time-mean fields. By contrast, synoptic time series generated via FFSM faithfully represent such behavior as varying diurnally, so it does not corrupt time-mean structure. Tracer structure mapped via the above algorithm was being applied to determine the Brewer-Dobson circulation of the stratosphere through Lagrangian considerations. Preliminary results are encouraging: They reveal diabatic ascent in the tropics, balanced by diabatic descent over the winter hemisphere. This Lagrangian description of transport derived from UARS tracer measurements would be potentially valuable for diagnosing interannual changes of transport, specifically, in relation to interannual changes of total ozone. Regrettably, support was reduced to 50% of the approved funding level, leaving this component incomplete.

Salby, Murry L.↗

Cometary rays - Magnetically channeled outflow

The present model of the structure and evolution of cometary rays views this radiation as the channeled outflow of cometary ions from the ionosphere. While ions in the central tail ray flow down the channel formed by the tailward closing of the draped interplanetary field, ions in the side rays flow through magnetic flux ropes that become embedded in the cometary ionosphere through Kelvin-Helmholtz instability of the cometary ionopause. The flux ropes tunnel through the ionosphere, penetrate the ionopause at two points, and extend out into the solar wind on both sides of the comet. Differences in observed cometary ray structure are attributed to differences in the radial distance and ecliptic latitude of different comets at the time of observation.

Wolff, R. S.↗

Photoevaporation of disks around massive stars and application to ultracompact H II regions

Young massive stars produce sufficient Lyman continuum photon luminosity Phi(sub i) to significantly affect the structure and evolution of the accretion disks surrounding them. A nearly static, ionized, isothermal 10(exp 4) K atmosphere forms above the neutral disk for disk radii r less than r(sub g) = 10(exp 15) M(sub 1) cm, where M(sub *) = 10 solar mass M(sub 1) is the stellar mass. For r approximately greater than r(sub g) the diffuse field created by hydrogen recombinations to the ground state in the photoionized gas above the disk produces a steady evaporation at the surface of the disk, and this H II gas flows freely out to the ISM (the 'disk wind'). The detailed structure depends on the mass-loss rate dot-M(sub w) of the fast, approximately greater than 1000 km/sec, stellar wind from the massive star. A critical mass-loss rate dot-M(sub cr) is defined such that the ram pressure of the stellar wind equals the thermal pressure of the H II atmosphere at r(sub g). In the weak stellar wind solution, dot-M(sub w) less than dot-M(sub cr), the diffuse photons from the atmosphere above r(sub g) produce a photoevaporative mass-loss rate from the disk at r approximately greater than r(sub g) of order 1 x 10(exp -5)(Phi(sub 49))(exp 1/2)(M(sub 1))(exp 1/2) solar mass/year, where Phi(sub i) = 10(exp 49) Phi(sub 49)/sec. The resulting slow (10 to 50 km/sec) ionized outflow, which persists for approximately greater than 10(exp 5) year for disk masses M(sub d) approximately 0.3 M(sub *), may explain the observational characteri stics of unresolved, ultracompact H II regions. In the strong stellar wind solution, dot-M(sub w) greater than dot-M(sub cr), the ram pressure of the stellar wind blows down the atmosphere for r less than r(sub g) and allows the stellar photons to penetrate to greater radii and smaller heights. A slow, ionized outflow produced mainly by diffuse photons is again created for r less than r(sub g); however, it is now dominated by the flow at r(sub w)(greater than r(sub g)), the radius at which the stellar wind ram pressure equals the thermal pressure in the evaporating flow. The mass-loss rate from the disk is of order 6 x 10(exp -5)dot-M(sub w-6) v(sub w8)(Phi (sub 49))(exp -1/2) solar mass/year, where dot-M(sub w-6) = M(sub w)/10(exp -6) solar mass/year and v(sub w8) = v(sub w)/1000 km/sec is the stellar wind velocity. The resulting outflow, which also persists for approximately greater than 10(exp 5) year may explain many of the more extended (r approximately greater than 10(exp 16) cm) ultracompact H II regions. Both the weak-wind and the strong-wind models depend entirely on stellar parameters Phi(sub i), M(sub *), dot-M(sub w)) and are independent of disk parameters as long as an extended r much greater than (r(sub g)), neutral disk exists. We compare both weak-wind and strong-wind model results to the observed radio free-free spectra and luminosities of ultracompact H II regions and to the interesting source MWC 349.

Hollenbach, David↗

Perception of interstellar structure - Facing complexity

Some orthodox notions concerning the structure and evolution of star-forming regions are challenged; it is proposed that they arise largely by a dual process in which conceptual models are fashioned after categories which are in great part reflections of observational limitations, and the models are projected onto interpretations of data. Several examples are discussed. The need for internal support of molecular clouds is questioned. It is suggested that the inverse density-size relation often claimed for clouds and accounted for by several theoretical models is an artifact caused by limited dynamic range column density detectability, selection bias, distance uncertainties, and internal density gradients, and is contradicted by several unbiased surveys. Column density structures mapped with a large spatial and column density dynamic range are found to be dominated by irregular, connected, and nested forms on all sides. These features and a technique for the quantification of the complex structure are illustrated with a densely-sampled column density image of the Taurus region constructed from IRAS data.

Scalo, John↗

Equatorial Density Irregularity Structures at Intermediate Scales and Their Temporal Evolution

We examine high resolution measurements of ion density in the equatorial ionosphere from the AE-E satellite during the years 1977-1981. Structure over spatial scales from 18 km to 200 m is characterized by the spectrum of irregularities at larger and smaller scales and at altitudes above 350 km and below 300 km. In the low-altitude region, only small amplitude large-scale (lambda greater than 5 km) density modulations are often observed, and thus the power spectrum of these density structures exhibits a steep spectral slope at kilometer scales. In the high-altitude region, sinusoidal density fluctuations, characterized by enhanced power near 1-km scale, are frequently observed during 2000-0200 LT. However, such fluctuations are confined to regions at the edges of larger bubble structures where the average background density is high. Small amplitude irregularity structures, observed at early local time hours, grow rapidly to high-intensity structures in about 90 min. Fully developed structures, which are observed at late local time hours, decay very slowly producing only-small differences in spectral characteristics even 4 hours later. The local time evolution of irregularity structure is investigated by using average statistics for low-(1% less than sigma less than 5%) and high-intensity (sigma greater than 10%) structures. At lower altitudes, little chance in the spectral slope is seen as a function of local time, while at higher attitudes the growth and maintenance of structures near 1 km scales dramatically affects the spectral slope.

Kil, Hyosub↗

Empirical Modeling of the Storm Time Innermost Magnetosphere Using Van Allen Probes and THEMIS Data: Eastward and Banana Currents

The structure of storm time currents in the inner magnetosphere, including its innermost region inside 4R(sub E), is studied for the first time using a modification of the empirical geomagnetic field model TS07D and new data from Van Allen Probes and Time History of Events and Macroscale Interactions during Substorms missions. It is shown that the model, which uses basis-function expansions instead of ad hoc current modules to approximate the magnetic field, consistently improves its resolution and magnetic field reconstruction with the increase of the number of basis functions and resolves the spatial structure and evolution of the innermost eastward current. This includes a connection between the westward ring current flowing largely at R > or approx. 3R(sub E) and the eastward ring current concentrated at R < or approx. 3R(sub E) resulting in a vortex current pattern. A similar pattern coined 'banana current' was previously inferred from the pressure distributions based on the energetic neutral atom imaging and first-principles ring current simulations. The morphology of the equatorial currents is dependent on storm phase. During the main phase, it is complex, with several asymmetries forming banana currents. Near SYM-H minimum, the banana current is strongest, is localized in the evening-midnight sector, and is more structured compared to the main phase. It then weakens during the recovery phase resulting in the equatorial currents to become mostly azimuthally symmetric.

Stephens, G. K.↗

Galaxy evolution and large-scale structure in the far-infrared. I - IRAS pointed observations

Redshifts for 66 galaxies were obtained from a sample of 93 60-micron sources detected serendipitously in 22 IRAS deep pointed observations, covering a total area of 18.4 sq deg. The flux density limit of this survey is 150 mJy, 4 times fainter than the IRAS Point Source Catalog (PSC). The luminosity function is similar in shape with those previously published for samples selected from the PSC, with a median redshift of 0.048 for the fainter sample, but shifted to higher space densities. There is evidence that some of the excess number counts in the deeper sample can be explained in terms of a large-scale density enhancement beyond the Pavo-Indus supercluster. In addition, the faintest counts in the new sample confirm the result of Hacking et al. (1989) that faint IRAS 60-micron source counts lie significantly in excess of an extrapolation of the PSC counts assuming no luminosity or density evolution.

Lonsdale, Carol J.↗

Evolution of basin and range structure in the Ruby Mountains and vicinity, Nevada

Results from various age dating techniques, seismic reflection profiling hydrocarbon maturation studies, and structural analysis were used to evaluate the Cenozoic deformation in the Ruby Mountains and adjoining ranges (pinyon Range and Cortez Range) in Elko and Eureka Counties, Nevada. Age dating techniques used include potassium-argon ages of biotites from granites published by Kistler et al. (1981) and fission track ages from apatite and zircon. Fission track ages from apatite reflect a closing temperature of 100 plus or minus 20 deg C. Zircon fission track ages reflect a closing temperature of 175 plus or minus 25 deg C and potassium-argon ages from brotite reflect a closing temperature of 250 plus or minus 30 deg C. Thus these results allow a reasonably precise tracking of the evolution of the ranges during the Cenozoic. Seismic reflection data are available from Huntington Valley. Access to seismic reflection data directly to the west of the Harrison Pass Pluton in the central Ruby Mountains was obtained. In addition results are available from several deep exploration holes in Huntington Valley.

Blackwell, D. D.↗

Early Spectral Evolution of Classical Novae: Consistent Evidence for Multiple Distinct Outflows

The physical mechanism driving mass ejection during a nova eruption is still poorly understood. Possibilities include ejection in a single ballistic event, a common-envelope interaction, a continuous wind, or some combination of these processes. Here, we present a study of 12 Galactic novae, for which we have premaximum high-resolution spectroscopy. All 12 novae show the same spectral evolution. Before optical peak, they show a slow P Cygni component. After peak, a fast component quickly arises, while the slow absorption remains superimposed on top of it, implying the presence of at least two physically distinct flows. For novae with high-cadence monitoring, a third, intermediate-velocity component is also observed. These observations are consistent with a scenario where the slow component is associated with the initial ejection of the accreted material and the fast component with a radiation-driven wind from the white dwarf. When these flows interact, the slow flow is swept up by the fast flow, producing the intermediate component. These colliding flows may produce theγ-ray emission observed in some novae. Our spectra also show that the transient heavy-element absorption lines seen in some novae have the same velocity structure and evolution as the other lines in the spectrum, implying an association with the nova ejecta rather than a preexisting circumbinary reservoir of gas or material ablated from the secondary. While this basic scenario appears to qualitatively reproduce multiwavelength observations of classical novae, substantial theoretical and observational work is still needed to untangle the rich diversity of nova properties.

E. Aydi↗

The Response of Disks to Oscillatory Modes in Galaxies

Recent studies suggest that galaxies can oscillate in normal modes with essentially no damping over a Hubble time. These modes may play an important role in the structure and evolution of disk/halo systems. Motivated by the possibility that normal mode oscillations exist in real galaxies, we are investigating the response of galactic disks to halo oscillations. The goal of these investigations is to ascertain whether or not observational signatures exist for such oscillations. Our approach is to perform numerical experiments on the response of a self-gravitating disk to a time-varying halo potential. We assume that a significant fraction of the mass in a galaxy is in a dark halo. The halo oscillates and the luminous disk material responds to these oscillations. Preliminary results are reported for disks embedded in a radially oscillating gravitational potential. The equilibrium initial disk is represented by an exponential density profile. Considerable care was taken to build an initial disk model that was "stable" over long time scales. A control experiment was run with the disk in a static halo potential. The disk responds to the time-varying potential by developing a ring structure, which forms and disappears during each halo oscillation cycle. The density of stars becomes depressed in an annular region at the radius where the disk epicycle frequency is equal to the halo oscillation period. This pattern of response persists over time periods approaching a Hubble time. In the oscillating potential, a bar develops in the inner disk. This bar is absent when the halo remains static. Specific targets of this study include the implications for large-scale disk structure, the gas dynamical response of the interstellar medium in such systems, and the inflow of material into the central regions of the galaxy.

Smith, Bruce F.↗

Recent Advances in Structural Carbon Nanotube Composites

Carbon nanotubes (CNT) have been a subject of investigation for more than two decades. Part of the interest in CNTs was driven by CNT tensile mechanical properties at the nanoscale that suggest their potential to outperform state-of-the-art (SOA) lightweight structural materials. Early work on the development of CNT nanocomposites involved the dispersion of CNT powder into polymer matrices. Although low levels of CNT doping resulted in meaningful increases in electrical conductivity to enable applications such as electrostatic charge dissipation, the highest concentrations of CNTs in these lightly doped nanocomposites did not yield mechanical properties that could compete with SOA carbon fiber composites. The maximum concentration of CNTs in lightly doped nanocomposites was limited to about 5% due to increases in viscosity that prevent the fabrication of good quality nanocomposites. At the lower concentrations, properties of the resulting nanocomposites were analogous to those of composites made with very short, chopped carbon fibers. About ten years ago, CNTs became available in continuous fiber form, allowing the investigation of processing methods to incorporate much higher CNT content into composites. Recent developments in CNT composite fabrication have yielded improved mechanical properties and indicate that while some standard carbon fiber composite processing methods can be used, CNTs have unique microstructural features that require techniques tailored to access the hierarchical features of CNTs. This presentation will discuss the evolution of structural CNT mechanical property advancements and collaboration mechanisms for research across federal labs, academia, and industry that accelerated the maturation of structural CNTs.

Carbon Nanotubes↗

Observations and Effects of Dipolarization Fronts Observed in Earth's Magnetotail

Dipolarization fronts in Earth's magnetotail are characterized by sharp jumps in magnetic field, a drop in density, and often follow earthward fast plasma flow. They are commonly detected near the equatorial plane of Earth s tail plasma sheet. Sometimes, but not always, dipolarization fronts are associated with global substorms and auroral brightenings. Both Cluster, THEMIS, and other spacecraft have detected dipolarization fronts in a variety of locations in the magnetotail. Using multi-spacecraft analyses together with simulations, we have investigated the propagation and evolution of some dipolarization events. We have also investigated the acceleration of electrons and ions that results from such magnetic-field changes. In some situations, the velocities of fast earthward flows are comparable to the Alfven speed, indicating that the flow bursts might have been generated by bursty reconnection that occurred tailward of the spacecraft. Based on multi-spacecraft timing analysis, dipolarization fronts are found to propagate mainly earthward at 160-335 km/s and have thicknesses of ~900-1500 km, which corresponds to the ion inertial length or gyroradius scale. Following the passage of dipolarization fronts, significant fluctuations are observed in the x and y components of the magnetic field. These peaks in the magnetic field come approximately 1-2 minutes after passage of the dipolarization front. These Bx and By fluctuations propagate primarily dawnward and earthward. Field-aligned electron beams are observed coincident with those magnetic field fluctuations. Non-Maxwellian electron and ion distributions are observed that are associated with the dipolarization that may be unstable to a range of electrostatic and/or whistler instabilities. Enhanced electrostatic broadband noise at frequencies below and near the lower-hybrid frequency is also observed at or very close to these fronts. This broadband noise is thought to play a role in further energizing the particles. Such studies provide insights into the particle acceleration mechanisms associated with substorm dipolarization, and, in turn, the effects of those acceleration mechanisms on the structure and evolution of dipolarization fronts.

Goldstein, Melvyn L.↗