Search NASASearch

SEARCH · Search NASA

Results for “Parallel Flow Path modeling”

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.

Stage-by-Stage and Parallel Flow Path Compressor Modeling for a Variable Cycle Engine

This paper covers the development of stage-by-stage and parallel flow path compressor modeling approaches for a Variable Cycle Engine. The stage-by-stage compressor modeling approach is an extension of a technique for lumped volume dynamics and performance characteristic modeling. It was developed to improve the accuracy of axial compressor dynamics over lumped volume dynamics modeling. The stage-by-stage compressor model presented here is formulated into a parallel flow path model that includes both axial and rotational dynamics. This is done to enable the study of compressor and propulsion system dynamic performance under flow distortion conditions. The approaches utilized here are generic and should be applicable for the modeling of any axial flow compressor design.

Parallel Flow Path modeling

Stage-by-Stage and Parallel Flow Path Compressor Modeling for a Variable Cycle Engine, NASA Advanced Air Vehicles Program - Commercial Supersonic Technology Project - AeroServoElasticity

This paper covers the development of stage-by-stage and parallel flow path compressor modeling approaches for a Variable Cycle Engine. The stage-by-stage compressor modeling approach is an extension of a technique for lumped volume dynamics and performance characteristic modeling. It was developed to improve the accuracy of axial compressor dynamics over lumped volume dynamics modeling. The stage-by-stage compressor model presented here is formulated into a parallel flow path model that includes both axial and rotational dynamics. This is done to enable the study of compressor and propulsion system dynamic performance under flow distortion conditions. The approaches utilized here are generic and should be applicable for the modeling of any axial flow compressor design accurate time domain simulations. The objective of this work is as follows. Given the parameters describing the conditions of atmospheric disturbances, and utilizing the derived formulations, directly compute the transfer function poles and zeros describing these disturbances for acoustic velocity, temperature, pressure, and density. Time domain simulations of representative atmospheric turbulence can then be developed by utilizing these computed transfer functions together with the disturbance frequencies of interest.

Compressor Modeling Modeling

Dependence of velocity fluctuations on solar wind speeds: A simple analysis with IPS method

A number of theoretical works have suggested that MHD plasma fluctuations in solar winds should play an important role particularly in the acceleration of high speed winds inside or near 0.1 AU from the sun. Since velocity fluctuations in solar winds are expected to be caused by the MHD plasma fluctuations, measurements of the velocity fluctuations give clues to reveal the acceleration process of solar winds. We made interplanetary scintillation (IPS) observations at the region out of 0.1 AU to investigate dependence of velocity fluctuations on flow speeds. For evaluating the velocity fluctuation of a flow, we selected the IPS data-set acquired at 2 separate antennas which located in the projected flow direction onto the baseline plane, and tried to compare skewness of the observed cross correlation function(CCF) with skewness of modeled CCFs in which velocity fluctuations were parametrized. The integration effect of IPS along a ray path was also taken into account in the estimation of modeled CCFs. Although this analysis method is significant to derive only parallel fluctuation components to the flow directions, preliminary analyses show following results: (1) High speed winds (Vsw greater than or equal to 500 km/s out of 0.3 AU) indicate enhancement of velocity fluctuations near 0.1 AU; and (2) Low speed winds (Vsw less than or equal to 400 Km/s out of 0.3 AU) indicate small velocity fluctuations at any distances.

Misawa, H.

Flow of GE90 Turbofan Engine Simulated

The objective of this task was to create and validate a three-dimensional model of the GE90 turbofan engine (General Electric) using the APNASA (average passage) flow code. This was a joint effort between GE Aircraft Engines and the NASA Lewis Research Center. The goal was to perform an aerodynamic analysis of the engine primary flow path, in under 24 hours of CPU time, on a parallel distributed workstation system. Enhancements were made to the APNASA Navier-Stokes code to make it faster and more robust and to allow for the analysis of more arbitrary geometry. The resulting simulation exploited the use of parallel computations by using two levels of parallelism, with extremely high efficiency.The primary flow path of the GE90 turbofan consists of a nacelle and inlet, 49 blade rows of turbomachinery, and an exhaust nozzle. Secondary flows entering and exiting the primary flow path-such as bleed, purge, and cooling flows-were modeled macroscopically as source terms to accurately simulate the engine. The information on these source terms came from detailed descriptions of the cooling flow and from thermodynamic cycle system simulations. These provided boundary condition data to the three-dimensional analysis. A simplified combustor was used to feed boundary conditions to the turbomachinery. Flow simulations of the fan, high-pressure compressor, and high- and low-pressure turbines were completed with the APNASA code.

