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 145 records · Page 8

Numerical analysis of the three-dimensional flow in the main injector assembly of the Space Shuttle Main Engine

The Main Injector Assembly of the Space Shuttle Main Engine supplies the propellants to the Main Combustion Chamber through a large number of vertical injector elements. The gas flow around these elements (LOX posts) is three-dimensional, turbulent and compressible. This paper presents the results of numerical modeling of the hot and cool sections of the Main Injector Assembly and shows that the shields on the outermost row of injector elements affect the flow and the pressure drop in the hot section significantly. It is suggested that appropriate experiments are undertaken to verify the predictions.

Mukerjee, T.↗

A stochastic regulator for integrated communication and control systems. I - Formulation of control law. II - Numerical analysis and simulation

A state feedback control law for integrated communication and control systems (ICCS) is formulated by using the dynamic programming and optimality principle on a finite-time horizon. The control law is derived on the basis of a stochastic model of the plant which is augmented in state space to allow for the effects of randomly varying delays in the feedback loop. A numerical procedure for synthesizing the control parameters is then presented, and the performance of the control law is evaluated by simulating the flight dynamics model of an advanced aircraft. Finally, recommendations for future work are made.

Liou, Luen-Woei↗

Computer-aided quantization and numerical analysis of superconducting circuits

The development of new superconducting circuits and the improvement of existing ones rely on the accurate modeling of spectral properties which are key to achieving the needed advances in qubit performance. Systematic circuit analysis at the lumped-element level, starting from a circuit network and culminating in a Hamiltonian appropriately describing the quantum properties of the circuit, is a well-established procedure, yet cumbersome to carry out manually for larger circuits. We present work utilizing symbolic computer algebra and numerical diagonalization routines versatile enough to tackle a variety of circuits. Results from this work are accessible through a newly released module of the scqubits package.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A numerical analysis of the British Experimental Rotor Program blade

Two Computational Fluid Dynamic codes which solve the compressible full-potential and the Reynolds-Averaged Thin-Layer Navier-Stokes equations were used to analyze the nonrotating aerodynamic characteristics of the British Experimental Rotor Program (BERP) helicopter blade at three flow regimes: low angle of attack, high angle of attack and transonic. Excellent agreement was found between the numerical results and experiment. In the low angle of attack regime, the BERP had less induced drag than a comparable aspect ratio rectangular planform wing. At high angle of attack, the blade attained high-lift by maintaining attached flow at the outermost spanwise locations. In the transonic regime, the BERP design reduces the shock strength at the outer spanwise locations which affects wave drag and shock-induced separation. Overall, the BERP blade exhibited many favorable aerodynamic characteristics in comparison to conventional helicopter rotor blades.

Duque, Earl P. N.↗

Numerical analysis of supersonic flow through oscillating cascade sections by using a deforming grid

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small-perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small-perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a section of the rotor on a supersonic throughflow compressor designed at NASA Lewis Research Center. Preliminary results indicate that two-dimensional, flat plate analysis predicts conservative flutter boundaries.

Huff, Dennis L.↗

Numerical analysis of supersonic flow through oscillating cascade sections by using a deforming grid

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small-perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small-perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a section of the rotor on a supersonic throughflow compressor designed at NASA Lewis Research Center. Preliminary results indicate that two-dimensional, flat plate analysis predicts conservative flutter boundaries.

Huff, Dennis L.↗

Numerical Analysis of Film Cooling at High Blowing Ratio

Computational Fluid Dynamics is used in the analysis of a film cooling jet in crossflow. Predictions of film effectiveness are compared with experimental results for a circular jet at blowing ratios ranging from 0.5 to 2.0. Film effectiveness is a surface quantity which alone is insufficient in understanding the source and finding a remedy for shortcomings of the numerical model. Therefore, in addition, comparisons are made to flow field measurements of temperature along the jet centerline. These comparisons show that the CFD model is accurately predicting the extent and trajectory of the film cooling jet; however, there is a lack of agreement in the near-wall region downstream of the film hole. The effects of main stream turbulence conditions, boundary layer thickness, turbulence modeling, and numerical artificial dissipation are evaluated and found to have an insufficient impact in the wake region of separated films (i.e. cannot account for the discrepancy between measured and predicted centerline fluid temperatures). Analyses of low and moderate blowing ratio cases are carried out and results are in good agreement with data.

El-Gabry, Lamyaa↗

Numerical analysis of the asymptotic two-point boundary value solution for N-body trajectories.

