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,189 records · Page 66

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Numerically Exact Computer Simulations of Light Scattering by Densely Packed, Random Particulate Media

Direct computer simulations of electromagnetic scattering by discrete random media have become an active area of research. In this progress review, we summarize and analyze our main results obtained by means of numerically exact computer solutions of the macroscopic Maxwell equations. We consider finite scattering volumes with size parameters in the range, composed of varying numbers of randomly distributed particles with different refractive indices. The main objective of our analysis is to examine whether all backscattering effects predicted by the low-density theory of coherent backscattering (CB) also take place in the case of densely packed media. Based on our extensive numerical data we arrive at the following conclusions: (i) all backscattering effects predicted by the asymptotic theory of CB can also take place in the case of densely packed media; (ii) in the case of very large particle packing density, scattering characteristics of discrete random media can exhibit behavior not predicted by the low-density theories of CB and radiative transfer; (iii) increasing the absorptivity of the constituent particles can either enhance or suppress typical manifestations of CB depending on the particle packing density and the real part of the refractive index. Our numerical data strongly suggest that spectacular backscattering effects identified in laboratory experiments and observed for a class of high-albedo Solar System objects are caused by CB.

Dlugach, Janna M.↗

NASA NTRS2013

Ablator Response Model Development: From Flight Data Back to Fundamental Experiments

The successful Mars atmospheric entry by the Mars Science Laboratory (MSL-Curiosity) combined with the success of the Earth atmospheric entry by the Stardust capsule have established PICA as a major Thermal Protection Systems (TPS) material. We expect that this class of materials will be on the short list selected by NASA for any atmospheric entry missions and that it will be the lead of that list of materials in any planning, feasibility studies or flight readiness studies. In addition to NASAs successes, the Dragon capsule, the successful commercial space vehicle built by SpaceX, uses PICA-X, while the European Space Agency is considering ASTERM for its exploration missions that involve atmospheric entries, both of these materials are of the same family as PICA. In the talk, a high-fidelity model will be detailed and discussed. The model tracks the chemical composition of the gases produced during pyrolysis. As in the conventional models, it uses equilibrium chemistry to determine the recession rate at high temperatures but switches to in-volume finite-rate ablation for lower temperatures. It also tracks the time evolution of the porosity of the material. Progress in implementing this high-fidelity model in a code will be presented. In addition, a set of basic experimental data being supported for model validation will be summarized. The validation process for the model development will be discussed. Preliminary results will be presented for a case where detailed pyrolysis product chemistry is computed. Finally, a wish list for a set of validation experiments will be outlined and discussed.

Mansour, Nagi N.↗

NASA NTRS2013

MHD Simulations of the Plasma Flow in the Magnetic Nozzle

The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.

Smith, T. E. R.↗

NASA NTRS2016

Modeling and Simulation of Radiative Compressible Flows in Aerodynamic Heating Arc-Jet Facility

Numerical simulations of an arc heated flow inside NASA's 20 [MW] Aerodynamics heating facility (AHF) are performed in order to investigate the three-dimensional swirling flow and the current distribution inside the wind tunnel. The plasma is considered in Local Thermodynamics Equilibrium(LTE) and is composed of Air-Argon gas mixture. The governing equations are the Navier-Stokes equations that include source terms corresponding to Joule heating and radiative cooling. The former is obtained by solving an electric potential equation, while the latter is calculated using an innovative massively parallel ray-tracing algorithm. The fully coupled system is closed by the thermodynamics relations and transport properties which are obtained from Chapman-Enskog method. A novel strategy was developed in order to enable the flow solver and the radiation calculation to be preformed independently and simultaneously using a different number of processors. Drastic reduction in the computational cost was achieved using this strategy. Details on the numerical methods used for space discretization, time integration and ray-tracing algorithm will be presented. The effect of the radiative cooling on the dynamics of the flow will be investigated. The complete set of equations were implemented within the COOLFluiD Framework. Fig. 1 shows the geometry of the Anode and part of the constrictor of the Aerodynamics heating facility (AHF). Fig. 2 shows the velocity field distribution along (x-y) plane and the streamline in (z-y) plane.

Pasma flows↗

NASA NTRSDate not supplied

Energy of Fluid (EOF) Simulation of Large-Scale Densification and Solidification of Hydrogen

