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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Wing-Body Aeroelasticity Using Finite-Difference Fluid/Finite-Element Structural Equations on Parallel Computers

This paper presents a procedure for computing the aeroelasticity of wing-body configurations on multiple-instruction, multiple-data (MIMD) parallel computers. In this procedure, fluids are modeled using Euler equations discretized by a finite difference method, and structures are modeled using finite element equations. The procedure is designed in such a way that each discipline can be developed and maintained independently by using a domain decomposition approach. A parallel integration scheme is used to compute aeroelastic responses by solving the coupled fluid and structural equations concurrently while keeping modularity of each discipline. The present procedure is validated by computing the aeroelastic response of a wing and comparing with experiment. Aeroelastic computations are illustrated for a High Speed Civil Transport type wing-body configuration.

Byun, Chansup↗

NASA NTRS1997

Continued Development of the Discontinuous Galerkin Method for Computational Aeroacoustic Applications

The formulation and the implementation of boundary conditions within the context of the quadrature-free form of the discontinuous Galerkin method are presented for several types of boundary conditions for the Euler equations. An important feature of the discontinuous Galerkin method is that the interior point algorithm is well behaved in the neighborhood of the boundary and requires no modifications. This feature leads to a simple and accurate treatment for wall boundary conditions and simple inflow and outflow boundary conditions. Curved walls are accurately treated with only minor changes to the implementation described in earlier work. The 'perfectly matched layer' approach to nonreflecting boundary conditions is easily applied to the discontinuous Galerkin. The compactness of the discontinuous Galerkin method makes it better suited for buffer-zone-type methods than high-order finite-difference methods. Results are presented for wall, characteristic inflow and outflow, and nonreflecting boundary conditions.

Atkins, H. L.↗

NASA NTRS1996

Wing-Body Aeroelasticity on Parallel Computers

This article presents a procedure for computing the aeroelasticity of wing-body configurations on multiple-instruction, multiple-data parallel computers. In this procedure, fluids are modeled using Euler equations discretized by a finite difference method, and structures are modeled using finite element equations. The procedure is designed in such a way that each discipline can be developed and maintained independently by using a domain decomposition approach. A parallel integration scheme is used to compute aeroelastic responses by solving the coupled fluid and structural equations concurrently while keeping modularity of each discipline. The present procedure is validated by computing the aeroelastic response of a wing and comparing with experiment. Aeroelastic computations are illustrated for a high speed civil transport type wing-body configuration.

Guruswamy, Guru P.↗

NASA NTRS2017

A Non-Intrusive Algorithm for Sensitivity Analysis of Chaotic Flow Simulations

We demonstrate a novel algorithm for computing the sensitivity of statistics in chaotic flow simulations to parameter perturbations. The algorithm is non-intrusive but requires exposing an interface. Based on the principle of shadowing in dynamical systems, this algorithm is designed to reduce the effect of the sampling error in computing sensitivity of statistics in chaotic simulations. We compare the effectiveness of this method to that of the conventional finite difference method.

Blonigan, Patrick J.↗

DOE OSTI2017

Orbital Advection with Magnetohydrodynamics and Vector Potential

Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear terms, and we find that we do not need an advective gauge since solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited to non-magnetic problems. Benchmarked results of (hydrodynamical) planet–disk interaction and of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high-frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficient for straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1978

Computational advances in fluid dynamics

An office size digital computer with capabilities exceeding a CDC 7600 can be constructed for less than 300 K and spectral methods offer a factor ten improvement over finite difference methods. Integration of these facts can produce two orders of magnitude reduction in cost of fluid dynamic computations and permit solution of 3-D problems which were not previously feasible.

Taylor, T. D.↗

NASA NTRS1981

Pseudospectral method for problems with nonperiodic boundary conditions

The capability of the pseudospectral technique and the time integration scheme for the time dependent nonperiodic problem is demonstrated. Significant improvements are shown in the accuracy for calculation of the two dimensional time dependent advection diffusion equation. The pseudospectral method is a viable alternative to the finite difference methods for the equation of the spatial deriviations.

Lee, H. N.↗

DOE OSTI2017

Numerical Modeling of Ultra Wideband Combined Antennas

With the help of a program we developed, based on the finite difference method in the time domain, we have investigated the characteristics of ultra wideband combined antennas in detail. The antennas were developed to radiate bipolar pulses with durations in the range 0.5–3 ns. Data obtained by numerical modeling are compared with the data of experimental studies on antennas and have been used in the synthesis of electromagnetic pulses with maximum field strength.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

The Interface Influence in TiN/SiN{sub x} Multilayer Nanocomposite Under Irradiation

The paper focuses on studying the kinetics of radiation-induced point defects formed in TiN/SiN{sub x} multilayer nanocomposites with account of their generation, diffusion recombination, and the influence of sinks functioning as interfaces. In order to describe the kinetics in nanocrystalline TiN and amorphous SiN{sub x} phases, a finite-difference method is used to solve the system of balance kinetic equations for absolute defect concentrations depending on the spatiotemporal variables. A model of the disclination-dislocation interface structure is used to study the absorption of radiation-induced point defects on the boundaries in created stress fields. It is shown that the interface effectively absorbs point defects in these phases of TiN/SiN{sub x} multilayer nanocomposite, thereby reducing their amount within the space between phases. This behavior of point defects partially explains a mechanism of the radiation resistance in this type of nanocomposites.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2019

Miniband formation engineering in GaN/AlN superlattices with constant total effective length

