Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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 253 records · Page 14

Numerical simulation of steady transonic flow about airfoils

A computer code has been developed that couples a fast transonic full-potential AF2 solver with both an efficient integral boundary-layer method and a viscous wedge approximation of the shock/boundary-layer interaction. The efficiency of the coupled analysis methods and the method of coupling has resulted in a uniquely efficient analysis tool. The airfoil geometry is modified by the displacement thickness before the shock and the displacement thickness plus the viscous wedge thickness after the shock by considering the viscous effects as an equivalent transpiration boundary condition. The flow about conventional and supercritical airfoils under moderately strong shock situations has been calculated. Comparisons with experimental data indicate that this viscous correction method has improved the accuracy of the full-potential analysis. Furthermore, the computer time required to obtain a converged solution has been reduced.

Lee, S. C.↗

Efficiency of silicon solar cells as a function of base layer resistivity

This paper reports on a theoretical study of the limitations on silicon solar-cell efficiency for both n(+)-p and n(+)-p-p(+) type cells. Detailed calculations have been made of solar-cell operation using a general computer analysis program for semiconductor devices. The computer program, which simultaneously solves Poisson's equation and the electron and hole quasi-Fermi level equations, provides an accurate numerical solution of solar-cell operation without limiting assumptions or approximations. It is found that minority-carrier lifetime and heavy doping effects in the n(+) surface region present serious limitations to efficiency in low-resistivity silicon solar cells.

Dunbar, P. M.↗

Accuracy of relative oceanic geoid computations

An error analysis based on a model which requires potential coefficients and gravity anomaly information for the computation of oceanic geoidal undulations is developed. The rigorous models are presented and abandoned owing to the numerical difficulties associated with them. An approximate error model is developed in the frequency domain through which an idea is obtained as to the size and behavior of the errors involved. The use of the method is demonstrated by obtaining data requirements for the realization of a 10-cm relative oceanic geoid. Various data sets could result in the desired accuracy. One such set, indicative of the strictness of the requirements, involves gravity profile spacings of approximately 3 km with observational noise not exceeding 0.5 mGal inside detailed data caps of 30 deg and potential coefficients of prescribed accuracy available to degree and order 70.

Christodoulidis, D. C.↗

Low frequency acoustic and electromagnetic scattering

This paper deals with two classes of problems arising from acoustics and electromagnetics scattering in the low frequency stations. The first class of problem is solving Helmholtz equation with Dirichlet boundary conditions on an arbitrary two dimensional body while the second one is an interior-exterior interface problem with Helmholtz equation in the exterior. Low frequency analysis show that there are two intermediate problems which solve the above problems accurate to 0(k(2) log k) where k is the frequency. These solutions greatly differ from the zero frequency approximations. For the Dirichlet problem numerical examples are shown to verify the theoretical estimates.

Hariharan, S. I.↗

High Speed Networking and Large-scale Simulation in Geodynamics

Large-scale numerical simulation has been one of the most important approaches for understanding global geodynamical processes. In this approach, peta-scale floating point operations (pflops) are often required to carry out a single physically-meaningful numerical experiment. For example, to model convective flow in the Earth's core and generation of the geomagnetic field (geodynamo), simulation for one magnetic free-decay time (approximately 15000 years) with a modest resolution of 150 in three spatial dimensions would require approximately 0.2 pflops. If such a numerical model is used to predict geomagnetic secular variation over decades and longer, with e.g. an ensemble Kalman filter assimilation approach, approximately 30 (and perhaps more) independent simulations of similar scales would be needed for one data assimilation analysis. Obviously, such a simulation would require an enormous computing resource that exceeds the capacity of a single facility currently available at our disposal. One solution is to utilize a very fast network (e.g. 10Gb optical networks) and available middleware (e.g. Globus Toolkit) to allocate available but often heterogeneous resources for such large-scale computing efforts. At NASA GSFC, we are experimenting with such an approach by networking several clusters for geomagnetic data assimilation research. We shall present our initial testing results in the meeting.

Kuang, Weijia↗

Approximate analysis of thermal convection in a crystal-growth cell for Spacelab 3

The transient and steady thermal convection in microgravity is described. The approach is applicable to many three dimensional flows in containers of various shapes with various thermal gradients imposed. The method employs known analytical solutions to two dimensional thermal flows in simpler geometries, and does not require recourse to numerical calculations by computer.

Dressler, R. F.↗

