Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational principles”

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 847 records · Page 47

An adaptive grid/Navier-Stokes methodology for the calculation of nozzle afterbody base flows with a supersonic freestream

This paper describes an adaptive grid method for base flows in a supersonic freestream. The method is based on the direct finite-difference statement of the equidistribution principle. The weighting factor is a combination of the Mach number, density, and velocity first-derivative gradients in the radial direction. Two key ideas of the method are to smooth the weighting factor by using a type of implicit smoothing and to allow boundary points to move in the grid adaptation process. An AGARD nozzle afterbody base flow configuration is used to demonstrate the performance of the adaptive grid methodology. Computed base pressures are compared to experimental data. The adapted grid solutions offer a dramatic improvement in base pressure prediction compared to solutions computed on a nonadapted grid. A total-variation-diminishing (TVD) Navier-Stokes scheme is used to solve the governing flow equations.

Williams, Morgan↗

Compressing Hamiltonians with ab initio downfolding for simulating strongly-correlated materials on quantum computers

The accurate first-principles description of strongly correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving compressed many-body Hamiltonians that maintain the essential physics of strongly correlated materials. The solution of these material-specific models is still exponentially difficult to generate on classical computers, but quantum algorithms allow for a significant speed-up in obtaining the ground states of these compressed Hamiltonians. Here, we demonstrate that using quantum algorithms to obtain the properties of downfolded Hamiltonians can indeed yield high-fidelity solutions. By combining ab initio downfolding and variational quantum eigensolvers, we correctly predict the antiferromagnetic state of one-dimensional cuprate Ca 2 Cu O 3 , the excitonic ground state of monolayer W Te 2 , and the charge-ordered state of correlated metal Sr VO 3 . Numerical simulations using a classical tensor network implementation of variational quantum eigensolvers allow us to simulate large models with up to 54 qubits and encompassing up to four bands in the correlated subspace, which is indicative of the complexity that our framework can address. Through these methods we demonstrate the potential of classical preoptimization and downfolding techniques for enabling efficient materials simulation using quantum algorithms.

Alvertis, Antonios M. [NASA, Ames; LBNL, Berkeley]↗

Lunar Brightness Variations with Phase at 4.3-MM Wave Length

The lunar radio observations and interpretations of Piddington and Minnett (1949) and Gibson (1958) show that the lunar brightness variation with phase at millimeter wave lengths can be used to determine the physical properties of the lunar surface. They found that the millimeter-wave brightness lagged the optical phase, and their interpretation was that the millimeter radiation originates below a surface layer that is a very good thermal insulator. The thickness of this layer could not be determined from observations at one frequency. Observations at different frequencies give different results because of the wave-length dependence of the radio absorption by the surface material. The attenuation in the material increases with decreasing wave length, and therefore it is possible, in principle, to determine surface layer thickness from radio observations at several wave lengths. For this reason, observations of lunar radiation were started at the Naval Research Laboratory at a wave length of 4.3 mm. This is half the wave length used by Gibson (1958) in his earlier studies. The radio telescope used for these observations has been described in detail in a previous publication (Coates 1958). The antenna is a parabolic reflector 10 feet in diameter, and it has a beam width of 6.7 minutes of arc at the wave length of 4.3 mm. This is about one-fifth the diameter of the moon. The receiver was a Dicke-type radiometer.

Coates, Robert J.↗

Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory

The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density-matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies, relative to reference values from configuration-interaction-based calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radiation transfer through a model sunspot

A cylindrically symmetric model for a sunspot atmosphere using the similarity principle of Schlueter and Temesvary for the magnetic field configuration is presented. The equations of magnetostatic equilibrium are used, augmented by a radial Evershed flow. The LTE radiative transfer equations for the Stokes vector were solved under a variety of conditions for a ray emerging from a typical penumbral point. The contribution from isolated lines to the broadband circular polarization in sunspot penumbrae is evaluated using a more realistic model sunspot atmosphere than has hitherto been considered. Results indicate that the inclusion of a velocity field along the magnetic field vector is unable to give a net circular polarization of sufficient magnitude, although the variation with the angle between the line-of-sight and the magnetic field vector is in qualitative agreement with observations. The corresponding results for the net linear polarization are satisfactory.

