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 487 records · Page 27

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗

Active Control of F/A-18 Vertical Tail Buffeting using Piezoelectric Actuators

Vertical tail buffeting is a serious multidisciplinary problem that limits the performance of twin-tail fighter aircraft. The buffet problem occurs at high angles of attack when the vortical flow breaks down ahead of the vertical tails resulting in unsteady and unbalanced pressure loads on the vertical tails. This paper describes a multidisciplinary computational investigation for buffet load alleviation of full F/A-18 aircraft using distributed piezoelectric actuators. The inboard and outboard surfaces of the vertical tail are equipped with piezoelectric actuators to control the buffet responses in the first bending and torsion modes. The electrodynamics of the smart structure are expressed with a three-dimensional finite element model. A single-input-single-output controller is designed to drive the active piezoelectric actuators. High-fidelity multidisciplinary analysis modules for the fluid dynamics, structure dynamics, electrodynamics of the piezoelectric actuators, fluid-structure interfacing, and grid motion are integrated into a multidisciplinary computing environment that controls the temporal synchronization of the analysis modules. Peak values of the power spectral density of tail tip acceleration are reduced by as much as 22% in the first bending mode and by as much as 82% in the first torsion mode. RMS values of tip acceleration are reduced by as much as 12%.

Sheta, Essam F.↗

Surface Energy Heat Fluxes Using Remotely Sensed Parameters

Realistic estimates of surface energy heat fluxes are needed for the study of water and energy interactions between the land and atmosphere. The primary objective of this work is to study the estimation of surface heat energy fluxes using remote sensing derived parameters under different spatial and temporal conditions. Surface energy fluxes and remote sensing derived data from two sources were analyzed. First, we used surface heat flux, remote sensing, and ancillary data from the International Satellite Land Surface Climatology Project (ISLSCP), mapped at a 1 deg. x 1 deg. grid. Second, we used NOAA AVHRR (1 km), weather station, and ancillary data to derive estimates of surface latent and sensible heat energy fluxes over a 100 sq kilometers area for three test sites: 1) First ISLSCP Field Experiment (FIFE) grassland site, Konza Prairie, Kansas; 2) Howland, Maine Forest Ecosystem Dynamics Site; and 3) Walnut Gulch, scrubland site, surrounding Tombstone, Arizona. Satellite derived estimates of land surface temperature, surface albedo, and spectral vegetation index are used in selected models to provide estimates of surface heat fluxes. Analysis of results from the 1 deg. x 1 deg. grid for North America indicated there were similar, overall correlations between sensible and latent heat energy fluxes versus remotely sensed vegetation index and ground temperature during dry and wet year conditions. However, there were significant differences in correlations between years when stratified by land cover class. Analysis of 100 km x 100 km data (1 km resolution) indicated partitioning the areas in to primary versus secondary cover, with the secondary cover comprising less than 5% of the area, significantly improved surface heat energy flux estimates.

Toll, David L.↗

Tracking Magnetic Perturbations, dB/dt and Geomagnetic Indices for Geospace Storms on a Routine Basis at the CCMC

The first comprehensive assessment of geospace model skill to specify magnetic perturbations (delta-B) on the ground and their time derivative (dB/dt) was performed at the Community Coordinated Modeling Center starting in 2010. This study resulted in the addition of the Space Weather Modling Framework to the suite of operational models run by the NOAA Space Weather Prediction Center (SWPC). Since then, magnetic pertubations have been made available on a larger scale for Run-on-Request simulations in geospace at the CCMC. We will demonstrate recent additions to the suite of analysis tools and model results including magnetic perturbations at more stations and on a grid of positions from both, original model outputs (SWMF using preset run configurations) and post-processed calculations using CalcDeltaB and their analysis and visualization using the open-source Kamodo data access and analysis suite and the Comprehensive Assessment of Models and Events using Library tools (CAMEL) application.