The current research develops a finite volume-based Computation Fluid Dynamic (CFD) model utilizing an Energy of Fluid (EOF) approach to generate a simulation of the densification of liquid hydrogen (LH2) in the Integrated Refrigeration and Storage (IRAS) experimental tank for the Ground Operations Demonstration Unit for Liquid Hydrogen Project (GODU-LH2) at the Kennedy Space Center (KSC). The computational code will integrate a commercial software pressure-based model with User Defined Functions (UDF) generated and added to the model. The modified model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. This method will decrease computational time when compared to the built-in commercially available Enthalpy solvers. Currently the research focused on utilizing User Defined Scalars (UDS) compared to the enthalpy model and a lumped node analysis to show the written UDS is as accurate as the built-in solvers, and the models will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

C. Mahony↗

NASA NTRSDate not supplied

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

NASA NTRS2023

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

Nodal Model↗

NASA NTRS2023

Simulations of Zero Boiloff, Densification, and Solidification in a Large-Scale Liquid Hydrogen Tank

The current research utilizes an Energy of Fluid (EOF) approach to develop a finite volume-based Computation Fluid Dynamic (CFD) model to create a benchmark simulation of the densification of liquid hydrogen (LH2) in an experimental Integrated Refrigeration and Storage (IRAS) tank at the Kennedy Space Center (KSC). The computational code will incorporate the FLUENT pressure-based model with User Defined Functions (UDF) to implement the EOF model. The enhanced model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. Specifically, a convection term will be added to an already existing code for a more accurate simulation of LH2 densification. However, currently the research focuses on utilizing User Defined Scalars (UDS) with FLUENT to show the built-in solvers are as accurate as the written UDS. This will reduce the computational time down significantly by utilizing an already built in solver for the momentum equation. A custom UDF has been written to solve transient and steady state conduction, and results obtained using the UDF agree with results obtained using the built-in FLUENT solver. Currently, a convection UDF is being added and will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

Colin Philip Mahony↗

NASA NTRS1975

PLANS: A finite element program for nonlinear analysis of structures. Volume 1: Theoretical manual

The PLANS system is described which is a finite element program for nonlinear analysis. The system represents a collection of special purpose computer programs each associated with a distinct physical problem class. Modules of PLANS specifically referenced and described in detail include: (1) REVBY, for the plastic analysis of bodies of revolution; (2) OUT-OF-PLANE, for the plastic analysis of 3-D built-up structures where membrane effects are predominant; (3) BEND, for the plastic analysis of built-up structures where bending and membrane effects are significant; (4) HEX, for the 3-D elastic-plastic analysis of general solids; and (5) OUT-OF-PLANE-MG, for material and geometrically nonlinear analysis of built-up structures. The SATELLITE program for data debugging and plotting of input geometries is also described. The theoretical foundations upon which the analysis is based are presented. Discussed are the form of the governing equations, the methods of solution, plasticity theories available, a general system description and flow of the programs, and the elements available for use.

Pifko, A.↗

NASA NTRS1999

High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.

Barth, Timothy↗

NASA NTRS1977

PLANS; a finite element program for nonlinear analysis of structures. Volume 2: User's manual

The PLANS system, rather than being one comprehensive computer program, is a collection of finite element programs used for the nonlinear analysis of structures. This collection of programs evolved and is based on the organizational philosophy in which classes of analyses are treated individually based on the physical problem class to be analyzed. Each of the independent finite element computer programs of PLANS, with an associated element library, can be individually loaded and used to solve the problem class of interest. A number of programs have been developed for material nonlinear behavior alone and for combined geometric and material nonlinear behavior. The usage, capabilities, and element libraries of the current programs include: (1) plastic analysis of built-up structures where bending and membrane effects are significant, (2) three dimensional elastic-plastic analysis, (3) plastic analysis of bodies of revolution, and (4) material and geometric nonlinear analysis of built-up structures.

Pifko, A.↗

NASA NTRS2011

Material Characterization and Geometric Segmentation of a Composite Structure Using Microfocus X-Ray Computed Tomography Image-Based Finite Element Modeling

