Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Backpropagated frame transformation theory: A reformulation

In this work, the energy-dependent frame transformation theory of Gao and Greene [H. Gao and C. H. Greene, Phys. Rev. A 42, 6946 (1990)] is extended to yield a quantitatively accurate description of the dissociative recombination process. Evidence is presented to show that direct application of the original theory leads to inaccurate cross sections. A major revision, based on an interaction-free backpropagation of the Born-Oppenheimer solutions, markedly improves the frame transformation theory, reducing its average error by orders of magnitude. The original theory and its extension are tested on the previously explored two-dimensional (2D) model that is tailored to describe the singlet ungerade states of molecular hydrogen. The 2D model can be solved exactly (within the numerical accuracy) without implementing the Born-Oppenheimer approximation. These exact results then serve as a benchmark for the frame transformation theory developed in this paper.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A Cooperative Game Theory-based Approach to Under-frequency Load Shedding Control

This paper proposes a cooperative game theory-based under-frequency load shedding (UFLS) approach for frequency stability and control in power systems. UFLS is a crucial factor for frequency stability and control especially in power grids with high penetration of renewable energy sources and restructured power systems. Conventional UFLS methods, most of which are off-line, usually shed fixed amounts of predetermined loads based on a predetermined schedule which can lead to over or under curtailment of load. This paper presents a co-operative game theory-based two-stage strategy to effectively and precisely determine locations and amounts of loads to be shed for UFLS control. In the first stage, the total amount of loads to be shed, also referred to as deficit in generation or the disturbance power, is computed using the initial rate of change of frequency (ROCOF) referred to the equivalent inertial center. In the second stage, the Shapley value, one of the solution concepts of cooperative game theory, is used to determine load shedding amounts and locations. The proposed method is implemented on the reduced 9-bus 3-machine Western Electricity Coordinating Council (WECC) system and simulated on Real-time Digital Simulators (RTDS). The results show that the proposed UFLS approach can effectively return the system to normal state after disturbances.

42 ENGINEERING↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A Cooperative Game Theory-based Secondary Frequency Regulation in Distribution Systems

Participation of distribution systems in frequency regulation has become an important factor for grid operation after the integration of distributed energy resources (DERs). While available reserves from a single distribution system may not be sufficient for frequency regulation, aggregated reserves from several distribution systems can provide frequency regulation at grid-scale. This paper proposes a cooperative game theory-based approach for secondary frequency regulation in distribution systems consisting of distributed energy resources (DERs). A two-stage strategy is proposed to effectively and precisely determine a change in active power set points of DERs following a command from system operators. In the first stage, the value or worth of each DER and their coalitions are determined using initial rates of change of frequency (ROCOF) for each DER/coalition. In the second stage, the Shapley value, one of the solution concepts of cooperative game theory, is used to determine changes in active power set points of DERs. The proposed method is implemented on several distribution systems including a modified IEEE 13-node system and modified 33-node distribution system.

Gautam, Mukesh↗
DOE OSTI2021

Alluvial fan morphology: A self-similar free boundary problem description

In this work, we examine approximate geometrically self-similar solutions to a parabolic free boundary value problem applied to alluvial fan surface morphology and growth. Alluvial fans are fan- or cone-shaped sedimentary deposits caused by the rapid deposition of sediment from a canyon discharging onto a flatter plain. Longitudinal, topographic profiles of fans can be readily described by a seemingly time independent dimensionless profile (DeChant et al., 1999). However, because an alluvial fan can be expected to grow over time, it is clear that this “steady” profile is certainly time dependent and can be described using a space-time self-similar solution. In an experimental and theory-based study, Guerit et al. (2014) developed a self-similar (or as they describe it a self-affine) linear solution based upon an approximate first order small parameter expansion solution for a 1-d homogeneous nonlinear diffusion equation. Direct substitution of this result into a linear diffusion equation suggests that this first order expression may not fully satisfy the associated governing equation. In contrast, we develop a more complete solution based upon a modeled approximation for the axi-symmetric formulation such that the associated temporal behavior is consistent with a 1/3 time power-law as described by Reitz and Jerolmack (2014). The resulting expression is an exact solution to a linear heat equation. Furthermore, we emphasize that a small parameter is not inherent to the resulting profile result and is not included in our model development. Though developed using rather different approaches, the formal solution developed here is in good agreement with the simple polynomial described by DeChant et al. (1999) suggesting that this self-similar solution is a suitable time dependent representation of alluvial fan longitudinal profile form and improves on earlier work.

