Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow dynamics”

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 1,135 records · Page 63

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗

Similarity laws of lunar and terrestrial volcanic flows

A mathematical model of a one dimensional, steady duct flow of a mixture of a gas and small solid particles (rock) was analyzed and applied to the lunar and the terrestrial volcanic flows under geometrically and dynamically similar conditions. Numerical results for the equilibrium two phase flows of lunar and terrestrial volcanoes under similar conditions are presented. The study indicates that: (1) the lunar crater is much larger than the corresponding terrestrial crater; (2) the exit velocity from the lunar volcanic flow may be higher than the lunar escape velocity but the exit velocity of terrestrial volcanic flow is much less than that of the lunar case; and (3) the thermal effects on the lunar volcanic flow are much larger than those of the terrestrial case.

Pai, S. I.↗

Characterization of medium and high flow series 1500 LARC-TPI

An imidized version of LARC-TPI was developed on NASA's request for use by both aerospace and electronic industries. Two batches of LARC-TPI 1500 Series, consisting of medium-flow and high-flow pellets with different inherent viscosities, are characterized for their dynamic and isothermal melt flow properties at 340 C, mechanical properties of extrudates, extrudate porosity, and degree of crystallinity. The Series shows significantly better melt flow properties than previosly examimed forms of this polymer.

Burks, H. D.↗

NASA. Marshall Space Flight Center Hydrostatic Bearing Activities

The basic approach for analyzing hydrostatic bearing flows at the Marshall Space Flight Center (MSFC) is briefly discussed. The Hydrostatic Bearing Team has responsibility for assessing and evaluating flow codes; evaluating friction, ignition, and galling effects; evaluating wear; and performing tests. The Office of Aerospace and Exploration Technology Turbomachinery Seals Tasks consist of tests and analysis. The MSFC in-house analyses utilize one-dimensional bulk-flow codes. Computational fluid dynamics (CFD) analysis is used to enhance understanding of bearing flow physics or to perform parametric analysis that are outside the bulk flow database. As long as the bulk flow codes are accurate enough for most needs, they will be utilized accordingly and will be supported by CFD analysis on an as-needed basis.

Benjamin, Theodore G.↗

Nuclear Thermal Propulsion Ground Test History

