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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

An experimental and analytical investigation of a hypersonic cruise vehicle at Mach 6

Flow field properties over a hypersonic vehicle are measured and compared with results from the numerical analysis of inviscid supersonic flow theory. Good agreement between theory and experiment at zero angle of attack, and significant discrepancy at angles of attack greater than or equal to 5 deg are found. The experiments indicate that turbulent boundary layer without flow separation exists over the model surface at zero angle of attack. Vortex lift-off, flow separation, and shock boundary layer interaction occur over the leeward surface, and near the wing-fuselage junction with angles of attack greater than or equal to 5 deg. Improvements of the existing numerical method to compute hypersonic cruise vehicle flow field at large angle of attack are presented.

Wang, C. R.↗
NASA NTRS2011

Consequences of the Breakout Model for Particle Acceleration in CMEs and Flares

The largest and most efficient particle accelerators in the solar system are the giant events consisting of a fast coronal mass ejection (CME) and an intense X-class solar flare. Both flares and CMEs can produce l0(exp 32) ergs or more in nonthermal particles. Two general processes are believed to be responsible: particle acceleration at the strong shock ahead of the CME, and reconnection-driven acceleration in the flare current sheet. Although shock acceleration is relatively well understood, the mechanism by which flare reconnection produces nonthermal particles is still an issue of great debate. We address the question of CME/flare particle acceleration in the context of the breakout model using 2.5D MHD simulations with adaptive mesh refinement (AMR). The AMR capability allows us to achieve ultra-high numerical resolution and, thereby, determine the detailed structure and dynamics of the flare reconnection region. Furthermore, we employ newly developed numerical analysis tools for identifying and characterizing magnetic nulls, so that we can quantify accurately the number and location of magnetic islands during reconnection. Our calculations show that flare reconnection is dominated by the formation of magnetic islands. In agreement with many other studies, we find that the number of islands scales with the effective Lundquist number. This result supports the recent work by Drake and co-workers that postulates particle acceleration by magnetic islands. On the other hand, our calculations also show that the flare reconnection region is populated by numerous shocks and other indicators of strong turbulence, which can also accelerate particles. We discuss the implications of our calculations for the flare particle acceleration mechanism and for observational tests of the models.

Antiochos, S. K.↗
NASA NTRS1978

Nonlinear axisymmetric liquid currents in spherical annuli

A numerical analysis of non-linear axisymmetric viscous flows in spherical annuli of different gap sizes is presented. Only inner sphere was supposed to rotate at a constant angular velocity. The streamlines, lines of constant angular velocity, kinetic energy spectra, and spectra of velocity components are obtained. A total kinetic energy and torque needed to rotate the inner sphere are calculated as functions of Re for different gap sizes. In small-gap annulus nonuniqueness of steady solutions of Navier-Stokes equations is established and regions of different flow regime existences are found. Numerical solutions in a wide-gap annulus and experimental results are used in conclusions about flow stability in the considered range of Re. The comparison of experimental and numerical results shows close qualitative and quantitative agreement.

Astafyeva, N. M.↗
NASA NTRS1988

Multi-airfoil Navier-Stokes simulations of turbine rotor-stator interaction

An accurate numerical analysis of the flows associated with rotor-stator configurations in turbomachinery can be extremely helpful in optimizing the performance of turbomachinery. In this study the unsteady, thin-layer, Navier-Stokes equations in two spatial dimensions are solved on a system of patched and overlaid grids for a rotor-stator configuration from an axial turbine. The governing equations are solved using a finite-difference, upwind algorithm that is set in an iterative, implicit framework. Results in the form of pressure contours, time-averaged pressures, unsteady pressure amplitudes and phase are presented. The numerical results are compared with experimental data and the agreement is found to be good. The results are also compared with those of an earlier study which used only one rotor and one stator. The current study uses multiple rotors and stators and a pitch ratio that is much closer to the experimental ratio. Consequently the results of this study are found to be closer to the experimental data.

Rai, Man Mohan↗
NASA NTRS1985

Navier-Stokes simulations of rotor-stator interaction using patched and overlaid grids

An accurate numerical analysis of the flows associated with rotor-stator configurations in turbomachinery can be very helpful in optimizing the performance of turbomachinery. In this study the unsteady, thin-layer, Navier-Stokes equations are solved using a system of patched and overlaid grids for a rotor-stator configuration of an axial turbine. The equations necessary for an accurate transfer of information between the several grids are briefly described within the framework of an iterative, implicit algorithm. Results in the form of Mach number contours, time-averaged pressures, unsteady pressures, amplitudes, and phase are presented. The numerical results are also compared with experimental data and the agreement is found to be good.

Rai, M. M.↗
NASA NTRS1981

