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 559 records · Page 31

Transient responses of phosphoric acid fuel cell power plant system

An analytical and computerized study of the steady state and transient response of a phosphoric acid fuel cell (PAFC) system was completed. Parametric studies and sensitivity analyses of the PAFC system's operation were accomplished. Four non-linear dynamic models of the fuel cell stack, reformer, shift converters, and heat exchangers were developed based on nonhomogeneous non-linear partial differential equations, which include the material, component, energy balance, and electrochemical kinetic features. Due to a lack of experimental data for the dynamic response of the components only the steady state results were compared with data from other sources, indicating reasonably good agreement. A steady state simulation of the entire system was developed using, nonlinear ordinary differential equations. The finite difference method and trial-and-error procedures were used to obtain a solution. Using the model, a PAFC system, that was developed under NASA Grant, NCC3-17, was improved through the optimization of the heat exchanger network. Three types of cooling configurations for cell plates were evaluated to obtain the best current density and temperature distributions. The steady state solutions were used as the initial conditions in the dynamic model. The transient response of a simplified PAFC system, which included all of the major components, subjected to a load change was obtained. Due to the length of the computation time for the transient response calculations, analysis on a real-time computer was not possible. A simulation of the real-time calculations was developed on a batch type computer. The transient response characteristics are needed for the optimization of the design and control of the whole PAFC system. All of the models, procedures and simulations were programmed in Fortran and run on IBM 370 computers at Cleveland State University and the NASA Lewis Research Center.

Lu, Cheng-Yi↗

Application of a new finite difference algorithm for computational aeroacoustics

Acoustic problems have become extremely important in recent years because of research efforts such as the High Speed Civil Transport program. Computational aeroacoustics (CAA) requires a faithful representation of wave propagation over long distances, and needs algorithms that are accurate and boundary conditions that are unobtrusive. This paper applies a new finite difference method and boundary algorithm to the Linearized Euler Equations (LEE). The results demonstrate the ability of a new fourth order propagation algorithm to accurately simulate the genuinely multidimensional wave dynamics of acoustic propagation in two space dimensions with the LEE. The results also show the ability of a new outflow boundary condition and fourth order algorithm to pass the evolving solution from the computational domain with no perceptible degradation of the solution remaining within the domain.

Goodrich, John W.↗

Numerical simulation of shock/turbulent boundary layer interaction

Most flows of aerodynamic interest are compressible and turbulent. However, our present knowledge on the structures and mechanisms of turbulence is mostly based on incompressible flows. In the present work, compressibility effects in turbulent, high-speed, boundary layer flows are systematically investigated using the Direct Numerical Simulation (DNS) approach. Three-dimensional, time-dependent, fully nonlinear, compressible Navier-Stokes equations were numerically integrated by high-order finite-difference methods; no modeling for turbulence is used during the solution because the available resolution is sufficient to capture the relevant scales. The boundary layer problem deals with fully-turbulent compressible flows over flat geometries. Apart from its practical relevance to technological flows, turbulent compressible boundary layer flow is the simplest experimentally realizable turbulent compressible flow. Still, measuring difficulties prohibit a detailed experimental description of the flow, especially in the near-wall region. DNS studies provide a viable means to probe the physics of compressible turbulence in this region. The focus of this work is to explore the paths of energy transfer through which compressible turbulence is sustained. The structural similarities and differences between the incompressible and compressible turbulence are also investigated. The energy flow patterns or energy cascades are found to be directly related to the evolution of vortical structures which are generated in the near-wall region. Near-wall structures, and mechanisms which are not readily accessible through physical experiments are analyzed and their critical role on the evolution and the behavior of the flow is documented extensively.

Biringen, Sedat↗

Solution of the Chapman-Ferraro problem with an arbitrary magnetopause