Nuclear Thermal Propulsion (NTP) was started in ~1955 under the Atomic Energy Commission as project Rover and was assigned to Los Alamos National Laboratory. The Nevada Test Site was selected in 1956 and facility construction began in 1957. The KIWI-A was tested on July 1, 1959 for 5 minutes at 70MW. KIWI-A1 was tested on July 8, 1960 for 6 minutes at 85MW. KIWI-A3 was tested on October 10, 1960 for 5 minutes at 100MW. The National Aeronautics and Space Administration (NASA) was formed in 1958. On August 31, 1960 the AEC and NASA established the Space Nuclear Propulsion Office and named Harold Finger as Director. Immediately following the formation of SNPO, contracts were awarded for the Reactor In Flight Test (RIFT), master plan for the Nuclear Rocket Engine Development Station (NRDS), and the Nuclear Engine for Rocket Vehicle Application (NERVA). From December 7, 1961 to November 30, 1962, the KIWI-B1A, KIWI-B1B, and KIWI-B4A were tested at test cell A. The last two engines were only tested for several seconds before noticeable failure of the fuel elements. Harold Finger called a stop to any further hot fire testing until the problem was well understood. The KIWI-B4A cold flow test showed the problem to be related to fluid dynamics of hydrogen interstitial flow causing fuel element vibrations. President Kennedy visited the NTS one week after the KIWI-B4A failure and got to see the engine starting to be disassembled in the maintenance facility. The KIWI-B4D and KIWI-B4E were modified to not have the vibration problems and were tested in test cell C. The NERVA NRX program started testing in early 1964 with NRX-A1 cold flow test series (unfueled graphite core), NRX-A2 and NRX-A3 power test series up to 1122 MW for 13 minutes. In March 1966, the NRX-EST (Engine System Test) was the first breadboard using flight functional relationship and total operating time of 116 minutes. The NRX-EST demonstrated the feasibility of a hot bleed cycle. The NRX-A5 had multiple start-ups in May-June 1966 with 30.75 minutes accumulative operating time at or above 1GW. The NRX-A6 was tested in December 1969 and ran for 62 minutes at 1100 MW. Each engine had post-test examination and found various structure anomalies which were identified for correction and the fuel element corrosion rate was reduced. The Phoebus series of research reactors began testing at test cell C, in June 1965 with Phoebus 1A. Phoebus 1A operated for 10.5 minutes at 1100 MW before unexpected loss of propellant and leading to an engine breakdown. Phoebus 1B ran for 30 minutes in February of 1967. Phoebus 2A was the highest steady state reactor built at 5GW. Phoebus 2A ran for 12 minutes at 4100 MW demonstrating sufficient power is available. The Peewee test bed reactor was tested November- December 1968 in test cell C for 40 minutes at 500MW with overall performance close to pre-run predictions. The XE' engine was the only engine tested with close to a flight configuration and fired downward into a diffuser at the Engine Test Stand (ETS) in 1969. The XE' was 1100 MW and had ~28 start-ups. The nuclear furnace NF-1 was operated at 44 MW with multiple test runs at 90 minutes in the summer of 1972. The NF-1 was the last NTP reactor tested. The Rover/NERVA program was cancelled in 1973. However, before cancellation, a lot of other engineering work was conducted by Aerojet on a 75, 000 lbf prototype flight engine and by Los Alamos on a ~16,000 lbf "Small Engine" nuclear rocket design. The ground test history of NTP at the NRDS also offers many lessons learned on how best to setup, operate, emergency shutdown, and post-test examine NTP engines. The reactor and engine maintenance and disassembly facilities were used for assembly and inspection of radioactive engines after testing. Most reactor/ engines were run at test cell A or test cell C with open air exhaust. The Rover/NERVA program became aware of a new environmental regulation that would restrict the amount of radioactive particulates allowed to be release in open air and successfully demonstrated a scrubber concept with the NF-1. The ETS stand was the only one with a high altitude test chamber used for XE'. The ETS and other test cells showed the effects the engine's radiation had on the facility materials and instrumentation as well as side effects the ground test facility has back on the engine operation. The breakdown of Phoebus 1A at test cell C showed how the site was cleaned up and back to operation for five more engines before the program was cancelled.

Gerrish, Harold P.↗

Bubble Generation in a Continuous Liquid Flow Under Reduced Gravity Conditions

The present work reports a study of bubble generation under reduced gravity conditions for both co-flow and cross-flow configurations. Experiments were performed aboard the DC-9 Reduced Gravity Aircraft at NASA Glenn Research Center, using an air-water system. Three different flow tube diameters were used: 1.27, 1.9, and 2.54 cm. Two different ratios of air injection nozzle to tube diameters were considered: 0.1 and 0.2. Gas and liquid volumetric flow rates were varied from 10 to 200 ml/s. It was experimentally observed that with increasing superficial liquid velocity, the bubbles generated decreased in size. The bubble diameter was shown to increase with increasing air injection nozzle diameters. As the tube diameter was increased, the size of the detached bubbles increased. Likewise, as the superficial liquid velocity was increased, the frequency of bubble formation increased and thus the time to detach forming bubbles decreased. Independent of the flow configuration (for either single nozzle or multiple nozzle gas injection), void fraction and hence flow regime transition can be controlled in a somewhat precise manner by solely varying the gas and liquid volumetric flow rates. On the other hand, it is observed that uniformity of bubble size can be controlled more accurately by using single nozzle gas injection than by using multiple port injection, since this latter system gives rise to unpredictable coalescence of adjacent bubbles. A theoretical model, based on an overall force balance, is employed to study single bubble generation in the dynamic and bubbly flow regime. Under conditions of reduced gravity, the gas momentum flux enhances bubble detachment; however, the surface tension forces at the nozzle tip inhibits bubble detachment. Liquid drag and inertia can act either as attaching or detaching force, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with performed experiments. However, at higher superficial,liquid velocities, the bubble neck length begins to significantly deviate from the value of the air injection nozzle diameter and thus the theory no longer predicts the experiment behavior. Effects of fluid properties, injection geometry and flow conditions on generated bubble size are investigated using the theoretical model. It is shown that bubble diameter is larger in a reduced gravity environment than in a normal gravity environment at similar flow condition and flow geometry.