Veres, Joseph P.

Column Number Density Expressions Through M = 0 and M = 1 Point Source Plumes Along Any Straight Path

Providers of payloads carried aboard the International Space Station must conduct analyses to demonstrate that any planned gaseous venting events generate no more than a certain level of material that may interfere with optical measurements from other experiments or payloads located nearby. This requirement is expressed in terms of a maximum column number density (CND). Depending on the level of rarefaction, such venting may be characterized by effusion for low flow rates, or by a sonic distribution at higher levels. Since the relative locations of other sensitive payloads are often unknown because they may refer to future projects, this requirement becomes a search for the maximum CND along any path.In another application, certain astronomical observations make use of CND to estimate light attenuation from a distant star through gaseous plumes, such as the Fermi Bubbles emanating from the vicinity of the black hole at the center of our Milky Way galaxy, in order to infer the amount of material being expelled via those plumes.This paper presents analytical CND expressions developed for general straight paths based upon a free molecule point source model for steady effusive flow and for a distribution fitted to model flows from a sonic orifice. Among other things, in this Mach number range it is demonstrated that the maximum CND from a distant location occurs along the path parallel to the source plane that intersects the plume axis. For effusive flows this value is exactly twice the CND found along the ray originating from that point of intersection and extending to infinity along the plumes axis. For sonic plumes this ratio is reduced to about 43.

Venting

Aircraft Engine Systems

This report outlines the detailed simulation of Aircraft Turbofan Engine. The objectives were to develop a detailed flow model of a full turbofan engine that runs on parallel workstation clusters overnight and to develop an integrated system of codes for combustor design and analysis to enable significant reduction in design time and cost. The model will initially simulate the 3-D flow in the primary flow path including the flow and chemistry in the combustor, and ultimately result in a multidisciplinary model of the engine. The overnight 3-D simulation capability of the primary flow path in a complete engine will enable significant reduction in the design and development time of gas turbine engines. In addition, the NPSS (Numerical Propulsion System Simulation) multidisciplinary integration and analysis are discussed.

Veres, Joseph

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen

Modeling and Visualizing Flow of Chemical Agents Across Complex Terrain

Release of chemical agents across complex terrain presents a real threat to homeland security. Modeling and visualization tools are being developed that capture flow fluid terrain interaction as well as point dispersal downstream flow paths. These analytic tools when coupled with UAV atmospheric observations provide predictive capabilities to allow for rapid emergency response as well as developing a comprehensive preemptive counter-threat evacuation plan. The visualization tools involve high-end computing and massive parallel processing combined with texture mapping. We demonstrate our approach across a mountainous portion of North California under two contrasting meteorological conditions. Animations depicting flow over this geographical location provide immediate assistance in decision support and crisis management.

Kao, David

Study of Critical Heat Flux and Two-Phase Pressure Drop Under Reduced Gravity

The design of the two-phase flow systems which are anticipated to be utilized in future spacecraft thermal management systems requires a knowledge of two-phase flow and heat transfer phenomena in reduced gravities. This program was funded by NASA headquarters in response to NRA-91-OSSA-17 and was managed by Lewis Research Center. The main objective of this program was to design and construct a two-phase test loop, and perform a series of normal gravity and aircraft trajectory experiments to study the effect of gravity on the Critical Heat Flux (CHF) and onset of instability. The test loop was packaged on two aircraft racks and was also instrumented to generate data for two-phase pressure drop. The normal gravity tests were performed with vertical up and downflow configurations to bound the effect of gravity on the test parameters. One set of aircraft trajectory tests was performed aboard the NASA DC-9 aircraft. These tests were mainly intended to evaluate the test loop and its operational performance under actual reduced gravity conditions, and to produce preliminary data for the test parameters. The test results were used to demonstrate the applicability of the normal gravity models for prediction of the two-phase friction pressure drop. It was shown that the two-phase friction multipliers for vertical upflow and reduced gravity conditions can be successfully predicted by the appropriate normal gravity models. Limited critical heat flux data showed that the measured CHF under reduced gravities are of the same order of magnitude as the test results with vertical upflow configuration. A simplified correlation was only successful in predicting the measured CHF for low flow rates. Instability tests with vertical upflow showed that flow becomes unstable and critical heat flux occurs at smaller powers when a parallel flow path exists. However, downflow tests and a single reduced gravity instability experiment indicated that the system actually became more stable with a parallel single-phase flow path. Several design modifications have been identified which will improve the system performance for generating reduced gravity data. The modified test loop can provide two-phase flow data for a range of operating conditions and can serve as a test bed for component evaluation.

