Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel planes”

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 973 records · Page 54

Computational heat transfer analysis for oscillatory channel flows

An accurate finite-difference scheme has been utilized to investigate oscillatory, laminar and incompressible flow between two-parallel-plates and in circular tubes. The two-parallel-plates simulate the regenerator of a free-piston Stirling engine (foil type regenerator) and the channel wall was included in the analysis (conjugate heat transfer problem). The circular tubes simulate the cooler and heater of the engine with an isothermal wall. The study conducted covered a wide range for the maximum Reynolds number (from 75 to 60,000), Valensi number (from 2.5 to 700), and relative amplitude of fluid displacement (0.714 and 1.34). The computational results indicate a complex nature of the heat flux distribution with time and axial location in the channel. At the channel mid-plane we observed two thermal cycles (out of phase with the flow) per each flow cycle. At this axial location the wall heat flux mean value, amplitude and phase shift with the flow are dependent upon the maximum Reynolds number, Valensi number and relative amplitude of fluid displacement. At other axial locations, the wall heat flux distribution is more complex.

Ibrahim, Mounir↗

Progress on the Plasmoid Thruster Experiment (PTX)

A plasmoid, also called a compact toroid, is a compact plasma structure with an integral magnetic field, that may be categorized according to the relative strength of the poloidal and toroidal magnetic field (Bp and Bt, respectively). An object with Bp/Bt much greater than 1 is called a Field Reverse Configuration (FRC); if Bp = Bt, it is called a Spheromak. The plasmoid thruster is a pulsed inductive device which operates by repetitively producing plasmoids that are accelerated and ejected at high velocity. As the process is inductive, this thruster avoids the problem of electrode erosion. Also, the magnetic structure of the plasmoid should suppress thermal and mass losses to the wall, and improve detachment of the plasma exhaust from the thruster. This concept should be capable of producing an Isp of 5,000 seconds and greater, with thrust densities of order 10(exp 5) N/sq m. The plasmoid thruster consists chiefly of a conical theta-pinch coil. Propellant is introduced onto a bias magnetic field, produced by an auxiliary coil, and is then pre-ionized, freezing in the magnetic field. The theta-pinch coil is then energized producing a field aligned anti- parallel to the bias field. The reversed field reconnects with the bias field to form the plasmoid. The magnetic pressure of the reversed field accelerates the plasmoid out of the thruster . A series of experiments have been conducted on the PTX device, which consisted of a single turn conical theta-pinch coil, driven by a 560 nF, 35 kV capacitor bank, which rang at a frequency of 500 kHz, and served all three functions required for formation: pre-ionization, bias field loading, and field reversal. Initial ionization was found to occur in an annular region at the exit plane of the coil, and was found to be reproducible with a variety of gases, including H2, D2, Ar, and an H2/N2 mixture (75% / 25%). A fast gas valve for injecting propellant has been tested, as well as a ringing pre-ionization circuit (operating at 5 Mhz) to better control the plasmoid formation. Velocities as high as 20 km/s have been measured for a plasmoid propagating into a static back-fill of neutral gas. Magnetic field and density measurements indicate that the plasmoid retains its structure well into the down-stream exhaust region, and that it is capable of pushing a dense cloud of neutral gas in front of it. Further experiments will focus on the effects of coil geometry and bias-field strength on the performance of the device. A high-power solid-state switching system is being developed in order to replace the spark-gap switch.

Martin, Adam↗

Spin bipolaron in the framework of emery model for high-T(sub c) copper oxide superconductors

