Search NASA⌕ Search

SEARCH · Search NASA

Results for “LU factorization”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

Numerical simulation of viscous flows over transonic aircraft configurations

Numerical simulation of compressible viscous flow fields is performed for a transonic transport configuration. A single structured grid system is constructed using analytical transformations such as conformal mapping, shearing/twisting/rotating/clustering/stretching transformations. The Reynolds-averaged, thin-layer Navier-Stokes equations are solved on a supercomputer, FACOM VP-400, using the LU-ADI factorization method.

Takanashi, Susumu↗
NASA NTRS1991

Navier-Stokes Simulation Of Wind-Tunnel Flow

Report describes three-dimensional Navier-Stokes simulation of transonic flow past model wing in wind tunnel and in free air. Numerical integration performed by LU-ADI factorization algorithm. Simulation shows effects of walls of tunnel.

Obayashi, S.↗
NASA NTRS1991

Navier-Stokes computation of wing/rotor interaction for a tilt rotor in hover

A method has been developed to analyze the wing/rotor interaction of tilt rotor aircraft in hover. The unsteady, thin-layer compressible Navier-Stokes equations are solved using an implicit, finite difference scheme that employs LU-ADI factorization. The rotor is modeled as an actuator disk which imparts a radial and azimuthal distribution of pressure rise and swirl to the flowfield. The 'chimera' approach of grid point blanking is used to update the rotor boundary conditions. Results are presented for both a rotor alone and for wing/rotor interaction where the thrust coefficient is 0.0164 and wing flap deflection is 67 degrees. Many of the complex flow features are captured including the fountain effect, leading and trailing edge separation, and the unsteady wake beneath the wing. Wing surface pressures compare fairly well with experimental data although the time-averaged download is about twenty percent higher than the measured value. This discrepancy is due to a combination of factors that are discussed.

Fejtek, Ian↗
NASA NTRS1993

Computational analysis of methods for reduction of induced drag

The purpose of this effort was to perform a computational flow analysis of a design concept centered around induced drag reduction and tip-vortex energy recovery. The flow model solves the unsteady three-dimensional Euler equations, discretized as a finite-volume method, utilizing a high-resolution approximate Riemann solver for cell interface flux definitions. The numerical scheme is an approximately-factored block LU implicit Newton iterative-refinement method. Multiblock domain decomposition is used to partition the field into an ordered arrangement of blocks. Three configurations are analyzed: a baseline fuselage-wing, a fuselage-wing-nacelle, and a fuselage-wing-nacelle-propfan. Aerodynamic force coefficients, propfan performance coefficients, and flowfield maps are used to qualitatively access design efficacy. Where appropriate, comparisons are made with available experimental data.

Janus, J. M.↗
NASA NTRS1993

Navier-Stokes flowfield computation of wing/rotor interaction for a tilt rotor aircraft in hover

The download on the wing produced by the rotor-induced downwash of a tilt rotor aircraft in hover is of major concern because of its severe impact on payload-carrying capability. A method has been developed to help gain a better understanding of the fundamental fluid dynamics that causes this download, and to help find ways to reduce it. In particular, the method is employed in this work to analyze the effect of a tangential leading edge circulation-control jet on download reduction. Because of the complexities associated with modeling the complete configuration, this work focuses specifically on the wing/rotor interaction of a tilt rotor aircraft in hover. The three-dimensional, unsteady, thin-layer compressible Navier-Stokes equations are solved using a time-accurate, implicit, finite difference scheme that employs LU-ADI factorization. The rotor is modeled as an actuator disk which imparts both a radical and an azimuthal distribution of pressure rise and swirl to the flowfield. A momentum theory blade element analysis of the rotor is incorporated into the Navier-Stokes solution method. Solution blanking at interior points of the mesh has been shown here to be an effective technique in introducing the effects of the rotor and tangential leading edge jet. Results are presented both for a rotor alone and for wing/rotor interaction. The overall mean characteristics of the rotor flowfield are computed including the flow acceleration through the rotor disk, the axial and swirl velocities in the rotor downwash, and the slipstream contraction. Many of the complex tilt rotor flow features are captured including the highly three-dimensional flow over the wing, the recirculation fountain at the plane of symmetry, wing leading and trailing edge separation, and the large region of separated flow beneath the wing. Mean wing surface pressures compare fairly well with available experimental data, but the time-averaged download/thrust ratio is 20-30 percent higher than the measured value. The discrepancy is due to a combination of factors that are discussed. Leading edge tangential blowing, of constant strength along the wing span, is shown to be effective in reducing download. The jet serves primarily to reduce the pressure on the wing upper surface. The computation clearly shows that, because of the three-dimensionality of the flowfield, optimum blowing would involve a spanwise variation in blowing strength.

