Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical algorithms”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Hydrodynamic Modeling of Oxidizer-Rich Staged Combustion Injector Flow

The main objective of this work is to determine the unsteady hydrodynamic characteristics of coaxial swirl atomizers of interest in oxidizer-rich staged combustion (ORSC) liquid rocket engines. To this end, the pseudo-density (homogeneous flow) treatment combined with the Marker-and-Cell (MAC) numerical algorithm has been used to develop an axisymmetric with swirl, two-phase, unsteady model. The numerical model is capable of assessing the time-dependent orifice exit conditions and internal mixing for arbitrary fuel and oxidizer gas injection conditions. Parametric studies have been conducted to determine the effect of geometry, gas properties, and liquid properties on the exit massflow rate and velocity. It has been found that the frequency at which the liquid film oscillates increases as the density ratio and thickness increase, decreases as film thickness and liquid swirl velocity increase, and is unaffected by the mixing length. Additionally, it has been determined that the variation in the massflow rate increases as the liquid swirl velocity and liquid film thickness increase, and decreases as the density ratio, collar thickness, and mixing length increase.

Canino, J. V.↗
NASA NTRS1997

Efficient Parallel Algorithm For Direct Numerical Simulation of Turbulent Flows

A distributed algorithm for a high-order-accurate finite-difference approach to the direct numerical simulation (DNS) of transition and turbulence in compressible flows is described. This work has two major objectives. The first objective is to demonstrate that parallel and distributed-memory machines can be successfully and efficiently used to solve computationally intensive and input/output intensive algorithms of the DNS class. The second objective is to show that the computational complexity involved in solving the tridiagonal systems inherent in the DNS algorithm can be reduced by algorithm innovations that obviate the need to use a parallelized tridiagonal solver.

Moitra, Stuti↗
NASA NTRS1995

Development of a linearized unsteady Euler analysis for turbomachinery blade rows

A linearized unsteady aerodynamic analysis for axial-flow turbomachinery blading is described in this report. The linearization is based on the Euler equations of fluid motion and is motivated by the need for an efficient aerodynamic analysis that can be used in predicting the aeroelastic and aeroacoustic responses of blade rows. The field equations and surface conditions required for inviscid, nonlinear and linearized, unsteady aerodynamic analyses of three-dimensional flow through a single, blade row operating within a cylindrical duct, are derived. An existing numerical algorithm for determining time-accurate solutions of the nonlinear unsteady flow problem is described, and a numerical model, based upon this nonlinear flow solver, is formulated for the first-harmonic linear unsteady problem. The linearized aerodynamic and numerical models have been implemented into a first-harmonic unsteady flow code, called LINFLUX. At present this code applies only to two-dimensional flows, but an extension to three-dimensions is planned as future work. The three-dimensional aerodynamic and numerical formulations are described in this report. Numerical results for two-dimensional unsteady cascade flows, excited by prescribed blade motions and prescribed aerodynamic disturbances at inlet and exit, are also provided to illustrate the present capabilities of the LINFLUX analysis.

Verdon, Joseph M.↗
NASA NTRS1990

Flowfield computation of entry vehicles

The equations governing the multidimensional flow of a reacting mixture of thermally perfect gasses were derived. The modeling procedures for the various terms of the conservation laws are discussed. A numerical algorithm, based on the finite-volume approach, to solve these conservation equations was developed. The advantages and disadvantages of the present numerical scheme are discussed from the point of view of accuracy, computer time, and memory requirements. A simple one-dimensional model problem was solved to prove the feasibility and accuracy of the algorithm. A computer code implementing the above algorithm was developed and is presently being applied to simple geometries and conditions. Once the code is completely debugged and validated, it will be used to compute the complete unsteady flow field around the Aeroassist Flight Experiment (AFE) body.

Prabhu, Dinesh K.↗
NASA NTRS2001

High-Performance Java Codes for Computational Fluid Dynamics

The computational science community is reluctant to write large-scale computationally -intensive applications in Java due to concerns over Java's poor performance, despite the claimed software engineering advantages of its object-oriented features. Naive Java implementations of numerical algorithms can perform poorly compared to corresponding Fortran or C implementations. To achieve high performance, Java applications must be designed with good performance as a primary goal. This paper presents the object-oriented design and implementation of two real-world applications from the field of Computational Fluid Dynamics (CFD): a finite-volume fluid flow solver (LAURA, from NASA Langley Research Center), and an unstructured mesh adaptation algorithm (2D_TAG, from NASA Ames Research Center). This work builds on our previous experience with the design of high-performance numerical libraries in Java. We examine the performance of the applications using the currently available Java infrastructure and show that the Java version of the flow solver LAURA performs almost within a factor of 2 of the original procedural version. Our Java version of the mesh adaptation algorithm 2D_TAG performs within a factor of 1.5 of its original procedural version on certain platforms. Our results demonstrate that object-oriented software design principles are not necessarily inimical to high performance.