58 GEOSCIENCES↗
DOE OSTI2022

Eigenvalue-Based Micromagnetic Analysis of Switching in Spin-Torque-Driven Structures

We present an eigenvalue-based approach for studying the magnetization dynamics in magnetic nanostructures driven by spintronic excitations, such as spin-transfer torque and spin-orbit torque. The approach represents the system dynamics in terms of normal oscillation modes (eigenstates) with corresponding complex eigenfrequencies. The dynamics is driven by a small number of active eigenstates and often considering just a single eigenstate is sufficient. We develop a perturbation theory that provides semianalytical dynamic solutions by using eigenstates for the case in the absence of damping and spintronic excitations as a basis. The approach provides useful insights into dynamics in such systems and allows solving several difficulties in their modeling, such as extracting the switching current in magnetic random-access memories and understanding switching mechanisms. We show that the presented approach directly predicts the critical switching current, i.e., switching current for an infinite time. The approach also provides solutions for the switching dynamics allowing the switching current to be obtained for a finite switching time, provided that the system symmetry is broken, e.g., by tilting the polarizer, so that switching by a finite pulse is possible.

42 ENGINEERING↗
DOE OSTI2022

Xenon mobility in γ-uranium and uranium–molybdenum alloys

Diffusion in bcc uranium and U–Mo alloys is of great interest because fission gas and other fission products impact the performance of nuclear fuels. We investigate the mobility of xenon and molybdenum in bcc uranium (⁠$γ$-U) and metallic U–Mo alloys by calculating the migration energies of xenon and molybdenum for various local compositions using density functional theory. We also calculate the solute–vacancy binding energies of different solutes to vacancies in bcc uranium. We find that the solute–vacancy binding energy in bcc uranium is significantly higher than it is in other bcc metals (e.g., Fe and W). We also find that the migration energy of molybdenum is substantially higher than the migration energy of xenon, indicating that xenon is much more mobile than molybdenum in bcc uranium. The presence of molybdenum in the nearest-neighbor shell around a xenon atom typically increases the migration energy of xenon, which indicates a reduction of xenon mobility in U–Mo alloys compared to pure bcc uranium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Theory-Guided Inelastic Neutron Scattering of Crystalline Alkaline Aluminate Salts Bearing Principal Motifs of Solution-State Species

Aluminate salts precipitated from caustic alkaline solutions exhibit a correlation between the anionic speciation and the identity of the alkali cation in the precipitate, with the aluminate ions occurring either in monomeric (Al(OH) 4 – ) or dimeric (Al 2 O(OH) 6 2– ) forms. The origin of this correlation is poorly understood as are the roles that oligomeric aluminate species play in determining the solution structure, prenucleation clusters, and precipitation pathways. Characterization of aluminate solution speciation with vibrational spectroscopy results in spectra that are difficult to interpret because the ions access a diverse and dynamic configurational space. To investigate the Al(OH) 4 – and Al 2 O(OH) 6 2– anions within a well-defined crystal lattice, inelastic neutron scattering (INS) and Raman spectroscopic data were collected and simulated by density functional theory for K 2 [Al 2 O(OH) 6 ], Rb 2 [Al 2 O(OH) 6 ], and Cs[Al(OH) 4 ]·2H 2 O. These structures capture archetypal solution aluminate species: the first two salts contain dimeric Al 2 O(OH) 6 2– anions, while the third contains the monomeric Al(OH) 4 – anion. Here, comparisons were made to the INS and Raman spectra of sodium aluminate solutions frozen in a glassy state. In contrast to solution systems, the crystal lattice of the salts results in well-defined vibrations and associated resolved bands in the INS spectra. The use of a theory-guided analysis of the INS of this solid alkaline aluminate series revealed that differences were related to the nature of the hydrogen-bonding network and showed that INS is a sensitive probe of the degree of completeness and strength of the bond network in hydrogen-bonded materials. Results suggest that the ionic size may explain cation-specific differences in crystallization pathways in alkaline aluminate salts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Ion-dependent protein–surface interactions from intrinsic solvent response

