SEARCH · Engineering Papers

Results for “finite volume method”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Numerical Modeling of Flow Distribution in Micro-Fluidics Systems

This paper describes an application of a general purpose computer program, GFSSP (Generalized Fluid System Simulation Program) for calculating flow distribution in a network of micro-channels. GFSSP employs a finite volume formulation of mass and momentum conservation equations in a network consisting of nodes and branches. Mass conservation equation is solved for pressures at the nodes while the momentum conservation equation is solved at the branches to calculate flowrate. The system of equations describing the fluid network is solved by a numerical method that is a combination of the Newton-Raphson and successive substitution methods. The numerical results have been compared with test data and detailed CFD (computational Fluid Dynamics) calculations. The agreement between test data and predictions is satisfactory. The discrepancies between the predictions and test data can be attributed to the frictional correlation which does not include the effect of surface tension or electro-kinetic effect.

Majumdar, Alok↗

NASA NTRS2007

Lunar Lander Structural Design Studies at NASA Langley

The National Aeronautics and Space Administration is currently developing mission architectures, vehicle concepts and flight hardware to support the planned human return to the Moon. During Phase II of the 2006 Lunar Lander Preparatory Study, a team from the Langley Research Center was tasked with developing and refining two proposed Lander concepts. The Descent-Assisted, Split Habitat Lander concept uses a disposable braking stage to perform the lunar orbit insertion maneuver and most of the descent from lunar orbit to the surface. The second concept, the Cargo Star Horizontal Lander, carries ascent loads along its longitudinal axis, and is then rotated in flight so that its main engines (mounted perpendicular to the vehicle longitudinal axis) are correctly oriented for lunar orbit insertion and a horizontal landing. Both Landers have separate crew transport volumes and habitats for surface operations, and allow placement of large cargo elements very close to the lunar surface. As part of this study, lightweight, efficient structural configurations for these spacecraft were proposed and evaluated. Vehicle structural configurations were first developed, and preliminary structural sizing was then performed using finite element-based methods. Results of selected structural design and trade studies performed during this activity are presented and discussed.

Wu, K. Chauncey↗

NASA NTRS2002

Analysis and Design of Crew Sleep Station for ISS

This paper details the analysis and design of the Temporary Sleep Station (TeSS) environmental control system for International Space Station (ISS). The TeSS will provide crewmembers with a private and personal space, to accommodate sleeping, donning and doffing of clothing, personal communication and performance of recreational activities. The need for privacy to accommodate these activities requires adequate ventilation inside the TeSS. This study considers whether temperature, carbon dioxide, and humidity within the TeSS remain within crew comfort and safety levels for various expected operating scenarios. Evaluation of these scenarios required the use and integration of various simulation codes. An approach was adapted for this study, whereby results from a particular code were integrated with other codes when necessary. Computational Fluid Dynamics (CFD) methods were used to evaluate the flow field inside the TeSS, from which local gradients for temperature, velocity, and species concentration such as CO (sub 2) could be determined. A model of the TeSS, containing a human, as well as equipment such as a laptop computer, was developed in FLUENT, a finite-volume code. Other factors, such as detailed analysis of the heat transfer through the structure, radiation, and air circulation from the TeSS to the US Laboratory Aisle, where the TeSS is housed, were considered in the model. A complementary model was developed in G189A, a code which has been used by NASA/JSC for environmental control systems analyses since the Apollo program. Boundary conditions were exchanged between the FLUENT and G189A TeSS models. G189A provides human respiration rates to the FLUENT model, while the FLUENT model provides local convective heat transfer coefficients to G189A model. An additional benefit from using an approach with both a systems simulation and CFD model, is the capability to verify the results of each model by comparison to the results of the other model. The G189A and FLUENT models were used to evaluate various ventilation designs for the TeSS over a range of operating conditions with varying crew metabolic load, equipment operating modes, ventilation flow rates, and with the TeSS doors open and closed. Results from the study were instrumental in the optimization of a design for the TeSS ventilation hardware. A special case was considered where failure of the TeSS ventilation system occurred. In this case, a study was conducted in order to determine the time required for the CO (sub 2) concentration inside the TeSS to increase to ISS limit values under transient conditions. A lumped-capacitance code, SINDA-FLUINT was used in this case to provide accurate predictions of the human reaction to the TeSS cabin conditions including core and skin temperatures and body heat storage. A simple two-dimensional CFD model of a crewmember inside the TeSS was developed in FLUENT in order to determine the volume envelope of the respired air from the human, which maintained a minimum velocity profile. This volume was then used in the SINDA-FLUINT model to facilitate the calculations of CO (sub 2) concentrations, dry bulb temperatures and humidity levels inside the TeSS.

