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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Rapid finite-difference computation of subsonic and transonic aerodynamic flows

Rapid iterative (or semidirect) computation methods are developed for the finite-difference solution of the nonlinear equations of subsonic and transonic aerodynamics. At each iteration, a fast, direct elliptic algorithm solves the entire computation field. In an application to subsonic flow over a lifting airfoil, the full nonlinear stream-function equation is solved. Finally, a direct Cauchy-Riemann solver is used for the nonlinear transonic small-disturbance equations for a biconvex airfoil. At M = 0.7, t/c = 0.1 (subcritical), three iterations on a 39 x 32 mesh (totaling 2.45 sec on an IBM 360/67 computer) obtain convergence within 0.1%. A slightly supercritical case requires seven iterations (6.75 sec) for convergence within 1%.

Martin, E. D.↗
NASA NTRS2015

Estimation of Surface Temperature and Heat Flux by Inverse Heat Transfer Methods Using Internal Temperatures Measured While Radiantly Heating a Carbon/Carbon Specimen up to 1920 F

The ability to solve the heat conduction equation is needed when designing materials to be used on vehicles exposed to extremely high temperatures; e.g. vehicles used for atmospheric entry or hypersonic flight. When using test and flight data, computational methods such as finite difference schemes may be used to solve for both the direct heat conduction problem, i.e., solving between internal temperature measurements, and the inverse heat conduction problem, i.e., using the direct solution to march forward in space to the surface of the material to estimate both surface temperature and heat flux. The completed research first discusses the methods used in developing a computational code to solve both the direct and inverse heat transfer problems using one dimensional, centered, implicit finite volume schemes and one dimensional, centered, explicit space marching techniques. The developed code assumed the boundary conditions to be specified time varying temperatures and also considered temperature dependent thermal properties. The completed research then discusses the results of analyzing temperature data measured while radiantly heating a carbon/carbon specimen up to 1920 F. The temperature was measured using thermocouple (TC) plugs (small carbon/carbon material specimens) with four embedded TC plugs inserted into the larger carbon/carbon specimen. The purpose of analyzing the test data was to estimate the surface heat flux and temperature values from the internal temperature measurements using direct and inverse heat transfer methods, thus aiding in the thermal and structural design and analysis of high temperature vehicles.

Pizzo, Michelle↗
NASA NTRS1986

A mixed pseudospectral/finite difference method for the axisymmetric flow in a heated, rotating spherical shell

For a Spacelab flight, a model experiment of the earth's atmospheric circulation has been proposed. This experiment is known as the Atmospheric General Circulation Experiment (AGCE). In the experiment concentric spheres will rotate as a solid body, while a dielectric fluid is confined in a portion of the gap between the spheres. A zero gravity environment will be required in the context of the simulation of the gravitational body force on the atmosphere. The present study is concerned with the development of pseudospectral/finite difference (PS/FD) model and its subsequent application to physical cases relevant to the AGCE. The model is based on a hybrid scheme involving a pseudospectral latitudinal formulation, and finite difference radial and time discretization. The advantages of the use of the hybrid PS/FD method compared to a pure second-order accurate finite difference (FD) method are discussed, taking into account the higher accuracy and efficiency of the PS/FD method.

Macaraeg, M. G.↗
NASA NTRS1978

Comparative study of the convergence rates of two numerical techniques

The paper examines the applicability of the three-step Stetter (1968) method to the problem of hypersonic viscous flow over a blunt axisymmetric body used for planetary entry probes at zero angle of attack. The flow-field results using the two-step finite-difference MacCormack (1969) method are reported by Kumar and Graves (1977). Only the computational efficiency of Stetter's method is compared with that of MacCormack's in terms of the iterative time steps and computing time required for the steady-state solution. Advantages of Stetter's method over MacCormack's are established.

Kumar, A.↗
DOE OSTI2022

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Computing Edge States without Hard Truncation

We present a numerical method which accurately computes the discrete spectrum and associated bound states of semi-infinite Hamiltonians which model electronic “edge” states localized at boundaries of one- and two-dimensional crystalline materials. The problem is nontrivial since arbitrarily large finite “hard” (Dirichlet) truncations of the Hamiltonian in the infinite bulk direction tend to produce spurious bound states partially supported at the truncation. Our method, which overcomes this difficulty, is to compute the Green's function of the semi-infinite Hamiltonian by imposing an appropriate boundary condition in the bulk direction; then, the spectral data is recovered via Riesz projection. We demonstrate our method's effectiveness by studies of edge states at a graphene zig-zag edge in the presence of defects modeled both by a discrete tight-binding model and a continuum PDE model under finite difference discretization. Our method may also be used to study states localized at domain wall-type edges in one- and two-dimensional materials where the edge Hamiltonian is infinite in both directions; we demonstrate this for the case of a tight-binding model of distinct honeycomb structures joined along a zig-zag edge. Here, we expect our method to be useful for designing novel devices based on precise wave-guiding by edge states.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Cell averaging Chebyshev methods for hyperbolic problems

