Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grid 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 1,279 records · Page 71

Finite element approximations in transient analysis

Alternative techniques for finite element models for transient analyses are explored, using as an example a pinned-pinned beam. Attention is given to the relationships between desired accuracy and the available computer resources. Stiffness, mass, damping, and loading components approximations are examined with note taken of optimization methods. Generation of regular grid results indicates a linear relationship between the number of elements or degrees of freedom and the number of eigenvalues for a given accuracy. Comparisons of the finite element models in terms of accuracy/element or accuracy/degree of freedom were found to be inaccurate for determining model efficiency. Higher accuracies were available from lumped mass and damping models than assumptions of consistent loading, which was preferred to statically equivalent loading.

Melosh, R. J.↗

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗

Analysis of an infinite array of rectangular microstrip patches with idealized probe feeds

A solution is presented to the problem of an infinite array of microstrip patches fed by idealized current probes. The input reflection coefficient is calculated versus scan angle in an arbitrary scan plane, and the effects of substrate parameters and grid spacing are considered. It is pointed out that even when a Galerkin method is used the impedance matrix is not symmetric due to phasing through a unit cell, as required for scanning. The mechanism by which scan blindness can occur is discussed. Measurement results are presented for the reflection coefficient magnitude variation with angle for E-plane, H-plane, and D-plane scans, for various substrate parameters. Measured results from waveguide simulators are also presented, and the scan blindness phenomenon is observed and discussed in terms of forced surface waves and a modified grating lobe diagram.

Pozar, D. M.↗

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗

Machine-Learning-Based Adaptive Thinning of CrIS Radiances to Improve Global Tropical Cyclone Analysis and Forecasts

This work is focused on optimizing the assimilation of hyperspectral infrared (IR) radiances from the Cross-track Infrared Sounder (CrIS) with the goal of improving the representation of tropical cyclones (TCs) in global analyses and forecasts. Current operational assimilation systems rely on subsampling IR radiances on a regular thinning grid. A new and improved adaptive methodology based on machine learning (ML) recognizes TCs from geostationary satellite imagery and is implemented in the Goddard Earth Observing System (GEOS) model and data assimilation framework. The ML methodology is extensively trained on existing TC data sets and creates for each TC a dynamic mask, based on the evolving shape and life cycle of that specific event. Once a TC mask is created, a switch is then activated in the data assimilation system to alter the thinning, ingesting more CrIS radiances within the moving mask, thus increasing the TC sampling. After the TC dissipates, the assimilation of CrIS radiances reverts to normal data density. Results of TC segmentation provided by a state-of-the-art generative machine learning model known as the Denoising Diffusion Probabilistic Model (DDPM) are compared to the previously used U-Net model. The new approach surpasses the performance of the previously developed one. The methodology is applied to both clear-sky and cloud-cleared radiances. Benefits from the latter methodology, particularly in improving the structure of TCs and the intensity forecasts, are presented.

Oreste Reale↗

Status and Mission Applicability of NASA's In-Space Propulsion Technology Project

The In-Space Propulsion Technology (ISPT) project develops propulsion technologies that will enable or enhance NASA robotic science missions. Since 2001, the ISPT project developed and delivered products to assist technology infusion and quantify mission applicability and benefits through mission analysis and tools. These in-space propulsion technologies are applicable, and potentially enabling for flagship destinations currently under evaluation, as well as having broad applicability to future Discovery and New Frontiers mission solicitations. This paper provides status of the technology development, near-term mission benefits, applicability, and availability of in-space propulsion technologies in the areas of advanced chemical thrusters, electric propulsion, aerocapture, and systems analysis tools. The current chemical propulsion investment is on the high-temperature Advanced Material Bipropellant Rocket (AMBR) engine providing higher performance for lower cost. Investments in electric propulsion technologies focused on completing NASA's Evolutionary Xenon Thruster (NEXT) ion propulsion system, a 0.6-7 kW throttle-able gridded ion system, and the High Voltage Hall Accelerator (HiVHAC) thruster, which is a mid-term product specifically designed for a low-cost electric propulsion option. Aerocapture investments developed a family of thermal protections system materials and structures; guidance, navigation, and control models of blunt-body rigid aeroshells; atmospheric models for Earth, Titan, Mars and Venus; and models for aerothermal effects. In 2009 ISPT started the development of propulsion technologies that would enable future sample return missions. The paper describes the ISPT project's future focus on propulsion for sample return missions. The future technology development areas for ISPT is: Planetary Ascent Vehicles (PAV), with a Mars Ascent Vehicle (MAV) being the initial development focus; multi-mission technologies for Earth Entry Vehicles (MMEEV) needed for sample return missions from many different destinations; propulsion for Earth Return Vehicles (ERV), transfer stages to the destination, and Electric Propulsion for sample return and low cost missions; and Systems/Mission Analysis focused on sample return propulsion. The ISPT project is funded by NASA's Science Mission Directorate (SMD).

