Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 685 records · Page 38

Computational simulation of intermingled-fiber hybrid composite behavior

Three-dimensional finite-element analysis and a micromechanics based computer code ICAN (Integrated Composite Analyzer) are used to predict the composite properties and microstresses of a unidirectional graphite/epoxy primary composite with varying percentages of S-glass fibers used as hydridizing fibers at a total fiber volume of 0.54. The three-dimensional finite-element model used in the analyses consists of a group of nine fibers, all unidirectional, in a three-by-three unit cell array. There is generally good agreement between the composite properties and microstresses obtained from both methods. The results indicate that the finite-element methods and the micromechanics equations embedded in the ICAN computer code can be used to obtain the properties of intermingled fiber hybrid composites needed for the analysis/design of hybrid composite structures. However, the finite-element model should be big enough to be able to simulate the conditions assumed in the micromechanics equations.

Mital, Subodh K.↗

User's manual for the one-dimensional hypersonic experimental aero-thermodynamic (1DHEAT) data reduction code

A FORTRAN computer code for the reduction and analysis of experimental heat transfer data has been developed. This code can be utilized to determine heat transfer rates from surface temperature measurements made using either thin-film resistance gages or coaxial surface thermocouples. Both an analytical and a numerical finite-volume heat transfer model are implemented in this code. The analytical solution is based on a one-dimensional, semi-infinite wall thickness model with the approximation of constant substrate thermal properties, which is empirically corrected for the effects of variable thermal properties. The finite-volume solution is based on a one-dimensional, implicit discretization. The finite-volume model directly incorporates the effects of variable substrate thermal properties and does not require the semi-finite wall thickness approximation used in the analytical model. This model also includes the option of a multiple-layer substrate. Fast, accurate results can be obtained using either method. This code has been used to reduce several sets of aerodynamic heating data, of which samples are included in this report.

Hollis, Brian R.↗

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗

Informatics and High Throughput Screening of Thermophysical Properties

The combination of computer-aided experiments with computational modeling enables a new class of powerful tools for materials research. A non-contact method for measuring density, thermal expansion, and creep of undercooled and high-temperature materials has been developed, using electrostatic levitation and optical diagnostics, including digital video. These experiments were designed to take advantage of the large volume of data (many gigabytes/experiment, terabytes/campaign) to gain additional information about the samples. For example, using sub-pixel interpolation to measure about 1000 vectors per image of the sample's surface allows the density of an axisymmetric sample to be determined to an accuracy of about 200 ppm (0.02%). A similar analysis applied to the surface shape of a rapidly rotating sample is combined with finite element modeling to determine the stress-dependence of creep in the sample in a single test. Details of the methods for both the computer-aided experiments and computational models will be discussed.

Hyers, Robert W.↗

Thermal Pollution Mathematical Model. Volume 3: User's Manual for One-Dimensional Numerical Model for the Seasonal Thermocline

A user's manual for a one dimensional thermal model to predict the temperature profiles of a deep body of water for any number of annual cycles is presented. The model is essentially a set of partial differential equations which are solved by finite difference methods using a high speed digital computer. The model features the effects of area change with depth, nonlinear interaction of wind generated turbulence and buoyancy, adsorption of radiative heat flux below the surface, thermal discharges, and the effects of vertical convection caused by discharge. The main assumption in the formulation is horizontal homogeneity. The environmental impact of thermal discharges from power plants is emphasized. Although the model is applicable to most lakes, a specific site (Lake Keowee, S.C.) application is described in detail. The programs are written in FORTRAN 5.

Lee, S. S.↗

Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations

One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. These algorithms are now used in the practical engineering design of devices such as cars and airplanes as well as theoretical studies of complex phenomena such as fluid turbulence. In past years, limited computer resources have motivated the development of efficient numerical methods in computational fluid dynamics (CFD) utilizing structured meshes. These meshes are comprised of systematic arrays of quadrilateral or hexahedral cells. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Structured meshes also permit the use of highly efficient solution techniques such as alternating direction implicit (ADI) iteration schemes or multigrid. Following the dramatic improvement in computing speed in recent years, emphasis has shifted towards the design of algorithms capable of treating complex geometries. The automatic generation of structured grids about complex geometries is problematic. Unstructured grids offer one promising alternative technique for treating these general geometries. Unstructured meshes have irregular connectivity and usually contain triangles and/or quadrilaterals in two dimensions and tetrahedra and/or hexahedra in three dimensions. The generation and use of unstructured grids poses new challenges in computational fluid dynamics. This is true for both grid generation as well as for the design of algorithms for flow solution. The purpose of these notes is to present recent developments in the unstructured grid generation and flow solution technology.

Barth, T. J.↗

Zonal techniques for flowfield simulation about aircraft

A technique for performing conservative flowfield calculations on zonal meshes is described. The underlying flow solver is an implicit, upwind finite volume scheme which can incorporate either a perfect gas or an equilibrium air equation of state. Two different approaches which yield identical results, in terms of performing a conservative flux calculation on a zonal interface, are described and compared in terms of numerical efficiency. The capability of the method to handle relatively complex geometries is demonstrated by considering the flowfield about a model SR71 aircraft.

Walters, Robert W.↗

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗

Euler solutions for high-speed flow about complex three-dimensional configurations

A numerical algorithm based on a finite-volume explicit scheme with Runge-Kutta time integration of the Euler equations is presented for calculating high-speed three-dimensional flow about complex aerospace configurations. The use of enhancing factors such as artificial dissipative terms, enthalpy damping and local time-stepping are described. An algebraic method for generating quasi-three-dimensional computational grids for realistic aerospace configurations is presented. Computed results for various three-dimensional bodies at different Mach numbers and angles of attack have been obtained using the methods for grid-generation and flow simulation. Comparison of computed and experimental data for an advanced tactical aircraft-like configuration is presented, and a reasonable agreement of the data is noticed.

