SEARCH · Engineering Papers

Results for “Finite Volume”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Determination of U, V, and W from single station Doppler radar radial velocities

The ST/MST (stratosphere troposphere/mesosphere stratosphere troposphere) clear air Doppler radar, or wind profiler, is an important tool in observational meteorology because of its capability to remote observe dynamic parameters of the atmosphere. There are difficulties in transforming the observed radial velocities into meteorological wind components. How this problem has been treated in the past is reviewed, and some of the analysis is recast to a form more suited to the high diagnostic abilities of a number of fixed beam configurations with reference to a linear wind field. The results, in conjunction with other works which treats problems such as the effects of finite sample volumes in the presence of nonhomogeneous atmospheric reflectivity, have implications important to the design of both individual MST/ST radars and MST/ST radar networks. The key parameters to uncoupling terms in the scaling equations are w sub x and w sub y. Whenever the stratiform condition, which states that these two parameters are negligible, is satisfied, a five beam ST radar may determine unbiased values of u, v, and w for sample volumes directly above the radar. The divergence and partial deformation of the flow may also be determined. Three beam systems can determine w and w sub z, but are unable to obtain u and v wind components uncontaminated by vertical sheer terms, even when the stratiform condition is satisfied.

Clark, W. L.↗

NASA NTRS1989

Control-volume based Navier-Stokes equation solver valid at all flow velocities

A control-volume based finite difference method to solve the Reynolds averaged Navier-Stokes equations is presented. A pressure correction equation valid at all flow velocities and a pressure staggered grid layout are used in the method. Example problems presented herein include: a developing laminar channel flow, developing laminar pipe flow, a lid-driven square cavity flow, a laminar flow through a 90-degree bent channel, a laminar polar cavity flow, and a turbulent supersonic flow over a compression ramp. A k-epsilon turbulence model supplemented with a near-wall turbulence model was used to solve the turbulent flow. It is shown that the method yields accurate computational results even when highly skewed, unequally spaced, curved grids are used. It is also shown that the method is strongly convergent for high Reynolds number flows.

Kim, S.-W.↗

NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗

NASA NTRS1989

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Numerical calculations of turbulent reattaching shear layers in a divergent channel are presented. The turbulence is described by a multiple-time-scale turbulence model. The turbulent flow equations are solved by a control-volume based finite difference method. The computational results are compared with those obtained using k-epsilon turbulence models and algebraic Reynolds stress turbulence models. It is shown that the multiple-time-scale turbulence model yields significantly improved computational results than the other turbulence models in the region where the turbulence is in a strongly inequilibrium state.

Kim, S.-W.↗

NASA NTRS1989

Numerical investigation of an internal layer in turbulent flow over a curved hill

The development of an internal layer in a turbulent boundary layer flow over a curved hill is investigated numerically. The turbulence field of the boundary layer flow over the curved hill is compared with that of a turbulent flow over a symmetric airfoil (which has the same geometry as the curved hill except that the leading and trailing edge plates were removed) to study the influence of the strongly curved surface on the turbulence field. The turbulent flow equations are solved by a control-volume based finite difference method. The turbulence is described by a multiple-time-scale turbulence model supplemented with a near-wall turbulence model. Computational results for the mean flow field (pressure distributions on the walls, wall shearing stresses and mean velocity profiles), the turbulence structure (Reynolds stress and turbulent kinetic energy profiles), and the integral parameters (displacement and momentum thicknesses) compared favorably with the measured data. Computational results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. Development of the internal layer was more obviously observed in the Reynolds stress profiles and in the turbulent kinetic energy profiles than in the mean velocity profiles. In this regard, the internal layers is significantly different from wall-bounded simple shear layers in which the mean velocity profile characterizes the boundary layer most distinguishably. Development of such an internal layer, characterized by an intense turbulence field, is attributed to the enormous mean flow strain rate caused by the streamline curvature and the strong pressure gradient. In the turbulent flow over the curved hill, the internal layer begin to form near the forward corner of the hill, merges with the external boundary layer, and develops into a new fully turbulent boundary layer as the fluid flows in the downstream direction. For the flow over the symmetric airfoil, the boundary layer began to form from almost the same location as that of the curved hill, grew in its strength, and formed a fully turbulent boundary layer from mid-part of the airfoil and in the downstream region. Computational results also show that the detailed turbulence structure in the region very close to the wall of the curved hill is almost the same as that of the airfoil in most of the curved regions except near the leading edge. Thus the internal layer of the curved hill and the boundary layer of the airfoil were also almost the same. Development of the wall shearing stress and separation of the boundary layer at the rear end of the curved hill mostly depends on the internal layer and is only slightly influenced by the external boundary layer flow.

