Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence 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 937 records · Page 52

ALPS: The Age-Layered Population Structure for Reducing the Problem of Premature Convergence

To reduce the problem of premature convergence we define a new attribute of an individual, its age, and propose the Age-Layered Population Structure (ALPS), in which age is used to restrict competition and breeding between members of the population. ALPS differs from a typical EA by segregating individuals into different age-layers by their age - a measure of how long the genetic material has been in the population - and by regularly replacing all individuals in the bottom layer with randomly generated ones. The introduction of new, randomly generated individuals at regular intervals results in an EA that is never completely converged and is always looking at new parts of the fitness landscape. By using age to restrict competition and breeding search is able to develop promising young individuals without them being dominated by older ones. We demonstrate the effectiveness of the ALPS algorithm on an antenna design problem in which evolution with ALPS produces antennas more than twice as good as does evolution with two other types of EAs. Further analysis shows that the ALPS model does allow the offspring of newly generated individuals to move the population out of mediocre local-optima to better parts of the fitness landscape.

Hornby, Gregory S.↗

The Impact of High-Resolution Sea Surface Temperatures on the Simulated Nocturnal Florida Marine Boundary Layer

High- and low-resolution sea surface temperature (SST) analysis products are used to initialize the Weather Research and Forecasting (WRF) Model for May 2004 for short-term forecasts over Florida and surrounding waters. Initial and boundary conditions for the simulations were provided by a combination of observations, large-scale model output, and analysis products. The impact of using a 1-km Moderate Resolution Imaging Spectroradiometer (MODIS) SST composite on subsequent evolution of the marine atmospheric boundary layer (MABL) is assessed through simulation comparisons and limited validation. Model results are presented for individual simulations, as well as for aggregates of easterly- and westerly-dominated low-level flows. The simulation comparisons show that the use of MODIS SST composites results in enhanced convergence zones. earlier and more intense horizontal convective rolls. and an increase in precipitation as well as a change in precipitation location. Validation of 10-m winds with buoys shows a slight improvement in wind speed. The most significant results of this study are that 1) vertical wind stress divergence and pressure gradient accelerations across the Florida Current region vary in importance as a function of flow direction and stability and 2) the warmer Florida Current in the MODIS product transports heat vertically and downwind of this heat source, modifying the thermal structure and the MABL wind field primarily through pressure gradient adjustments.

LaCasse, Katherine M.↗

Low-Level Jets and the Convergence of Mars Data Assimilation Algorithms

Data assimilation is an increasingly popular technique in Mars atmospheric science, but its effect on the mean states of the underlying atmosphere models has not been thoroughly examined. The robustness of results to the choice of model and assimilation algorithm also warrants further study. We investigate these issues using two Mars general circulation models (MGCMs), with particular emphasis on zonal wind and temperature fields. When temperature retrievals from the Mars Global Surveyor Thermal Emission Spectrometer (TES) are assimilated into the U.K.-Laboratoire de M ́et ́eorologie Dynamique (UK-LMD) MGCM to create the Mars Analysis Correction Data Assimilation (MACDA) reanalysis, low-level zonal jets in the winter northern hemisphere shift equatorward and weaken relative to a free-running control simulation from the same MGCM. The Ensemble Mars Atmosphere Reanalysis System (EMARS) reanalysis, which is also based on TES temperature retrievals, also shows jet weakening (but less if any shifting) relative to a control simulation performed with the underlying Geophysical Fluid Dynamics Laboratory (GFDL) MGCM. Examining higher levels of the atmosphere, monthly mean three-dimensional temperature and zonal wind fields are in generally better agreement between the two reanalyses than between the two control simulations. In conjunction with information about the MGCMs’ physical parametrizations, intercomparisons between the various reanalyses and control simulations suggest that overall the EMARS control run is plausibly less biased (relative to the true state of the Martian atmosphere) than the MACDA control run. Implications for future observational studies are discussed. Plain Language Summary An increasingly popular way to study Martian weather and climate is to combine atmospheric temperature observations with a computer model (specifically, a global climate model). The process of combining model and observations is called “data assimilation”, and the resulting merged data set is called a “reanalysis”. One advantage of reanalyses is that they include variables (such as wind) that are not directly observed. For scientific and practical applications we want these variables to be reasonably accurate— however, it is not clear how well data assimilation algorithms compute them. Our study investigates this issue using two Mars reanalyses and two model simulations that do not assimilate temperature data. We focus on slowly-varying atmospheric phenomena (timescales from 10 Mars days to a season). Assimilating temperature data into two different global climate models changes the strength and/or spatial pattern of east-west winds at low altitudes. Furthermore, monthly mean three-dimensional temperature and east-west wind fields agree better between reanalyses than between non-assimilating model simulations. This suggests that the data assimilation process is basically successful. One non-assimilating model simulation has less realistic representations of atmospheric physical processes than the other—we argue that this plausibly gives it larger biases relative to the true state of the atmosphere.