We present a global model of the magnetic field of the magnetosphere that includes the effects of the Chapman-Ferraro currents at the magnetopause. In contrast to ealier models, the magnetopause shape is arbitrary, thus allowing the use of more realistic geometries. The internal magnetospheric field model of Hilmer and Voigt (1993), is completely shielded within the magnetopause by solving the Laplace equation with Neumann boundary conditions using a finite difference method on a non-orthogonal, curvilinear grid. The resulting model magnetosphere is perfectly closed although the method can also be applied with more general boundary conditions, to generate a set of open models based on the approach of Toffoletto and Hill (1989, 1993). The purpose of this paper is to demonstrate the feasibility of a purely numerical approach to solving the Chapman-Ferraro problem with arbitrary magnetopause shape and boundary conditions.

Toffoletto, F. R.↗

Adaptively-refined overlapping grids for the numerical solution of systems of hyperbolic conservation laws

Adaptive mesh refinement (AMR) in conjunction with higher-order upwind finite-difference methods have been used effectively on a variety of problems in two and three dimensions. In this paper we introduce an approach for resolving problems that involve complex geometries in which resolution of boundary geometry is important. The complex geometry is represented by using the method of overlapping grids, while local resolution is obtained by refining each component grid with the AMR algorithm, appropriately generalized for this situation. The CMPGRD algorithm introduced by Chesshire and Henshaw is used to automatically generate the overlapping grid structure for the underlying mesh.

Brislawn, Kristi D.↗

Local multiplicative Schwarz algorithms for convection-diffusion equations

We develop a new class of overlapping Schwarz type algorithms for solving scalar convection-diffusion equations discretized by finite element or finite difference methods. The preconditioners consist of two components, namely, the usual two-level additive Schwarz preconditioner and the sum of some quadratic terms constructed by using products of ordered neighboring subdomain preconditioners. The ordering of the subdomain preconditioners is determined by considering the direction of the flow. We prove that the algorithms are optimal in the sense that the convergence rates are independent of the mesh size, as well as the number of subdomains. We show by numerical examples that the new algorithms are less sensitive to the direction of the flow than either the classical multiplicative Schwarz algorithms, and converge faster than the additive Schwarz algorithms. Thus, the new algorithms are more suitable for fluid flow applications than the classical additive or multiplicative Schwarz algorithms.

Cai, Xiao-Chuan↗

A fast Poisson solver for unsteady incompressible Navier-Stokes equations on the half-staggered grid

In this paper, a fast Poisson solver for unsteady, incompressible Navier-Stokes equations with finite difference methods on the non-uniform, half-staggered grid is presented. To achieve this, new algorithms for diagonalizing a semi-definite pair are developed. Our fast solver can also be extended to the three dimensional case. The motivation and related issues in using this second kind of staggered grid are also discussed. Numerical testing has indicated the effectiveness of this algorithm.

Golub, G. H.↗

Free-Surface and Contact Line Motion of Liquid in Microgravity

This project involves fundamental studies of the role of nonlinearity in determining the motion of liquid masses under the principal influences of surface tension, viscosity and inertia. Issues to be explored are relevant to aspects of terrestrial processes, as well as being immediately applicable to fluid management in a low-gravity environment. Specific issues include: (1) the mechanic's of liquid masses in large-amplitude motions, (2) the influence of bounding surfaces on the motion, and (3) the ability of such surfaces to control liquid motion by wetting forces, especially when they are augmented by various surface treatments. Mathematical techniques include asymptotic analysis of the governing equations, for problem simplification, and numerical simulation, using both boundary-element and finite-difference methods. The flow problem is divided into an 'outer' or inviscid potential-flow region and one or more inner, or viscous dominated, regions. Relevant to one inner region, the vicinity of the contact line, we discuss time-dependent simulation of slow droplet motion, on a surface of variable wettability, using the lubrication approximation. The simulation uses a disjoining pressure model and reproduces realistic wetting-dewetting behavior.

Schwartz, Leonard W.↗

A 2.5D Single Passage CFD Model for Centrifugal Pumps

This paper describes the single passage model based on CFD to analyze the flow in blade passages of a centrifugal pump. The model consists of the flow passage between two impeller blades and the spaces in the inlet eye as well as in the volute. The incompressible Navier-Stokes equations in the conservation form are solved by a finite difference method. The code is designed to investigate the velocity and pressure distributions and intended to investigate how the pump design affects fluid flow through the rotor as well as the pump performance. An early part of the paper investigates the behavior of the model as well as validity of the assumptions made in the model. Then, applications to a rotodynamic heart pump are presented.

