Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion 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 325 records · Page 18

Variational formulation of hybrid problems for fully 3-D transonic flow with shocks in rotor

Based on previous research, the unified variable domain variational theory of hybrid problems for rotor flow is extended to fully 3-D transonic rotor flow with shocks, unifying and generalizing the direct and inverse problems. Three variational principles (VP) families were established. All unknown boundaries and flow discontinuities (such as shocks, free trailing vortex sheets) are successfully handled via functional variations with variable domain, converting almost all boundary and interface conditions, including the Rankine Hugoniot shock relations, into natural ones. This theory provides a series of novel ways for blade design or modification and a rigorous theoretical basis for finite element applications and also constitutes an important part of the optimal design theory of rotor bladings. Numerical solutions to subsonic flow by finite elements with self-adapting nodes given in Refs., show good agreement with experimental results.

Liu, Gao-Lian↗

Maximum likelihood techniques applied to quasi-elastic light scattering

There is a necessity of having an automatic procedure for reliable estimation of the quality of the measurement of particle size from QELS (Quasi-Elastic Light Scattering). Getting the measurement itself, before any error estimates can be made, is a problem because it is obtained by a very indirect measurement of a signal derived from the motion of particles in the system and requires the solution of an inverse problem. The eigenvalue structure of the transform that generates the signal is such that an arbitrarily small amount of noise can obliterate parts of any practical inversion spectrum. This project uses the Maximum Likelihood Estimation (MLE) as a framework to generate a theory and a functioning set of software to oversee the measurement process and extract the particle size information, while at the same time providing error estimates for those measurements. The theory involved verifying a correct form of the covariance matrix for the noise on the measurement and then estimating particle size parameters using a modified histogram approach.

Edwards, Robert V.↗

Integrated Approach to the Dynamics and Control of Maneuvering Flexible Aircraft

This work uses a fundamental approach to the problem of simulating the flight of flexible aircraft. To this end, it integrates into a single formulation the pertinent disciplines, namely, analytical dynamics, structural dynamics, aerodynamics, and controls. It considers both the rigid body motions of the aircraft, three translations (forward motion, sideslip and plunge) and three rotations (roll, pitch and yaw), and the elastic deformations of every point of the aircraft, as well as the aerodynamic, propulsion, gravity and control forces. The equations of motion are expressed in a form ideally suited for computer processing. A perturbation approach yields a flight dynamics problem for the motions of a quasi-rigid aircraft and an 'extended aeroelasticity' problem for the elastic deformations and perturbations in the rigid body motions, with the solution of the first problem entering as an input into the second problem. The control forces for the flight dynamics problem are obtained by an 'inverse' process and the feedback controls for the extended aeroservoelasticity problem are determined by the LQG theory. A numerical example presents time simulations of rigid body perturbations and elastic deformations about 1) a steady level flight and 2) a level steady turn maneuver.

Waszak, Martin R.↗

Graph Theory Roots of Spatial Operators for Kinematics and Dynamics

Spatial operators have been used to analyze the dynamics of robotic multibody systems and to develop novel computational dynamics algorithms. Mass matrix factorization, inversion, diagonalization, and linearization are among several new insights obtained using such operators. While initially developed for serial rigid body manipulators, the spatial operators and the related mathematical analysis have been shown to extend very broadly including to tree and closed topology systems, to systems with flexible joints, links, etc. This work uses concepts from graph theory to explore the mathematical foundations of spatial operators. The goal is to study and characterize the properties of the spatial operators at an abstract level so that they can be applied to a broader range of dynamics problems. The rich mathematical properties of the kinematics and dynamics of robotic multibody systems has been an area of strong research interest for several decades. These properties are important to understand the inherent physical behavior of systems, for stability and control analysis, for the development of computational algorithms, and for model development of faithful models. Recurring patterns in spatial operators leads one to ask the more abstract question about the properties and characteristics of spatial operators that make them so broadly applicable. The idea is to step back from the specific application systems, and understand more deeply the generic requirements and properties of spatial operators, so that the insights and techniques are readily available across different kinematics and dynamics problems. In this work, techniques from graph theory were used to explore the abstract basis for the spatial operators. The close relationship between the mathematical properties of adjacency matrices for graphs and those of spatial operators and their kernels were established. The connections hold across very basic requirements on the system topology, the nature of the component bodies, the indexing schemes, etc. The relationship of the underlying structure is intimately connected with efficient, recursive computational algorithms. The results provide the foundational groundwork for a much broader look at the key problems in kinematics and dynamics. The properties of general graphs and trees of nodes and edge were examined, as well as the properties of adjacency matrices that are used to describe graph connectivity. The nilpotency property of such matrices for directed trees was reviewed, and the adjacency matrices were generalized to the notion of block weighted adjacency matrices that support block matrix elements. This leads us to the development of the notion of Spatial Kernel Operator SKO kernels. These kernels provide the basis for the development of SKO resolvent operators.

