Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Theory”

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 1,441 records · Page 80

LSZ in action: extracting form factors from correlators nonperturbatively in 2d ϕ4 theory

Abstract In this paper, we compute multiparticle form factors of local operators in 2dϕ 4 theory using a recently proposed method [1] for efficiently implementing the LSZ prescription with Hamiltonian Truncation methods, and we adopt Lightcone Conformal Truncation (LCT) in particular for our calculations. We perform various checks of our results at weak and strong coupling, and elucidate the parametric behavior of truncation errors. This opens up the possibility to compute S-matrix in various strongly coupled models using the LSZ method in LCT.

Physics↗

Systematic Study of the Self-Renormalized Nucleon Gluon PDF in Large-Momentum Effective Theory

We present a systematic study of the nucleon gluon parton distribution function (PDF) using the self-renormalized large-momentum effective theory (LaMET) approach in lattice QCD. This work extends previous gluon-PDF extractions by performing a detailed analysis of key systematic effects, including gauge-link smearing, lattice spacing, pion mass, and nucleon boost momentum. The self-renormalization framework mitigates ultraviolet divergences associated with Wilson-line self-energy and renormalon contributions by combining lattice matrix elements with perturbative short-distance information, thereby preserving the correct infrared structure. Calculations are performed on $N_f=2+1+1$ HISQ ensembles generated by the MILC Collaboration at three lattice spacings and two pion masses, with boosted nucleon states reaching momenta up to 2.2~GeV. We determine renormalization factors from zero-momentum matrix elements and apply hybrid renormalization to suppress discretization artifacts. After extrapolating large-separation behavior and performing Fourier transforms, we reconstruct quasi-PDFs and match them to lightcone PDFs using next-to-leading order Wilson coefficients. Our results demonstrate that smearing and lattice-spacing effects are under control, and pion-mass and lattice-spacing dependence is mild relative to the current $O(10^6)$ statistics; however, momentum dependence remains a significant source of uncertainty. Future work including even larger boost momenta will be essential to reduce systematics in lattice determinations of the gluon PDF and to advance toward precision QCD phenomenology at the LHC and the future Electron-Ion Collider.

FOS: Physical sciences↗

Vibrational dependence of pressure induced spectral linewidths and line shifts - Application of the infinite order sudden scattering approximation

The infinite order sudden (IOS) approximation to molecular rotation is applied to simplify the theory of linewidths and shifts in vibration-rotation spectra. This approximation is expected to be most accurate for hard, short-range collisions and is therefore complementary to Anderson theory which is best for weak, glancing collisions. The IOS approximation predicts identical linewidths and shifts for P- and R-branch transitions with the same line number. It also predicts zero line shifts for pure rotational spectra. The dependence of linewidths and shifts on vibrational band is seen to be due mainly to variation in diagonal vibrational matrix elements of the intermolecular potential. Calculations are performed for the 0-0, 0-1, and 0-2 bands of CO perturbed by He, using a theoretical interaction potential with no semiempirical or adjustable parameters; results are in satisfactory accord with experimental data.

Green, S.↗

Incompressible lifting-surface aerodynamics for a rotor-stator combination

Current literature on the three dimensional flow through compressor cascades deals with a row of rotor blades in isolation. Since the distance between the rotor and stator is usually 10 to 20 percent of the blade chord, the aerodynamic interference between them has to be considered for a proper evaluation of the aerothermodynamic performance of the stage. A unified approach to the aerodynamics of the incompressible flow through a stage is presented that uses the lifting surface theory for a compressor cascade of arbitrary camber and thickness distribution. The effects of rotor stator interference are represented as a linear function of the rotor and stator flows separately. The loading distribution on the rotor and stator flows separately. The loading distribution on the rotor and stator blades and the interference factor are determined concurrently through a matrix iteration process.

Ramachandra, S. M.↗

Radar scene generation for tactical decision aids