Moitra, A.↗

The efficient calculation of chemically reacting flow

A semi-implicit finite volume formulation is used to study flows with chemical reactions. In this formulation the source terms resulting from the chemical reactions are treated implicitly and the resulting system of partial differential equations is solved using two time-stepping schemes. The first is based on the Runge-Kutta method while the second is based on an Adams predictor-corrector method. Results show that improvements in computational efficiency depend to a large extent on the manner in which the source term is treated. Further, analysis and computation indicate that the Runge-Kutta method is more efficient than the Adams methods. Finally, an adaptive time stepping scheme is developed to study problems involving shock ignition. Calculations for a hydrogen-air system agree well with other methods.

Eklund, D. R.↗

Numerical solution of the Euler equation for compressible inviscid fluids

An effort is made to develop a satisfactory numerical method for the calculation of steady solutions of the Euler equations for inviscid compressible gas flows. The intended application is the prediction of the aerodynamic properties of aircraft flying at transonic speeds. Particular consideration is given to the modification of the equations to improve convergence to a steady state; finite volume formulation; adaptive dissipation; schemes designed to improve the resolution of shock waves; hybrid multistage time stepping schemes; residual averaging; a multigrid scheme; and trials with Burgers' equation. Results clearly demonstrate that the convergence of a time-dependent hyperbolic system to a steady state can be substantially accelerated by the introduction of multiple grids.

Jameson, A.↗

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

Chapter 4: Physically informed deep learning networks for simulating microstructure evolution of 3D polycrystals

As discussed in the previous chapter, high energy diffraction microscopy (HEDM) is used to study the micromechanical evolution of a material during in situ loading. HEDM experiments have been used to verify crystal plasticity (CP) simulations [119, 91, 90, 120], for experimental planning, material design, and to further analyze experimental results. However, Fast Fourier transform-based CP (CP-FFT) or finite element-based CP (CP-FE) methods are often too slow to be used in real-time during an experiment. CP-FFT is faster than CP-FE simulations due to the absence of meshing, but can still take hours to simulate the response of a single volume depending on the size and number of strain steps [127]. Reducing computation time would create a larger exploration space in planning and design, and enable faster analysis of experimental results and real-time feedback during an experiment. This research expands upon previous works to develop a workflow for predicting the full-field evolution of a 3D polycrystal. The workflow is simplified from previous works to predict only orientation and elastic strain tensors (from which stress tensors are calculated). The network is physically informed through loss functions and network architecture for a more robust model. The orientation predictions are informed about the cubic crystal symmetry of the material by incorporating disorientation and misorientation information into the network architecture and loss. The Von Mises stress is used to enforce the correct stress-strain trends in the strain tensor predictions. Additional total strain steps from the elastic and elastoplastic region are included to better capture the stress-strain evolution at smaller total strain steps. Material and hardening parameters are additional inputs into the networks to further inform the network and to study the network’s ability to predict different materials other than those used for training.

36 MATERIALS SCIENCE↗

Subplane Decusping for BWRs in MPACT

Control blade cusping can introduce significant error in boiling water reactor (BWR) calculations with MPACT when blade tips fall partway within an axial method of characteristics (MOC) plane, requiring homogenization of controlled and uncontrolled regions. This work implements subplane decusping for BWRs in MPACT by enabling BWR-compatible subplane coarse mesh finite difference (CMFD) and extending the existing decusping framework to represent between-assembly control blades that insert from the bottom of the core. The method resolves axial heterogeneity on a refined subplane mesh in the low-order solve and uses the resulting subplane flux shape to form flux-volume homogenized transport cross sections for the partially rodded MOC plane. The capability is evaluated using a single physics General Electric (GE)-14 assembly and a multiphysics Peach Bottom Unit 2 Type 1 assembly (PB2T1A) with thermal hydraulic feedback. In both cases, coarse axial meshes with and without subplane decusping are compared against fine mesh reference solutions over the full range of blade withdrawal positions. Subplane decusping reduces maximum/average eigenvalue errors from 3,369/220 pcm to 303/36 pcm for GE-14 and from 6,689/585 pcm to 1,582/109 pcm for PB2T1A. Additionally, it reduces pin power root mean square errors from 5.2%/0.8% to 1.8%/0.2% for GE-14 and from 6.9%/0.6% to 2.3%/0.1% for PB2T1A. These results demonstrate an effective, practical correction for BWR blade cusping in VERA-MPACT.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical boundary condition procedures and multigrid methods; Proceedings of the Symposium, NASA Ames Research Center, Moffett Field, CA, October 19-22, 1981

Papers presented in this volume provide an overview of recent work on numerical boundary condition procedures and multigrid methods. The topics discussed include implicit boundary conditions for the solution of the parabolized Navier-Stokes equations for supersonic flows; far field boundary conditions for compressible flows; and influence of boundary approximations and conditions on finite-difference solutions. Papers are also presented on fully implicit shock tracking and on the stability of two-dimensional hyperbolic initial boundary value problems for explicit and implicit schemes.

Source record↗

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

Improved Boundary Conditions for Cell-centered Difference Schemes

Cell-centered finite-volume (CCFV) schemes have certain attractive properties for the solution of the equations governing compressible fluid flow. Among others, they provide a natural vehicle for specifying flux conditions at the boundaries of the physical domain. Unfortunately, they lead to slow convergence for numerical programs utilizing them. In this report a method for investigating and improving the convergence of CCFV schemes is presented, which focuses on the effect of the numerical boundary conditions. The key to the method is the computation of the spectral radius of the iteration matrix of the entire demoralized system of equations, not just of the interior point scheme or the boundary conditions.

VanderWijngaart, Rob F.↗