Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “precondition”

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 361 records · Page 20

Simulation of Synthetic Jets in Quiescent Air Using Unsteady Reynolds Averaged Navier-Stokes Equations

We report research experience in applying an Unsteady Reynolds-Averaged Navier-Stokes (URANS) solver for the prediction of time-dependent flows in the presence of an active flow control device. The configuration under consideration is a synthetic jet created by a single diaphragm piezoelectric actuator in quiescent air. Time-averaged and instantaneous data for this case were obtained at Langley Research Center, using multiple measurement techniques. Computational results for this case using one-equation Spalart-Allmaras and two-equation Menter s turbulence models are presented here along with comparisons with the experimental data. The effect of grid refinement, preconditioning and time-step variation are also examined.

Vatsa, Veer N.↗

Evaluation of Isolated Fuselage and Rotor-Fuselage Interaction Using CFD

The US Army Aeroflightdynamics Directorate (AFDD), the French Office National d'Etudes et de Recherches Aerospatiales (ONERA) and the Georgia Institute of Technology (GIT) are working under the United States/France Memorandum of Agreement on Helicopter Aeromechanics to study rotorcraft aeromechanics issues of interest to both nations. As a task under this agreement, a comparative study of the Dauphin 365N helicopter has been undertaken to analyze the capabilities and weaknesses of state-of-the-art computational fluid dynamics (CFD) codes, with the aim of fuselage performance prediction and investigation of rotor-fuselage interaction. Three CFD flow solvers applied on three meshes provide similar results in terms of pressure coefficient. Force predictions vary somewhat. This paper presents details on the grid sensitivity and the low Mach number preconditioning influence. The importance of taking into account the wind tunnel strut and the rotor hub is shown. The pressure coefficients along top and bottom centerlines of the fuselage are in good agreement with the experiment except in the area aft of the hub. There remains a discrepancy between the computed forces and the experimental data due in part to modeling inaccuracies. Rotor- fuselage interactions are performed using uniform and non-uniform actuator disk models in order to simulate the rotor downwash.

Renaud, Thomas↗

Impact of Variations on 1-D Flow in Gas Turbine Engines via Monte Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented using a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, to determine the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions motivate this research: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance? This paper presents an initial attempt to model the flow behavior in terms of their eigenfrequencies subject to the assumption of the uncertainty or variation (perturbation). The flow behavior is explored using simulation outputs from a customer-deck model obtained from Pratt & Whitney. Variations of the main variables (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine are modeled. Flow behavior is analyzed for the high-pressure compressor and combustion chamber utilizing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high-pressure and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000 sample simulation. Furthermore, the characteristic waves and wave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures in the case of the combustion chamber and high-pressure turbine. Finally, a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗

Vectran Fiber Time-Dependent Behavior and Additional Static Loading Properties

Vectran HS appears from literature and testing to date to be an ideal upgrade from Kevlar braided cords for many long-term, static-loading applications such as tie-downs on solar arrays. Vectran is a liquid crystalline polymer and exhibits excellent tensile properties. The material has been touted as a zero creep product. Testing discussed in this report does not support this statement, though the creep is on the order of four times slower than with similar Kevlar 49 products. Previous work with Kevlar and new analysis of Vectran testing has led to a simple predictive model for Vectran at ambient conditions. The mean coefficient of thermal expansion (negative in this case) is similar to Kevlar 49, but is not linear. A positive transition in the curve occurs near 100 C. Out-gassing tests show that the material performs well within parameters for most space flight applications. Vectran also offers increased abrasion resistance, minimal moisture regain, and similar UV degradation. The effects of material construction appear to have a dramatic effect in stress relaxation for braided Vectran. To achieve the improved relaxation rate, upgrades must also examine alternate construction or preconditioning methods. This report recommends Vectran HS as a greatly improved replacement material for applications where time-dependent relaxation is a major factor.

Fette, Russell B.↗

SWIFT: Multistage Turbomachinery Analysis Code Developed

A new turbomachinery analysis code called SWIFT has been developed at the NASA Lewis Research Center. SWIFT solves thin-layer Navier-Stokes equations with the Baldwin-Lomax turbulence model and an explicit finite-difference scheme. Preconditioning allows the code to be used for all speed ranges, from incompressible to supersonic flows. Multiblock capability allows three types of grids to be patched together to simulate many types of turbomachinery geometries, including hub and tip clearances. In addition, code users can analyze multistage turbomachinery by using a steady averaging-plane approach. This approach uses linearized characteristic boundary conditions to pass information accurately between the stages.

Chima, Rodrick V.↗

Convergence Acceleration and Documentation of CFD Codes for Turbomachinery Applications

The development and analysis of turbomachinery components for industrial and aerospace applications has been greatly enhanced in recent years through the advent of computational fluid dynamics (CFD) codes and techniques. Although the use of this technology has greatly reduced the time required to perform analysis and design, there still remains much room for improvement in the process. In particular, there is a steep learning curve associated with most turbomachinery CFD codes, and the computation times need to be reduced in order to facilitate their integration into standard work processes. Two turbomachinery codes have recently been developed by Dr. Daniel Dorney (MSFC) and Dr. Douglas Sondak (Boston University). These codes are entitled Aardvark (for 2-D and quasi 3-D simulations) and Phantom (for 3-D simulations). The codes utilize the General Equation Set (GES), structured grid methodology, and overset O- and H-grids. The codes have been used with success by Drs. Dorney and Sondak, as well as others within the turbomachinery community, to analyze engine components and other geometries. One of the primary objectives of this study was to establish a set of parametric input values which will enhance convergence rates for steady state simulations, as well as reduce the runtime required for unsteady cases. The goal is to reduce the turnaround time for CFD simulations, thus permitting more design parametrics to be run within a given time period. In addition, other code enhancements to reduce runtimes were investigated and implemented. The other primary goal of the study was to develop enhanced users manuals for Aardvark and Phantom. These manuals are intended to answer most questions for new users, as well as provide valuable detailed information for the experienced user. The existence of detailed user s manuals will enable new users to become proficient with the codes, as well as reducing the dependency of new users on the code authors. In order to achieve the objectives listed, the following tasks were accomplished: 1) Parametric Study Of Preconditioning Parameters And Other Code Inputs; 2) Code Modifications To Reduce Runtimes; 3) Investigation Of Compiler Options To Reduce Code Runtime; and 4) Development/Enhancement of Users Manuals for Aardvark and Phantom

