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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

On the accuracy of the pseudocompressibility method in solving the incompressible Navier-Stokes equations

The method of pseudocompressibility is tested for its accuracy in solving the incompressible Navier-Stokes equations. An implicit, finite difference computer code is used to solve the equations in a three-dimensional, curvilinear coordinate system. The code employs artificial compressibiity for solving the pressure field, coupled with an implicit, approximate factorization scheme. This coupling is known as the pseudocompressibility method and introduces pressure waves of finite speed into the fluid medium that would otherwise have an infinite sound speed. The waves die out as the solution converges, and the steady state solution approaches a divergence-free condition. These waves, however, limit the time accuracy of the computations. The effects of these waves are analysed and criteria are set for choosing the pseudocompressibility parameter that governs the pressure wave speed; test cases are presented that verify these criteria. The code is tested by computing laminar flow over a two-dimensional, backward-facing step and over a two-dimensional, impulsively started circular cylinder.

Rogers, S. E.↗
NASA NTRS1988

Power flow in a beam using a 5-accelerometer probe

A method to determine structure-borne power flow is examined which uses central finite differences to approximate the field variables required to calculate the shear and bending components of flexural vibration. An experimental study using a beam driven by a shaker was performed with the aim of comparing estimates using this proposed five-accelerometer method and the conventional two-accelerometer, free-field method. The advantages and disadvantages of the new method are pointed out.

Downing, J. Micah↗
DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗
NASA NTRS1996

On DGS Relaxation: The Stokes Problem

Multigrid methods have proven to be efficient methods for solving partial differential equations (especially those of elliptic type). There is also growing experience with multigrid solvers for fluids problems, e.g., the Stokes and Navier-Stokes equations (using both finite element and finite difference discretizations). It is also well known that at the heart of any multigrid method is the smoother. In this work we look at a smoother introduced by Brandt and Dinar (DGS relaxation), and we examine some of its properties and consider some possible modifications to it. It is well known that multigrid performance using DGS relaxation is sensitive to the treatment of boundaries; this issue is addressed.

Meir, A. J.↗
NASA NTRS1976

Polynominal Interpolation Methods for Viscous Flow Calculations

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and by Hermitian (Taylor series) finite-difference discretization. The similarities and special features of these different developments are discussed. The governing systems apply for both uniform and variable meshes. Hybrid schemes resulting from two different polynomial approximations for the first and second derivatives lead to a nonuniform mesh extension of the so-called compact or Pad? difference technique (Hermite 4). A variety of fourth-order methods are described and the Hermitian approach is extended to sixth-order (Hermite 6). The appropriate spline boundary conditions are derived for all procedures. For central finite differences, this leads to a two-point, second-order accurate generalization of the commonly used three-point end-difference formula. Solutions with several spline and Hermite procedures are presented for the boundary layer equations, with and without mass transfer, and for the incompressible viscous flow in a driven cavity. Divergence and nondivergence equations are considered for the cavity. Among the fourth-order techniques, it is shown that spline 4 has the smallest truncation error. The spline 4 procedure generally requires one-quarter the number of mesh points in a given coordinate direction as a central finite-difference calculation of equal accuracy. The Hermite 6 procedure leads to remarkably accurate boundary layer solutions.

Rubin, S. G.↗
DOE OSTI2022

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Examination of Mean Stress Calculation Approaches in Rock Mechanics

The mean stress, as a fundamental statistical property of a group of stress data, is essential for stress variability characterisation. However, currently in rock mechanics, the mean stress is customarily and erroneously calculated by separately averaging the principal stress magnitudes and orientations. In order to draw the attention of our community to the appropriate approach for stress variability characterisation, here we compare the customary scalar/vector mean with that obtained by the mathematically rigorous tensorial approach—the Euclidean mean. Calculation of mean stress using both a small group of actual in situ stress measurement results and a large group of simulated stress data (obtained using the combined finite– discrete element method, FEMDEM) demonstrates that the two approaches yield different results. Further investigation of these results shows that the scalar/vector approach may yield non-unique and non-orthogonal mean principal stresses, and these may deviate significantly from those of the Euclidean mean. Our calculations and comparisons reveal that the scalar/ vector approach is deficient because it processes the principal stress magnitudes and orientations separately as independent quantities and ignores the connection between them. Conversely, the tensorial approach appropriately averages the tensors that simultaneously carry not only the information of stress magnitudes and orientations, but also the inherent relations between them. Therefore, arbitrarily using scalar/vector mean stress of in situ stress measurements as input in further rock engineering analyses may yield significantly erroneous results. Finally, we advise that stress data should be statistically processed in a tensorial manner using tensors referred to a common Cartesian coordinate system.