Keener, John F.↗

NASA NTRS1987

Three-dimensional unsteady Euler solutions for propfans and counter-rotating propfans in transonic flow

An Euler code designed for computing the unsteady, three-dimensional, transonic flow about single-rotating and counter-rotating propfans using dynamic blocked-grids is presented. The algorithm is a finite volume, flux-split, upwind, implicit scheme and solves the equations which have been written in a time-dependent curvilinear coordinate system. Relative motion of the blades for counter-rotating configurations is handled by requiring that grid lines be aligned after each discrete rotation of fore and aft rotor grid blocks. The method by which information is passed across block interfaces, as well as how downstream characteristic outflow boundary conditions which enforce simple radial equilibrium are implemented, is discussed. Comparisons of computed flow-field parameters and propfan performance with experimental data indicate good overall agreement between predictions and measurements.

Whitfield, D. L.↗

NASA NTRS1991

Unified Aeroacoustics Analysis for High Speed Turboprop Aerodynamics and Noise: Development of Theory for Blade Loading, Wakes, and Noise - Volume 1

A unified theory for the aerodynamics and noise of advanced turboprops are presented. Aerodynamic topics include calculation of performance, blade load distribution, and non-uniform wake flow fields. Blade loading can be steady or unsteady due to fixed distortion, counter-rotating wakes, or blade vibration. The aerodynamic theory is based on the pressure potential method and is therefore basically linear. However, nonlinear effects associated with finite axial induction and blade vortex flow are included via approximate methods. Acoustic topics include radiation of noise caused by blade thickness, steady loading (including vortex lift), and unsteady loading. Shielding of the fuselage by its boundary layer and the wing are treated in separate analyses that are compatible but not integrated with the aeroacoustic theory for rotating blades.

Hanson, D. B.↗

NASA NTRS1982

Transonic potential flow and coordinate generation for bodies in a wind tunnel

An accurate method has been developed for computation of transonic potential flow about a 2-D lifting airfoil or an axisymmetric body in a wind tunnel. The computational mesh is nearly orthogonal everywhere and is generated by a sequence of Schwarz-Christoffel transformations and shearings to obtain an '0' grid near the body. A conservative finite-volume scheme for the full-potential equation and exact boundary conditions is used together with 'retarded' density to solve efficiently transonic flow with embedded shocks and large regions of supersonic flow. Supersonic free-stream flows are also solved with captured bow shocks and embedded subsonic regions.

Doria, M. L.↗

NASA NTRS1990

A user guide for the EMTAC-MZ CFD code