Abdollahian, Davood

Column Number Density Expressions Through M = 0 and M = 1 Point Source Plumes Along Any Straight Path

Analytical expressions for column number density (CND) are developed for optical line of sight paths through a variety of steady free molecule point source models including directionally-constrained effusion (Mach number M = 0) and flow from a sonic orifice (M = 1). Sonic orifice solutions are approximate, developed using a fair simulacrum fitted to the free molecule solution. Expressions are also developed for a spherically-symmetric thermal expansion (M = 0). CND solutions are found for the most general paths relative to these sources and briefly explored. It is determined that the maximum CND from a distant location through directed effusion and sonic orifice cases occurs along the path parallel to the source plane that intersects the plume axis. For the effusive case this value is exactly twice the CND found along the ray originating from that point of intersection and extending to infinity along the plume's axis. For sonic plumes this ratio is reduced to about 4/3. For high Mach number cases the maximum CND will be found along the axial centerline path. Keywords: column number density, plume flows, outgassing, free molecule flow.

Venting

Column Number Density Expressions Through M = 0 and M = 1 Point Source Plumes Along Any Straight Path

Analytical expressions for column number density (CND) are developed for optical line of sight paths through a variety of steady free molecule point source models including directionally-constrained effusion (Mach number M = 0) and flow from a sonic orifice (M 1). Sonic orifice solutions are approximate, developed using a fair simulacrum fitted to the free molecule solution. Expressions are also developed for a spherically-symmetric thermal expansion (M = 0). CND solutions are found for the most general paths relative to these sources and briefly explored. It is determined that the maximum CND from a distant location through directed effusion and sonic orifice cases occurs along the path parallel to the source plane that intersects the plume axis. For the effusive case this value is exactly twice the CND found along the ray originating from that point of intersection and extending to infinity along the plumes axis. For sonic plumes this ratio is reduced to about 43. For high Mach number cases the maximum CND will be found along the axial centerline path.

Column Density

Performance Capability of Single-Cavity Vortex Gaseous Nuclear Rockets

An analysis was made to determine the maximum powerplant thrust-to-weight ratio possible with a single-cavity vortex gaseous reactor in which all the hydrogen propellant must diffuse through a fuel-rich region. An assumed radial temperature profile was used to represent conduction, convection, and radiation heat-transfer effects. The effect of hydrogen property changes due to dissociation and ionization was taken into account in a hydrodynamic computer program. It is shown that, even for extremely optimistic assumptions of reactor criticality and operating conditions, such a system is limited to reactor thrust-to-weight ratios of about 1.2 x 10(exp -3) for laminar flow. For turbulent flow, the maximum thrust-to-weight ratio is less than 10(exp -3). These low thrusts result from the fact that the hydrogen flow rate is limited by the diffusion process. The performance of a gas-core system with a specific impulse of 3000 seconds and a powerplant thrust-to-weight ratio of 10(exp -2) is shown to be equivalent to that of a 1000-second advanced solid-core system. It is therefore concluded that a single-cavity vortex gaseous reactor in which all the hydrogen must diffuse through the nuclear fuel is a low-thrust device and offers no improvement over a solid-core nuclear-rocket engine. To achieve higher thrust, additional hydrogen flow must be introduced in such a manner that it will by-pass the nuclear fuel. Obviously, such flow must be heated by thermal radiation. An illustrative model of a single-cavity vortex system employing supplementary flow of hydrogen through the core region is briefly examined. Such a system appears capable of thrust-to-weight ratios of approximately 1 to 10. For a high-impulse engine, this capability would be a considerable improvement over solid-core performance. Limits imposed by thermal radiation heat transfer to cavity walls are acknowledged but not evaluated. Alternate vortex concepts that employ many parallel vortices to achieve higher hydrogen flow rates offer the possibility of sufficiently high thrust-to-weight ratios, if they are not limited by short thermal-radiation path lengths.

