Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power law process”

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 361 records · Page 20

Fermi Large Area Telescope Observations of the Supernova Remnant GS.7-0.1

We present a detailed analysis of the GeV gamma-ray emission toward the supernova remnant (SNR) G8.7-0.1 with the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope. An investigation of the relationship among G8.7-0.1 and the TeV unidentified source HESS J1804-216 provides us with an important clue on diffusion process of cosmic rays if particle acceleration operates in the SNR. The GeV gamma-ray emission is extended with most of the emission in positional coincidence with the SNR G8.7-0.1 and a lesser part located outside the western boundary of G8.7-0.1. The region of the gamma-ray emission overlaps spatially-connected molecular clouds, implying a physical connection for the gamma-ray structure. The total gamma-ray spectrum measured with LAT from 200 MeV-100 GeV can be described by a broken power-law function with a break of 2.4 +/- 0.6 (stat) +/- 1.2 (sys) GeV, and photon indices of2.10 +/- 0.06 (stat) +/- 0.10 (sys) below the break and 2.70 +/- 0.12 (stat) +/- 0.14 (sys) above the break. Given the spatial association among the gamma rays, the radio emission ofG8.7-0.1, and the molecular clouds, the decay of pions produced by particles accelerated in the SNR and hitting the molecular clouds naturally explains the GeV gamma-ray spectrum. We also find that the GeV morphology is not well represented by the TeV emission from HESS Jl804-2l6 and that the spectrum in the Ge V band is not consistent with the extrapolation of the TeV gamma-ray spectrum. The spectral index of the TeV emission is consistent with the particle spectral index predicted by a theory that assumes energy-dependent diffusion of particles accelerated in an SNR. We discuss the possibility that the TeV-spectrum originates from the interaction of particles accelerated in G8.7-0.l with molecular clouds, and we constrain the diffusion coefficient of the particles.

Ferrara, E. C.↗

The Role of Viscous Deformation in the Morphology of the Martian North Polar Cap

The first detailed topographic measurements of Mars' north polar cap have recently been collected by the Mars Orbiter Laser Altimeter (MOLA), an instrument on Mars Global Surveyor (MGS). These new data provide the opportunity to address in a quantitative way the processes responsible for the origin and evolution of the martian polar caps. On the basis of imaging and spectral observations from the Mariner 9 and Viking orbiters, it was recognized that a number of physical mechanisms. including radiative transfer, wind erosion, and viscous flow, represent possible or probable contributors to the morphology of the polar cap. Here we review rheological data for H2O and CO2 at conditions relevant to the martian polar regions. We then use MOLA topographic profiles to perform a preliminary assessment of the role of power law flow in contributing to the regional-scale planform of the north polar cap.

Zuber, M. T.↗

The 'Supercritical Pile' GRB Model: Afterglows and GRB, XRR, XRF Unification

We present the general notions and observational consequences of the "Supercritical Pile" GRB model; the fundamental feature of this model is a detailed process for the conversion of the energy stored in relativistic protons in the GRB Relativistic Blast Waves (RBW) into relativistic electrons and then into radiation. The conversion is effected through the $p \, \gamma \rightarrow p \, e circumflex + e circumflex -$ reaction, whose kinematic threshold is imprinted on the GRB spectra to provide a peak of their emitted luminosity at energy \Ep $\sim 1$ MeV in the lab frame. We extend this model to include, in addition to the (quasi--)thermal relativistic post-shock protons an accelerated component of power law form. This component guarantees the production of $e circumflex +e circumflex- - $pairs even after the RBW has slowed down to the point that its (quasi-) thermal protons cannot fulfill the threshold of the above reaction. We suggest that this last condition marks the transition from the prompt to the afterglow GRB phase. We also discuss conditions under which this transition is accompanied by a significant drop in the flux and could thus account for several puzzling, recent observations. Finally, we indicate that the same mechanism applied to the late stages of the GRB evolution leads to a decrease in \Ep $\propto \Gamma circumflex 2(t)\propto t circumflex {-3/4}$, a feature amenable to future observational tests.