The computer code (EMTAC-MZ) was applied to investigate the flow field over a variety of very complex three-dimensional (3-D) configurations across the Mach number range (subsonic, transonic, supersonic, and hypersonic flow). In the code, a finite volume, multizone implementation of high accuracy, total variation diminishing (TVD) formulation (based on Roe's scheme) is used to solve the unsteady Euler equations. In the supersonic regions of the flow, an infinitely large time step and a space-marching scheme is employed. A finite time step and a relaxation or 3-D approximate factorization method is used in subsonic flow regions. The multizone technique allows very complicated configurations to be modeled without geometry modifications, and can easily handle combined internal and external flow problems. An elliptic grid generation package is built into the EMTAC-MZ code. To generate the computational grid, only the surface geometry data are required. Results obtained for a variety of configurations, such as fighter-like configurations (F-14, AVSTOL), flow through inlet, multi-bodies (shuttle with external tank and SRBs), are reported and shown to be in good agreement with available experimental data.

Szema, Kuo-Yen↗

NASA NTRS1992

An incremental strategy for calculating consistent discrete CFD sensitivity derivatives

In this preliminary study involving advanced computational fluid dynamic (CFD) codes, an incremental formulation, also known as the 'delta' or 'correction' form, is presented for solving the very large sparse systems of linear equations which are associated with aerodynamic sensitivity analysis. For typical problems in 2D, a direct solution method can be applied to these linear equations which are associated with aerodynamic sensitivity analysis. For typical problems in 2D, a direct solution method can be applied to these linear equations in either the standard or the incremental form, in which case the two are equivalent. Iterative methods appear to be needed for future 3D applications; however, because direct solver methods require much more computer memory than is currently available. Iterative methods for solving these equations in the standard form result in certain difficulties, such as ill-conditioning of the coefficient matrix, which can be overcome when these equations are cast in the incremental form; these and other benefits are discussed. The methodology is successfully implemented and tested in 2D using an upwind, cell-centered, finite volume formulation applied to the thin-layer Navier-Stokes equations. Results are presented for two laminar sample problems: (1) transonic flow through a double-throat nozzle; and (2) flow over an isolated airfoil.

Korivi, Vamshi Mohan↗

NASA NTRS2020

Implicit Large-Eddy Simulations of Hot and Cold Supersonic Jets in Loci-CHEM

Turbulence is a critical factor in most aero- and propulsion-related flows. For many years, the industry-standard approach to incorporating the effect of turbulence in computational fluid dynamics (CFD)calculations has been the Reynolds Averaged Navier Stokes (RANS) modeling approach. RANS implicitly time-averages turbulent motion, and models its effect on the mixing of species, momentum, and energy in a flow field. The approach is economical and often allows a CFD calculation to proceed to steady-state. However, it makes several simplifying assumptions, and the results will only be as good as the RANS turbulence model. Meanwhile, the last two decades have witnessed rapid and explosive growth in computational power, the development of CFD codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. CFD engineering calculations of 108– 109 finite-volume cells have become routine. Thus, large eddy simulation (LES), which attempts to directly resolve the unsteady motion of the largest scales of turbulence, is increasingly of interest for many fluids engineering problems. Additionally, hybrid RANS-LES approaches such as Spalart’s Detached Eddy Simulation (DES), which attempt to marry the best strengths of both methods, are in active use and development. Because it seeks to accurately resolve the unsteady motion of the largest scales of turbulence, LES places certain requirements on the numerical method used to discretize and advance the Navier Stokes equations. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide range of length scales.

C I Morris↗

NASA NTRS1991

Vapor condensation at the free surface of an axisymmetric liquid mixed by a laminar jet

This paper presents numerical solutions of jet-induced mixing in a partially full cryogenic tank. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to slove the nondimensional form of steady-state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of previously stated parameters on the condensation Nusselt and Stanton numbers that characterize the steady-state interface condensation process are investigated. Detailed analysis is performed to gain a better understanding of the fundamentals of fluid mixing and interface condensation.

Lin, Chin-Shun↗

NASA NTRS1993

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

In this study involving advanced fluid flow codes, an incremental iterative formulation (also known as the delta or correction form) together with the well-known spatially-split approximate factorization algorithm, is presented for solving the very large sparse systems of linear equations which are associated with aerodynamic sensitivity analysis. For smaller 2D problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. Iterative methods are needed for larger 2D and future 3D applications, however, because direct methods require much more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioning of the coefficient matrix; this problem can be overcome when these equations are cast in the incremental form. These and other benefits are discussed. The methodology is successfully implemented and tested in 2D using an upwind, cell-centered, finite volume formulation applied to the thin-layer Navier-Stokes equations. Results are presented for two sample airfoil problems: (1) subsonic low Reynolds number laminar flow; and (2) transonic high Reynolds number turbulent flow.

Taylor, Arthur C., III↗

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 NTRS2001

DTFM Modeling and Analysis Method for Gossamer Structures

Gossamer systems are mostly composed of support structures formed by highly flexible, long tubular elements and pre-tensioned thin-film membranes. These systems offer order-of-magnitude reductions in mass and launch volume and will revolutionize the architecture and design of space flight systems that require large in-orbit configurations and apertures. A great interest has been generated in recent years to fly gossamer systems on near-term and future space missions. Modeling and analysis requirements for gossamer structures are unique. Simulation of in-space performance issues of gossamer structures, such as inflation deployment of flexible booms, formation and effects of wrinkle in tensioned membranes, synthesis of tubular and membrane elements into a complete structural system, usually cannot be accomplished by using the general-purpose finite-element structural analysis codes. This has led to the need of structural modeling and analysis capabilities specifically suitable for gossamer structures. The Distributed Transfer Function Method (DTFM) can potentially meet this urgent need. Additional information is contained in the original extended abstract.

Fang, Hou-Fei↗

NASA NTRS1989

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Numerical solutions of jet-induced mixing in a partially full cryogenic tank are presented. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear-free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to solve the nondimensional form of steady state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of above stated parameters on the condensation Nusselt and Stanton numbers which characterize the steady-state interface condensation process are investigated. Detailed analysis to gain a better understanding of the fundamentals of fluid mixing and interface condensation is performed.

Lin, Chin-Shun↗

NASA NTRS1989

A multigrid and upwind viscous flow solver on 3-D embedded and overlapped grids

A numerically efficient method is presented for solving the three-dimensional governing equations of the viscous compressible flow about complex configurations with topologically different components. The physical domain is decomposed into regions for which the grid generation is relatively simple and virtually with no significant restrictions. The Navier-Stokes equations are solved by an implicit, approximately factored, upwind, finite-volume scheme. The block inversions and the diagonalized scalar inversions of the coefficient matrices are modified to allow the holes created in the computational domain by the embedded and overlapped grids. The convergence is accelerated by a multigrid algorithm despite the existence of such holes. The solution for s supersonic flow past a blunt-nose-cylinder at high angle-of-attack is obtained using a C-O grid embedded in a global Cartesian grid.

Baysal, Oktay↗

NASA NTRS1989

Numerical studies of the effects of jet-induced mixing on liquid-vapor interface condensation

Lin, Chin-Shun↗

NASA NTRS1993

Automatic partitioning of unstructured meshes for the parallel solution of problems in computational mechanics

Most of the recently proposed computational methods for solving partial differential equations on multiprocessor architectures stem from the 'divide and conquer' paradigm and involve some form of domain decomposition. For those methods which also require grids of points or patches of elements, it is often necessary to explicitly partition the underlying mesh, especially when working with local memory parallel processors. In this paper, a family of cost-effective algorithms for the automatic partitioning of arbitrary two- and three-dimensional finite element and finite difference meshes is presented and discussed in view of a domain decomposed solution procedure and parallel processing. The influence of the algorithmic aspects of a solution method (implicit/explicit computations), and the architectural specifics of a multiprocessor (SIMD/MIMD, startup/transmission time), on the design of a mesh partitioning algorithm are discussed. The impact of the partitioning strategy on load balancing, operation count, operator conditioning, rate of convergence and processor mapping is also addressed. Finally, the proposed mesh decomposition algorithms are demonstrated with realistic examples of finite element, finite volume, and finite difference meshes associated with the parallel solution of solid and fluid mechanics problems on the iPSC/2 and iPSC/860 multiprocessors.

Farhat, Charbel↗

NASA NTRS1993

Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach

Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.

Zuffada, C.↗