Kim, S-W.↗

NASA NTRS1989

A numerical analysis of Stefan problems for generalized multi-dimensional phase-change structures using the enthalpy transforming model

An enthalpy transforming scheme is proposed to convert the energy equation into a nonlinear equation with the enthalpy, E, being the single dependent variable. The existing control-volume finite-difference approach is modified so it can be applied to the numerical performance of Stefan problems. The model is tested by applying it to a three-dimensional freezing problem. The numerical results are in agreement with those existing in the literature. The model and its algorithm are further applied to a three-dimensional moving heat source problem showing that the methodology is capable of handling complicated phase-change problems with fixed grids.

Cao, Yiding↗

NASA NTRS1990

Calculation of reattaching shear layers in divergent channel with a multiple-time-scale turbulence model

Kim, S.-W.↗

NASA NTRS1990

Thermal buoyancy and Marangoni convection in a two fluid layered system - A numerical study

Thermal buoyancy and surface tension driven convection is numerically investigated in a system with two, immiscible fluids. The geometry investigated, has an open cavity configuration with the lighter fluid situated on top of the heavier fluid, forming a stable layered system. The upper fluid meniscus and the interface are assumed to be flat and undeformable in the calculations. The governing equations and boundary/interface conditions are solved by a control volume based finite difference scheme for two pairs of immiscible fluids; the water-hexadecane system and a so called generic system. The steady state calculations show, that dramatically different flows are predicted when the interfacial tension effects are included or excluded from the system model. These differences are particularly appreciable in surface tension dominated flows, that are typical of microgravity situations. Complex flow patterns, with induced secondary flows are noticed in both the fluids. In general, the overall system heat transfer is found to increase with increases in the thermal buoyancy and surface tension effects, but the behavior of the system flow and thermal fields is not easy to characterize, when different combinations of these forces are considered.

Ramachandran, N.↗

NASA NTRS1990

Small domes on Venus - Characteristics and origin

The areal and size-frequency distribution, abundance, individual characteristics, and geologic and regional associations of small domes in the Venera data set are examined. It is noted that, because of their numbers and widespread occurence, regardless of origin, these structures will be important in the local and global geographic interpretation of the surface. It is concluded that small domes of 20 km in diameter occur in numbers of order 10,000 on the northern quarter of the surface of Venus, and that small domes are interpreted to be predominantly low shield volcanos and to represent multiple centralized effusive eruptions of discrete volumes over finite periods. It is noted that small domes occur throughout the northern latitudes, but attain maximum regional concentrations in two areas, one located north of Thetis Regio in Akkruva Colles, and another located on the northeast flanks of Beta Regio in Guinevere Planitia. The presence of two global concentrations implies that magma production and eruption were enhanced within these two areas.

Aubele, Jayne C.↗

NASA NTRS1991

Slave finite elements for nonlinear analysis of engine structures, volume 1

A 336 degrees of freedom slave finite element processing capability to analyze engine structures under severe thermomechanical loading is presented. Description of the theoretical development and demonstration of that element is presented in this volume.

Gellin, S.↗

NASA NTRS1992

Computation of non-equilibrium flow downstream of a plasma torch

Numerical solutions of the Navier-Stokes equations for compressible reacting flow in an axisymmetric geometry are presented for a nitrogen plasma torch with both thermal and chemical nonequilibrium. The Navier-Stokes equations are solved using a new axisymmetric finite element/finite volume formulation in which the convective flux is treated by an extension of the approximate Riemann solver due to Osher. The numerical scheme is validated by comparison with a previous solution of the same problem using a different computational scheme. Results obtained using two different models of nonthermal dissociation rates are compared to experimental data.

Harle, Christophe↗

NASA NTRS1992