Jain, Abhinandan↗

Design of refractory multi-principal-element alloys for high-temperature applications

Abstract Refractory multi-principal-element alloys (RMPEAs) exhibit high specific strength at elevated temperatures ( T ). However, current RMPEAs lack a balance of room-temperature (RT) ductility, high- T strength, and high- T creep resistance. Using density-functional theory methods, we scanned composition space using four criteria: (1) formation energies for operational stability: $$-150\le {E}_{{\rm {f}}}$$ − 150 ≤ E f ≤ +70 meV per atom; (2) higher strength found via interstitial electron density with Young’s moduli E > 250 GPa; (3) inverse Pugh ratio for ductility: G / B < 0.57; and (4) high melting points: T m > 2500 °C. Using rapid bulk alloy synthesis and characterization, we validated theory and down-selected promising alloy compositions and discovered Mo 72.3 W 12.8 Ta 10.0 Ti 2.5 Zr 2.5 having well-balanced RT and high- T mechanical properties. This alloy has comparable high- T compressive strength to well-known MoNbTaW but is more ductile and more creep resistant. It is also superior to a commercial Mo-based refractory alloy and a nickel-based superalloy (Haynes-282) with improved high- T tensile strength and creep resistance.

36 MATERIALS SCIENCE↗

A finite-difference method for transonic airfoil design.

This paper describes an inverse method for designing transonic airfoil sections or for modifying existing profiles. Mixed finite-difference procedures are applied to the equations of transonic small disturbance theory to determine the airfoil shape corresponding to a given surface pressure distribution. The equations are solved for the velocity components in the physical domain and flows with embedded shock waves can be calculated. To facilitate airfoil design, the method allows alternating between inverse and direct calculations to obtain a profile shape that satisfies given geometric constraints. Examples are shown of the application of the technique to improve the performance of several lifting airfoil sections. The extension of the method to three dimensions for designing supercritical wings is also indicated.

Steger, J. L.↗

Absence of High-Pressure Ground-State Reentrant Ferroelectricity in PbTiO 3

We study ferroelectricity in the classic perovskite ferroelectric PbTiO 3 to high pressures with density functional theory (DFT) and experimental diamond-anvil techniques. We use second harmonic generation spectroscopy to detect lack of inversion symmetry. Consistent with early understanding and experiments, we find that ferroelectricity disappears at moderate pressures. However, DFT computations show that the disappearance arises from the overtaking of zone boundary instabilities, and not from the squeezing out of the off-centering ferroelectric displacements with pressure, as previously thought. Moreover, at high pressures the distorted perovskite phases are metastable with respect to a new dense centrosymmetric postperovskite phase with P⁢2 1 /m symmetry and 8-coordinated Ti, which becomes stable at about 70 GPa.

36 MATERIALS SCIENCE↗

First observation of increased DT yield over prediction due to addition of hydrogen

