Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-Difference Methods”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

An implicit finite-difference solution to the viscous shock layer, including the effects of radiation and strong blowing

An implicit finite-difference scheme is developed for the fully coupled solution of the viscous, radiating stagnation-streamline equations, including strong blowing. Solutions are presented for both air injection and injection of carbon-phenolic ablation products into air at conditions near the peak radiative heating point in an earth entry trajectory from interplanetary return missions. A detailed radiative-transport code that accounts for the important radiative exchange processes for gaseous mixtures in local thermodynamic and chemical equilibrium is utilized in the study. With minimum number of assumptions for the initially unknown parameters and profile distributions, convergent solutions to the full stagnation-line equations are rapidly obtained by a method of successive approximations. Damping of selected profiles is required to aid convergence of the solutions for massive blowing. It is shown that certain finite-difference approximations to the governing differential equations stabilize and improve the solutions. Detailed comparisons are made with the numerical results of previous investigations. Results of the present study indicate lower radiative heat fluxes at the wall for carbonphenolic ablation than previously predicted.

Garrett, L. B.↗
NASA NTRS1974

Fast direct numerical solution of the nonhomogeneous Cauchy-Riemann equations

A fast direct (noniterative) 'Cauchy-Riemann Solver' is developed for solving the finite-difference equations representing systems of first-order elliptic partial differential equations in the form of the nonhomogeneous Cauchy-Riemann equations. The method is second-order accurate and requires approximately the same computer time as a fast cyclic-reduction Poisson solver. The accuracy and efficiency of the direct solver are demonstrated in an application to solving an example problem in aerodynamics: subsonic inviscid flow over a biconvex airfoil. The analytical small-perturbation solution contains singularities, which are captured well by the computational technique. The algorithm is expected to be useful in nonlinear subsonic and transonic aerodynamics.

Lomax, H.↗
NASA NTRS1974

Development of solution techniques for nonlinear structural analysis

Nonlinear structural solution methods in the current research literature are classified according to order of the solution scheme, and it is shown that the analytical tools for these methods are uniformly derivable by perturbation techniques. A new perturbation formulation is developed for treating an arbitrary nonlinear material, in terms of a finite-difference generated stress-strain expansion. Nonlinear geometric effects are included in an explicit manner by appropriate definition of an applicable strain tensor. A new finite-element pilot computer program PANES (Program for Analysis of Nonlinear Equilibrium and Stability) is presented for treatment of problems involving material and geometric nonlinearities, as well as certain forms on nonconservative loading.

Vos, R. G.↗
NASA NTRS1956

Extrapolation techniques applied to matrix methods in neutron diffusion problems

A general matrix method is developed for the solution of characteristic-value problems of the type arising in many physical applications. The scheme employed is essentially that of Gauss and Seidel with appropriate modifications needed to make it applicable to characteristic-value problems. An iterative procedure produces a sequence of estimates to the answer; and extrapolation techniques, based upon previous behavior of iterants, are utilized in speeding convergence. Theoretically sound limits are placed on the magnitude of the extrapolation that may be tolerated. This matrix method is applied to the problem of finding criticality and neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron-diffusion equations is treated. Results for this example are indicated.

Mccready, Robert R↗
NASA NTRS1973

The transonic aerofoil problem with embedded shocks.

The integral equation approach to the mixed flow problem of infinite wings at high subsonic speeds is adopted for non-circulatory and circulatory (lifting) flows. The solutions are determined from a system of non-linear algebraic equations and, to ensure always unique solutions, the method of differentiation with respect to a parameter has been applied. The resulting Cauchy problem is then solved with the linearised flow solution as the initial value vector. For the case of embedded shocks in the flow field, the method of steepest descent has been added to the calculation scheme. Results for subcritical and supercritical flows past aerofoils are given and compared with solutions obtained by finite-difference techniques.

Norstrud, H.↗
NASA NTRS1973

Preliminary study of a wall stabilized constructed arc

An iterative, implicit, finite-difference numerical technique is described which is suitable for obtaining solutions to the governing equations for a gas flowing in an axially symmetric constricted-arc heater. The method is shown to provide adequate solutions for three cases of simple pipe flows found in the literature, and for flow in a constricted-arc heater by direct comparison with experimental data. The comparison with arc-heater data includes static pressure, arc voltage, and wall heat flux, all as a function of axial location, and a radial temperature profile at one axial station. The arc-heater data were taken with air as the test gas at a heater inlet pressure approximately 0.40 atm and two flow rates of 2.2 and 4.8 g/sec. The arc currents investigated were between 377 and 584 amperes.

Graves, R. A.↗
NASA NTRS1972

Thermochemical nonequilibrium in atomic hydrogen at elevated temperatures