Anderson, David J.↗

Using RotCFD to Predict Isolated XV-15 Rotor Performance

Experimental techniques to measure rotorcraft aerodynamic performance are widely used. However, the need exists to understand the limitations of ground based testing by augmenting the analysis of experimental test results with Computational Fluid Dynamics (CFD) modeling. The immediate objective of the present research is to develop an XV-15 Tilt Rotor Research Aircraft rotor model for investigation of wind tunnel wall interference. The predicted performance of the XV-15 during various flight modes is compared to theoretical and experimental data. This research is performed to support wind tunnel tests scheduled for 2016. A mid-fidelity RANS solver, RotCFD, is used with an unsteady, incompressible flow model and a realizable k- turbulence model. The rotor is modeled using an actuator disk model or blade element model with a momentum source approach. In RotCFD the setup, grid generation and running of cases is faster than many CFD codes which makes it a useful engineering tool. Performance predictions need not be as accurate as high-fidelity CFD codes, as long as wall effects can be properly simulated. Being able to accurately predict unsteady rotorcraft performance on desktop-class computers provides a quicker analysis of highly complex flows during the initial design phase.

CFD↗

Seasonal Variations of the Earth's Gravitational Field: An Analysis of Atmospheric Pressure, Ocean Tidal, and Surface Water Excitation

Monthly mean gravitational field parameters (denoted here as C(sub even)) that represent linear combinations of the primarily even degree zonal spherical harmonic coefficients of the Earth's gravitational field have been recovered using LAGEOS I data and are compared with those derived from gridded global surface pressure data of the National meteorological center (NMC) spanning 1983-1992. The effect of equilibrium ocean tides and surface water variations are also considered. Atmospheric pressure and surface water fluctuations are shown to be the dominant cause of observed annual C(sub even) variations. Closure with observations is seen at the 1sigma level when atmospheric pressure, ocean tide and surface water effects are include. Equilibrium ocean tides are shown to be the main source of excitation at the semiannual period with closure at the 1sigma level seen when both atmospheric pressure and ocean tide effects are included. The inverted barometer (IB) case is shown to give the best agreement with the observation series. The potential of the observed C(sub even) variations for monitoring mass variations in the polar regions of the Earth and the effect of the land-ocean mask in the IB calculation are discussed.

Dong, D,↗

Application of a Two-dimensional Unsteady Viscous Analysis Code to a Supersonic Throughflow Fan Stage

The Rai ROTOR1 code for two-dimensional, unsteady viscous flow analysis was applied to a supersonic throughflow fan stage design. The axial Mach number for this fan design increases from 2.0 at the inlet to 2.9 at the outlet. The Rai code uses overlapped O- and H-grids that are appropriately packed. The Rai code was run on a Cray XMP computer; then data postprocessing and graphics were performed to obtain detailed insight into the stage flow. The large rotor wakes uniformly traversed the rotor-stator interface and dispersed as they passed through the stator passage. Only weak blade shock losses were computerd, which supports the design goals. High viscous effects caused large blade wakes and a low fan efficiency. Rai code flow predictions were essentially steady for the rotor, and they compared well with Chima rotor viscous code predictions based on a C-grid of similar density.

