Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 55 records · Page 3

Extension of a streamwise upwind algorithm to a moving grid system

A new streamwise upwind algorithm was derived to compute unsteady flow fields with the use of a moving-grid system. The temporally nonconservative LU-ADI (lower-upper-factored, alternating-direction-implicit) method was applied for time marching computations. A comparison of the temporally nonconservative method with a time-conservative implicit upwind method indicates that the solutions are insensitive to the conservative properties of the implicit solvers when practical time steps are used. Using this new method, computations were made for an oscillating wing at a transonic Mach number. The computed results confirm that the present upwind scheme captures the shock motion better than the central-difference scheme based on the beam-warming algorithm. The new upwind option of the code allows larger time-steps and thus is more efficient, even though it requires slightly more computational time per time step than the central-difference option.

Obayashi, Shigeru↗

Convergence acceleration of an aeroelastic Navier-Stokes solver

New capabilities have been added to a Navier-Stokes solver to perform steady-state simulations more efficiently. The flow solver for solving the Navier-Stokes equations is completely rewritten with a combination of the LU-SGS (Lower-Upper factored Symmetric Gauss-Seidel) implicit method and the modified HLLE (Harten-Lax-van Leer-Einfeldt) upwind scheme. A pseudo-time marching method is used for the directly coupled structural equations to improve overall convergence rates for static aeroelastic analysis. Results are demonstrated for transonic flows over rigid and flexible wings.

Obayashi, S.↗

Numerical Simulations of an Unsteady Rocket Launch from the AH-64D Apache Longbow Helicopter

Rocket and missile firings from attack helicopters can cause main engine compressor stall. Studies of this phenomenon suggest that the main engine ingests either the plume from the rockets or the rocket blast waves. This creates surges at the inlet face, causing a loss of power in the main engine. The objective of this project is to set-up a computational fluid dynamics (CFD) simulation of the AH-64D Apache Longbow helicopter during a rocket launch, in order to qualitatively study the fluid dynamics of the problem. This project presents a progression of three unsteady Navier-Stokes solutions. The first unsteady solution involves only a rocket launch from its launch canister. The second solution is a launch from a canister mounted on the Apache's wing-pylon assembly. The last solution includes the Apache main engine and fuselage. The computations use a series of structured, overset grid systems, which allow for a rocket moving in a prescribed path. The method implements a Roe upwind scheme with LU-SGS (lower-upper factored symmetric Gauss-Seidel). A rotor pressure disk model approximates the helicopter rotor, while the rocket engine exit properties are applied as a prescribed boundary condition. Although the project is only at the half-way point, the first and second CFD simulations suggest the possibility of pressure wave interference. Sudden surges in pressure occur from two sources: at rocket start-up, and as the rocket leaves the canister. Wave patterns set-up by these sources appear to propagate to the location of the engine inlet. However the simplified geometry simulation with the main engine needs to be performed before coming to a conclusion.

Okamoto, Kevin↗

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗

An implict LU scheme for the Euler equations applied to arbitrary cascades

An implicit scheme for solving the Euler equations is derived and demonstrated. The alternating-direction implicit (ADI) technique is modified, using two implicit-operator factors corresponding to lower-block-diagonal (L) or upper-block-diagonal (U) algebraic systems which can be easily inverted. The resulting LU scheme is implemented in finite-volume mode and applied to 2D subsonic and transonic cascade flows with differing degrees of geometric complexity. The results are presented graphically and found to be in good agreement with those of other numerical and analytical approaches. The LU method is also 2.0-3.4 times faster than ADI, suggesting its value in calculating 3D problems.

Buratynski, E. K.↗

Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a tubine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of independence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Choo, Yung K.↗

Application of a lower-upper implicit scheme and an interactive grid generation for turbomachinery flow field simulations

A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a turbine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of dependence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Choo, Yung K.↗

Transient Gene and miRNA Expression Profile Changes of Confluent Human Fibroblast Cells in Space

Microgravity or an altered gravity environment from the static 1 gravitational constant has been shown to influence global gene expression patterns and protein levels in cultured cells. However, most of the reported studies conducted in space or using simulated microgravity on the ground have focused on the growth or differentiation of the cells. Whether non-dividing cultured cells will sense the presence of microgravity in space has not been specifically addressed. In an experiment conducted on the International Space Station, confluent human fibroblast cells were fixed after being cultured in space for 3 and 14 days for investigations of gene and miRNA (microRNA) expression profile changes in these cells. A fibroblast is a type of cell that synthesizes the extracellular matrix and collagen, the structural framework for tissues, and plays a critical role in wound healing and other functions. Results of the experiment showed that on Day 3, both the flown and ground cells were still proliferating slowly even though they were confluent, as measured by the expression of the protein Ki-67 positive cells, and the cells in space grew slightly faster. Gene and miRNA expression data indicated activation of NF(sub kappa)B (nuclear factor kappa-light-chain-enhancer of activated B cells) and other growth related pathways involving HGF and VEGF in the flown cells. On Day 14 when the cells were mostly non-dividing, the gene and miRNA expression profiles between the flight and ground samples were indistinguishable. Comparison of gene and miRNA expressions in the Day 3 samples in respect to Day 14 revealed that most of the changes observed on Day 3 were related to cell growth for both the flown and ground cells. Analysis of cytoskeleton changes by immunohistochemistry staining of the cells with antibodies for alpha-tubulin showed no difference between the flight and ground samples. Results of our study suggest that in true non-dividing human fibroblast cells, microgravity in space has little effect on the gene and miRNA expression. Gene and miRNA expression changes were observed in cells that were confluent, but still proliferating slowly. The faster growth in the flown cells was associated with the activation of NF(sub kappa)B pathways which triggers the expression of several growth factors and the suppression of the cell cycle checkpoint.

