Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW MODELS”

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 1,747 records · Page 97

Idealized tip-to-tail waverider model

The flow field of an idealized cone-derived waverider is axisymmetric. This forebody feature is preserved for the rest of the vehicle, including the inlet, cowl, combustor, and nozzle. There is thus an inviscid, tip-to-tail model in which both the external and internal flows are axisymmetric. The assumption of axial symmetry provides a major simplification for the analysis and allows for a systematic integration of the propulsion unit with the aerodynamics. The code is an initial formulation that provides only the most basic engineering data, such as lift, thrust, drag, and fuel consumption for a point-designed vehicle that may be cruising at a low hypersonic Mach number. The user may specify flight altitude and Mach number, a multiple shock configuration for the inlet, a few basic geometric parameter, H2 or CH4 as fuel, the fuel/air ratio, etc. A new design concept is used for the nozzle that avoids shock waves, minimizes the nozzle length, and may maximize its thrust. A general description of the model is provided with emphasis on the design of the nozzle and fins. Preliminary results are presented that compare the cruise flight range using H2 or CH4 as the fuel.

Emanuel, G.↗

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(y) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(y) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(v) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(v) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗

A Flow Rate Dependent 1D Model for Thermally Stratified Hot-Water Energy Storage

Stratified tank models are used to simulate thermal storage in applications such as residential or commercial hot-water storage tanks, chilled-water storage tanks, and solar thermal systems. The energy efficiency of these applications relates to the system components and the level of stratification maintained during various flow events in the tank. One-dimensional (1D) models are used in building energy simulations because of the short computation time but often do not include flow-rate dependent mixing. The accuracy of 1D models for plug flow, plug flow with axial conduction, and two convection eddy-diffusivity models were compared with experimental data sets for discharging a 50-gal residential tank and recharging the tank with hot water from an external hot-water source. A minimum and maximum relationship for the eddy diffusivity factor were found at Re <2100 and >10,000 for recirculation of hot water to the top of the tank and vertical tubes inletting cold water at the bottom. The root mean square error decreased from >4 °C to near 2 °C when considering flow-based mixing models during heating, while the exponential decay of the eddy diffusion results in a root mean square error reduction of 1 °C for cone-shaped diffusers that begin to relaminarize flow at the inlet.

1D flow model↗

Modeling of supersonic combustor flows using parallel computing

While current 3D CFD codes and modeling techniques have been shown capable of furnishing engineering data for complex scramjet flowfields, the usefulness of such efforts is primarily limited by solutions' CPU time requirements, and secondarily by memory requirements. Attention is presently given to the use of parallel computing capabilities for engineering CFD tools for the analysis of supersonic reacting flows, and to an illustrative incompressible CFD problem using up to 16 iPSC/2 processors with single-domain decomposition.

Riggins, D.↗

Application of a new K-tau model to near wall turbulent flows

A recently developed K-tau model for near wall turbulent flows is applied to two severe test cases. The turbulent flows considered include the incompressible flat plate boundary layer with the adverse pressure gradients and incompressible flow past a backward facing step. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.

Thangam, S.↗

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗

Finite element modeling of biomass hopper flow

Hoppers are widely used biomass handling devices that channel bulk biomass from storage to subsequent handling equipment. Jenike’s longstanding approach, based on the Mohr-Coulomb model, has been successfully used to design hoppers handling cohesionless granular materials such as grains and other agricultural produces. However, designing a hopper to ensure reliable biomass flow is found to be challenging due to cohesion, irregular particle shape, and bulk material elastoplasticity. This study aims to address the biomass handling engineering challenge with alternative constitutive material models concerning the flow behavior of bulk solids. Finite element modeling is an approach that allows for implementation of different material models, whose underlying constitutive theories assist in investigating the origin and manifestation of bulk mechanical behavior of granular materials. This study focuses on the incipient gravity hopper flow of two types of biomass feedstocks, i.e., ground corn stover and Douglas fir wood. Three widely used constitutive material models, i.e., Mohr-Coulomb model, modified Cam-Clay model, and Drucker-Prager/Cap model, are implemented. Using the flow pattern represented by the volume of biomass exhibiting more than 7% of axial strain (Kamath and Puri, 1999), the finite element model predicts that the bulk corn stover particulate material forms an arch, which represents a hampered transition from the static state to the dynamic flow-state out of the hopper, whereas bulk Douglas fir wood particulate material develops a reliable mass flow pattern. A laboratory scale hopper was used to experimentally determine the biomass flow conditions, which were subsequently compared with the predicted onset of flow by a finite element model (FEM). The developed FEM was found to correctly predict the initiation of mass flow for the milled Douglas fir wood, whereas corn stover was predicted to establish a strong core flow suggesting an unreliable handling characteristic. This observation aligns with the reported poor handling of milled corn stover.