Marquart, Jed E.↗

Matching Multistage Schemes to Viscous Flow

Multistage, explicit time stepping can be tailored to accelerate convergence for scalar advection-diffusion problems by using optimized multistage coefficients that vary with the local cell Reynolds number. And, when combined with local preconditioning, variable-coefficient multistage schemes for computational fluid dynamics codes can also provide an order of magnitude faster convergence, relative to standard, fixed-coefficient schemes, for the Navier-Stokes system of equations.

Kleb, William L.↗

Accuracy Quantification of the Loci-CHEM Code for Chamber Wall Heat Transfer in a GO2/GH2 Single Element Model Problem

All solutions with Loci-CHEM achieved demonstrated steady state and mesh convergence. Preconditioning had no effect on solution accuracy and typically yields a 3-5times solution speed-up. The SST turbulence model has superior performance, relative to the data in the head end region, for the rise rate and peak heat flux. It was slightly worse than the others in the downstream region where all over-predicted the data by 30-100%.There was systematic mesh refinement in the unstructured volume and structured boundary layer areas produced only minor solution differences. Mesh convergence was achieved. Overall, Loci-CHEM satisfactorily predicts heat flux rise rate and peak heat flux and significantly over predicts the downstream heat flux.

West, Jeff↗

Convergence Acceleration for Multistage Time-Stepping Schemes

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 could be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. Numerical dissipation operators (based on the Roe scheme, a matrix formulation, and the CUSP scheme) as well as the number of RK stages are considered in evaluating the RK/implicit scheme. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. In two dimensions, turbulent flows over an airfoil at subsonic and transonic conditions are computed. The effects of mesh cell aspect ratio on convergence are investigated for Reynolds numbers between 5.7 x 10(exp 6) and 100.0 x 10(exp 6). Results are also obtained for a transonic wing flow. For both 2-D and 3-D problems, the computational time of a well-tuned standard RK scheme is reduced at least a factor of four.

Swanson, R. C.↗

Climate Data Assimilation on a Massively Parallel Supercomputer

We have designed and implemented a set of highly efficient and highly scalable algorithms for an unstructured computational package, the PSAS data assimilation package, as demonstrated by detailed performance analysis of systematic runs on up to 512-nodes of an Intel Paragon. The preconditioned Conjugate Gradient solver achieves a sustained 18 Gflops performance. Consequently, we achieve an unprecedented 100-fold reduction in time to solution on the Intel Paragon over a single head of a Cray C90. This not only exceeds the daily performance requirement of the Data Assimilation Office at NASA's Goddard Space Flight Center, but also makes it possible to explore much larger and challenging data assimilation problems which are unthinkable on a traditional computer platform such as the Cray C90.

supercomputer↗

A General Purpose Sparse Matrix Parallel Solvers Package