Todd A Mooring↗

Oceanic geoid and tides derived from GEOS 3 satellite data in the Northwestern Atlantic Ocean

Two sets of GEOS 3 altimeter data which fall within about a 2.5-deg width are analyzed for ocean geoid and tides. One set covers a path from Newfoundland to Cuba, and the other a path from Puerto Rico to the North Carolina coast. Forty different analyses using various parameters are performed in order to investigate convergence. Profiles of the geoid and four tides, M2, O1, S2, and K1, are derived along the two strips. While the analyses produced convergent solutions for all 40 cases, the uncertainty caused by the linear orbital bias error of the satellite is too large to claim that the solutions represent the true ocean tides in the area. A spot check of the result with the Mode deep-sea tide gauge data shows poor agreement. A positive conclusion of this study is that despite the uncertain orbital error the oceanic geoid obtained through this analysis can improve significantly the short-wavelength structure over existing spherical harmonic geoid models.

Won, I. J.↗

Analysis of airfoil transitional separation bubbles

A previously developed local inviscid-viscous interaction technique for the analysis of airfoil transitional separation bubbles, ALESEP (Airfoil Leading Edge Separation), has been modified to utilize a more accurate windward finite difference procedure in the reversed flow region, and a natural transition turbulence model has been incorporated for the prediction of transition within the separation bubble. Numerous calculations and experimental comparisons are presented to demonstrate the effects of the windward differencing scheme and the natural transition turbulence model. Grid sensitivity and convergence capabilities of this inviscid-viscous interaction technique are briefly addressed. A major conclusion of this paper is that a second, counter-rotating eddy has been found to exist in the wall layer of the primary separation bubble with the use of windward differencing.

Davis, R. L.↗

Propulsive-jet flow field analysis using the three-dimensional Navier-Stokes equations

A three-dimensional Navier-Stokes code has been applied to the analysis of flow fields containing propulsive jets. Specifically, the application was made to a flow field containing a supersonic jet injected at an angle of 90 degrees to a subsonic free stream. Although wind tunnel data were available, the computational results were not readily comparable to the experimental data because of significant differences between the two plume trajectories. Reasons for the differences are suggested in the report and include: (1) incomplete convergence, (2) inadequate grid resolution in the high gradient regions, and (3) use of a low-order turbulence closure model.

Reed, C. L.↗

The cloud and surface contributions to the net radiation balance during the 1982-83 El Nino

This paper investigates the interannual fluctuations in cloud and surface forcing associated with the unusually intense ENSO event of 1982-1983. The analysis was performed over a six-year period for which global radiation budget and cloudiness measurements from the Nimbus-7 experiments were available simultaneously. It is shown that major anomalies in cloudiness, due mostly to shifts of a major convective center in the Pacific Ocean associated with inactive or active ENSO phase, are most notable over the equatorial central Pacific, the maritime continent over the western Pacific, the South Pacific convergence zone, and the Indian Ocean.

Sohn, Byung-Ju↗

Modal truncation, Ritz vectors, and derivatives of closed-loop damping ratios

