Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization methods”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗
NASA NTRS1990

High-performance equation solvers and their impact on finite element analysis

The role of equation solvers in modern structural analysis software is described. Direct and iterative equation solvers which exploit vectorization on modern high-performance computer systems are described and compared. The direct solvers are two Cholesky factorization methods. The first method utilizes a novel variable-band data storage format to achieve very high computation rates and the second method uses a sparse data storage format designed to reduce the number of operations. The iterative solvers are preconditioned conjugate gradient methods. Two different preconditioners are included; the first uses a diagonal matrix storage scheme to achieve high computation rates and the second requires a sparse data storage scheme and converges to the solution in fewer iterations that the first. The impact of using all of the equation solvers in a common structural analysis software system is demonstrated by solving several representative structural analysis problems.

Poole, Eugene L.↗
NASA NTRS1991

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form were iterated in time until a steady solution was reached. Evidence was given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at the domain boundaries was proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects were accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. The flow in a developing S-duct was then solved in the laminar regime in a Reynolds number (Re) of 790 and in the turbulent regime at Re equals 40,000 by using the Baldwin-Lomax model. The Stanitz elbow was then solved by using an invicid version of the same code at M sub inlet equals 0.4. Grid dependence and convergence rate were investigated, showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re equals 2.5 times 10(exp 6) was solved with the Baldwin-Lomax and the q-omega models. Both approaches show satisfactory agreement with experiments, although the q-omega model was slightly more accurate.

Michelassi, V.↗
NASA NTRS1990

A finite element solver for 3-D compressible viscous flows

Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.

Reddy, K. C.↗
NASA NTRS1990

Simplified procedures for designing composite bolted joints

Simplified procedures are described to design and analyze single and multi-bolt composite joints. Numerical examples illustrate the use of these methods. Factors affecting composite bolted joints are summarized. References are cited where more detailed discussion is presented on specific aspects of composite bolted joints. Design variables associated with these joints are summarized in the appendix.

Chamis, Christos C.↗
NASA NTRS1991

Elliptic flow computation by low Reynolds number two-equation turbulence models

A detailed comparison of ten low-Reynolds-number k-epsilon models is carried out. The flow solver, based on an implicit approximate factorization method, is designed for incompressible, steady two-dimensional flows. The conservation of mass is enforced by the artificial compressibility approach and the computational domain is discretized using centered finite differences. The turbulence model predictions of the flow past a hill are compared with experiments at Re = 10 exp 6. The effects of the grid spacing together with the numerical efficiency of the various formulations are investigated. The results show that the models provide a satisfactory prediction of the flow field in the presence of a favorable pressure gradient, while the accuracy rapidly deteriorates when a strong adverse pressure gradient is encountered. A newly proposed model form that does not explicitly depend on the wall distance seems promising for application to complex geometries.

Michelassi, V.↗
NASA NTRS1990

A user guide for the EMTAC-MZ CFD code