The Mueller matrix and polarization covariance matrix for polarimetric radar systems was studied. The clutter is modeled by a layer of random permittivity, described by a three-dimensional correlation function, with variance, and horizontal and vertical correlation lengths. A general mixing formula was derived for discrete scatters immersed in a host medium. The results are applicable to general multiphase mixtures. The strong fluctuation theory was used to derive the backscattering cross sections, and was further extended to include higher order co-polarized and cross-polarized moments. A two-layer anisotropic random medium model was developed for the active and passive remote sensing of ice fields. A three-layer random medium model was adopted to study the volume scattering effects for the active and passive microwave remote sensing of snow-covered ice fields. The snow layer was simulated by an isotropic random medium and the ice layer by an anisotropic random medium. The vegetation canopy and snow-covered ice field were studied with a three-layer model, an isotropic random medium layer overlying an anisotropic random medium. The dyadic Green's function of the three-layer medium and the scattered electromagnetic intensities with Born approximation were calculated.

Kong, J. A.↗

The Surfaces of Pluto and Charon

Much of the surface of Pluto consists of high-albedo regions covered to an unknown depth by Beta-N2, contaminated with CH4, CO, and other molecules. A portion of the exposed surface appears to consist of solid H2O. The remainder is covered by lower albedo material of unknown composition. The N2 ice may occur as polar caps of large extent, leaving ices and other solids of lower volatility in the equatorial regions. The low-albedo material found primarily in the equatorial regions may consist in part of solid hydrocarbons and nitriles produced from N2 and CH4 in the atmosphere or in the surface ices. Alternatively, it may arise from deposition from impacting bodies and/or the chemistry of the impact process itself. Charon's surface is probably more compositionally uniform than that of Pluto, and is covered by H2O ice with possible contaminants or exposures of other materials that are as yet unidentified. The molecular ices discovered on Pluto and Charon have been identified from near-infrared spectra obtained with Earth-based telescopes. The quantitative interpretation of those data has been achieved through the computation of synthetic spectra using the Hapke scattering theory and the optical constants of various ices observed in the laboratory. Despite limitations imposed by the availability of laboratory data on ices in various mixtures, certain specific results have been obtained. It appears that CH4 and CO are trace constituents, and that some fraction of the CH4 (and probably the CO) on Pluto is dissolved in the matrix of solid N2. Pure CH4 probably also occurs on Pluto's surface, allowing direct access to the atmosphere. Study of the nitrogen absorption band at 2.148 micrometers shows that the temperature of the N2 in the present epoch is 40 +/-2 K. The global temperature regime of Pluto can be modeled from observations of the thermal flux at far-infrared and millimeter wavelengths. The low-albedo equatorial regions must be significantly warmer than the polar regions covered by N2 (at T = 40 K) to account for the total thermal flux measured. At the present season, the diurnal skin depth of the insolation-driven thermal wave is small, and the observed mm-wave fluxes may arise from a greater depth. Alternatively, the mm-wave flux may arise from the cool, sublimation source region. The surface microstructure in the regions covered by N2 ice is likely governed by the sintering properties of this highly volatile material. The observed nitrogen infrared band strength requires that expanses of the surface be covered with cm-sized crystals of N2. Grains of H2O ice on Charon, in contrast, are probably of order 50 micrometers in size, and do not metamorphose into larger grains at a significant rate. Because of the similarities in size, density, atmosphere and surface composition between Pluto and Neptune's satellite Triton, the surface structures observed by Voyager on Triton serve as a plausible paradigm for what might be expected on Pluto. Such crater forms, tectonic structures, aeolian features, cryovolcanic structures, and sublimation-degraded topography as are eventually observed on Pluto and Charon by spacecraft will give information on their interior compositions and structures, as well as on the temperature and wind regimes over the planet's extreme seasonal cycle.

Cruikshank, Dale P.↗

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

