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.
Nonsteady flow conditions in subsonic section of Laval nozzle to investigate nonlinear transverse combustion instability in liquid propellant rocket motors
Ignition characteristics of AN-F-32 and two AN-F-58a fuels were studied in a single can-type turbojet combustor under air-flow conditions representing engine speeds of 1600, 2500, and 4000 rpm, altitudes from sea level to 30,000 feet, ambient temperatures at sea level from 90 degrees to minus 36 degrees F, and flight Mach numbers of 0 and 0.6. Critical fuel-flow rates for ignition increased with increase in preignition engine speed, with increase in altitude, or with decrease in sea-level ambient temperature. This flow rate appears to increase in a direct relation to decrease in fuel volatility as indicated by the 10-percent-evaporated temperature.
An analysis of the flow state on a trapezoidal wing model from the NASA 3-D High Lift Flow Physics Experiment is presented. The objective of the experiment was to characterize the flow over a non-proprietary semi-span three-element high-lift configuration to aid in assessing the state of the art in the computation of three-dimensional high-lift flows. Surface pressures and hot-film sensors are used to determine the flow conditions on the slat, main, and flap. The locations of the attachments lines and the values of the attachment line Reynolds number are estimated based on the model surface pressures. Data from the hot-films are used to determine if the flow is laminar, transitional, or turbulent by examining the hot-film time histories, statistics, and frequency spectra.
Time/frequency analysis of fluctuations measured by pressure taps and strain gauges in the experimental studies of the flexible semispan model of a high-speed civil transport wing configuration is performed. The interest is in determining the coupling between the aerodynamic loads and structural motions that led to the hard flutter conditions and loss of the model. The results show that, away from the hard flutter point, the aerodynamic loads at all pressure taps near the wing tip and the structural motions contained the same frequency components. On the other hand, in the flow conditions leading to the hard flutter, the frequency content of the pressure fluctuations near the leading and trailing edges varied significantly. This led to contribution to the structural motions over two frequency ranges. The ratio of these ranges was near 2:1, which suggests the possibility of nonlinear structural coupling.
A novel approach of using counterflowing jets positioned strategically on the aircraft and exploiting its long penetration mode (LPM) of interaction towards sonic-boom mitigation forms the motivation for this study. Given that most previous studies on the counterflowing LPM jet have all been on blunt bodies and at high supersonic or hypersonic flow conditions, exploring the feasibility to obtain a LPM jet issuing from a slender body against low supersonic freestream conditions is the main focus of this study. Computational fluid dynamics computations of axisymmetric models (cone-cylinder and quartic geometry), of relevance to NASA's High Speed project, are carried out using the space-time conservation element solution element viscous flow solver with unstructured meshes. A systematic parametric study is conducted to determine the optimum combination of counterflowing jet size, mass flow rate, and nozzle geometry for obtaining LPM jets. Details from these computations will be used to assess the potential of the LPM counterflowing supersonic jet as a means of active flow control for enabling supersonic flight over land and to establish the knowledge base for possible future implementation of such technologies.
This report documents a test of aircraft exhaust nozzle designs applicable to commercial supersonic aircraft as envisioned in the 2025-2035 time frame. The test is primarily intended to provide data about the noise such aircraft would produce during landing and takeoff (LTO) operation, specifically the jet noise component of this noise. The tests were conducted at the NASA Glenn Research Center’s Aero-Acoustic Propulsion Laboratory in March 2020. The tests used the High-Flow Jet Exhaust Rig in the Nozzle Acoustic Test Rig which simulates the flow from a dual-stream turbofan engine in a flight stream. Seven exhaust configurations were tested with various combinations of internal mixer and external plug, and covered engine cycles having nozzle pressures from 1.5 to 2.3 at a flight speed of Mach 0.3. A few unheated flow conditions were also tested to tie results to other baseline test results. Data acquired and included in this report are 1) far-field acoustic spectra, 2) phased array measurements of noise source distributions, and 3) background-oriented schlieren of the flow near the nozzles. These data, along with the nozzle geometry in CAD format, are given as Appendices and are available online.
Compact heat exchangers are being examined for use in advanced, high-temperature nuclear or concentrated solar power plants to improve thermal efficiencies and reduce capital costs. At the time of this work, no experimental investigation of liquid sodium in the popular diffusion-bonded heat exchanger design had been performed. Herein this paper describes the thermal-hydraulic and overall operations performance observed during testing of a 316/L stainless steel “zig-zag” microchannel heat exchanger operating between a (0.5 kg/s) sodium loop and a (0.3 kg/s) recuperated supercritical CO 2 system with outlet conditions up to 550 °C at 16 MPa. Observed pressure drop is presented in the form of an effective friction factor correlation which aligns closely with recent experimental data. The thermal performance of this heat exchanger, measured by its overall conductance under various flow conditions, was decomposed to examine the independent convection behaviors of both the CO 2 and sodium streams. These results are presented in the form of Nusselt correlations fitted to the experimental data. The sodium heat transfer correlation aligns closely with the widely accepted Lubarsky-Kaufman correlation, while the CO 2 heat transfer correlation lies between previous correlations developed by computational and experimental methods. This was the first long-duration use of sodium in a microchannels (~1.8 mm hydraulic diameter) heat exchanger. The use of an oxide control system enabled trouble-free operation of the heat exchanger for over 150 h of high-temperature operation. This first-of-a-kind demonstration using liquid sodium and supercritical carbon dioxide shows great promise for further application of this technology in advanced nuclear or solar power plants.
Despite the extensive research on bubble collapse near rigid walls, the bubble collapse dynamics in the presence of shear flow near a rigid wall is poorly understood. We conduct direct simulations of the Navier–Stokes equations to explore the bubble dynamics and pressures during bubble collapse near a rigid, flat wall under linear shear flow conditions. We examine the dependence of the bubble collapse morphology and wall pressures on the initial bubble location and shear rate. We find that shear distorts the bubble, generating two re-entrant jets – one developing from the side opposite to the mean flow and the other from the far end toward the wall. Upon impact of the jet on the opposite side of the bubble, water-hammer shocks are produced, which propagate outward and interact with the convoluted bubble shape. The shock stretches the bubble towards the wall, resulting in a closer impact location for the jet originating from the far end compared with the case with no shear flow. The water-hammer pressure location can be approximated as the theoretical distance travelled by a particle initialised at the bubble centre with the corresponding constant shear flow velocity. The maximum wall pressures can thus be predicted by considering the distance between the far jet impingement location and the wall along the wall-normal direction. As the shear rate is increased, the maximum wall pressure increases, although only marginally. We determine the critical initial stand-off distance from the wall at which the bubble morphology is shear dominated, i.e. characterised by converging re-entrant jets.
A methodology for simulation of molecular mixing and the resulting velocity and temperature fields has been developed. The ideas are applied to the flow conditions present in the NASA Lewis Planar Reacting Shear Layer (PRSL) facility, and results compared to experimental data. A gaussian transverse turbulent velocity distribution is used in conjunction with a linearly increasing time scale to describe the mixing of different regions of the flow. Equilibrium reaction calculations are then performed on the mix to arrive at a new species composition and temperature. Velocities are determined through summation of momentum contributions. The analysis indicates a combustion efficiency of the order of 80 percent for the reacting mixing layer, and a turbulent Schmidt number of 2/3. The success of the model is attributed to the simulation of large-scale transport of fluid. The favorable comparison shows that a relatively quick and simple PC calculation is capable of simulating the basic flow structure in the reacting and non-reacting shear layer present in the facility given basic assumptions about turbulence properties.
The buffet response of the flexible twin-tail configuration of the F/A-18 and a generic F-111 aircraft are computationally simulated and experimentally validated. The problem is a multidisciplinary one which requires the sequential solution of three sets of equations on a multi-block grid structure. The first set is the unsteady, compressible, full Navier-Stokes equations. The second set is the aeroelastic equations for bending and torsional twin-tail responses. The third set is the grid-displacement equations which are used to update the grid coordinates due to the tail deflections. The computational models consist of a 76 deg. swept back, sharp edged delta wing of aspect ratio of one and a swept-back F/A-18 or F-117 twin-tail. The configuration is pitched at 30 deg. angle-of-attack. The problem is solved for the initial flow conditions with the twin tails kept rigid. Next, the aeroelastic equations of the tails are turned on along with the grid-displacement equations to solve for the bending and torsional tails responses due to the unsteady loads produced by the vortex breakdown flow of the leading-edge vortex cores of the delta wing. Several spanwise locations of the twin tails are investigated. The computational results are validated using several existing experimental data.
This study supports the NASA Langley sponsored project aimed at determining the viability of using Euler technology for preliminary design use. The primary objective of this study was to assess the accuracy and efficiency of the Boeing, St. Louis unstructured grid flow field analysis system, consisting of the MACGS grid generation and NASTD flow solver codes. Euler solutions about the Aero Configuration/Weapons Fighter Technology (ACWFT) 1204 aircraft configuration were generated. Several variations of the geometry were investigated including a standard wing, cambered wing, deflected elevon, and deflected body flap. A wide range of flow conditions, most of which were in the non-linear regimes of the flight envelope, including variations in speed (subsonic, transonic, supersonic), angles of attack, and sideslip were investigated. Several flowfield non-linearities were present in these solutions including shock waves, vortical flows and the resulting interactions. The accuracy of this method was evaluated by comparing solutions with test data and Navier-Stokes solutions. The ability to accurately predict lateral-directional characteristics and control effectiveness was investigated by computing solutions with sideslip, and with deflected control surfaces. Problem set up times and computational resource requirements were documented and used to evaluate the efficiency of this approach for use in the fast paced preliminary design environment.
Concerns regarding noise, propulsive efficiency, and fuel burn are inspiring aircraft designs wherein the propulsive turbomachines are partially (or fully) embedded within the airframe; such designs present serious concerns with regard to aerodynamic and aeromechanic performance of the compression system in response to inlet distortion. Previously, a preliminary design of a forward-swept high-speed fan exhibited flutter concerns in clean-inlet flows, and the present author then studied this fan further in the presence of off-design distorted in-flows. Continuing this research, a three-dimensional, unsteady, Navier-Stokes computational fluid dynamics code is again applied to analyze and corroborate fan performance with clean inlet flow and now with a simplified, sinusoidal distortion of total pressure at the aerodynamic interface plane. This code, already validated in its application to assess aerodynamic damping of vibrating blades at various flow conditions using a one-way coupled energy-exchange approach, is modified to include a two-way coupled timemarching aeroelastic simulation capability. The two coupling methods are compared in their evaluation of flutter stability in the presence of distorted in-flows.
A computational framework for the simulation of dynamic stall suppression with combustion-powered actuation (COMPACT) is validated against wind tunnel experimental results on a VR-12 airfoil. COMPACT slots are located at 10% chord from the leading edge of the airfoil and directed tangentially along the suction-side surface. Helicopter rotor-relevant flow conditions are used in the study. A computationally efficient two-dimensional approach, based on unsteady Reynolds-averaged Navier-Stokes (RANS), is compared in detail against the baseline and the modified airfoil with COMPACT, using aerodynamic forces, pressure profiles, and flow-field data. The two-dimensional RANS approach predicts baseline static and dynamic stall very well. Most of the differences between the computational and experimental results are within two standard deviations of the experimental data. The current framework demonstrates an ability to predict COMPACT efficacy across the experimental dataset. Enhanced aerodynamic lift on the downstroke of the pitching cycle due to COMPACT is well predicted, and the cycleaveraged lift enhancement computed is within 3% of the test data. Differences with experimental data are discussed with a focus on three-dimensional features not included in the simulations and the limited computational model for COMPACT.
Shock-wave boundary-layer interactions (SWBLI) pose unique challenges for the slender airframes and sharp control surfaces characteristic of high-speed airbreathing aircraft. Spikes in pressure and thermal loading have been shown to occur in SWBLI regions; thus, predicting their location and behavior proves critical for assuring aircraft structural integrity and fatigue resilience. SWBLIs frequently lead to severe turbulent boundary-layer separation and corresponding unsteady motions characterized by a wide range of frequencies. Low-frequency motion has been shown to couple with the structure of a compliant panel, and turbulent boundary layers thus present a particular danger for thin control surfaces. Previous studies have shown that the thermal loading on a compliant panel in a ramp-induced SWBLI may have a signi cant impact on the panel vibrations. A multitude of studies has investigated the heating on rigid structures, but a charactrization of the heating effects on compliant structures is necessary to fully understand the fluid-thermal-structural interactions (FTSI) occurring in ramp-induced SWBLIs. This studuses infrared thermography and high-frequency Kulite pressure transducers to examine the heating loads on rigid and compliant ramps at various angles and flow conditions at Mach 6.
This report documents the progress in preparing the ANSYS/Fluent ® and MELCOR models to perform characterization studies of aerosols dispersal and deposition within a select spent fuel cask system. Steady state thermal-hydraulic cask response is modeled with both codes, at present, while the MELCOR source code is being modified to allow imposed thermal-hydraulic conditions with aerosol physics calculations. This will allow the MELCOR model to assume the thermal-hydraulic calculation from the ANSYS/Fluent ® , while only computing the aerosol physics. Detailed results are presented on the thermal-hydraulic analysis of the MAGNASTOR ® cask for the current ANSYS/Fluent ® model, with convergent conditions observed over two fidelities. While the MELCOR computation computes the steady state conditions, they differ sufficiently from that computed with ANSYS/Fluent ® . The MELCOR analysis include a set of UO 2 sources to investigate system response for interim reporting. The airborne concentrations and evolving distributions are presented. Model development is anticipated to continue as additional components impacting flow conditions or aerosol deposition within the MAGNASTOR ® cask and Westinghouse fuel assemblies are identified.
A distributed parameter model for computing unsteady flow conditions in liquid systems is presented. The model is adapted for use in conjunction with a digital computer. The analytical method employed is essentially a solution synthesized from the effects of incremental step pressure pulses. The pressure pulses are generated because of incremental flow-rate changes that originate in a fluid system from a variety of sources, including the mechanical motion of the system structure. The pressure pulses propagate throughout the system at sonic velocity and are partially transmitted and reflected at each discontinuity. The velocity change caused by each pressure pulse is obtained from the relation for the characteristic acoustic impedance. Pressure and velocity time histories at any point in the system are obtained by a timewise summation of the contributions of the incremental pressure pulses passing that point. The analysis is presented in a form general enough to be applied to a variety of liquid-filled fluid systems. To illustrate the application of the method to a specific system, the response of a long, straight propellant line to a sinusoidal inlet flow and pressure perturbation is computed. Both a constant-cross-section line and a tapered line are analyzed in the example.
Experimental data are presented for the growth of thrombi (blood clots) in a stagnation point flow of fresh blood. Thrombus shape, size and structure are shown to depend on local flow conditions. The evolution of a thrombus is described in terms of a physical model that includes platelet diffusion, a platelet aggregation mechanism, and diffusion and convection of the chemical species responsible for aggregation. Diffusion-controlled and convection-controlled regimes are defined by flow parameters and thrombus location, and the characteristic growth pattern in each regime is explained. Quantitative comparisons with an approximate theoretical model are presented, and a more general model is formulated.
The influence of wall slip and catalytic atom-recombination on the flow field and wall heat flux are calculated for high altitude flight and arc jet flow conditions. Boundary equations, which include velocity slip, temperature jump, and wall catalytic atom recombination, are coupled to the viscous reacting multicomponent Navier-Stokes equations. These equations are solved using a time-dependent finite difference technique applied to spheres in an arc jet flow (Reynolds number of 550) and a high altitude flight case representative of the Space Shuttle Orbiter (Reynolds number of 450). The results indicate that catalysis strongly influences the temperature jump, but not the velocity slip. Slip increases the atom fraction and temperature at both the wall and the flow field. Likewise, the shock stand-off distance, the wall heat flux, and friction coefficient are increased over the nonslip cases. The reacting gas calculations confirm the chemically frozen nature of the shock layer in arc jet flows.