Ragsdale, Robert G.

Field-aligned dynamics of chemically induced perturbations to the ionosphere

Ionospheric modification using photochemically reactive vapors is studied with a one-dimensional, multiion, fluid model of plasma flow along magnetic field lines. The magnitudes of ion and electron density changes are determined by considering both chemical processes (i.e., photoionization, ion-molecule reactions, dissociative recombination, electron attachment) and transport processes (i.e., multispecies diffusion, electric currents, ambipolar electric fields). The numerical treatment in the model is not specific to any type of release or any interaction chemistry. It has been used to simulate releases of Ba, CO2, and CF3Br in the ionosphere, but generalization to other species may be easily accomplished. The results of the calculations are found to be in good agreement with experimental observations. The feasibility of modifying the parallel current paths in the auroral F region is examined.

Bernhardt, P. A.

The role of optimization in structural model refinement

To evaluate the role that optimization can play in structural model refinement, it is necessary to examine the existing environment for the structural design/structural modification process. The traditional approach to design, analysis, and modification is illustrated. Typically, a cyclical path is followed in evaluating and refining a structural system, with parallel paths existing between the real system and the analytical model of the system. The major failing of the existing approach is the rather weak link of communication between the cycle for the real system and the cycle for the analytical model. Only at the expense of much human effort can data sharing and comparative evaluation be enhanced for the two parallel cycles. Much of the difficulty can be traced to the lack of a user-friendly, rapidly reconfigurable engineering software environment for facilitating data and information exchange. Until this type of software environment becomes readily available to the majority of the engineering community, the role of optimization will not be able to reach its full potential and engineering productivity will continue to suffer. A key issue in current engineering design, analysis, and test is the definition and development of an integrated engineering software support capability. The data and solution flow for this type of integrated engineering analysis/refinement system is shown.

Lehman, L. L.

Boundary layer receptivity and control

Receptivity processes initiate natural instabilities in a boundary layer. The instabilities grow and eventually break down to turbulence. Consequently, receptivity questions are a critical element of the analysis of the transition process. Success in modeling the physics of receptivity processes thus has a direct bearing on technological issues of drag reduction. The means by which transitional flows can be controlled is also a major concern: questions of control are tied inevitably to those of receptivity. Adjoint systems provide a highly effective mathematical method for approaching many of the questions associated with both receptivity and control. The long term objective is to develop adjoint methods to handle increasingly complex receptivity questions, and to find systematic procedures for deducing effective control strategies. The most elementary receptivity problem is that in which a parallel boundary layer is forced by time-harmonic sources of various types. The characteristics of the response to such forcing form the building blocks for more complex receptivity mechanisms. The first objective of this year's research effort was to investigate how a parallel Blasius boundary layer responds to general direct forcing. Acoustic disturbances in the freestream can be scattered by flow non-uniformities to produce Tollmien-Schlichting waves. For example, scattering by surface roughness is known to provide an efficient receptivity path. The present effort is directed towards finding a solution by a simple adjoint analysis, because adjoint methods can be extended to more complex problems. In practice, flows are non-parallel and often three-dimensional. Compressibility may also be significant in some cases. Recent developments in the use of Parabolized Stability Equations (PSE) offer a promising possibility. By formulating and solving a set of adjoint parabolized equations, a method for mapping the efficiency with which external forcing excites the three-dimensional motions of a non-parallel boundary layer was developed. The method makes use of the same computationally efficient formulation that makes the PSE currently so appealing. In the area of flow control, adjoint systems offer a powerful insight into the effect of control forces. One of the simplest control strategies for boundary layers involves the application of localized mean wall suction.

Hill, D. C.

Measurement and Computation of Supersonic Flow in a Lobed Diffuser-Mixer for Trapped Vortex Combustors