Significance Hard–soft interfaces between inorganic surfaces and biomolecules promote self-assembly processes with broad implications in biogeochemistry, energy sciences, nanomedicine, and origins of life. Yet, detailed molecular-scale understanding of inorganic–biomolecule interactions and their dependence on solution conditions is missing. We present a theory for the initial stages of inorganic–biomolecule assembly based on the far-field response of water, using experimentally characterized interactions between muscovite surfaces and mica-binding proteins as model systems. Our work connects molecular details of the solution to assembly outcomes and suggests the initial driving forces for assembly are dominated by long-range, ion-specific interactions. The connections made between interfacial structure and long-range surface–biomolecule interactions provide insights toward a predictive understanding of biomolecular self-assembly on mineral surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Engineering Biomolecular Self–Assembly at Solid–Liquid Interfaces

Biomolecular self–assembly is a key process used by life to build functional materials from the “bottom up.” In the last few decades, bioengineering and bionanotechnology have borrowed this strategy to design and synthesize numerous biomolecular and hybrid materials with diverse architectures and properties. However, engineering biomolecular self–assembly at solid–liquid interfaces into predesigned architectures lags the progress made in bulk solution both in practice and theory. Here, recent achievements in programming self–assembly of peptides, proteins, and peptoids at solid–liquid interfaces are summarized and corresponding applications are described. Recent advances in the physical understandings of self–assembly pathways obtained using in situ atomic force microscopy are also discussed. Furthermore, these advances will lead to novel strategies for designing biomaterials organized at and interfaced with inorganic surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Through the interface: New insights of the hydrogen evolution and oxidation reactions in aqueous solutions

Abstract Hydrogen evolution and oxidation reactions (HER/HOR) are the most fundamental reactions in electrocatalysis. Despite the practical significance, the mechanisms of HER/HOR in aqueous solutions are still elusive. Various theories have been proposed to rationalize the pH effect, cation effect, and structure effect of HER/HOR but none of them can explain all observations. In this review, we discuss four schools of thought for the HER/HOR, focusing on the strengths and shortcomings of each hypothesis and highlighting the magnitude of electrochemical interface structure in hydrogen electrocatalysis.

Liu, Ershuai↗
DOE OSTI2021

Statistical learning for predicting density–matrix-based electron dynamics

In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

New heavenly double copies

The double copy relates scattering amplitudes and classical solutions in Yang-Mills theory, gravity, and related field theories. Previous work has shown that this has an explicit realisation in self-dual YM theory, where the equation of motion can be written in a form that maps directly to Plebański’s heavenly equation for self-dual gravity. The self-dual YM equation involves an area-preserving diffeomorphism algebra, two copies of which appear in the heavenly equation. In this paper, we show that this construction is a special case of a wider family of heavenly-type examples, by (i) performing Moyal deformations, and (ii) replacing the area-preserving diffeomorphisms with a less restricted algebra. As a result, we obtain a double-copy interpretation for hyper-Hermitian manifolds, extending the previously known hyper-Kähler case. We also introduce a double-Moyal deformation of the heavenly equation. The examples where the construction of Lax pairs is possible are manifestly consistent with Ward’s conjecture, and suggest that the classical integrability of the gravity-type theory may be guaranteed in general by the integrability of at least one of two gauge-theory-type single copies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Unsteady aerodynamic loads on pitching aerofoils represented by Gaussian body force distributions