Impulsive phase of solar flares. I - Characteristics of high energy electrons

The models describing the operative physical processes occurring during the impulsive phase of a solar flare can be broadly classified as thermal or nonthermal models. The current investigation is part of a series of studies which analyze the general characteristics of the 'nonthermal' models. The results are presented of an analysis regarding the steady state distribution of energetic electrons within the flare plasma. The time independent Fokker-Planck equation is used to determine the steady state distribution of the electrons in the ambient plasma. The model investigated encompasses thick-target, thin-target, and trap aspects of nonthermal hard X-ray models. High energy electrons with some initial energy spectrum and pitch angle distribution are injected at a height above the photosphere either into a closed magnetic loop structure or into an open field structure. Some approximate analytic solutions are discussed and the results of a numerical analysis are presented.

Leach, J.↗
NASA NTRS2003

Analysis of ETMS Data Quality for Traffic Flow Management Decisions

The data needed for air traffic flow management decision support tools is provided by the Enhanced Traffic Management System (ETMS). This includes both the tools that are in current use and the ones being developed for future deployment. Since the quality of decision support provided by all these tools will be influenced by the quality of the input ETMS data, an assessment of ETMS data quality is needed. Motivated by this desire, ETMS data quality is examined in this paper in terms of the unavailability of flight plans, deviation from the filed flight plans, departure delays, altitude errors and track data drops. Although many of these data quality issues are not new, little is known about their extent. A goal of this paper is to document the magnitude of data quality issues supported by numerical analysis of ETMS data. Guided by this goal, ETMS data for a 24-hour period were processed to determine the number of aircraft with missing flight plan messages at any given instant of time. Results are presented for aircraft above 18,000 feet altitude and also at all altitudes. Since deviation from filed flight plan is also a major cause of trajectory-modeling errors, statistics of deviations are presented. Errors in proposed departure times and ETMS-generated vertical profiles are also shown. A method for conditioning the vertical profiles for improving demand prediction accuracy is described. Graphs of actual sector counts obtained using these vertical profiles are compared with those obtained using the Host data for sectors in the Fort Worth Center to demonstrate the benefit of preprocessing. Finally, results are presented to quantify the extent of data drops. A method for propagating track positions during ETMS data drops is also described.

Chatterji, Gano B.↗
NASA NTRS1974

Calculation of turbulent shear stress in supersonic boundary layer flows

An analysis of turbulent boundary layer flow characteristics and the computational procedure used are discussed. The integrated mass and momentum flux profiles and differentials of the integral quantities are used in the computations so that local evaluation of the streamwise velocity gradient is not necessary. The computed results are compared with measured shear stress data obtained by using hot wire anemometer and laser velocimeter techniques. The flow measurements were made upstream and downstream of an adiabatic unseparated interaction of an oblique shock wave with the turbulent boundary layer on the flat wall of a two dimensional wind tunnel. A comparison of the numerical analysis and actual measurements is made and the effects of small differences in mean flow profiles on the computed shear stress distributions are discussed.

Sun, C. C.↗
NASA NTRS1992

Validating a large geophysical data set: Experiences with satellite-derived cloud parameters

We are validating the global cloud parameters derived from the satellite-borne HIRS2 and MSU atmospheric sounding instrument measurements, and are using the analysis of these data as one prototype for studying large geophysical data sets in general. The HIRS2/MSU data set contains a total of 40 physical parameters, filling 25 MB/day; raw HIRS2/MSU data are available for a period exceeding 10 years. Validation involves developing a quantitative sense for the physical meaning of the derived parameters over the range of environmental conditions sampled. This is accomplished by comparing the spatial and temporal distributions of the derived quantities with similar measurements made using other techniques, and with model results. The data handling needed for this work is possible only with the help of a suite of interactive graphical and numerical analysis tools. Level 3 (gridded) data is the common form in which large data sets of this type are distributed for scientific analysis. We find that Level 3 data is inadequate for the data comparisons required for validation. Level 2 data (individual measurements in geophysical units) is needed. A sampling problem arises when individual measurements, which are not uniformly distributed in space or time, are used for the comparisons. Standard 'interpolation' methods involve fitting the measurements for each data set to surfaces, which are then compared. We are experimenting with formal criteria for selecting geographical regions, based upon the spatial frequency and variability of measurements, that allow us to quantify the uncertainty due to sampling. As part of this project, we are also dealing with ways to keep track of constraints placed on the output by assumptions made in the computer code. The need to work with Level 2 data introduces a number of other data handling issues, such as accessing data files across machine types, meeting large data storage requirements, accessing other validated data sets, processing speed and throughput for interactive graphical work, and problems relating to graphical interfaces.

Kahn, Ralph↗
NASA NTRS1972