Pais, Salvatore Cezar↗

High-response on-line gas analysis system for hydrogen-reaction combustion products

The results of testing an on-line quadrupole gas analyzer system are reported. Gas samples were drawn from the exhaust of a hydrogen-oxygen-nitrogen rocket which simulated the flow composition and dynamics at the combustor exit of a supersonic combustion ramjet engine. System response time of less than 50 milliseconds was demonstrated, with analytical accuracy estimated to be + or - 5 percent. For more complex chemical systems with interfering atom patterns, analysis would be more difficult. A cooled-gas pyrometer probe was evaluated as a total temperature indicator and as the primary mass flow measuring element for the total sample flow rate.

Metzler, A. J.↗

Spacelab 2 Upper Atmospheric Modification Experiment over Arecibo. II - Plasma dynamics

Results are presented from an experiment performed on Spacelab 2 over Arecibo to study the neutral gas dynamics of supersonic flows in a rarefield atmosphere and to modify the plasma density by releasing chemically reactive vapors. Exhaust vapor was released at an altitude of 317 km, where the plasma density was 300,000/cu cm. Observations were made with high resolution incoherent scatter radar. A localized depletion formed in the ionosphere. The depletion fell and eventually disappeared within the bottomside F-region ionosphere. The dynamics of the evolution of the depletion are discussed. Optical and radar data are compared, setting an upper limit of 3 percent for the branching ratio to produce O(D-1) from dissociative recombination of CO(2+) and electrons.

Bernhardt, P. A.↗

Unsteady aerodynamic analysis of space shuttle vehicles. Part 1: Summary report

An analysis of the unsteady aerodynamics of space shuttle vehicles was performed. The results show that slender wing theory can be modified to give the potential flow static and dynamic characteristics over a large Mach number range from M = 0 to M 1. A semi-empirical analytic approximation is derived for the loads induced by the leading edge vortex; and it is shown that the developed analytic technique gives good prediction of experimentally determined steady and unsteady delta wing aerodynamics, including the effects of leading edge roundness. At supersonic speeds, attached leading edge flow is established and shock-induced flow separation effects become of concern. Analysis of experimental results for a variety of boost configurations led to a definition of the main features of the flow interference effects between orbiter (delta wing) and booster. The effects of control deflection on the unsteady aerodynamics of the delta-wing orbiter were also evaluated.

Ericsson, L. E.↗

Analytical modeling requirements for tilting proprotor aircraft dynamics

Proprotor and cantilever wing aeroelastic behavior is applied to a gimballed rotor and a hingeless rotor to develop an analytical model for prediction of tilting proprotor aircraft dynamics. Particular attention is given to: the influence of coupled flap/lag bending modes; the influence of rotor blade torsion degrees of freedom on proprotor dynamics; and, to a constant coefficient approximation representing the dynamics in nonaxial flow through the rotor. The following are also examined: the number of blade bending and torsion modes required; the influence of the rotor aerodynamic model; the influence of the blade trim bending deflection; the importance of the rotor rotational speed degree of freedom; and the effect of the wing aerodynamic forces. The origin of the significant influence of the blade pitch motion on the proprotor dynamics is discussed.