Two-dimensional thermal modeling of power monolithic microwave integrated circuits (MMIC's)

Numerical simulations of the two-dimensional temperature distributions for a typical GaAs MMIC circuit are conducted, aiming at understanding the heat conduction process of the circuit chip and providing temperature information for device reliability analysis. The method used is to solve the two-dimensional heat conduction equation with a control-volume-based finite difference scheme. In particular, the effects of the power dissipation and the ambient temperature are examined, and the criterion for the worst operating environment is discussed in terms of the allowed highest device junction temperature.

Fan, Mark S.↗

NASA NTRS1993

Verification of a two-dimensional infiltration model for the resin transfer molding process

A two-dimensional finite element model for the infiltration of a dry textile preform by an injected resin was verified. The model, which is based on the finite element/control volume technique, determines the total infiltration time and the pressure increase at the mold inlet associated with the RTM process. Important input data for the model are the compaction and permeability behavior of the preform along with the kinetic and rheological behavior of the resin. The compaction behavior for several textile preforms was determined by experimental methods. A power law regression model was used to relate fiber volume fraction to the applied compaction pressure. Results showed a large increase in fiber volume fraction with the initial application of pressure. However, as the maximum fiber volume fraction was approached, the amount of compaction pressure required to decrease the porosity of the preform rapidly increased. Similarly, a power law regression model was used to relate permeability to the fiber volume fraction of the preform. Two methods were used to measure the permeability of the textile preform. The first, known as the steady state method, measures the permeability of a saturated preform under constant flow rate conditions. The second, denoted the advancing front method, determines the permeability of a dry preform to an infiltrating fluid. Water, corn oil, and an epoxy resin, Epon 815, were used to determine the effect of fluid type and viscosity on the steady state permeability behavior of the preform. Permeability values measured with the different fluids showed that fluid viscosity had no influence on the permeability behavior of 162 E-glass and TTI IM7/8HS preforms. Permeabilities measured from steady state and advancing front experiments for the warp direction of 162 E-glass fabric were similar. This behavior was noticed for tests conducted with corn oil and Epon 815. Comparable behavior was observed for the warp direction of the TTI IM7/8HS preform and corn oil. Mold filling and flow visualization experiments were performed to verify the analytical computer model. Frequency dependent electromagnetic sensors were used to monitor the resin flow front as a function of time. For the flow visualization tests, a video camera and high resolution tape recorder were used to record the experimental flow fronts. Comparisons between experimental and model predicted flow fronts agreed well for all tests. For the mold filling tests conducted at constant flow rate injection, the model was able to accurately predict the pressure increase at the mold inlet during the infiltration process. A kinetics model developed to predict the degree of cure as a function of time for the injected resin accurately calculated the increase in the degree of cure during the subsequent cure cycle.

Hammond, Vincent H.↗

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 NTRS1996

Minimum Heating Re-Entry Trajectories for Advanced Hypersonic Launch Vehicles

Optimal re-entry trajectories are generated for reusable launch vehicles which minimize: (1) the heat absorbed at the vehicle surface, (2) the lower surface temperature, and (3) the heat absorbed by the internal structure. The approach uses the energy state approximation technique and a finite control volume heat transfer code coupled to a flight path integration code. These trajectories are compared to the optimal re-entry trajectory minimizing the integrated convective heat rate to determine which trajectory produces the minimum internal structural temperatures for a given thermal protection system. Three different thermal protection systems are considered: tile, blanket, and metallic.

Robert Windhorst↗

NASA NTRS2000

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗

NASA NTRS2011

Extension of Generalized Fluid System Simulation Program's Fluid Property Database

This internship focused on the development of additional capabilities for the General Fluid Systems Simulation Program (GFSSP). GFSSP is a thermo-fluid code used to evaluate system performance by a finite volume-based network analysis method. The program was developed primarily to analyze the complex internal flow of propulsion systems and is capable of solving many problems related to thermodynamics and fluid mechanics. GFSSP is integrated with thermodynamic programs that provide fluid properties for sub-cooled, superheated, and saturation states. For fluids that are not included in the thermodynamic property program, look-up property tables can be provided. The look-up property tables of the current release version can only handle sub-cooled and superheated states. The primary purpose of the internship was to extend the look-up tables to handle saturated states. This involves a) generation of a property table using REFPROP, a thermodynamic property program that is widely used, and b) modifications of the Fortran source code to read in an additional property table containing saturation data for both saturated liquid and saturated vapor states. Also, a method was implemented to calculate the thermodynamic properties of user-fluids within the saturation region, given values of pressure and enthalpy. These additions required new code to be written, and older code had to be adjusted to accommodate the new capabilities. Ultimately, the changes will lead to the incorporation of this new capability in future versions of GFSSP. This paper describes the development and validation of the new capability.

Patel, Kishan↗

NASA NTRS2011

Coherent Backscattering by Polydisperse Discrete Random Media: Exact T-Matrix Results

The numerically exact superposition T-matrix method is used to compute, for the first time to our knowledge, electromagnetic scattering by finite spherical volumes composed of polydisperse mixtures of spherical particles with different size parameters or different refractive indices. The backscattering patterns calculated in the far-field zone of the polydisperse multiparticle volumes reveal unequivocally the classical manifestations of the effect of weak localization of electromagnetic waves in discrete random media, thereby corroborating the universal interference nature of coherent backscattering. The polarization opposition effect is shown to be the least robust manifestation of weak localization fading away with increasing particle size parameter.

Mishchenko, Michael I.↗