A general purpose solver package for constructing and solving a range of sparse linear systems arising from discretization of PDEs on unstructured meshes is developed. Once the sparse symmetric complex matrix is constructed, it can be solved by either a preconditioned bi-conjugate gradient solver, a two-stage Cholesky LDLT factorization solver, or a hybrid solver combining the above two methods.

solver sparse matrix solvers package PDE↗

A General Purpose Sparse Matrix Parallel Solvers Package

A general purpose solver package for constructing and solving a range of sparse linear systems arising from discretization of PDEs on unstructured meshes is developed. Once the sparse symmetric complex matrix is constructed, it can be solved by either a preconditioned bi-conjugate gradient solver, a two-stage Cholesky LDLT factorization solver, or a hybrid solver combining the above two methods. (More detailed than 95-0127).

solver sparse matrix solvers package PDE↗

On large outflows of Arctic sea ice into the Barents Sea

Winter outflows of Arctic sea ice into the Barents Sea are estimated using a 10-year record of satellite ice motion and thickness. The mean winter volume export through the Svalbard/Franz Josef Land passage is 40 km3, and ranges from -280 km3 to 340 km3. A large outflow in 2003 is preconditioned by an unusually high concentration of thick perennial ice over the Nansen Basin at the end of the 2002 summer. With a deep atmospheric low situated over the eastern Barents Sea in winter, the result is an increased export of Arctic ice. The Oct-Mar ice area flux, at 110 x 10 to the third power km3, is not only unusual in magnitude but also remarkable in that >70% of the area is multiyear ice; the ice volume flux at~340 km3 is almost one-fifth of the ice flux through the Fram Strait. Another large outflow of Arctic sea ice through this passage, comparable to that in 2003, is found in 1996. This southward flux of sea ice represents one of two major sources of freshwater in the Barents Sea; the other is the eastward flux of water via the Norwegian Coastal Current. The possible consequences of variable freshwater input on the Barents Sea hydrography and its impact on transformation of Atlantic Water en route to the Arctic Ocean are examined with a 25-year coupled ice-ocean model.

Artic ice↗

Modeling of Heat and Mass Transfer in a TEC-Driven Lyophilizer

Dewatering of wet waste during space exploration missions is important for crew safety as it stabilizes the waste. It may also be used to recover water and serve as a preconditioning step for waste compaction. A thermoelectric cooler (TEC)-driven lyophilizer is under development at NASA Ames Research Center for this purpose. It has three major components: (i) an evaporator section where water vapor sublimes from the frozen waste, (ii) a condenser section where this water vapor deposits as ice, and (iii) a TEC section which serves as a heat pump to transfer heat from the condenser to the evaporator. This paper analyses the heat and mass transfer processes in the lyophilizer in an effort to understand the ice formation behavior in the condenser. The analysis is supported by experimental observations of ice formation patterns in two different condenser units.

Yuan, Zeng-Guang↗

Simulations of Dynamics and Transport during the September 2002 Antarctic Major Warming