Fejtek, Ian G.↗
NASA NTRS1993

CFD analyses of combustor and nozzle flowfields

The objectives of the research are to improve design capabilities for low thrust rocket engines through understanding of the detailed mixing and combustion processes. A Computational Fluid Dynamic (CFD) technique is employed to model the flowfields within the combustor, nozzle, and near plume field. The computational modeling of the rocket engine flowfields requires the application of the complete Navier-Stokes equations, coupled with species diffusion equations. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an ongoing experimental program at NASA LeRC. The numerical procedure is performed on both time-marching and time-accurate algorithms, using an LU approximate factorization in time, flux split upwinding differencing in space. The integrity of fuel film cooling along the wall, its effectiveness in the mixing with the core flow including unsteady large scale effects, the resultant impact on performance and the assessment of the near plume flow expansion to finite pressure altitude chamber are addressed.

Tsuei, Hsin-Hua↗
NASA NTRS1994

Numerical study of the trailing vortex of a wing with wing-tip blowing

Trailing vortices generated by lifting surfaces such as helicopter rotor blades, ship propellers, fixed wings, and canard control surfaces are known to be the source of noise, vibration, cavitation, degradation of performance, and other hazardous problems. Controlling these vortices is, therefore, of practical interest. The formation and behavior of the trailing vortices are studied in the present research. In addition, wing-tip blowing concepts employing axial blowing and spanwise blowing are studied to determine their effectiveness in controlling these vortices and their effects on the performance of the wing. The 3D, unsteady, thin-layer compressible Navier-Stokes equations are solved using a time-accurate, implicit, finite difference scheme that employs LU-ADI factorization. The wing-tip blowing is simulated using the actuator plane concept, thereby, not requiring resolution of the jet slot geometry. Furthermore, the solution blanking feature of the chimera scheme is used to simplify the parametric study procedure for the wing-tip blowing. Computed results are shown to compare favorably with experimental measurements. It is found that axial wing-tip blowing, although delaying the rolling-up of the trailing vortices and the near-field behavior of the flowfield, does not dissipate the circulation strength of the trailing vortex farther downstream. Spanwise wing-tip blowing has the effect of displacing the trailing vortices outboard and upward. The increased 'wing-span' due to the spanwise wing-tip blowing has the effect of lift augmentation on the wing and the strengthening of the trailing vortices. Secondary trailing vortices are created at high spanwise wing-tip blowing intensities.

Lim, Hock-Bin↗
NASA NTRS1992

Electromagnetic scattering calculations on the Intel Touchstone Delta

During the first year's operation of the Intel Touchstone Delta system, software which solves the electric field integral equations for fields scattered from arbitrarily shaped objects has been transferred to the Delta. To fully realize the Delta's resources, an out-of-core dense matrix solution algorithm that utilizes some or all of the 90 Gbyte of concurrent file system (CFS) has been used. The largest calculation completed to date computes the fields scattered from a perfectly conducting sphere modeled by 48,672 unknown functions, resulting in a complex valued dense matrix needing 37.9 Gbyte of storage. The out-of-core LU matrix factorization algorithm was executed in 8.25 h at a rate of 10.35 Gflops. Total time to complete the calculation was 19.7 h-the additional time was used to compute the 48,672 x 48,672 matrix entries, solve the system for a given excitation, and compute observable quantities. The calculation was performed in 64-b precision.