Previously published asymptotic solutions for lunar and interplanetary trajectories have been modified and combined to formulate a general analytical boundary value solution applicable to a broad class of trajectory problems. In addition, the earlier first-order solutions have been extended to second-order to determine if improved accuracy is possible. Comparisons between the asymptotic solution and numerical integration for several lunar and interplanetary trajectories show that the asymptotic solution is generally quite accurate. Also, since no iterations are required, a solution to the boundary value problem is obtained in a fraction of the time required for numerically integrated solutions.

Lancaster, J. E.↗

Numerical Analysis of a Radiant Heat Flux Calibration System

A radiant heat flux gage calibration system exists in the Flight Loads Laboratory at NASA's Dryden Flight Research Center. This calibration system must be well understood if the heat flux gages calibrated in it are to provide useful data during radiant heating ground tests or flight tests of high speed aerospace vehicles. A part of the calibration system characterization process is to develop a numerical model of the flat plate heater element and heat flux gage, which will help identify errors due to convection, heater element erosion, and other factors. A 2-dimensional mathematical model of the gage-plate system has been developed to simulate the combined problem involving convection, radiation and mass loss by chemical reaction. A fourth order finite difference scheme is used to solve the steady state governing equations and determine the temperature distribution in the gage and plate, incident heat flux on the gage face, and flat plate erosion. Initial gage heat flux predictions from the model are found to be within 17% of experimental results.

Jiang, Shanjuan↗

Numerical Analysis of Radiation Effects on Fiber Optic Sensors

Optical fiber sensors (OFS) are a potential candidate for monitoring physical parameters in nuclear environments. However, under an irradiation field the optical response of the OFS is modified via three primary mechanisms: (i) radiation-induced attenuation (RIA), (ii) radiation-induced emission (RIE), and (iii) radiation-induced compaction (RIC). For resonance-based sensors, RIC plays a significant role in modifying their performance characteristics. In this paper, we numerically investigate independently the effects of RIC and RIA on three types of OFS widely considered for radiation environments: fiber Bragg grating (FBG), long-period grating (LPG), and Fabry-Perot (F-P) sensors. In our RIC modeling, experimentally calculated refractive index (RI) changes due to low-dose radiation are extrapolated using a power law to calculate density changes at high doses. The changes in RI and length are subsequently calculated using the Lorentz–Lorenz relation and an established empirical equation, respectively. The effects of both the change in the RI and length contraction on OFS are modeled for both low and high doses using FIMMWAVE, a commercially available vectorial mode solver. An in-depth understanding of how radiation affects OFS may reveal various potential OFS applications in several types of radiation environments, such as nuclear reactors or in space.

42 ENGINEERING↗

Numerical Analysis of Base Flowfield for a Four-Engine Clustered Nozzle Configuration

Excessive base heating has been a problem for many launch vehicles. For certain designs such as the direct dump of turbine exhaust inside and at the lip of the nozzle, the potential burning of the turbine exhaust in the base region can be of great concern. Accurate prediction of the base environment at altitudes is therefore very important during the vehicle design phase. Otherwise, undesirable consequences may occur. In this study, the turbulent base flowfield of a cold flow experimental investigation for a four-engine clustered nozzle was numerically benchmarked using a pressure-based computational fluid dynamics (CFD) method. This is a necessary step before the benchmarking of hot flow and combustion flow tests can be considered. Since the medium was unheated air, reasonable prediction of the base pressure distribution at high altitude was the main goal. Several physical phenomena pertaining to the multiengine clustered nozzle base flow physics were deduced from the analysis.

Wang, Ten-See↗

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗

Numerical Analysis of an Impinging Jet Reactor for the CVD and Gas-Phase Nucleation of Titania

We model a cold-wall atmospheric pressure impinging jet reactor to study the CVD and gas-phase nucleation of TiO2 from a titanium tetra-iso-propoxide (TTIP)/oxygen dilute source gas mixture in nitrogen. The mathematical model uses the computational code FIDAP and complements our recent asymptotic theory for high activation energy gas-phase reactions in thin chemically reacting sublayers. The numerical predictions highlight deviations from ideality in various regions inside the experimental reactor. Model predictions of deposition rates and the onset of gas-phase nucleation compare favorably with experiments. Although variable property effects on deposition rates are not significant (approximately 11 percent at 1000 K), the reduction rates due to Soret transport is substantial (approximately 75 percent at 1000 K).