A mechanistic model simulation initialized on 14 September 2002, forced by 100-hPa geopotential heights from Met Office analyses, reproduced the dynamical features of the 2002 Antarctic major warming. The vortex split on approx.25 September; recovery after the warming, westward and equatorward tilting vortices, and strong baroclinic zones in temperature associated with a dipole pattern of upward and downward vertical velocities were all captured in the simulation. Model results and analyses show a pattern of strong upward wave propagation throughout the warming, with zonal wind deceleration throughout the stratosphere at high latitudes before the vortex split, continuing in the middle and upper stratosphere and spreading to lower latitudes after the split. Three-dimensional Eliassen-Palm fluxes show the largest upward and poleward wave propagation in the 0(deg)-90(deg)E sector prior to the vortex split (coincident with the location of strongest cyclogenesis at the model's lower boundary), with an additional region of strong upward propagation developing near 180(deg)-270(deg)E. These characteristics are similar to those of Arctic wave-2 major warmings, except that during this warming, the vortex did not split below approx.600 K. The effects of poleward transport and mixing dominate modeled trace gas evolution through most of the mid- to high-latitude stratosphere, with a core region in the lower-stratospheric vortex where enhanced descent dominates and the vortex remains isolated. Strongly tilted vortices led to low-latitude air overlying vortex air, resulting in highly unusual trace gas profiles. Simulations driven with several meteorological datasets reproduced the major warming, but in others, stronger latitudinal gradients at high latitudes at the model boundary resulted in simulations without a complete vortex split in the midstratosphere. Numerous tests indicate very high sensitivity to the boundary fields, especially the wave-2 amplitude. Major warmings occurred for initial fields with stronger winds and larger vortices, but not smaller vortices, consistent with the initiation of wind-deceleration by upward-propagating waves near the poleward edge of the region where wave 2 can propagate above the jet core. Thus, given the observed 100-hPa boundary forcing, stratospheric preconditioning is not needed to reproduce a major warming similar to that observed. The anomalously strong forcing in the lower stratosphere can be viewed as the primary direct cause of the major warming.

winds↗

Annual Cycles of Multiyear Sea Ice Coverage of the Arctic Ocean: 1999-2003

For the years 1999-2003, we estimate the time-varying perennial ice zone (PIZ) coverage and construct the annual cycles of multiyear (MY, including second year) ice coverage of the Arctic Ocean using QuikSCAT backscatter, MY fractions from RADARSAT, and the record of ice export from satellite passive microwave observations. An area balance approach extends the winter MY coverage from QuikSCAT to the remainder of the year. From these estimates, the coverage of MY ice at the beginning of each year is 3774 x 10(exp 3) sq km (2000), 3896 x 10(exp 3) sq km (2001), 4475 x 10(exp 3) sq km (2002), and 4122 x 10(exp 3) sq km (2003). Uncertainties in coverage are approx.150 x 10(exp 3) sq km. In the mean, on 1 January, MY ice covers approx.60% of the Arctic Ocean. Ice export reduces this coverage to approx.55% by 1 May. From the multiple annual cycles, the area of first-year (FY) ice that survives the intervening summers are 1192 x 10(exp 3) sq km (2000), 1509 x 10(exp 3) sq km (2001), and 582 x 10(exp 3) sq km (2002). In order for the MY coverage to remain constant from year to year, these replenishment areas must balance the overall area export and melt during the summer. The effect of the record minimum in Arctic sea ice area during the summer of 2002 is seen in the lowest area of surviving FY ice of the three summers. In addition to the spatial coverage, the location of the PIZ is important. One consequence of the unusual location of the PIZ at the end of the summer of 2002 is the preconditioning for enhanced export of MY ice into the Barents and Kara seas. Differences between the minimums in summer sea ice coverage from our estimates and passive microwave observations are discussed.

Arctic Ocean↗

Using the GeoFEST Faulted Region Simulation System

GeoFEST (the Geophysical Finite Element Simulation Tool) simulates stress evolution, fault slip and plastic/elastic processes in realistic materials, and so is suitable for earthquake cycle studies in regions such as Southern California. Many new capabilities and means of access for GeoFEST are now supported. New abilities include MPI-based cluster parallel computing using automatic PYRAMID/Parmetis-based mesh partitioning, automatic mesh generation for layered media with rectangular faults, and results visualization that is integrated with remote sensing data. The parallel GeoFEST application has been successfully run on over a half-dozen computers, including Intel Xeon clusters, Itanium II and Altix machines, and the Apple G5 cluster. It is not separately optimized for different machines, but relies on good domain partitioning for load-balance and low communication, and careful writing of the parallel diagonally preconditioned conjugate gradient solver to keep communication overhead low. Demonstrated thousand-step solutions for over a million finite elements on 64 processors require under three hours, and scaling tests show high efficiency when using more than (order of) 4000 elements per processor. The source code and documentation for GeoFEST is available at no cost from Open Channel Foundation. In addition GeoFEST may be used through a browser-based portal environment available to approved users. That environment includes semi-automated geometry creation and mesh generation tools, GeoFEST, and RIVA-based visualization tools that include the ability to generate a flyover animation showing deformations and topography. Work is in progress to support simulation of a region with several faults using 16 million elements, using a strain energy metric to adapt the mesh to faithfully represent the solution in a region of widely varying strain.

Geophyical Finite Element Simulation Tool (GeoFEST↗

Effects of Gravity on Start-up of a Miniature Loop Heat Pipe with Multiple Evaporators and Multiple Condensers

This paper describes an experimental study on the effect of gravity on the start-up of a miniature loop heat pipe (MLHP) with two evaporators and two condensers. Each evaporator has an outer diameter of 9 mm and has its own integral compensation chamber (CC). The MLHP was placed under five configurations where the relative elevation and tilt among the loop components were varied. The four well-known initial conditions between the evaporator and CC prior to the LHP start-up were created in this experimental study through combinations of: 1) the test configuration; 2) the method of preconditioning the loop prior to start-up, and 3) the heat load distribution among the evaporators. A total of 165 start-ups were conducted under the five test configurations. All of these start-ups were successful. However, the effect of gravity on start-up transients was clearly seen under otherwise the same heat load distribution and sink temperatures. An analytical model was used to simulate the MLHP transient behaviors, and the model predictions agreed very well with the experimental results.

Ku, Jentung↗