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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Using Whispering-Gallery-Mode Resonators for Refractometry

A method of determining the refractive and absorptive properties of optically transparent materials involves a combination of theoretical and experimental analysis of electromagnetic responses of whispering-gallery-mode (WGM) resonator disks made of those materials. The method was conceived especially for use in studying transparent photorefractive materials, for which purpose this method affords unprecedented levels of sensitivity and accuracy. The method is expected to be particularly useful for measuring temporally varying refractive and absorptive properties of photorefractive materials at infrared wavelengths. Still more particularly, the method is expected to be useful for measuring drifts in these properties that are so slow that, heretofore, the properties were assumed to be constant. The basic idea of the method is to attempt to infer values of the photorefractive properties of a material by seeking to match (1) theoretical predictions of the spectral responses (or selected features thereof) of a WGM of known dimensions made of the material with (2) the actual spectral responses (or selected features thereof). Spectral features that are useful for this purpose include resonance frequencies, free spectral ranges (differences between resonance frequencies of adjacently numbered modes), and resonance quality factors (Q values). The method has been demonstrated in several experiments, one of which was performed on a WGM resonator made from a disk of LiNbO3 doped with 5 percent of MgO. The free spectral range of the resonator was approximately equal to 3.42 GHz at wavelengths in the vicinity of 780 nm, the smallest full width at half maximum of a mode was approximately equal to 50 MHz, and the thickness of the resonator in the area of mode localization was 30 microns. In the experiment, laser power of 9 mW was coupled into the resonator with an efficiency of 75 percent, and the laser was scanned over a frequency band 9 GHz wide at a nominal wavelength of approximately equal to 780 nm. Resonance frequencies were measured as functions of time during several hours of exposure to the laser light. The results of these measurements, plotted in the figure, show a pronounced collective frequency drift of the resonator modes. The size of the drift has been estimated to correspond to a change of 8.5 x 10(exp -5) in the effective ordinary index of refraction of the resonator material.

Matsko, Andrey↗
DOE OSTI2024

Time-dependent density functional theory with the orthogonal projector augmented wave method

The projector augmented wave (PAW) method of Blöchl linearly maps smooth pseudo wavefunctions to the highly oscillatory all-electron DFT orbitals. Compared to norm-conserving pseudopotentials (NCPP), PAW has the advantage of lower kinetic energy cutoffs and larger grid spacing at the cost of having to solve for non-orthogonal wavefunctions. We earlier developed orthogonal PAW (OPAW) to allow the use of PAW when orthogonal wavefunctions are required. In OPAW, the pseudo wavefunctions are transformed through the efficient application of powers of the PAW overlap operator with essentially no extra cost compared to NCPP methods. Previously, we applied OPAW to DFT. Here, we take the first step to make OPAW viable for post-DFT methods by implementing it in real-time time-dependent (TD) DFT. Using fourth-order Runge–Kutta for the time-propagation, we compare calculations of absorption spectra for various organic and biological molecules and show that very large grid spacings are sufficient, 0.6–0.7 bohr in OPAW-TDDFT rather than the 0.4–0.5 bohr used in traditional NCPP-TDDFT calculations. This reduces the memory and propagation costs by around a factor of 3. Our method would be directly applicable to any post-DFT methods that require time-dependent propagations such as the GW approximation and the Bethe–Salpeter equation.

Chemistry↗
DOE OSTI2025

Fermilab Automation of Coil Winding System

The automated coil winder project creates a repeatable, scalable method for producing superconducting magnet coils with increased accuracy and efficiency. The design uses an ortho cyclic winding method to achieve high fill factors by combining computer numerical control (CNC), microcontroller-based controls, and sophisticated optimization techniques. The method optimizes coil turns, core geometry, and current consumption by utilizing analytical models based on magnetic moment, Ohm’s law, and power constraints. For particle accelerator applications, prototype testing using a 3D-printed coil winder has shown enhanced winding repeatability, thermal management, and manufacturability.

Dzida, Rafal [Northern Illinois U.]↗
DOE OSTI2024

Gravitational form factors of charmonia