The high-T(sub c) oxide compounds discovered recently exhibit a number of interesting physical properties. Two-dimensional antiferromagnetic spin order has been observed in these materials at the oxygen deficiency. This fact can be explained by strong correlation of the spins, situated on Cu sites in the conducting planes of the oxide superconductors. The doping or the oxygen deficiency lead to the occurrence of holes, occupying the oxygen p-orbitals according to the Emery model. At the small hole concentration they can move along the antiferromagnetic lattice of spins, localized on Cu sites. Researchers consider the two holes situation and describe in what way their behavior depends on the antiferromagnetic exchange interation J. It is known that in the framework of Hubbard model with strong on-site Coulomb repulsion, a single hole can form a spin polaron of the large radius. It is reasonable to admit that two holes with parallel spins (triplet) form the spin bipolaron complex owing to the hole excitations' capability to polarize Cu spin surroundings. Such an excitation was considered in the phenomenological way. Here the problem is discussed on the basis of the microscopic approach in the framework of the variational principle. A special kind of wave function is used for such a purpose. The wave function is constructed by generalizing the trial functions proposed in over two holes excitation situation (triplet) and then the region of spin bipolaron existance in the framework of Emery model is studied. In this model the Hamiltonian can be easily rewritten by forming the oxygen states transforming as the irreducible representations of the group D(sub 4).

Golub, A. A.↗

Lower-Hybrid Drift Waves and Their Interaction With Plasmas in A 3D Symmetric Reconnection Simulation With Zero Guide Field

We investigate lower-hybrid drift waves (LHDW) in symmetric magnetic reconnection with zero guide field using three-dimensional particle-in-cell simulations. The long-wavelength mode with develops in the bifurcated electron current layer around the X-line within the width of the electron meandering motion from the mid-plane, where is the ion (electron) gyroradius. The short-wavelength mode with develops in the separatrix region downstream of the electron outflow jet, producing electron vortices in the background flow frame. Electrons follow the E × B drift with corrections from the diamagnetic drift and are heated inside the vortices with diverging electric fields. In the vortices, ions have comparable E × B and inertia drifts, which together mostly cancel the diamagnetic drift. Toward the center of diverging field vortices, ions are decelerated, leading to a decrease in the perpendicular temperature, while the loss of low-energy ions results in an increase in the parallel temperature. Parallel electric fields exist as a combination of the LHDW wave field projected to the magnetic field direction and the penetration of whistler waves that are mainly outside of the LHDW layer. The magnetic flux tube is twisted in the vortices. The twist may potentially lead to slippage reconnection, as indicated by the non-uniform parallel potential variation across field lines, while the periodic variations of the twisting directions are a limiting factor.

Shan Wang↗

Bit-Serial Adder Based on Quantum Dots

A proposed integrated circuit based on quantum-dot cellular automata (QCA) would function as a bit-serial adder. This circuit would serve as a prototype building block for demonstrating the feasibility of quantum-dots computing and for the further development of increasingly complex and increasingly capable quantum-dots computing circuits. QCA-based bit-serial adders would be especially useful in that they would enable the development of highly parallel and systolic processors for implementing fast Fourier, cosine, Hartley, and wavelet transforms. The proposed circuit would complement the QCA-based circuits described in "Implementing Permutation Matrices by Use of Quantum Dots" (NPO-20801), NASA Tech Briefs, Vol. 25, No. 10 (October 2001), page 42 and "Compact Interconnection Networks Based on Quantum Dots" (NPO-20855), which appears elsewhere in this issue. Those articles described the limitations of very-large-scale-integrated (VLSI) circuitry and the major potential advantage afforded by QCA. To recapitulate: In a VLSI circuit, signal paths that are required not to interact with each other must not cross in the same plane. In contrast, for reasons too complex to describe in the limited space available for this article, suitably designed and operated QCA-based signal paths that are required not to interact with each other can nevertheless be allowed to cross each other in the same plane without adverse effect. In principle, this characteristic could be exploited to design compact, coplanar, simple (relative to VLSI) QCA-based networks to implement complex, advanced interconnection schemes. To enable a meaningful description of the proposed bit-serial adder, it is necessary to further recapitulate the description of a quantum-dot cellular automation from the first-mentioned prior article: A quantum-dot cellular automaton contains four quantum dots positioned at the corners of a square cell. The cell contains two extra mobile electrons that can tunnel (in the quantum-mechanical sense) between neighboring dots within the cell. The Coulomb repulsion between the two electrons tends to make them occupy antipodal dots in the cell. For an isolated cell, there are two energetically equivalent arrangements (denoted polarization states) of the extra electrons. The cell polarization is used to encode binary information. Because the polarization of a nonisolated cell depends on Coulomb-repulsion interactions with neighboring cells, universal logic gates and binary wires could be constructed, in principle, by arraying QCA of suitable design in suitable patterns. Again, for reasons too complex to describe here, in order to ensure accuracy and timeliness of the output of a QCA array, it is necessary to resort to an adiabatic switching scheme in which the QCA array is divided into subarrays, each controlled by a different phase of a multiphase clock signal. In this scheme, each subarray is given time to perform its computation, then its state is frozen by raising its inter-dot potential barriers and its output is fed as the input to the successor subarray. The successor subarray is kept in an unpolarized state so it does not influence the calculation of preceding subarray. Such a clocking scheme is consistent with pipeline computation in the sense that each different subarray can perform a different part of an overall computation. In other words, QCA arrays are inherently suitable for pipeline and, moreover, systolic computations. This sequential or pipeline aspect of QCA would be utilized in the proposed bit-serial adders.