Zhang, Ye↗

Modeling the interaction of ultrasound with pores

Factors that affect ultrasonic velocity sensing of density during consolidation of metal powders are examined. A comparison is made between experimental results obtained during the final stage of densification and the predictions of models that assume either a spherical or a spheroidal pore shape. It is found that for measurements made at low frequencies during the final stage of densification, relative density (pore fraction) and pore shape are the two most important factors determining the ultrasonic velocity, the effect of pore size is negligible.

Lu, Yichi↗

Using Strassen's algorithm to accelerate the solution of linear systems

Strassen's algorithm for fast matrix-matrix multiplication has been implemented for matrices of arbitrary shapes on the CRAY-2 and CRAY Y-MP supercomputers. Several techniques have been used to reduce the scratch space requirement for this algorithm while simultaneously preserving a high level of performance. When the resulting Strassen-based matrix multiply routine is combined with some routines from the new LAPACK library, LU decomposition can be performed with rates significantly higher than those achieved by conventional means. We succeeded in factoring a 2048 x 2048 matrix on the CRAY Y-MP at a rate equivalent to 325 MFLOPS.

Bailey, David H.↗

Observation of Solar Wind Charge Exchange Emission From Exospheric Material in and Outside Earth's Magnetosheath 2008 September 25

A long XMM-Newton exposure is used to observe solar wind charge exchange (SWCX) emission from exospheric material in and outside Earth's magnetosheath. The light curve of the O vii (0.5-0.62 keV) band is compared with a model for the expected emission, and while the emission is faint and the light curve has considerable scatter, the correlation is significant to better than 99.9%. This result demonstrates the validity of the geocoronal SWCX emission model for predicting a contribution to astrophysical observations to a scale factor of order unity (1.5). In addition, an average value of the SWCX O vii emission from the magnetosheath over the observation of 2.6 +/- 0.5 LU is derived. The results also demonstrate the potential utility of using X-ray observations to study global phenomena of the magnetosheath which currently are only investigated using in situ measurements.

Snowden, S. L.↗

Formation of the muonic helium atom /alpha particle-muon-electron/ and observation of its Larmor precession

Experiments are described in which it proved possible to form the muonic helium atom by stopping polarized negative muons in a helium gas with a 2% xenon admixture at a pressure of 14 atm. The observed Larmor precession amplitudes are plotted against the gyromagnetic ratio for both muons and antimuons stopped in He + 2% Xe. In addition, a non-zero residual polarization of 0.06 plus or minus 0.01 was measured for muons stopped in pure helium gas, which corresponds to a depolarization factor of 18 plus or minus 3.

Souder, P. A.↗

TGF-beta1 release from biodegradable polymer microparticles: its effects on marrow stromal osteoblast function

BACKGROUND: Controlled release of transforming growth factor-beta1 (TGF-beta1) to a bone defect may be beneficial for the induction of a bone regeneration cascade. The objectives of this work were to assess the feasibility of using biodegradable polymer microparticles as carriers for controlled TGF-beta1 delivery and the effects of released TGF-beta1 on the proliferation and differentiation of marrow stromal cells in vitro. METHODS: Recombinant human TGF-beta1 was incorporated into microparticles of blends of poly(DL-lactic-co-glycolic acid) (PLGA) and poly(ethylene glycol) (PEG). Fluorescein isothiocynate-labeled bovine serum albumin (FITC-BSA) was co-encapsulated as a porogen. The effects of PEG content (0, 1, or 5% by weight [wt%]) and buffer pH (3, 5, or 7.4) on the protein release kinetics and the degradation of PLGA were determined in vitro for as long as 28 days. Rat marrow stromal cells were seeded on a biodegradable poly(propylene fumarate) (PPF) substrate. The dose response and biological activity of released TGF-beta1 was determined after 3 days in culture. The effects of TGF-beta1 released from PLGA/PEG microparticles on marrow stromal cell proliferation and osteoblastic differentiation were assessed during a 21-day period. RESULTS: TGF-beta1 was encapsulated along with FITC-BSA into PLGA/PEG blend microparticles and released in a multiphasic fashion including an initial burst for as long as 28 days in vitro. Increasing the initial PEG content resulted in a decreased cumulative mass of released proteins. Aggregation of FITC-BSA occurred at lower buffer pH, which led to decreased release rates of both proteins. The degradation of PLGA was increased at higher PEG content and significantly accelerated at acidic pH conditions. Rat marrow stromal cells cultured on PPF substrates showed a dose response to TGF-beta1 released from the microparticles similar to that of added TGF-beta1, indicating that the activity of TGF-beta1 was retained during microparticle fabrication and after growth factor release. At an optimal TGF-beta1 dosage of 1.0 ng/ml after 3 days, the released TGF-beta1 enhanced the proliferation and osteoblastic differentiation of marrow stromal cells over 21 days of culture, with increased total cell number, alkaline phosphatase activity, and osteocalcin production. CONCLUSIONS: PLGA/PEG blend microparticles can serve as delivery vehicles for controlled release of TGF-beta1, and the released growth factor enhances marrow stromal cell proliferation and osteoblastic differentiation in vitro. CLINICAL RELEVANCE: Controlled release of TGF-beta1 from PLGA/PEG microparticles is representative of emerging tissue engineering technologies that may modulate cellular responses to encourage bone regeneration at a skeletal defect site.