We investigate the gravitational form factors of charmonium. Our method is based on a Hamiltonian formalism on the light front known as basis light-front quantization. The charmonium mass spectrum and light-front wave functions were obtained from diagonalizing an effective Hamiltonian that incorporates confinement from holographic QCD and one-gluon exchange interaction from light-front QCD. We proposed a quantum many-body approach to construct the hadronic matrix elements of the energy momentum tensor T + + and T + − , which are used to extract the gravitational form factors A ( Q 2 ) and D ( Q 2 ) . The obtained form factors satisfy the known constraints, e.g., the von Laue condition. From these quantities, we also extract the energy, pressure and light-front energy distributions of the system. We find that hadrons are multilayer systems. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
NASA NTRS1991

Efficient multitasking of Choleski matrix factorization on CRAY supercomputers

A Choleski method is described and used to solve linear systems of equations that arise in large scale structural analysis. The method uses a novel variable-band storage scheme and is structured to exploit fast local memory caches while minimizing data access delays between main memory and vector registers. Several parallel implementations of this method are described for the CRAY-2 and CRAY Y-MP computers demonstrating the use of microtasking and autotasking directives. A portable parallel language, FORCE, is used for comparison with the microtasked and autotasked implementations. Results are presented comparing the matrix factorization times for three representative structural analysis problems from runs made in both dedicated and multi-user modes on both computers. CPU and wall clock timings are given for the parallel implementations and are compared to single processor timings of the same algorithm.

Overman, Andrea L.↗
NASA NTRS1978

Three-dimensional radiative exchange factors for nongray, nondiffuse surfaces

A computer program using a variation of the Monte Carlo method was developed to predict radiative exchange factors in three-dimensional configurations with curved surfaces. The model allows for arbitrary emission, absorption, and reflection characteristics and for radiation exchange between surfaces as well as external irradiation through openings. It was shown that exchange factors can be determined efficiently by using the normal Monte Carlo method for closed configurations and/or high surface emissivities, while energy partitioning is vastly superior in the case of open configurations and/or low surface emissivities. Comparison with some experiments performed on the Space Shuttle heat rejection system demonstrates the model's accuracy as well as its superiority over a simplified diffuse-surface analysis, which requires similar amounts of computer time.

Modest, M. F.↗
NASA NTRS1977

On the logarithmic-singularity correction in the kernel function method of subsonic lifting-surface theory

A logarithmic-singularity correction factor is derived for use in kernel function methods associated with Multhopp's subsonic lifting-surface theory. Because of the form of the factor, a relation was formulated between the numbers of chordwise and spanwise control points needed for good accuracy. This formulation is developed and discussed. Numerical results are given to show the improvement of the computation with the new correction factor.

Lan, C. E.↗
NASA NTRS1996

Multigrid and Krylov Subspace Methods for the Discrete Stokes Equations

Discretization of the Stokes equations produces a symmetric indefinite system of linear equations. For stable discretizations, a variety of numerical methods have been proposed that have rates of convergence independent of the mesh size used in the discretization. In this paper, we compare the performance of four such methods: variants of the Uzawa, preconditioned conjugate gradient, preconditioned conjugate residual, and multigrid methods, for solving several two-dimensional model problems. The results indicate that where it is applicable, multigrid with smoothing based on incomplete factorization is more efficient than the other methods, but typically by no more than a factor of two. The conjugate residual method has the advantage of being both independent of iteration parameters and widely applicable.

Elman, Howard C.↗
NASA NTRS1985

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗
NASA NTRS1974

The radiative effect of aerosols in the earth's atmosphere

A modified two-flux approximation is employed to compute the transfer of radiation in a finite, inhomogeneous, turbid atmosphere. A perturbation technique is developed to allow the treatment of nongray gaseous absorption with multiple scattering. The perturbation method, which employs a backscatter factor as a parameter, can be used with anisotropic particle scattering as well as Rayleigh scattering. This method is used to study the effect of aerosols on radiative solar heating and infrared cooling as well as the radiative-convective temperature distribution in the earth's atmosphere. It is found that the effect of aerosols in the infrared cannot be neglected; while in the visible, the effect can be of the same order as that due to absorption by water vapor. For a high surface albedo (greater than 0.30) heating of the earth-atmosphere system results due to the presence of aerosols. The aerosols also reduce the amount of convection needed to maintain a stable atmosphere. For the case of a dense haze a temperature inversion is found to exist close to the ground.