The actuator line model (ALM) is an approach commonly used to represent lifting and dragging devices like wings and blades in large-eddy simulations (LES). The crux of the ALM is the projection of the actuator point forces onto the LES grid by means of a Gaussian regularisation kernel. The minimum width of the kernel is constrained by the grid size; however, for most practical applications like LES of wind turbines, this value is an order of magnitude larger than the optimal value that maximises accuracy. This discrepancy motivated the development of corrections for the actuator line, which, however, neglect the effect of unsteady spanwise shed vorticity. In this work we develop a model for the impact of spanwise shed vorticity on the unsteady loading of an aerofoil modelled as a Gaussian body force distribution, where the model is applicable within the regime of unsteady attached flow. The model solution is derived both in the time and frequency domain and features an explicit dependence on the Gaussian kernel width. We verify the model with ALM-LES for both pitch steps and periodic pitching. The model solution is compared with Theodorsen theory and validated with both computational fluid dynamics using body fitted grids and experiment. It is concluded that the optimal kernel width for unsteady aerodynamics is approximately 40 % of the chord. The ALM is able to predict the magnitude of the unsteady loading up to a reduced frequency of 𝑘 ≈ 0.2.

17 WIND ENERGY↗
DOE OSTI2020

Rietveld Refinement for Macromolecular Powder Diffraction

Powder X-ray diffraction (PXRD) has been employed extensively for the structural characterization of materials, quantitative analysis of multicomponent mixtures, phase identification, texture, and microstructural analysis. At the heart of all those lies the Rietveld method, which has revolutionized the use of powder diffraction for materials studies. Initially applied to minerals and inorganics, the Rietveld method has progressively been used for more complex materials, for instance, zeolites and pharmaceuticals. A major advance of the method came some 20 years ago, when the first protein structure was successfully refined. Because of the sheer complexity of macromolecules, several new approaches and algorithms had to be pioneered or adapted from macromolecular single-crystal diffraction experiments, thus constituting macromolecular PXRD a quite unique field of study. This review aims to provide necessary elements of theory and application of structure solution and refinement via the Rietveld method for macromolecular PXRD data. Here, practical explanations and highlighted case studies are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Spin polarization and correlation of quarks from the glasma

We investigate the interaction of strong color fields in the glasma stage of high-energy nuclear collisions with the spins of quarks and antiquarks. We employ the perturbative solution of the quantum kinetic theory for the spin transport of (massive) quarks in a background color field governed by the linearized Yang-Mills equation and derive expressions for the quark-spin polarization and quark-antiquark spin correlation at small momentum in terms of field correlators. For the Golec-Biernat–Wusthoff dipole distribution the quark-spin polarization vanishes, but the out-of-plane spin correlation of quarks and antiquarks is nonzero. Our order-of-magnitude estimate of the correlation far exceeds that caused by vorticity effects, but does not fully explain the data for vector meson alignment. We identify possible mechanisms that could further increase the predicted spin correlation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Thermal conductivity of the quasi-one-dimensional materials TaSe 3 and ZrTe 3

The high breakdown current densities and resilience to scaling of the metallic transition metal trichalcogenides TaSe 3 and ZrTe 3 make them of interest for possible interconnect applications, and it motivates this study of their thermal conductivities and phonon properties. These crystals consist of planes of strongly bonded one-dimensional chains more weakly bonded to neighboring chains. Phonon dispersions and the thermal conductivity tensors are calculated using density functional theory combined with an iterative solution of the phonon Boltzmann transport equation. The phonon velocities and the thermal conductivities of TaSe 3 are considerably more anisotropic than those of ZrTe 3 . The maximum LA velocity in ZrTe 3 occurs in the cross-chain direction, and this is consistent with the strong cross-chain bonding that gives rise to large Fermi velocities in that direction. The thermal conductivities are similar to those of other metallic two-dimensional transition metal dichalcogenides. At room temperature, a significant portion of the heat is carried by the optical modes. In the low frequency range, the phonon lifetimes and mean free paths in TaSe 3 are considerably shorter than those in ZrTe 3 . Here, we find that the shorter lifetimes in TaSe 3 are consistent with the presence of lower frequency optical branches and zone-folding features in the acoustic branches that arise due to the doubling of the TaSe 3 unit cell within the plane.

36 MATERIALS SCIENCE↗