Gokoglu, Suleyman A.↗

Earth zonal harmonics from rapid numerical analysis of long satellite arcs

A zonal geopotential is presented to degree 21 from evaluation of mean elements for 21 satellites including 2 of low inclination. Each satellite is represented by an arc of at least one apsidal rotation. The lengths range from 200 to 800 days. Differential correction of the initial elements in all of the arcs, together with radiation pressure and atmospheric drag coefficients, was accomplished simultaneously with the correction for the zonal harmonics. The satellite orbits and their variations are generated by numerical integration of the Lagrange equations for mean elements. Disturbances due to precession and nutation of the earth's pole, atmospheric drag, radiation pressure and luni-solar gravity are added at from 1- to 8-day intervals in the integrated orbits. The results agree well with recent solutions from other authors using different methods and different satellite sets.

Wagner, C. A.↗

Analytical and Numerical Analysis of Lithium Plating Onset in Single and Bilayer Graphite Electrodes during Fast Charging

Lithium plating is commonly observed in anodes charged at fast rates, and can lead to capacity loss and battery safety issues. The increased risk of plating has been attributed to transport limitations, and architectured electrodes may reduce plating risk. However, while theoretical studies have shown that reaction non-uniformity arises due to interplay of transport limitations, anode open circuit voltage behavior and reaction kinetics, its effect on lithium plating has not been studied. We use analytic and numerical simulations to predict onset of plating in graphite anode half-cells at high C-rates and demonstrate how anodes with layered porosities can delay plating. Simplified analytical models identify trends for plating onset and predictions are calibrated against numerical models. A calibrated numerical model of graphite demonstrates qualitative agreement with analytical model predictions. This reaction inhomogeneity mechanism occurs in the absence of lithium ion depletion, indicating that these mechanisms may contribute to capacity loss independently or simultaneously. A bilayer model of graphite exhibits delayed plating onset, and an optimization procedure is presented. This theoretical work presents quantitative and mechanistic insight on how reaction inhomogenity affects lithium metal plating onset and can be used as a guide to engineer anodes resistant to lithium plating.

Mijailovic, Aleksandar S.↗

Numerical Analysis of Incipient Separation on 53 Deg Swept Diamond Wing

A systematic analysis of incipient separation and subsequent vortex formation from moderately swept blunt leading edges is presented for a 53 deg swept diamond wing. This work contributes to a collective body of knowledge generated within the NATO/STO AVT-183 Task Group titled 'Reliable Prediction of Separated Flow Onset and Progression for Air and Sea Vehicles'. The objective is to extract insights from the experimentally measured and numerically computed flow fields that might enable turbulence experts to further improve their models for predicting swept blunt leading-edge flow separation. Details of vortex formation are inferred from numerical solutions after establishing a good correlation of the global flow field and surface pressure distributions between wind tunnel measurements and computed flow solutions. From this, significant and sometimes surprising insights into the nature of incipient separation and part-span vortex formation are derived from the wealth of information available in the computational solutions.

Frink, Neal T.↗

Hurricane Debby - Analysis and numerical forecasts using VAS soundings

The utility of VISSR Atmospheric Sounder (VAS) temperature and moisture soundings in defining the storm and its surroundings at subsynoptic scales has been examined using a numerical analysis and prognosis system. In particular, VAS temperature and moisture soundings and cloud and water vapor motion winds have been used in numerical analysis for three time periods. It is shown that the VAS temperature and moisture data which specify temperature and moisture well in cloud free regions are complemented by cloud and water vapor wind data which provide horizontal gradient information for the cloudy areas. The loss of analysis integrity due to the reduction of VAS data density in the cloudy regions associated with synoptic activity is ameliorated by using cloud and water vapor motion winds. The improvement in numerical forecasts resulting from the addition of these data to the numerical data base is also recorded.

Le Marshall, J. F.↗

A numerical analysis of high-temperature heat pipe startup from the frozen state

Continuum and rarefied vapor flows co-exist along the heat pipe length for most of the startup period. A two-region model is proposed in which the vapor flow in the continuum region is modeled by the compressible Navier-Stokes equations, and the vapor flow in the rarefied region is simulated by a self-diffusion model. The two vapor regions are linked with appropriate boundary conditions, and heat pipe wail, wick, and vapor flow are solved as a conjugate problem. The numerical solutions for the entire heat pipe startup process from the frozen state are compared with the corresponding experimental data with good agreement.

Cao, Y.↗