Lutz Rastaetter↗

Comparison and analysis of aircraft measurements and mesoscale atmospheric chemistry model simulations of tropospheric ozone

The Regional Acid Deposition Model (RADM) has been applied to several of the field experiments which were part of the Acid Models Operational and Diagnostic Evaluation Study (Acid MODES). The experiment which was of particular interest with regards to ozone photochemistry involved horizontal zig-zag flight patterns (ZIPPER) over an area from the eastern Ohio River valley to the Adirondacks of New York. Model simulations by both the standard resolution RADM (delta x = 80 km) and the nested grid RADM (delta x = 26.7 km) compare well to measurements in the low emission regions in central Pennsylvania and upstate New York, but underestimate in the high emission upper Ohio River valley. The nested simulation does considerably better, however, than the coarse grid simulation in terms of horizontal pattern and concentration magnitudes. Analysis of NO(x) and HO(x) concentrations and photochemical products rates of ozone show that the model's response to large point source emissions is very unsystematic both spatially and temporally. This is due to the models instability to realistically simulate the small scale (subgrid) gradients in precursor concentrations in and around large point source plumes.

Pleim, Jonathan E.↗

An Active Smart Material Control System for F/A-18 Buffet Alleviation

The vertical tail buffet problem of fighter aircraft occurs at high angles of attack when the vortical flow breaks down ahead of the vertical tails resulting in unsteady and unbalanced pressure loads on the vertical tails. The buffet loads imposed upon the vertical tails resulted in a premature fatigue failure of the tails, and consequently limits the performance and super maneuverability of twin-tail fighter aircraft. An active smart material control system using distributed piezoelectric actuators has been developed for buffet alleviation and is presented. The inboard and outboard surfaces of the vertical tail are equipped with piezoelectric actuators to control the buffet responses in the first bending and torsion modes. The electrodynamics of the piezoelectric actuators are expressed with a three-dimensional finite-element model. A single-input-single-output controller is designed to drive the active piezoelectric actuators. High-fidelity multidisciplinary analysis modules for the fluid dynamics, structure dynamics, electrodynamics of the piezoelectric actuators, control law, fluid structure interfacing, and grid motion are integrated into a multidisciplinary computing environment that controls the temporal synchronization of the analysis modules. At 30 degree angle of attack, RMS values of tip acceleration are reduced by as much as 12%. The peak values of the power spectral density of tail-tip acceleration are reduced by as much as 22% in the first bending mode and by as much as 82% in the first torsion mode. The actively controlled piezoelectric actuators were also effective in adding damping at wide range of angles of attack.

Sheta, Essam F.↗

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗

Multi-resolution analysis for ENO schemes

Given an function, u(x), which is represented by its cell-averages in cells which are formed by some unstructured grid, we show how to decompose the function into various scales of variation. This is done by considering a set of nested grids in which the given grid is the finest, and identifying in each locality the coarsest grid in the set from which u(x) can be recovered to a prescribed accuracy. This multi-resolution analysis was applied to essentially non-oscillatory (ENO) schemes in order to advance the solution by one time-step. This is accomplished by decomposing the numerical solution at the beginning of each time-step into levels of resolution, and performing the computation in each locality at the appropriate coarser grid. An efficient algorithm for implementing this program in the 1-D case is presented; this algorithm can be extended to the multi-dimensional case with Cartesian grids.

Harten, Ami↗

Modification of a Successive Corrections Objective Analysis for Improved Derivative Calculations

The use of objectively analyzed fields of meteorological data for complex diagnostic studies and for the initialization of numerical prediction models places the requirements upon the objective method that derivatives of the gridded fields be accurate and free from interpolation error. A modification of an objective analysis developed by Barnes provides improvements in analyses of both the field and its derivatives. Theoretical comparisons, comparisons between analyses of analytical monochromatic waves, and comparisons between analyses of actual weather data are used to show the potential of the new method. The new method restores more of the amplitudes of desired wavelengths while simultaneously filtering more of the amplitudes of undesired wavelengths. These results also hold for the first and second derivatives calculated from the gridded fields. Greatest improvements were for the Laplacian of the height field; the new method reduced the variance of undesirable very short wavelengths by 72 pct. Other improvements were found in the divergence of the gridded wind field and near the boundaries of the field data.