Here, in a number of reported instances, implosions utilizing fuel mixtures have resulted in anomalously low fusion yields below those predicted by radiation-hydrodynamics simulations. Inter-species ion diffusion has been suggested as a possible cause of the observed yield degradation in fuel mixture implosions. An experimental platform utilizing hydro-equivalent deuterium–tritium (DT), deuterium–tritium–hydrogen (DTH), and deuterium-tritium-helium3 (DT 3 He) capsule implosions was developed to determine whether the inter-species ion diffusion theory may describe the resulting fuel mixture implosion behavior. The implosion experiments were performed at the Omega laser facility. X-ray images and shell areal density diagnostics results show that the hydro-equivalent three capsules (DT, DTH, and DT 3 He) have similar compression behavior. However, nuclear yield deviation was observed from the scaling determined using a fusion yield formula. In the DT 3 He mixture, a reduced yield of a factor of 0.65 ± 0.13 was observed, which is similar to a yield reduction observed in D 3 He mixture by Rygg et al. (i.e., Rygg effect). In contrast, in the DTH mixture, a factor of 1.17 ±0.15 yield increase was observed, which we named the inverse Rygg effect. The yield increase observed in the DTH mixture is consistent with the inter-species ion diffusion theory where lighter H diffuses away from the core and concentrated DT in the core produces higher yield. An inter-species ion diffusion model, the Zimmerman–Paquette–Kagan–Zhdanov model, implemented in a Lagrangian radiation-hydrodynamics fluid code, was also used to analyze the present data, without the need to assume hydrodynamic equivalence of the capsules, but it does not completely explain the DTH or DT 3 He capsules although its effects are in the correct direction. Simulation-based Bayesian inference was used in the latter analysis to quantify the uncertainty in the numerical simulations. The simulation-based analysis resulted in an inferred Rygg-effect yield decrease factor of 0.91 ± 0.02 for the DT 3 He mixture, and an inferred inverse-Rygg yield increase factor of 1.21 ± 0.04 for the DTH mixture, based on simulations ignoring ion diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spin-orbit-proximitized ferromagnetic metal by monolayer transition metal dichalcogenide: Atlas of spectral functions, spin textures, and spin-orbit torques in Co/MoSe 2 , Co/WSe 2 and Co/TaSe 2 heterostructures

In this work, the heterostructures composed of an ultrathin ferromagnetic metal (FM) and a material hosting strong spin-orbit (SO) coupling are the principal resource for SO torque and spin-to-charge conversion nonequilibrium effects in spintronics. The key quantity in theoretical description of these effects is nonequilibrium spin density, which can appear on any monolayer of the heterostructure through which the current is flowing as long as the monolayer bands are affected by the native or proximity induced SO coupling. Here we demonstrate how hybridization of wave functions of Co layer and a monolayer of transition metal dichalcogenides (TMDs)—such as semiconducting MoSe 2 and WSe 2 or metallic TaSe 2 —can lead to dramatic transmutation of electronic and spin structure of Co within some distance away from its interface with TMD, when compared to the bulk of Co or its surface in contact with a vacuum. This is due to proximity induced SO splitting of Co bands encoded in the spectral functions and spin textures on its monolayers, which we obtain using noncollinear density functional theory (ncDFT) combined with equilibrium Green function (GF) calculations. In fact, SO splitting is present due to structural inversion asymmetry of the bilayer—i.e., just the presence of the interface—even if SO coupling within TMD monolayer is artificially switched off in ncDFT calculations. However, switching it on makes the effects associated with proximity SO coupling within Co layer about five times larger. Injecting spin-unpolarized charge current through SO-proximitized monolayers of Co generates nonequilibrium spin density over them, so that its cross product with the magnetization of Co determines SO torque. The SO torque computed via first-principles quantum transport methodology, which combines ncDFT with nonequilibrium GF calculations, can be used as the screening parameter to identify optimal combination of materials and their interfaces for applications in spintronics. In particular, we identify heterostructure two-monolayer-Co/monolayer-WSe 2 as the most optimal one, at least in the clean limit.

2-dimensional systems↗

Magnetic spectral signatures in the Earth's magnetosheath and plasma depletion layer

