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.
Results of an experimental and theoretical study of factors affecting self-sustaining hydrodynamic oscillations in boiling-water loops are reported. Data on flow variables, and the effects of geometry, subcooling and pressure on the development of oscillatory behavior in a natural-circulation two-phase-flow boiling loop are included.
In the present work a solver for the Reynolds averaged compressible Navier-Stokes equations, to compute high speed turbulent flows characterized by interacting shock waves and viscous layers, is presented. A k-epsilon turbulence model that accounts for compressibility effects is developed. The numerical algorithm is based on a finite volume multistage Runge Kutta technique that is explicit for the solution of the mean flow variables, and implicit for the solution of the k-epsilon equations. The model is validated by extensive comparison with experimental results of flows over compression ramps characterized by interacting shock waves/boundary layers.
One-dimensional isentropic flow variables of partially ionized argon - tables and data analysis
A FORTRAN 4 computer program has been developed that obtains a subsonic or shock-free transonic flow solution on the hub-shroud mid-channel flow surface of a turbomachine. The blade row may be fixed or rotating, and may be twisted and leaned. Flow may be axial or mixed, up to 45 deg from axial. Upstream and downstream flow variables may vary from hub to shroud, and provision is made to correct for loss of stagnation pressure. The results include velocities, streamlines, and flow angles on the flow surface; and approximate blade surface velocities. Subsonic solutions are obtained by a finite-difference stream-function solution. Transonic solutions are obtained by a velocity-gradient method, using information from a finite-difference stream-function solution at a reduced mass flow.
A FORTRAN-IV computer program, MERIDL, has been developed that obtains a subsonic or shock-free transonic flow solution on the hub-shroud mid-channel flow surface of a turbomachine. The blade row may be fixed or rotating and may be twisted and leaned. Flow may be axial or mixed, up to 45 deg from axial. Upstream and downstream flow variables can vary from hub to shroud, and provision is made to correct for loss of stagnation pressure. The results include velocities, streamlines, and flow angles on the flow surface and approximate blade surface velocities. Subsonic solutions are obtained by a finite-difference stream-function solution. Transonic solutions are obtained by a velocity-gradient method, using information from a finite-difference stream-function solution at a reduced mass flow.
A FORTRAN 4 computer program was developed that obtains a detailed subsonic or shock-free transonic flow solution on the hub-shroud midchannel stream surface of a turbomachine. The blade row may be fixed or rotating, and the blades may be twisted and leaned. Flow may be axial, mixed, or radial. Upstream and downstream flow variables may vary from hub to shroud, and provision is made to correct for loss of stagnation pressure. The results include velocities, streamlines, and flow angles on the stream surface as well as approximate blade surface velocities. Subsonic solutions are obtained by a finite-difference, stream-function solution. Transonic solutions are obtained by a velocity-gradient method that uses information from a finite-difference, stream-function solution at a reduced mass flow.
A FORTRAN IV computer program has been developed that obtains a detailed subsonic or shock free transonic flow solution on the hub-shroud midchannel stream surface of a turbomachine. The blade row may be fixed or rotating, and the blades may be twisted and leaned. Flow may be axial, mixed, or radial. Upstream and downstream flow variables may vary from hub to shroud, and provisions are made to correct for loss of stagnation pressure. The results include velocities, streamlines, and flow angles on the stream surface and approximate blade surface velocities.
The NASA Redox energy storage system is under active technology development. The hardware undergoing laboratory testing is either 310 sq. cm. or 929 sq. cm. (0.33 sq. ft. or 1.0 sq. ft. per cell active area with up to 40 individual cells connected to make up a modular cell stack. This size of hardware allows rather accurate projections to be made of the shunt power/pump power tradeoffs. The modeling studies that were completed on the system concept are reviewed along with the approach of mapping the performance of Redox cells over a wide range of flow rates and depths of discharge of the Redox solutions. Methods are outlined for estimating the pumping and shunt current losses for any type of cell and stack combination. These methods are applicable to a variety of pumping options that are present with Redox systems. The results show that a fully developed Redox system has acceptable parasitic losses when using a fixed flow rate adequate to meet the worst conditions of current density and depth of discharge. These losses are reduced by about 65 percent if variable flow schedules are used. The exact value of the overall parasitics will depend on the specific system requirements of current density, voltage limits, charge, discharge time, etc.
The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.
Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.
The effect of shock impingement on an axial corner simulating the inlet of a hypersonic vehicle engine is modeled using a finite-difference procedure. A three-dimensional dynamic grid adaptation procedure is utilized to move the grids to regions with strong flow gradients. The adaptation procedure uses a grid relocation stencil that is valid at both the interior and boundary points of the finite-difference grid. A linear combination of spatial derivatives of specific flow variables, calculated with finite-element interpolation functions, are used as adaptation measures. This computational procedure is used to study laminar and turbulent Mach 6 flows in the axial corner. The description of flow physics and qualitative measures of heat transfer distributions on cowl and strut surfaces obtained from the analysis are compared with experimental observations. Conclusions are drawn regarding the capability of the numerical scheme for enhanced modeling of high-speed compressible flows.
The design and techno-economic performance of a compact additively manufactured (AM) molten salt (MS)-to-supercritical carbon di-oxide (sCO2) primary heat exchanger (PHE) for solar thermal application is described. The PHE design consists of sCO2 flow through an array of microscale pin fins while the MS flows through mm-scale rectangular channels. Constraints imposed by AM using laser powder bed fusion method are considered in the design. Structural and fluid flow simulations are performed to arrive at a viable design of the core and headers. A simplified one-dimensional steady state model for the PHE is developed including the impact of surface roughness from the AM process. A process-based cost model is used to determine the tradeoff between thermofluidic design and manufacturing cost. A parametric study is performed using the thermo-fluidic and cost models to determine the set of geometrical and flow variables that result in high power density and low cost, while restricting the pressure drop on the sCO2 side to less than 2% of line pressure. Flow rates of MS and sCO2 were varied over heat capacity rate ratios ranging from 0.2 to 1. Results indicate that it is possible to design a low-pressure drop AM PHE with an effectiveness of 90% and a power density in excess of 10 MW/m3 (including headers). Fabrication of representative nickel superalloy specimens are shown to demonstrate that low-porosity parts with the requisite dimensional tolerance of PHE core can be generated.
A technology base for the thermal acoustic shield concept as a noise suppression device for single stream exhaust nozzles was developed. Acoustic data for 314 test points for 9 scale model nozzle configurations were obtained. Five of these configurations employed an unsuppressed annular plug core jet and the remaining four nozzles employed a 32 chute suppressor core nozzle. Influence of simulated flight and selected geometric and aerodynamic flow variables on the acoustic behavior of the thermal acoustic shield was determined. Laser velocimeter and aerodynamic measurements were employed to yield valuable diagnostic information regarding the flow field characteristics of these nozzles. An existing theoretical aeroacoustic prediction method was modified to predict the acoustic characteristics of partial thermal acoustic shields.
Equations for the mean and turbulent quantities for compressible turbulent flows are derived. Both the conventional Reynolds average and the mass-weighted, Favre average were employed to decompose the flow variable into a mean and a turbulent quality. These equations are to be used later in developing second order Reynolds stress models for high speed compressible flows. A few recent advances in modeling some of the terms in the equations due to compressibility effects are also summarized.
Equations for the mean and the turbulence quantities of compressible turbulent flows are derived in this report. Both the conventional Reynolds average and the mass-weighted Favre average were employed to decompose the flow variable into mean and turbulent quantities. These equations are to be used later in developing second-order Reynolds stress models for high-speed compressible flows. A few recent advances in modeling some of the terms in the equation due to compressibility effects are also summarized.
We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.
A mathematical model of a multi-stage compressor with variable molecular weight flow medium is derived. The modeled system consists of a five stage, six cylinder, double acting, piston type compressor. Each stage is followed by a water cooled heat exchanger which serves to transfer the heat of compression from the gas. A high molecular weight gas (CFC-12) mixed with air in varying proportions is introduced to the suction of the compressor. Condensation of the heavy gas may occur in the upper stage heat exchangers. The state equations for the system are integrated using the Advanced Continuous Simulation Language (ACSL) for determining the system's dynamic and steady state characteristics under varying operating conditions.
A detailed study of two strong adverse pressure gradient flows, one with a free-stream velocity of 35 m/sec, at throat (producing a Re sub theta of 27000 at detachment) and another with free-stream velocity of 22 m/sec, at throat (producing a Re sub theta of 19000 at detachment) is presented. In these examples flows separate slowly and reattach very rapidly over a very short distance in a streamwise direction. In the backflow region, there appears to be a semi-logarithmically flat region in the streamwise fluctuating velocity component, u', which spreads over a definite range of y/delta. In power spectra, the flow variables phi sub upsilon upsilon (kappa sub 1 delta)/ -uv bar sub max vs. kappa sub 1 delta forms a unique set of scaling parameters for adverse pressure gradient flows. Experimental results show good agreement with previous studies.