Non-NASA Center↗

Preconditioned conjugate gradient methods for the Navier-Stokes equations

A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulation. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a conventional implicit solver, namely line Gauss-Seidel relaxation (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU time on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow case are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner.

Ajmani, Kumud↗

Numerical Procedures for Inlet/Diffuser/Nozzle Flows

Two primitive variable, pressure based, flux-split, RNS/NS solution procedures for viscous flows are presented. Both methods are uniformly valid across the full Mach number range, Le., from the incompressible limit to high supersonic speeds. The first method is an 'optimized' version of a previously developed global pressure relaxation RNS procedure. Considerable reduction in the number of relatively expensive matrix inversion, and thereby in the computational time, has been achieved with this procedure. CPU times are reduced by a factor of 15 for predominantly elliptic flows (incompressible and low subsonic). The second method is a time-marching, 'linearized' convection RNS/NS procedure. The key to the efficiency of this procedure is the reduction to a single LU inversion at the inflow cross-plane. The remainder of the algorithm simply requires back-substitution with this LU and the corresponding residual vector at any cross-plane location. This method is not time-consistent, but has a convective-type CFL stability limitation. Both formulations are robust and provide accurate solutions for a variety of internal viscous flows to be provided herein.

Rubin, Stanley G.↗

Relative Accuracy Assessment of LANDSAT-4 MSS and TM Data for Level 1 Land Cover Inventory

Digital data for the Washington, DC scene simultaneously acquired by the LANDSAT-4 Multispectral Scanner (MSS) and the LANDSAT-4 thematic mapper (TM) was compared. Classification success for the TM and MSS data sets was determined by a per pixel comparison with digitized ground verification data (GVD). These GVD were comprised of Level 7 land cover (developed, agriculture, forest, water, wetlands, and barren) for four USGS 7.5 minute topographic quadrangle maps. The relative improvement in classification success for TM was between 11% and 14%, or about a factor of 1.3, for these data. This represents a meaningful improvement in accuracy for Level 7 land cover categorization for TM relative to MSS, particularly when errors of omission and commission were considered.

Middleton, E. M.↗

Emission Factors for Crop Residue and Prescribed Fires in the Eastern US during FIREX-AQ

Agricultural and prescribed burning activities emit large amounts of trace gases and aerosols on regional to global scales. We present a compilation of emission factors (EFs) and emission ratios (ERs) from the eastern portion of the Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) campaign in 2019 in the United States, which sampled burning of crop residues and other prescribed fire fuels. FIREX-AQ provided comprehensive chemical characterization of 53 crop residue and 22 prescribed fires. Crop residues burned at different modified combustion efficiencies (MCE), with corn residue burning at higher MCE than other fuel types. Prescribed fires burned at lower MCE (<0.90) which is typical, while grasslands burned at lower MCE (0.90) than normally observed due to moist, green, growing season fuels. Most non-methane volatile organic compounds (NMVOCs) were significantly anticorrelated with MCE except for ethanol and NMVOCs that were measured with less certainty. We identified 23 species where crop residue fires differed by more than 50% from prescribed fires at the same MCE. Crop residue EFs were greater for species related to agricultural chemical use and fuel composition as well as oxygenated NMVOCs possibly due to the presence of metals such as potassium. Prescribed EFs were greater for monoterpenes (5×). FIREX-AQ crop residue average EFs generally agreed with the previous agricultural fire study in the US but had large disagreements with global compilations. FIREX-AQ observations show the importance of regionally-specific and fuel-specific EFs as first steps to reduce uncertainty in modeling the air quality impacts of fire emissions.

air quality↗

Orientation of Space Station Freedom electrical power system in environmental effects assessment

The orientation effects of six Space Station Freedom Electrical Power System (EPS) components are evaluated for three environmental interactions: aerodynamic drag, atomic oxygen erosion, and orbital debris impact. Designers can directly apply these orientation factors to estimate the magnitude of the examined environment and the environmental effects for the EPS component of interest. The six EPS components are the solar array, photovoltaic module radiator, integrated equipment assembly, solar dynamic concentrator, solar dynamic radiator, and beta gimbal.

Lu, Cheng-Yi↗