Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference 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 1,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A fast marching algorithm for the factored eikonal equation

The eikonal equation is instrumental in many applications in several fields ranging from computer vision to geoscience. This equation can be efficiently solved using the iterative Fast Sweeping (FS) methods and the direct Fast Marching (FM) methods. However, when used for a point source, the original eikonal equation is known to yield inaccurate numerical solutions, because of a singularity at the source. In this case, the factored eikonal equation is often preferred, and is known to yield a more accurate numerical solution. One application that requires the solution of the eikonal equation for point sources is travel time tomography. This inverse problem may be formulated using the eikonal equation as a forward problem. While this problem has been solved using FS in the past, the more recent choice for applying it involves FM methods because of the efficiency in which sensitivities can be obtained using them. However, while several FS methods are available for solving the factored equation, the FM method is available only for the original eikonal equation. In this paper we develop a Fast Marching algorithm for the factored eikonal equation, using both first and second order finite-difference schemes. Our algorithm follows the same lines as the original FM algorithm and requires the same computational effort. In addition, we show how to obtain sensitivities using this FM method and apply travel time tomography, formulated as an inverse factored eikonal equation. Numerical results in two and three dimensions show that our algorithm solves the factored eikonal equation efficiently, and demonstrate the achieved accuracy for computing the travel time. We also demonstrate a recovery of a 2D and 3D heterogeneous medium by travel time tomography using the eikonal equation for forward modeling and inversion by Gauss–Newton.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

FEHM V3.0 - V3.X

The numerical background of the FEHM computer code can be traced to the early 1970s when it was used to simulate geothermal and hot dry rock reservoirs. The primary use over a number of years was to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the potential Yucca Mountain repository. Today FEHM is used to simulate groundwater and contaminant flow and transport in deep and shallow, fractured and un-fractured porous media throughout the US DOE complex. FEHM has proved to be a valuable asset on a variety of projects of national interest including Environmental Remediation of the Nevada Test Site, the LANL Groundwater Protection Program, geologic CO2 sequestration, Enhanced Geothermal Energy (EGS) programs, Oil and Gas production, Nuclear Waste Isolation, and Arctic Permafrost. Subsurface physics has ranged from single fluid/single phase fluid flow when simulating basin scale groundwater aquifers to complex multifluid/ multi-phase fluid flow that includes phase change with boiling and condensing in applications such as unsaturated zone surrounding nuclear waste storage facility or leakage of CO2/brine through faults or wellbores. The numerical method used in FEHM is the control volume method (CV) for fluid flow and heat transfer equations which allows FEHM to exactly enforce energy/mass conservation; while an option is available to use the finite element (FE) method for displacement equations to obtain more accurate stress calculations. In addition to these standard methods, an option to use FE for flow is available, as well as a simple Finite Difference scheme.

Pawar, Rajesh↗
DOE OSTI2019

Enhancements to the Subplane Capability in MPACT

In recent years, much of the development effort for the Virtual Environment for Reactor Applications (VERA) has been devoted to improving performance and targeting industry-class computers. This includes reducing both runtime (improving efficiency) and core-count requirements to enable simulation on smaller machines. In MPACT, in particular, many aspects of the neutronic iteration strategy have been evaluated and improved, such as the subgroup and method of characteristics (MOC) runtimes, and many efforts have focused on improving the coarse mesh finite difference (CMFD) runtime, which has become a dominant component of the runtime as other aspects have been improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Irradiation of Advanced Cladding Specimens in the High Flux Isotope Reactor: Capsule Designs and Test Matrix

The Advanced Fuels Campaign (AFC) has initiated the Advanced Reactor Cladding (ARC) irradiation campaign to generate irradiation performance data for candidate fuel cladding concepts. The campaign includes a diverse set of ferritic/martensitic steels, oxide dispersion strengthened (ODS) alloys, FeCrAlbased alloys, coated materials, and welded cladding specimens produced through multiple US Department of Energy (DOE) programs and international collaborations. Three complementary experimental thrusts comprise the campaign: tensile testing (ARC Tensile) to rapidly screen candidate alloys, fracture toughness testing (ARC Fracture) to evaluate irradiation effects on crack resistance, and tubular weld testing (ARC Weld) to quantify irradiation-induced changes in the mechanical performance of end cap welds. This report documents the irradiation campaign design, including the selected materials, specimen types, irradiation matrix, and capsule designs for irradiation within the High Flux Isotope Reactor (HFIR). A total of 14 irradiation capsules were developed to achieve target irradiation temperatures between 300°C and 600°C and doses up to 30 dpa. Thermal analyses were performed using finite element methods to establish capsule geometries capable of achieving the desired specimen temperatures while accommodating differences in specimen geometry and material properties. The resulting capsule designs provide the basis for irradiation of the AFC-ARC experimental matrix and subsequent post-irradiation examination to assess the effects of neutron irradiation on advanced cladding materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1976

Multidimensional difference schemes with fourth-order accuracy

An explicit finite-difference algorithm is presented for the solution of quasilinear divergence free multidimensional hyperbolic systems. The method consists of four steps per time level. The resulting scheme is fourth-order accurate in both space and time, though the intermediate steps are only first-order accurate. The family of schemes introduced is dissipative, and hence, suitable for both smooth flows and flows containing shocks. This method is compared, in several numerical examples, with both second-order schemes and others that are fourth order in space, but second order in time.

Turkel, E.↗
NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗
NASA NTRS1978

Computational transonic flow about realistic aircraft configurations

