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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

L{sup ∞} error bound of conservative compact difference scheme for the generalized symmetric regularized long-wave (GSRLW) equations

In this paper, we design a compact finite difference scheme which preserves the original conservative properties to solve the generalized symmetric regularized long-wave equations. The existence of the difference solution is proved by the Brouwer fixed-point theorem. Applying the discrete energy method, the convergence and stability of the difference scheme is obtained, and its numerical convergence order is O(τ{sup 2}+h{sup 4}) in the L{sup ∞}-norm for u and ρ. For computing the nonlinear algebraic system generated by the compact scheme, a decoupled iterative algorithm is constructed and proved to be convergent. Numerical experiment results show that the theory is accurate and the method is efficient and reliable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Hexagonal Geometries in MPACT

The MPACT code is a high-fidelity light-water reactor analysis code using whole-core pin-resolved neutron transport calculations on modern parallel-computing hardware. MPACT uses the 2D/1D method to solve 3D neutron transport problems by decomposing the problem into a stack of 2D slices, each of which is solved independently using the method of characteristics (MOC). The slices are then coupled axially using the P3 nodal expansion method (NEM-P3) for the 1D axial calculations. MPACT also employs the coarse mesh finite difference (CMFD) method to accelerate calculations. This manuscript details work supporting advanced reactor designs using hexagonal pins and hexagonal assemblies such as the VVER-1000. If performed correctly, MOC is geometry agnostic. However, MPACT previously had optimizations in place for Cartesian geometries, specifically in the modularization and current calculations. Sections 2 and 3 detail the changes made to MPACT to support MOC and CMFD calculations on hexagonal geometries. Section 4 reports results demonstrating solution consistency for problems run with and without CMFD acceleration, results demonstrating solution consistency when run in serial and parallel, and pincell results using the Monte Carlo code, McCard’s benchmark results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2018

Lebedev Scheme for Ultrasound Simulation in Composites

The growing use of composite materials for aerospace applications has resulted in a need for quantitative nondestructive evaluation (NDE) methods appropriate for characterizing damage in composite components. NDE simulation tools, such as ultrasound models, can aid in enabling optimized inspection methods and establishing confidence in inspection capabilities. In this paper, a mathematical approach using the Lebedev Finite Difference (LFD) method is presented for ultrasonic wave simulation in composites. Boundary condition equations for implementing stress-free boundaries (necessary for simulation of NDE scenarios) are also presented. Quantitative comparisons between LFD guided wave ultrasound simulation results, experimental guided wave data, and dispersion curves are described. Additionally, stability tests are performed to establish the LFD code behavior in the presence of stress-free boundaries and low-symmetry anisotropy. Results show that LFD is an appropriate approach for simulating ultrasound in anisotropic composite materials and that the method is stable in the presence of low-symmetry anisotropy and stress-free boundaries. Studies presented in this paper include guided wave simulation in hexagonal, monoclinic, triclinic and layered composite laminates.

Guided wave↗
NASA NTRS1989

A review of high-speed, convective, heat-transfer computation methods

Formulations useful for engineering formulations are sought, with a view to a clearer physical understanding, of the phenomena that govern convective aerodynamic heating at the elevated speeds encountered in atmospheric missile trajectories. After discussing cold-wall laminar boundary layer heating, the complex boundary layer transition phenomenon and cold-wall turbulent boundary layer heating are treated. The current understanding of separated flow-region and shock-interaction heating is presented, together with an evaluation of heat-protection methods and a characterization of the influence of mass addition on laminar and turbulent boundary layers. Finite-difference method-based CFD code results are evaluated.

Tauber, Michael E.↗
DOE OSTI2021

Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems

Multipole moments are the first-order responses of the energy to spatial derivatives of the electric field strength. The quality of density functional theory prediction of molecular multipole moments thus characterizes errors in modeling the electron density itself, as well as the performance in describing molecules interacting with external electric fields. However, only the lowest non-zero moment is translationally invariant, making the higher-order moments origin-dependent. Therefore, instead of using the 3 × 3 quadrupole moment matrix, we utilize the translationally invariant 3 × 3 matrix of second cumulants (or spatial variances) of the electron density as the quantity of interest (denoted by K). The principal components of K are the square of the spatial extent of the electron density along each axis. A benchmark dataset of the principal components of K for 100 small molecules at the coupled cluster singles and doubles with perturbative triples at the complete basis set limit is developed, resulting in 213 independent K components. The performance of 47 popular and recent density functionals is assessed against this Var213 dataset. Several functionals, especially double hybrids, and also SCAN and SCAN0 predict reliable second cumulants, although some modern, empirically parameterized functionals yield more disappointing performance. The H, Li, and Be atoms, in particular, are challenging for nearly all methods, indicating that future functional development could benefit from the inclusion of their density information in training or testing protocols.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1987