09 BIOMASS FUELS↗

Numerical study of three-dimensional turbulent flow interactions between blockage models and wind tunnels including longitudinally slotted test sections

A spatial forward-marching approach is applied to compute three-dimensional turbulent flows for several blockage models in free flight, in a solid wall wind tunnel, and in a wind tunnel with longitudinal slots in the test section. The effects of area blockage in the tunnel, model growth and tunnel wall boundary layers, and of the slots are included. The large blockage models are found to have significant wall interference effects which can be reduced by the slots. The effects of the latter are confined to the region near the tunnel wall. Model/wall interference effects are not limited to the effects of area blockage; in particular, boundary layer profile shapes for a wind tunnel model in a tunnel are different from shapes for a model in free flight even when slots are used. This indicates that the flow responds differently in these two cases with the same pressure gradient.

Lin, S.-J.↗

A Multi-Dimensional Heat Transfer Model of a Tie-Tube and Hexagonal Fuel Element for Nuclear Thermal Propulsion

The Space Capable Cryogenic Thermal Engine (SCCTE) effort considers a nuclear thermal rocket design based around a Low-Enriched Uranium (LEU) design fission reactor. The reactor core is comprised of bundled hexagonal fuel elements that directly heat hydrogen for expansion in a thrust chamber and hexagonal tie-tubes that house zirconium hydride moderator mass for the purpose of thermalizing fast neutrons resulting from fission events. Created 3D steady state Hex fuel rod model with 1D flow channels. Hand Calculation were used to set up initial conditions for fluid flow. The Hex Fuel rod uses 1D flow paths to model the channels using empirical correlations for heat transfer in a pipe. Created a 2-D axisymmetric transient to steady state model using the CFD turbulent flow and Heat Transfer module in COMSOL. This model was developed to find and understand the hydrogen flow that might effect the thermal gradients axially and at the end of the tie tube where the flow turns and enters an annulus. The Hex fuel rod and Tie tube models were made based on requirements given to us by CSNR and the SCCTE team. The models helped simplify and understand the physics and assumptions. Using pipe correlations reduced the complexity of the 3-D fuel rod model and is numerically more stable and computationally more time-efficient compared to the CFD approach. The 2-D axisymmetric tie tube model can be used as a reference "Virtual test model" for comparing and improving 3-D Models.

Gomez, C. F.↗

Comparison of liquid rocket engine base region heat flux computations using three turbulence models

The flow in the base region of launch vehicles is characterized by flow separation, flow reversals, and reattachment. Computation of the convective heat flux in the base region and on the nozzle external surface of Space Shuttle Main Engine and Space Transportation Main Engine (STME) is an important part of defining base region thermal environments. Several turbulence models were incorporated in a CFD code and validated for flow and heat transfer computations in the separated and reattaching regions associated with subsonic and supersonic flows over backward facing steps. Heat flux computations in the base region of a single STME engine and a single S1C engine were performed using three different wall functions as well as a renormalization-group based k-epsilon model. With the very limited data available, the computed values are seen to be of the right order of magnitude. Based on the validation comparisons, it is concluded that all the turbulence models studied have predicted the reattachment location and the velocity profiles at various axial stations downstream of the step very well.

Kumar, Ganesh N.↗

Turbulence modeling in aerodynamic shear flows - Status and problems

This paper briefly summarizes the status and problems of turbulence modeling for aerodynamical applications. For complex flows the 'approach of choice' is (increasingly) full second-order (Reynolds stress equation) closure. These closures have not yet developed to anywhere near their full potential, significant further research is required especially regarding length-scale equations, representation of pressure-strain correlations, and wall region treatments. Recent developments in computer capability, algorithms, numerical simulations, theory and quantitative flow visualization should assist in and hasten this research. Several problem areas such as shock interaction and discrete dynamic instabilities of turbulent flows may require mega-to-large eddy simulation or theoretical adjuncts.

Bushnell, D. M.↗

Transition of oscillatory flow in tubes - An empirical model for application to Stirling engines