Fijany, Amir↗

On the Presence and Thermalization of Cold Ions in the Exhaust of Antiparallel Symmetric Reconnection

Using fully kinetic 2.5 dimensional particle-in-cell simulations of anti-parallel symmetric magnetic reconnection, we investigate how initially cold ions are captured by the reconnection process, and how they evolve and behave in the exhaust. We find that initially cold ions can remain cold deep inside the exhaust. Cold ions that enter the exhaust downstream of active separatrices, closer to the dipolarization front, appear as cold counter-streaming beams behind the front. In the off-equatorial region, these cold ions generate ion-acoustic waves that aid in the thermalization both of the incoming and outgoing populations. Closest to the front, due to the stronger magnetization, the ions can remain relatively cold during the neutral plane crossing. In the intermediate exhaust, the weaker magnetization leads to enhanced pitch angle scattering and reflection. Cold ions that enter the exhaust closer to the X line, at active separatrices, evolve into a thermalized exhaust. Here, the cold populations are heated through a combination of thermalization at the separatrices and pitch angle scattering in the curved magnetic field around the neutral plane. Depending on where the ions enter the exhaust, and how long time they have spent there, they are accelerated to different energies. The superposition of separately thermalized ion populations that have been accelerated to different energies form the hot exhaust population.

Cecilia Norgren↗

Comprehensive Modeling and Analysis of Rotorcraft Variable Speed Propulsion System With Coupled Engine/Transmission/Rotor Dynamics

This project develops comprehensive modeling and simulation tools for analysis of variable rotor speed helicopter propulsion system dynamics. The Comprehensive Variable-Speed Rotorcraft Propulsion Modeling (CVSRPM) tool developed in this research is used to investigate coupled rotor/engine/fuel control/gearbox/shaft/clutch/flight control system dynamic interactions for several variable rotor speed mission scenarios. In this investigation, a prototypical two-speed Dual-Clutch Transmission (DCT) is proposed and designed to achieve 50 percent rotor speed variation. The comprehensive modeling tool developed in this study is utilized to analyze the two-speed shift response of both a conventional single rotor helicopter and a tiltrotor drive system. In the tiltrotor system, both a Parallel Shift Control (PSC) strategy and a Sequential Shift Control (SSC) strategy for constant and variable forward speed mission profiles are analyzed. Under the PSC strategy, selecting clutch shift-rate results in a design tradeoff between transient engine surge margins and clutch frictional power dissipation. In the case of SSC, clutch power dissipation is drastically reduced in exchange for the necessity to disengage one engine at a time which requires a multi-DCT drive system topology. In addition to comprehensive simulations, several sections are dedicated to detailed analysis of driveline subsystem components under variable speed operation. In particular an aeroelastic simulation of a stiff in-plane rotor using nonlinear quasi-steady blade element theory was conducted to investigate variable speed rotor dynamics. It was found that 2/rev and 4/rev flap and lag vibrations were significant during resonance crossings with 4/rev lagwise loads being directly transferred into drive-system torque disturbances. To capture the clutch engagement dynamics, a nonlinear stick-slip clutch torque model is developed. Also, a transient gas-turbine engine model based on first principles mean-line compressor and turbine approximations is developed. Finally an analysis of high frequency gear dynamics including the effect of tooth mesh stiffness variation under variable speed operation is conducted including experimental validation. Through exploring the interactions between the various subsystems, this investigation provides important insights into the continuing development of variable-speed rotorcraft propulsion systems.