The effect of modal truncation on the damping ratio and their derivatives with respect to an added mass is investigated for a simply supported, multispan beam with a linear quadratic Gaussian control system. It is found that both the damping ratios and derivatives converge slowly, but the derivatives converge more slowly than the damping ratios. However, it is shown that when Ritz vectors corresponding to static displacements due to actuator forces are added to the reduced model, the convergence of both the damping ratios and their derivatives is accelerated. It is also shown that the accuracy of the damping ratio predicted by a reduced-model control design can be improved significantly if the Ritz vectors are included in the design of the control system. Thus, it appears that Ritz vectors added to the reduced model of flexible structures can improve greatly the accuracy of both the design and analysis of the control system.

Sandridge, Chris A.↗

Applications of computer algebra to distributed parameter systems

In the analysis of vibrations of continuous elastic systems, one often encounters complicated transcendental equations with roots directly related to the system's natural frequencies. Typically, these equations contain system parameters whose values must be specified before a numerical solution can be obtained. The present paper presents a method whereby the fundamental frequency can be obtained in analytical form to any desired degree of accuracy. The method is based upon truncation of rapidly converging series involving inverse powers of the system natural frequencies. A straightforward method to developing these series and summing them in closed form is presented. It is demonstrated how Computer Algebra can be exploited to perform the intricate analytical procedures which otherwise would render the technique difficult to apply in practice. We illustrate the method by developing two analytical approximations to the fundamental frequency of a vibrating cantilever carrying a rigid tip body. The results are compared to the numerical solution of the exact (transcendental) frequency equation over a range of system parameters.

Storch, Joel A.↗

First-Order Model Management With Variable-Fidelity Physics Applied to Multi-Element Airfoil Optimization

First-order approximation and model management is a methodology for a systematic use of variable-fidelity models or approximations in optimization. The intent of model management is to attain convergence to high-fidelity solutions with minimal expense in high-fidelity computations. The savings in terms of computationally intensive evaluations depends on the ability of the available lower-fidelity model or a suite of models to predict the improvement trends for the high-fidelity problem, Variable-fidelity models can be represented by data-fitting approximations, variable-resolution models. variable-convergence models. or variable physical fidelity models. The present work considers the use of variable-fidelity physics models. We demonstrate the performance of model management on an aerodynamic optimization of a multi-element airfoil designed to operate in the transonic regime. Reynolds-averaged Navier-Stokes equations represent the high-fidelity model, while the Euler equations represent the low-fidelity model. An unstructured mesh-based analysis code FUN2D evaluates functions and sensitivity derivatives for both models. Model management for the present demonstration problem yields fivefold savings in terms of high-fidelity evaluations compared to optimization done with high-fidelity computations alone.

Alexandrov, N. M.↗

Tropospheric Ozone Over a Tropical Atlantic Station in the Northern Hemisphere: Paramaribo, Surinam (6 deg N, 55 deg W)

We present an analysis of 2.5 years of weekly ozone soundings conducted at a new monitoring station in Paramaribo, Surinam (6 deg N,55 deg W). This is currently one of only three ozone sounding stations in the northern hemisphere (NH) tropics, and the only one in the equatorial Atlantic region. Paramaribo is part of the Southern Hemisphere ADditional Ozone Sounding program (SHADOZ). Due to its position close to the equator, the Inter Tropical Convergence Zone (ITCZ) passes over Paramaribo twice per year, which results in a semi-annual seasonality of many parameters including relative humidity and ozone. The dataset from Paramaribo is used to: (1) evaluate ozone variability relative to precipitation, atmospheric circulation patterns and biomass burning; (2) contrast ozone at the NH equatorial Atlantic with that at nearby southern hemisphere (SH) stations Natal (6 deg S,35 deg W) and Ascension (8 deg S,14 deg W); (3) compare the seasonality of tropospheric ozone with a satellite-derived ozone product: Tropical Tropospheric Ozone Columns from the Modified Residual method (MR-TTOC). We find that Paramaribo is a distinctly Atlantic station. Despite its position north of the equator, it resembles nearby SH stations during most of the year. Transport patterns in the lower and middle troposphere during February and March differ from SH stations, which leads to a seasonality of ozone with two maxima. MR-TTOC over Paramaribo does not match the observed seasonality of ozone due to the use of a SH ozone sonde climatology in the MR method. The Paramaribo ozone record is used to suggest an improvement for northern hemisphere MR-TTOC retrievals. We conclude that station Paramaribo shows unique features in the region, and clearly adds new information to the existing SHADOZ record.

