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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Analysis of Combustion Instability in Liquid Fuel Rocket Motors

The development of a technique to be used in the solution of nonlinear velocity-sensitive combustion instability problems is described. The orthogonal collocation method was investigated. It found that the results are heavily dependent on the location of the collocation points and characteristics of the equations, so the method was rejected as unreliabile. The Galerkin method, which has proved to be very successful in analysis of the pressure sensitive combustion instability was found to work very well. It was found that the pressure wave forms exhibit a strong second harmonic distortion and a variety of behaviors are possible depending on the nature of the combustion process and the parametric values involved. A one-dimensional model provides further insight into the problem by allowing a comparison of Galerkin solutions with more exact finite-difference computations.

Wong, K. W.↗
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 NTRS1978

Analysis of turbulent free jet hydrogen-air diffusion flames with finite chemical reaction rates

The nonequilibrium flow field resulting from the turbulent mixing and combustion of a supersonic axisymmetric hydrogen jet in a supersonic parallel coflowing air stream is analyzed. Effective turbulent transport properties are determined using the (K-epsilon) model. The finite-rate chemistry model considers eight reactions between six chemical species, H, O, H2O, OH, O2, and H2. The governing set of nonlinear partial differential equations is solved by an implicit finite-difference procedure. Radial distributions are obtained at two downstream locations of variables such as turbulent kinetic energy, turbulent dissipation rate, turbulent scale length, and viscosity. The results show that these variables attain peak values at the axis of symmetry. Computed distributions of velocity, temperature, and mass fraction are also given. A direct analytical approach to account for the effect of species concentration fluctuations on the mean production rate of species (the phenomenon of unmixedness) is also presented. However, the use of the method does not seem justified in view of the excessive computer time required to solve the resulting system of equations.

Sislian, J. P.↗
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.↗
NASA NTRS1975

A finite-difference program for stresses in anisotropic, layered plates in bending

The interlaminar stresses induced in a layered laminate that is bent into a cylindrical surface are studied. The laminate is modeled as a continuum, and the resulting elasticity equations are solved using the finite difference method. The report sets forth the mathematical framework, presents some preliminary results, and provides a listing and explanation of the computer program. Significant among the results are apparent symmetry relationships that will reduce the numerical size of certain problems and an interlaminar stress behavior having a sharp rise at the free edges.

Salamon, N. J.↗
NASA NTRS1975

An investigation of several numerical procedures for time-asymptotic compressible Navier-Stokes solutions

The status of an investigation of four numerical techniques for the time-dependent compressible Navier-Stokes equations is presented. Results for free shear layer calculations in the Reynolds number range from 1000 to 81000 indicate that a sequential alternating-direction implicit (ADI) finite-difference procedure requires longer computing times to reach steady state than a low-storage hopscotch finite-difference procedure. A finite-element method with cubic approximating functions was found to require excessive computer storage and computation times. A fourth method, an alternating-direction cubic spline technique which is still being tested, is also described.

Rudy, D. H.↗
NASA NTRS1978

Numerical solution of the Navier-Stokes equations for super-sonic flows with strong shocks

The numerical solution of the full Navier-Stokes Equations for viscous flows with high Mach numbers and a strong detached bow shock was obtained. Two dimensional flows around a circular cylinder, and a circular cylinder with an aft-body in the form of a fairing, were considered. The solution of the compressible N.S. equations was accomplished by the method of finite differences. An implicit scheme of solution, the S.O.R., was used with the optimum acceleration parameters determined by trial and error. The tensor notation was used in writing the N-S Equations transformed into general curvilinear coordinates. The equations for the generation of the coordinate system were solved, followed by the solution of the N.S. equations, at the end of a set of given number of time steps. "Wiggles", constituted the one major problem that needed to be overcome. These oscillations give rise to quantities such as negative temperatures, which ultimately caused the computational program to break down. Certain dissipative finite-difference schemes damped these oscillations.

Devarayalu, K.↗
NASA NTRS1978

Numerical method for boundary layers with blowing - The exponential box scheme

The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.

El-Mistikawy, T. M.↗