Steinke, Ronald J.↗

Euler and Navier-Stokes solutions for hypersonic flows

An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.

Thareja, Rajiv R.↗

Modelling a single phase voltage controlled rectifier using Laplace transforms

The development of a 20 kHz, AC power system by NASA for large space projects has spurred a need to develop models for the equipment which will be used on these single phase systems. To date, models for the AC source (i.e., inverters) have been developed. It is the intent of this paper to develop a method to model the single phase voltage controlled rectifiers which will be attached to the AC power grid as an interface for connected loads. A modified version of EPRI's HARMFLO program is used as the shell for these models. The results obtained from the model developed in this paper are quite adequate for the analysis of problems such as voltage resonance. The unique technique presented in this paper uses the Laplace transforms to determine the harmonic content of the load current of the rectifier rather than a curve fitting technique. Laplace transforms yield the coefficient of the differential equations which model the line current to the rectifier directly.

Kraft, L. Alan↗

Effects of habitat map generalization in biodiversity assessment

Species richness is being mapped as part of an inventory of biological diversity in California (i.e., gap analysis). Species distributions are modeled with a GIS on the basis of maps of each species' preferred habitats. Species richness is then tallied in equal-area sampling units. A GIS sensitivity analysis examined the effects of the level of generalization of the habitat map on the predicted distribution of species richness in the southern Sierra Nevada. As the habitat map was generalized, the number of habitat types mapped within grid cells tended to decrease with a corresponding decline in numbers of species predicted. Further, the ranking of grid cells in order of predicted numbers of species changed dramatically between levels of generalization. Areas predicted to be of greatest conservation value on the basis of species richness may therefore be sensitive to GIS data resolution.

Stoms, David M.↗

A 3D-CFD code for accurate prediction of fluid flows and fluid forces in seals

Current and future turbomachinery requires advanced seal configurations to control leakage, inhibit mixing of incompatible fluids and to control the rotodynamic response. In recognition of a deficiency in the existing predictive methodology for seals, a seven year effort was established in 1990 by NASA's Office of Aeronautics Exploration and Technology, under the Earth-to-Orbit Propulsion program, to develop validated Computational Fluid Dynamics (CFD) concepts, codes and analyses for seals. The effort will provide NASA and the U.S. Aerospace Industry with advanced CFD scientific codes and industrial codes for analyzing and designing turbomachinery seals. An advanced 3D CFD cylindrical seal code has been developed, incorporating state-of-the-art computational methodology for flow analysis in straight, tapered and stepped seals. Relevant computational features of the code include: stationary/rotating coordinates, cylindrical and general Body Fitted Coordinates (BFC) systems, high order differencing schemes, colocated variable arrangement, advanced turbulence models, incompressible/compressible flows, and moving grids. This paper presents the current status of code development, code demonstration for predicting rotordynamic coefficients, numerical parametric study of entrance loss coefficients for generic annular seals, and plans for code extensions to labyrinth, damping, and other seal configurations.

Athavale, M. M.↗

Oceanic Excitation of Polar Motion During 1992-1994

Two global ocean general circulation models have been used to compute the angular momentum of nontidal oceanic current and sea level height variations: (1) the Princeton Modular Ocean Model (MOM) having 22 vertical layers and a rigid lid, and (2) the Miami Isopycnic-Coordinate Ocean Model (MICOM) having 11 vertical layers with a mixed layer and a free surface. Both models were run on the same 2 degree longitude by 1 degree latitude grid spanning 80 S to 80 N latitude. Following a 10-year spin-up with climatological air-sea fluxes, both models were forced during 1992-1994 with daily wind and heat flux from the NCEP operational analysis and sea surface salinity restoring to Levitus climatology. After correcting for the effects of mass non-conservation in these models, they are shown to predict similar effects on polar motion excitation of sea level height variations, but the effect of currents predicted by MOM is much larger than that predicted by MICOM.

Gross, R.↗

Modeling of Supersonic Film Cooling on the J-2X Nozzle Extension