The modified equation approach to the stability and accuracy analysis of finite-difference methods

The stability and accuracy of finite-difference approximations to simple linear partial differential equations are analyzed by studying the modified partial differential equation. Aside from round-off error, the modified equation represents the actual partial differential equation solved when a numerical solution is computed using a finite-difference equation. The modified equation is derived by first expanding each term of a difference scheme in a Taylor series and then eliminating time derivatives higher than first order by certain algebraic manipulations. The connection between 'heuristic' stability theory based on the modified equation approach and the von Neumann (Fourier) method is established. In addition to the determination of necessary and sufficient conditions for computational stability, a truncated version of the modified equation can be used to gain insight into the nature of both dissipative and dispersive errors.

Warming, R. F.↗

Improved Euler analysis of advanced turboprop propeller flows

An implicit approximate factorization scheme in conjunction with a new boundary treatment was used to compute the inviscid flow field about an advanced high-speed propeller. The method of characteristics was used to apply impermeable boundary conditions. The convergence history of numerical calculations shows substantial decrease in the residual error decay. Euler solutions were computed for SR-3 propeller geometry and the results were compared with the experimental data and previous numerical results.

Yamamoto, O.↗

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗

Numerical analysis of discrete fractional integrodifferential structural dampers

This paper develops solution algorithms enabling the handling of the dynamic response of nonlinear structures contained discretely attached dampers modeled by fractional integrodifferential operators of the Grunwald-Liouville-Riemann type. The development consists of two levels of formulation, namely: (1) numerical approximations of fractional operators and, (2) the establishment of global level implicit schemes enabling the solution to nonlinear structural formulations. To generalize the overall results, error estimates are derived for the fractional operator approximation algorithm. These enable an ongoing optimization of solution efficiency for a given error tolerance. To benchmark the scheme, the results of several numerical experiments are presented. These illustrate the numerical characteristics of the overall formulation.

Padovan, Joseph↗

On the relationship between acoustic energy density flux near the jet and far field acoustic intensity

The relationship between the distribution of the outflow of acoustic energy over the jet boundary and the far-field directivity and intensity distribution is established by measurement and analysis. The numerical and experimental procedures involved have been checked out by using a known source. The results indicate that the acoustic power output per unit length of the jet, in the region from which the sound emanates, peaks at approximately 9 diameters downstream. The acoustic emission for a jet Strouhal number of about 0.3 exceeds the emission for all other Strouhal numbers nearly everywhere along the measurement plane. However, the far-field peak intensity distribution obtained from the contribution of each station was found to depend on the spatial extent of the region where sound emanates from the jet, which, in turn, depends more on the far-field angle than on the Strouhal number. The implications of these results for sound suppression techniques are discussed.

Maestrello, L.↗

Structural tailoring of advanced turboprops

A computer program has been developed for the performance of numerical optimizations of highly swept propfan blades by minimizing an objective function that is defined either as direct operating cost or the aeroelastic difference between a blade and its scaled model. Three component analysis categories are employed: an optimization algorithm, approximate analysis procedures for objective function and constraint evaluation, and refined analysis procedures for optimum design validation. The analyses conducted by the program encompass aerodynamic efficiency evaluation, finite element stress and vibration analysis, acoustics, flutter, and forced response life prediction.

Brown, K. W.↗

Comparison of Mesospheric Winds From a High-Altitude Meteorological Analysis System and Meteor Radar Observations During the Boreal Winters of 2009-2010 and 2012-2013

We present a study of horizontal winds in the mesosphere and lower thermosphere (MLT) during the boreal winters of 2009-2010 and 2012-2013 produced with a new high-altitude numerical weather prediction (NWP) system. This system is based on a modified version of the Navy Global Environmental Model (NAVGEM) with an extended vertical domain up to approximately 116 km altitude coupled with a hybrid four-dimensional variational (4DVAR) data assimilation system that assimilates both standard operational meteorological observations in the troposphere and satellite-based observations of temperature, ozone and water vapor in the stratosphere and mesosphere. NAVGEM-based MLT analyzed winds are validated using independent meteor radar wind observations from nine different sites ranging from 69 deg N-67 deg S latitude. Time-averaged NAVGEM zonal and meridional wind profiles between 75 and 95 km altitude show good qualitative and quantitative agreement with corresponding meteor radar wind profiles. Wavelet analysis finds that the 3-hourly NAVGEM and 1-hourly radar winds both exhibit semi-diurnal, diurnal, and quasi-diurnal variations whose vertical profiles of amplitude and phase are also in good agreement. Wavelet analysis also reveals common time-frequency behavior in both NAVGEM and radar winds throughout the Northern extra tropics around the times of major stratospheric sudden warmings (SSWs) in January 2010 and January 2013, with a reduction in semi-diurnal amplitudes beginning around the time of a mesospheric wind reversal at 60 deg N that precedes the SSW, followed by an amplification of semi-diurnal amplitudes that peaks 10-14 days following the onset of the mesospheric wind reversal. The initial results presented in this study demonstrate that the wind analyses produced by the high altitude NAVGEM system accurately capture key features in the observed MLT winds during these two boreal winter periods.