The computer code (EMTAC-MZ) was applied to investigate the flow field over a variety of very complex three-dimensional (3-D) configurations across the Mach number range (subsonic, transonic, supersonic, and hypersonic flow). In the code, a finite volume, multizone implementation of high accuracy, total variation diminishing (TVD) formulation (based on Roe's scheme) is used to solve the unsteady Euler equations. In the supersonic regions of the flow, an infinitely large time step and a space-marching scheme is employed. A finite time step and a relaxation or 3-D approximate factorization method is used in subsonic flow regions. The multizone technique allows very complicated configurations to be modeled without geometry modifications, and can easily handle combined internal and external flow problems. An elliptic grid generation package is built into the EMTAC-MZ code. To generate the computational grid, only the surface geometry data are required. Results obtained for a variety of configurations, such as fighter-like configurations (F-14, AVSTOL), flow through inlet, multi-bodies (shuttle with external tank and SRBs), are reported and shown to be in good agreement with available experimental data.

Szema, Kuo-Yen↗
NASA NTRS1992

High-performance equation solvers and their impact on finite element analysis

The role of equation solvers in modern structural analysis software is described. Direct and iterative equation solvers which exploit vectorization on modern high-performance computer systems are described and compared. The direct solvers are two Cholesky factorization methods. The first method utilizes a novel variable-band data storage format to achieve very high computation rates and the second method uses a sparse data storage format designed to reduce the number od operations. The iterative solvers are preconditioned conjugate gradient methods. Two different preconditioners are included; the first uses a diagonal matrix storage scheme to achieve high computation rates and the second requires a sparse data storage scheme and converges to the solution in fewer iterations that the first. The impact of using all of the equation solvers in a common structural analysis software system is demonstrated by solving several representative structural analysis problems.

Poole, Eugene L.↗
NASA NTRS1992

Multivariable frequency domain identification via 2-norm minimization

The author develops a computational approach to multivariable frequency domain identification, based on 2-norm minimization. In particular, a Gauss-Newton (GN) iteration is developed to minimize the 2-norm of the error between frequency domain data and a matrix fraction transfer function estimate. To improve the global performance of the optimization algorithm, the GN iteration is initialized using the solution to a particular sequentially reweighted least squares problem, denoted as the SK iteration. The least squares problems which arise from both the SK and GN iterations are shown to involve sparse matrices with identical block structure. A sparse matrix QR factorization method is developed to exploit the special block structure, and to efficiently compute the least squares solution. A numerical example involving the identification of a multiple-input multiple-output (MIMO) plant having 286 unknown parameters is given to illustrate the effectiveness of the algorithm.

Bayard, David S.↗
NASA NTRS1997

Massively Parallel and Scalable Implicit Time Integration Algorithms for Structural Dynamics

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because of the following additional facts: (a) explicit schemes are easier to parallelize than implicit ones, and (b) explicit schemes induce short range interprocessor communications that are relatively inexpensive, while the factorization methods used in most implicit schemes induce long range interprocessor communications that often ruin the sought-after speed-up. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet be offset by the speed of the currently available parallel hardware. Therefore, it is essential to develop efficient alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating the low-frequency dynamics of aerospace structures.

Farhat, Charbel↗
NASA NTRS2004

Improvements in Electron-Probe Microanalysis: Applications to Terrestrial, Extraterrestrial, and Space-Grown Materials

Improvement in the accuracy of electron-probe microanalysis (EPMA) has been accomplished by critical assessment of standards, correction algorithms, and mass absorption coefficient data sets. Experimental measurement of relative x-ray intensities at multiple accelerating potential highlights errors in the absorption coefficient. The factor method has been applied to the evaluation of systematic errors in the analysis of semiconductor and silicate minds. Accurate EPMA of Martian soil stimulant is necessary in studies that build on Martian rover data in anticipation of missions to Mars.

Carpenter, Paul↗
NASA NTRS2011

Assessing Uncertainties in Boundary Layer Transition Predictions for HIFiRE-1 at Non-zero Angles of Attack

Boundary layer stability was analyzed for the HIFiRE-1 flight vehicle geometry for ground tests conducted at the CUBRC LENS I hypersonic shock test facility and the Langley Research Center (LaRC) 20- inch Mach 6 Tunnel. Boundary layer stability results were compared to transition onset location obtained from discrete heat transfer measurements from thin film gauges during the CUBRC test and spatially continuous heat transfer measurements from thermal phosphor paint data during the LaRC test. The focus of this analysis was on conditions at non-zero angles of attack as stability analysis has already been performed at zero degrees angle of attack. Also, the transition onset data obtained during flight testing was at nonzero angles of attack, so this analysis could be expanded in the future to include the results of the flight test data. Stability analysis was performed using the 2D parabolized stability software suite STABL (Stability and Transition Analysis for Hypersonic Boundary Layers) developed at the University of Minnesota and the mean flow solutions were computed using the DPLR finite volume Navier-Stokes computational fluid dynamics (CFD) solver. A center line slice of the 3D mean flow solution was used for the stability analysis to incorporate the angle of attack effects while still taking advantage of the 2D STABL software suite. The N-factors at transition onset and the value of Re(sub theta)/M(sub e), commonly used to predict boundary layer transition onset, were compared for all conditions analyzed. Ground test data was analyzed at Mach 7.2 and Mach 6.0 and angles of attack of 1deg, 3deg and 5deg. At these conditions, the flow was found to be second mode dominant for the HIFiRE-1 slender cone geometry. On the leeward side of the vehicle, a strong trend of transition onset location with angle of attack was observed as the boundary layer on the leeward side of the vehicle developed inflection points at streamwise positions on the vehicle that correlated to angle of attack. Inflection points are a strong instability mechanism that lead to rapid breakdown and transition to turbulence. The transition onset location on the windward side of the vehicle displayed no trend with angle of attack or freestream Reynolds number and transition was observed farther down the vehicle than observed on the leeward side of the vehicle. In analysis of both windward and leeward sides of the vehicle, use of the N factor methodology to develop trends to predict boundary layer transition onset showed improvements over the Re(sub theta)/M(sub e) empirical correlation methodology. Stronger correlations and less scatter in the data were observed when using the N factor method for these cases.

Marek, Lindsay C.↗
NASA NTRS2017

Metabolomic and Genomic Markers of Atherosclerosis as Related to Oxidative Stress, Inflammation, and Vascular Function in Twin Astronauts

BACKGROUND: Future human space travel will consist primarily of long-duration missions onboard the International Space Station (ISS) or exploration-class missions to Mars, its moons, or nearby asteroids. Astronauts participating in long-duration missions may be at an increased risk of oxidative stress and inflammatory damage due to radiation, psychological stress, altered physical activity, nutritional insufficiency, and hyperoxia during extravehicular activity. By studying one identical twin during his 1-year ISS mission and his ground-based twin, this work extends a current NASA-funded investigation to determine whether these spaceflight factors contribute to an accelerated progression of atherosclerosis. This study of twins affords a unique opportunity to examine spaceflight-related atherosclerosis risk that is independent of the confounding factors associated with different genotypes. PURPOSE: The purpose of this investigation was to determine whether biomarkers of oxidative and inflammatory stress are elevated during and after long-duration spaceflight and determine if a relation exists between levels of these biomarkers and structural and functional indices of atherosclerotic risk measured in the carotid and brachial arteries. These physiological and biochemical data will be extended by using an exploratory approach to investigate the relationship between intermediate phenotypes and risk factors for atherosclerosis and the metabolomic signature from plasma and urine samples. Since metabolites are often the indirect products of gene expression, we simultaneously assessed gene expression and DNA methylation in leukocytes. HYPOTHESIS: We predict that, compared to the ground-based twin, the space-flown twin will experience elevated biomarkers of oxidative stress and inflammatory damage, altered arterial structure and function, accelerated telomere shortening, dysregulation of genes associated with oxidative stress and inflammation, and a metabolic profile shift that is associated with elevated atherosclerosis risk factors. METHODS: In the space-flown twin, a panel of biomarkers of oxidative and inflammatory stress were measured in venous blood samples and in 24-h (in-flight) and 48-h (pre- and post-flight) urine pools collected twice before flight, six times during the mission (~FD15, 75, 180, 240, 300, 335), and early in the post-flight recovery phase (3-5 days after landing). We also measured metabolomic (targeted and untargeted approaches) and genomic markers (DNA methylation, mRNA gene expression, telomere length) in these samples. Arterial structure, assessed from measures of intima-media thickness, also were measured using standard clinical ultrasound at the same time points. Arterial function was assessed using brachial flow-mediated dilation, a well-validated measure used to assess endothelium-dependent vasodilation and a sensitive predictor of atherosclerotic risk, only before and after spaceflight. All of the same measures were obtained in the ground-based twin, but less frequently. DISCUSSION: All data collection has been completed for both the space-flown twin and the ground-based twin. Vascular structure and function measures have been analyzed, blood and urine samples have been batch-processed. Results from these individuals will be compared to each other, to data from other Twin Study investigations, and to the larger complement of subjects participating in the companion study currently ongoing in ISS astronauts.

Lee, Stuart M. C.↗
NASA NTRS2016

On the Distance to the North Polar Spur and the Local CO-H2 Factor

Aims. Most models identify the X-ray bright North Polar Spur (NPS) with a hot interstellar (IS) bubble in the Sco-Cen star-forming region at approx. = to 130 pc. An opposite view considers the NPS as a distant structure associated with Galactic nuclear outflows. Constraints on the NPS distance can be obtained by comparing the foreground IS gas column inferred from X-ray absorption to the distribution of gas and dust along the line of sight. Absorbing columns towards shadowing molecular clouds simultaneously constrain the CO-H2 conversion factor. Methods. We derived the columns of X-ray absorbing matter N(sub Habs) from spectral fitting of dedicated XMM-Newton observations towards the NPS southern terminus (l(sup II) approx. = to 29 degrees, b(sup II) approx. = to +5 to +11 degrees). The distribution of the IS matter was obtained from absorption lines in stellar spectra, 3D dust maps and emission data, including high spatial resolution CO measurements recorded for this purpose. Results. N(sub Habs) varies from approx. = to 4.3 to approx. = to 1.3 x 10(exp 21) per sq cm along the 19 fields. Relationships between X-ray brightness, absorbing column and hardness ratio demonstrate a brightness decrease with latitude governed by increasing absorption. The comparison with absorption data, local and large-scale dust maps rules out a NPS near side closer than 300 pc. The correlation between N(sub Habs) and the reddening increases with the sightline length from 300 pc to 4 kpc and is the tightest with Planck tau(sub 353GHz) -based reddening, suggesting a much larger distance. N(H)/E(B-V)(sub tao) approx. = to 4.1 x 10(exp 21) per sq cm mag, close to Fermi-Planck determinations. N(sub Habs) absolute values are compatible with HI-CO clouds at -5 less than or equal to V(sub LSR) less than or equal to +25 to +45 km/s and a NPS potentially far beyond the Local Arm. A shadow cast by a b=+9 degree molecular cloud constrains X(sub CO) in that direction to less than or equal to 1.0 x 10(exp 20) sec/((sq cm)K km). The average X(sub CO) over the fields is less than or equal to 0.75 x 10(exp 20) sec/((sq cm)K km).

X-rays: ISM / radio lines: ISM / local insterstell↗
NASA NTRS2020

Automated Noise Calibration System (VT-1000)

This paper details an automated Noise Source calibration system in development at Jet Propulsion Laboratory, California Institute of Technology (JPL). The paper begins with a discussion on Noise Figure and Excess Noise Ratio (ENR) theory, fundamentals and governing equations. As part of the fundamentals is a discussion of the system’s use of the Y-factor method to obtain accurate measurements of the Unit Under Test (UUT), and how these measurements are compared against a known ENR standard to obtain the UUT’s ENR values. There is also an in-depth discussion on uncertainty quantification for Noise Source system calibrations. The architecture of the automated calibration system is provided, which includes both the system’s hardware and software configuration. The software is written in Python 3, and provides the user detailed instruction on how to proceed, including step-by-step connection requirements. This system automates much of the measurement process, including real-time uncertainty quantification and report generation, as well as real-time feedback to the user to allow intervention if necessary. The system takes advantage of a database of results from previous measurements to compare calibration history of the ENR measurements. The automated system presented here operates over a frequency range from 10 MHz to 50 GHz, and has shown substantial time savings over traditional manual methods of performing this calibration.

Timpe, Scott↗
NASA NTRS1985

Technique for measuring gas conversion factors

A method for determining hydrocarbon conversion factors for a flowmeter. A mixture of air, O2 and C sub x H sub y is burned and the partial paressure of O2 in the resulting gas is forced to equal the partial pressure of O2 in air. The flowrate of O2 flowing into the mixture is measured by flowmeter and the flowrate of C sub x H sub y flowing into the mixture is measured by the flowmeter conversion factor is to be determined. These measured values are used to calculate the conversion factor.

Singh, J. J.↗
NASA NTRS2018

Method and Apparatus for Spacecraft Gyroscope Scale Factor Calibration

A method and apparatus for estimating gyro scale factor during normal spacecraft operations, using small attitude motions that are compliant with mission pointing accuracy and stability requirements and a signal processing method that specifically detects the intentionally induced motions. This process increases operational availability by avoiding the need to take the spacecraft offline for large calibration maneuvers.

Reth, Alan D↗
NASA NTRS1974

Geodesic synchrotron radiation in the Kerr geometry by the method of asymptotically factorized Green's functions

The scalar, electromagnetic, and gravitational geodesic-synchrotron-radiation (GSR) spectra are determined for the case of a test particle moving on a highly relativistic circular orbit about a rotating (Kerr) black hole. It is found that the spectral shape depends only weakly on the value of the angular-momentum parameter (a/M) of the black hole, but the total radiated power drops unexpectedly for a value of at least 0.95 and vanishes as the value approaches unity. A spin-dependent factor (involving the inner product of the polarization of a radiated quantum with the source) is isolated to explain the dependence of the spectral shape on the spin of the radiated field. Although the scalar wave equation is solved by separation of variables, this procedure is avoided for the vector and tensor cases by postulating a sum-over-states expansion for the Green's function similar to that found to hold in the scalar case. The terms in this sum, significant for GSR, can then be evaluated in the geometric-optics approximation without requiring the use of vector or tensor spherical harmonics.

Chrzanowski, P. L.↗