For scattering calculations involving nonspherical particles such as ice crystals, we show that the transverse wave condition is not applicable to the refracted electromagnetic wave in the context of geometric optics when absorption is involved. Either the TM wave condition (i.e., where the magnetic field of the refracted wave is transverse with respect to the wave direction) or the TE wave condition (i.e., where the electric field is transverse with respect to the propagating direction of the wave) may be assumed for the refracted wave in an absorbing medium to locally satisfy the electromagnetic boundary condition in the ray tracing calculation. The wave mode assumed for the refracted wave affects both the reflection and refraction coefficients. As a result, a nonunique solution for these coefficients is derived from the electromagnetic boundary condition. In this study we have identified the appropriate solution for the Fresnel reflection/refraction coefficients in light scattering calculation based on the ray tracing technique. We present the 3 x 2 refraction or transmission matrix that completely accounts for the inhomogeneity of the refracted wave in an absorbing medium. Using the Fresnel coefficients for an absorbing medium, we derive an asymptotic solution in an analytical format for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations. The asymptotic theory can produce reasonable accuracy in the phase function calculations in the infrared window region (wavelengths near 10 micron) if the particle size (in diameter) is on the order of 40 micron or larger. However, since strong absorption is assumed in the computation of the single-scattering albedo in the asymptotic theory, the single scattering albedo does not change with variation of the particle size. As a result, the asymptotic theory can lead to substantial errors in the computation of single-scattering albedo for small and moderate particle sizes. However, from comparison of the asymptotic results with the FDTD solution, it is expected that a convergence between the FDTD results and the asymptotic theory results can be reached when the particle size approaches 200 micron. We show that the phase function at side-scattering and backscattering angles is insensitive to particle shape if the random orientation condition is assumed. However, if preferred orientations are assumed for particles, the phase function has a strong dependence on scattering azimuthal angle. The single-scattering albedo also shows very strong dependence on the inclination angle of incident radiation with respect to the rotating axis for the preferred particle orientations.

Yang, Ping↗

Comparison of Micromechanics Theories for Modeling Discontinuous Long Fiber Reinforced Composites

Two micromechanics theories, the high fidelity generalized method of cells and Carrera Unified Formulation, are used to model a repeating unit cell of a discrete long fiber composite. The models are used to predict the effective properties of a single fiber with varying fiber-to-fiber end gap lengths, as well as the local stresses in the fiber and matrix. Both models produced similar results and demonstrated that they could capture the complex local load transfer due to the shear lag effect.

discrete long fiber composites↗

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The j and k dependencies of the line coupling and line mixing effects: Theoretical studies of the relaxation matrices of N 2 -broadened CH 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. Throughout this work, it is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Securely Aggregated Coded Matrix Inversion

Coded computing is a method for mitigating straggling workers in a centralized computing network, by using erasure-coding techniques. Federated learning is a decentralized model for training data distributed across client devices. In this work we propose approximating the inverse of an aggregated data matrix, where the data is generated by clients; similar to the federated learning paradigm, while also being resilient to stragglers. To do so, we propose a coded computing method based on gradient coding. We modify this method so that the coordinator does not access the local data at any point; while the clients access the aggregated matrix in order to complete their tasks. Here, the network we consider is not centrally administrated, and the communications which take place are secure against potential eavesdroppers.

97 MATHEMATICS AND COMPUTING↗

On the theory of quantum measurement

Many so called paradoxes of quantum mechanics are clarified when the measurement equipment is treated as a quantized system. Every measurement involves nonlinear processes. Self consistent formulations of nonlinear quantum optics are relatively simple. Hence optical measurements, such as the quantum nondemolition (QND) measurement of photon number, are particularly well suited for such a treatment. It shows that the so called 'collapse of the wave function' is not needed for the interpretation of the measurement process. Coherence of the density matrix of the signal is progressively reduced with increasing accuracy of the photon number determination. If the QND measurement is incorporated into the double slit experiment, the contrast ratio of the fringes is found to decrease with increasing information on the photon number in one of the two paths.

Haus, Hermann A.↗

Praseodymium in the formal +5 oxidation state

