SEARCH · Engineering Papers

Results for “numerical solution”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Axially grooved heat pipe study

A technology evaluation study on axially grooved heat pipes is presented. The state-of-the-art is reviewed and present and future requirements are identified. Analytical models, the Groove Analysis Program (GAP) and a closed form solution, were developed to facilitate parametric performance evaluations. GAP provides a numerical solution of the differential equations which govern the hydrodynamic flow. The model accounts for liquid recession, liquid/vapor shear interaction, puddle flow as well as laminar and turbulent vapor flow conditions. The closed form solution was developed to reduce computation time and complexity in parametric evaluations. It is applicable to laminar and ideal charge conditions, liquid/vapor shear interaction, and an empirical liquid flow factor which accounts for groove geometry and liquid recession effects. The validity of the closed form solution is verified by comparison with GAP predictions and measured data.

Source record↗

NASA NTRS1980

Grid generation using differential systems techniques

The errors in approximating the derivatives of a function by traditional central differences at grid points of a curvilinear coordinate system were examined. The implications concerning the accuracy of the numerical solution of a partial differential equation are explained by considering several numerical examples. Although this study only considers the two dimensional case, the techniques and implications are equally valid for three dimensional grids. An interesting feature of the error analysis is its simplicity. Most of the results follow by merely working with the truncation terms of some power series expansion. These series expansions also give rise to higher order difference approximations which can significantly reduce error when the grid spacing changes rapidly, as might be the case in problems with shock waves or thin boundary layers.

Thompson, J. F.↗

NASA NTRS1980

Transient and diffusion analysis of Hg Cd Te

The solution of the one dimensional planar interface solidification problem is presented. Numerical solutions are applied to experimental solute profiles to determine the effective diffusion constant for Hg Co Te. Progress on the design and construction of a high gradient furnace is reported.

Clayton, J. C.↗

NASA NTRS1983

A preconditioned formulation of the Cauchy-Riemann equations

A preconditioning of the Cauchy-Riemann equations which results in a second-order system is described. This system is shown to have a unique solution if the boundary conditions are chosen carefully. This choice of boundary condition enables the solution of the first-order system to be retrieved. A numerical solution of the preconditioned equations is obtained by the multigrid method.

Phillips, T. N.↗

NASA NTRS1985

PIFEX: An advanced programmable pipelined-image processor

PIFEX is a pipelined-image processor being built in the JPL Robotics Lab. It will operate on digitized raster-scanned images (at 60 frames per second for images up to about 300 by 400 and at lesser rates for larger images), performing a variety of operations simultaneously under program control. It thus is a powerful, flexible tool for image processing and low-level computer vision. It also has applications in other two-dimensional problems such as route planning for obstacle avoidance and the numerical solution of two-dimensional partial differential equations (although its low numerical precision limits its use in the latter field). The concept and design of PIFEX are described herein, and some examples of its use are given.

Gennery, D. B.↗

NASA NTRS1988

Analysis of the spectral vanishing viscosity method for periodic conservation laws

The convergence of the spectral vanishing method for both the spectral and pseudospectral discretizations of the inviscid Burgers' equation is analyzed. It is proven that this kind of vanishing viscosity is responsible for a spectral decay of those Fourier coefficients located toward the end of the computed spectrum; consequently, the discretization error is shown to be spectrally small independent of whether the underlying solution is smooth or not. This in turn implies that the numerical solution remains uniformly bounded and convergence follows by compensated compactness arguments.

Maday, Yvon↗

NASA NTRS1988

Computational issues in control-structure interaction analysis

This paper surveys computational issues for large-scale simulation of dynamics and controls of space structures, which involve structural elements that are capable of large combined rigid and flexible motions, accurate and efficient treatment of constraints, robust integration of both translational and large rotational motions in the equations of motion, modular interface with active control synthesis packages, and a capability from wave motions to slowly varying transient responses. A particular feature of the present survey is a partitioned solution (or divide-and-conquer) procedure that can handle the numerical solution of multidisciplinary simulation problems by relying on individual modular solution packages that treat each aspect of simulation requirements.

Park, K. C.↗

NASA NTRS1990

Asymptotic analysis of dissipative waves with applications to their numerical simulation

Various problems involving the interplay of asymptotics and numerics in the analysis of wave propagation in dissipative systems are studied. A general approach to the asymptotic analysis of linear, dissipative waves is developed. It was applied to the derivation of asymptotic boundary conditions for numerical solutions on unbounded domains. Applications include the Navier-Stokes equations. Multidimensional traveling wave solutions to reaction-diffusion equations are also considered. A preliminary numerical investigation of a thermo-diffusive model of flame propagation in a channel with heat loss at the walls is presented.

Hagstrom, Thomas↗

NASA NTRS1991

On the nonlinear development of the most unstable Goertler vortex mode

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

NASA NTRS1992

Boundary layer study on nozzle wall at hypersonic velocities