Riley, Christopher↗
DOE OSTI2022

Estimation of the Water Table Position in Unconfined Aquifers with MODFLOW 6

Abstract The numerical estimation of the position of the water table in unconfined aquifers is important for many practical applications. Its determination through observations or analytical methods is restricted to a few cases. Therefore, it is often estimated through numerical simulations, which may be affected by numerical artifacts and/or poor stability. We use MODFLOW to estimate the position of the water table for a seemingly simple example problem and demonstrate difficulties that can be faced when performing this kind of numerical simulation. We explain the causes for the numerical challenges that originate from the properties of the mathematical equations that must be solved. Based on the results of more than 600 steady‐state simulations, we show how the stability of the numerical solution can be affected by the values of physical parameters that define the problem (e.g., recharge rate, anisotropy ratio, and other parameters that control the numerical algorithm such as settings of the linear and nonlinear solution methods). Finally, we comment on some best practices to apply numerical simulations to estimate the water table position.

Geology↗
NASA NTRS1983

Three-dimensional Euler solutions for long-duct nacelles

A three-dimensional Euler-equation computational technique has been developed to solve for the transonic flow past flow-through nacelles. The technique employs an approximately-factored alternating-direction implicit numerical algorithm and a radiation treatment of the outflow boundary. Studies are presented which show that the radiation treatment gives better numerical convergence than the condition of specifying the pressure at the outflow boundary. Calculations made with the technique are presented for a long-duct turbofan engine nacelle at a Mach number of 0.80 and angles of attack of 0 deg and 4 deg. Good agreement is shown between the computational results and wind-tunnel data. Problem areas are identified and recommendations are made for further numerical studies.

Compton, W. B., III↗
NASA NTRS1978

Supersonic flow over ablated nosetips using an unsteady, implicit numerical procedure

The axisymmetric supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical region, containing multiple zones of embedded subsonic flow into a rectangular computational domain at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing experimental data for very mildly and severely indented blunt nosetips.

Kutler, P.↗
NASA NTRS1980

Supersonic flow over three-dimensional ablated nosetips using an unsteady implicit numerical procedure

The three-dimensional supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical domain, containing multiple zones of embedded subsonic flow and separated flow regions into a rectangular computational volume at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing numerical solutions and experimental data for simple spheres which validates the program. Results are also presented for analytically defined indented bodies for both laminar and turbulent flow conditions that demonstrate the program's capability for computing such flows.

Kutler, P.↗
NASA NTRS1988

Implementation and analysis of a Navier-Stokes algorithm on parallel computers

The results of the implementation of a Navier-Stokes algorithm on three parallel/vector computers are presented. The object of this research is to determine how well, or poorly, a single numerical algorithm would map onto three different architectures. The algorithm is a compact difference scheme for the solution of the incompressible, two-dimensional, time-dependent Navier-Stokes equations. The computers were chosen so as to encompass a variety of architectures. They are the following: the MPP, an SIMD machine with 16K bit serial processors; Flex/32, an MIMD machine with 20 processors; and Cray/2. The implementation of the algorithm is discussed in relation to these architectures and measures of the performance on each machine are given. The basic comparison is among SIMD instruction parallelism on the MPP, MIMD process parallelism on the Flex/32, and vectorization of a serial code on the Cray/2. Simple performance models are used to describe the performance. These models highlight the bottlenecks and limiting factors for this algorithm on these architectures. Finally, conclusions are presented.

Fatoohi, Raad A.↗
NASA NTRSDate not supplied

Investigating the Feasibility and Stability for Modeling Broadband Acoustic Wave Scattering using a Time-Domain Boundary Integral Equation

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. Acoustic liners are an effective tool for reducing aircraft noise and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model is discussed in which the liner impedance is specified along a wide range of frequencies. An impedance boundary condition is derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate. It is assumed the flat plate has surfaces treated with acoustic liners. The stability of the numerical algorithm is assessed using eigenvalue analysis.

mathmatics↗
NASA NTRS1990

Numerical approach of collision avoidance and optimal control on robotic manipulators