Peters, W.↗

Coupled-Flow Simulation of HP-LP Turbines Has Resulted in Significant Fuel Savings

Our objective was to create a high-fidelity Navier-Stokes computer simulation of the flow through the turbines of a modern high-bypass-ratio turbofan engine. The simulation would have to capture the aerodynamic interactions between closely coupled high- and low-pressure turbines. A computer simulation of the flow in the GE90 turbofan engine's high-pressure (HP) and low-pressure (LP) turbines was created at GE Aircraft Engines under contract with the NASA Glenn Research Center. The three-dimensional steady-state computer simulation was performed using Glenn's average-passage approach named APNASA. The areas upstream and downstream of each blade row mutually interact with each other during engine operation. The embedded blade row operating conditions are modeled since the average passage equations in APNASA actively include the effects of the adjacent blade rows. The turbine airfoils, platforms, and casing are actively cooled by compressor bleed air. Hot gas leaks around the tips of rotors through labyrinth seals. The flow exiting the high work HP turbines is partially transonic and, therefore, has a strong shock system in the transition region. The simulation was done using 121 processors of a Silicon Graphics Origin 2000 (NAS 02K) cluster at the NASA Ames Research Center, with a parallel efficiency of 87 percent in 15 hr. The typical average-passage analysis mesh size per blade row was 280 by 45 by 55, or approx.700,000 grid points. The total number of blade rows was 18 for a combined HP and LP turbine system including the struts in the transition duct and exit guide vane, which contain 12.6 million grid points. Design cycle turnaround time requirements ran typically from 24 to 48 hr of wall clock time. The number of iterations for convergence was 10,000 at 8.03x10(exp -5) sec/iteration/grid point (NAS O2K). Parallel processing by up to 40 processors is required to meet the design cycle time constraints. This is the first-ever flow simulation of an HP and LP turbine. In addition, it includes the struts in the transition duct and exit guide vanes.

Veres, Joseph P.↗

Multistage Schemes with Multigrid for Euler and Navier-Strokes Equations: Components and Analysis

A class of explicit multistage time-stepping schemes with centered spatial differencing and multigrids are considered for the compressible Euler and Navier-Stokes equations. These schemes are the basis for a family of computer programs (flow codes with multigrid (FLOMG) series) currently used to solve a wide range of fluid dynamics problems, including internal and external flows. In this paper, the components of these multistage time-stepping schemes are defined, discussed, and in many cases analyzed to provide additional insight into their behavior. Special emphasis is given to numerical dissipation, stability of Runge-Kutta schemes, and the convergence acceleration techniques of multigrid and implicit residual smoothing. Both the Baldwin and Lomax algebraic equilibrium model and the Johnson and King one-half equation nonequilibrium model are used to establish turbulence closure. Implementation of these models is described.

Swanson, R. C.↗

Viscid/inviscid interaction analysis of thrust augmenting ejectors

A method was developed for calculating the static performance of thrust augmenting ejectors by matching a viscous solution for the flow through the ejector to an inviscid solution for the flow outside the ejector. A two dimensional analysis utilizing a turbulence kinetic energy model is used to calculate the rate of entrainment by the jets. Vortex panel methods are then used with the requirement that the ejector shroud must be a streamline of the flow induced by the jets to determine the strength of circulation generated around the shroud. In effect, the ejector shroud is considered to be flying in the velocity field of the jets. The solution is converged by iterating between the rate of entrainment and the strength of the circulation. This approach offers the advantage of including external influences on the flow through the ejector. Comparisons with data are presented for an ejector having a single central nozzle and Coanda jet on the walls. The accuracy of the matched solution is found to be especially sensitive to the jet flap effect of the flow just downstream of the ejector exit.

