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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Single Phase Passive Rectification Versus Active Rectification Applied to High Power Stirling Engines

Stirling engine converters are being considered as potential candidates for high power energy conversion systems required by future NASA explorations missions. These types of engines typically contain two major moving parts, the displacer and the piston, in which a linear alternator is attached to the piston to produce a single phase sinusoidal waveform at a specific electric frequency. Since all Stirling engines perform at low electrical frequencies (less or equal to 100 Hz), space explorations missions that will employ these engines will be required to use DC power management and distribution (PMAD) system instead of an AC PMAD system to save on space and weight. Therefore, to supply such DC power an AC to DC converter is connected to the Stirling engine. There are two types of AC to DC converters that can be employed, a passive full bridge diode rectifier and an active switching full bridge rectifier. Due to the inherent line inductance of the Stirling Engine-Linear Alternator (SE-LA), their sinusoidal voltage and current will be phase shifted producing a power factor below 1. In order to keep power the factor close to unity, both AC to DC converters topologies will implement power factor correction. This paper discusses these power factor correction methods as well as their impact on overall mass for exploration applications. Simulation results on both AC to DC converters topologies with power factor correction as a function of output power and SE-LA line inductance impedance are presented and compared.

Santiago, Walter↗
NASA NTRS1978

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗
DOE OSTI2022

Case studies in bias reduction and inference for electronic health record data with selection bias and phenotype misclassification

Electronic health records (EHR) are not designed for population‐based research, but they provide easy and quick access to longitudinal health information for a large number of individuals. Many statistical methods have been proposed to account for selection bias, missing data, phenotyping errors, or other problems that arise in EHR data analysis. However, addressing multiple sources of bias simultaneously is challenging. We developed a methodological framework (R package, SAMBA ) for jointly handling both selection bias and phenotype misclassification in the EHR setting that leverages external data sources. These methods assume factors related to selection and misclassification are fully observed, but these factors may be poorly understood and partially observed in practice. As a follow‐up to the methodological work, we demonstrate how to apply these methods for two real‐world case studies, and we evaluate their performance. In both examples, we use individual patient‐level data collected through the University of Michigan Health System and various external population‐based data sources. In case study (a), we explore the impact of these methods on estimated associations between gender and cancer diagnosis. In case study (b), we compare corrected associations between previously identified genetic loci and age‐related macular degeneration with gold standard external summary estimates. These case studies illustrate how to utilize diverse auxiliary information to achieve less biased inference in EHR‐based research.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the 'loop unrolling' technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large-scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the loop unrolling technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗
DOE OSTI2023

Liquid solution centrifugation for safe, scalable, and efficient isotope separation

A general method of separating isotopes by centrifuging dissolved chemical compounds in a liquid is introduced. This technique can be applied to almost all elements and leads to large separation factors. The method has been demonstrated in several isotopic systems including Ca, Mo, O, and Li with single-stage selectivities of 1.046 to 1.067 per neutron mass difference (e.g., 1.43 in 40 Ca/ 48 Ca), which are beyond the capabilities of various conventional methods. Equations are derived to model the process, and the results agree with those of the experiments. The scalability of the technique has been demonstrated by a three-stage enrichment of 48 Ca with a total 40 Ca/ 48 Ca selectivity of 2.43, and the scalability is more broadly supported through analogies to gas centrifuge, whereby countercurrent centrifugation can further multiply the separation factor by 5 to 10 times per stage in a continuous process. Optimal centrifuge conditions and solutions can achieve both high-throughput and highly efficient isotope separation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Analysis of Hydrocyclone Geometry via Rapid Optimization Based on Computational Fluid Dynamics

Abstract Hydrocyclones exploit density gradients for the centrifugal separation of dispersions in a continuous liquid. Selection of the geometrics for optimal separation is case specific, like the media characteristics. The existing optimization method based on computational fluid dynamics (CFD) provides a powerful analytical tool but requires long computational times. The most common praxis for CFD optimization is via the single‐factor optimization method (SFOM). In this study, a novel approach is presented as an improved rapid optimization method that implements a dynamic‐mesh and user‐defined function optimization method (DUOM). The DUOM adapts the dynamic‐mesh approach from other applications to the optimization analysis of hydrocyclones. The DUOM reduced the computational time by 31.1 %, compared to the SFOM.