The finite-difference time-domain (FD-TD) method for electromagnetic scattering and interaction problems

The formulation and recent applications of the finite-difference time-domain (FD-TD) method for the numerical modeling of electromagnetic scattering and interaction problems are considered. It is shown that improvements in FD-TD modeling concepts and software implementation often make it a preferable choice for structures which cannot be easily treated by conventional integral equations and asymptotic approaches. Recent FD-TD modeling validations in research areas including coupling to wires and wire bundles in free space and cavities, scattering from surfaces in relativistic motion, inverse scattering, and radiation condition theory, are reviewed. Finally, the advantages and disadvantages of FD-TD, and guidelines concerning when FD-TD should and should not be used in high-frequency electromagnetic modeling problems, are summarized.

Taflove, A.↗
NASA NTRS1986

A full potential rotor analysis with wake influence using an inner-outer domain technique

A three-dimensional, quasi-steady, full potential flow solver was adapted to include realistic rotor wake influence for the aerodynamic analysis of helicopter rotors. The method uses an inner/outer domain technique to accommodate wake effects. Nonlinear flow is computed in the inner domain using a finite difference solution method. The wake is modeled using prescribed wake techniques to allow for the inclusion of realistic wake geometries. Portions of the wake passing inside the inner domain are treated using an embedded vortex technique. The procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. Correlation with measured lifting transonic data in hover and forward flight is shown which demonstrates the merits of the approach.

Egolf, T. A.↗
NASA NTRS1985

Selecting step sizes in sensitivity analysis by finite differences

This paper deals with methods for obtaining near-optimum step sizes for finite difference approximations to first derivatives with particular application to sensitivity analysis. A technique denoted the finite difference (FD) algorithm, previously described in the literature and applicable to one derivative at a time, is extended to the calculation of several simultaneously. Both the original and extended FD algorithms are applied to sensitivity analysis for a data-fitting problem in which derivatives of the coefficients of an interpolation polynomial are calculated with respect to uncertainties in the data. The methods are also applied to sensitivity analysis of the structural response of a finite-element-modeled swept wing. In a previous study, this sensitivity analysis of the swept wing required a time-consuming trial-and-error effort to obtain a suitable step size, but it proved to be a routine application for the extended FD algorithm herein.

Iott, J.↗
NASA NTRS1972

Large deflection elastic-plastic dynamic response of stiffened shells of revolution

The formulation and check out porblems for a computer code DYNAPLAS, which analyzes the large deflection elastic-plastic dynamic response of stiffened shells of revolution, are presented. The formulation for special discretization is by the finite element method with finite differences being used for the evaluation of the pseudo forces due to material and geometric nonlinearities. Time integration is by the Houbolt method. The stiffeners may be due to concentrated or distributed eccentric rings and spring supports at arbitrary angles around the circumference of the elements. Check out porblems include the comparison of solutions from DYNAPLAS with experimental and other computer solutions for rings, conical and cylindrical shells and a curved panel. A hypothetical submarine including stiffeners and missile tube is studied under a combination of hydrostatic and dynamically applied asymmetrical pressure loadings.

Stricklin, J. A.↗
NASA NTRS1974

Large deflection elastic-plastic dynamic response of stiffened shells of revolution

This paper presents the formulation and check-out problems for a computer code DYNAPLAS, which analyzes the large deflection elastic-plastic dynamic response of stiffened shells of revolution. The formulation for spacial discretization is by the finite element method with finite differences being used for the evaluation of the pseudo forces due to material and geometric nonlinearities. Time integration is by the Houbolt method or central differences. The stiffeners may be due to concentrated or distributed eccentric rings and spring supports at arbitrary angles around the circumference of the elements. Check-out problems include the comparison of solutions from DYNAPLAS with experimental and other computer solutions for rings and conical and cylindrical shells. A hypothetical submarine including stiffeners and missile tube is studied under a combination of hydrostatic and dynamically applied asymmetrical pressure loadings.