Johnson, W.↗

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗

A fully coupled flow simulation around spacecraft in low earth orbit

The primary objective of this investigation is to provide a full flow simulation of a spacecraft in low earth orbit (LEO). Due to the nature of the environment, the simulation includes the highly coupled effects of neutral particle flow, free stream plasma flow, nonequilibrium gas dynamics effects, spacecraft charging and electromagnetic field effects. Emphasis is placed on the near wake phenomenon and will be verified in space by the Wake Shield Facility (WSF) and developed for application to Space Station conditions as well as for other spacecraft. The WSF is a metallic disk-type structure that will provide a controlled space platform for highly accurate measurements. Preliminary results are presented for a full flow around a metallic disk.

Justiz, C. R.↗

The Effect of Acoustic Disturbances on the Operation of the Space Shuttle Main Engine Fuel Flowmeter

The Space Shuttle Main Engine (SSME) uses a turbine fuel flowmeter (FFM) in its Low Pressure Fuel Duct (LPFD) to measure liquid hydrogen flowrates during engine operation. The flowmeter is required to provide accurate and robust measurements of flow rates ranging from 10000 to 18000 GPM in an environment contaminated by duct vibration and duct internal acoustic disturbances. Errors exceeding 0.5% can have a significant impact on engine operation and mission completion. The accuracy of each sensor is monitored during hot-fire engine tests on the ground. Flow meters which do not meet requirements are not flown. Among other parameters, the device is screened for a specific behavior in which a small shift in the flow rate reading is registered during a period in which the actual fuel flow as measured by a facility meter does not change. Such behavior has been observed over the years for specific builds of the FFM and must be avoided or limited in magnitude in flight. Various analyses of the recorded data have been made prior to this report in an effort to understand the cause of the phenomenon; however, no conclusive cause for the shift in the instrument behavior has been found. The present report proposes an explanation of the phenomenon based on interactions between acoustic pressure disturbances in the duct and the wakes produced by the FFM flow straightener. Physical insight into the effects of acoustic plane wave disturbances was obtained using a simple analytical model. Based on that model, a series of three-dimensional unsteady viscous flow computational fluid dynamics (CFD) simulations were performed using the MSFC PHANTOM turbomachinery code. The code was customized to allow the FFM rotor speed to change at every time step according to the instantaneous fluid forces on the rotor, that, in turn, are affected by acoustic plane pressure waves propagating through the device. The results of the simulations show the variation in the rotation rate of the flowmeter due to the interaction of the flow straightener wakes and the upstream propagating acoustic waves. A detailed analysis of the acoustic disturbance effects is presented along with an assessment of the impact on measurement accuracy.

Marcu, Bogdan↗

Investigation of Compressibility Effect for Aeropropulsive Shear Flows

Rocket Based Combined Cycle (RBCC) engines operate within a wide range of Mach numbers and altitudes. Fundamental fluid dynamic mechanisms involve complex choking, mass entrainment, stream mixing and wall interactions. The Propulsion Research Center at the University of Alabama in Huntsville is involved in an on- going experimental and numerical modeling study of non-axisymmetric ejector-based combined cycle propulsion systems. This paper attempts to address the modeling issues related to mixing, shear layer/wall interaction in a supersonic Strutjet/ejector flow field. Reynolds Averaged Navier-Stokes (RANS) solutions incorporating turbulence models are sought and compared to experimental measurements to characterize detailed flow dynamics. The effect of compressibility on fluids mixing and wall interactions were investigated using an existing CFD methodology. The compressibility correction to conventional incompressible two- equation models is found to be necessary for the supersonic mixing aspect of the ejector flows based on 2-D simulation results. 3-D strut-base flows involving flow separations were also investigated.

Balasubramanyam, M. S.↗

Automatic control of cryogenic wind tunnels