Praseodymium in the +5 oxidation state is a long-sought connection between lanthanide, early-transition and actinide metal redox chemistries. Unique among the lanthanide series, evidence for molecular pentavalent praseodymium species has been observed in the gas phase and noble gas matrix isolation conditions. Here we report the low-temperature synthesis and characterization of a molecular praseodymium complex in the formal +5 oxidation state, [Pr 5+ (NP t Bu 3 ) 4 ][X − ] (where t Bu = tert -butyl and X − = tetrakis(pentafluorophenyl)borate or hexafluorophosphate). Single-crystal X-ray diffraction, solution-state spectroscopic, solution magnetometric, density functional theory and multireference wavefunction-based methods indicate a highly multiconfigurational singlet ground state. Finally, an inverted ligand field drives this unique electronic structure, which establishes a critical link in understanding the bonding of high-valent metal complexes across the periodic table.

chemical bonding↗

Applications of Laplace transform methods to airfoil motion and stability calculations

This paper reviews the development of generalized unsteady aerodynamic theory and presents a derivation of the generalized Possio integral equation. Numerical calculations resolve questions concerning subsonic indicial lift functions and demonstrate the generation of Kutta waves at high values of reduced frequency, subsonic Mach number, or both. The use of rational function approximations of unsteady aerodynamic loads in aeroelastic stability calculations is reviewed, and a reformulation of the matrix Pade approximation technique is given. Numerical examples of flutter boundary calculations for a wing which is to be flight tested are given. Finally, a simplified aerodynamic model of transonic flow is used to study the stability of an airfoil exposed to supersonic and subsonic flow regions.

Edwards, J. W.↗

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗

Vibrational Entropy of Crystalline Solids from Covariance of Atomic Displacements

The vibrational entropy of a solid at finite temperature is investigated from the perspective of information theory. Ab initio molecular dynamics (AIMD) simulations generate ensembles of atomic configurations at finite temperature from which we obtain the N-body distribution of atomic displacements, ρN. We calculate the information-theoretic entropy from the expectation value of lnρN. At a first level of approximation, treating individual atomic displacements independently, our method may be applied using Debye–Waller B-factors, allowing diffraction experiments to obtain an upper bound on the thermodynamic entropy. At the next level of approximation we correct the overestimation through inclusion of displacement covariances. We apply this approach to elemental body-centered cubic sodium and face-centered cubic aluminum, showing good agreement with experimental values above the Debye temperatures of the metals. Below the Debye temperatures, we extract an effective vibrational density of states from eigenvalues of the covariance matrix, and then evaluate the entropy quantum mechanically, again yielding good agreement with experiment down to low temperatures. Our method readily generalizes to complex solids, as we demonstrate for a high entropy alloy. Further, our method applies in cases where the quasiharmonic approximation fails, as we demonstrate by calculating the HCP/BCC transition in Ti.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An SU(5) × U(1)' SUSY GUT with a “vector-like chiral” fourth family to fit all low energy data, including the muon g – 2

An additional generation of quarks and leptons and their SUSY counterparts, which are vector-like under the Standard Model gauge group but are chiral with respect to the new U(1) 3–4 gauge symmetry, are added to the Minimal Supersymmetric Standard Model (MSSM). We show that this model is a GUT and unifies the three SM gauge couplings and also the additional U(1) 3–4 coupling at a GUT scale of ≈ 5 × 10 16 GeV and explains the experimentally observed deviation of the muon g – 2. We also fit the quark flavor changing processes consistent with the latest experimental data and look at the effect of the new particles on the W boson mass without obviously conflicting with the observed masses of particles, CKM matrix elements, neutrino mixing angles, their mass differences, and the lepton-flavor violating bounds. This model predicts sparticle masses less than 25 TeV, with a gluino mass ≈ 2.3 – 3 TeV consistent with constraints, and one of the neutralinos as the LSP with a mass of ≈ 480 – 580 GeV, which is a potential dark matter candidate. The model is string theory motivated and predicts the VL quarks, leptons, a massive Z' and two Dirac neutrinos at the TeV scale and the branching ratios of μ → eγ, τ → μγ and τ → 3μ with BR(μ → eγ) within reach of future experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