A cell averaging method for the Chebyshev approximations of first order hyperbolic equations in conservation form is described. Formulas are presented for transforming between pointwise data at the collocation points and cell averaged quantities, and vice-versa. This step, trivial for the finite difference and Fourier methods, is nontrivial for the global polynomials used in spectral methods. The cell averaging methods presented are proven stable for linear scalar hyperbolic equations and present numerical simulations of shock-density wave interaction using the new cell averaging Chebyshev methods.

Wei, Cai↗
NASA NTRS1991

A comparison of numerical methods for the Rayleigh equation in unbounded domains

A second-order finite difference and two spectral methods, including a Chebyshev tau and a Chebyshev collocation method were implemented to determine the linear hydrodynamic stability of an unbounded shear flow. The velocity profile of the basic flow in the stability analysis mimicks that of a one-stream free mixing layer. Local and global eigenvalue solution methods are used to determine individual eigenvalues and the eigenvalue spectrum, respectively. The calculated eigenvalue spectrum includes a discrete mode, a continuous spectrum associated with the equation singularity and a continuous spectrum associated with the domain unboundedness. The efficiency and the accuracy of these discretization methods in the prediction of the eigensolutions of the discrete mode were evaluated by comparison with a conventional shooting procedure. Their capabilities in mapping out the continuous eigenvalue spectra are also discussed.

Liou, W. W.↗
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

A numerical study of the unsteady leading edge separation bubble on an oscillating airfoil

A preliminary study of the unsteady viscous flow in the region of an airfoil leading edge was performed, in which the interaction between the viscous and inviscid flow fields is neglected. The solution method uses the finite difference form of the governing equations throughout the separated flow field and incorporates a transition model based on the integral turbulence kinetic energy equation. The validity of the numerical procedure is verified by making comparisons with analytical solutions to several test problems, including unsteady flow over a plate oscillating in its own plane. The method was then applied to the problem of unsteady viscous flow over a NACA 0012 airfoil oscillating sinusoidally in pitch. The flow field characteristics were in qualitative agreement with experimental results. The bubble moved forward on the airfoil and decreased in size as incidence was increased. Viscous flow in the leading edge region was found to be quasi-steady, while bubble height varied inversely with Reynolds number.

Kreskovsky, J. P.↗
DOE OSTI2019

A design procedure for fan-out improvement in all-optical photonic crystal logic design

In this paper, a general method is proposed to improve the fan-out parameter for all-optical logic gates and functions. Two different types of devices are designed to increase the fan-out of a logic operation and its inverted logic. Nonlinear cavities are used to design driving units. Silicon nanocrystal is used as the nonlinear material to create the required frequency shift for different values of input power. Plane wave expansion and finite difference time domain methods are used to simulate and analyze the proposed structures. The propagation delay of the proposed structure is less than 1.5 ps and the maximum required power for fan-out of three is 15 W.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1991

Application of the boundary element method to transient heat conduction

An advanced boundary element method (BEM) is presented for the transient heat conduction analysis of engineering components. The numerical implementation necessarily includes higher-order conforming elements, self-adaptive integration and a multiregion capability. Planar, three-dimensional and axisymmetric analyses are all addressed with a consistent time-domain convolution approach, which completely eliminates the need for volume discretization for most practical analyses. The resulting general purpose algorithm establishes BEM as an attractive alternative to the more familiar finite difference and finite element methods for this class of problems. Several detailed numerical examples are included to emphasize the accuracy, stability and generality of the present BEM. Furthermore, a new efficient treatment is introduced for bodies with embedded holes. This development provides a powerful analytical tool for transient solutions of components, such as casting moulds and turbine blades, which are cumbersome to model when employing the conventional domain-based methods.

Dargush, G. F.↗
NASA NTRS1985

Numerical solution of the vertical structure equation in the normal mode method

In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.

Sasaki, Y. K.↗
NASA NTRS1981

Viscous-inviscid interaction for transonic wing-body configurations including wake effects

