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 existence of an accretion column in 32 Cyg is reported, and the presence of the accretion column in Zeta Aurigae is confirmed. Secular variations have been detected for Zeta Aur. Refinements are suggested for the model of Ahmad, Chapman, and Kondo (1983). Resonance line equivalent widths have been calculated and the implications for flow geometry, secular variability, and plasma density are discussed.
The initial research effort was an in-depth analysis of the shuttle main engine plumes in an effort to improve the flowfield model and to enhance shuttle base heating equipment predictions during ascent. A prediction methodology code was developed incorporating the improved plume model into a predictive tool which could consider different trajectoreis and engine perfromance variables. Various plume flow model improvement studies were ongoing at the time of the 51-L accident. Since that time, base heating and plume methodology improvements have continued as part of the overall emphasis on Shuttle design assurance before resuming flight schedule.
A solution-adaptive grid algorithm has been developed for use in two and three dimensions. The algorithm uses a transformation from the cartesian coordinate system to a general coordinate space, which will be defined as a parallelepiped. A weighting function for adaption of the grid is developed that will allow adaption to the gradients of any combination of dependent variables in the flow. The adaption is carried out in the parametric space and a simple inverse mapping to return the new parametric space to the physical space is derived. The concept used to relocate the grid-points in the parametric space is based on the center of mass of distributed weights. Solution-adaptive results are presented for various laminar flows in two dimensions and for mathematical weighting functions in three dimensions.
The overall objective is to establish theoretically and confirm experimentally the ultimate capabilities of continuous flow electrophoresis chambers operating in an environment essentially free of particle sedimentation and buoyancy. The efforts are devoted to: (1) studying the effects of particle concentration on sample conductivity and dielectric constant. The dielectric constant and conductivity were identified as playing crucial roles in the behavior of the sample and on the resolving power and throughput of continuous flow devices; and (2) improving the extant mathematical models to predict flow fields and particle trajectories in continuous flow electrophoresis. A dielectric spectrometer was designed and built to measure the complex dielectric constant of a colloidal dispersion as a function of frequency between 500 Hz and 200 kHz. The real part of the signal can be related to the sample's conductivity and the imaginary part to its dielectric constant. Measurements of the dielectric constants of several different dispersions disclosed that the dielectric constants of dilute systems of the sort encountered in particle electrophoresis are much larger than would be expected based on the extant theory. Experiments were carried out to show that, in many cases, this behavior is due to the presence of a filamentary structure of small hairs on the particle surface. A technique for producing electrokinetically ideal synthetic latex particles by heat treating was developed. Given the ubiquitous nature of hairy surfaces with both cells and synthetic particles, it was deemed necessary to develop a theory to explain their behavior. A theory for electrophoretic mobility of hairy particles was developed. Finally, the extant computer programs for predicting the structure of electro-osmotically driven flows were extended to encompass flow channels with variable wall mobilities.
Numerical boundary conditions for mass injection/suction at the wall are incorporated in the thin-layer Navier-Stokes code, F3D. The accuracy of the boundary conditions and the code is assessed by a detailed comparison of the predictions of velocity distributions and skin-friction coefficients with exact similarity solutions for laminar flow over a flat plate with variable blowing/suction, and measurements for turbulent flow past a flat plate with uniform blowing. In laminar flow, F3D predictions for friction coefficient compare well with exact similarity solution with and without suction, but produces large errors at moderate-to-large values of blowing. A slight Mach number dependence of skin-friction coefficient due to blowing in turbulent flow is computed by F3D code. Predicted surface pressures for turbulent flow past an airfoil with mass injection are in qualitative agreement with measurements for a flat plate.
The effect of the vertical scale is investigated by considering a simple kinematic two-scale model of fluid flow inducing a variable magnetic field. Depending on the time constant, the induced magnetic field displays a variety of behaviors and geometries. In the high-frequency case, for example, a strong magnetic field tangential to the core mantle boundary, and hidden in the Delta layer, can be generated. A detailed computation and description of this magnetic field are presented. Some possible features of the secular variation of the actual geomagnetic field are discussed in the light of the model proposed here.
This paper is concerned with the effect of a weak spanwise-variable mean-flow distortion on the growth of oblique instability waves in a Blasius boundary layer. The streamwise component of the distortion velocity initially grows linearly with increasing streamwise distance, reaches a maximum, and eventually decays through the action of viscosity. This decay occurs slowly and allows the distortion to destabilize the Blasius flow over a relatively large streamwise region. It is shown that even relatively weak distortions can cause certain oblique Rayleigh instability waves to grow much faster than the usual two-dimensional Tollmien-Schlichting waves that would be the dominant instability modes in the absence of the distortion. The oblique instability waves can then become large enough to interact nonlinearly within a common critical layer. It is shown that the resulting nonlinearity is weak and that the common amplitude of the interacting oblique waves is governed by the amplitude evolution equation derived in Goldstein & Choi (1989). The implications of these results for Klebanoff-type transition are discussed.
A numerical analysis of forebody tangential slot blowing as a means of generating side force and yawing moment is conducted using an aircraft geometry. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved using a partially flux-split, approximately-factored algorithm. An algebraic turbulence model is used to determine the turbulent eddy viscosity values. Solutions are obtained using both patched and overset grid systems. In the patched grid model, and actuator plane is used to introduce jet variables into the flow field. The overset grid model is used to model the physical slot geometry and facilitate modeling of the full aircraft configuration. A slot optimization study indicates that a short slot located close to the nose of the aircraft provided the most side force and yawing moment per unit blowing coefficient. Comparison of computed surface pressure with that obtained in full-scale wind tunnel tests produce good agreement, indicating the numerical method and grid system used in the study are valid. Full aircraft computations resolve the changes in vortex burst point due to blowing. A time-accurate full-aircraft solution shows the effect of blowing on the changes in the frequency of the aerodynamic loads over the vertical tails. A study of the effects of freestream Mach number and various jet parameters indicates blowing remains effective through the transonic Mach range. An investigation of the force onset time lag associated with forebody blowing shows the lag to be minimal. The knowledge obtained in this study may be applied to the design of a forebody tangential slot blowing system for use on flight aircraft.
An experimental investigation has been conducted in the 2-inch helium tunnel at the Langley Research Center at a Mach number of 19.4 to determine the pressure distributions and heat-transfer characteristics of a family of reentry nose shapes. The pressure and heat-transfer-rate distributions on the nose shapes are compared with theoretical predictions to ascertain the limitations and validity of the theories at hypersonic speeds. The experimental results were found to be adequately predicted by existing theories. Two of the nose shapes were tested with variable-length flow-separation spikes. The results obtained by previous investigators of spike-nose bodies were found to prevail at the higher Mach number of the present investigation.
We have implemented a three-dimensional compressible Navier-Stokes code on the Connection Machine CM-5. The code is set up for implicit time-stepping on single or multiple structured grids. For multiple grids and geometrically complex problems, we follow the 'chimera' approach, where flow data on one zone is interpolated onto another in the region of overlap. We will describe our design philosophy and give some timing results for the current code. A parallel machine like the CM-5 is well-suited for finite-difference methods on structured grids. The regular pattern of connections of a structured mesh maps well onto the architecture of the machine. So the first design choice, finite differences on a structured mesh, is natural. We use centered differences in space, with added artificial dissipation terms. When numerically solving the Navier-Stokes equations, there are liable to be some mesh cells near a solid body that are small in at least one direction. This mesh cell geometry can impose a very severe CFL (Courant-Friedrichs-Lewy) condition on the time step for explicit time-stepping methods. Thus, though explicit time-stepping is well-suited to the architecture of the machine, we have adopted implicit time-stepping. We have further taken the approximate factorization approach. This creates the need to solve large banded linear systems and creates the first possible barrier to an efficient algorithm. To overcome this first possible barrier we have considered two options. The first is just to solve the banded linear systems with data spread over the whole machine, using whatever fast method is available. This option is adequate for solving scalar tridiagonal systems, but for scalar pentadiagonal or block tridiagonal systems it is somewhat slower than desired. The second option is to 'transpose' the flow and geometry variables as part of the time-stepping process: Start with x-lines of data in-processor. Form explicit terms in x, then transpose so y-lines of data are in-processor. Form explicit terms in y, then transpose so z-lines are in processor. Form explicit terms in z, then solve linear systems in the z-direction. Transpose to the y-direction, then solve linear systems in the y-direction. Finally transpose to the x direction and solve linear systems in the x-direction. This strategy avoids inter-processor communication when differencing and solving linear systems, but requires a large amount of communication when doing the transposes. The transpose method is more efficient than the non-transpose strategy when dealing with scalar pentadiagonal or block tridiagonal systems. For handling geometrically complex problems the chimera strategy was adopted. For multiple zone cases we compute on each zone sequentially (using the whole parallel machine), then send the chimera interpolation data to a distributed data structure (array) laid out over the whole machine. This information transfer implies an irregular communication pattern, and is the second possible barrier to an efficient algorithm. We have implemented these ideas on the CM-5 using CMF (Connection Machine Fortran), a data parallel language which combines elements of Fortran 90 and certain extensions, and which bears a strong similarity to High Performance Fortran. We make use of the Connection Machine Scientific Software Library (CMSSL) for the linear solver and array transpose operations.
A robust optimization is demonstrated on a two-dimensional inviscid airfoil problem in subsonic flow. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables), an approximate first-order statistical moment method is employed to represent the Computational Fluid Dynamics (CFD) code outputs as expected values with variances. These output quantities are used to form the objective function and constraints. The constraints are cast in probabilistic terms; that is, the probability that a constraint is satisfied is greater than or equal to some desired target probability. Gradient-based robust optimization of this stochastic problem is accomplished through use of both first and second-order sensitivity derivatives. For each robust optimization, the effect of increasing both input standard deviations and target probability of constraint satisfaction are demonstrated. This method provides a means for incorporating uncertainty when considering small deviations from input mean values.
Two methods that exploit the availability of sensitivity derivatives are successfully employed to predict uncertainty propagation through Computational Fluid Dynamics (CFD) code for an inviscid airfoil problem. An approximate statistical second-moment method and a Sensitivity Derivative Enhanced Monte Carlo (SDEMC) method are successfully demonstrated on a two-dimensional problem. First- and second-order sensitivity derivatives of code output with respect to code input are obtained through an efficient incremental iterative approach. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables); these sensitivity derivatives enable one to formulate first- and second-order Taylor Series approximations for the mean and variance of CFD output quantities. Additionally, incorporation of the first-order sensitivity derivatives into the data reduction phase of a conventional Monte Carlo (MC) simulation allows for improved accuracy in determining the first moment of the CFD output. Both methods are compared to results generated using a conventional MC method. The methods that exploit the availability of sensitivity derivatives are found to be valid when considering small deviations from input mean values.
The solar wind is a direct manifestation of the coronal heating processes which continue to elude us. For over three decades, observations in interplanetary space have identified two types of wind: a slow component with highly variable physical properties also characterized by speeds typically beow 500 kn/s, and a much less variable fast wind flowing on average at 750 km/s1.
We have studied the application of the baseline Model Predictive Control (MPC) algorithm to the control of main fuel flow rate (WF36), variable bleed valve (AE24) and variable stator vane (STP25) control of a simulated high-bypass turbofan engine. Using reference trajectories for thrust and turbine inlet temperature (T41) generated by a simulated new engine, we have examined MPC for tracking these two reference outputs while controlling a deteriorated engine. We have examined the results of MPC control for six different transients: two idle-to-takeoff transients at sea level static (SLS) conditions, one takeoff-to-idle transient at SLS, a Bode power command and reverse Bode power command at 20,000 ft/Mach 0.5, and a reverse Bode transient at 35,000 ft/Mach 0.84. For all cases, our primary focus was on the computational effort required by MPC for varying MPC update rates, control horizons, and prediction horizons. We have also considered the effects of these MPC parameters on the performance of the control, with special emphasis on the thrust tracking error, the peak T41, and the sizes of violations of the constraints on the problem, primarily the booster stall margin limit, which for most cases is the lone constraint that is violated with any frequency.
The NASA Fundamental Aeronautics Hypersonics project is focused on technologies for combined cycle, airbreathing propulsions systems to enable reusable launch systems for access to space. Turbine Based Combined Cycle (TBCC) propulsion systems offer specific impulse (Isp) improvements over rocket-based propulsion systems in the subsonic takeoff and return mission segments and offer improved safety. The potential to realize more aircraft-like operations with expanded launch site capability and reduced system maintenance are additional benefits. The most critical TBCC enabling technologies as identified in the National Aeronautics Institute (NAI) study were: 1) mode transition from the low speed propulsion system to the high speed propulsion system, 2) high Mach turbine engine development, 3) transonic aero-propulsion performance, 4) low-Mach-number dual-mode scramjet operation, 5) innovative 3-D flowpath concepts and 6) innovative turbine based combined cycle integration. To address several of these key TBCC challenges, NASA s Hypersonics project (TBCC Discipline) initiated an experimental mode transition task that includes an analytic research endeavor to assess the state-of-the-art of propulsion system performance and design codes. This initiative includes inlet fluid and turbine performance codes and engineering-level algorithms. This effort has been focused on the Combined Cycle Engine Large-Scale Inlet Mode Transition Experiment (CCE LIMX) which is a fully integrated TBCC propulsion system with flow path sizing consistent with previous NASA and DoD proposed Hypersonic experimental flight test plans. This experiment is being tested in the NASA-GRC 10 x 10 Supersonic Wind Tunnel (SWT) Facility. The goal of this activity is to address key hypersonic combined-cycle-engine issues: (1) dual integrated inlet operability and performance issues unstart constraints, distortion constraints, bleed requirements, controls, and operability margins, (2) mode-transition constraints imposed by the turbine and the ramjet/scramjet flow paths (imposed variable geometry requirements), (3) turbine engine transients (and associated time scales) during transition, (4) high-altitude turbine engine re-light, and (5) the operating constraints of a Mach 3-7 combustor (specific to the TBCC). The model will be tested in several test phases to develop a unique TBCC database to assess and validate design and analysis tools and address operability, integration, and interaction issues for this class of advanced propulsion systems. The test article and all support equipment is complete and available at the facility. The test article installation and facility build-up in preparation for the inlet performance and operability characterization is near completion and testing is planned to commence in FY11.
Pump Design (PUMPDES) is a computer program for designing a rotating pump for liquid hydrogen, liquid oxygen, liquid nitrogen, water, methane, or ethane. Using realistic properties of these fluids provided by another program called GASPAK, this code performs a station-by-station, mean-line analysis along the pump flow path, obtaining thermodynamic properties of the pumped fluid at each station and evaluating hydraulic losses along the flow path. The variables at each station are obtained under constraints that are consistent with the underlying physical principles. The code evaluates the performance of each stage and the overall pump. In addition, by judiciously choosing the givens and the unknowns, the code can perform a geometric inverse design function: that is, it can compute a pump geometry that yields a closest approximation of given design point. The code contains two major parts: one for an axial-rotor/inducer and one for a multistage centrifugal pump. The inducer and the centrifugal pump are functionally integrated. The code can be used in designing and/or evaluating the inducer/centrifugal-pump combination or the centrifugal pump alone. The code is written in standard Fortran 77.
In support of NASA's Entry, Descent, and Landing technology development efforts, testing of Langley's Trim Tab Parametric Models was conducted in Test Section 2 of NASA Langley's Unitary Plan Wind Tunnel. The objectives of these tests were to generate quantitative aerodynamic data and qualitative surface pressure data for experimental and computational validation and aerodynamic database development. Six component force-and-moment data were measured on 38 unique, blunt body trim tab configurations at Mach numbers of 2.5, 3.5, and 4.5, angles of attack from -4deg to +20deg, and angles of sideslip from 0deg to +8deg. Configuration parameters investigated in this study were forebody shape, tab area, tab cant angle, and tab aspect ratio. Pressure Sensitive Paint was used to provide qualitative surface pressure mapping for a subset of these flow and configuration variables. Over the range of parameters tested, the effects of varying tab area and tab cant angle were found to be much more significant than varying tab aspect ratio relative to key aerodynamic performance requirements. Qualitative surface pressure data supported the integrated aerodynamic data and provided information to aid in future analyses of localized phenomena for trim tab configurations.
Advances in the understanding of multiphase flow characteristics under variable gravity conditions will ultimately lead to improved and as of yet unknown process designs for advanced space missions. Such novel processes will be of paramount importance to the success of future manned space exploration as we venture into our solar system and beyond. In addition, because of the ubiquitous nature and vital importance of biological and environmental processes involving airwater mixtures, knowledge gained about fundamental interactions and the governing properties of these mixtures will clearly benefit the quality of life here on our home planet. The techniques addressed in the current research involving multiphase transport in porous media and gas-liquid phase separation using capillary pressure gradients are also a logical candidate for a future International Space Station (ISS) flight experiment. Importantly, the novel and potentially very accurate Lattice-Boltzmann (LB) modeling of multiphase transport in porous media developed in this work offers significantly improved predictions of real world fluid physics phenomena, thereby promoting advanced process designs for both space and terrestrial applications.This 3-year research effort has culminated in the design and testing of a zero-g demonstration prototype. Both the hydrophilic (glass) and hydrophobic (Teflon) media Capillary Pressure Gradient (CPG) cartridges prepared during the second years work were evaluated. Results obtained from ground testing at 1-g were compared to those obtained at reduced gravities spanning Martian (13-g), Lunar (16-g) and zero-g. These comparisons clearly demonstrate the relative strength of the CPG phenomena and the efficacy of its application to meet NASAs unique gas-liquid separation (GLS) requirements in non-terrestrial environments.LB modeling software, developed concurrently with the zero-g test effort, was shown to accurately reproduce observed CPG driven gas-liquid separation phenomena. The design and fabrication of a micropost plate-lamina Hele-Shaw (HS) cell was performed which served as a computationally attainable geometric structure facilitating direct comparison between physical phenomena observed in our laboratory and the LB software predictions.