DeSmidt, Hans A.↗

Software for Use with Optoelectronic Measuring Tool

A computer program has been written to facilitate and accelerate the process of measurement by use of the apparatus described in "Optoelectronic Tool Adds Scale Marks to Photographic Images" (KSC-12201). The tool contains four laser diodes that generate parallel beams of light spaced apart at a known distance. The beams of light are used to project bright spots that serve as scale marks that become incorporated into photographic images (including film and electronic images). The sizes of objects depicted in the images can readily be measured by reference to the scale marks. The computer program is applicable to a scene that contains the laser spots and that has been imaged in a square pixel format that can be imported into a graphical user interface (GUI) generated by the program. It is assumed that the laser spots and the distance(s) to be measured all lie in the same plane and that the plane is perpendicular to the line of sight of the camera used to record the image

Ballard, Kim C.↗

Optical Flow-Field Techniques Used for Measurements in High-Speed Centrifugal Compressors

The overall performance of a centrifugal compressor depends on the performance of the impeller and diffuser as well as on the interactions occurring between these components. Accurate measurements of the flow fields in each component are needed to develop computational models that can be used in compressor design codes. These measurements must be made simultaneously over an area that covers both components so that researchers can understand the interactions occurring between the two components. Optical measurement techniques are being used at the NASA Lewis Research Center to measure the velocity fields present in both the impeller and diffuser of a 4:1 pressure ratio centrifugal compressor operating at several conditions ranging from design flow to surge. Laser Doppler Velocimetry (LDV) was used to measure the intrablade flows present in the impeller, and the results were compared with analyses obtained from two three-dimensional viscous codes. The development of a region of low throughflow velocity fluid within this high-speed impeller was examined and compared with a similar region first observed in a large low-speed centrifugal impeller at Lewis. Particle Image Velocimetry (PIV) is a relatively new technique that has been applied to measuring the diffuser flow fields. PIV can collect data rapidly in the diffuser while avoiding the light-reflection problems that are often encountered when LDV is used. The Particle Image Velocimeter employs a sheet of pulsed laser light that is introduced into the diffuser in a quasi-radial direction through an optical probe inserted near the diffuser discharge. The light sheet is positioned such that its centerline is parallel to the hub and shroud surfaces and such that it is parallel to the diffuser vane, thereby avoiding reflections from the solid surfaces. Seed particles small enough to follow the diffuser flow are introduced into the compressor at an upstream location. A high-speed charge-coupled discharge (CCD) camera is synchronized to the laser pulse rate; this allows it to capture images of seed particle position that are separated by a small increment in time. A crosscorrelation of a particle's position in two consecutive images provides an estimate of flow velocity and direction. Multiple image pairs obtained in rapid succession at a particular flow condition provide enough measurements for statistical significance. PIV provides simultaneous velocity measurements over the entire plane that is illuminated by the light sheet instead of at a single point, as is the case when LDV is used. PIV has a further advantage in that the laser light pulse can be triggered by an external source such as a high-response pressure transducer. This feature will allow PIV to synchronize flow imaging to physical phenomena such as rotating stall or stall precursor waves. We hope that this technique can be used to obtain images of the flow field during and just prior to stall.

Skoch, Gary J.↗