The boundary layer on the wall of the Hypervelocity Tunnel 9 was investigated with pitot pressure and total temperature measurements. Experimental results are presented for standard and supercooled Mach 14 runs. The boundary layer data at supercooled conditions are compared to numerical predictions made with a Navier-Stokes algorithm including vibrational nonequilibrium and intermolecular force effects. For standard tunnel conditions, the numerical solutions agree well with experimental data. For the supercooled cases, the numerical code predicts the total temperature but overpredicts the pitot pressure.

Jones, Kenneth M.↗

NASA NTRS1993

On the nonlinear development of the most unstable Goertler vortex mode TLSP - Final report

The nonlinear development of the most unstable Goertler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Goertler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Goertler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Goertler equations valid for O(1) Goertler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

NASA NTRS2018

Modeling the Effect of Damage on the Electrical Resistivity of Melt-Infiltrated SiC/SiC Composites

Electrical resistivity (ER) measurements are a possible health monitoring technique for CMC components in future aerospace applications. In order to use ER measurements to detect and identify damage, it is necessary to understand how each specific damage state will affect the ER response. In this paper, finite element models are developed and applied to quantify the effect of specific damage states on the electrical resistivity response in a melt-infiltrated silicon carbide (SiC) fiber-reinforced silicon carbide composite. The electrical resistivity of several damage states are calculated by simulating the electric current flow through the damaged micro-structure. This is achieved by performing the numerical solution of the steady-state conservation of charge density equation. Numerical results reveal that cracking of the inter-tow matrix has the most profound effect on the composite electrical resistivity. Also, fiber/matrix debonding at matrix cracks in the 0° tows (tows aligned with the loading direction) may cause a significant increase in the electrical resistivity, but only if the fiber coating resistivity is 1000 Ω-cm or less. Cracks in the 90° tows and the crack opening displacement have very little effect on the composite electrical resistivity.

finite element analysis↗

NASA NTRS2019

Numerical analysis of unsteady flow in a converging-diverging nozzle

This study investigates the physics of jet noise by the numerical solution of the three dimensional unsteady Navier-Stokes equations for a low-supersonic jet (Mach numbers on the order of 1.4) from a converging-diverging (C-D) nozzle and mixing with a subsonic surrounding flow. The numerical method used for the solution of these equations is the second order accurate MacCormack explicit predictor-corrector scheme. The numerical results are compared to the experimental results of Yamamoto for the same C-D nozzle geometry and similar operating conditions. The goal of the study is to use numerical simulation techniques to resolve accurately physical processes which contribute to jet noise, and to assess whether performing these simulations in a three dimensional geometry is worth the additional cost in computing resources as compared to a comparable simulation on a two dimensional planar or axisymmetric geometry. Scientific visualization is used to observe the spreading characteristics of the jet, as well as to examine the structure of the shocks within the jet core. The time average velocity and turbulence intensity level at various locations are computed and compared with experimental results. A two point space-time cross correlation is used to compute the mean convection speed through the shear layer, and these predictions are compared with results from theory.

James N Scott↗

NASA NTRS1974

Finite element solution theory for three-dimensional boundary flows

A finite element algorithm is derived for the numerical solution of a three-dimensional flow field described by a system of initial-valued, elliptic boundary value partial differential equations. The familiar three-dimensional boundary layer equations belong to this description when diffusional processes in only one coordinate direction are important. The finite element algorithm transforms the original description into large order systems of ordinary differential equations written for the dependent variables discretized at node points of an arbitrarily irregular computational lattice. The generalized elliptic boundary conditions is piecewise valid for each dependent variable on boundaries that need not explicitly coincide with coordinate surfaces. Solutions for sample problems in laminar and turbulent boundary flows illustrate favorable solution accuracy, convergence, and versatility.

Baker, A. J.↗

NASA NTRS1989

Solutions of differential equations with regular coefficients by the methods of Richmond and Runge-Kutta

Numerical solutions of the differential equation which describe the electric field within an inhomogeneous layer of permittivity, upon which a perpendicularly-polarized plane wave is incident, are considered. Richmond's method and the Runge-Kutta method are compared for linear and exponential profiles of permittivities. These two approximate solutions are also compared with the exact solutions.

Cockrell, C. R.↗

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 NTRS1990

Stability analysis of multigrid acceleration methods for the solution of partial differential equations

A calculation is made of the stability of various relaxation schemes for the numerical solution of partial differential equations. A multigrid acceleration method is introduced, and its effects on stability are explored. A detailed stability analysis of a simple case is carried out and verified by numerical experiment. It is shown that the use of multigrids can speed convergence by several orders of magnitude without adversely affecting stability.

Fay, John F.↗

NASA NTRS1981

Solution of the transonic full potential equation in conservative form using an implicit algorithm

The paper presents numerical solutions of the full potential equation in conservative form. The iteration scheme used is a fully implicit approximate factorization technique and provides a significant improvement in convergence speed relative to standard successive line overrelaxation algorithms. The spatial differencing algorithm is centrally differenced in both subsonic and supersonic regions to maintain stability. This effectively approximates rotated differencing, thereby greatly improving the reliability of the algorithm.

Holst, T.↗