Landman, D. A.↗

Anisotropic-strain-enhanced hole mobility in GaN by lattice matching to ZnGeN 2 and MgSiN 2

The key obstacle toward realizing integrated gallium nitride (GaN) electronics is its low hole mobility. Here, we explore the possibility of improving the hole mobility of GaN via epitaxial matching to II–IV nitride materials that have recently become available, namely, ZnGeN 2 and MgSiN 2 . We perform state-of-the-art calculations of the hole mobility of GaN using the ab initio Boltzmann transport equation. We show that effective uniaxial compressive strain of GaN along the [$1\bar{1}00$] by lattice matching to ZnGeN 2 and MgSiN 2 results in the inversion of the heavy hole band and split-off hole band, thereby lowering the effective hole mass in the compression direction. We find that lattice matching to ZnGeN 2 and MgSiN 2 induces an increase in the room-temperature hole mobility by 50% and 260% as compared to unstrained GaN, respectively. Further, examining the trends as a function of strain, we find that the variation in mobility is highly nonlinear; lattice matching to a hypothetical solid solution of Zn 0.75 Ge 0.75 Mg 0.25 Si 0.25 N 2 would already increase the hole mobility by 160%.

36 MATERIALS SCIENCE↗

Quantifying Large Lattice Relaxations in Photovoltaic Devices

Temporal variations of Cu(In,Ga)Se 2 photovoltaic device properties during light exposure at various temperatures and voltage biases for times up to 100 h are analyzed using the kinetic theory of large lattice relaxations. Open-circuit voltage and p-type doping increased with charge injection and decreased with temperature at low injection conditions. Furthermore, lattice relaxation can account for both trends and activation energies extracted from the data are approximately 0.9 and 1.2 eV for devices with lower and higher sodium content, respectively. In these devices, increased sodium content resulted in higher initial p-type doping with greater stability. First-principles calculations providing revised activation energies for the (V Se -V Cu ) complex suggest that this defect does not account for the metastability observed here.

14 SOLAR ENERGY↗

On the flight dynamics of aeroelastic vehicles

The nonlinear equations of motion for an elasic airplane are developed from first principles. Lagrange's Equation and the Principle of Virtual Work are used to generate the equations of motion and strip theory is then employed to obtain closed form integral expressions for the generalized forces. The inertial coupling is minimized by appropriate choice of the body reference axes and by making use of free vibration modes of the body. In addition, particular attention is paid to the simplifying assumptions used during the development of the equations of motion. A unique aspect of this modeling process is that since the generalized aerodynamic forces are determined from closed form, analytic expressions, this method can be used to gain insight into the effects of parameter variations not easily obtained from numerical models. A numerical example is also presented in which the modeling method is applied to a representative elastic aircraft. The model is used to address the effects of aerodynamic coupling which occurs between the rigid body degrees of freedom and the elastic degrees of freedom. Finally, model simplification is addressed and two methods are evaluated. The resulting frequency responses are compared.

Schmidt, D. K.↗

Age as a diagnostic of stratospheric transport

Estimates of stratospheric age from observations of long-lived trace gases with increasing tropospheric concentrations invoke the implicit assumption that an air parcel has been transported intact from the tropopical tropopause. However, because of rapid and irreversible mixing in the stratosphere, a particular air parcel cannot be identified with one that left the troposphere at some prior time. The parcel contains a mix of air with a range of transit times, and the mean value over this range is the most appropriate definition of age. The measured tracer concentration is also a mean over the parcel, but its value depends both on the transit time distribution and the past history of the tracer in the troposphere. In principle, only if the tropospheric concentration is increasing linearly can the age be directly inferred. We illustrate these points by employing both a one-dimensional diffusive analog of stratospheric transport, and the general circulation model (GCM) of the Goddard Institute for Space Studies (GISS). Within the limits of the GCM, we estimate the time over which tropospheric tracer concentrations must be approximately linear in order to determine stratospheric age unambiguously; the concentration of an exponentially increasing tracer is a function only of age if the growth time constant is greater than about 7 years, which is true for all the chlorofluorocarbons. More rapid source variations (for example, the annual cycle in CO2) have no such direct relationship with age.