Bevilacqua, P. M.↗

Sonic Boom Propagation Codes Validated by Flight Test

The sonic boom propagation codes reviewed in this study, SHOCKN and ZEPHYRUS, implement current theory on air absorption using different computational concepts. Review of the codes with a realistic atmosphere model confirm the agreement of propagation results reported by others for idealized propagation conditions. ZEPHYRUS offers greater flexibility in propagation conditions and is thus preferred for practical aircraft analysis. The ZEPHYRUS code was used to propagate sonic boom waveforms measured approximately 1000 feet away from an SR-71 aircraft flying at Mach 1.25 to 5000 feet away. These extrapolated signatures were compared to measurements at 5000 feet. Pressure values of the significant shocks (bow, canopy, inlet and tail) in the waveforms are consistent between extrapolation and measurement. Of particular interest is that four (independent) measurements taken under the aircraft centerline converge to the same extrapolated result despite differences in measurement conditions. Agreement between extrapolated and measured signature duration is prevented by measured duration of the 5000 foot signatures either much longer or shorter than would be expected. The duration anomalies may be due to signature probing not sufficiently parallel to the aircraft flight direction.

Poling, Hugh W.↗

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗

Rotor Airloads Prediction Using Loose Aerodynamic Structural Coupling

This work couples a computational fluid dynamics (CFD) code and rotorcraft computational structural dynamics (CSD) code to calculate helicopter rotor airloads across a range of flight conditions. An iterative loose (weak) coupling methodology is used to couple the CFD and CSD codes on a per revolution, periodic basis. The CFD uses a high fidelity, Navier-Stokes, overset grid methodology with first principles-based wake capturing. Modifications are made to the CFD code for aeroelastic analysis. For a UH-60A Blackhawk helicopter, four challenging level flight conditions are computed: 1) low speed (u = 0.15) with blade-vortex interaction, 2) high speed (u = 0.37) with advancing blade negative lift, 3) high thrust with dynamic stall (u = 0.24), and 4) hover. Results are compared with UH-60A Airloads Program fight test data. Most importantly, for all cases the loose coupling methodology is shown to be stable, convergent, and robust with full coupling of normal force, pitching moment, and chord force. In comparison with flight test data, normal force and pitching moment magnitudes are in good agreement. For the high speed and dynamic stall cases a phase lag in comparison with the data is seen, nonetheless, the shapes of the curves are very good. Overall, the results are noteworthy improvement over lifting line aerodynamics used in rotorcraft comprehensive codes.

Potsdam, Mark↗

Geodetic Strain Analysis Tool

A geodetic software analysis tool enables the user to analyze 2D crustal strain from geodetic ground motion, and create models of crustal deformation using a graphical interface. Users can use any geodetic measurements of ground motion and derive the 2D crustal strain interactively. This software also provides a forward-modeling tool that calculates a geodetic velocity and strain field for a given fault model, and lets the user compare the modeled strain field with the strain field obtained from the user s data. Users may change parameters on-the-fly and obtain a real-time recalculation of the resulting strain field. Four data products are computed: maximum shear, dilatation, shear angle, and principal components. The current view and data dependencies are processed first. The remaining data products and views are then computed in a round-robin fashion to anticipate view changes. When an analysis or display parameter is changed, the affected data products and views are invalidated and progressively re-displayed as available. This software is designed to facilitate the derivation of the strain fields from the GPS and strain meter data that sample it to facilitate the understanding of the strengths and weaknesses of the strain field derivation from continuous GPS (CGPS) and other geodetic data from a variety of tectonic settings, to converge on the "best practices" strain derivation strategy for the Solid Earth Science ESDR System (SESES) project given the CGPS station distribution in the western U.S., and to provide SESES users with a scientific and educational tool to explore the strain field on their own with user-defined parameters.

Kedar, Sharon↗