Collision-free optimal motion and trajectory planning for robotic manipulators are solved by a method of sequential gradient restoration algorithm. Numerical examples of a two degree-of-freedom (DOF) robotic manipulator are demonstrated to show the excellence of the optimization technique and obstacle avoidance scheme. The obstacle is put on the midway, or even further inward on purpose, of the previous no-obstacle optimal trajectory. For the minimum-time purpose, the trajectory grazes by the obstacle and the minimum-time motion successfully avoids the obstacle. The minimum-time is longer for the obstacle avoidance cases than the one without obstacle. The obstacle avoidance scheme can deal with multiple obstacles in any ellipsoid forms by using artificial potential fields as penalty functions via distance functions. The method is promising in solving collision-free optimal control problems for robotics and can be applied to any DOF robotic manipulators with any performance indices and mobile robots as well. Since this method generates optimum solution based on Pontryagin Extremum Principle, rather than based on assumptions, the results provide a benchmark against which any optimization techniques can be measured.

Wang, Jyhshing Jack↗
NASA NTRS1982

Analysis of a finite difference grid

Some means of assessing the suitability of a mesh network for a finite difference calculation are investigated in this study. This has been done by a study of the nonlinear truncation errors of the scheme. It turns out that the mesh can not be properly assessed a priori. The effect of the mesh on the numerical solution depends on several factors including the mesh itself, the numerical algorithm, and the solution. Several recommendations are made with regard to generating the mesh and to assessing its suitability for a particular numerical calculation.

Klopfer, G. H.↗
DOE OSTI2024

Development of High-Fidelity and Efficient Modeling Capabilities for Enabling Co-Optimization of Fuels and Multi-Mode Engines

The objective of this closely coordinated research project is to develop improved physical sub-models and innovative numerical algorithms to enable accurate and efficient predictions of multi-mode combustion to support EERE’s Co-Optima program. To this end, a research pro-gram is proposed that seeks • to develop accurate submodels for predicting multi-mode combustion regimes, wall-heat transfer, non-equilibrium plasma ignition, and combustion-mode transition, • to develop dynamic load-rebalancing algorithms and efficient semi-implicit time-integration schemes for enabling exascale computing, • to validate the computational submodels against experiments and direct numerical simulation (DNS) data in conjunction with sensitivity analysis, and • to support the direct transition of the resulting sub-models into other engine-simulation tools through the development of self-contained modules. All sub-models will be integrated into the open-source exascale engine-simulation platform Nek5000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1974

Eigenproblem solution of damped structural systems

This article presents an efficient numerical algorithm and a complete listing of the associated computer program, developed for the eigenproblem solution of discrete damped structures, including spinning ones. The numerically stable procedure is based on a combined Sturm sequence and inverse iteration technique, which fully exploits the banded form of the relevant matrices and proves to be most substantially economical when compared to similar existing softwares. Numerical results are presented for representative structures, solved by the present computer program written in FORTRAN V for the JPL UNIVAC 1108 computer.

Gupta, K. K.↗
NASA NTRS1979

First principles numerical model of avalanche-induced arc discharges in electron-irradiated dielectrics

The model consists of four phases: single electron dynamics, single electron avalanche, negative streamer development, and tree formation. Numerical algorithms and computer code implementations are presented for the first three phases. An approach to developing a code description of fourth phase is discussed. Numerical results are presented for a crude material model of Teflon.

Beers, B. L.↗
NASA NTRS1992

Navier-Stokes analyses of flows over slender airframes

An account is given of the current development status of Navier-Stokes algorithms for the analysis of 3D viscous flows over slender airframes, emphasizing the design and analysis of practical configurations. The pacing items of the numerical algorithms used encompass grid generation, spatial and temporal differencing, and transition and turbulence modeling. Attention is given to the high angle-of-attack prediction capability of the various methods whose results are compared with experimental data; it is in that regime that Navier-Stokes methods surpass the prediction capabilities of flow solvers based on simpler mathematical models.

Thomas, James L.↗
NASA NTRS1994

Dynamics of Numerics and CFD

This lecture attempts to illustrate the basic ideas of how the recent advances in nonlinear dynamical systems theory (dynamics) can provide new insights into the understanding of numerical algorithms used in solving nonlinear differential equations (DEs). Examples will be given of the use of dynamics to explain unusual phenomena that occur in numerics. The inadequacy of the use of linearized analysis for the understanding of long time behavior of nonlinear problems will be illustrated, and the role of dynamics in studying the nonlinear stability, accuracy, convergence property and efficiency of using time- dependent approaches to obtaining steady-state numerical solutions in computational fluid dynamics (CFD) will briefly be explained.

Yee, H. C.↗