Hall, Timothy M.↗

The imprint of proper motion of nonlinear structures on the cosmic microwave background

We investigate the imprint of nonlinear matter condensations on the cosmic microwave background (CMB) in an Omega = 1, cold dark matter (CDM) model universe. Temperature anisotropies are obtained by numerically evolving matter inhomogeneities and CMB photons from the beginning of decoupling until the present epoch. The underlying density field produced by the inhomogeneities is followed from the linear, through the weakly clustered, into the fully nonlinear regime. We concentrate on CMB temperature distortions arising from variations in the gravitational potentials of nonlinear structures. We find two sources of temperature fluctuations produced by time-varying potentials: (1) anisotropies due to intrinsic changes in the gravitational potentials of the inhomogeneities and (2) anisotropies generated by the peculiar, bulk motion of the structures across the microwave sky. Both effects generate CMB anisotropies in the range of 10(exp -7) approximately less than or equal to (Delta T/T) approximately less than or equal to 10(exp -6) on scales of approximately 1 deg. For isolated structures, anisotropies due to proper motion exhibit a dipole-like signature in the CMB sky that in principle could yield information on the transverse velocity of the structures.

Tuluie, Robin↗

Numerical studies of collinear laser-assisted injection from a foil for plasma wakefield accelerators

We present a laser-assisted electron injection scheme for beam-driven plasma wakefield acceleration. The laser is collinear with the driver and triggers the injection of hot electrons into the plasma wake by interaction with a thin solid target. We present a baseline case using the AWAKE Run 2 parameters and then perform variations on key parameters to explore the scheme. It is found that the trapped witness electron charge may be tuned by altering laser parameters, with a strong dependence on the phase of the wake upon injection. Normalized emittance settles at the order of micrometres and varies with witness charge. The scheme is robust to misalignment, with a 1/10th plasma skin-depth offset ( 20 μ m for the AWAKE case) having a negligible effect on the final beam. The final beam quality is better than similar existing schemes, and several avenues for further optimization are indicated. The constraints on the AWAKE experiment are very specific, but the general principles of this mechanism can be applied to future beam-driven plasma wakefield accelerator experiments. Published by the American Physical Society 2024

Physics↗

𝜂 and 𝜂′ Production in 𝐽/𝜓 Radiative Decays from Quantum Chromodynamics

We present a first principles calculation within quantum chromodynamics (QCD) of the radiative decays of the 𝐽/𝜓 into the light pseudoscalar mesons 𝜂 and 𝜂′. Within a lattice computation we obtain the transition form factors as a function of photon virtuality from the timelike region, accessible experimentally via the “Dalitz” decay 𝐽/𝜓 →𝑒 + ⁢𝑒 − ⁢𝜂 (′) , through to the real photon point corresponding to 𝐽/𝜓 →𝛾⁢𝜂 (′) . This is the first calculation in lattice QCD with two (heavier than physical) degenerate flavors of light quark and a heavier strange quark, in which the 𝜂′ appears as the first excited state with pseudoscalar isoscalar quantum numbers. We access it reliably by using variationally optimized operators , the use of which also improves the purity of the 𝐽/𝜓 and 𝜂 signals, reducing systematic uncertainties. High quality results at a large number of kinematic points are obtained in a typically noisy disconnected process by using a novel correlator averaging procedure. Our results show the expected enhanced production of the 𝜂′ over the 𝜂 in this process, and suggest that the demonstrated lattice technology is suitable for future calculations considering processes in which light meson resonances are produced.