58 GEOSCIENCES↗
DOE OSTI2020

A matrix-free hyperviscosity formulation for high-order ALE hydrodynamics

The numerical approximation of compressible hydrodynamics is at the core of high-energy density (HED) multiphysics simulations as shocks are the driving force in experiments like inertial confinement fusion (ICF). In this work, we describe our extension of the hyperviscosity technique, originally developed for shock treatment in finite difference simulations, for use in arbitrarily high-order finite element methods for Lagrangian hydrodynamics. Hyperviscosity enables shock capturing while preserving the high-order properties of the underlying discretization away from the shock region. Specifically, we compute a high-order term based on a product of the mesh length scale to a high power scaled by a hyper-Laplacian operator applied to a scalar field. We then form the total artificial viscosity by taking a non-linear blend of this term and a traditional artificial viscosity term. We also present a matrix-free formulation for computing the finite element based hyper-Laplacian operator. Such matrix-free methods have superior performance characteristics compared to traditional full matrix assembly approaches and offer advantages for GPU based HPC hardware. We demonstrate the numerical convergence of our method and its application to complex, multi-material ALE simulations on high-order (curved) meshes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Measurements of skin friction on the wall of a hypersonic nozzle.

Direct skin friction measurements obtained on the wall of a Mach 19.8 nozzle are presented. These measurements of skin friction and values deduced from the profile data are compared with predictions from a finite difference theory and also with several 'flat-plate' prediction methods for skin friction.

Harvey, W. D.↗
NASA NTRS1972

Numerical calculations of flow fields

The solutions to the equations of motion for inviscid fluid flow around a pointed elliptic cone at incidence are presented. The numerical method used, MacCormack's second order preferential predictor-corrector finite difference approximation, is applied to the fluid flow equations derived in conservation-law form. The entropy boundary condition, hitherto unused for elliptic cone problems, is investigated and compared to reflection boundary condition solutions. The stagnation streamline movement of the inclined elliptic cone is noted and surface pressure coefficients are plotted. Also presented are solutions for an elliptic cone and a circular cone at zero incidence and a circular cone at a small angle of attack. Comparisons are made between these present solutions and previously published theory.

Anderson, D. M.↗
NASA NTRS1975

Finite element computation of a viscous compressible free shear flow governed by the time dependent Navier-Stokes equations

A finite element algorithm for solution of fluid flow problems characterized by the two-dimensional compressible Navier-Stokes equations was developed. The program is intended for viscous compressible high speed flow; hence, primitive variables are utilized. The physical solution was approximated by trial functions which at a fixed time are piecewise cubic on triangular elements. The Galerkin technique was employed to determine the finite-element model equations. A leapfrog time integration is used for marching asymptotically from initial to steady state, with iterated integrals evaluated by numerical quadratures. The nonsymmetric linear systems of equations governing time transition from step-to-step are solved using a rather economical block iterative triangular decomposition scheme. The concept was applied to the numerical computation of a free shear flow. Numerical results of the finite-element method are in excellent agreement with those obtained from a finite difference solution of the same problem.

Cooke, C. H.↗
NASA NTRS1977

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. Part 2: Delamination behavior and fracture mechanics parameters

Based on theories of laminate anisotropic elasticity and interlaminar fracture, the complete solution structure associated with a composite delamination is determined. Fracture mechanics parameters characterizing the interlaminar crack behavior are defined from asymptotic stress solutions for delaminations with different crack-tip deformation configurations. A numerical method employing singular finite elements is developed to study delaminations in fiber composites with any arbitrary combinations of lamination, material, geometric, and crack variables. The special finite elements include the exact delamination stress singularity in its formulation. The method is shown to be computationally accurate and efficient, and operationally simple. To illustrate the basic nature of composite delamination, solutions are shown for edge-delaminated (0/-0/-0/0) and (+ or - 0/+ or - 0/90/90 deg) graphite-epoxy systems under uniform axial extenstion. Three-dimensional crack-tip stress intensity factors, associated energy release rates, and delamination crack-closure are determined for each individual case. The basic mechanics and mechanisms of composite delamination are studied, and fundamental characteristics unique to recently proposed tests for interlaminar fracture toughness of fiber composite laminates are examined.