This study utilizes microfocus x-ray computed tomography (CT) slice sets to model and characterize the damage locations and sizes in thermal protection system materials that underwent impact testing. ScanIP/FE software is used to visualize and process the slice sets, followed by mesh generation on the segmented volumetric rendering. Then, the local stress fields around several of the damaged regions are calculated for realistic mission profiles that subject the sample to extreme temperature and other severe environmental conditions. The resulting stress fields are used to quantify damage severity and make an assessment as to whether damage that did not penetrate to the base material can still result in catastrophic failure of the structure. It is expected that this study will demonstrate that finite element modeling based on an accurate three-dimensional rendered model from a series of CT slices is an essential tool to quantify the internal macroscopic defects and damage of a complex system made out of thermal protection material. Results obtained showing details of segmented images; three-dimensional volume-rendered models, finite element meshes generated, and the resulting thermomechanical stress state due to impact loading for the material are presented and discussed. Further, this study is conducted to exhibit certain high-caliber capabilities that the nondestructive evaluation (NDE) group at NASA Glenn Research Center can offer to assist in assessing the structural durability of such highly specialized materials so improvements in their performance and capacities to handle harsh operating conditions can be made.

Abdul-Aziz, Ali↗

NASA NTRS1991

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗

NASA NTRS1977

Multiple-compartment venting

Program simulates venting characteristics of multireservoir multiconductor configurations. Reservoirs may be finite or infinite volume.

Donehoo, L. K.↗

NASA NTRS1977

An automated data generator for NASTRAN

A modeling system based on construction-in-context of geometry and physical properties data is adapted for NASTRAN finite element data generation. This system was originally developed for generating three-dimensional finite data for composite structures in PATCHES-III using parametric cubic models. These finite line, surface and volume models constructed using simple data directives will be subdivided into NASTRAN finite elements and their supporting data in the new data generator. Any geometry can be constructed as illustrated by a compressor fan blade with camber, twist and changing airfoil section and the shape then subdivided using a uniform or non-uniform mesh. Special emphasis is given to the adaption of construction-in-context to composite material property modeling. A PMAT series of directives is described for any one, two or three-dimensionally reinforced composite that allows the phase properties, as well as composite properties, to be generated.

Stanton, E. L.↗

NASA NTRS1989

A parametric heat transfer study for cryogenic ball bearings in SSME HPOTP

A numerical modeling is to examine the effects of coolant convective heat transfer coefficient and frictional heating on the local temperature characteristics of a ball element in Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP) bearing. The present modeling uses a control-volume based, finite-difference method to solve the non-dimensionalized heat conduction equation in spherical coordinate system. The dimensionless temperature is found as a function of Biot number, heat flux ratio between the two race contacts, and location in the ball. The current results show that, for a given cooling capability, the ball temperature generally increases almost linearly with the heat input from the race-contacts. This increase is always very high at one of the two contacts. An increase in heat transfer coefficient generally reduces the ball temperature and alleviates the temperature gradient, except for the regions very close to the race contacts. For a 10-fold increase of heat transfer coefficient, temperature decrease is 35 percent for the average over entire ball, and 10 percent at the inner-race contact. The corresponding change of temperature gradient displays opposing trends between the regions immediately adjacent to the contacts and the remaining portion of the ball. The average temperature gradient in the vicinity of both contacts increases approximately 70 to 100 percent. A higher temperature gradient produces excessive thermal stress locally which may be detrimental to the material integrity. This, however, is the only unfavorable issue for an increase of heat transfer coefficient.

Chyu, Mingking K.↗

NASA NTRS1990

Numerical simulation of solidification in a horizontal cylindrical annulus charged with an aqueous salt solution

A computational study of solidification of a binary Na2CO3 solution in a horizontal cylindrical annulus is performed using a continuum formulation with a control-volume based, finite-difference scheme. The initial conditions were selected to facilitate the study of counter thermal and solutal convection, accompanied by extensive mushy region growth. Numerical results are compared with experimental data with mixed success. Qualitative agreement is obtained for the overall solidification process and associated physical phenomena. However, the plume thickness calculated for the solutally-driven convective upflow is substantially smaller than the observed value. Evolution of double-diffusive layers is predicted, but over a time scale much smaller than that observed experimentally. Good agreement is obtained between predicted and measured results for solid growth, but the mushy region thickness is significantly overpredicted.

Neilson, D. G.↗

NASA NTRS1991

Polymers And Riblets Reduce Hydrodynamic Skin Friction

Polymers injected into riblet grooves dramatically reduce polymer flow rate required for drag reduction. Polymer solution injected into valleys of grooves through array of holes or slots angled downstream to keep injected streams within grooves. Injection repeated some distance downstream because volumes of grooves finite and polymer becomes depleted as slowly pulled from groove by turbulence. Potentially useful for oil tankers as means of markedly reducing cost of fuel and used extensively on submarines, other ships, and other marine vehicles.

Bushnell, Dennis M.↗