An investigation to improve the Menhaden fishery prediction and detection model through the application of ERTS-A data

The author has identified the following significant results. Preliminary analyses indicate that several important relationships have been observed utilizing ERTS-1 imagery. Of most significance is that in the Mississippi Sound, as elsewhere, considerable detail exists as to turbidity patterns in the water column. Simple analysis is complicated by the apparent interaction between actual turbidity, turbidity induced by shoal water, and actual imaging of the bottom in extreme shoal water. A statistical approach is being explored which shows promise of at least partially separating these effects so that partitioning of true turbid plumes can be accomplished. This partitioning is of great importance to this program in that supportive data seem to indicate that menhaden occur more frequently in turbid areas. In this connection four individual captures have been associated with a major turbid feature imaged on 6 August. If a significant relationship between imaged turbid features and catch distribution can be established, for example by graphic and/or numeric analysis, it will represent a major advancement for short term prediction of commercially accessible menhaden.

Maughan, P. M.↗
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
NASA NTRSDate not supplied

Development of a C-ELS Specimen-Based Numerical Benchmark for Mode II Delamination and Assessment of Two VCCT-Based Propagation Strategies

A finite element (FE) benchmark example inspired by the calibrated end-loaded split (C-ELS) specimen is developed and used to assess the performance of delamination propagation capabilities based on linear elastic fracture mechanics (LEFM). The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the virtual crack closure technique (VCCT) delamination propagation capabilities in the Abaqus/Standard®1 FE code and the recently developed Progressive Release eXplicit-VCCT (PRX-VCCT) method are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination encourages stable delamination propagation. Detailed results from three-dimensional analyses with aligned and misaligned meshes and two levels of mesh refinement are provided. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur for propagation analyses using Abaqus/Standard VCCT. In comparison, continuous, smooth, delamination fronts are observed for PRX-VCCT. The use of the benchmark case to assess different VCCT-based propagation strategies illustrates the value of establishing benchmark cases.

Composite Materials↗
NASA NTRS1986

Sediment-transport experiments in zero-gravity

One of the important parameters in the analysis of sediment entrainment and transport is gravitational attraction. The availability of a laboratory in Earth orbit would afford an opportunity to conduct experiments in zero and variable gravity environments. Elimination of gravitational attraction as a factor in such experiments would enable other critical parameters (such as particle cohesion and aerodynamic forces) to be evaluated much more accurately. A Carousel Wind Tunnel (CWT) is proposed for use in conducting experiments concerning sediment particle entrainment and transport in a space station. In order to test the concept of this wind tunnel design a one third scale model CWT was constructed and calibrated. Experiments were conducted in the prototype to determine the feasibility of studying various aeolian processes and the results were compared with various numerical analysis. Several types of experiments appear to be feasible utilizing the proposed apparatus.

Iversen, J. D.↗
NASA NTRS1987

Sediment-transport experiments in zero-gravity

One of the important parameters in the analysis of sediment entrainment and transport is gravitational attraction. The availability of a laboratory in earth orbit would afford an opportunity to conduct experiments in zero and variable gravity environments. Elimination of gravitational attraction as a factor in such experiments would enable other critical parameters (such as particle cohesion and aerodynamic forces) to be evaluated much more accurately. A Carousel Wind Tunnel (CWT) is proposed for use in conducting experiments concerning sediment particle entrainment and transport in a space station. In order to test the concept of this wind tunnel design a one third scale model CWT was constructed and calibrated. Experiments were conducted in the prototype to determine the feasibility of studying various aeolian processes and the results were compared with various numerical analysis. Several types of experiments appear to be feasible utilizing the proposed apparatus.

Iversen, James D.↗
NASA NTRS1976

On a numerical solution of the supersonic panel flutter eigenproblem

An automated digital computer procedure is presented in this paper which enables efficient solution of the eigenvalue problem associated with the supersonic panel flutter phenomena. The step-by-step incremental solution procedure is based on an inverse iteration technique which effectively utilizes solution results from the previous step in determining such results during the current solution step. Also, the computations are limited to the determination of a few specific roots only, which are expected to contain the flutter mode, and this is achieved at each step without having to compute any other root. The structural discretization achieved by the finite-element method yields highly banded stiffness, mass, and aerodynamic matrices; the aerodynamic matrix evaluated by the linearized piston theory is real but unsymmetric in nature. The solution algorithm presented in this paper fully exploits the banded form of the associated matrices, and the resulting computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be most efficient and economical when compared to existing procedures of such analysis. Numerical results are presented for a two-dimensional panel flutter problem.

Gupta, K. K.↗
DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗
DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗
NASA NTRS1999

Computation of Steady and Unsteady Laminar Flames: Theory

In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Hagstrom, Thomas↗