Cwik, Tom↗
NASA NTRS1994

Multi-dimensional combustor flowfield analyses in gas-gas rocket engine

The objectives of the present research are to improve design capabilities for low thrust rocket engines through understanding of the detailed mixing and combustions processes. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an on-going experimental program at NASA LeRC. Detailed computational modeling requires the application of the full three-dimensional Navier Stokes equations, coupled with species diffusion equations. The numerical procedure is performed on both time-marching and time-accurate algorithms and using an LU approximate factorization in time, flux split upwinding differencing in space. The emphasis in this paper is focused on using numerical analysis to understand detailed combustor flowfields, including the shear layer dynamics created between fuel film cooling and the core gas in the vicinity on the nearby combustor wall; the integrity and effectiveness of the coolant film; three-dimensional fuel jets injection/mixing/combustion characteristics; and their impacts on global engine performance.

Tsuei, Hsin-Hua↗
NASA NTRS1995

Spectral ordering techniques for incomplete LU preconditoners for CG methods

The effectiveness of an incomplete LU (ILU) factorization as a preconditioner for the conjugate gradient method can be highly dependent on the ordering of the matrix rows during its creation. Detailed justification for two heuristics commonly used in matrix ordering for anisotropic problems is given. The bandwidth reduction and weak connection following heuristics are implemented through an ordering method based on eigenvector computations. This spectral ordering is shown to be a good representation of the heuristics. Analysis and test cases in two and three dimensional diffusion problems demonstrate when ordering is important, and when an ILU decomposition will be ordering insensitive. The applicability of the heuristics is thus evaluated and placed on a more rigorous footing.

Clift, Simon S.↗
NASA NTRS1996

NASA flow fields analysis

The objectives of the present research are to improve design capabilities for low thrust rocket engines through understanding the detailed mixing and combustion processes in a representative combustor. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an on-going experimental program at NASA LERC. Detailed computational modeling involves the solution of both the two- and three-dimensional Navier Stokes equations, coupled with chemical reactions and the species diffusion equations. Computations of interest include both steady state and time-accurate flowfields and are obtained by means of LU approximate factorization in time and flux split upwinding differencing in space. The emphasis in the research is focused on using numerical analysis to understand detailed combustor flowfields, including the shear layer dynamics created between fuel film cooling and the core gas in the vicinity on the nearby combustor wall; the integrity and effectiveness of the coolant film; and three-dimensional fuel and oxidizer jet injection/mixing/combustion characteristics in the primary combustor along with their joint impacts on global engine performance.

Merkle, Charles L.↗
NASA NTRS1988

The ijk forms of factorization methods. II - Parallel systems

The paper considers the 'ijk forms' of LU and Cholesky factorization on certain parallel computers. This extends an earlier analysis for vector computers. Attention is restricted to local memory systems with processors that may or may not have vector capability. Special attention is given to bus architectures but qualitative analyses are given for other interconnection systems.

Ortega, J. M.↗
NASA NTRS1988

The ijk forms of factorization methods. I - Vector computers

This paper gives a detailed exposition of the 'ijk forms' of LU and Choleski factorization. Several aspects of these different organizations are discussed and their properties on vector computers are compared. Extensions of the ijk formalism to other algorithms is also given.

Ortega, J. M.↗
NASA NTRS1990

Characteristic-based algorithms for flows in thermo-chemical nonequilibrium