Supersonic film cooling (SSFC) of nozzles has been used in several liquid rocket engine designs, and is being applied to the nozzle extension (NE) of the J-2X upper stage engine currently under development. Turbine exhaust gas (TEG) is injected tangentially from a manifold along the NE, and provides a thermal barrier from the core nozzle flow for the NE. As the TEG stream mixes with the nozzle flow, the effectiveness of the thermal barrier is reduced. This paper documents computational fluid dynamics (CFD) analysis work performed by NASA Marshall Space Flight Center (MSFC) to model the flow of the TEG through the manifold, into the nozzle, and the subsequent mixing of the TEG stream with the core flow. The geometry and grid of the TEG manifold, structural support ribs, and the NE wall will be shown, and the CFD boundary conditions described. The Loci-CHEM CFD code used in this work will also be briefly described. A unique approach to modeling the combined TEG manifold/thrust chamber assembly (TCA) was employed, as it was not practical to model the entire 360 circumferential range in one simulation. Prior CFD validation work modeling Calspan SSFC experiments in the early 1990s, documented in a previous AIAA paper, will also be briefly discussed. The fluid dynamics of the TEG flow through the manifold, into and between the structural support ribs, and into the nozzlette that feeds the TCA will be described. Significant swirl and non-uniformities are present, which along with the wakes from the ribs, act to degrade the film cooling effectiveness compared to idealized injection of TEG gas. The effect of these flow characteristics on the adiabatic wall temperature profile on the NE will be discussed.

Ruf, Joseph H.↗

Discharge Chamber Performance of the NEXIS Ion Thruster

The Nuclear-Electric Xenon Ion System (NEXIS) thruster was designed to produce greater than or equal to 70% efficiency at ISPs in excess of 6500 sec and total power levels in excess of 15 kW. In order to achieve this performance, the thruster requires a large area plasma generator capable of high propellant utilimtion efficiency and low discharge loss while producing a very flat, uniform beam profile. Fortunately, larger thrusters can be made more uniform and efficient due to the higher volume to surface ratio, provided that the magnetic cusp confinement is designed properly and the thruster length to diameter ratio is adequate. This paper describes the discharge chamber performance of the NEXIS Laboratory Model (LM) thruster. The LM discharge chamber is 65 cm in diameter at the grid plane and uses 6 ring-cusps to provide magnetic confinement of the plasma. The thruster was tested with flat carbon-carbon composite grids with the hole pattern masked to 57 cm in diameter and a conventional Type-B "1/2" diameter hollow cathode. During the preliminary "discharge only" tests, the LM thruster demonstrated profile factors of 0.84 and a discharge loss of about 160 eV/ion at 25 V discharge voltage and over 90% propellant utilization efficiency in simulated beam extraction experiments at 3.9 A of beam current. Analysis of the data from these tests used the discharge-only model developed by Brophy. Subsequent beam extraction experiments validated the key variables used in the model to predict the performance from the discharge-only data, and demonstrated 3.9 A of beam current at over 90% propellant utilization efficiency with a flatness parameter of better than 0.8 and a discharge loss of about 185 eV/ion. The slightly higher discharge loss measured during beam extractions was found to be due to a lower screen transparency in the as-manufactured LM grid set. Plasma measurements with a scanning probe internal to the thruster near the screen grid showed plasma densities over l x 10(exp 11) per cubic centimeter and electron temperatures of 3.5 to 5.5 eV depending on the operation parameters. The performance of the NEXIS discharge chamber contributed to the over 78% thruster efficiency measured during beam extraction at 7500 sec ISP and 25 kW of power, and over 81% thruster efficiency measured at 8500 sec ISP.

plasma generation↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

An objective analysis technique for extrapolating tidal fields

An interpolation technique which allows accurate extrapolation of tidal height fields in the ocean basins by making use of selected satellite altimetry measurements and/or conventional gauge measurements was developed and tested. A normal mode solution for the Atlantic and Indian Oceans was obtained by means of a finite difference grid. Normal mode amplitude maps are presented.

Sanchez, B. V.↗