Parallel arrays of Josephson junctions for submillimeter local oscillators

In this paper we discuss the influence of the DC biasing circuit on operation of parallel biased quasioptical Josephson junction oscillator arrays. Because of nonuniform distribution of the DC biasing current along the length of the bias lines, there is a nonuniform distribution of magnetic flux in superconducting loops connecting every two junctions of the array. These DC self-field effects determine the state of the array. We present analysis and time-domain numerical simulations of these states for four biasing configurations. We find conditions for the in-phase states with maximum power output. We compare arrays with small and large inductances and determine the low inductance limit for nearly-in-phase array operation. We show how arrays can be steered in H-plane using the externally applied DC magnetic field.

Pance, Aleksandar↗

Structural anomalies in undoped Gallium Arsenide observed in high resolution diffraction imaging with monochromatic synchrotron radiation

Novel, streak-like disruption features restricted to the plane of diffraction have recently been observed in images obtained by synchrotron radiation diffraction from undoped, semi-insulating gallium arsenide crystals. These features were identified as ensembles of very thin platelets or interfaces lying in (110) planes, and a structural model consisting of antiphase domain boundaries was proposed. We report here the other principal features observed in high resolution monochromatic synchrotron radiation diffraction images: (quasi) cellular structure; linear, very low-angle subgrain boundaries in (110) directions, and surface stripes in a (110) direction. In addition, we report systematic differences in the acceptance angle for images involving various diffraction vectors. When these observations are considered together, a unifying picture emerges. The presence of ensembles of thin (110) antiphase platelet regions or boundaries is generally consistent not only with the streak-like diffraction features but with the other features reported here as well. For the formation of such regions we propose two mechanisms, operating in parallel, that appear to be consistent with the various defect features observed by a variety of techniques.

Steiner, B.↗

Structural anomalies in undoped gallium arsenide observed in high-resolution diffraction imaging with monochromatic synchrotron radiation

Novel, streak-like disruption features restricted to the plane of diffraction have recently been observed in images obtained by synchrotron radiation diffraction from undoped, semi-insulating gallium arsenide crystals. These features were identified as ensembles of very thin platelets or interfaces lying in (110) planes, and a structural model consisting of antiphase domain boundaries was proposed. We report here the other principal features observed in high resolution monochromatic synchrotron radiation diffraction images: (quasi) cellular structure; linear, very low-angle subgrain boundaries in (110) directions, and surface stripes in a (110) direction. In addition, we report systematic differences in the acceptance angle for images involving various diffraction vectors. When these observations are considered together, a unifying picture emerges. The presence of ensembles of thin (110) antiphase platelet regions or boundaries is generally consistent not only with the streak-like diffraction features but with the other features reported here as well. For the formation of such regions we propose two mechanisms, operating in parallel, that appear to be consistent with the various defect features observed by a variety of techniques.

Steiner, B.↗

Chaos and Beyond in a Water Filled Ultrasonic Resonance System

Finite amplitude ultrasonic wave resonances in a one-dimensional liquid-filled cavity, formed by a narrow band transducer and a plane reflector, are reported. The resonances are observed to include not only the expected harmonic and subharmonic signals (1,2) but chaotic signals as well. The generation mechanism requires attaining a threshold value of the driving amplitude that the liquid-filled cavity system becomes sufficiently nonlinear in response. The nonlinear features of the system were recently investigated via the construction of an ultrasonic interferometer having optical precision. The transducers were compressional, undamped quartz and lithium niobate crystals having the frequency range 1-10 MHz, driven by a high power amplifier. Both an optical diffraction system to characterize the diffraction pattern of laser light normally incident to the cavity and a receiving transducer attached to an aligned reflector with lapped flat and parallel surfaces were used to assess the generated resonance response in the cavity. At least 5 regions of excitation are identified.

