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.
The boundary condition on the acoustic perturbation velocity at an impermeable surface in a flow is considered for the cases in which the surface generates a sound field by vibration or is acoustically deformed by an incident sound field. It is shown that in general the condition is not equivalent to the requirement of continuity of acoustic particle displacement in the direction normal to the unperturbed surface.
The aerodynamic characteristics of an 11 percent thick symmetrical supercritical airfoil were determined in an 8 foot transonic pressure tunnel over a Mach number range of 0.30 to 0.85. The Reynolds number of the tests, based on the airfoil chord, varied with Mach number over a range of 360 million to 774 million. The geometric angle of attack varied from -0.5 to 10.5 deg. The abrupt drag rise for the supercritical airfoil at zero normal force conditions occurs at a Mach number just above 0.82. The corresponding drag rise Mach number for a conventional NACA 0012 airfoil is approximately 0.70. At zero normal force conditions, the level of supervelocity over the supercritical airfoil is considerably reduced from that for the NACA 0012 airfoil. Also, the shock wave for the supercritical airfoil is substantially weaker than that for the NACA 0012 airfoil. For a Mach number of 0.82 and zero normal force, the flow over the present airfoil is supercritical; however, there is no discernible shock wave in the flow, indicating near isentropic recompression.
Thermal convection phenomena in a vertically oriented Bridgman-Stockbarger apparatus were modeled by computer simulations for different gravity conditions, ranging from earth conditions to extremely low gravity, approximate space conditions. The modeling results were obtained by the application of a state-of-the art, transient, multi-dimensional, completely densimetrically coupled, discrete-element computational model which was specifically developed for the simulation of flow, temperature, and species concentration conditions in two-phase (solid-liquid) systems. The computational model was applied to the simulation of the flow and the thermal conditions associated with the convection phenomena in a modified Germanium-Silicon charge enclosed in a stationary fused-silica ampoule. The results clearly indicated that the gravitational field strength influences the characteristics of the coherent vortical flow patterns, interface shape and position, maximum melt velocity, and interfacial normal temperature gradient.
The linearized equation of motion is given in a Lagrangian representation for a rotating plasma with anisotropic pressure. A WKB theory is developed for large-n ballooning modes in an axisymmetric configuration with field-aligned and rigid toroidal flows. In the presence of field-aligned flows, it is shown that a resonance occurs which is strongly suggestive of a generalized mirror instability. In the presence of toroidal rotation, a possible stabilizing effect is identified for P(normal) greater than P(parallel). Finally, as a special case of the theory, the necessary and sufficient conditions for stability in a static, anisotropic plasma are obtained.
Modifications have been proposed to effect further improvement of the device described in Improved Piezoelectrically Actuated Microvalve (NPO-30158), NASA Tech Briefs, Vol. 26, No. 1 (January 2002), page 29. To recapitulate: What is being developed is a prototype of valves for microfluidic systems and other microelectromechanical systems (MEMS). The version of the valve reported in the cited previous article included a base (which contained a seat, an inlet, and an outlet), a diaphragm, and a linear actuator. With the exception of the actuator, the parts were micromachined from silicon. The linear actuator consisted of a stack of piezoelectric disks in a rigid housing. To make the diaphragm apply a large sealing force on the inlet and outlet, the piezoelectric stack was compressed into a slightly contracted condition during assembly of the valve. Application of a voltage across the stack caused the stack to contract into an even more compressed condition, lifting the diaphragm away from the seat, thereby creating a narrow channel between the inlet and outlet. The positions of the inlet and outlet, relative to the diaphragm and seat, were such that the inlet flow and pressure contributed to sealing and thus to a desired normally-closed mode of operation.
Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.
The Flow Boiling and Condensation Experiment (FBCE) Transfer Line (TL) test module is being developed to study fluid physics of cryogenic transfer lines in microgravity. The TL module will reside in the FBCE facility located within the Fluids Integration Rack (FIR) aboard the International Space Station (ISS). TL will use the same normal perfluoro hexane (nPFH) used in the previous FBCE test sections. To simulate cryogenic line transfer in the visualization section, a sapphire tube coated with indium tin oxide (ITO) is heated via Joule heated while a bypass line is used to flow liquid nPFH at the desired inlet conditions. Once the visualization line reaches steady-state temperature and pressure conditions, flow is diverted through the heated line and chilldown begins. All flow regimes observed within a cryogenic TL are created and studied using visible imaging with diffuse backlighting to image inverted annular flow, transition, and nucleate boiling. A visible laser light sheet provided by FIR will provide a means to observe and track droplets in the dispersed flow regime. IR imaging of the visualization tube outer surface with the emissivity of the ITO coating provides the surface chilldown temperature profile. Preliminary imaging system development includes some of the components used in the Flow Boiling Module (FBM) and some new components that are easily integrated into the PC-based system but will require some refinement to operate properly on the FIR PCIe/104 imaging system. Optical design and packaging within the limited volume of the module enclosure is also presented. This paper presents FBCE-TL visible and infrared laboratory-based imaging studies and anticipated visualization capabilities of the TL module under preliminary breadboard development and under the utilization constraints within the FBCE/FIR host facility.
Nuclei are nearly transparent to each other when they collide at high energy, but the collisions do produce high-energy density matter in the central rapidity region where most experimental measurements are made. What happens to the receding nuclear fireballs? We calculate the energy loss of the nuclei using the color glass condensate model. We then use a simple space-time picture of the collision to calculate the baryon and energy densities of the receding fireballs. For central collisions of large nuclei at the BNL Relativistic Heavy Ion Collider and the CERN Large Hadron Collider we find baryon densities more than ten times that of normal nuclear matter. Furthermore, these results provide initial conditions for subsequent hydrodynamic evolution and could test the equation of state at very high baryon densities.
High-fidelity analysis has been utilized in the design of beam target options for an accelerator driven subcritical system. Designs featuring stacks of plates with square cross section have been investigated for both tungsten and uranium target materials. The presented work includes the first thermal-hydraulic simulations of the full, detailed target geometry. The innovative target cooling manifold design features many regions with complex flow features, including 90 bends and merging jets, which necessitate three-dimensional fluid simulations. These were performed using the commercial computational fluid dynamics code STAR-CCM+. Conjugate heat transfer was modeled between the plates, cladding, manifold structure, and fluid. Steady-state simulations were performed but lacked good residual convergence. Unsteady simulations were then performed, which converged well and demonstrated that flow instability existed in the lower portion of the manifold. It was established that the flow instability had little effect on the peak plate temperatures, which were well below the melting point. The estimated plate surface temperatures and target region pressure were shown to provide sufficient margin to subcooled boiling for standard operating conditions. This demonstrated the safety of both potential target configurations during normal operation.
The objective of this program was to investigate analytically and experimentally the performance of heat pipes with composite wicks--specifically, those having pedestal arteries and screwthread circumferential grooves. An analytical model was developed to describe the effects of screwthreads and screen secondary wicks on the transport capability of the artery. The model describes the hydrodynamics of the circumferential flow in triangular grooves with azimuthally varying capillary menisci and liquid cross-sections. Normalized results were obtained which give the influence of evaporator heat flux on the axial heat transport capability of the arterial wick. In order to evaluate the priming behavior of composite wicks under actual load conditions, an 'inverted' glass heat pipe was designed and constructed. The results obtained from the analysis and from the tests with the glass heat pipe were applied to the OAO-C Level 5 heat pipe, and an improved correlation between predicted and measured evaporator and transport performance were obtained.
The interaction between a shock wave and a supersonic streamwise vortex is a fundamental fluid-dynamics problem with numerous practical applications. This paper describes an experimental study of this phenomenon. In particular, supersonic streamwise vortices of varying strength and Mach number were generated and measured using five-hole and total-temperature probes. In addition, the interactions between a vortex and either an oblique or a normal shock wave were visualized using schlieren and planar-laser-scattering techniques. The mean-flow measurements show both similarities and differences between the supersonic streamwise vortex and its incompressible counterpart, while the flow-visualization results show that the shock/vortex interaction is always unsteady and that, under certain conditions, the vortex can burst. The conditions necessary for supersonic vortex breakdown are presented.
The objective of this work was to enhance the predictive capability of widely used computational fluid dynamic (CFD) codes through the use of solution adaptive gridding. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. In order to study the accuracy and efficiency improvements due to the grid adaptation, it is necessary to quantify grid size and distribution requirements as well as computational times of non-adapted solutions. Flow fields about launch vehicles of practical interest often involve supersonic freestream conditions at angle of attack exhibiting large scale separate vortical flow, vortex-vortex and vortex-surface interactions, separated shear layers and multiple shocks of different intensity. In this work, a weight function and an associated mesh redistribution procedure is presented which detects and resolves these features without user intervention. Particular emphasis has been placed upon accurate resolution of expansion regions and boundary layers. Flow past a wedge at Mach=2.0 is used to illustrate the enhanced detection capabilities of this newly developed weight function.
A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.
Volatile Removal Assembly (VRA) is a subsystem of the Closed Environment Life Support System (CELSS) installed in the International Space Station. It is used for removing contaminants (volatile organics) in the wastewater produced by the space station crews. The major contaminants are formic acid, ethanol, and propylene glycol. The VRA contains a slim packbed reactor (3.5 cm diameter and four 28 cm long tubes in series) to perform catalyst oxidation of wastewater at elevated pressure and temperature under microgravity conditions. In the reactor, the contaminants are burned with oxygen gas (O2) to form water and carbon dioxide (CO2) that dissolves in the water stream. Optimal design of the reactor requires a thorough understanding about how the reactor performs under microgravity conditions. The objective of this study was to develop a mathematical model to interpret experimental data obtained from normal and microgravity conditions, and to predict the performance of VRA reactor under microgravity conditions. Catalyst oxidation kinetics and the total oxygen-water contact area control the efficiency of catalyst oxidation for mass transfer, which depends on oxygen gas holdup and distribution in the reactor. The process involves bubbly flow in porous media with chemical reactions in microgravity environment. This presents a unique problem in fluid dynamics that has not been studied. Guo et al. (2004) developed a mathematical model that predicts oxygen holdup in the VRA reactor. No mathematical model has been found in the literature that can be used to predict the efficiency of catalyst oxidation under microgravity conditions.
A new technique is reported for calculating the entire flow field on spherically blunted cones at high angles of attack and high laminar Reynolds numbers. An approximate system of parabolic equations obtained from the steady Navier-Stokes equations by assuming the viscous, streamwise derivative terms are small compared to the viscous normal and circumferential derivatives is the basis of the calculations. These equations are valid for both the inviscid and viscous regions, including the circumferential separation zone that develops on the leeward side at high angles of attack. Two different methods are used to obtain the initial conditions for these equations at the sphere cone tangency plane. For small nose Reynolds numbers, an axisymmetric merged layer solution around a sphere is rotated to provide a three-dimensional initial plane of data. For large nose Reynolds numbers, the nose region is solved using an inviscid, three dimensional time dependent solution combined with a boundary layer solution for the viscous flow. The computed flowfield including the leeward separation region is described and compared with data for a 7 deg half angle cone at 10 deg angle of attack, and a blunt 15 deg half angle cone at 15 deg angle of attack.
To support High-Temperature Gas-Cooled Reactor designs, the effects of corrosion and wear on a SS316L telescopic control rod prototype were investigated in helium at 775 °C with up to 5 ppm oxygen for durations up to 500 h. Three deployment conditions were tested to assess their impact on wear and oxide formation: severe cycles, daily cycles, and gravity-driven deployments. Here, a chromium-rich oxide layer developed on regions exposed to the helium flow, reaching up to 2 µm in the most affected areas. No carburization or decarburization was observed. Oxide breakthrough occurred mainly at grooves and contact points during severe cycling of the telescopic rod, while in the other tests it was limited to contact points. Friction forces increased due to oxide growth and thermal expansion, the latter being the dominant factor, but never enough to hinder rod deployment significantly. Coaxial misalignment of stages was identified as the main potential risk to reliable operation under normal conditions.
The panel flutter phenomenon is studied assuming free-molecule flow. This kind of analysis is relevant in the case of hypersonic flight vehicles traveling at high altitudes, especially in the leeward portion of the vehicle. In these conditions the aerodynamic shear can be expected to be considerably larger than the pressure at a given point, so that the effects of such a loading are incorporated into the structural model. Both the pressure and shear loadings are functions of the panel temperature, which can lead to great variations on the location of the stability boundaries for parametric studies. Different locations can, however, be 'collapsed' onto one another by using as ordinate an appropriately normalized dynamic pressure parameter. This procedure works better for higher values of the panel temperature for a fixed undisturbed flow temperature. Finally, the behavior of the system is studied when the panel has some initial curvature. This leads to the conclusion that it may be unrealistic to try to distinguish between a parabolic or sinusoidal initial shape.
Measurements of mass flux are obtained in a vitiated supersonic ground test facility using a sensor based on line-of-sight (LOS) diode laser absorption of water vapor. Mass flux is determined from the product of measured velocity and density. The relative Doppler shift of an absorption transition for beams directed upstream and downstream in the flow is used to measure velocity. Temperature is determined from the ratio of absorption signals of two transitions (lambda(sub 1)=1349 nm and lambda(sub 2)=1341.5 nm) and is coupled with a facility pressure measurement to obtain density. The sensor exploits wavelength-modulation spectroscopy with second-harmonic detection (WMS-2f) for large signal-to-noise ratios and normalization with the 1f signal for rejection of non-absorption related transmission fluctuations. The sensor line-of-sight is translated both vertically and horizontally across the test section for spatially-resolved measurements. Time-resolved measurements of mass flux are used to assess the stability of flow conditions produced by the facility. Measurements of mass flux are within 1.5% of the value obtained using a facility predictive code. The distortion of the WMS lineshape caused by boundary layers along the laser line-of-sight is examined and the subsequent effect on the measured velocity is discussed. A method for correcting measured velocities for flow non-uniformities is introduced and application of this correction brings measured velocities within 4 m/s of the predicted value in a 1630 m/s flow.