A numerical method has been developed to compute transonic flows about realistic wing-fuselage configurations. The finite difference scheme employs an improved small disturbance flow equation. A unique grid embedding technique, which was heretofore applied to airfoils and wings, has been extended to include the treatment of both body and wing-body shapes. The resulting high-density mesh is shown to be a valuable asset in resolving details of the three-dimensional flow. A mathematical modeling system is used to process arbitrary fuselage geometries for body boundary conditions. Correlations with experimental data for simple isolated bodies, an isolated fuselage, and wing-fuselage combinations are included.

Boppe, C. W.↗
NASA NTRS1983

Numerical simulation of an electrothermal deicer pad

A numerical simulation is developed to investigate the removal of ice from composite aircraft blades by means of electrothermal deicing. The model considers one dimensional, unsteady state heat transfer in the composite blade-ice body. The heat conduction equations are approximated by using the Crank-Nicolson finite difference scheme, and the phase change in the ice layer is handled using the Enthalpy method. To solve the system of equations which result, Gauss-Seidel iteration is used. The simulation computes the temperature profile in the composite blade-ice body, as well as the movement of the ice-water interface, as a function of time. This information can be used to evaluate deicer performance. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Marano, J. J.↗
NASA NTRS1983

Computation of potential flows with embedded vortex rings and applications to helicopter rotor wakes

A finite difference scheme for solving the motion of a number of vortex rings is developed. The method is an adaptation of the 'cloud-in-cell' technique to axisymmetric flows, and is thus a combined Eulerian-Lagrangian technique. A straightforward adaptation of the cloud-in-cell scheme to an axisymmetric flow field is shown to introduce a grid dependent self-induced velocity to each vortex ring. To correct this behavior the potential is considered to consist of two parts, a local and a global field. An improved difference formula is derived, allowing the accurate calculation of the potential at points near vortex locations. The local potential is then subtracted before calculating the velocity, leaving only the influences of the remaining vortices. The correct self-induced velocity is then explicitly added to the vortex velocity. Calculations of the motion of one and two vortex rings are performed, demonstrating the ability of the new method to eliminate the grid dependence of the self-induced velocity. The application of the method to the calculation of helicopter rotor flows in hover is attempted.

Roberts, T. W.↗
NASA NTRS1982

Constitutive model development for isotropic materials

The objective is to develop a unified constitutive model for finite-element structural analysis of turbine engine hot section components. This effort constitutes a different approach for nonlinear finite-element computer codes which were heretofore based on classical inelastic methods. A unified constitutive theory will avoid the simplifying assumptions of classical theory and should more accurately represent the behavior of superalloy materials under cyclic loading conditions and high temperature environments. Model development will be directed toward isotropic, cast nickel-base alloys used for aircooled turbine blades and vanes. The contractor will select a base material for model development and an alternate material for verification purposes from a list of three alloys specified by NASA. The candidate alloys represent a cross-section of turbine blade and vane materials of interest to both large and small size engine manufacturers. Material stock for the base and alternate materials will be supplied to the Contractor by the government.

Kaufman, A.↗
NASA NTRS1988

FASTWO - A 2-D interactive algebraic grid generator

This paper presents a very simple and effective computational procedure, FASTWO, for generating patched composite finite difference grids in 2-D for any geometry. Major components of the interactive graphics based method that is closely akin to and borrows many tools from transfinite interpolation are highlighted. Several grids produced by FASTWO are shown to illustrate its powerful capability. Comments about extending the methodology to 3-D are also given.

Luh, Raymond Ching-Chung↗
NASA NTRS1989

Adaptive methods, rolling contact, and nonclassical friction laws

Results and methods on three different areas of contemporary research are outlined. These include adaptive methods, the rolling contact problem for finite deformation of a hyperelastic or viscoelastic cylinder, and non-classical friction laws for modeling dynamic friction phenomena.

Oden, J. T.↗
NASA NTRS1992

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗
NASA NTRS1993

High Temperature Composite Analyzer (HITCAN) theoretical manual, version 1.0

This manual outlines some of the theoretical aspects embedded in the computer code, HITCAN (High Temperature Composite ANalyzer). HITCAN is a general purpose computer program for predicting nonlinear global structural and local stress-strain response of arbitrarily oriented, multilayered high temperature metal matrix composite structures. This code combines composite mechanics and laminate theory with an internal data base of the constituents (matrix, fiber and interface) material properties. The thermal and mechanical properties of the constituents are considered to be nonlinearly dependent on several parameters including temperature, stress and stress rate. The computational procedure for the analysis of the composite structure uses the finite element method. HITCAN consists of three modules: COBSTRAN, METCAN and MHOST. COBSTRAN generates the geometry (pre-processor) and defines the layup of the different plies. METCAN computes the material behavior of the composite and of the constituents. Finally, MHOST is a finite element program based on the mixed iterative solution technique. MHOST has a library for 2D and 3D isoparametric elements. HITCAN is written in FORTRAN 77 computer language and has been configured and executed on the NASA Lewis Research Center CRAY XMP and YMP computers.

Lackney, J. J.↗
NASA NTRS1995

Aeroelasticity of wing and wing-body configurations on parallel computers

The objective of this research is to develop computationally efficient methods for solving aeroelasticity problems on parallel computers. Both uncoupled and coupled methods are studied in this research. For the uncoupled approach, the conventional U-g method is used to determine the flutter boundary. The generalized aerodynamic forces required are obtained by the pulse transfer-function analysis method. For the coupled approach, the fluid-structure interaction is obtained by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗
NASA NTRS1994

Low-Density Nozzle Flow by the Direct Simulation Monte Carlo and Continuum Methods

Two different approaches, the direct simulation Monte Carlo (DSMC) method based on molecular gasdynamics, and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gasdynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental data, in which density and rotational temperature variations along the centerline and at various locations inside a low-density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data, both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity-temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chang-Hong↗