Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass conservation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Computer Modeling Of Atomization

Improved mathematical models based on fundamental principles of conservation of mass, energy, and momentum developed for use in computer simulation of atomization of jets of liquid fuel in rocket engines. Models also used to study atomization in terrestrial applications; prove especially useful in designing improved industrial sprays - humidifier water sprays, chemical process sprays, and sprays of molten metal. Because present improved mathematical models based on first principles, they are minimally dependent on empirical correlations and better able to represent hot-flow conditions that prevail in rocket engines and are too severe to be accessible for detailed experimentation.

Giridharan, M.↗
NASA NTRS1992

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗
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 NTRS1995

Analysis of gas absorption to a thin liquid film in the presence of a zero-order chemical reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Rajagopalan, S.↗
NASA NTRS1995

Transient studies of G-induced capillary dryout and rewet

A transient, one-dimensional numerical code is developed to model the liquid motion in an axial groove with square cross section. Axial variation in liquid level, shear stress and heat transfer between the groove wall and the liquid, evaporation and transient body forces are accounted for in the model. Dryout and rewet of the groove are allowed; the front location is determined numerically using conservation of mass and linear extrapolation. Several numerical test results are presented and discussed.

Reagan, M. K.↗
NASA NTRS1994

Quasi-periodic atmosphere-regolith-cap CO2 redistribution in the martian past

Our earlier Mars regolith-atmosphere-cap CO2 distribution model (Fanale et al., 1982, Icarus 50, 381-407) has been improved, revised, and extended back over Mars' mid to late history. The present model takes into account four new factors: (1) a more realistic long-term obliquity cycle, (2) thermal conduction as it affects the surface energy balance, (3) the changing solar constant, and (4) atmospheric erosion 3.5 byr ago to the present. Solar insolation and temperatures are computed for the full range of obliquities, latitudes, and epochs, and a CO2 adsorption relation is used, together with a conservation of mass constraint, to calculate atmospheric pressures and exchangeable CO2 mass as functions of obliquity and epoch for the regolith, atmosphere, and polar caps for two assumed thicknesses of a basalt regolith. It is found that the heat conduction term in the surface boundary condition has an important effect in reducing the range of atmospheric pressures over the obliquity cycle at all epochs.

Fanale, Fraser P.↗
NASA NTRS1995

Resolution of the buoyancy in the 8-foot high temperature tunnel combustor

Currently, the 8-Foot High Temperature Tunnel (8-Ft. HTT) combustor produces a good profile at only one point (2000 psia and 3650 R with oxygen enrichment). Air is enriched with oxygen (liquid) so that the combustor product gas will contain the volumetric amount of oxygen normally found in air. The oxygen enriched air has a large fraction that is not reacted and flows through the outer periphery of the fuel injector. This ring of cold air in addition to the relatively cold walls of the combustor set up buoyancy forces that produce a segregation of relatively cool gases at the bottom of the combustor exit. The basic problem is to produce a test gas that has uniform properties at all combustor conditions. The combustor temperature may be as high as 3700 R or as low as 2000 R. Combustor pressures can be as high as 3500 psia (no oxygen enrichment) and as low as 600 psia. The segregation is most severe with oxygen enriched air, since its temperature is lower and its density is high. The combustor is lined with nickel 201 and can be operated at about 1600 R maximum. A global mixing process is desired that produces an acceptable profile of temperature, species, and velocity at the exit of the combustor. The ultimate goal is a temperature profile with about 100 R variance and about 2 percent variance in oxygen. The exit total temperature must not be lowered significantly by the mixing apparatus or mechanisms employed. If immersed bodies are used, they must also be kept very hot. All combustor wall modifications must be able to survive the heat and structural conditions of the varied operating conditions. Our approach to resolving this issue is being conducted in three stages: (1) Consider mixing exclusively, (2) Resolve the heat transfer concerns resulting from the chosen mixing strategy, and (3) Solve the material and structural problems resulting from stages (1) and (2). Since the 8-Ft. HTT is unavailable for experimentation, the study is conducted exclusively with computational fluid dynamic (CPD) codes (Fluent/Uns and Rampant 3.1) using unstructured grid through body fitted coordinates. Both CFD codes are general purpose Navier-Stokes solution packages that can solve integral conservation equations for conservation of mass, momentum and energy. The governing equations are discretized using a control-volume finite-element method on unstructured triangular 2-D grids. In the interest of time, a 3-D tetrahedral grid was used to check the 2-D results on one mixing strategy and the 2-D results were confirmed. Preliminary results indicate that excellent mixing can be achieved with a body placed in the center of the flowing hot fluid with a minor modification to the combustor wall similar to a model positioned in the test section of a wind tunnel. The concept here, is to create longitudinal vortices strong enough to bring both fluids into intimate contact with each other near the wall where the cold fluid resides. However, there may be a trade-off in the choice of mixing strategy, heat transfer and structural requirements. Currently, we are examining the more promising geometries for heat transfer concerns and developing strategies for the material of construction of the center body.

