SEARCH · Engineering Papers
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.
Report presents computational study of accretion of ice on wings and on horizontal tail surfaces of simplified version of Boeing 737 airplane. Calculations performed by use of LEWICE3D computer code, which incorporates panel-method mathematical model of three-dimensional potential flow and mathematical models of trajectories of drops of water entrained in flow, impingement of drops on surfaces, and transfer of heat.
Report presents mathematical model for use in designing rocket-launching stand. Predicts pressure and thermal environment, as well as thermal responses of structures to impinging rocket-exhaust plumes. Enables relatively inexperienced analyst to determine time-varying distributions and absolute levels of pressure and heat loads on structures.
An experimental set-up to measure low-energy (below 1 keV) sputtering of materials is described. The materials to be bombarded represent ion thruster components as well as insulators used in the stationary plasma thruster. The sputtering takes place in a 9 inch diameter spherical vacuum chamber. Ions of argon, krypton and xenon are used to bombard the target materials. The sputtered neutral atoms are detected by a secondary neutral mass spectrometer (SNMS). Samples of copper, nickel, aluminum, silver and molybdenum are being sputtered initially to calibrate the spectrometer. The base pressure of the chamber is approximately 2 x 10(exp -9) Torr. the primary ion beam is generated by an ion gun which is capable of delivering ion currents in the range of 20 to 500 nA. The ion beam can be focused to a size approximately 1 mm in diameter. The mass spectrometer is positioned 10 mm from the target and at 90 deg angle to the primary ion beam direction. The ion beam impinges on the target at 45 deg. For sputtering of insulators, charge neutralization is performed by flooding the sample with electrons generated from an electron gun. Preliminary sputtering results, methods of calculating the instrument response function of the spectrometer and the relative sensitivity factors of the sputtered elements will be discussed.
A system for and method of detecting and measuring concentrations of an ultrasonically-reflective microsphere contrast agent involving detecting non-linear sum and difference beat frequencies produced by the microspheres when two impinging signals with non-identical frequencies are combined by mixing. These beat frequencies can be used for a variety of applications such as detecting the presence of and measuring the flow rates of biological fluids and industrial liquids, including determining the concentration level of microspheres in the myocardium.
A noncontact pyroelectric infrared detector is described. A pyroelectric film that also has piezoelectric properties is held in place so that it is free to vibrate. It is electrically stimulated to vibrate at a resonance frequency. The vibrating film forms part of a balanced bridge circuit. As thermal radiation impinges on the film the pyroelectric effect causes the resonance frequency to change, thereby unbalancing the bridge circuit. A differential amplifier tracks the change in voltage across the bridge. The resulting voltage signal is further processed by a bandpass filter and a precision rectifier. The device allows for DC or static temperature measurements without the use of a mechanical chopping device.
The computational fluid dynamics code, PARC3D, is tested to see if its use of non-physical artificial dissipation affects the accuracy of its results. This is accomplished by simulating a shock-laminar boundary layer interaction and several hypersonic flight conditions of the Pegasus(TM) launch vehicle using full artificial dissipation, low artificial dissipation, and the Engquist filter. Before the filter is applied to the PARC3D code, it is validated in one-dimensional and two-dimensional form in a MacCormack scheme against the Riemann and convergent duct problem. For this explicit scheme, the filter shows great improvements in accuracy and computational time as opposed to the nonfiltered solutions. However, for the implicit PARC3D code it is found that the best estimate of the Pegasus experimental heat fluxes and surface pressures is the simulation utilizing low artificial dissipation and no filter. The filter does improve accuracy over the artificially dissipative case but at a computational expense greater than that achieved by the low artificial dissipation case which has no computational time penalty and shows better results. For the shock-boundary layer simulation, the filter does well in terms of accuracy for a strong impingement shock but not as well for weaker shock strengths. Furthermore, for the latter problem the filter reduces the required computational time to convergence by 18.7 percent.
The results of further spectroscopic studies on the plume from a 3 cm ion source operated on an argon propellant is reported on. In particular, it is shown that it should be possible to use the spectroscopic technique to measure the plasma density of the ion plume close to the grids, where it is difficult to use electrical probe measurements. How the technique, along with electrical probe measurements in the far downstream region of the plume, can be used to characterize the operation of a three-grid, 15 cm diameter thruster from NASA JPL is outlined. Pumping speed measurements on the Vacuum Research Facility have shown that this facility should be adequate for testing the JPL thruster at pressures in the low 10(exp -5) Torr range. Finally, we describe a simple analytical model which can be used to calculate the grid impingement current which results from charge-exchange collisions in the ion plume.
A predictive tool for hypergolic bipropellant spray combustion and flow evolution in the OMV VTE (orbital maneuvering vehicle variable thrust engine) is described. It encompasses a computational technique for the gas phase governing equations, a discrete particle method for liquid bipropellant sprays, and constitutive models for combustion chemistry, interphase exchanges, and unlike impinging liquid hypergolic stream interactions. Emphasis is placed on the phenomenological modelling of the hypergolic liquid bipropellant gasification processes. An application to the OMV VTE combustion chamber is given in order to show some of the capabilities and inadequacies of this tool.
One source contributing to wind tunnel background noise is microphone self-noise. An experiment was conducted to investigate the flow-induced acoustic oscillations of Bruel & Kjaer (B&K) in-flow microphones. The results strongly suggest the B&K microphone cavity behaves more like an open cavity. Their cavity acoustic oscillations are likely caused by strong interactions between the cavity shear layer and the cavity trailing edge. But the results also suggest that cavity shear layer oscillations could be coupled with cavity acoustic resonance to generate tones. Detailed flow velocity measurements over the cavity screen have shown inflection points in the mean velocity profiles and high disturbance and spectral intensities in the vicinity of the cavity trailing edge. These results are the evidence for strong interactions between cavity shear layer oscillations and the cavity trailing edge. They also suggest that beside acoustic signals, the microphone inside the cavity has likely recorded hydrodynamic pressure oscillations, too. The results also suggest that the forebody shape does not have a direct effect on cavity oscillations. For the FITE (Flow Induced Tone Eliminator) microphone, it is probably the forebody length and the resulting boundary layer turbulence that have made it work. Turbulence might have thickened the boundary layer at the separation point, weakened the shear layer vortices, or lifted them to miss impinging on the cavity trailing edge. In addition, the study shows that the cavity screen can modulate the oscillation frequency but not the cavity acoustic oscillation mechanisms.
The pressure-sensitive paint method was used in the test of a high-sweep oblique wing model, conducted in the NASA Ames 9- by 7-ft Supersonic Wind Tunnel. Surface pressure data was acquired from both the luminescent paint and conventional pressure taps at Mach numbers between M = 1.6 and 2.0. In addition, schlieren photographs of the outer flow were used to determine the location of shock waves impinging on the model. The results show that the luminescent pressure-sensitive paint can capture both global and fine features of the static surface pressure field. Comparison with conventional pressure tap data shows good agreement between the two techniques, and that the luminescent paint data can be used to make quantitative measurements of the pressure changes over the model surface. The experiment also demonstrates the practical considerations and limitations that arise in the application of this technique under supersonic flow conditions in large-scale facilities, as well as the directions in which future research is necessary in order to make this technique a more practical wind-tunnel testing tool.
An extensive experimental program was initiated for the purpose of understanding the mechanisms leading to bubble generation during fluid handling procedures in a microgravity environment. Several key fluid handling procedures typical for PCG experiments were identified for analysis in that program. Experiments were designed to specifically understand how such procedures can lead to bubble formation. The experiments were then conducted aboard the NASA KC-135 aircraft which is capable of simulating a low gravity environment by executing a parabolic flight attitude. However, such a flight attitude can only provide a low gravity environment of approximately 10-2go for a maximum period of 30 seconds. Thus all of the tests conducted for these experiments were designed to last no longer than 20 seconds. Several experiments were designed to simulate some of the more relevant fluid handling procedures during protein crystal growth experiments. These include submerged liquid jet cavitation, filling of a cubical vessel, submerged surface scratch, attached drop growth, liquid jet impingement, and geysering experiments. To date, four separate KC-135 flight campaigns were undertaken specifically for performing these experiments. However, different experiments were performed on different flights.
The final figuring step in the fabrication of an optical component involves imparting a specified contour onto the surface. This can be expensive and time consuming step. The recent development of ion beam figuring provides a method for performing the figuring process with advantages over standard mechanical methods. Ion figuring has proven effective in figuring large optical components. The process of ion beam figuring removes material by transferring kinetic energy from impinging neutral particles. The process utilizes a Kaufman type ion source, where a plasma is generated in a discharge chamber by controlled electric potentials. Charged grids extract and accelerate ions from the chamber. The accelerated ions form a directional beam. A neutralizer outside the accelerator grids supplies electrons to the positive ion beam. It is necessary to neutralize the beam to prevent charging workpieces and to avoid bending the beam with extraneous electro-magnetic fields. When the directed beam strikes the workpiece, material sputters in a predicable manner. The amount and distribution of material sputtered is a function of the energy of the beam, material of the component, distance from the workpiece, and angle of incidence of the beam. The figuring method described here assumes a constant beam removal, so that the process can be represented by a convolution operation. A fixed beam energy maintains a constant sputtering rate. This temporally and spatially stable beam is held perpendicular to the workpiece at a fixed distance. For non-constant removal, corrections would be required to model the process as a convolution operation. Specific figures (contours) are achieved by rastering the beam over the workpiece at varying velocities. A unique deconvolution is performed, using series-derivative solution developed for the system, to determine these velocities.
An experimental study is conducted to determine the dependence of vortex generator geometry and impinging flow conditions on shed vortex circulation and crossplane peak vorticity for one type of vortex generator. The vortex generator is a symmetric airfoil having a NACA 0012 cross-sectional profile. The geometry and flow parameters varied include angle-of-attack alfa, chordlength c, span h, and Mach number M. The vortex generators are mounted either in isolation or in a symmetric counter-rotating array configuration on the inside surface of a straight pipe. The turbulent boundary layer thickness to pipe radius ratio is delta/R = 0. 17. Circulation and peak vorticity data are derived from crossplane velocity measurements conducted at or about 1 chord downstream of the vortex generator trailing edge. Shed vortex circulation is observed to be proportional to M, alfa, and h/delta. With these parameters held constant, circulation is observed to fall off in monotonic fashion with increasing airfoil aspect ratio AR. Shed vortex peak vorticity is also observed to be proportional to M, alfa, and h/delta. Unlike circulation, however, peak vorticity is observed to increase with increasing aspect ratio, reaching a peak value at AR approx. 2.0 before falling off.
A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.
During the past year, various computational fluid dynamic (CFD) analyses were performed at Marshall Space Flight Center to support the Reusable Solid Rocket Motor program. The successful completion of these analyses involved application of the CFD codes FDNS and CELMINT. The topics addressed by the analyses were: (1) the design and prediction of slag pool accumulation within the five inch test motor, (2) prediction of slag pool behavior and its response to lateral accelerations, (3) the clogging of potential insulation debonds within the nozzle by slag accumulation, (4) the behavior of jets within small voids inside nozzle joint gaps, (5) The effect of increased inhibitor stiffness on motor acoustics, and (6) the effect of a nozzle defect on particle impingement enhanced erosion. The emphasis of this presentation will be to further discuss the work in topics 3, 4, and 5.
Recent trends towards an increase in the stiffness of the acrylonitrile butadiene rubber (NBR) insulation material used in the construction of the redesigned solid rocket motor (RSRM) propellant inhibitors prompted questions about possible effects on RSRM performance. The specific objectives of the computational fluid dynamics (CFD) task included: (1) the definition of pressure loads to calculate the deformed shape of stiffer inhibitors, (2) the calculation of higher port velocities over the inhibitors to determine shifts in the vortex shedding or edge tone frequencies, and (3) the quantification of higher slag impingement and collection rates on the inhibitors and in the submerged nose nozzle cavity.
This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.
Pressure sensitive paint is a newly-developed optical measurement technique with which one can get a continuous pressure distribution in much shorter time and lower cost than a conventional pressure tap measurement. However, most of the current pressure sensitive paint applications are restricted to steady pressure measurement at high speeds because of the small signal-to-noise ratio at low speed and a slow response to pressure changes. In the present study, three phases of work have been completed to extend the application of the pressure sensitive paint technique to low-speed testing and to investigate the applicability of the paint technique to unsteady flow. First the measurement system using a commercially available PtOEP/GP-197 pressure sensitive paint was established and applied to impinging jet measurements. An in-situ calibration using only five pressure tap data points was applied and the results showed good repeatability and good agreement with conventional pressure tap measurements on the whole painted area. The overall measurement accuracy in these experiments was found to be within 0.1 psi. The pressure sensitive paint technique was then applied to low-speed wind tunnel tests using a 60 deg delta wing model with leading edge blowing slots. The technical problems encountered in low-speed testing were resolved by using a high grade CCD camera and applying corrections to improve the measurement accuracy. Even at 35 m/s, the paint data not only agreed well with conventional pressure tap measurements but also clearly showed the suction region generated by the leading edge vortices. The vortex breakdown was also detected at alpha=30 deg. It was found that a pressure difference of 0.2 psi was required for a quantitative pressure measurement in this experiment and that temperature control or a parallel temperature measurement is necessary if thermal uniformity does not hold on the model. Finally, the pressure sensitive paint was applied to a periodically changing pressure field with a 12.8s time period. A simple first-order pole model was applied to deal with the phase lag of the paint. The unsteady pressure estimated from the time-changing pressure sensitive paint data agreed well with the pressure transducer data in regions of higher pressure and showed the possibility of extending the technique to unsteady pressure measurements. However, the model still needs further refinement based on the physics of the oxygen diffusion into the paint layer and the oxygen quenching on the paint luminescence.