Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical 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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗
NASA NTRS1990

Elimination of Gibbs' phenomena from error analysis of finite element results

This paper is one of a series on error analysis and correction of finite element solutions for plates and shells. The error analysis in the earlier papers used half-range double Fourier sine series for numerical harmonic analysis. The half-range formulas are simple to apply, but they can be inaccurate near the ends of the ranges of the independent variables. The Gibbs' phenomenon exhibited by half-range sine series in one independent variable has a two-dimensional analog; a classic example is the Navier solution in a double half-range sine series for the simply supported plate under a uniform load. A simple change of variables is introduced in the paper to improve the accuracy of the double Fourier sine series without adding complexity to the numerical analysis. The change of variables is applied to the problem of approximating a transverse load that is tabulated on a rectangular grid. A solution based on the change of variables is compared with results from the Navier solution for the simply supported plate problem and finite element results for the same problem.

Thurston, Gaylen A.↗
NASA NTRS1986

Formulation of numerical procedures for dynamic analysis of spinning structures

The paper presents the descriptions of recently developed numerical algorithms that prove to be useful for the solution of the free vibration problem of spinning structures. First, a generalized procedure for the computation of nodal centrifugal forces in a finite element owing to any specified spin rate is derived in detail. This is followed by a description of an improved eigenproblem solution procedure that proves to be economical for the free vibration analysis of spinning structures. Numerical results are also presented which indicate the efficacy of the currently developed procedures.

Gupta, K. K.↗
NASA NTRS1995

An efficient numerical procedure for thermohydrodynamic analysis of cavitating bearings

An efficient and accurate numerical procedure to determine the thermo-hydrodynamic performance of cavitating bearings is described. This procedure is based on the earlier development of Elrod for lubricating films, in which the properties across the film thickness are determined at Lobatto points and their distributions are expressed by collocated polynomials. The cavitated regions and their boundaries are rigorously treated. Thermal boundary conditions at the surfaces, including heat dissipation through the metal to the ambient, are incorporated. Numerical examples are presented comparing the predictions using this procedure with earlier theoretical predictions and experimental data. With a few points across the film thickness and across the journal and the bearing in the radial direction, the temperature profile is very well predicted.

Vijayaraghavan, D.↗
NASA NTRS1990

Filtering analysis of a direct numerical simulation of the turbulent Rayleigh-Benard problem

A filtering analysis of a turbulent flow was developed which provides details of the path of the kinetic energy of the flow from its creation via thermal production to its dissipation. A low-pass spatial filter is used to split the velocity and the temperature field into a filtered component (composed mainly of scales larger than a specific size, nominally the filter width) and a fluctuation component (scales smaller than a specific size). Variables derived from these fields can fall into one of the above two ranges or be composed of a mixture of scales dominated by scales near the specific size. The filter is used to split the kinetic energy equation into three equations corresponding to the three scale ranges described above. The data from a direct simulation of the Rayleigh-Benard problem for conditions where the flow is turbulent are used to calculate the individual terms in the three kinetic energy equations. This is done for a range of filter widths. These results are used to study the spatial location and the scale range of the thermal energy production, the cascading of kinetic energy, the diffusion of kinetic energy, and the energy dissipation. These results are used to evaluate two subgrid models typically used in large-eddy simulations of turbulence. Subgrid models attempt to model the energy below the filter width that is removed by a low-pass filter.

Eidson, T. M.↗
DOE OSTI2024

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

Here, we present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface-based controlled-source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three-component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical-mesh-based electromagnetic modelling code and 3D tetrahedral-mesh-based electromagnetic modelling code. The former is particularly well-suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel-cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite-element electromagnetic code. This multi-step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel-fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple-source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep-localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

58 GEOSCIENCES↗
NASA NTRS2004

Superposed Epoch Analysis of Ring Current Geoeffectiveness Related to Solar Wind and Plasma Sheet Drivers

The goal of that proposal was to examine the relationship between solar wind drivers and ring current dynamics through data analysis and numerical simulations. The data analysis study was a statistical examination (via superposed epoch analyses) of a solar cycle's worth of storm data. Solar wind data, geophysical indices, and geosynchronous plasma data were collected for every time period with Dst< -50 nT from 1989 through 2002, and the storm list now exceeds 400 entries. This work was first conducted by a summer undergraduate student, Mr. John Vann (University of Kansas), with funding from the NSF Research Experience for Undergraduates program. It was then continued by a University of Michigan graduate student, Mr. Jichun Zhang. Mr. Zhang is now in his fourth year at U-M and is progressing very well toward a PhD in space science. His dissertation will be based on his data analysis and modeling efforts using this geomagnetic storm database. The results of the data analysis study have been the focus of several conference presentations, and the first manuscript has just been published. Two additional papers are presently being prepared, one on average (superposed) solar wind features for various storm subsets (e.g., intense storms at solar maximum), and another on geosynchronous plasma features for these same storm subsets. The latter result was highlighted by the TR&T program director in his presentation at the COSPAR meeting this summer.

Liemohm, M. W.↗
NASA NTRS1979

Sensitivity of Forecast Skill to Different Objective Analysis Schemes

Numerical weather forecasts are characterized by rapidly declining skill in the first 48 to 72 h. Recent estimates of the sources of forecast error indicate that the inaccurate specification of the initial conditions contributes substantially to this error. The sensitivity of the forecast skill to the initial conditions is examined by comparing a set of real-data experiments whose initial data were obtained with two different analysis schemes. Results are presented to emphasize the importance of the objective analysis techniques used in the assimilation of observational data.