Loney, Norman W.↗
NASA NTRS1995

An estimate of solar wind velocity profiles in an coronal hole and a coronal streamer area (6-40 solar radius)

Using the total electron content data obtained by the Ulysses Solar Corona Experiment during the superior solar conjunction in summer 1991, we selected two data sets, one associated with a coronal hole and the other one with coronal streamer crossings. By doing this data splitting, we find two entirely different density profiles varying as r(exp -2.7) and r(exp -2.3) for the coronal hole and coronal streamers, respectively. Assuming mass flux conservation from the inner corona to one AU, an estimate for the velocity profiles or acceleration in these two different regions can be determined. The more negative exponent of the coronal hole density profile indicates a more extended heating and acceleration region or more flaring, or both. Various possible explanations will be discussed.

Paetzold, M.↗
NASA NTRS1996

Determination of Uncertainties for the New SSME Model

This report discusses the uncertainty analysis performed in support of a new test analysis and performance prediction model for the Space Shuttle Main Engine. The new model utilizes uncertainty estimates for experimental data and for the analytical model to obtain the most plausible operating condition for the engine system. This report discusses the development of the data sets and uncertainty estimates to be used in the development of the new model. It also presents the application of uncertainty analysis to analytical models and the uncertainty analysis for the conservation of mass and energy balance relations is presented. A new methodology for the assessment of the uncertainty associated with linear regressions is presented.

Coleman, Hugh W.↗
NASA NTRS1996

Modeling of Thermal Performance of Multiphase Nuclear Fuel Cell Under Variable Gravity Conditions

A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.

Ding, Z.↗
NASA NTRS1999

A One Dimensional, Time Dependent Inlet/Engine Numerical Simulation for Aircraft Propulsion Systems

The NASA Lewis Research Center (LeRC) and the Arnold Engineering Development Center (AEDC) have developed a closely coupled computer simulation system that provides a one dimensional, high frequency inlet/engine numerical simulation for aircraft propulsion systems. The simulation system, operating under the LeRC-developed Application Portable Parallel Library (APPL), closely coupled a supersonic inlet with a gas turbine engine. The supersonic inlet was modeled using the Large Perturbation Inlet (LAPIN) computer code, and the gas turbine engine was modeled using the Aerodynamic Turbine Engine Code (ATEC). Both LAPIN and ATEC provide a one dimensional, compressible, time dependent flow solution by solving the one dimensional Euler equations for the conservation of mass, momentum, and energy. Source terms are used to model features such as bleed flows, turbomachinery component characteristics, and inlet subsonic spillage while unstarted. High frequency events, such as compressor surge and inlet unstart, can be simulated with a high degree of fidelity. The simulation system was exercised using a supersonic inlet with sixty percent of the supersonic area contraction occurring internally, and a GE J85-13 turbojet engine.

Doug Garrard↗
NASA NTRS2000

Parametric Study of Reactive Melt Infiltration

Reactive melt infiltration is viewed as a promising means of achieving near-net shape manufacturing with quick processing time and at low cost. Since the reactants and products are, in general, of varying density, overall conservation of mass dictates that there is a force related to chemical conversion which can directly influence infiltration behavior. In effect, the driving pressure forces may compete with the forces from chemical conversion, affecting the advancement of the front. We have developed a two-dimensional numerical code to examine these effects, using reaction-formed silicon carbide as a model system for this process. We have examined a range of initial porosities, pore radii, and reaction rates in order to investigate their effects on infiltration dynamics.

Nelson, Emily S.↗
NASA NTRS1995

Analysis of Aerothermal Heating Over Reusable Single-Stage to Orbit Vehicles