An existing three-dimensional compressible integral boundary-layer method was modified to account for mean dilatation effects, to model transition properly, and to provide better numerical stability near computational boundaries. Results of this method were compared against those from a three-dimensional finite-difference boundary-layer method on a difficult test case. An interaction procedure was developed to couple this integral method with a number of wing-alone and wing-body transonic potential codes to account for viscous effects. A strip wake model, including thickness and curvature effects, was developed and incorporated into this interaction procedure. Results from this procedure were compared against experimental data and results from previous procedures, on test cases where viscous effects were large.

Streett, C. L.↗
DOE OSTI2019

Nanoscale plasmonic detector of wave intensity difference and uni-directional waveguide

In this paper, a new nanoscale plasmonic differential detector of wave intensities is proposed. This structure can be useful for calibration applications and for finding out if a wave intensity is high enough to pass a threshold level or not. Next, a nanoscale plasmonic uni-directional waveguide is presented. The distinction of this waveguide with similar cases is that the one-way wave-guiding is provided by adding a properly located controlling source, without employing any non-reciprocal material. This feature can be used in many applications that only uni-directional wave propagation is desirable. Both structures have been designed in metal–insulator–metal structure including cavity resonators. The operations are based on the coupling of the wave from input waveguides to the cavity resonators and decoupling to output waveguides in the cavity resonance wavelengths. The simulation results of performance of the structures, as elements applicable to photonic integrated circuits, have been obtained using the numerical method of the finite difference time domain.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2017

Large scale Brownian dynamics of confined suspensions of rigid particles

We introduce methods for large-scale Brownian Dynamics (BD) simulation of many rigid particles of arbitrary shape suspended in a fluctuating fluid. Our method adds Brownian motion to the rigid multiblob method [F. Balboa Usabiaga et al., Commun. Appl. Math. Comput. Sci. 11(2), 217–296 (2016)] at a cost comparable to the cost of deterministic simulations. We demonstrate that we can efficiently generate deterministic and random displacements for many particles using preconditioned Krylov iterative methods, if kernel methods to efficiently compute the action of the Rotne-Prager-Yamakawa (RPY) mobility matrix and its “square” root are available for the given boundary conditions. These kernel operations can be computed with near linear scaling for periodic domains using the positively split Ewald method. Here we study particles partially confined by gravity above a no-slip bottom wall using a graphical processing unit implementation of the mobility matrix-vector product, combined with a preconditioned Lanczos iteration for generating Brownian displacements. We address a major challenge in large-scale BD simulations, capturing the stochastic drift term that arises because of the configuration-dependent mobility. Unlike the widely used Fixman midpoint scheme, our methods utilize random finite differences and do not require the solution of resistance problems or the computation of the action of the inverse square root of the RPY mobility matrix. We construct two temporal schemes which are viable for large-scale simulations, an Euler-Maruyama traction scheme and a trapezoidal slip scheme, which minimize the number of mobility problems to be solved per time step while capturing the required stochastic drift terms. We validate and compare these schemes numerically by modeling suspensions of boomerang-shaped particles sedimented near a bottom wall. Using the trapezoidal scheme, we investigate the steady-state active motion in dense suspensions of confined microrollers, whose height above the wall is set by a combination of thermal noise and active flows. We find the existence of two populations of active particles, slower ones closer to the bottom and faster ones above them, and demonstrate that our method provides quantitative accuracy even with relatively coarse resolutions of the particle geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Enriched immersed finite element and isogeometric analysis: algorithms and data structures

Immersed finite element methods provide a convenient analysis framework for problems involving geometrically complex domains, such as those found in topology optimization and microstructures for engineered materials. However, their implementation remains a major challenge due to, among other things, the need to apply nontrivial stabilization schemes and generate custom quadrature rules. This article introduces the robust and computationally efficient algorithms and data structures comprising an immersed finite element preprocessing framework. The input to the preprocessor consists of a background mesh and one or more geometries defined on its domain. The output is structured into groups of elements with custom quadrature rules formatted such that common finite element assembly routines may be used without or with only minimal modifications. The key to the preprocessing framework is the construction of material topology information, concurrently with the generation of a quadrature rule, which is then used to perform enrichment and generate stabilization rules. While the algorithmic framework applies to a wide range of immersed finite element methods using different types of meshes, integration, and stabilization schemes, the preprocessor is presented within the context of the extended isogeometric analysis. This method utilizes a structured B-spline mesh, a generalized Heaviside enrichment strategy considering the material layout within individual basis functions’ supports, and face-oriented ghost stabilization. Using a set of examples, the effectiveness of the enrichment and stabilization strategies is demonstrated alongside the preprocessor’s robustness in geometric edge cases. Additionally, the performance and parallel scalability of the implementation are evaluated.

Computer implementation↗