In this work, we study the miniband formation procedure in GaN/AlN constant total effective length superlattices by means of the subband energy calculations. We calculate the transmission coefficients and miniband structures of the systems by means of the transfer matrix and finite difference methods, respectively. The miniband structures obtained by using these methods confirm each other. Here, we observe a nonlinear miniband gap behavior as the number of wells changes. Now, we can tune the width of each miniband and minigap by means of the number of wells and the total effective length of the system. By using the strategy of fixing the total length of the system and optimizing the number of layers inside it, we can find the same miniband widths as the usual superlattice (fixing the well width, not the total system length). However, in the earlier case (our proposal), fabrication of the device may need fewer amounts of the material and may lead to smaller device sizes, because we fixed the total system length.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model

We present a semi-numerical solution to a fractal telegraphic dual-porosity fluid flow model. It combines Laplace transform and finite difference schemes. The Laplace transform handles the time variable whereas the finite difference method deals with the spatial coordinate. This semi-numerical scheme is not restricted by space discretization and allows the computation of a solution at any time without compromising numerical stability or the mass conservation principle. Our formulation results in a non-analytically-solvable second-order differential equation whose spatial discretization yields a tridiagonal linear algebraic system. Moreover, we describe comparisons between semi-numerical and semi-analytical solutions for particular cases. Results agree well with those from semi-analytic solutions. Furthermore, we expose a parametric analysis from the coupled model in order to show the effects of relevant parameters on pressure profiles and flow rates for the case where neither analytic nor semi-analytic solutions are available.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for high Reynolds number laminar flows

A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.

Kim, Sang-Wook↗

NASA NTRS1987

Development of computational methods for unsteady aerodynamics at the NASA Langley Research Center

The current scope, recent progress, and plans for research and development of computational methods for unsteady aerodynamics at the NASA Langley Research Center are reviewed. Both integral equations and finite difference methods for inviscid and viscous flows are discussed. Although the great bulk of the effort has focused on finite difference solution of the transonic small perturbation equation, the integral equation program is given primary emphasis here because it is less well known.

Yates, E. Carson, Jr.↗

DOE OSTI2023

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗

NASA NTRS1984

Improved finite-difference vibration analysis of pretwisted, tapered beams

An improved finite difference procedure based upon second order central differences is developed. Several difficulties encountered in earlier works with fictitious stations that arise in using second order central differences, are eliminated by developing certain recursive relations. The need for forward or backward differences at the beam boundaries or other similar procedures is eliminated in the present theory. By using this improved theory, the vibration characteristics of pretwisted and tapered blades are calculated. Results of the second order theory are compared with published theoretical and experimental results and are found to be in good agreement. The present method generally produces close lower bound solutions and shows fast convergence. Thus, extrapolation procedures that are customary with first order finite-difference methods are unnecessary. Furthermore, the computational time and effort needed for this improved method are almost the same as required for the conventional first order finite-difference approach.

Subrahmanyam, K. B.↗

DOE OSTI2018

Numerical investigation on two-fluid model (micropolar-Newtonian) for pulsatile flow of blood in a tapered arterial stenosis with radially variable magnetic field and core fluid viscosity

In the present study, an unsteady two-fluid model of blood through a tapered arterial stenosis with variable viscosity in the presence of variable magnetic field has been investigated. In this model, blood in the core region is assumed to be micropolar and plasma in the peripheral layer as Newtonian. Finite difference method is employed in solving the governing equations. The solutions for velocity, flow rate, wall shear stress and flow resistance are computed numerically. A comparison between the velocity profiles obtained by the present study and the experimental data is made and a good agreement between them is found. The model is used to study the effect of parameters such as taper angle, radially variable viscosity, hematocrit, the coupling number, the micropolar parameter, magnetic field and plasma layer thickness on physiologically important parameters such as wall shear stress distribution in the stenotic region and flow resistance. The results in the case of constant magnetic field and variable magnetic field are compared to study the effects of magnetic field on the flow of blood. It is found that the magnitudes of wall shear stress and flow resistance are higher in the case of variable magnetic field. It is important to note that the flow resistance is higher for magneto-micropolar fluid than the micropolar fluid. The wall shear stress decreases with increasing hematocrit parameter whereas flow resistance increases with hematocrit when the applied pressure gradient is held fixed. The model clearly shows the situation of a patient with a tapered arterial stenosis under feverish condition. Due care has been taken to compare the present theoretical results with the existing ones including experimental results and good agreement between them has been observed both qualitatively and quantitatively.

42 ENGINEERING↗

NASA NTRS2009

Ablative Thermal Response Analysis Using the Finite Element Method

A review of the classic techniques used to solve ablative thermal response problems is presented. The advantages and disadvantages of both the finite element and finite difference methods are described. As a first step in developing a three dimensional finite element based ablative thermal response capability, a one dimensional computer tool has been developed. The finite element method is used to discretize the governing differential equations and Galerkin's method of weighted residuals is used to derive the element equations. A code to code comparison between the current 1-D tool and the 1-D Fully Implicit Ablation and Thermal Response Program (FIAT) has been performed.

Dec John A.↗

NASA NTRS1993

A total variation diminishing finite difference algorithm for sonic boom propagation models

It is difficult to accurately model the rise phases of sonic boom waveforms with traditional finite difference algorithms because of finite difference phase dispersion. This paper introduces the concept of a total variation diminishing (TVD) finite difference method as a tool for accurately modeling the rise phases of sonic booms. A standard second order finite difference algorithm and its TVD modified counterpart are both applied to the one-way propagation of a square pulse. The TVD method clearly outperforms the non-TVD method, showing great potential as a new computational tool in the analysis of sonic boom propagation.

Sparrow, Victor W.↗