The trapped vortex combustor (TVC) pioneered by Air Force Research Laboratories (AFRL) is under consideration as an alternative to conventional gas turbine combustors. The TVC has demonstrated excellent operational characteristics such as high combustion efficiency, low NO(x) emissions, effective flame stabilization, excellent high-altitude relight capability, and operation in the lean-burn or rich burn-quick quench-lean burn (RQL) modes of combustion. It also has excellent potential for lowering the engine combustor weight. This performance at low to moderate combustor mach numbers has stimulated interest in its ability to operate at higher combustion mach number, and for aerospace, this implies potentially higher flight mach numbers. To this end, a lobed diffuser-mixer that enhances the fuel-air mixing in the TVC combustor core was designed and evaluated, with special attention paid to the potential shock system entering the combustor core. For the present investigation, the lobed diffuser-mixer combustor rig is in a full annular configuration featuring sixfold symmetry among the lobes, symmetry within each lobe, and plain parallel, symmetric incident flow. During hardware cold-flow testing, significant discrepancies were found between computed and measured values for the pitot-probe-averaged static pressure profiles at the lobe exit plane. Computational fluid dynamics (CFD) simulations were initiated to determine whether the static pressure probe was causing high local flow-field disturbances in the supersonic flow exiting the diffuser-mixer and whether shock wave impingement on the pitot probe tip, pressure ports, or surface was the cause of the discrepancies. Simulations were performed with and without the pitot probe present in the modeling. A comparison of static pressure profiles without the probe showed that static pressure was off by nearly a factor of 2 over much of the radial profile, even when taking into account potential axial displacement of the probe by up to 0.25 in. (0.64 cm). Including the pitot probe in the CFD modeling and data interpretation lead to good agreement between measurement and prediction. Graphical inspection of the results showed that the shock waves impinging on the probe surface were highly nonuniform, with static pressure varying circumferentially among the pressure ports by over 10 percent in some cases. As part of the measurement methodology, such measurements should be routinely supplemented with CFD analyses that include the pitot probe as part of the flow-path geometry.

Brankovic, Andreja

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery

Defining Top-of-Atmosphere Flux Reference Level for Earth Radiation Budget Studies

To estimate the earth's radiation budget at the top of the atmosphere (TOA) from satellite-measured radiances, it is necessary to account for the finite geometry of the earth and recognize that the earth is a solid body surrounded by a translucent atmosphere of finite thickness that attenuates solar radiation differently at different heights. As a result, in order to account for all of the reflected solar and emitted thermal radiation from the planet by direct integration of satellite-measured radiances, the measurement viewing geometry must be defined at a reference level well above the earth s surface (e.g., 100 km). This ensures that all radiation contributions, including radiation escaping the planet along slant paths above the earth s tangent point, are accounted for. By using a field-of- view (FOV) reference level that is too low (such as the surface reference level), TOA fluxes for most scene types are systematically underestimated by 1-2 W/sq m. In addition, since TOA flux represents a flow of radiant energy per unit area, and varies with distance from the earth according to the inverse-square law, a reference level is also needed to define satellite-based TOA fluxes. From theoretical radiative transfer calculations using a model that accounts for spherical geometry, the optimal reference level for defining TOA fluxes in radiation budget studies for the earth is estimated to be approximately 20 km. At this reference level, there is no need to explicitly account for horizontal transmission of solar radiation through the atmosphere in the earth radiation budget calculation. In this context, therefore, the 20-km reference level corresponds to the effective radiative top of atmosphere for the planet. Although the optimal flux reference level depends slightly on scene type due to differences in effective transmission of solar radiation with cloud height, the difference in flux caused by neglecting the scene-type dependence is less than 0.1%. If an inappropriate TOA flux reference level is used to define satellite TOA fluxes, and horizontal transmission of solar radiation through the planet is not accounted for in the radiation budget equation, systematic errors in net flux of up to 8 W/sq m can result. Since climate models generally use a plane-parallel model approximation to estimate TOA fluxes and the earth radiation budget, they implicitly assume zero horizontal transmission of solar radiation in the radiation budget equation, and do not need to specify a flux reference level. By defining satellite-based TOA flux estimates at a 20-km flux reference level, comparisons with plane-parallel climate model calculations are simplified since there is no need to explicitly correct plane-parallel climate model fluxes for horizontal transmission of solar radiation through a finite earth.

Loeb, N. G.