Stricklin, J. A.↗
DOE OSTI2021

Approximating and incorporating model uncertainty in an inversion for seismic source functions: Preliminary results

We present preliminary work on propagating model uncertainty into the estimation of the time domain source time functions of the seismic source. Our method is based on an estimated model covariance function, which we estimate from the data. The model covariance function is then used to construct a suite of surrogate Greens functions which we use in a Monte Carlo type inversion scheme. The result is a probability density function of the six independent source time functions, each of which corresponds to an individual component of the seismic moment tensor. We compare the results of our method with those obtained using a computationally expensive finite difference Monte Carlo method and find that our new method produces results that are deficient in low frequencies. The advantage of our new method, which we term the Karhunen-Loeve Monte Carlo (KLMC) method, is that is several orders of magnitude faster than our current method, which uses a finite difference scheme to produce the suite of forward models.

42 ENGINEERING↗
DOE OSTI2026

A sweeping positivity-preserving high-order finite difference WENO scheme for Euler equations

We develop a simple, high-order, conservative and robust positivity-preserving sweeping procedure for the density and the nonlinear pressure function in the compressible Euler equations. Using the scaling limiter in Zhang and Shu (J Comput Phys 229:3091–3120, 2010), we obtain a non-trivial extension of the scalar sweeping technique in Liu et al. (J Sci Comput 73:1028–1071, 2017) for the positivity of pressure. The sweeping procedure developed in this paper is a post-processing technique, which can be applied to any concave functions of the conserved variables in hyperbolic conservation law systems. Thus, it has applications beyond the Euler equations. This procedure preserves positivity and conservation of physical quantities without destroying the accuracy of the underlying scheme. The algorithm works for general schemes including finite difference, finite volume and discontinuous Galerkin methods; however, in this paper we focus on finite difference weighted essentially non-oscillatory (WENO) methods. As a result, we provide numerical tests of the fifth-order finite difference WENO scheme to demonstrate the accuracy and robustness of the technique.

Compressible Euler equations↗
NASA NTRS1974

Estimation of deviation angle for axial-flow compressor blade sections using inviscid-flow solutions

Development of a method of estimating deviation angles by analytical procedures was begun. Solutions for inviscid, irrotational flow in the blade-to-blade plane were obtained with a finite-difference calculation method. Deviation angles for a plane cascade with a rounded trailing edge were estimated by using the inviscid-flow solutions and three trailing-edge hypotheses. The estimated deviation angles were compared with existing experimental data over a range of incidence angles at inlet flow angles of 30 deg and 60 deg. The results indicate that deviation angles can be estimated accurately (within 1 deg) by using one of the three trailing-edge hypotheses, but only when pressure losses are low. A new trailing-edge hypotheses is presented which is suitable (for the cascade considered) for both low- and high-loss operating points.

Miller, M. J.↗
NASA NTRS1974

Deviation angle estimation for axial-flow compressors using inviscid flow solutions

Development of a more satisfactory method of estimating deviation angles was attempted using analytical procedures. Solutions for inviscid, irrotational flow in the blade-to-blade plane were obtained with a finite-difference calculation method. Deviation angles for three cascades with rounded trailing edges were estimated using the inviscid flow solutions and four existing trailing edge hypotheses. The estimated deviation angles were compared with experimental values for a range of incidence angle and inlet flow angle. None of the hypotheses were satisfactory at all incidence angles for all the cascades considered. An alternative hypothesis is presented which is suitable for one of the cascade profiles over a range of incidence and blade setting angle.

Miller, M. J.↗
NASA NTRS1993

Finite-Difference Algorithms For Computing Sound Waves

Governing equations considered as matrix system. Method variant of method described in "Scheme for Finite-Difference Computations of Waves" (ARC-12970). Present method begins with matrix-vector formulation of fundamental equations, involving first-order partial derivatives of primitive variables with respect to space and time. Particular matrix formulation places time and spatial coordinates on equal footing, so governing equations considered as matrix system and treated as unit. Spatial and temporal discretizations not treated separately as in other finite-difference methods, instead treated together by linking spatial-grid interval and time step via common scale factor related to speed of sound.

Davis, Sanford↗
NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