Liu, Lin↗
DOE OSTI2020

Determination of acoustic nonlinearity parameters using thermal modulation of ultrasonic waves

This study presents a test method and its theoretical framework to determine the acoustic nonlinearity parameters (α, β, δ) of material using thermal modulation of ultrasonic waves. Temperature change-induced thermal strain excites the nonlinear response of the material and modulates the ultrasonic wave propagating in it. Experimental results showed a strong correlation between the relative wave velocity change and the temperature change. With a quadratic polynomial model, the acoustic nonlinearity parameters were obtained from the polynomial coefficients by curve fitting the experimental data. Their effects on thermal-induced velocity change were discussed. The parameters α, β, and δ govern the hysteretic gap, average slope, and curvature of the correlation curve, respectively. The proposed theory was validated on aluminum, steel, intact and damaged concrete samples. The obtained nonlinear parameters show reasonable agreement with values reported in the literature. Compared to other nonlinear acoustic methods using vibration or acoustic excitation, the thermal modulation method generates more uniform, slow changing, and larger strain field in the test sample. Employing the thermal effect as the driving force for nonlinearity instead of an undesired influencing factor, this method can measure the absolute values of α, β, and δ with good accuracy using a simple ultrasonic test setup.

Physics↗
NASA NTRS1976

Data reduction and evaluation procedures

The computational procedures that are involved in exhaust emissions data reduction and the use of these computational procedures for determining the quality of the data that is obtained from exhaust measurements were considered. Four problem areas were calculated: (1) the various methods for performing the carbon balance, (2) the method for calculating water correction factors, (3) the method for calculating the exhaust molecular weight, and (4) assessing the quality of the data.

Mirsky, W.↗
NASA NTRS1992

Active strut placement using integer programming for the CSI Revolutionary Model

A method for determining the most effective locations for active struts on large space structures is developed and tested on the NASA CSI Evolutionary Model. Depending on the choice of weighting factors, the method can be used to maximize the maximum modal damping ratio or decay rate, or to maximize a mission-oriented measure of performance. Placement of 8, 16, and 32 active struts out of 1507 candidate truss elements is demonstrated. Preliminary estimates of damping enhancement are reported pending refined structural models and dynamic test results. The method can handle complicated FEM models wih a large number of truss elements and many target modes. It can treat each mode equally or it can emphasize the importance of selected modes. The method can eliminate some combinations of actuator locations based on topological constraints.

Padula, Sharon L.↗
NASA NTRS1996

Biases In A Magnitude Limited Versus A Distant Limited Planet Search

All methods utilized to conduct a search or survey inevitably have some built-in biases. These biases are often traced to some limitation of the instrument used or some inherent character of the signal being recorded. We address these limitations for various methods used or proposed for planet detection: spectroscopy, astrometry, interferometry, and photometry. For spectroscopy, the turbulence in the photosphere limits the minimum measurable dossier velocity to 3 m/s and hence the minimum planet to star mass ratio, thereby favoring massive close-in planets. Limited available observing time will necessarily introduce additional selection biases in the targets observed and telescope aperture will limit the faintest magnitude stars to be measured; For astrometry, the angular resolution of the instrument along with motions in the photometric center of the star limit the furthest distance for which giant planets can be detected to about 10 pc and favor massive outer orbit planets around low mass non-solar like stars; For imaging interferometry, the minimum angular size of the central null limits both the distance to the star and closeness of the planet to the host star to about 1 AU at 10 pc and thus to the very few solar-like stars within 10 pc. Solar and extra-solar zodiacal emission will limit the minimum size of the detectable planet. For photometry, the inherent variability of the star does not limit the minimum planet size until earth-sized or smaller planets are considered around solar-like stars. The telescope aperture limits the faintest stars that can be monitored. As with spectroscopy, there is no inherent distant limit to the method. After addressing the limiting factors of each method, an estimate is made of the number of planets of various sizes that could be found for each stellar type based on the detection probability and the number of stars that can be searched for planets.