Wang, W.-C.↗
NASA NTRS1952

Radiant-interchange Configuration Factors

A study is presented of the geometric configuration factors required for computing radiant heat transfer between opaque surfaces separated by a nonabsorbing medium and various methods of determining the configuration factors are discussed. Configuration-factor solutions available in the literature have been checked and the more complicated equations are presented as families of curves. Cases for point, line, and finite-area sources are worked out over a wide range of geometric proportions. These cases include several new configurations involving rectangles, triangles, and cylinders of finite length which are integrated and tabulated. An analysis is presented, in which configuration factors are employed of the radiant heat transfer to the rotor blades of a typical gas turbine under different conditions of temperature and pressure. (author)

Hamilton, D C :↗
NASA NTRS1992

Organic crystal growth experiment facility (13-IML-1)

The interesting nature of metal-like organic compounds composed of charge transfer complexes has been recently realized. Crystals of these complexes can usually be grown by the solution crystallization method. It is difficult to grow such organic crystals on Earth, especially from the chemical reactions through diffusion controlled process in the solutions, because of gravitational disturbances, or sedimentation. The International Microgravity Lab. (IML-1) Organic Crystal Growth with G-Gitter Preventive Measure (OCGP) experiment is expected to grow a single crystal large enough to allow its intrinsic physical properties to be measured and its detailed crystal structure to be determined. This experiment also attempts to assess the experimental conditions including the microgravity environment for further study of the fundamental process of solution crystallization, nucleation, and growth from supersaturated phases including chemical reactions. Microgravity disturbances, G-jitter, may be an important environmental factor in the experimental method to assess. The vibration damping effects on organic crystal growth can be carefully studied.

Kanbayashi, Akio↗
NASA NTRS1975

Satellite animal tracking feasibility studies

A study was initiated in Tsavo National Park to determine movements and home ranges of individual elephants and their relations to overall distribution patterns and environmental factors such as rainfall. Methods used were radio tracking and observations of visually identifiable individuals. Aerial counts provided data on overall distribution. Two bulls and two cows were radio-tagged in Tsavo West and two bulls and four cows in Tsavo East, providing home range and movement data. The movements of individuals were useful in interpreting relatively major shifts in elephant distribution. Results point to the following preliminary conclusions: (1) elephants in the Tsavo area undertook long distance movements in fairly direct response to localized rainfall; (2) a subdivision of the overall population into locally distinct units may exist during the dry season but did not occur after significant rainfall; and (3) food appears to be the primary factor governing movements and distribution of elephants in the area.

Buechner, H. K.↗
DOE OSTI2020

Low synchronization Gram–Schmidt and generalized minimal residual algorithms

The Gram–Schmidt process uses orthogonal projection to construct the A = QR factorization of a matrix. When Q has linearly independent columns, the operator P = I - Q(QTQ)-1QT defines an orthogonal projection onto Q⊥. In finite precision, Q loses orthogonality as the factorization progresses. A family of approximate projections is derived with the form P = I - QTQT, with correction matrix T. When T = (QTQ)-1, and T is triangular, it is postulated that the best achievable orthogonality is $\mathcal{O}(ε)\mathcal{K}(A)$. We present new variants of modified (MGS) and classical Gram–Schmidt algorithms that require one global reduction step. An interesting form of the projector leads to a compact WY representation for MGS. In particular, the inverse compact WY MGS algorithm is equivalent to a lower triangular solve. Our main contribution is to introduce a backward normalization lag into the compact WY representation, resulting in a $\mathcal{O}(ε)\mathcal{K}[r_0, AV_m])$ stable Generalized Minimal Residual Method (GMRES) algorithm that requires only one global reduce per iteration. Finally, further improvements in performance are achieved by accelerating GMRES on GPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Westcott g factors extended to arbitrary neutron energy spectra

