SEARCH · Engineering Papers

Results for “fluid instabilities”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Cosmic bullets as particle accelerators and radio sources

We have simulated in two dimensions the dynamical evolution of dense gas clouds(`cosmic bullets') moving supersonically through a uniform low-density medium. The diffusive shock acceleration of relativistic protons (cosmic rays) and their dynamical feedback on the background flow are included by the two-fluid model for this process. The acceleration of relativistic electrons is approximated by a test-particle model, and a passive magnetic field is followed by a simple advection scheme. Strong bow shocks, with Mach numbers similar to that of a bullet's motion, are the most important particle accelerators in the flow, while tail shocks and shocks inside the bullets do not play generally significant roles in this regard. For our simulation parameters, approximately greater than 10% of the initial bullet kinetic energy is converted to a combination of internal energy of gas and cosmic-ray protons by the time the bullets begin to be disrupted. Characteristically, the cosmic rays gain several percent of the available kinetic energy. Bullet destruction on timescales only a little larger than the ram pressure bullet crushing time begins in response to Kelvin-Helmholtz and especially to Rayleigh-Taylor instabilities along the forward bullet surface. For dense bullets this happens before the bullet is stopped by ram pressure. According to our simple model for synchrotron emission from relativistic electrons accelerated and transported within the flows, that emission increases rapidly as the bullet begins to fragment, when it is strongly dominated by field enhancement in sheared flows. Synchrotron emission from the acceleration region within the bow shock is, by contrast, much weaker.

Jones, T. W.↗

NASA NTRS1989

Axisymmetric shapes and stability of charged drops in an external electric field

A highly conducting charged drop that is surrounded by a fluid insulator of another density can be levitated by suitably applying a uniform electric field. Axisymmetric equilibrium shapes and stability of the levitated drop are found by solving simultaneously the augmented Young-Laplace equation for surface shape and the Laplace equation for the elecric field, together with constraints of fixed drop volume, charge, and center of mass. The means are a method of subdomains, finite element basis functions, and Galerkin's method of weighted residuals, all facilitated by a large-scale computer. Shape families of fixed charge are treated systematically by first-order continuation. Previous analyses by Abbas et al. in 1967 and Abbas and Latham in 1969, in which the shapes of levitated drops are approximated as spheroids, are corrected. The new analysis shows that drops charged to less than the Rayleigh limit lose shape stability at turning points, with respect to external field strength, and that the instability seen in experiments of Doyle et al. in 1964 and others is not a bifurcation to a family of two-lobed shapes, but rather is a related imperfect bifurcation.

Basaran, O. A.↗

NASA NTRS1990

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

NASA NTRS1991

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I - The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics

Yee, H. C.↗

NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2-D, two-temperature, single fluid magnetohydrodynamic code which incorporates classical plasma transport coefficients and Hall effects was developed to predict steady-state, self-field magnetoplasmadynamic (MPD) thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects, the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries, assuming a fully ionized argon propellant. A maximum predicted specific impulse of 1680 s is obtained for a thruster with an anode radius of 2.5 cm, a cathode radius of 0.5 cm, and equal electrode lengths of 2.5 cm. A scaling relation is developed to predict, within limits, the onset of cylindrical, self-field thruster instability as a function of geometry and operating condition.

Lapointe, M.↗

NASA NTRS2011

Assessing Uncertainties in Boundary Layer Transition Predictions for HIFiRE-1 at Non-zero Angles of Attack

Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests conducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially continuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at nonzero angles of attack, so this analysis could be expanded in the future to include the results of the flight test data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and the mean flow solutions were computed using the DPLR finite volume Navier-Stokes computational fluid dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The N-factors at transition onset and the value of Re(sub theta)/M(sub e), commonly used to predict boundary layer transition onset, were compared for all conditions analyzed. Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1deg, 3deg and 5deg. At these conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary layer transition onset showed improvements over the Re(sub theta)/M(sub e) empirical correlation methodology. Stronger correlations and less scatter in the data were observed when using the N factor method for these cases.

Marek, Lindsay C.↗

NASA NTRS1995

Diffusion at the magnetopause: Hybrid simulations

Electromagnetic wave generation and resulting cross-field diffusion of plasma are considered at a tangential discontinuity, which characterizes the magnetopause for northward interplanetary magnetic field. Two-dimensional hybrid (particle ions, massless fluid electrons) simulations, in which the tangential discontinuity is generated self-consistently via a stream-stream interaction, are used to show that wave growth occurs when the ambient magnetic field is predominantly perpendicular to the direction of the density gradient. Low-frequency (much less than ion gyrofrequency) waves, with amplitudes delta B/B less than or equal to 0.2 and anticorrelated density fluctuations delta n/n less than or equal to 0.6, are generated at the discontinuity, resulting in cross-field diffusion that is comparable to the Bohm rate. Both the fluctuation level and the lack of ion heating in the calculations are consistent with observations at the magnetopause. The magnitude of the diffusion is considered in the presence of numerical effects and in the context of the inferred diffusion rate at the magnetopause. The relation of the low-frequency waves and their consequences to faster growing, short-wavelength waves due to the lower hybrid drift instability is also addressed. The overall conclusion of this initial study is that diffusion due to low frequency waves is not likely to be a major effect at the magnetopause.

Winske, D.↗

NASA NTRS2002

Lattice Boltzmann Method for 3-D Flows with Curved Boundary

In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.

Mei, Renwei↗

NASA NTRS2021

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗

NASA NTRS2020

Generation and Calibration of Linear Models of Aircraft with Highly Coupled Aeroelastic and Flight Dynamics

The lightweight structures and unconventional configurations being considered for the next generation of aircraft mean that any effort to predict or control the flight dynamics is impacted by the structural dynamics. One of the most severe forms of coupling between aeroelasticity and flight dynamics is an instability called body freedom flutter. The existing tools often assume a relatively weak effect of structural dynamics on the flight dynamics, and are therefore incapable of modeling strong interactions like body freedom flutter. A method of combining different sources of data traditionally used for aeroelasticity and flight dynamics is described by reconciling many of the differences between these models. By building upon past modeling efforts, a level of familiarity in the approach is achieved. Generally the differences from the traditional approaches are subtle but significant. The traditional frequency domain flutter model in a modal coordinate system is converted to a form consistent with a time domain flight dynamics model. The time domain rational function approximation about a non-inertial coordinate system and the unique constraints for the conversion between the inertial and non-inertial coordinate systems are discussed. A consistent transformation of the states of aeroelastic models to flight dynamics models is derived, which enables the integration of data from higher fidelity computational fluid dynamics models or wind-tunnel testing. The present method of integrating multidisciplinary data was used to create models that compare well with X-56A flight-test data, including conditions past the flutter speed.

Ouellette, Jeffrey↗

NASA NTRS1991

Turbulization of shear flows in astrophysics

The problem of instability in shear flows is treated in a number of works which present a global analysis of irrotational (potential) modes in a 2-D compressible shear flow, implying that, in the presence of reflecting boundaries, those modes can increase and that the characteristic time of the increase is several orders as high as the dynamic time in the medium. Due to the slowness of the increase of irrotational modes and to the problematic character of the theory's applicability for thin accretion disks, an attempt was made to develop an alternative theory of turbulization shear flows. It is quite contrary to the above theory. Instead of the global analysis, a local one is made far from the flow boundaries. Instead of irrotational perturbations, vortical ones are examined. Lastly, the medium is considered to be incompressible. The flow in the accretion disk can be modeled by a plane Couette flow, without regard for boundary conditions: it can be modeled by a plane shear flow in infinite incompressible fluid. According to this scenario, there is a critical perturbation level in the free shear flow. If the level of initial perturbations exceeds the critical one, the flow is turbulized.

Chagelishvili, G. D.↗

NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗

NASA NTRSDate not supplied

X-57 “Maxwell” High-Lift Propeller Testing and Model Development

NASA’s X-57 “Maxwell” distributed electric propulsion flight demonstrator has a high-lift system that includes 12 fixed-pitch high-lift propellers located upstream of the wing leading edge for lift augmentation at low speeds. These high-lift propellers are not required at higher speeds and are folded conformally along the nacelles to reduce drag when not in operation. Aircraft performance models and flight simulations incorporate propeller performance to predict thrust, moments, and power consumption and propeller model accuracy is important in identifying safe operating regimes for this aircraft. Current high-lift propeller performance models have been verified and calibrated against numerous computational fluid dynamics analyses under a variety of flight conditions. We performed a series of full-scale wind tunnel tests at the NASA Langley Research Center Low Speed Aeroacoustic Wind Tunnel to further validate these models, to identify any adverse operating conditions for these propellers, and to assess the accuracy of the facility’s new Propeller Test Stand. The results indicate that the models accurately predicted performance and that the right-handed propeller showed lower torque and thrust for a given propeller speed compared to the left-handed version. Acoustic data support these measurements, showing higher tonal noise for the left versus right propeller. We suspect this is due to slight differences in the material properties between the early production (right) and more recent (left) blade sets leading to greater detwisting in the right-handed blades under load. Both propellers demonstrated very stable operation throughout the test including during deployment, stowing, and windmilling. Furthermore, high-speed imagery of the spinning propeller indicated no blade oscillations or instabilities for any test case. We also leveraged these images to analyze the effects of propeller loading on blade detwisting and tip position. Unfortunately, the thrust measurements varied significantly during the test due to cross-loading of the Propeller Test Stand load cell, though trends in the thrust data were still observable. The torque measurement, however, remained very precise for all test points and was used for performance comparisons. As a result of this test, we validated and improved the X-57 high-lift propeller performance models, obtained detailed propeller acoustic data and high-speed imagery, observed the propeller operating safely during a variety of dynamic events, and identified two potentially adverse, low-propeller-speed conditions which will be avoided during normal aircraft operation.

X-57↗

NASA NTRS1993

LIF measurements of scalar mixing in turbulent shear layers

The structure of shear layer flows at high Reynolds numbers remains a very interesting problem. Straight mixing layers have been studied and yielded information on the probability density function (pdf) of a passive scalar across the layer. Konrad and Koochesfahani & Dimotakis measured the pdf of the mixture fraction for mixing layers of moderate Reynolds numbers, each about 25,000 (Re based on velocity difference and visual thickness). Their measurements showed a 'non-marching' pdf (central hump which is invariant from edge to edge across the layer), a result which is linked to the visualizations of the spanwise Kelvin-Helmholtz (K-H) instability mode, which is the primary instability for plane shear layer flows. A secondary instability mode, the Taylor-Gortler (T-G) instability, which is associated with streamwise vortical structures, has also been observed in shear layers. Image reconstruction by Jimenez et al. and volume renderings by Karasso & Mungal at low Re numbers have demonstrated that the K-H and the T-G instability modes occur simultaneously in a non-mutually destructive way, evidence that supports the quasi two-dimensional aspect of these flows and the non-marching character of the pdf at low Reynolds numbers. At higher Re numbers though, the interaction of these two instability modes is still unclear and may affect the mixing process. In this study, we perform measurements of the concentration pdf of plane mixing layers for different operating conditions. At a speed ratio of r = U(sub 1)/U(sub 2) = 4:1, we examine three Reynolds number cases: Re = 14,000, Re = 31,000, and Re = 62,000. Some other Re number cases' results, not presented in detail, are invoked to explain the behavior of the pdf of the concentration field. A case of r = 2.6:1 at Re = 20,000 is also considered. The planar laser-induced fluorescence technique is used to yield quantitative measurements. The different Re are obtained by changing the velocity magnitudes of the two streams. The question of resolution of these measurements is addressed. In order to investigate the effects of the initial conditions on the development and the structure of the mixing layer, the boundary layer on the high-speed side of the splitter plate is tripped. The average concentration and the average mixed fluid concentration are also calculated to further understand the changes in the shear layer for the different cases examined.

Karasso, Paris S.↗

NASA NTRS1995

SRM Internal Flow Tests and Computational Fluid Dynamic Analysis: RSRM Full Scale Motor Analyses - Volume 2

The purpose of the RSRM Nozzle Slag Ejection Precursor Test is to investigate the effect that slag ejection from the RSRM nozzle has on the chamber pressure and trust of the SRB's. In past firings of the Reusable Solid Rocket Motor (RSRM) both static test and flight motors have shown small pressure perturbations occurring primarily between 65 and 80 seconds. A joint NASA/Thiokol team investigation concluded that the cause of the pressure perturbations was the periodic ingestion and ejection of molten aluminum oxide slag from the cavity around the submerged nozzle nose which tends to trap and collect individual aluminum oxide droplets from the approach flow. The conclusions of the team were supported by numerous data and observations from special tests including high speed photographic films, real time radiography, plume calorimeters, accelerometers, strain gauges, nozzle TVC system force gauges, and motor pressure and thrust data. A simplistic slag ballistics model was formulated to relate a given pressure perturbation to a required slag quantity. Also, a cold flow model using air and water was developed to provide data on the relationship between the slag flow rate and the chamber pressure increase. Both the motor and the cold flow model exhibited low frequency oscillations in conjunction with periods of slag ejection. Motor and model frequencies were related to scaling parameters. The data indicate that there is a periodicity to the slag entrainment and ejection phenomena which is possibly related to organized oscillations from instabilities in the dividing streamline shear layer which impinges on the underneath surface of the nozzle.

Source record↗

NASA NTRS1989

Time-dependent computational studies of flames in microgravity

The research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of the NASA Microgravity Science and Applications Program is described. The primary focus was on investigating fundamental questions concerning the propagation and extinction of premixed flames in Earth gravity and in microgravity environments. The approach was to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. A combination of one-dimensional and two-dimensional simulations was used to investigate the effects of curvature and dilution on ignition and propagation of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the mechanism leading to cellular instability, and to study the effects of gravity on the transition to cellular structure. A flame in a microgravity environment can be extinguished without external losses, and the mechanism leading to cellular structure is not preferential diffusion but a thermo-diffusive instability. The simulations have also lead to a better understanding of the interactions between buoyancy forces and the processes leading to thermo-diffusive instability.

Oran, Elaine S.↗

NASA NTRS2005

The Dynamics of Miscible Interfaces: Simulations

This research project focuses on the dynamics of interfacial regions between miscible fluids. While much attention has focused on immiscible interfaces in the past, miscible interfaces have been explored to a much lesser degree, so that there are many open questions regarding their dynamics at this time. Among the more pressing issues is the role that nonconventional stresses can play in such interfacial regions. Such stresses are typically not accounted for in efforts to model the dynamics of miscible flows. Our research aims to clarify under which circumstances these stresses do have to be taken into account, and what quantitative approaches are most suitable in this regard. In order to address these issues, we have focused on conducting linear stability analyses and nonlinear simulations for capillary tube and Hele-Shaw flows, and to compare the results with corresponding experiments performed in the labs of our co-investigators Prof. Maxworthy at USC, and Dr. Balasubramaniam at NASA. Over the duration of the project we have, among other things, focused on the effects of variable diffusion coefficients in such flows, and specifically on their influence in the growth of instabilities. Furthermore, our three-dimensional spectral element simulations have made good progress, so that we have come to a point where we can conduct more detailed comparisons with experimental observations. We are currently focusing our efforts on reproducing the tip-splitting instability observed by Maxworthy. Finally, we have discovered a new core-annular flow instability in the Stokes flow regime during the last year. This represents a significant finding, as this instability does not have an immiscible counterpart.

Meiburg, Eckart↗

NASA NTRS2016

Predictive Evaluations of Oxygen-Rich Hydrocarbon Combustion Gas-Centered Swirl Coaxial Injectors using a Flamelet-Based 3-D CFD Simulation Approach

NASA Marshall Space Flight Center (MSFC) has embarked upon a joint project with the Air Force to improve the state-of-the-art of space application combustion device design and operational understanding. One goal of the project is to design, build and hot-fire test a 40,000 pound-thrust Oxygen/Rocket Propellant-2 (RP-2) Oxygen-Rich staged engine at MSFC. The overall project goals afford the opportunity to test multiple different injector designs and experimentally evaluate the any effect on the engine performance and combustion dynamics. To maximize the available test resources and benefits, pre-test, combusting flow, Computational Fluid Dynamics (CFD) analysis was performed on the individual injectors to guide the design. The results of the CFD analysis were used to design the injectors for specific, targeted fluid dynamic features and the analysis results also provided some predictive input for acoustic and thermal analysis of the main Thrust Chamber Assembly (TCA). MSFC has developed and demonstrated the ability to utilize a computationally efficient, flamelet-based combustion model to guide the pre-test design of single-element Gas Centered Swirl Coaxial (GCSC) injectors. Previous, Oxygen/RP-2 simulation models utilizing the Loci-STREAM flow solver, were validated using single injector test data from the EC-1 Air Force test facility. The simulation effort herein is an extension of the validated, CFD driven, single-injector design approach applied to single injectors which will be part of a larger engine array. Time-accurate, Three-Dimensional, CFD simulations were performed for five different classes of injector geometries. Simulations were performed to guide the design of the injector to achieve a variety of intended performance goals. For example, two GCSC injectors were designed to achieve stable hydrodynamic behavior of the propellant circuits while providing the largest thermal margin possible within the design envelope. While another injector was designed to purposefully create a hydrodynamic instability in the fuel supply circuit as predicted by the CFD analysis. Future multi-injector analysis and testing will indicate what if any changes occur in the predicted behavior for the single-element injector when the same injector geometry is placed in a multi-element array.

Richardson, Brian R.↗