Achtemeier, Gary L.↗

Propulsion Powertrain Real-Time Simulation Using Hardware-in-the-Loop (HIL) for Aircraft Electric Propulsion System

It is essential to design a propulsion powertrain real-time simulator using the hardware-in-the-loop (HIL) system that emulates an electrified aircraft propulsion (EAP) systems power grid. This simulator would enable us to facilitate in-depth understanding of the system principles, to validate system model analysis and performance prediction, and to demonstrate the proof-of-concept of the EAP electrical system. This paper describes how subscale electrical machines with their controllers can mimic the power components in an EAP powertrain. In particular, three powertrain emulations are presented to mimic 1) a gas turbo-=shaft engine driving a generator, consisting of two permanent magnet (PM) motors with brushless motor drives, coupled by a shaft, 2) a motor driving a propulsive fan, and 3) a turbo-shaft engine driven fan (turbofan engine) operation. As a first step towards the demonstration, experimental dynamic characterization of the two motor drive systems, coupled by a mechanical shaft, were performed. The previously developed analytical motor models1 were then replaced with the experimental motor models to perform the real-time demonstration in the predefined flight path profiles. This technique can convert the plain motor system into a unique EAP power grid emulator that enables rapid analysis and real-time simulation performance using hardware-in-the-loop (HIL).

turbo-electric propulsion↗

Addressing the Big-Earth-Data Variety Challenge with the Hierarchical Triangular Mesh

We have implemented an updated Hierarchical Triangular Mesh (HTM) as the basis for a unified data model and an indexing scheme for geoscience data to address the variety challenge of Big Earth Data. We observe that, in the absence of variety, the volume challenge of Big Data is relatively easily addressable with parallel processing. The more important challenge in achieving optimal value with a Big Data solution for Earth Science (ES) data analysis, however, is being able to achieve good scalability with variety. With HTM unifying at least the three popular data models, i.e. Grid, Swath, and Point, used by current ES data products, data preparation time for integrative analysis of diverse datasets can be drastically reduced and better variety scaling can be achieved. In addition, since HTM is also an indexing scheme, when it is used to index all ES datasets, data placement alignment (or co-location) on the shared nothing architecture, which most Big Data systems are based on, is guaranteed and better performance is ensured. Moreover, our updated HTM encoding turns most geospatial set operations into integer interval operations, gaining further performance advantages.

SciDB↗

SEASAT-A SASS wind processing

The techniques used for ingesting SEASAT-A SASS wind retrievals into the existing operational software are described. The intent is to assess the impact of SEASAT data in he marine wind fields produced by the global marine wind/sea level pressure analysis. This analysis is performed on a 21/2 deg latitude/longitude global grid which executes at three hourly time increments. Wind fields with and without SASS winds are being compared. The problems of data volume reduction and aliased wind retrieval ambiquity are treated.

Langland, R. A.↗

Calculation of transonic steady and oscillatory pressures on a low aspect ratio model and comparison with experiment

Pressure data measured by the British Royal Aircraft Establishment for the AGARD SMP tailplane are compared with results calculated using the transonic small perturbation code XTRAN3S. A brief description of the analysis is given and a recently developed finite difference grid is described. Results are presented for five steady and nine harmonically oscillating cases near zero angle of attack and for a range of subsonic and transonic Mach numbers.

Bennett, R. M.↗

Calculation of transonic steady and oscillatory pressures on a low aspect ratio model and comparison with experiment