Westcott 𝑔 factors are used in Neutron Activation Analysis (NAA) and Prompt Gamma-ray Activation Analysis (PGAA) to evaluate the impact of non-1∕𝑣 behavior in the neutron-capture cross sections of certain nuclei on activation product yields. This non-1∕𝑣 behavior arises from the presence of neutron resonances in the neutron- capture cross sections that overlap with the source neutron spectrum at low (< 5 eV) energies. Historically, Westcott 𝑔 factors that have been cataloged for NAA and PGAA applications are the result of calculations that assume a Maxwellian neutron flux distribution with a given temperature. In this work, we use this approach with updated neutron-capture cross sections from the Evaluated Nuclear Data File, version VIII.1 (ENDF/B-VIII.1) to tabulate Westcott 𝑔 factor values for a broad range of Maxwellian distribution temperatures, comparing the results against currently-available 𝑔 factors from International Atomic Energy Agency tables and other sources. Here, it was discovered during this analysis that the use of guided thermal and cold-neutron beams at certain facilities necessitates an approach for evaluating Westcott 𝑔 factors based on arbitrary non-Maxwellian spectra. In this paper, we present an approach for calculating 𝑔 factors with user-specified neutron spectra, and we demonstrate these methods to obtain Westcott 𝑔-factors for guided- and cold-neutron beams at the Budapest Research Reactor and the Forschungsreaktor München II reactor. As part of this work, open-source software has been developed that can be used to perform these calculations for applications in PGAA and NAA experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Development of high throughput light-sheet fluorescence lifetime imaging microscopy for 3D functional imaging of metabolic pathways in plant and microorganisms (Final Technical Report)

This research program will enable new biochemical contrast in the nanosecond lifetime domain through use of the recently demonstrated electro-optic fluorescence lifetime imaging technique (EO-FLIM) for wide-field lifetime imaging. The Stanford/Stanford Linear Accelerator Center multidisciplinary collaboration -- physics, applied physics, and structural biology -- will develop a light-sheet fluorescence lifetime imaging microscope for functional studies of microbial and plant metabolic pathways and dynamic interactions between plants and microorganisms in the rhizosphere. The proposed approach overcomes the imaging time bottleneck associated with existing fluorescence lifetime imaging methods. Initial demonstrations have shown a factor of 100,000 improvement in photon throughput compared to existing methods. High photon efficiency allowed the first wide-field fluorescence lifetime imaging of single molecules. Recent work has improved the technique’s repetition rate to enable compatibility with mode-locked lasers and demonstrated the combination of wide-field fluorescence lifetime imaging with super-resolution localization microscopy, observations of single molecule dynamics, and observation of donor lifetime quenching in single-molecule imaging. These results were achieved on standard camera sensors and would not have been possible with other wide-field approaches. The throughput and photon economy of the EO-FLIM method enables new BER-relevant imaging opportunities. In particular, scanned single- and two-photon light-sheet excitation will be used to achieve volumetric imaging with time-domain contrast.

47 OTHER INSTRUMENTATION↗
NASA NTRS1987

Numerical study of laminar and turbulent flows inside a turnaround duct with and without guide vanes

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside a turnaround duct with and without guide vanes and to investigate the effects of vanes on the flow characteristics. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. Method of segmental sweeping is developed to overcome the multiple-zone problem due to imposition of guide vanes. In the present study, several test cases have been computed. These test cases include laminar and turbulent flows inside a turnaround duct without and with two or three guide vanes. The study reveals that: (1) there exists large recirculation zones inside the duct if no vanes are present; (2) properly shaped and positional guide vanes are effective in eliminating flow separation; and (3) laminar and turbulent flows have similar flow features. But turbulent flow has less total pressure drop.

Lin, S.-J.↗
NASA NTRS1987

A Navier-Stokes flow simulation of the Space Shuttle Main Engine Hot Gas Manifold

Incompressible viscous flow inside the turnaround duct, the fuel bowl, the transfer duct and the racetrack of the Space Shuttle Main Engine (SSME) Hot Gas Manifold (HGM) has been computed using the method of pseudo-compressibility together with an implicit, approximate-factorization algorithm. A multiple-zone method is used to make solution of flows in complex geometries easy. A model which predicts the pressure loading for the shield and the injector post arrangement without solving the complex flow field in the main injector region is proposed. The computed results show good qualitative agreement with experimental data.

Yang, Ruey-Jen↗