Nakamura S.↗

High Order Numerical Simulation of Sound Generated by the Kirchhoff Vortex

An improved high order finite difference method for low Mach number computational aeroacoustics (CAA) is described. The improvements involve the conditioning of the Euler equations in perturbation form to minimize numerical cancellation error, and the use of a stable non-dissipative sixth-order central spatial differencing for the interior points and third-order at the boundary points. The spatial difference operator satisfies the summation-by-parts property to guarantee strict stability for linear hyperbolic systems. Spurious high frequency oscillations are damped by a third-order characteristic-based filter. The objective of this paper is to apply these improvements in the simulation of sound generated by the Kirchhoff vortex.

Mueller, Bernhard↗

Accurate Evaluation of Quantum Integrals

Combining an appropriate finite difference method with Richardson's extrapolation results in a simple, highly accurate numerical method for solving a Schr\"{o}dinger's equation. Important results are that error estimates are provided, and that one can extrapolate expectation values rather than the wavefunctions to obtain highly accurate expectation values. We discuss the eigenvalues, the error growth in repeated Richardson's extrapolation, and show that the expectation values calculated on a crude mesh can be extrapolated to obtain expectation values of high accuracy.

Galant, David C.↗

Accurate Evaluation of Quantum Integrals

Combining an appropriate finite difference method with Richardson's extrapolation results in a simple, highly accurate numerical method for solving a Schrodinger's equation. Important results are that error estimates are provided, and that one can extrapolate expectation values rather than the wavefunctions to obtain highly accurate expectation values. We discuss the eigenvalues, the error growth in repeated Richardson's extrapolation, and show that the expectation values calculated on a crude mesh can be extrapolated to obtain expectation values of high accuracy.

Galant, D. C.↗

Interaction of Porosity with a Planar Solid/Liquid Interface

In this paper, an investigation of the interaction between gas porosity and a planar solid/liquid (SL) interface is reported. A two-dimensional numerical model able to accurately track sharp SL interfaces during solidification of pure metals and alloys is proposed. The finite difference method and a rectangular undeformed grid are used for computation. The SL interface is described through the points of intersection with the grid lines. Its motion is determined by the thermal and solute gradients at each particular point. Changes of the interface temperature because of capillarity or solute redistribution as well as any perturbation of the thermal and solute field produced by the presence of non-metallic inclusions can be computed. To validate the model, the dynamics of the interaction between a gas pore and a solidification front in metal alloys was observed using a state of the art X-ray Transmission Microscope. The experiments included observation of the distortion of the SL interface near a pore, real-time measurements of the growth rate and the change in shape of the porosity during interaction with an advancing SL interface in pure Al and Al-0.25 wt% Au alloy. In addition, porosity induced solute segregation patterns surrounding a pore were also quantified.

Catalina, Adrian V.↗

Transverse Mode Dynamics of VCSELs Through Space-Time Domain Simulation

Modeling and simulation are important to understand laser operation and to optimize and design device functions. Numerical simulation of VCSEL (Vertical Cavity Surface Emitting Lasers) has been largely based on solving time-independent Helmholtz equation or time dependent coupled mode equations. There are various advantages for choosing these approaches. However, the disadvantages are also apparent. The former cannot handle dynamical mode competition seen in VCSELs, while the latter assumes a given type and number of modes a priori. Furthermore, the microscopic physics of heterstructures and electron-hole plasma is very often represented by a few parameters such as linear gain coefficients and the linewidth enhancement factor. These are over simplification of space and frequency (wavelength) dependent gain and refractive index functions. When the space-time dynamical operation of VCSELs becomes important, these simple approximations become questionable. In this paper, we apply a recently developed model for edge-emitting lasers to a gain guided VCSEL for space-time domain simulation. This model takes into account the actual nonlinear dependence of gain and refractive index on frequency and carrier density within the frame work of the effective Bloch equations. The corresponding partial differential equations are solved directly by finite difference methods. Laser behavior with increasing pumping current is investigated in detail. Special attention is paid to the dynamical competition of the transverse modes.