Pressure data measured by the British Royal Aircraft Establishment for the AGARD SMP tailplane are compared with results calculated using the transonic small perturbation code XTRAN3S. A brief description of the analysis is given and a recently developed finite difference grid is described. Results are presented for five steady and nine harmonically oscillating cases near zero angle of attack and for a range of subsonic and transonic Mach numbers.

Bennett, R. M.↗

Development of liquid xenon detectors for gamma ray astronomy

The application of liquid xenon in high-resolution detectors for gamma-ray astronomy is being investigated. Initial results from a pulse-shape analysis of ionization signals in a liquid-xenon gridded chamber indicate that it is possible to achieve the necessary liquid purity for the transport of free electrons with simple techniques. The energy resolution has been measured as a function of applied electric field, using electrons and gamma-rays from a 207Bi source. At a field of 12 kV/cm the noise-substracted energy resolution of the dominant 569-keV gamma-ray line is 34 keV FWHM (full width at half maximum). This value is mostly determined by recombination of electron-ion pairs on delta-electron tracks.

Aprile, Elena↗

Conical Euler simulation and active suppression of delta wing rocking motion

A conical Euler code was developed to study unsteady vortex-dominated flows about rolling highly-swept delta wings, undergoing either forced or free-to-roll motions including active roll suppression. The flow solver of the code involves a multistage Runge-Kutta time-stepping scheme which uses a finite volume spatial discretization of the Euler equations on an unstructured grid of triangles. The code allows for the additional analysis of the free-to-roll case, by including the rigid-body equation of motion for its simultaneous time integration with the governing flow equations. Results are presented for a 75 deg swept sharp leading edge delta wing at a freestream Mach number of 1.2 and at alpha equal to 10 and 30 deg angle of attack. A forced harmonic analysis indicates that the rolling moment coefficient provides: (1) a positive damping at the lower angle of attack equal to 10 deg, which is verified in a free-to-roll calculation; (2) a negative damping at the higher angle of attack equal to 30 deg at the small roll amplitudes. A free-to-roll calculation for the latter case produces an initially divergent response, but as the amplitude of motion grows with time, the response transitions to a wing-rock type of limit cycle oscillation. The wing rocking motion may be actively suppressed, however, through the use of a rate-feedback control law and antisymmetrically deflected leading edge flaps. The descriptions of the conical Euler flow solver and the free-to-roll analysis are presented. Results are also presented which give insight into the flow physics associated with unsteady vortical flows about forced and free-to-roll delta wings, including the active roll suppression of this wing-rock phenomenon.

Lee, Elizabeth M.↗

Aerodynamic shape optimization via sensitivity analysis on decomposed computational domains

Direct and iterative method considered to be most applicable to large systems of linear equations arising in discrete sensitivity analysis are assessed. Based on a single-domain grid, computations are performed using a banded matrix solver and an iterative solver, the generalized minimum residual (GMRES) method. The banded matrix solver is found to be generally the most economical method for those applications where the number of right-hand sides is large (i.e., a large number of design variables or a large number of adjoint vectors). For systems of equations that are too large to be solved by direct methods, an approach is proposed whereby the computational domain is divided into small subdomains, and each subdomain is solved separately.

Eleshaky, Mohamed E.↗

Future requirements in surface modeling and grid generation

The past ten years have seen steady progress in surface modeling procedures, and wholesale changes in grid generation technology. Today, it seems fair to state that a satisfactory grid can be developed to model nearly any configuration of interest. The issues at present focus on operational concerns such as cost and quality. Continuing evolution of the engineering process is placing new demands on the technologies of surface modeling and grid generation. In the evolution toward a multidisciplinary analysis-bascd design environment, methods developed for Computational Fluid Dynamics are finding acceptance in many additional applications. These two trends, the normal evolution of the process and a watershed shift toward concurrent and multidisciplinary analysis, will be considered in assessing current capabilities and needed technological improvements.

Cosner, Raymond R.↗