Lazlo, Adler↗

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: 1) a k-epsilon turbulence model computation on a 6.8 million point grid using wall functions, 2) a k-epsilon turbulence model computation on a 14 million point grid integrating to the wall, and 3) a k-omega turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-epsilon and k-omega computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-epsilon computations which extended from 70% to 100% span at the exit rating plane, whereas the k-omega computation had reversed flow from 95% to 100% span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-epsilon case resulted in an under-prediction in adiabatic efficiency of 8.3 points, whereas the k-omega case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: (1) a k-ε turbulence model computation on a 6.8 million point grid using wall functions, (2) a k-ε turbulence model computation on a 14 million point grid integrating to the wall, and (3) a k-ω turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-ε and k-ω computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-ε computations which extended from 70 to 100 percent span at the exit rating plane, whereas the k-ω computation had reversed flow from 95 to 100 percent span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-ε case resulted in an underprediction in adiabatic efficiency of 8.3 points, whereas the k-ω case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗

SOFIA Far-infrared Imaging Polarimetry of M82 and NGC 253: Exploring the Supergalactic Wind

We present far-infrared polarimetry observations of M82 at 53 and 154μm and NGC 253 at 89μm, which were taken with High-resolution Airborne Wideband Camera-plus(HAWC+)in polarimetry mode on the Stratospheric Observatory for Infrared Astronomy. The polarization of M82 at 53μm clearly shows a magnetic field geometry perpendicular to the disk in the hot dust emission. For M82 the polarization at 154μm shows a combination of field geometry perpendicular to the disk in the nuclear region, but closer to parallel to the disk away from the nucleus. The fractional polarization at 53μm(154μm)ranges from 7%(3%)off nucleus to 0.5%(0.3%)near the nucleus. A simple interpretation of the observations of M82 invokes a massive polar outflow, dragging the field along, from a region∼700pc in diameter that has entrained some of the gas and dust, creating a vertical field geometry seen mostly in the hotter(53μm)dust emission. This outflow sits within a larger disk with a more typical planar geometry that more strongly contributes to the cooler(154μm)dust emission. For NGC 253, the polarization at 89μm is dominated by a planar geometry in the tilted disk, with weak indication of a vertical geometry above and below the plane from the nucleus. The polarization observations of NGC 253 at 53μm were of a insufficient signal-to-noise ratio for a detailed analysis.

Galaxies↗

Fluid/chemistry modeling for hypersonic flight analysis

Design studies are underway for a variety of hypersonic flight vehicles. The National Aero-Space Plane will provide a reusable, single-stage-to-orbit capability for routine access to low earth orbit. Flight-capable satellites will dip into the atmosphere to maneuver to new orbits, while planetary probes will decelerate at their destination by atmospheric aerobraking. To supplement limited experimental capabilities in the hypersonic regime, computational fluid dynamics is being used to analyze the flow about these configurations. The governing equations include fluid dynamic as well as chemical species equations, which are being solved with new, robust numerical algorithms. Examples of CFD applications to hypersonic vehicles suggest an important role this technology will play in the development of future aerospace systems. The computational resources needed to obtain solutions are large, but solution-adaptive grids, convergence acceleration, and parallel processing may make run times manageable.

Edwards, Thomas A.↗

CFD analysis of hypersonic, chemically reacting flow fields

Design studies are underway for a variety of hypersonic flight vehicles. The National Aero-Space Plane will provide a reusable, single-stage-to-orbit capability for routine access to low earth orbit. Flight-capable satellites will dip into the atmosphere to maneuver to new orbits, while planetary probes will decelerate at their destination by atmospheric aerobraking. To supplement limited experimental capabilities in the hypersonic regime, computational fluid dynamics (CFD) is being used to analyze the flow about these configurations. The governing equations include fluid dynamic as well as chemical species equations, which are being solved with new, robust numerical algorithms. Examples of CFD applications to hypersonic vehicles suggest an important role this technology will play in the development of future aerospace systems. The computational resources needed to obtain solutions are large, but solution adaptive grids, convergence acceleration, and parallel processing may make run times manageable.

Edwards, T. A.↗