Correlations between plasma properties and magnetic fluctuations in the sub-solar magnetosheath downstream of a quasi-perpendicular shock have been found and indicate that mirror and ion cyclotronlike fluctuations correlate with the magnetosheath proper and plasma depletion layer, respectively (Anderson and Fueselier, 1993). We explore the entire range of magnetic spectral signatures observed from the Active Magnetospheric Particle Tracer Explorers/Charge Composition Explorer (AMPTE/CCE)spacecraft in the magnetosheath downstream of a quasi-perpendicular shock. The magnetic spectral signatures typically progress from predominantly compressional fluctuations,delta B(sub parallel)/delta B perpendicular to approximately 3, with F/F (sub p) less than 0.2 (F and F (sub p) are the wave frequency and proton gyrofrequency, respectively) to predominantly transverse fluctuations, delta B(sub parallel)/delta B perpendicular to approximately 0.3, extending up to F(sub p). The compressional fluctuations are characterized by anticorrelation between the field magnitude and electron density, n(sub e), and by a small compressibility, C(sub e) identically equal to (delta n(sub e)/n(sub e)) (exp 2) (B/delta B(sub parallel)) (exp 2) approximately 0.13, indicative of mirror waves. The spectral characteristics of the transverse fluctuations are in agreement with predictions of linear Vlasov theory for the H(+) and He(2+) cyclotron modes. The power spectra and local plasma parameters are found to vary in concert: mirror waves occur for beta(s ub parallel p) (beta (sub parallel p) identically = 2 mu(sub zero) n(sub p) kT (sub parallel p) / B(exp 2) approximately = 2, A(sub p) indentically = T(sub perpendicular to p)/T(sub parallel p) - 1 approximately = 0.4, whereas cyclotron waves occur for beta (sub parallel p) approximately = 0.2 and A(sub p) approximately = 2. The transition from mirror to cyclotron modes is predicted by linear theory. The spectral characteristics overlap for intermediate plasma parameters. The plasma observations are described by A(sub p) = 0.85 beta(sub parallel P) (exp - 0.48) with a log regression coefficient of -0.74. This inverse A(sub p) - beta(sub parallel p) correlation corresponds closely to the isocontours of maximum ion anisotropy instability growth, gamma (sub m)/omega(sub p) = 0.01, for the mirror and cyclotron modes. The agreement of observed properties and predictions of local theory suggests that the spectral signatures reflect the local plasma environment and that the anisotropy instabilities regulate A(sub p). We suggest that the spectral characteristics may provide a useful basis for ordering observations in the magnetosheath and that the A(sub p) - beta(sub parallel p) inverse correlation may be used as a beta-dependent upper limit on the proton anisotropy to represent kinetic effects.

Anderson, Brian J.↗

Level density within a micro-macroscopic approach

Statistical level density $ρ$($E, A$) is derived for nucleonic system with a given energy $E$, particle number $A$ and other integrals of motion in the micro-macroscopic approximation beyond the standard saddle-point method of the Fermi gas model. This level density reaches the two limits; the well-known Fermi gas grand-canonical ensemble limit for a large entropy $S$ related to large excitation energies, and the finite micro-canonical limit for a small combinatorical entropy $S$ at low excitation energies. In conclusion, the inverse level density parameter $K$ as function of the particle number $A$ in the semiclassical periodic orbit theory, taking into account the extended Thomas-Fermi and Strutinsky shell corrections, is calculated and compared with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Atomistic origin of the entropy of melting from inelastic neutron scattering and machine learned molecular dynamics

The latent heat, L, is central to melting, but its atomic origin remains elusive. It is proportional to the entropy of fusion, ΔS fus = L/T m (T m is the melting temperature), which depends on changes of atom configurations, atom vibrations, and thermal electron excitations. Here, we combine inelastic neutron scattering and machine-learned molecular dynamics to separate ΔS fus into these components for Ge, Si, Bi, Sn, Pb, and Li. When the vibrational entropy of melting, ΔS vib , is zero, ΔS fus ≃ 1.2 kB per atom. This result provides a baseline for ΔSconfig and nearly coincides with “Richard’s Rule” of melting. The ΔSfus deviates from this value for most elements, however, and we show that this deviation originates with extra ΔS vib and extra ΔS config . These two components are correlated for positive and negative deviations from Richard’s rule – the extra ΔS config is consistently ~ 80% of ΔS vib . Our results, interpreted with potential energy landscape theory, imply a correlation between the change in the number of basins and the change in the inverse of their curvature for the melting of pure elements.

36 MATERIALS SCIENCE↗

On the degeneracy of spin ice graphs, and its estimate via the Bethe permanent

The concept of spin ice can be extended to a general graph. We study the degeneracy of spin ice graph on arbitrary interaction structures via graph theory. We map spin ice graphs to the Ising model on a graph and clarify whether the inverse mapping is possible via a modified Krausz construction. From the gauge freedom of frustrated Ising systems, we derive exact, general results about frustration and degeneracy. We demonstrate for the first time that every spin ice graph, with the exception of the one-dimensional Ising model, is degenerate. We then study how degeneracy scales in size, using the mapping between Eulerian trails and spin ice manifolds, and a permanental identity for the number of Eulerian orientations. Furthermore, we show that the Bethe permanent technique provides both an estimate and a lower bound to the frustration of spin ices on arbitrary graphs of even degree. While such a technique can also be used to obtain an upper bound, we find that in all finite degree examples we studied, another upper bound based on Schrijver inequality is tighter.

97 MATHEMATICS AND COMPUTING↗

Coherent phonon pairs and rotational symmetry breaking of charge density wave order in the kagome superconductor CsV 3 Sb 5

Here, in this work we perform ultrafast time-resolved reflectivity measurements to study the symmetry breaking in the charge-density wave (CDW) phase of CsV 3 ⁢Sb 5 . By extracting the coherent phonon spectrum in the CDW phase of CsV 3 ⁢Sb 5 , we discover close phonon pairs near 1.3 THz and 3.1 THz, as well as a new mode at 1.84 THz. The 1.3 THz phonon pair and the 1.84 THz mode are observed up to the CDW transition temperature. Combining density-functional theory calculations, we show that these phonon pairs arise from the coexistence of Star of David and inverse Star of David distortions combined with sixfold rotational symmetry breaking. An anisotropy in the magnitude of transient reflectivity change is also revealed at the onset of CDW order. Our results thus indicate broken sixfold rotational symmetry in the charge-density wave state of CsV 3⁢ Sb 5 , along with the absence of nematic fluctuation above 𝑇 CDW . Meanwhile, the measured coherent phonon spectrum in the CDW phase of CsV 3 ⁢Sb 5−𝑥 ⁢Sn 𝑥 with 𝑥 = 0.03–0.04 matches the staggered inverse Star of David with interlayer 𝜋 phase shift. This CDW structure contrasts with undoped CsV3⁢Sb5 and explains the evolution from a phonon pair to a single mode at 1.3 THz by 𝑥 = 0.03–0.04 Sn doping.

Deng, Qinwen [Univ. of Pennsylvania, Philadelphia,↗

A two-dimensional model for the transport of pollutants in an urban basin

The distribution of an inert chemical species, carbon monoxide (CO), emitted mainly from automobile exhaust is investigated. The transport model, in analogy with the shallow water theory in fluid dynamics, considers variation of all physical quantities in the horizontal direction below the temperature inversion layer. Pollutants are found to be carried primarily by the wind; turbulent diffusion in a normal day plays only a minor role. The concentration of CO predicted by the present model for the entire Los Angeles basin is compared with observed data at nine stations. Accuracy of four different numerical schemes, the effect of turbulent diffusivity, and the source strengths are examined.

Liu, C. Y.↗

A two-dimensional model for the transport of pollutants in an urban basin

The distribution of carbon monoxide emitted mainly from an automobile exhaust is investigated. Carbon monoxide is assumed to be chemically inert. The transport model, in analogy with the shallow-water theory in fluid dynamics, considers variation of all physical quantities in the horizontal direction below the temperature inversion layer. Pollutants are found to be carried primarily by the wind; turbulent diffusion in a normal day plays only a minor role. The concentration of CO predicted by the present model for the entire Los Angeles basin is compared with observed data at nine stations. Accuracy of four different numerical schemes, the effect of turbulent diffusivity, and the source strengths are examined.

Liu, C. Y.↗

Ultrasonic characterization of changes in viscoelastic properties of epoxy during cure

A technique for using the longitudinal velocity (LV) of an ultrasonic wave to monitor the extent of cross linking (CL) during the cure of thermosetting resins is described. The method was developed by monitoring the rate of change in LV during the cure of a bisphenol-A epoxy resin with an amine adduct. The experiment included variations in the temperature and stoichiometry in order to express the rate of change in terms of the reaction kinetics. The pulse-echo method was used with a single transducer operating at 20 MHz. Numerical models were defined to account for the acoustic response of a single layer, the attenuation and the reflection coefficient. A linear relationship was observed between the inverse of the temperature and the log of the rate of change in the velocity, supporting the theory that the velocity could be used to monitor the extent of the cross-linking reaction. An activation energy of 11.9 kcal/mole was calculated for the mixture being investigated.

Winfree, W. P.↗

Unsteady thermocapillary migration of bubbles

Upon the introduction of a gas bubble into a liquid possessing a uniform thermal gradient, an unsteady thermo-capillary flow begins. Ultimately, the bubble attains a constant velocity. This theoretical analysis focuses upon the transient period for a bubble in a microgravity environment and is restricted to situations wherein the flow is sufficiently slow such that inertial terms in the Navier-Stokes equation and convective terms in the energy equation may be safely neglected (i.e., both Reynolds and Marangoni numbers are small). The resulting linear equations were solved analytically in the Laplace domain with the Prandtl number of the liquid as a parameter; inversion was accomplished numerically using a standard IMSL routine. In the asymptotic long-time limit, the theory agrees with the steady-state theory of Young, Goldstein, and Block. The theory predicts that more than 90 percent of the terminal steady velocity is achieved when the smallest dimensionless time, i.e., the one based upon the largest time scale-viscous or thermal-equals unity.

Dill, Loren H.↗