Baker, W. E.↗
NASA NTRS1995

Evaluation of numerical schemes for the analysis of sound generation by blade-gust interaction

In this investigation three different numerical algorithms have been utilized to compute the flow about a flat plate in the presence of a transverse gust described by a sinusoidal disturbance. The three schemes include the MacCormack explicit finite difference scheme which is second order accurate in both time and space, the Gottlieb and Turkel modification of MacCormack's scheme which is fourth order accurate in space and second order accurate in time, (referred to as the 2-4 scheme), and a two step scheme developed by Bayliss et al. which has second order temporal accuracy and sixth order spatial accuracy (a 2-6 scheme). The flow field results are obtained with these schemes by using the same code with the only difference being the implementation of the respective solution algorithms. The problem is set up so that the sinusoidal disturbance is imposed at the surface of the flat plate as a surface boundary condition. Thus the problem is treated as scattering problem. The computed results include the time average of the acoustic pressure squared along grid lines five points away from the boundaries; distribution throughout the computational domain is monitored at various times. The numerical results are compared with an exact solution obtained by Atassi, Dusey, and Davis.

Scott, J. N.↗
NASA NTRS1988

Analysis Of Flow Across Cylinders

Report reviews various experimental, theoretical, and numerical methods of analysis of interactions between fluid and one or more cylinders across which it flows. Begins with discussion of mechanisms inducing vibrations in nests of cylinders. Major portion of review devoted to methods for numerical analysis of flow/cylinder interactions.

Kim, Sang-Wook↗
NASA NTRS1977

Computation of unsteady transonic flows through rotating and stationary cascades. 1: Method of analysis

A numerical method of solution of the inviscid, compressible, two-dimensional unsteady flow on a blade-to-blade stream surface through a stage (rotor and stator) or a single blade row of an axial flow compressor or fan is described. A cyclic procedure has been developed for representation of adjacent blade-to-blade passages which asymptotically achieves the correct phase between all passages of a stage. A shock-capturing finite difference method is employed in the interior of the passage, and a method of characteristics technique is used at the boundaries. The blade slipstreams form two of the passage boundaries and are treated as moving contact surfaces capable of supporting jumps in entropy and tangential velocity. The Kutta condition is imposed by requiring the slipstreams to originate at the trailing edges, which are assumed to be sharp. Results are presented for several transonic fan rotors and compared with available experimental data, consisting of holographic observations of shock structure and pressure contour maps. A subcritical stator solution is also compared with results from a relaxation method. Finally, a periodic solution for a stage consisting of 44 rotor blades and 46 stator blades is discussed.

Erdos, J. I.↗
NASA NTRS1986

An LU implicity scheme for high speed inlet analysis

A numerical method is developed to analyze the inviscid flowfield of a high speed inlet by the solution of the Euler equations. The lower-upper implicit scheme in conjunction with adaptive dissipation proves to be an efficient and robust nonoscillatory shock capturing technique for high Mach number flows as well as for transonic flows.

Yoon, S.↗
NASA NTRS1986

An LU implicit scheme for high speed inlet analysis

A numerical method is developed to analyze the inviscid flowfield of a high speed inlet by the solution of the Euler equations. The LU implicit scheme in conjunction with adaptive dissipation proves to be an efficient and robust nonoscillatory shock capturing technique for high Mach number flows as well as for transonic flows.

Yoon, S.↗
NASA NTRS1988

Liquid rocket boosters for the Space Shuttle - An outline of trade issues and analysis techniques

Numerous liquid rocket booster configurations were studied to replace the Solid Rocket Boosters on the current Space Shuttle. While improved safety, reliability, performance and economy of operation were sought, the designs also sought to minimize volumetric growth to reduce aerodynamic loads and structural modifications. Aero loads could be significant as a function of booster diameter if extrapolated directly from solid rocket characteristics. Procedures to meet these goals with analytical booster sizing and simple dynamic simulations are shown. Engine chamber pressure proved to be a good independent variable for implications of pressure-fed and pump-fed engines, fuel and material selection. For some configurations, ascent trajectory modifications were suggested: longer vertical rise, lofting, heads-up orientation, and higher first and second stage target altitudes. Since engine manufacture would dominate booster costs, booster simplification of reuse are discussed based on cost estimates.

Kelly, W. D.↗
NASA NTRS1989

Turbine rotor/stator flowfield analysis

A numerical study of the unsteady aerodynamic and thermal environment associated with axial turbine stages is presented. Computations are performed using a modification of the ROTOR1 rotor/stator interaction code. Two different turbine states are analyzed: the first state of the United Technologies Research Center large scale rotating rig and the first state of the Space Shuttle main engine (SSME) high pressure fuel turbopump. Time-averaged blade midspan pressure and heat transfer profiles are calculated using the following different surface boundary conditions: adiabatic wall, prescribed wall temperature, and prescribed heat flux. Numerical solutions for the large scale rotating rig are compared with experimental data. Unsteady pressure envelopes are also presented for each geometry. In addition, instantaneous contours are plotted for the SSME configuration.

Griffin, Lisa↗
NASA NTRS1974

Numerical methods for the design and analysis of wings at supersonic speeds

Numerical methods for the design and analysis of arbitrary-planform wings at supersonic speeds are reviewed. Certain deficiencies are revealed, particularly in application to wings with slightly subsonic leading edges. Recently devised numerical techniques which overcome the major part of these deficiencies are presented. The original development as well as the more recent revisions are subjected to a thorough review.

Carlson, H. W.↗