A generalized finite-rate chemistry algorithm with Steger-Warming, Van Leer, and Roe characteristic-based flux splittings is presented in three-dimensional generalized coordinates for the Navier-Stokes equations. Attention is placed on convergence to steady-state solutions with fully coupled chemistry. Time integration schemes including explicit m-stage Runge-Kutta, implicit approximate-factorization, relaxation and LU decomposition are investigated and compared in terms of residual reduction per unit of CPU time. Practical issues such as code vectorization and memory usage on modern supercomputers are discussed.

Walters, Robert W.↗
DOE OSTI2021

High-throughput characterization of Lu-doped zirconia

Lu-doped zirconia (LuDZ) was reported in the 1960s to be a poor oxygen ion conductor. However, we identified this system as worthy of reexamination based on the results of a machine learning-based screen for promising oxygen ion conductors, further supported by a physical argument: because Lu is the heaviest and smallest of the lanthanides, it should offer the possibility of high conductivity due to its high polarizability and size compatibility with the zirconia structure. Here the relationship between dopant concentration and conductivity in the LuDZ system was determined using a high throughput experimental methodology in which spatially resolved conductivity was measured across a compositionally graded film, spanning from Zr 0.93 Lu 0.07 O 1.965 to Zr 0.82 Lu 0.18 O 1.91 . The conductivity broadly peaks with composition at 15 cation % Lu, coinciding with the approximate phase boundary of the fully stabilized cubic phase. The composition-dependent conductivity measurements reveal a surprising increase in activation energy for ion transport with increasing Lu content with no change in this trend at the composition of the conductivity maximum. The pre-exponential factor increases sharply with Lu content, but displays a distinct plateau once the conductivity maximum composition is attained. At 500 °C, the conductivity of Zr 0.85 Lu 0.15 O 1.925 is 2 × 10 –3 S cm –1 , approaching that of optimally doped Sc-stabilized zirconia, which has the highest reported conductivity among zirconia compounds. In comparison to other oxide ion conductors, Lu-stabilized zirconia offers the advantage over ceria of negligible electronic conductivity and is easier to process than (La,Sr)(Ga,Mg)O 3-δ .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A continuous analogue of the tensor-train decomposition

We develop new approximation algorithms and data structures for representing and computing with multivariate functions using the functional tensor-train (FT), a continuous extension of the tensor-train (TT) decomposition. The FT represents functions using a tensor-train ansatz by replacing the three-dimensional TT cores with univariate matrix-valued functions. The main contribution of this paper is a framework to compute the FT that employs adaptive approximations of univariate fibers, and that is not connected to any tensorized discretization. The algorithm can be coupled with any univariate linear or nonlinear approximation procedure. We demonstrate that this approach can generate multivariate function approximations that are several orders of magnitude more accurate, for the same cost, than those based on the conventional approach of compressing the coefficient tensor of a tensor-product basis. Our approach is in the spirit of other continuous computation packages such as Chebfun, and yields an algorithm which requires the computation of “continuous” matrix factorizations such as the LU and QR decompositions of vector-valued functions. To support these developments, we describe continuous versions of an approximate maximum-volume cross approximation algorithm and of a rounding algorithm that re-approximates an FT by one of lower ranks. We demonstrate that our technique improves accuracy and robustness, compared to TT and quantics-TT approaches with fixed parameterizations, of high-dimensional integration, differentiation, and approximation of functions with local features such as discontinuities and other nonlinearities.

42 ENGINEERING↗
DOE OSTI2018

Lifetime-ever pneumonia among pre-school children across China – Associations with pre-natal and post-natal early life environmental factors

Highlights: • Childhood pneumonia can be influenced by prenatal and perinatal indoor exposure. • Mould and window condensation at birth is a risk factor for childhood pneumonia. • Prenatal exposure to redecoration and new furniture can increase the risk. • Cockroaches, mosquitos and flies are other risk factors in the home environment. • Farm exposure, breastfeeding and suburban or rural environment can be protective. Biomass burning is associated with childhood pneumonia but few studies exist on other indoor factors.

54 ENVIRONMENTAL SCIENCES↗