form factors↗

Fundamental Aerodynamic Investigations for Development of Arrow-Stabilized Projectiles

The numerous patent applications on arrow-stabilized projectiles indicate that the idea of projectiles without spin is not new, but has appeared in various proposals throughout the last decades. As far as projectiles for subsonic speeds are concerned, suitable shapes have been developed for sometime, for example, numerous grenades. Most of the patent applications, though, are not practicable particularly for projectiles with supersonic speed. This is because the inventor usually does not have any knowledge of aerodynamic flow around the projectile nor any particular understanding of the practical solution. The lack of wind tunnels for the development of projectiles made it necessary to use firing tests for development. These are obviously extremely tedious or expensive and lead almost always to failures. The often expressed opinion that arrow-stabilized projectiles cannot fly supersonically can be traced to this condition. That this is not the case has been shown for the first time by Roechling on long projectiles with foldable fins. Since no aerodynamic investigations were made for the development of these projectiles, only tedious series of firing tests with systematic variation of the fins could lead to satisfactory results. These particular projectiles though have a disadvantage which lies in the nature cf foldable fins. They occasionally do not open uniformly in flight, thus causing unsymmetry in flow and greater scatter. The junctions of fins and body are very bad aerodynamically and increase the drag. It must be possible to develop high-performance arrow-stabilized projectiles based on the aerodynamic research conducted during the last few years at Peenemuende and new construction ideas. Thus the final shape, ready for operational use, could be developed in the wind tunnel without loss of expensive time in firing tests. The principle of arrow-stabilized performance has been applied to a large number of caliburs which were stabilized by various means Most promising was the development of a subcaliber wing-stabilized projectile with driving disc (Treibspiegel) where rigid control surfaces extend beyond the caliber of the projectile into the free stream. The stabilized projectiles of full-caliber, wing-stabilized projectiles with fins within the caliber is considerably more difficult. A completely satisfactory solution for the latter has not been found yet.

Kurzweg, Hermann↗

Modeling Wave Driven Non-linear Flow Oscillations: The Terrestrial QBO and a Solar Analog

The Quasi Biennial Oscillation (QBO) of the zonal circulation observed in the terrestrial atmosphere at low latitudes is driven by wave mean flow interaction as was demonstrated first by Lindzen and Holton (1968), shown in a laboratory experiment by Plumb and McEwan (1978), and modeled by others (e.g., Plumb, Dunkerton). Although influenced by the seasonal cycle of solar forcing, the QBO, in principle, represents a nonlinear flow oscillation that can be maintained by a steady source of upward propagating waves. The wave driven non-linearity is of third or odd order in the flow velocity, which regenerates the fundamental harmonic itself to keep the oscillation going - the fluid dynamical analog of the displacement mechanism in the mechanical clock. Applying Hines' Doppler Spread Parameterization (DSP) for gravity waves (GW), we discuss with a global-scale spectral model numerical experiments that elucidate some properties of the QBO and its possible effects on the climatology of the atmosphere. Depending on the period of the QBO, wave filtering can cause interaction with the seasonal variations to produce pronounced oscillations with beat periods around 10 years. Since the seasonal cycle and its variability influence the period of the QBO, it may also be a potent conduit of solar activity variations to lower altitudes. Analogous to the terrestrial QBO, we propose that a flow oscillation may account for the 22-year periodicity of the solar magnetic cycle, potentially answering Dicke (1978) who asked, "Is there a chronometer hidden deep inside the Sun?" The oscillation would occur below the convection region, where gravity waves can propagate. Employing a simplified, analytic model, Hines' DSP is applied to estimate the flow oscillation. Depending on the adopted horizontal wavelengths of GW's, wave amplitudes less than 10 m/s can be made to produce oscillating zonal flows of about 20 m/s that should be large enough to generate a significant oscillation in the magnetic field. For the large length scales of the Sun, the flow cycle period tends to be very long. The period, however, can be made to be 22 years, provided the buoyancy frequency (stability) is sufficiently small, thus placing the proposed flow near the base of the convection zone where a dynamo is now believed to operate.

