Divergent-flow contact-ionization electrostatic thrustor for satellite attitude control and station keeping
Divergent-flow contact-ionization electrostatic thrustor for satellite attitude control and stationkeeping
SEARCH · Search NASA
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.
Divergent-flow contact-ionization electrostatic thrustor for satellite attitude control and stationkeeping
Experiments have been made at Stanford University to determine the performance characteristics of plane-wall, two-dimensional diffusers which were so proportioned as to insure reasonable approximation of two-dimensional flow. All of the diffusers had identical entrance cross sections and discharged directly into a large plenum chamber; the test program included wide variations of divergence angle and length. During all tests a dynamic pressure of 60 pounds per square foOt was maintained at the diffuser entrance and the boundary layer there was thin and fully turbulent. The most interesting flow characteristics observed were the occasional appearance of steady, unseparated, asymmetric flow - which was correlated with the boundary-layer coalescence - and the rapid deterioration of flow steadiness - which occurred as soon as the divergence angle for maximum static pressure recovery was exceeded. Pressure efficiency was found to be controlled almost exclusively by divergence angle, whereas static pressure recovery was markedly influenced by area ratio (or length) as well as divergence angle. Volumetric efficiency. diminished as area ratio increased, and at a greater rate with small lengths than with large ones. Large values of the static-pressure-recovery coefficient were attained only with long diffusers of large area ratio; under these conditions pressure efficiency was high and. volumetric efficiency low. Auxiliary tests with asymmetric diffusers demonstrated that longitudinal pressure gradient, rather than wall divergence angle, controlled flow separation. Others showed that the addition of even a short exit duct of uniform section augmented pressure recovery. Finally, it was found that the installation of a thin, central, longitudinal partition suppressed flow separation in short diffusers and thereby improved pressure recovery
In a previous work, embedded ensemble propagation was proposed to improve the efficiency of sampling-based uncertainty quantification methods of computational models on emerging computational architectures. It consists of simultaneously evaluating the model for a subset of samples together, instead of evaluating them individually. A first approach introduced to solve parametric linear systems with ensemble propagation is ensemble reduction. In Krylov methods for example, this reduction consists in coupling the samples together using an inner product that sums the sample contributions. Ensemble reduction has the advantages of being able to use optimized implementations of BLAS functions and having a stopping criterion which involves only one scalar. However, the reduction potentially decreases the rate of convergence due to the gathering of the spectra of the samples. In this paper, we investigate a second approach: ensemble propagation without ensemble reduction in the case of GMRES. This second approach solves each sample simultaneously but independently to improve the convergence compared to ensemble reduction. This raises two new issues which are solved in this paper: the fact that optimized implementations of BLAS functions cannot be used anymore and that ensemble divergence, whereby individual samples within an ensemble must follow different code execution paths, can occur. We tackle those issues by implementing a high-performing ensemble GEMV and by using masks. The proposed ensemble GEMV leads to a similar cost per GMRES iteration for both approaches, i.e. with and without reduction. For illustration, we study the performances of the new linear solver in the context of a mesh tying problem. Furthermore, this example demonstrates improved ensemble propagation speed-up without reduction.
An experimental investigation has been conducted at the NASA Langley 16-Foot Transonic Tunnel Static Test Facility to determine the concept feasibility of using fluidics to achieve multiplane thrust vector control in a 2D convergent-divergent (2D-CD) fixed aperture nozzle. Pitch thrust vector control is achieved by injection of flow through a slot in the divergent flap into the primary nozzle flow stream. Yaw vector control results from secondary air delivered tangentially to vertical Coanda flaps. These flaps are offset laterally and aligned parallel to the primary nozzle side walls. All tests were conducted at static (no external flow) conditions. Flow visualization was conducted using a paint flow technique and Focus Schlieren. Significant levels of pitch deflection angles (19 deg) were achieved at low pressure ratios and practical levels (14 deg) resulted at typical intermediate power settings. The ability of the Coanda surface blowing concept to produce yaw deflection was limited to NPR not greater than 4.
The problem of determining whether or not a point lies inside a given polygon occurs in many applications. In air traffic management concepts, a correct solution to the point-in-polygon problem is critical to geofencing systems for Unmanned Aerial Vehicles and in weather avoidance applications. Many mathematical methods can be used to solve the point-in-polygon problem. Unfortunately, a straightforward floating- point implementation of these methods can lead to incorrect results due to round-off errors. In particular, these errors may cause the control flow of the program to diverge with respect to the ideal real-number algorithm. This divergence potentially results in an incorrect point-in- polygon determination even when the point is far from the edges of the polygon. This paper presents a provably correct implementation of a point-in-polygon method that is based on the computation of the winding number. This implementation is mechanically generated from a source- to-source transformation of the ideal real-number specification of the algorithm. The correctness of this implementation is formally verified within the Frama-C analyzer, where the proof obligations are discharged using the Prototype Verification System (PVS).
In floating-point programs, test instability occurs when the control flow of a conditional statement diverges from its ideal execution under real arithmetic. This phenomenon is caused by the presence of round-off errors in floating-point computations. Writing programs that correctly handle test instability often require expertise on finite precision computations and rounding errors. This paper presents a fully automatic tool chain that generates and formally verifies a test-stable floating-point C program from its functional specification in real arithmetic. The generated program is instrumented to soundly detect when unstable tests may occur and, in these cases, to issue a warning. The proposed approach combines the PRECiSA floating-point static analyzer, the Frama-C software verification suite, and the PVS theorem prover.
An investigation was conducted in the model preparation area of the Langley 16-Foot Transonic Tunnel to study a passive cavity concept for improving the off-design performance of fixed-geometry exhaust nozzles. Passive cavity ventilation (through a porous surface) was applied to divergent flap surfaces and tested at static conditions in a sub-scale, nonaxisymmetric, convergent-divergent nozzle. As part of a comprehensive investigation, force, moment and pressure measurements were taken and focusing schlieren flow visualization was obtained for a baseline configuration and D passive cavity configurations. All tests were conducted with no external flow and high-pressure air was used to simulate jet-exhaust flow at nozzle pressure ratios from 1.25 to approximately 9.50. Results indicate that baseline nozzle performance was dominated by unstable shock-induced boundary-layer separation at off-design conditions, which came about through the natural tendency of overexpanded exhaust flow to satisfy conservation requirements by detaching from the nozzle divergent flaps. Passive cavity ventilation added the ability to control off-design separation in the nozzle by either alleviating separation or encouraging stable separation of the exhaust flow. Separation alleviation offers potential for installed nozzle performance benefits by reducing drag at forward flight speeds, even though it may reduce off-design static thrust efficiency as much as 3.2 percent. Encouraging stable separation of the exhaust flow offers significant performance improvements at static, low NPR and low Mach number flight conditions by improving off-design static thrust efficiency as much as 2.8 percent. By designing a fixed-geometry nozzle with fully porous divergent flaps, where both cavity location and percent open porosity of the flaps could be varied, passive flow control would make it possible to improve off-design nozzle performance across a wide operating range. In addition, the ability to encourage separation on one flap while alleviating it on the other makes it possible to generate thrust vectoring in the nozzle through passive flow control.
M = 0.83 Laminar Flow Control (LFC) transports, carrying large percentage payloads over a range of 20000 kilometers at cruise L/D's of 39 appear feasible with large space externally braced wings, external fuel pods, active controls, and 70 percent laminar flow on wing and tail surfaces, engine nacelles and struts, and a turbulent fuselage. A combination of a swept-forward inboard and a swept-back outer wing appears superior overall, especially for laminar flow and eliminating leading edge contamination probably caused by flyspecks and ice crystals. Wing divergence appears controllable by a combination of various methods. Wind-mounted superfans with extensive laminar flow on their nacelles appear practical. Their dominant tone noise is below the frequency range of the most strongly amplified TS-waves.
Results are presented of an experimental program that investigated the use of a secondary air stream to control the amount of flow through a convergent-divergent nozzle. These static tests utilized high pressure, ambient temperature air that was injected at the throat of the nozzle through an annular slot. Multiple injection slot sizes and injection angles were tested. The introduction of secondary flow was made in an opposing direction to the primary flow and the resulting flow field caused the primary stream to react as though the physical throat size had been reduced. The percentage reduction in primary flow rate was generally about twice the injected flow rate. The most effective throttling was achieved by injecting through the smallest slot in an orientation most nearly opposed to the approaching primary flow. Thrust edliciency, as measured by changes in nozzle thrust coefficient, was highest at high nozzle pressure ratios, NPR. The static test results agreed with predictions obtained prior from PABSD, a fully viscous computational fluid dynamics program. Since use of such an injection system on gas turbine engine exhaust nozzles would be primarily at high NPRs, it was concluded that fluidic control holds promise for reducing nozzle weight and complexity on future systems.
The NASA DAIDALUS library provides formal definitions for Detect-and-Avoid avionics concepts such as when an aircraft is well-clear with respect to the surrounding air traffic, i.e., it does not operate in such proximity to create a collision hazard. While several properties are proven correct for DAIDALUS assuming ideal real number arithmetic, an actual implementation that uses floating-point numbers may behave unexpectedly because of round-off errors and run-time exceptions. This paper presents an experience report on the application of a formal methods toolchain to extract and verify floating-point C code from a real-valued specification of the well-clear module of DAIDALUS. This toolchain comprises the PVS theorem prover, the PRECiSA floating-point analyzer and code generator, and the Frama-C analysis suite. The generated code is automatically instrumented to detect when the control flow of the floating-point program may diverge from the ideal real number specification, and it is annotated with contracts that state the maximum accumulated round-off error. The absence of overflows is also formally verified for the generated code. In order to apply the toolchain to an industrial case study such as DAIDALUS, a formally verified pre-processing of the input specification is performed, which includes a program slicing and several semantic-preserving simplifications.
The required type and amount of numerical dissipation/filter to accurately resolve all relevant multiscales of complex MHD unsteady high-speed shock/shear/turbulence/combustion problems are not only physical problem dependent, but also vary from one flow region to another. In addition, proper and efficient control of the divergence of the magnetic field (Div(B)) numerical error for high order shock-capturing methods poses extra requirements for the considered type of CPU intensive computations. The goal is to extend our adaptive numerical dissipation control in high order filter schemes and our new divergence-free methods for ideal MHD to non-ideal MHD that include viscosity and resistivity. The key idea consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. These scheme-independent detectors are capable of distinguishing shocks/shears, flame sheets, turbulent fluctuations and spurious high-frequency oscillations. The detection algorithm is based on an artificial compression method (ACM) (for shocks/shears), and redundant multiresolution wavelets (WAV) (for the above types of flow feature). These filters also provide a natural and efficient way for the minimization of Div(B) numerical error.
Electrostatic ion thrustor employing cesium on tungsten contact ionization for satellite attitude control and stationkeeping missions
While Monte Carlo Neutron Transport (MCNT) is near-embarrasingly parallel, the effectively unpredictable lifetime of neutrons can lead to divergence when MCNT is evaluated on GPUs. Divergence is the phenomenon of adjacent threads in a warp executing different control flow paths; on GPUS, it reduces performance because each work group may only execute one path at a time. The process of Thread Data Remapping (TDR) resolves these discrepancies by moving data across hardware such that data in the same warp will be processed through similar paths. A common issue among prior implementations of TDR is the synchronous nature of its remapping and processing cycles, which exhaustively sort data produced by prior processing passes and exhaustively evaluate the sorted data. In another work, we defined a method of remapping data through an asynchronous scheduler which allows for work to be stored in shared memory and deferred arbitrarily until that work is a viable option for low-divergence evaluation. This article surveys a wider set of cases, with the goal of characterizing performance trends across a more comprehensive set of parameters. These parameters include cross sections of scattering/capturing/fission, use of implicit capture, source neutron counts, simulation time spans, and tuned memory allocations. Across these cases, we have recorded minimum and average execution times, as well as a heuristically tuned near-optimal memory allocation size for both synchronous and asynchronous scheduling. Across the collected data, it is shown that the asynchronous method is faster and more memory efficient in the majority of cases, and that it requires less tuning to achieve competitive performance.
A method and apparatus is disclosed for levitating a specimen of material in a containerless environment at a stable nodal position independent of gravity. An elongated levitation tube has a contoured interior in the form of convergent section, constriction, and a divergent section in which the levitation node is created. A gas flow control means prevents separation of flow from the interior walls in the region of a specimen. The apparatus provides for levitating and heating the specimen simultaneously by combustion of a suitable gas mixture combined with an inert gas.
A three-Dimensional computational model is used to simulate flow in a non-axisymmetric, convergent-divergent nozzle incorporating porous cavities for shock-boundary layer interaction control. The nozzle has an expansion ratio (exit area/throat area) of 1.797 and a design nozzle pressure ratio of 8.78. Flow fields for the baseline nozzle (no porosity) and for the nozzle with porous surfaces of 10% openness are computed for Nozzle Pressure Ratio (NPR) varying from 1.29 to 9.54. The three dimensional computational results indicate that baseline (no porosity) nozzle performance is dominated by unstable, shock-induced, boundary-layer separation at over-expanded conditions. For NPR less than or equal to 1.8, the separation is three dimensional, somewhat unsteady, and confined to a bubble (with partial reattachment over the nozzle flap). For NPR greater than or equal to 2.0, separation is steady and fully detached, and becomes more two dimensional as NPR increased. Numerical simulation of porous configurations indicates that a porous patch is capable of controlling off design separation in the nozzle by either alleviating separation or by encouraging stable separation of the exhaust flow. In the present paper, computational simulation results, wall centerline pressure, mach contours, and thrust efficiency ratio are presented, discussed and compared with experimental data. Results indicate that comparisons are in good agreement with experimental data. The three-dimensional simulation improves the comparisons for over-expanded flow conditions as compared with two-dimensional assumptions.
A new porous condition has been implemented in the PAB3D solver for simulating the flow over porous surfaces. The newly-added boundary condition is utilized to compute the flow field of a non-axisymmetric, convergent-divergent nozzle incorporating porous cavities for shock-boundary layer interaction control. The nozzle has an expansion ratio (exit area/throat area) of 1.797 and a design nozzle pressure ratio of 8.78. The flow fields for a baseline nozzle (no porosity) and for a nozzle with porous surfaces (10% porosity ratio) are computed for NPR varying from 2.01 to 9.54. Computational model results indicate that the over-expanded nozzle flow was dominated by shock-induced boundary-layer separation. Porous configurations were capable of controlling off-design separation in the nozzle by encouraging stable separation of the exhaust flow. Computational simulation results, wall centerline pressure, mach contours, and thrust efficiency ratio are presented and discussed. Computed results are in excellent agreement with experimental data.
A new porous condition has been implemented in the PAB3D solver for simulating the flow over porous surfaces. The newly-added boundary condition is utilized to compute the flow field of a non-axisymmetric, convergent-divergent nozzle incorporating porous cavities for shock-boundary layer interaction control. The nozzle has an expansion ratio (exit area/throat area) of 1.797 and a design nozzle pressure ratio of 8.78. The flow fields for a baseline nozzle (no porosity) and for a nozzle with porous surfaces (10% porosity ratio) are computed for NPR varying from 2.01 to 9.54. Computational model results indicate that the over-expanded nozzle flow is dominated by shock-induced boundary-layer separation. Porous configurations are capable of controlling off-design separation in the nozzle by encouraging stable separation of the exhaust flow. Computational simulation results, wall centerline pressure, mach contours, and thrust efficiency ratio are presented and discussed. Computed results are in excellent agreement with experimental data.
Low thrust divergent flow cesium-on-tungsten contact ionization electrostatic thrustor for satellite attitude control and stationkeeping missions