Search NASA⌕ Search

SEARCH · Search NASA

Results for “preconditioning”

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 433 records · Page 24

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↗

Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations

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 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. 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. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.

Swanson, Roy C., Jr.↗

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

We apply an unsteady Reynolds-averaged Navier-Stokes (URANS) solver for the simulation of a synthetic jet created by a single diaphragm piezoelectric actuator in quiescent air. This configuration was designated as Case 1 for the CFDVAL2004 workshop held at Williamsburg, Virginia, in March 2004. Time-averaged and instantaneous data for this case were obtained at NASA 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 along with the experimental data. The effect of grid refinement, preconditioning and time-step variation are also examined in this paper.

Vatsa, Veer N.↗

Anticipating Cycle 24 Minimum and its Consequences: An Update

This Technical Publication updates estimates for cycle 24 minimum and discusses consequences associated with cycle 23 being a longer than average period cycle and cycle 24 having parametric minimum values smaller (or larger for the case of spotless days) than long term medians. Through December 2007, cycle 23 has persisted 140 mo from its 12-mo moving average (12-mma) minimum monthly mean sunspot number occurrence date (May 1996). Longer than average period cycles of the modern era (since cycle 12) have minimum-to-minimum periods of about 139.0+/-6.3 mo (the 90-percent prediction interval), inferring that cycle 24 s minimum monthly mean sunspot number should be expected before July 2008. The major consequence of this is that, unless cycle 24 is a statistical outlier (like cycle 21), its maximum amplitude (RM) likely will be smaller than previously forecast. If, however, in the course of its rise cycle 24 s 12-mma of the weighted mean latitude (L) of spot groups exceeds 24 deg, then one expects RM >131, and if its 12-mma of highest latitude (H) spot groups exceeds 38 deg, then one expects RM >127. High-latitude new cycle spot groups, while first reported in January 2008, have not, as yet, become the dominant form of spot groups. Minimum values in L and H were observed in mid 2007 and values are now slowly increasing, a precondition for the imminent onset of the new sunspot cycle.

Wilson, Robert M.↗

Computational Analyses of Pressurization in Cryogenic Tanks

A) Advanced Gas/Liquid Framework with Real Fluids Property Routines: I. A multi-fluid formulation in the preconditioned CRUNCH CFD(Registered TradeMark) code developed where a mixture of liquid and gases can be specified: a) Various options for Equation of state specification available (from simplified ideal fluid mixtures, to real fluid EOS such as SRK or BWR models). b) Vaporization of liquids driven by pressure value relative to vapor pressure and combustion of vapors allowed. c) Extensive validation has been undertaken. II. Currently working on developing primary break-up models and surface tension effects for more rigorous phase-change modeling and interfacial dynamics B) Framework Applied to Run-time Tanks at Ground Test Facilities C) Framework Used For J-2 Upper Stage Tank Modeling: 1) NASA MSFC tank pressurization: a) Hydrogen and oxygen tank pre-press, repress and draining being modeled at NASA MSFC. 2) NASA AMES tank safety effort a) liquid hydrogen and oxygen are separated by a baffle in the J-2 tank. We are modeling pressure rise and possible combustion if a hole develops in the baffle and liquid hydrogen leaks into the oxygen tank. Tank pressure rise rates simulated and risk of combustion evaluated.

Ahuja, Vineet↗

Automatic Testcase Generation for Flight Software

The TacSat3 project is applying Integrated Systems Health Management (ISHM) technologies to an Air Force spacecraft for operational evaluation in space. The experiment will demonstrate the effectiveness and cost of ISHM and vehicle systems management (VSM) technologies through onboard operation for extended periods. We present two approaches to automatic testcase generation for ISHM: 1) A blackbox approach that views the system as a blackbox, and uses a grammar-based specification of the system's inputs to automatically generate *all* inputs that satisfy the specifications (up to prespecified limits); these inputs are then used to exercise the system. 2) A whitebox approach that performs analysis and testcase generation directly on a representation of the internal behaviour of the system under test. The enabling technologies for both these approaches are model checking and symbolic execution, as implemented in the Ames' Java PathFinder (JPF) tool suite. Model checking is an automated technique for software verification. Unlike simulation and testing which check only some of the system executions and therefore may miss errors, model checking exhaustively explores all possible executions. Symbolic execution evaluates programs with symbolic rather than concrete values and represents variable values as symbolic expressions. We are applying the blackbox approach to generating input scripts for the Spacecraft Command Language (SCL) from Interface and Control Systems. SCL is an embedded interpreter for controlling spacecraft systems. TacSat3 will be using SCL as the controller for its ISHM systems. We translated the SCL grammar into a program that outputs scripts conforming to the grammars. Running JPF on this program generates all legal input scripts up to a prespecified size. Script generation can also be targeted to specific parts of the grammar of interest to the developers. These scripts are then fed to the SCL Executive. ICS's in-house coverage tools will be run to measure code coverage. Because the scripts exercise all parts of the grammar, we expect them to provide high code coverage. This blackbox approach is suitable for systems for which we do not have access to the source code. We are applying whitebox test generation to the Spacecraft Health INference Engine (SHINE) that is part of the ISHM system. In TacSat3, SHINE will execute an on-board knowledge base for fault detection and diagnosis. SHINE converts its knowledge base into optimized C code which runs onboard TacSat3. SHINE can translate its rules into an intermediate representation (Java) suitable for analysis with JPF. JPF will analyze SHINE's Java output using symbolic execution, producing testcases that can provide either complete or directed coverage of the code. Automatically generated test suites can provide full code coverage and be quickly regenerated when code changes. Because our tools analyze executable code, they fully cover the delivered code, not just models of the code. This approach also provides a way to generate tests that exercise specific sections of code under specific preconditions. This capability gives us more focused testing of specific sections of code.

Bushnell, David Henry↗

Towards an "All Speed" Unstructured Upwind Scheme

In the authors previous studies [1], a time-accurate, upwind finite volume method (ETAU scheme) for computing compressible flows on unstructured grids was proposed. The scheme is second order accurate in space and time and yields high resolution in the presence of discontinuities. The scheme features a multidimensional limiter and multidimensional numerical dissipation. These help to stabilize the numerical process and to overcome the annoying pathological behaviors of upwind schemes. In the present paper, it will be further shown that such multidimensional treatments also lead to a nearly all-speed or Mach number insensitive upwind scheme. For flows at very high Mach number, e.g., 10, local numerical instabilities or the pathological behaviors are suppressed, while for flows at very low Mach number, e.g., 0.02, computation can be directly carried out without invoking preconditioning. For flows in different Mach number regimes, i.e., low, medium, and high Mach numbers, one only needs to adjust one or two parameters in the scheme. Several examples with low and high Mach numbers are demonstrated in this paper. Thus, the ETAU scheme is applicable to a broad spectrum of flow regimes ranging from high supersonic to low subsonic, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics).

Loh, Ching Y.↗

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