Mayr, Hans G.↗

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗

Machine learning predictions for local electronic properties of disordered correlated electron systems

We present a scalable machine learning (ML) model to predict local electronic properties such as on-site electron number and double occupation for disordered correlated electron systems. Our approach is based on the locality principle, or the nearsightedness nature, of many-electron systems, which means local electronic properties depend mainly on the immediate environment. A ML model is developed to encode this complex dependence of local quantities on the neighborhood. We demonstrate our approach using the square-lattice Anderson-Hubbard model, which is a paradigmatic system for studying the interplay between Mott transition and Anderson localization. We develop a lattice descriptor based on the group-theoretical method to represent the on-site random potentials within a finite region. The resultant feature variables are used as input to a multilayer fully connected neural network, which is trained from data sets of variational Monte Carlo (VMC) simulations on small systems. We show that the ML predictions agree reasonably well with the VMC data. Our work underscores the promising potential of ML methods for multiscale modeling of correlated electron systems.

36 MATERIALS SCIENCE↗

D/H ratios and H2O contents of mantle-derived amphibole megacrysts from Dish Hill, California

D/H ratios are, in principle, useful in characterizing reservoirs of mantle hydrogen and as tracers of volatile transfer processes in Earth's interior. In practice, however, interpretation of isotopic measurements on mantle derived H is complicated by surface processes such as contamination and degassing which may alter the primary D/H ratio. Although there are indications that water associated with subduction zones and certain chemically enriched basalts is enriched in D relative to 'typical' upper mantle water, the extent of isotopic heterogeneity of mantle H remains uncertain. Kaersutitic amphibole megacrysts in alkaline basalts are one of the most widespread sources of mantle water and are therefore potentially useful for large-scale regional studies of D/H variation. However, D/H ratios of these amphiboles vary widely (from plus 8 to minus 113 percent), even in samples from the same locality, so that this potential has yet to be realized. In order to investigate the origin of this variability, and to explore the possibility that primary mantle D/H ratios may be deduced from these amphiboles, we analyzed the D/H ratios and chemical compositions of a suite of 17 kaersutitic amphiboles from Dish Hill, California. This work contrasts with previous studies in which sampling is widespread, but representatives from any given locality are few. Samples were collected from a restricted area on the southern flank of the volcanic center and are associated with the basal volcanic breccia. Fourteen of the samples were large single crystals or crystal fragments (megacrysts, 0.4 to 30 grams), believed to derive from pegmatitic veins crystallized from melts in the mantle. Two were coarse-grained intergrowths of amphibole with olivine and spinel, and one was a thin (2 mm) selvage on a peridotite xenolith.

Bell, David R.↗

Nonlinear Analysis of Bonded Composite Tubular Lap Joints

The present study describes a semi-analytical solution method for predicting the geometrically nonlinear response of a bonded composite tubular single-lap joint subjected to general loading conditions. The transverse shear and normal stresses in the adhesive as well as membrane stress resultants and bending moments in the adherends are determined using this method. The method utilizes the principle of virtual work in conjunction with nonlinear thin-shell theory to model the adherends and a cylindrical shear lag model to represent the kinematics of the thin adhesive layer between the adherends. The kinematic boundary conditions are imposed by employing the Lagrange multiplier method. In the solution procedure, the displacement components for the tubular joint are approximated in terms of non-periodic and periodic B-Spline functions in the longitudinal and circumferential directions, respectively. The approach presented herein represents a rapid-solution alternative to the finite element method. The solution method was validated by comparison against a previously considered tubular single-lap joint. The steep variation of both peeling and shearing stresses near the adhesive edges was successfully captured. The applicability of the present method was also demonstrated by considering tubular bonded lap-joints subjected to pure bending and torsion.

E. Oterkus↗