A numerical study of the nonequilibrium flow of atomic hydrogen in a cascade arc was performed to obtain insight into the physics of the hydrogen cascade arc. A rigorous mathematical model of the flow problem was formulated, incorporating the important nonequilibrium transport phenomena and atomic processes which occur in atomic hydrogen. Realistic boundary conditions, including consideration of the wall electrostatic sheath phenomenon, were included in the model. The governing equations of the asymptotic region of the cascade arc were obtained by writing conservation of mass and energy equations for the electron subgas, an energy conservation equation for heavy particles and an equation of state. Finite-difference operators for variable grid spacing were applied to the governing equations and the resulting system of strongly coupled, stiff equations were solved numerically by the Newton-Raphson method.

Scott, R. K.↗
NASA NTRS1974

Numerical calculations of two dimensional, unsteady transonic flows with circulation

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data.

Beam, R. M.↗
DOE OSTI2025

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1950

Direct method of design and stress analysis of rotating disks with temperature gradient

A method is presented for the determination of the contour of disks, typified by those of aircraft gas turbines, to incorporate arbitrary elastic-stress distributions resulting from either centrifugal or combined centrifugal and thermal effects. The specified stress may be radial, tangential, or any combination of the two. Use is made of the finite-difference approach in solving the stress equations, the amount of computation necessary in the evolution of a design being greatly reduced by the judicious selection of point stations by the aid of a design chart. Use of the charts and of a preselected schedule of point stations is also applied to the direct problem of finding the elastic and plastic stress distribution in disks of a given design, thereby effecting a great reduction in the amount of calculation. Illustrative examples are presented to show computational procedures in the determination of a new design and in analyzing an existing design for elastic stress and for stresses resulting from plastic flow.

Manson, S S↗
NASA NTRS1973

Nonlinear transient responses of structures by the spatial finite-element method.

Based upon the Principle of Virtual Work and D'Alembert's Principle, the assumed-displacement version of the spatial finite-element method is developed to predict the large deflection transient responses of structures including elastic-plastic, strain hardening, and strain-rate material behavior. The formulations are developed in detail for curved beamlike structures undergoing planar (1) Bernouilli-Euler-type or (2) Timoshenko-type deformation behavior. The resulting equations of motion are solved timewise by a finite-difference numerical procedure. The present predictions are evaluated via several beam and ring examples for which experimental measurements and independent finite-difference predictions in both space and time are available; very good agreement is noted. The consequences of employing several types of timewise finite-difference operators are examined. Also, some comparisons between finite-element predictions and finite-difference predictions are shown to illustrate 'typical comparisons' of efficiency for a given prediction accuracy.

Wu, R. W.-H.↗
DOE OSTI2025

HFBTHO-AD: Differentiation of a nuclear energy density functional code

The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.

Algorithmic differentiation↗
NASA NTRS1972

Numerical solutions of reactive fluid flows during postignition transients in hybrid rocket systems.

A computational method has been developed for the study of the post-ignition transients in hybrid rocket systems. The particular system chosen consisted of a gaseous oxidizer flowing within a tube of solid fuel, resulting in heterogeneous combustion. With the appropriate assumptions, two-dimensional, time-dependent conservation equations were derived for the reacting gas phase, and for the solid phase, in a cylindrical coordinate system. These were then programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and the Crank-Nicolson scheme for the solid phase. Appropriate initial and boundary conditions were represented, including heat and mass conservation at the interface between gas and solid. Initially, no attempt was made to relate the recession rate at the surface to the surface temperature, or to include heat transfer by radiation. A simple case was selected for preliminary calculations, with aluminum and oxygen as fuel and oxidizer, and aluminum oxide as the product.

Hung, W. S. Y.↗
DOE OSTI2026

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Formation and optical properties of indium nanoparticle arrays for deep-UV plasmonics

We utilize a combined computational-experimental approach to examine the influence of indium nanoparticle (NP) array distributions on deep-ultraviolet (UV) plasmon resonances. For photon energies < 5.7 eV, analysis of ellipsometric spectra reveals an increase in silicon reflectance induced by indium NP arrays on silicon. For various energies in the range 5.7–7.0 eV, a decrease in reflectance is induced by the NP arrays. Similar trends in reflectance are predicted from finite-difference time-domain (FDTD) simulations using NP size distributions extracted from atomic-force micrographs as input. In addition, in the energy range of 7.4–9.2 eV, the FDTD simulations reveal reflectance minima, characteristic of localized surface plasmon resonances. Here, electron energy-loss spectroscopy collected from individual indium NPs reveals the presence of LSPR at ≈ 8 eV, further supporting the promise of indium NP arrays on silicon for deep-UV plasmonics.

36 MATERIALS SCIENCE↗