Kazanas, D.↗

MPD thruster technology

MPD thrusters have demonstrated between 2000 and 7000 sec specific impulse at efficiencies approaching 40 percent, and have been operated continuously at power levels over 500 kW. These demonstrated capabilities, combined with the simplicity and robustness of the thruster, make them attractive candidates for application to both unmanned and manned orbit raising, lunar, and planetary missions. This work reviews the present status of MPD thruster research, including developments in the measured performance levels and electrode erosion rates, and theoretical studies of the thruster dynamics. Significant progress has been made in establishing empirical scaling laws, performance and lifetime limitations, and in the development of numerical codes to simulate the flowfield and the electrode processes.

Myers, Roger M.↗

ASCA and Radio/RXTE Observations of GX 339-4

We have analyzed three separate archival Advanced Satellite for Cosmology and Astrophysics (ASCA) observations and eight separate Rossi X-ray Timing Explorer (RXTE) observations of the black hole candidate GX 339-4 in its low luminosity, spectrally hard state. Three of the RXTE observations were strictly simultaneous with 843 MHz and 8.3- 9.1 GHz radio observations. All data sets show evidence for an approximately 6.4 keV Fe line with equivalent widths approximately 20-100 eV. 'Reflection models' show a hardening of the RXTE spectra with decreasing X-ray flux; however, these models do not exhibit evidence of a correlation between the photon index of the incident power law flux and the solid angle subtended by the reflector. None of the models fit to the X-ray data, however, simultaneously explain the observed radio properties. We argue that the spatial extent of the observed radio emission is at least 0(10(exp 7 GM/c squared). Timing analysis reveals that all observations save one show evidence of a persistent f(qpo approximately equals 0.3 Hz quasi-periodic oscillations(quasi-periodic oscillations)). The broad band (10-3-102 Hz) power appears to be dominated by two independent processes that can be modeled as very broad Lorentzians with Q approximately less than 1. Similar to Cyg X-1, the hard photon variability is seen to lag the soft photon vaxiability with the lag time increasing with decreasing Fourier frequency. The magnitude of this time lag is seen to be positively correlated with the flux of GX 339-4.

Nowak, M. A.↗

Chandra Observation of the Core of the Galaxy Cluster AWM 7

We present results from a Chandra observation of the core region of the nearby X-ray bright galaxy cluster AWM 7. There are blob-like substructures, which are seen in the energy band 2-10 keV, within 10 kpc (20") of the cD galaxy NGC 1129, and the brightest sub-peak has a spatial extent more than 4 kpc. We also notice that the central soft X-ray peak is offset from the optical center by 1.3 kpc. These structures have no correlated features in optical, infrared, or radio band. Energy spectrum of the hard sub-peak indicates a temperature higher than 3 keV with a metallicity less than 0.3 solar, or a power-law spectrum with photon index approximately 1.2. A hardness ratio map and a narrow Fe-K band image jointly indicate two Fe-rich blobs symmetrically located around the cD galaxy, with the direction perpendicular to the sub-peak direction. In larger scales (r less than 60 kpc), the temperature gradually drops from 4 keV to 2 keV toward the cluster center and the metal abundance rises steeply to a peak of 1.5 solar at r approximately equal to 7 kpc. These results indicate that a dynamical process is going on in the central region of AWM 7, which probably creates heated gas blobs and drives metal injection.

Furusho, T.↗

Corrosion Response of Ni-19Cr-5Fe in Molten NaCl-MgCl 2 Salt

Ni-based alloys are strong candidates for use in high-temperature molten salt reactors due to their superior corrosion resistance and mechanical stability compared to stainless steels. In this study, the static corrosion behavior of a Ni-19Cr-5Fe model alloy was systematically evaluated in a purified molten NaCl-MgCl 2 salt at 700°C for 30, 240, and 500 h. Post-exposure analyses were conducted to assess microstructural evolution, corrosion depth, and elemental depletion profiles. Corrosion rates, quantified by chromium depletion depth, followed an inverse power-law trend with increasing exposure time, indicating diffusion-limited kinetics. This trend is attributed to limitations in mass transport in the static salt and to the progressive local depletion of reactive chromium species. Coupled electron backscatter diffraction and energy dispersive X-ray spectroscopy analysis further revealed that the grain boundary character significantly influences corrosion susceptibility: high-angle grain boundaries exhibited pronounced Cr depletion and pitting, while low-angle and Σ3 boundaries remained comparatively resistant. These results offer valuable insight into the role of microstructure in corrosion processes and reinforce the importance of time-dependent material evaluation in molten salt environments relevant to advanced reactor designs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Absolute spectral energy distribution of quasi-stellar objects from 0.3 to 10 microns

Energy distributions of 24 quasrs have been measured from 0.3 to 10 microns. In addition, optical emission-line strengths have been obtained for 35 quasars. The sample includes quasars discovered from both radio and optical surveys. Over the observed range of wavelengths the energy distributions are roughly characterized by a power law of slope -1, but there are vast differences in the shapes of individual spectra. In the rest frame the average slope between 1 and 3 microns is significantly steeper than the average slope between 0.3 and 1 micron; the power emitted at 0.3 and 3 microns is significantly greater than that emitted around 1 micron. On the basis of arguments of variability and the shape of the energy distribution between optical and radio wavelengths, it is concluded that the bulk of the radiation in the observed range results from nonthermal processes. Some contribution from thermal emission from dust cannot be ruled out and may explain the break in the spectra at 1 micron seen in many quasars.

Neugebauer, G.↗

Analysis of the fetch-related bias in altimeter wind speed measurements

The present interpretation of the radar cross section sigma exp 0 measured by satellite altimeters implies that the rms wave slope gamma is controlled solely by the local wind. However, parameters of wave spectra, including the exponent in the power law for the equilibrium range, depend on sea maturity. The latter is characterized by the nondimensional fetch, x = gX/U-squared. Consequently, gamma and sigma exp 0 are controlled by both U and the wind fetch X. Geosat data for one year are used jointly with in-situ wind and wave observations to assess the fetch-related error trend in altimeter wind speeds. This trend results in overestimated winds in the regions and seasons characterized by a high x, and vice versa. A procedure for wind speed retrieval based on processing sigma exp 0 jointly with the significant wave height information contained in the altimeter wave forms is proposed.

Glazman, Roman E.↗

An estimate of the second law thermodynamic efficiency of the various units comprising an Environmental Control and Life Support System (ECLSS)

The objective of this paper is to present an estimate of the second law thermodynamic efficiency of the various units comprising an Environmental Control and Life Support System (ECLSS). The technique adopted here is based on an evaluation of the 'lost work' within each functional unit of the subsystem. Pertinent information for our analysis is obtained from a user interactive integrated model of an ECLSS. The model was developed using ASPEN. A potential benefit of this analysis is the identification of subsystems with high entropy generation as the most likely candidates for engineering improvements. This work has been motivated by the fact that the design objective for a long term mission should be the evaluation of existing ECLSS technologies not only the basis of the quantity of work needed for or obtained from each subsystem but also on the quality of work. In a previous study Brandhorst showed that the power consumption for partially closed and completely closed regenerable life support systems was estimated as 3.5 kw/individual and 10-12 kw/individual respectively. With the increasing cost and scarcity of energy resources, our attention is drawn to evaluate the existing ECLSS technologies on the basis of their energy efficiency. In general the first law efficiency of a system is usually greater than 50 percent. From literature, the second law efficiency is usually about 10 percent. The estimation of second law efficiency of the system indicates the percentage of energy degraded as irreversibilities within the process. This estimate offers more room for improvement in the design of equipment. From another perspective, our objective is to keep the total entropy production of a life support system as low as possible and still ensure a positive entropy gradient between the system and the surroundings. The reason for doing so is as the entropy production of the system increases, the entropy gradient between the system and the surroundings decreases, and the system will gradually approach equilibrium with the surroundings until it reaches the point where the entropy gradient is zero. At this point no work can be extracted from the system. This is called the 'dead state' of the system.

Chatterjee, Sharmista↗

State of the AGN: Progress Toward Understanding Black Hole Accretion Processes

Accretion of plasma onto black holes power some of the most powerful systems in the cosmos. Supermassive black holes at the centers of galaxies represent the high-mass limit of these objects and so account for the most luminous accretors. As a result, their influence spans vast spatial and temporal scales of cosmic phenomena: intracluster heating, intergalactic media, galactic feedback and star formation, kiloparsec-scale jets/outflows, variability over time scales of minutes to centuries, and luminous multi-wavelength electromagnetic emission extending all the way down to its event horizon. Their intrigue is heightened by the fact that they lie at the intersection of various physical laws---e.g., general relativistic gravity, magnetohydrodynamics, radiation, high-energy particle physics, thermodynamics, and photo-ionization---which all must be reconciled to arrive at a fundamental understanding and probe for new physics, like tests of general relativity. These physics ingredients must be incorporated into simulations performed somehow on dynamical scales ranging from that of the event horizon to parsec-scales. Fortunately, new computational and theoretical techniques---such as GPU computing, radiation transport, and and novel gridding techniques---are enabling progress to larger scales, more degrees of freedom, and even to binary systems. Some of the topics we will survey include recent progress on simulating the relationship between the disk-jet interaction, how tilted black holes behave, radiation-dominated flow, and how binary AGN affect the standard picture of black hole accretion. Along the way, we will highlight how new technologies have enabled these scientific rewards. Future directions and open questions will be provided to inspire discussion and interaction during the session.

LISA↗

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical “slaving onset” stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Environmental fatigue of an Al-Li-Cu alloy. Part 1: Intrinsic crack propagation kinetics in hydrogenous environments

Deleterious environmental effects on steady-state, intrinsic fatigue crack propagation (FCP) rates (da/dN) in peak aged Al-Li-Cu alloy 2090 are established by electrical potential monitoring of short cracks with programmed constant delta K and K(sub max) loading. The da/dN are equally unaffected by vacuum, purified helium, and oxygen but are accelerated in order of decreasing effectiveness by aqueous 1 percent NaCl with anodic polarization, pure water vapor, moist air, and NaCl with cathodic polarization. While da/dN depends on delta K(sup 4.0) for the inert gases, water vapor and chloride induced multiple power-laws, and a transition growth rate 'plateau'. Environmental effects are strongest at low delta K. Crack tip damage is ascribed to hydrogen embrittlement because of the following: (1) accelerated da/dN due to part-per-million levels of H2O without condensation; (2) impeded molecular flow model predictions of the measured water vapor pressure dependence of da/dN as affected by mean crack opening; (3) the lack of an effect of film-forming O2; (4) the likelihood for crack tip hydrogen production in NaCl, and (5) the environmental and delta K-process zone volume dependencies of the microscopic cracking modes. For NaCl, growth rates decrease with decreasing loading frequency, with the addition of passivating Li2CO3, and upon cathodic polarization. These variables increase crack surface film stability to reduce hydrogen entry efficiency. The hydrogen environmental FCP resistance of 2090 is similar to other 2000 series alloys and is better than 7075.

Piascik, Robert S.↗

Environmental fatigue of an Al-Li-Cu alloy. I - Intrinsic crack propagation kinetics in hydrogenous environments

Deleterious environmental effects on steady-state, intrinsic fatigue crack propagation (FCP) rates (da/dN) in peak aged Al-Li-Cu alloy 2090 are established by electrical potential monitoring of short cracks with programmed constant delta K and K(sub max) loading. The da/dN are equally unaffected by vacuum, purified helium, and oxygen but are accelerated in order of decreasing effectiveness of aqueous 1 percent NaCl with anodic polarization, pure water vapor, moist air, and NaCl with cathodic polarization. While da/dN depends on delta K(sup 4.0) for the inert gases, water vapor and chloride induced multiple power-laws, and a transition growth rate 'plateau'. Environmental effects are strongest at low delta K. Crack tip damage is ascribed to hydrogen embrittlement because of the following: (1) accelerated da/dN due to part-per-million levels of H2O without condensation; (2) impeded molecular flow model predictions of the measured water vapor pressure dependence of da/dN as affected by mean crack opening; (3) the lack of an effect of film-forming O2; (4) the likelihood for crack tip hydrogen production in NaCl; and (5) the environmental and delta K-process zone volume dependencies of the microscopic cracking modes. For NaCl, growth rates decrease with decreasing loading frequency, with the addition of passivating Li2CO3, and upon cathodic polarization. These variables increase crack surface film stability to reduce hydrogen entry efficiency. The hydrogen environmental FCP resistance of 2090 is similar to other 2000 series alloys and is better than 7075.

Piascik, Robert S.↗

Inverse Compton X-ray emission from the superluminal quasar 3C 345

In quasars with strong radio cores, the inverse-Compton process is believed to be the dominant source X-ray emission. For objects with parsec-scale radio jets, simple models have predicted that components in the jet emerging from the quasar nucleus generate the observed X-ray emission. We have tested this hypothesis in detail for the quasar 3C 345 using a ROSAT X-ray observation in 1990 July, together with quasi-simultaneous very long base interferometry (VLBI) imaging of the parsec-scale jet at five frequencies. The ROSAT spectrum is well fitted by a power law with index alpha = -0.96 +/- -0.13, consistent with models in which the X-ray emission results from inverse-Compton scattering of radio radiation from high-energy electrons in compact components. We show that the radio properties of brightest `knot' in the jet (`C5') can be fitted with a homogeneous sphere model whose parameters require bulk relativistic motion of the emitting material; otherwise the predicted model whose parameters require bulk relativistic motion of the emitting material; otherwise the predicted inverse-Compton X-ray emission exceeds the observed flux. If C5 is the origin of the X-ray emission, then it has a Doppler factor delta = 7.5((sup +3 sub -2)). If the nucleus or other components contribute to the X-ray emission, then this becomes a firm lower limit to delta. The inhomogeneous jet model of Koenigl is a good fit both to the barely resolved (less than 1 pc) flat-spectrum nucleus in the radio, and also to the ROSAT X-ray spectrum. The synchrotron and inverse-Compton emitting fluid moves down a narrow cone (opening angle 2 phi approximately 5 deg) nucleus relativistically, with delta approximately 4.6. Doppler factors for the nucleus and C5, derived from our ROSAT observation, provide evidence for bulk relativistic motion in the jet. By combining these constraints with well-known superluminal motion of jet components, we can deduce geometry. For epoch 1990.5 we infer the Lorentz factor gamma = 7.5 ((sup +1.0 sub -1.5)) and angle to the line of sight theta = 8((sup +2 deg sub -3 deg)) for H(sub 0) = 100 km/s/Mpc. These values are the most reliable yet derived using this method, because of the near-simultaneity of our X-ray and VLBI observations and the quality of the multifrequency of VLBI images and component radio spectra.

Unwin, S. C.↗

STITCH: A Subgrid-Scale Model for Energy Buildup in the Solar Corona

The solar corona routinely exhibits explosive activity, in particular coronal mass ejections and their accompanying eruptive ares, that have global-scale consequences. These events and their smaller counterparts, coronal jets, originate in narrow, sinuous lament channels. The key processes that form and evolve the channels operate on still smaller spatial scales and much longer time scales, culminating in a vast separation of characteristic lengths and times that govern these explosive phenomena. In this article, we describe implementation and tests of an efficient subgrid-scale model for generating eruptive structures in magnetohydrodynamics (MHD) coronal simulations. STITCH {STatistical InjecTion of Condensed Helicity { is a physics-based, reduced representation of helicity condensation: a process wherein small-scale vortical surface convection forms ubiquitous current sheets, and pervasive reconnection across the sheets mediates an inverse cascade of magnetic helicity and free energy, thereby forming the lament channels. We have developed a formalism, STITCH, that abstracts these complex processes into a single term in Ohm's law and the induction equation that directly injects tangential magnetic flux into the low corona. We show that our approach is in very good agreement with a full helicity-condensation calculation that treats all of the dynamics explicitly, while enabling substantial reductions in temporal duration and spatial resolution. In addition, we illustrate the flexibility of STITCH at forming localized lament channels and at energizing complex surface flux distributions that have sinuous boundaries. STITCH is simple to implement and computationally ecient, making it a powerful technique for physics-based modeling of solar eruptive events.

J T Dahlin↗

Development of Nanocrystalline Graphite from Lignin Sources

Carbon composites are attractive to a variety of high-impact applications, such as carbon fibers, batteries, and vehicle parts, due to their multifunctional properties. The properties of carbon are highly dependent on the allotrope the carbon takes and the functionality, impurities, and defects contained within the structure. The increase in demand for sustainable carbon sources in energy storage devices motivates interest in understanding synthesis parameters of lignin value-added products. Also, as the dependence on oil for fuel decreases, alternative sources for carbon in many applications will be needed. In this work, the thermochemical conversion of lignin powders from different feedstocks was evaluated via small and wide-angle X-ray scattering techniques to resolve the amorphous, disordered, and crystalline domains present in the lignin carbons. Scattering analyses indicated an evolution of hierarchical structures along with an increase in ordered domains as a function of carbonization temperature. Qualitative and quantitative methods were used to describe isotropic scattering intensity profiles at multiple length scales. The use of power law models in the mesoscopic region served as the basis to describe morphological changes related to structural features, for example, graphene stacking, degree of roughness, and surface fractals. Kraft softwood and switchgrass produced carbon powder with the most crystalline domains and the least surface roughness. Softwoods reached the highest degree of crystallinity followed by switchgrass samples and had less variability in particle sizes. These results suggest lignin carbons extracted from softwoods and switchgrass are viable substitutes for graphite. Interpretation of X-ray scattering data from lignin carbon powders elucidates feedstock- and processing-dependent morphological features across multiple length scales providing a straightforward framework to evaluate the feasibility of leveraging lignin carbons for producing tunable application-specific materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Developing and flight testing the HL-10 lifting body: A precursor to the Space Shuttle

The origins of the lifting-body idea are traced back to the mid-1950's, when the concept of a manned satellite reentering the Earth's atmosphere in the form of a wingless lifting body was first proposed. The advantages of low reentry deceleration loads, range capability, and horizontal landing of a lifting reentry vehicle (as compared with the high deceleration loads and parachute landing of a capsule) are presented. The evolution of the hypersonic HL-10 lifting body is reviewed from the theoretical design and development process to its selection as one of two low-speed flight vehicles for fabrication and piloted flight testing. The design, development, and flight testing of the low-speed, air-launched, rocket-powered HL-10 was part of an unprecedented NASA and contractor effort. NASA Langley Research Center conceived and developed the vehicle shape and conducted numerous theoretical, experimental, and wind-tunnel studies. NASA Flight Research Center (now NASA Dryden Flight Research Center) was responsible for final low-speed (Mach numbers less than 2.0) aerodynamic analysis, piloted simulation, control law development, and flight tests. The prime contractor, Northrop Corp., was responsible for hardware design, fabrication, and integration. Interesting and unusual events in the flight testing are presented with a review of significant problems encountered in the first flight and how they were solved. Impressions by the pilots who flew the HL-10 are included. The HL-10 completed a successful 37-flight program, achieved the highest Mach number and altitude of this class vehicle, and contributed to the technology base used to develop the space shuttle and future generations of lifting bodies.

Kempel, Robert W.↗