Inadequate Reynolds number similarity in testing of scaled models affects the quality of aerodynamic data from wind tunnels. This is due to scale effects of boundary-layer shock wave interaction which is likely to be severe at transonic speeds. The idea of operation of wind tunnels using test gas cooled to cryogenic temperatures has yielded a quantrum jump in the ability to realize full scale Reynolds number flow similarity in small transonic tunnels. In such tunnels, the basic flow control problem consists of obtaining and maintaining the desired test section flow parameters. Mach number, Reynolds number, and dynamic pressure are the three flow parameters that are usually required to be kept constant during the period of model aerodynamic data acquisition. The series of activity involved in modeling, control law development, mechanization of the control laws on a microcomputer, and the performance of a globally stable automatic control system for the 0.3-m Transonic Cryogenic Tunnel (TCT) are discussed. A lumped multi-variable nonlinear dynamic model of the cryogenic tunnel, generation of a set of linear control laws for small perturbation, and nonlinear control strategy for large set point changes including tunnel trajectory control are described. The details of mechanization of the control laws on a 16 bit microcomputer system, the software features, operator interface, the display and safety are discussed. The controller is shown to provide globally stable and reliable temperature control to + or - 0.2 K, pressure to + or - 0.07 psi and Mach number to + or - 0.002 of the set point value. This performance is obtained both during large set point commands as for a tunnel cooldown, and during aerodynamic data acquisition with intrusive activity like geometrical changes in the test section such as angle of attack changes, drag rake movements, wall adaptation and sidewall boundary-layer removal. Feasibility of the use of an automatic Reynolds number control mode with fixed Mach number control is demonstrated.

Balakrishna, S.↗

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗

Large eddy simulation of a boundary layer with concave streamwise curvature

One of the most exciting recent developments in the field of large eddy simulation (LES) is the dynamic subgrid-scale model. The dynamic model concept is a general procedure for evaluating model constants by sampling a band of the smallest scales actually resolved in the simulation. To date, the procedure has been used primarily in conjunction with the Smagorinsky model. The dynamic procedure has the advantage that the value of the model constant need not be specified a priori, but rather is calculated as a function of space and time as the simulation progresses. This feature makes the dynamic model especially attractive for flows in complex geometries where it is difficult or impossible to calibrate model constants. The dynamic model was highly successful in benchmark tests involving homogeneous and channel flows. Having demonstrated the potential of the dynamic model in these simple flows, the overall direction of the LES effort at CTR shifted toward an evaluation of the model in more complex situations. The current test cases are basic engineering-type flows for which Reynolds averaged approaches were unable to model the turbulence to within engineering accuracy. Flows currently under investigation include a backward-facing step, wake behind a circular cylinder, airfoil at high angles of attack, separated flow in a diffuser, and boundary layer over a concave surface. Preliminary results from the backward-facing step and cylinder wake simulations are encouraging. Progress on the LES of a boundary layer on a concave surface is discussed. Although the geometry of a concave wall is not very complex, the boundary layer that develops on its surface is difficult to model due to the presence of streamwise Taylor-Gortler vortices. These vortices arise as a result of a centrifugal instability associated with the convex curvature.

Lund, Thomas S.↗

Machine learning methods for particle stress development in suspension Poiseuille flows

Numerical simulations are used to study the dynamics of a developing suspension Poiseuille flow with monodispersed and bidispersed neutrally buoyant particles in a planar channel, and machine learning is applied to learn the evolving stresses of the developing suspension. The particle stresses and pressure develop on a slower time scale than the volume fraction, indicating that once the particles reach a steady volume fraction profile, they rearrange to minimize the contact pressure on each particle. Here we consider how the stress development leads to particle migration, time scales for stress development, and present a new physics-informed Galerkin neural network that allows for learning the particle stresses when direct measurements are not possible. The particle fluxes are compared with the Suspension Balance Model with good agreement. We show that when stress measurements are possible, the MOR-physics operator learning method can also capture the particle stresses.

97 MATHEMATICS AND COMPUTING↗