Wang, S. S.↗
NASA NTRS1983

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 2nd, California State University, Long Beach, CA, January 17-20, 1983, Proceedings

The present conference covers topics concerning the measurement and calculation of interactive flows, together with problems posed by subsonic and transonic wings, missiles, and ships. Discussions are presented on the time-dependent finite difference simulation of unsteady interactive flows, Navier-Stokes equation methods, numerical solutions for spatially periodic boundary layers, the application of unsteady laminar tripple deck theory to viscous-inviscid interaction, the coupling of boundary layer and Euler equation solutions, and viscous-inviscid flow interactions. Also discussed are leading and trailing edge flows, three-dimensional wing flows, small disturbance calculations including entropy corrections, an inviscid computational method for tactical missiles, and boundary layer and flow separation characteristics of bodies of revolution at incidence.

Source record↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. II - The delamination behavior and fracture mechanics parameters

Based on theories of laminate anisotropic elasticity and interlaminar fracture, the complete solution structure associated with a composite delamination is determined. Fracture mechanics parameters characterizing the interlaminar crack behavior are defined from asymptotic stress solutions for delaminations with different crack-tip deformation configurations. A numerical method employing singular finite elements is developed to study delaminations in fiber composites with any arbitrary combinations of lamination, material, geometric, and crack variables. The special finite elements include the exact delamination stress singularity in its formulation. The method is shown to be computationally accurate and efficient, and operationally simple. To illustrate the basic nature of composite delamination, solutions are shown for edge-delaminated (0/-0/-0/0) and (+ or - 0/+ or - 0/90/90 deg) graphite-epoxy systems under uniform axial extension. Three-dimensional crack-tip stress intensity factors, associated energy release rates, and delamination crack-closure are determined for each individual case. The basic mechanics and mechanisms of composite delamination are studied, and fundamental characteristics unique to recently proposed tests for interlaminar fracture toughness of fiber composite laminates are examined. Previously announced in STAR as N84-13222

Wang, S. S.↗
NASA NTRS1987

Lectures series in computational fluid dynamics

The lecture notes cover the basic principles of computational fluid dynamics (CFD). They are oriented more toward practical applications than theory, and are intended to serve as a unified source for basic material in the CFD field as well as an introduction to more specialized topics in artificial viscosity and boundary conditions. Each chapter in the test is associated with a videotaped lecture. The basic properties of conservation laws, wave equations, and shock waves are described. The duality of the conservation law and wave representations is investigated, and shock waves are examined in some detail. Finite difference techniques are introduced for the solution of wave equations and conservation laws. Stability analysis for finite difference approximations are presented. A consistent description of artificial viscosity methods are provided. Finally, the problem of nonreflecting boundary conditions are treated.

Thompson, Kevin W.↗
NASA NTRS1987

Wave envelope technique for multimode wave guide problems

A numerical solution technique for two-dimensional EM waveguide problems is developed analytically and demonstrated. The frequency-domain waveguide problem is transformed to a time-dependent problem and treated following the limiting-amplitude principle; the wave-envelope technique is applied; and the resulting problem is solved numerically using a finite-difference scheme based on a variation of the finite-area method (Baumeister, 1977 and 1984). Results for several sample problems, obtained in 10 min or less of execution time on PCs, are presented in tables.

Hariharan, S. I.↗
NASA NTRS1989

Finite element solvers for incompressible fluid flows and heat transfer

Two different finite-element solvers for incompressible viscous flow, i.e., the mixed interpolation method and the SIMPLE-type iterative method, are compared and tested with some benchmark problems. The advantages of the SIMPLE-type iterative method are the decoupling of the governing equations and the use of equal-order interpolation functions for both velocity and pressure. Even though there is a significant difference between the two methods in terms of the pressure field, similar solutions are obtained for the velocity field.

Sohn, Jeong L.↗