SEARCH · Engineering Papers

Results for “matrix approximation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Efficient Implementation of the Invariant Imbedding T-Matrix Method and the Separation of Variables Method Applied to Large Nonspherical Inhomogeneous Particles

Three terms, ''Waterman's T-matrix method'', ''extended boundary condition method (EBCM)'', and ''null field method'', have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIMþSOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIMþSOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIMþSOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIMþSOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging.

Invariant imbedding method↗

DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Disordered hydrogen adsorption at the three-fold site of W(110)

Hydrogen often forms disordered phases on metals; however, deciphering the atomic-scale behavior requires experimental and modeling techniques that directly probe the short-range order between neighboring hydrogen adsorbates. Here, to address this challenge, we applied direct recoil spectroscopy (DRS) to investigate hydrogen adsorption on the W(110) surface. We show that the recoiled hydrogen flux measured during DRS is sensitive to the short range order of the adsorbed hydrogen. Ordered hydrogen neighbors are much more efficient at dechanneling incident ions along low-index surface channels, strongly suppressing the recoiled hydrogen flux along these directions. By modeling the DRS measurements with molecular dynamics for ordered and disordered hydrogen phases, we find that at a hydrogen surface coverage of approximately Θ = 0.6, the hydrogen adsorbates are bound to three-fold sites in a disordered phase at room temperature. This finding is consistent with transfer-matrix scaling model calculations that suggest the transition to an ordered hydrogen phase occurs at a greater surface coverage of Θ = 0.83.

08 HYDROGEN↗

DOE OSTI2026

Quantum Monte Carlo calculation of {delta}C in the superallowed beta decay of 10C

We perform an ab initio quantum Monte Carlo calculation of the isospin-symmetry-breaking correction δC to the superallowed β decay of ¹⁰C. Using both phenomenological and chiral nuclear interactions, we evaluate the Fermi matrix element and quantify its deviation from the canonical √2 value. The resulting δC values lie in the range ≈ 0.15–0.25% and are consistent, within sizable uncertainties (approximately 34%–65% relative), across Hamiltonians, indicating no statistically significant dependence on the choice of nuclear interaction. The extracted values of Vud are also found to be compatible with current determinations within these uncertainties.

Piarulli, M↗

DOE OSTI2023

Online Convex Optimization of Programmable Quantum Computers to Simulate Time-Varying Quantum Channels

Simulating quantum channels is a fundamental primitive in quantum computing, since quantum channels define general (trace-preserving) quantum operations. An arbitrary quantum channel cannot be exactly simulated using a finite-dimensional programmable quantum processor, making it important to develop optimal approximate simulation techniques. In this paper, we study the challenging setting in which the channel to be simulated varies adversarially with time. We propose the use of matrix exponentiated gradient descent (MEGD), an online convex optimization method, and analytically show that it achieves a sublinear regret in time. Through experiments, we validate the main results for time-varying dephasing channels using a programmable generalized teleportation processor.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1975

Multidimensional Hermite-Gaussian quadrature formulae and their application to nonlinear estimation

A simplified technique is proposed for calculating multidimensional Hermite-Gaussian quadratures that involves taking the square root of a matrix by the Cholesky algorithm rather than computation of the eigenvectors of the matrix. Ways of reducing the dimension, number, and order of the quadratures are set forth. If the function f(x) under the integral sign is not well approximated by a low-order algebraic expression, the order of the quadrature may be reduced by factoring f(x) into an expression that is nearly algebraic and one that is Gaussian.

Mcreynolds, S. R.↗

NASA NTRS1977

Comment on the relation between the nonadiabatic coupling and the complex intersection of potential energy curves

Simple relations are discussed that provide a correspondence between the complex intersection of two potential surfaces and the nonadiabatic coupling matrix element between those surfaces. These are key quantities in semiclassical and quantum-mechanical theories of collision-induced electronic transitions. Within the two-state approximation, the complex intersection is shown to be directly related to the location and magnitude of the peak in the nonadiabatic coupling. Two cases have been considered: (1) the avoided crossing between two potential surfaces, and (2) the spin-orbit interaction due to a 2P halogen atom. Comparisons are made between the results of the two-state model and the results of ab initio quantum chemical calculations.

Jaffe, R. L.↗

NASA NTRS1978

Is photon angular momentum important in molecular collision processes occurring in a laser field

The importance of the rigorous treatment of photon angular momentum in molecular-collision processes occurring in the presence of intense radiation is investigated. An alternate approximate treatment, which essentially neglects the angular momentum coupling between the photon and the molecular degrees of freedom by averaging over the angular dependence of the interaction matrix elements, is presented and applied to a model calculation. The degeneracy-averaged results of this calculation compare remarkably well with the results of a rigorous calculation, from which we conclude (with reservation) that the explicit consideration of photoangular momentum coupling in molecular-collision problems is unnecessary.

Devries, P. L.↗

NASA NTRS1982

On the application and extension of Harten's high resolution scheme

Extensions of a second order high resolution explicit method for the numerical computation of weak solutions of one dimensonal hyperbolic conservation laws are discussed. The main objectives were (1) to examine the shock resoluton of Harten's method for a two dimensional shock reflection problem, (2) to study the use of a high resolution scheme as a post-processor to an approximate steady state solution, and (3) to construct an implicit in the delta-form using Harten's scheme for the explicit operator and a simplified iteration matrix for the implicit operator.

Yee, H. C.↗

NASA NTRS1988

