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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Viscous flow solutions with a cubic spline approximation

A cubic spline approximation is used for the solution of several problems in fluid mechanics. This procedure provides a high degree of accuracy even with a nonuniform mesh, and leads to a more accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several typical integration schemes are presented. For two-dimensional flows a spline-alternating-direction-implicit (SADI) method is evaluated. The spline procedure is assessed and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.

Rubin, S. G.↗
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 NTRS1975

Subgrid or Reynolds stress-modeling for three-dimensional turbulence computations

A review is given of recent advances in two distinct computational methods for evaluating turbulence fields, namely, statistical Reynolds stress modeling and turbulence simulation, where large eddies are followed in time. It is shown that evaluation of the mean Reynolds stresses, rather than use of a scalar eddy viscosity, permits an explanation of streamline curvature effects found in several experiments. Turbulence simulation, with a new volume averaging technique and third-order accurate finite-difference computing is shown to predict the decay of isotropic turbulence in incompressible flow with rather modest computer storage requirements, even at Reynolds numbers of aerodynamic interest.

Rubesin, M. W.↗
NASA NTRS1976

Improved interactive calculation procedure for supersonic flows

An interactive numerical procedure has been developed for supersonic viscous flows (either two-dimensional or axisymmetric configurations). The flow field is divided into two regions: (1) an inner region which is highly viscous and mostly subsonic, and (2) an outer region where the flow is supersonic and in which viscous effects are small, but not negligible. This paper presents a detailed description of: I. Outer Region - numerical solution obtained by applying the method of characteristics to a system of equations which includes viscous and conduction transport terms only normal to the streamlines; II. Inner Region - treated by a system of equations of the boundary-layer type that includes higher order effects, such as longitudinal and transverse curvature and normal pressure gradients (equations are coupled and solved simultaneously in physical coordinates, using an implicit finite-difference scheme); III. Interactive Procedure - in the interaction mode, the two regions are coupled iteratively along a matching line, where the Mach number is of the order of 1.2.

Tassa, Y.↗
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.↗
NASA NTRS1976

A numerical study of viscous flow around an airfoil

An integrodifferential method, previously formulated in terms of velocity and vorticity vectors, is reformulated in terms of stream function and vorticity for two-dimensional incompressible viscous flows. The reformulated integrodifferential method is shown to retain the distinguishing feature of the previous formulation in permitting the confinement of the solution field to the viscous region of the flow and consequently offers great computational advantages. The application of this procedure in a study of an incompressible flow around an impulsively started 9% thick symmetric Joukowski airfoil at an angle of attack of 15 deg and a Reynolds number of 1000 is discussed. Numerical results are presented and compared with available finite-difference results.

Wu, J. C.↗
NASA NTRS1976

A two-dimensional finite difference solution for the transient thermal behavior of tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique was analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗
NASA NTRS1977

A perturbation solution for interlaminar stresses in bidirectional laminates

The paper considers a thin elastic bidirectional symmetric composite laminate subjected to a uniform axial strain. All interlaminar stresses are determined by a zeroth-order perturbation analysis of the governing equations. The method of solution considers the laminate to be composed of an interior region (removed from the free edge) and a boundary layer region. The interior region is analyzed by dropping terms related to powers of the thickness-to-width ratio in the dimensionless elasticity equations. The boundary layer region is studied by introducing a stretching transformation to the governing equations. A uniform solution is then formed by satisfying Prandtl's matching principle of perturbation theory. Results for a four-layer bidirectional graphite-epoxy laminate are compared with existing finite-difference results. It is shown that the perturbation analysis provides more insight into the stress distribution near the free edge.

Hsu, P. W.↗
NASA NTRS1974

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗
NASA NTRS1974

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗
NASA NTRS1975

Time-dependent transonic flow solutions for axial turbomachinery

Three-dimensional unsteady transonic flow through an axial turbomachine stage is described in terms of a pair of two-dimensional formulations pertaining to orthogonal surfaces, namely, a blade-to-blade surface and a hub-to-casing surface. The resulting systems of nonlinear, inviscid, compressible equations of motion are solved by an explicit finite-difference technique. The blade-to-blade program includes the periodic interaction between rotor and stator blade rows. Treatment of the boundary conditions and of the blade slipstream motion by a characteristic type procedure is discussed in detail. Harmonic analysis of the acoustic far field produced by the blade row interaction, including an arbitrary initial transient, is outlined. Results from the blade-to-blade program are compared with experimental measurements of the rotating pressure field at the tip of a high-speed fan. The hub-to-casing program determines circumferentially averaged flow properties on a meridional plane. Blade row interactions are neglected in this formulation, but the force distributions over the entire blade surface for both the rotor and stator are obtained. Results from the hub-to-casing program are compared with a relaxation method solution for a subsonic rotor. Results are also presented for a quiet fan stage which includes transonic flow in both the rotor and stator and a normal shock in the stator.

Erdos, J.↗
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↗
NASA NTRS1972

Finite element analysis of large transient elastic-plastic deformations of simple structures, with application to the engine rotor fragment containment/deflection problem

Assumed-displacement versions of the finite-element method are developed to predict large-deformation elastic-plastic transient deformations of structures. Both the conventional and a new improved finite-element variational formulation are derived. These formulations are then developed in detail for straight-beam and curved-beam elements undergoing (1) Bernoulli-Euler-Kirchhoff or (2) Timoshenko deformation behavior, in one plane. For each of these categories, several types of assumed-displacement finite elements are developed, and transient response predictions are compared with available exact solutions for small-deflection, linear-elastic transient responses. The present finite-element predictions for large-deflection elastic-plastic transient responses are evaluated via several beam and ring examples for which experimental measurements of transient strains and large transient deformations and independent finite-difference predictions are available.

Wu, R. W.↗
NASA NTRS1974

Liquefaction dynamics of n-octadecane in cylindrical coordinates

The study presents one method of mathematically modeling the effects of gravity-induced free convection upon the melting phenomena of a finite paraffin slab. The theoretical model is based on the numerical computer solution of the solid phase and liquid phase energy equations, coupled with an imposed ideal-viscous flow approximation. An implicit alternating-direction technique is used as the finite-difference approximation for the basic energy equations in the numerical solution. The theoretical predictions are compared to experimental data from an annular test cell, in which n-octadecane, C18H38, was used as the test material. The results show that the model is a good first approximation to the solution of solid-liquid phase change involving gravity-induced free convection.

Bain, R. L.↗