McCormack, J.↗

NASA-UCLA Workshop on Computational Techniques in Identification and Control of Flexible Flight Structures, Lake Arrowhead, CA, Nov. 2-4, 1989, Proceedings

The present conference on the identification and control of flexible flight structures with computational techniques encompasses existing and planned testbeds, modeling techniques, adaptive control, and numerical computation and simulation. Specific issues addressed include the NASA/MFSC ground test facility, an active control test on the NASA Minimast, the Astrex testbed for large precision space structures, the dynamic analysis of truss-beam modeling, combined structure-controls-integrated optimization, distributed parameter modeling applicable to the control of flexible flight structures, and approaches to nonclassically damped dynamic systems. Also addressed are adaptive control concepts, techniques for vibration and decentralized slew maneuver control, the discretization of second-order-in-time differential equations which are applicable to nonlinear wave problems, the finite-dimensional approximations of Riccati equations, and numerical simulations of actively controlled space structures.

Balakrishnan, A. V.↗

Reflection and transmission of acoustic waves from a moving layer

The refraction of acoustic waves by a moving medium layer is theoretically treated and the expressions for reflection and transmission coefficients are determined. The moving medium layer velocity is assumed to have a space dependence in one direction. A partitioning of the moving medium layer into constant-velocity sublayers is introduced and the number of sublayers is allowed to increase until the reflection and transmission coefficients converage to their respective values. Numerical results for several sublayer approximations of Poiseuille's flow are presented as functions of the moving layer velocity for several angles of incidence of the acoustic wave. The degenerate case of single constant-velocity layer is also treated, both theoretically and by a numerical analysis.

Steinmetz, G. G.↗

A perturbation solution for transverse pressure-sensitive combustion instability in an annulus

A method for analysis of combustion instability in rocket motors based on a combination of the Galerkin method and the two-variable (or multiple-scales) perturbation method is developed. The method is illustrated by applying it to the problem of pressure-sensitive combustion instability in a liquid-fuel annular combustion chamber. To the order of approximation inherent in the method it is found that the complete linear stability analysis and most of the nonlinear stability analysis can be carried out in closed form. The results thus obtained are used to clarify several aspects of the interpretation of numerical solutions obtained by previous investigators.

Googerdy, A.↗

Analysis and computation of three-dimensional flow in strongly curved ducts

Numerical solutions are presented for three-dimensional laminar and turbulent flow in curved ducts of rectangular cross section and significant curvature. The analysis is based on a primary-secondary velocity decomposition in a given coordinate system, and leads to approximate governing equations which correct an a priori inviscid solution for viscous effects, secondary flows, total pressure distortion, heat transfer, and internal flow blockage and losses. Solution of the correction equations is accomplished as an initial-value problem in space using an implicit forward-marching technique. The overall solution procedure requires significantly less computational effort than Navier-Stokes algorithms. The present solution procedure is effective even with the extreme local mesh resolution which is necessary to resolve near-wall sublayer regions in turbulent flow calculations. Computed solutions for both laminar and turbulent flow compare very favorably with available analytical and experimental results.

Kreskovsky, J. P.↗

Improved tangent-cone method for the Aerodynamic Preliminary Analysis System (APAS) version of the hypersonic arbitrary-body program

The Aerodynamic Preliminary Analysis System (APAS) utilizes a modified version of the Hypersonic Arbitrary-Body Program (HABP) Mark 3 code in its analysis rationale. Four methods are considered for incorporation into the code as the tangent-cone method. The combination of second-order slender body theory and the approximate solution of Hammitt and Murthy shows the best agreement with the exact numerical solutions and is thus included in the APAS production version of the HABP code.

Cruz, Christopher I.↗