Engineering applications of heuristic multilevel optimization methods

Some engineering applications of heuristic multilevel optimization methods are presented and the discussion focuses on the dependency matrix that indicates the relationship between problem functions and variables. Coordination of the subproblem optimizations is shown to be typically achieved through the use of exact or approximate sensitivity analysis. Areas for further development are identified.

Barthelemy, Jean-Francois M.↗

NASA NTRS1989

Engineering applications of heuristic multilevel optimization methods

Barthelemy, Jean-Francois M.↗

NASA NTRS1989

Engineering applications of heuristic multilevel optimization methods

Barthelemy, Jean-Francois M.↗

NASA NTRS1956

Extrapolation techniques applied to matrix methods in neutron diffusion problems

A general matrix method is developed for the solution of characteristic-value problems of the type arising in many physical applications. The scheme employed is essentially that of Gauss and Seidel with appropriate modifications needed to make it applicable to characteristic-value problems. An iterative procedure produces a sequence of estimates to the answer; and extrapolation techniques, based upon previous behavior of iterants, are utilized in speeding convergence. Theoretically sound limits are placed on the magnitude of the extrapolation that may be tolerated. This matrix method is applied to the problem of finding criticality and neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron fluxes in a nuclear reactor with control rods. The two-dimensional finite-difference approximation to the two-group neutron-diffusion equations is treated. Results for this example are indicated.

Mccready, Robert R↗

NASA NTRS1995

Kinetics of diffusional droplet growth in a liquid/liquid two-phase system

A new powerful experimental technique based on holographic observations, developed at the NASA Marshall Space Flight Center, now permits observation of small liquid droplets coarsening. This technique was developed and used for mixed-dimensional coarsening studies. Experiments were conducted on an isopycnic two-phase alloy of succinonitrile and water, annealed isothermally over a four-month period. The succinonitrile-rich droplets precipitate from a water-rich liquid matrix having a density very close to that of the droplets. The matrix and droplets, however, have different optical indices. The results of these experiments, along with the results of computer simulation based on the quasi-static diffusion approximation developed at Rensselaer are reported. These results were published recently. Copies of these papers are attached to this report.

Glickman, M. E.↗

NASA NTRS1995

A far-wing line shape theory which satisfies the detailed balance principle

A far-wing theory in which the validity of the detailed balance principle is maintained in each step of the derivation is presented. The role of the total density matrix including the initial correlations is analyzed rigorously. By factoring out the rapidly varying terms in the complex-time development operator in the interaction representation, better approximate expressions can be obtained. As a result, the spectral density can be expressed in terms of the line-coupling functions in which two coupled lines are arranged symmetrically and whose frequency detunings are omega - 1/2(omega(sub ji) + omega (sub j'i'). Using the approximate values omega - omega(sub ji) results in expressions that do not satisfy the detailed balance principle. However, this principle remains satisfied for the symmetrized spectral density in which not only the coupled lines are arranged symmetrically, but also the initial and final states belonging to the same lines are arranged symmetrically as well.

Ma, Q.↗

NASA NTRS2008

Non-Boltzmann Modeling for Air Shock-Layer Radiation at Lunar-Return Conditions

This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular electronic states present in lunar-return shock-layers. The Master Equation is derived for a general atom or molecule while accounting for a variety of excitation and de-excitation mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N2+, which are the main radiating species for most lunar-return shock-layers. Based on these new rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic and molecular states is developed. This new approach provides a simple and accurate method for calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix inversion procedure required for the detailed solution of the Master Equation. The radiative flux values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting approach are shown to agree within 5% for the Fire 1634 s case.

Johnston, Christopher O.↗

NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗

NASA NTRS2009

Application of Reduced Order Transonic Aerodynamic Influence Coefficient Matrix for Design Optimization

Supporting the Aeronautics Research Mission Directorate guidelines, the National Aeronautics and Space Administration [NASA] Dryden Flight Research Center is developing a multidisciplinary design, analysis, and optimization [MDAO] tool. This tool will leverage existing tools and practices, and allow the easy integration and adoption of new state-of-the-art software. Today s modern aircraft designs in transonic speed are a challenging task due to the computation time required for the unsteady aeroelastic analysis using a Computational Fluid Dynamics [CFD] code. Design approaches in this speed regime are mainly based on the manual trial and error. Because of the time required for unsteady CFD computations in time-domain, this will considerably slow down the whole design process. These analyses are usually performed repeatedly to optimize the final design. As a result, there is considerable motivation to be able to perform aeroelastic calculations more quickly and inexpensively. This paper will describe the development of unsteady transonic aeroelastic design methodology for design optimization using reduced modeling method and unsteady aerodynamic approximation. The method requires the unsteady transonic aerodynamics be represented in the frequency or Laplace domain. Dynamically linear assumption is used for creating Aerodynamic Influence Coefficient [AIC] matrices in transonic speed regime. Unsteady CFD computations are needed for the important columns of an AIC matrix which corresponded to the primary modes for the flutter. Order reduction techniques, such as Guyan reduction and improved reduction system, are used to reduce the size of problem transonic flutter can be found by the classic methods, such as Rational function approximation, p-k, p, root-locus etc. Such a methodology could be incorporated into MDAO tool for design optimization at a reasonable computational cost. The proposed technique is verified using the Aerostructures Test Wing 2 actually designed, built, and tested at NASA Dryden Flight Research Center. The results from the full order model and the approximate reduced order model are analyzed and compared.

Pak, Chan-gi↗