Using the General Aerodynamic Simulation Program (GASP) with finite-rate surface catalysis and radiative equilibrium surface temperature conditions, a computational study has been conducted to predict three-dimensional aero-convective heating over the proposed reusable single-stage to orbit vehicles. The solutions for three basic classes of usable single-stage-to-orbit vehicles, including wing body, lifting body and vertical lander, are presented and compared to provide information for vehicle design. The free stream conditions used in this study occur at an altitude of 57.5 km, ambient density of 4.129e-4 kg/cubic m, and velocity of 4,849 m/s. The convective heating over the vehicle surface in both laminar and turbulent flow regimes were computed. The Baldwin-Lomax algebraic model with a modified van Driest damping factor was applied in the turbulent flow calculations. A bifurcation model which guarantees the conservation of mass diffusive fluxes was adopted to account for the multi-species mass diffusion. The effects of turbulence and surface catalysis on the surface heating are discussed in detail The solutions indicate the degree of oxygen dissociation over the leeward side is sensitive to the shape of vehicle body. Predicted convective heating over the entire surface of lifting body is strongly affected by the presence of surface catalysis. For wing body and vertical lander, the surface catalysis effect does not appear to be important in the far downstream region, and turbulence is the dominant effect on the heating prediction in this region.

Chen, Y. K.↗
NASA NTRS1999

1-D Numerical Analysis of ABCC Engine Performance

ABCC engine combines air breathing and rocket engine into a single engine to increase the specific impulse over an entire flight trajectory. Except for the heat source, the basic operation of the ABCC is similar to the basic operation of the RBCC engine. The ABCC is intended to have a higher specific impulse than the RBCC for single stage Earth to orbit vehicle. Computational fluid dynamics (CFD) is a useful tool for the analysis of complex transport processes in various components in ABCC propulsion system. The objective of the present research was to develop a transient 1-D numerical model using conservation of mass, linear momentum, and energy equations that could be used to predict flow behavior throughout a generic ABCC engine following a flight path. At specific points during the development of the 1-D numerical model a myriad of tests were performed to prove the program produced consistent, realistic numbers that follow compressible flow theory for various inlet conditions.

Holden, Richard↗
NASA NTRS2003

A Study of the Impact of Variations on Aerodynamic Flow in Gas Turbine Engines via Monte-Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented with a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, under the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions that motivate this research are: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance. This paper presents an initial attempt to the flow behavior in terms (perturbation) using simulation outputs from a customer-deck model obtained from Pratt&Whitney, (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine. Flow behavior is analyzed for the high pressure compressor and combustion chamber employing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000-sample simulation. Furthermore, the characteristic waves and eave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures; which are clearly the case in the combustion chamber and high-pressure turbine. Finally a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗
NASA NTRS2003

A Reduced-Order Model for Efficient Simulation of Synthetic Jet Actuators

A new reduced-order model of multidimensional synthetic jet actuators that combines the accuracy and conservation properties of full numerical simulation methods with the efficiency of simplified zero-order models is proposed. The multidimensional actuator is simulated by solving the time-dependent compressible quasi-1-D Euler equations, while the diaphragm is modeled as a moving boundary. The governing equations are approximated with a fourth-order finite difference scheme on a moving mesh such that one of the mesh boundaries coincides with the diaphragm. The reduced-order model of the actuator has several advantages. In contrast to the 3-D models, this approach provides conservation of mass, momentum, and energy. Furthermore, the new method is computationally much more efficient than the multidimensional Navier-Stokes simulation of the actuator cavity flow, while providing practically the same accuracy in the exterior flowfield. The most distinctive feature of the present model is its ability to predict the resonance characteristics of synthetic jet actuators; this is not practical when using the 3-D models because of the computational cost involved. Numerical results demonstrating the accuracy of the new reduced-order model and its limitations are presented.

Yamaleev, Nail K.↗
NASA NTRS2004

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗
NASA NTRS2004

Impact of Variations on 1-D Flow in Gas Turbine Engines via Monte Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented using a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, to determine the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions motivate this research: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance? This paper presents an initial attempt to model the flow behavior in terms of their eigenfrequencies subject to the assumption of the uncertainty or variation (perturbation). The flow behavior is explored using simulation outputs from a customer-deck model obtained from Pratt & Whitney. Variations of the main variables (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine are modeled. Flow behavior is analyzed for the high-pressure compressor and combustion chamber utilizing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high-pressure and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000 sample simulation. Furthermore, the characteristic waves and wave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures in the case of the combustion chamber and high-pressure turbine. Finally, a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