Attention is given to an empirical model for transition to turbulence in oscillatory flows in straight tubes. Designed after a correlation for transition of a boundary layer on a flat plate, the model yields the laminar flow momentum thickness Reynolds number that must be met before transition to turbulence will occur. The transition point is located by comparing this to the actual momentum thickness Reynolds number. A scheme is proposed for estimating the momentum thickness Reynolds number in terms of the position within the cycle, the maximum value of the diameter Reynolds within the cycle, Re(max), and the dimensionless frequency, Valensi number. Results from an experimental study of oscillatory flow in a tube are employed to develop the model. When the flow is determined to be turbulent, it is proposed that a fully-developed, steady flow friction coefficient be applied. When the flow is laminar, the assumption of fully developed flow cannot be made; thus, a method is suggested for estimating the friction factor.

Simon, T. W.↗

A coupled thermo-hydro-mechanical model for simulating leakoff-dominated hydraulic fracturing with application to geologic carbon storage

A potential risk of injecting CO2 into storage reservoirs with marginal permeability (≲ 10 mD (1 mD = 10 –15 m 2 )) is that commercial injection rates could induce fracturing of the reservoir and/or the caprock. Such fracturing is essentially fluid-driven fracturing in the leakoff-dominated regime. Recent studies suggested that fracturing, if contained within the lower portion of the caprock complex, could substantially improve the injectivity without compromising the overall seal integrity. Modeling this phenomenon entails complex coupled interactions among the fluids, the fracture, the reservoir, and the caprock. Here, we develop a simple method to capture all these interplays in high fidelity by sequentially coupling a hydraulic fracturing module with a coupled thermal-hydrological-mechanical (THM) model for nonisothermal multiphase flow. The model was made numerically tractable by taking advantage of self-stabilizing features of leakoff-dominated fracturing. The model is validated against the PKN solution in the leakoff-dominated regime. Moreover, we employ the model to study thermo-poromechanical responses of a fluid-driven fracture in a field-scale carbon storage reservoir that is loosely based on the In Salah project's Krechba reservoir. The model reveals complex yet intriguing behaviors of the reservoir-caprock-fluid system with fracturing induced by cold CO 2 injection. We also study the effects of the in situ stress contrast between the reservoir and caprock and thermal contraction on the vertical containment of the fracture. The proposed model proves effective in simulating practical problems on length and time scales relevant to geological carbon storage.

58 GEOSCIENCES↗

Helicopter rotor noise due to ingestion of atmospheric turbulence

A theoretical study was conducted to develop an analytical prediction method for helicopter main rotor noise due to the ingestion of atmospheric turbulence. This study incorporates an atmospheric turbulence model, a rotor mean flow contraction model and a rapid distortion turbulence model which together determine the statistics of the non-isotropic turbulence at the rotor plane. Inputs to the combined mean inflow and turbulence models are controlled by atmospheric wind characteristics and helicopter operating conditions. A generalized acoustic source model was used to predict the far field noise generated by the non-isotropic flow incident on the rotor. Absolute levels for acoustic spectra and directivity patterns were calculated for full scale helicopters, without the use of empirical or adjustable constants. Comparisons between isotropic and non-isotropic turbulence at the rotor face demonstrated pronounced differences in acoustic spectra. Turning and contraction of the flow for hover and low speed vertical ascent cases result in a 3 dB increase in the acoustic spectrum energy and a 10 dB increase in tone levels. Compared to trailing edge noise, turbulence ingestion noise is the dominant noise mechanism below approximately 30 rotor harmonics, while above 100 harmonics, trailing edge noise levels exceed turbulence ingestion noise by 25 dB.

Simonich, J. C.↗

Scalar and momentum turbulent transport experiments with swirling and nonswirling flows

Combustor models for the aircraft gas turbine industry have been obtained because of the need to reduce the costs of developing improved performance and more durable engines. A few years ago, it became apparent that the mass concentration and velocity predictions provided by the computer codes were not representing the data measured in some confined recirculating flows. It is pointed out that errors in the mass concentration distribution are an especially serious problem because of their influence on the heat release, temperature, and reactant distributions. Combined mass and momentum turbulent transport experiments with swirling and nonswirling flow have been conducted with the objective to obtain an experimental data base which can be used to evaluate and improve the turbulent transport submodes employed in the aerothermal models. The present paper is mainly concerned with the overall characteristics of the mass turbulent transport processes in complex flows with recirculation and the deficiencies of the conventional models.

Johnson, B. V.↗