Goorjian, Peter M.↗

Closed-Loop Endoatmospheric Ascent Guidance

This paper will present a complete formulation of the optimal control problem for atmospheric ascent of rocket powered launch vehicles subject to usual load constraints and final condition constraints. We shall demonstrate that the classical finite difference method for two-point-boundary-value-problems (TPBVP) is suited for solving the ascent trajectory optimization problem in real time, therefore closed-loop optimal endoatmospheric ascent guidance becomes feasible. Numerical simulations with a the vehicle data of a reusable launch vehicle will be provided.

Lu, Ping↗

Interaction of Porosity with a Planar Solid/Liquid Interface

In this article, an investigation of the interaction between gas porosity and a planar solid/liquid (SL) interface is reported. A two-dimensional numerical model able to accurately track sharp SL interfaces during solidification of pure metals and alloys is proposed. The finite-difference method and a rectangular undeformed grid are used for computation. The SL interface is described through the points of intersection with the grid lines. Its motion is determined by the thermal and solute gradients at each particular point. Changes of the interface temperature because of capillarity or solute redistribution as well as any perturbation of the thermal and solute field produced by the presence of non-metallic inclusions can be computed. To validate the model, the dynamics of the interaction between a gas pore and a solidification front in metal alloys was observed using a state of the art X-ray transmission microscope (XTM). The experiments included observation of the distortion of the SL interface near a pore, real-time measurements of the growth rate, and the change in shape of the porosity during interaction with the SL interface in pure Al and Al-0.25 wt pct Au alloy. In addition, porosity-induced solute segregation patterns surrounding a pore were also quantified.

Catalina, Adrian V.↗

Scalability of Parallel Spatial Direct Numerical Simulations on Intel Hypercube and IBM SP1 and SP2

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube and IBM SP1 and SP2 parallel computers is documented. Spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows are computed with the PSDNS code. The feasibility of using the PSDNS to perform transition studies on these computers is examined. The results indicate that PSDNS approach can effectively be parallelized on a distributed-memory parallel machine by remapping the distributed data structure during the course of the calculation. Scalability information is provided to estimate computational costs to match the actual costs relative to changes in the number of grid points. By increasing the number of processors, slower than linear speedups are achieved with optimized (machine-dependent library) routines. This slower than linear speedup results because the computational cost is dominated by FFT routine, which yields less than ideal speedups. By using appropriate compile options and optimized library routines on the SP1, the serial code achieves 52-56 M ops on a single node of the SP1 (45 percent of theoretical peak performance). The actual performance of the PSDNS code on the SP1 is evaluated with a "real world" simulation that consists of 1.7 million grid points. One time step of this simulation is calculated on eight nodes of the SP1 in the same time as required by a Cray Y/MP supercomputer. For the same simulation, 32-nodes of the SP1 and SP2 are required to reach the performance of a Cray C-90. A 32 node SP1 (SP2) configuration is 2.9 (4.6) times faster than a Cray Y/MP for this simulation, while the hypercube is roughly 2 times slower than the Y/MP for this application. KEY WORDS: Spatial direct numerical simulations; incompressible viscous flows; spectral methods; finite differences; parallel computing.

Joslin, Ronald D.↗

Parallel Computation for Natural Convection in Cavities

Parallel computation for thermal convective flows in cavities with adiabatic horizontal boundaries and driven by differential heating of the two vertical end walls, is investigated using supercomputers. A parallel computation code has been implemented using a finite-difference method with a multigrid elliptic solver and a Dufort-Frankel scheme. The domain decomposition techniques are discussed in detail. The parallel code is numerically stable, computationally efficient, and portable to various parallel architectures which support either PVM or NX libraries for communications. Finally, numerical results for various Rayleigh numbers and Prandtl numbers are presented.

convection simulation parallel computation↗