Koch, David↗
NASA NTRS1988

A fast and accurate method for evaluating the nonrelativistic free-free Gaunt factor for hydrogenic ions

A fast method for obtaining a two-dimensional Chebyshev expansion of the nonrelativistic hydrogenic free-free Gaunt factor is presented which gives a maximum relative error of 0.7 percent over a wide range of temperatures and frequencies. The expansion derived by Karzas and Latter (1961) is used to derive primary values of the free-free Gaunt factor as a function of the initial and final electron energies. Thermal averages have been numerically computed and are presented in tabular form in intervals of 0.25 dex.

Hummer, D. G.↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Implicit methods in CFD

A class of implicit approximate factorization schemes is examined for stability and convergence characteristics. These schemes include Newton's method, factorization, and flux-vector splitting. Examples are used to show that all practical methods suffer from some limited stability or asymptotic convergence restriction. A three-dimensional factored scheme is shown which suffers from unconditional instability which can only be ameliorated by added artificial dissipation. An F3D + or - flux split scheme is described which avoids unconditional instability, but in the end has similar convergence characteristics.

Pulliam, Thomas H.↗
NASA NTRS1992

Viscous simulation method for unsteady flows past multicomponent configurations

The present algorithm for the numerical simulation of flows about complex configurations (whose multiple and nonsimilar components have arbitraty geometries) employs a hybridization of the domain decomposition techniques for grid generation as well as to reduce computer-memory requirements. A fully vectorized, finite-volume, upwind-biased, approximately factored multigrid method is used to solve 3D Reynolds-averaged unsteady and compressible Navier-Stokes equations simulating supersonic flows past an ogive-nose-cylinder near or within a cavity. The time-averaged surface pressures obtained compare favorably with the wind tunnel data.

Fouladi, Kamran↗
NASA NTRS1984

Supersonic nonlinear potential analysis

The NCOREL computer code was established to compute supersonic flow fields of wings and bodies. The method encompasses an implicit finite difference transonic relaxation method to solve the full potential equation in a spherical coordinate system. Two basic topic to broaden the applicability and usefulness of the present method which is encompassed within the computer code NCOREL for the treatment of supersonic flow problems were studied. The first topic is that of computing efficiency. Accelerated schemes are in use for transonic flow problems. One such scheme is the approximate factorization (AF) method and an AF scheme to the supersonic flow problem is developed. The second topic is the computation of wake flows. The proper modeling of wake flows is important for multicomponent configurations such as wing-body and multiple lifting surfaces where the wake of one lifting surface has a pronounced effect on a downstream body or other lifting surfaces.

Siclari, M. J.↗
NASA NTRS1986

Upwind relaxation algorithms for Euler/Navier-Stokes equations

A description of and results from a solution algorithm for the compressible Navier-Stokes equations are presented. The main features of the algorithm are second or third order accurate upwind discretization of the convection and pressure derivatives and a relaxation scheme for the unfactored implicit backward Euler time method, implemented in a finite-volume formulation. Upwind methods were successfully used to obtain solutions to the Euler equations for flows with strong shock waves. The particular upwind method being used is based on the flux vector splitting technique developed by Van Leer and both second and third order accurate discretizations were developed. Currently, the most widely used implicit solution technique for the Navier-Stokes equations use approximate factorization (AF) methods to treat multidimensional problems. The time integration scheme being used in the present algorithm corresponds to a line Gauss-Seidel relaxation method. This method produces good convergence rates for steady-state flows, and most of the algorithm was vectorized on the NASA Langley VPS 32 computer. The Navier-Stokes algorithm was tested for several two-dimensional flow problems. Solutions for the problems gave excellent results. The presented effort is directed toward the extension of the scheme to the full three-dimensional Navier-Stokes equations.

Thomas, J. L.↗
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