Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean”

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 73 records · Page 4

Mean-Field Ring Polymer Rates Using a Population Dividing Surface

Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cycle equivalence classes, orthogonal Weingarten calculus, and the mean field theory of memristive systems

Abstract It has been recently noted that for a class of dynamical systems with explicit conservation laws represented via projector operators, the dynamics can be understood in terms of lower dimensional equations. This is the case, for instance, of memristive circuits. Memristive systems are important classes of devices with wide-ranging applications in electronic circuits, artificial neural networks, and memory storage. We show that such mean-field theories can emerge from averages over the group of orthogonal matrices, interpreted as cycle-preserving transformations applied to the projector operator describing Kirchhoff’s laws. Our results provide insights into the fundamental principles underlying the behavior of resistive and memristive circuits and highlight the importance of conservation laws for their mean-field theories. In addition, we argue that our results shed light on the nature of the critical avalanches observed in quasi-two-dimensional nanowires as boundary phenomena.

97 MATHEMATICS AND COMPUTING↗

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optical conductivity of the two-dimensional Hubbard model: Vertex corrections, emergent Galilean invariance, and the accuracy of the single-site dynamical mean field approximation

We compute the frequency-dependent conductivity of the two-dimensional square lattice Hubbard model at zero temperature as a function of density to second order in the interaction strength, and compare the results to the predictions of single-site dynamical mean field theory computed at the same order. We find that despite the neglect of vertex corrections, the single-site dynamical mean field approximation produces semiquantitatively accurate results for most carrier concentrations, but fails qualitatively for the nearly empty or nearly filled band cases where the model exhibits an emergent Galilean invariance. The DMFT approximation also becomes qualitatively inaccurate very near half filling if nesting is important.

2-dimensional systems↗

Dynamics of fractionalized mean-field theories: Consequences for Kitaev materials

There have been substantial recent efforts, both experimentally and theoretically, to find a material realization of the Kitaev spin liquid—the ground state of the exactly solvable Kitaev model on the honeycomb lattice. Candidate materials are now plentiful, but the presence of non-Kitaev terms makes comparison between theory and experiment challenging. Here, we rederive time-dependent Majorana mean-field theory and extend it to include quantum phase information, allowing the direct computation of the experimentally relevant dynamical spin-spin correlator, which reproduces exact results for the unperturbed model. In contrast to previous work, we find that small perturbations do not substantially alter the exact result, implying that α-RuCl 3 is perhaps farther from the Kitaev phase than originally thought. Our approach generalizes to any correlator and to any model where Majorana mean-field theory is a valid starting point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multinucleon transfer mechanism in Cf 250 + Th 232 collisions using the quantal transport description based on the stochastic mean-field approach

Production cross sections of heavy neutron-rich isotopes are calculated by employing quantal transport description in 250 Cf + 232 Th collisions. This quantal transport description is based on the stochastic mean-field approach, and it provides a microscopic approach beyond time-dependent Hartree-Fock theory to include mean-field fluctuations. Deexcitation of primary fragments is determined by employing the statistical GEMINI++ code. Here, calculations provide predictions for production cross sections of neutron rich transfermium isotopes without any adjustable parameters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multinucleon transfer mechanism in Gd 160 + W 186 collisions in stochastic mean-field theory

The multinucleon transfer mechanism in 160 Gd + 186 W collisions is investigated in the framework of quantal transport description, based on the stochastic mean-field (SMF) theory. The SMF theory provides a microscopic approach for nuclear dynamics beyond the time-dependent Hartree-Fock approach by including mean-field fluctuations. Finally, cross sections for the primary fragment production are determined in the quantal transport description and compared with the available data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimal Electrification Using Renewable Energies: Microgrid Installation Model with Combined Mixture k-Means Clustering Algorithm, Mixed Integer Linear Programming, and Onsset Method

Optimal planning and design of microgrids are priorities in the electrification of off-grid areas. Indeed, in one of the Sustainable Development Goals (SDG 7), the UN recommends universal access to electricity for all at the lowest cost. Several optimization methods with different strategies have been proposed in the literature as ways to achieve this goal. This paper proposes a microgrid installation and planning model based on a combination of several techniques. The programming language Python 3.10 was used in conjunction with machine learning techniques such as unsupervised learning based on K-means clustering and deterministic optimization methods based on mixed linear programming. These methods were complemented by the open-source spatial method for optimal electrification planning: onsset. Four levels of study were carried out. The first level consisted of simulating the model obtained with a cluster, which is considered based on the elbow and k-means clustering method as a case study. The second level involved sizing the microgrid with a capacity of 40 kW and optimizing all the resources available on site. The example of the different resources in the Togo case was considered. At the third level, the work consisted of proposing an optimal connection model for the microgrid based on voltage stability constraints and considering, above all, the capacity limit of the source substation. Finally, the fourth level involved a planning study of electrification strategies based mainly on microgrids according to the study scenario. The results of the first level of study enabled us to obtain an optimal location for the centroid of the cluster under consideration, according to the different load positions of this cluster. Then, the results of the second level of study were used to highlight the optimal resources obtained and proposed by the optimization model formulated based on the various technology costs, such as investment, maintenance, and operating costs, which were based on the technical limits of the various technologies. In these results, solar systems account for 80% of the maximum load considered, compared to 7.5% for wind systems and 12.5% for battery systems. Next, an optimal microgrid connection model was proposed based on the constraints of a voltage stability limit estimated to be 10% of the maximum voltage drop. The results obtained for the third level of study enabled us to present selective results for load nodes in relation to the source station node. Finally, the last results made it possible to plan electrification using different network technologies and systems in the short and long term. The case study of Togo was taken into account. The various results obtained from the different techniques provide the necessary leads for a feasibility study for optimal electrification of off-grid areas using microgrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mean Field Theory of Self‐Organizing Memristive Connectomes

Abstract Biological neuronal networks are characterized by nonlinear interactions and complex connectivity. Given the growing impetus to build neuromorphic computers, understanding physical devices that exhibit structures and functionalities similar to biological neural networks is an important step toward this goal. Self‐organizing circuits of nanodevices are at the forefront of the research in neuromorphic computing, as their behavior mimics synaptic plasticity features of biological neuronal circuits. However, an effective theory to describe their behavior is lacking. This study provides for the first time an effective mean field theory for the emergent voltage‐induced polymorphism of circuits of a nanowire connectome, showing that the behavior of these circuits can be explained by a low‐dimensional dynamical equation. The equation can be derived from the microscopic dynamics of a single memristive junction in analytical form. The effective model is tested on experiments of nanowire networks and show that it fits both the potentiation and depression of these synapse‐mimicking circuits. It is shown that this theory applies beyond the case of nanowire networks by formulating a general mean‐field theory of conductance transitions in self‐organizing memristive connectomes.

Caravelli, Francesco↗

Number of sites-based solver for determining coverages from steady-state mean-field micro-kinetic models

Kinetic models parameterized by ab-initio calculations have led to significant improvements in understanding chemical reactions in heterogeneous catalysis. These studies have been facilitated by implementations which determine steady-state coverages and rates of mean-field micro-kinetic models. As implemented in the open-source kinetic modeling program, CatMAP, the conventional solution strategy is to use a root-finding algorithm to determine the coverage of all intermediates through the steady-state expressions, constraining all coverages to be non-negative and to properly sum to unity. Though intuitive, this root-finding strategy causes issues with convergence to solution due to these imposed constraints. In this work, we avoid explicitly imposing these constraints, solving the mean-field steady-state micro-kinetic model in the space of number of sites instead of solving it in the space of coverages. Here, we transform the constrained root-finding problem to an unconstrained least-squares minimization problem, leading to significantly improved convergence in solving micro-kinetic models and thus enabling the efficient study of more complex catalytic reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using observational mean-flow data to drive large-eddy simulations of a diurnal cycle at the SWiFT site

Reproducing realistic date- and site-specific unsteady wind conditions in large-eddy simulations is becoming increasingly useful in wind energy. How to run a large-eddy simulation to match observed conditions, however, remains an open research question. One approach that has received considerable attention is mesoscale-to-microscale coupling, in which information about the mesoscale weather, most commonly acquired from a mesoscale numerical weather model, is passed on to a microscale model. In this paper, we demonstrate how the recently developed profile-assimilation technique, a form of mesoscale-to-microscale coupling, can be used to drive large-eddy simulations solely based on observed mean-flow profiles at a single location, bypassing the need for auxiliary mesoscale simulations. The new approach is evaluated for a diurnal cycle at the Scaled Wind Farm Technology site. Observed mean-flow profiles from the ground up to a height of 2 km are reconstructed by aggregating measurements from multiple instruments, and gaps in the data are infilled with natural neighbor interpolation. We perform nine simulations using various forcing approaches to deal with data limitations. The results show that it is indeed possible to drive microscale large-eddy simulation with observations using the profile-assimilation technique, notwithstanding large gaps in virtual potential temperature measurements. However, profile assimilation with vertical smoothing of the error between the desired and actual profiles is required. Without that smoothing, the microscale simulations develop unrealistically high turbulence levels under many situations. Finally, we show that simulated mesoscale data can account for missing observations, although care is needed as both data sources are not necessarily compatible.

17 WIND ENERGY↗

Non-Boolean quantum amplitude amplification and quantum mean estimation

This paper generalizes the quantum amplitude amplification and amplitude estimation algorithms to work with non-Boolean oracles. The action of a non-Boolean oracle $U_\varphi $ on an eigenstate $\mathinner {|{x}\rangle }$ is to apply a state-dependent phase-shift $\varphi (x)$. Unlike Boolean oracles, the eigenvalues $\exp (i\varphi (x))$ of a non-Boolean oracle are not restricted to be $\pm 1$. Two new oracular algorithms based on such non-Boolean oracles are introduced. The first is the non-Boolean amplitude amplification algorithm, which preferentially amplifies the amplitudes of the eigenstates based on the value of $\varphi (x)$. Starting from a given initial superposition state $\mathinner {|{\psi _0}\rangle }$, the basis states with lower values of $\cos (\varphi )$ are amplified at the expense of the basis states with higher values of $\cos (\varphi )$. The second algorithm is the quantum mean estimation algorithm, which uses quantum phase estimation to estimate the expectation $\mathinner {\langle {\psi _0|U_\varphi |\psi _0}\rangle }$, i.e., the expected value of $\exp (i\varphi (x))$ for a random x sampled by making a measurement on $\mathinner {|{\psi _0}\rangle }$. It is shown that the quantum mean estimation algorithm offers a quadratic speedup over the corresponding classical algorithm. Both algorithms are demonstrated using simulations for a toy example. Potential applications of the algorithms are briefly discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigating the kinetic effect of prenol on iso -octane auto-ignition by means of an experimental and modeling study

A detailed experimental and kinetic modeling study was dedicated to understand the reported octane hyperboosting effect of prenol, by means of the measurement of the ignition delay times of its blends with iso -octane, and measurement of the mole fraction profiles of the fuels and intermediates inside the ULille rapid compression machine. Here, these results show that prenol addition leads to a reduction of the first-stage ignition phenomena and negative temperature coefficient behavior, which is only qualitatively captured by the model and is consistent with knock resistance improvement. It is suggested that this behavior is caused by two different factors. The first originates from gas-phase reactivity of prenol, and spans from the formation of unreactive unsaturated species through resonance-stabilized radicals, thereby constituting a competitive pathway for the radical pool generated by iso -octane. The second is of catalytic nature and cannot be captured by means of gas-phase kinetic modeling, but could also play an important role in the behavior of prenol in internal combustion engines.

42 ENGINEERING↗

Active populations and growth of soil microorganisms are framed by mean annual precipitation in three California annual grasslands

Climate influences soil microbial composition and function, but the relative importance of a site's historic climate versus its more immediate environmental conditions is unclear. Using quantitative stable isotope probing (qSIP), we characterized actively growing soil microbial communities and soil properties in three California annual grasslands that span a rainfall gradient and have developed on similar parent material. The soils were assayed in the wet winter season, when environmental conditions are most similar across sites. Since growing populations might be expected to be most responsive to contemporary environmental conditions, we hypothesized that the structure of growing microbial communities would be more similar across the gradient than that of total communities (i.e., including non-growing populations). In addition, we hypothesized that population growth rates would be slowest in the driest site, reflecting a legacy effect of low soil moisture on microbial growth. Soils along the rainfall gradient differed in pH, texture, and cation exchange capacity, but not in total C, C:N or dominant minerals. The radiocarbon (14C) age of soil C (reflecting turnover time) increased with mean annual precipitation but soil respiration was uniformly modern, reflecting microbial reliance on recent C inputs across the sites. The structure of both total and growing microbial communities differed across sites. Across major microbial phyla, including the Actinobacteria, Acidobacteria, Bacteroidetes, Gemmatimonadetes and Proteobacteria, bacterial growth rates were consistently lower in the site with the lowest mean annual precipitation. Taxa that were growing at the dry site alone grew more slowly than taxa that grew at multiple sites. These results reflect the influence of climate history and point to the role of environmental filtering at the driest site in shaping its slower growing microbial community, possibly reflecting adaptation to repeated exposure to water stress. Lastly, across taxa, the growth rate of a taxon at one site was correlated with its growth rate in the other sites. Furthermore, this growth rate coherence is likely a consequence of genetically determined physiological traits and is consistent with the idea that evolutionary history constrains growth rate.

Environmental filtering↗

Non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor instability using the macroscopic forcing method

The importance of non-locality of mean scalar transport in two-dimensional Rayleigh–Taylor Instability (RTI) is investigated. The macroscopic forcing method is utilized to measure spatio-temporal moments of the eddy diffusivity kernel representing passive scalar transport in the ensemble averaged fields. Presented in this work are several studies assessing the importance of the higher-order moments of the eddy diffusivity, which contain information about non-locality, in models for RTI. First, it is demonstrated through a comparison of leading-order models that a purely local eddy diffusivity is insufficient to capture the mean field evolution of the mass fraction in RTI. Therefore, higher-order moments of the eddy diffusivity operator are not negligible. Models are then constructed by utilizing the measured higher-order moments. It is demonstrated that an explicit operator based on the Kramers–Moyal expansion of the eddy diffusivity kernel is insufficient. An implicit operator construction that matches the measured moments is shown to offer improvements relative to the local model in a converging fashion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scalable Implementation of Mean-Field and Correlation Methods Based on Lie-Algebraic Similarity Transformation of Spin Hamiltonians in the Jordan–Wigner Representation

Recent work has highlighted that the strong correlation inherent in spin Hamiltonians can be effectively reduced by mapping spins to Fermions via the Jordan−Wigner transformation (JW). The Hartree−Fock method is straightforward in the Fermionic domain and may provide a reasonable approximation to the ground state. Correlation with respect to the Fermionic mean field can be recovered based on Lie-algebraic similarity transformation (LAST) with two-body correlators. Specifically, a unitary LAST variant eliminates the dependence on site ordering, while a nonunitary LAST yields size-extensive correlation energies. Whereas the first recent demonstration of such methods was restricted to small spin systems, we present efficient implementations using analytical gradients for the optimization with respect to the mean-field reference and the LAST parameters, thereby enabling the treatment of larger clusters, including systems with local spins s > $\frac{1}{2}$.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Annual Mean Arctic Amplification 1970–2020: Observed and Simulated by CMIP6 Climate Models

While the annual mean Arctic Amplification (AA) index varied between two and three during the 1970–2000 period, it reached values exceeding four during the first two decades of the 21st century. The AA did not change in a continuous fashion but rather in two sharp increases around 1986 and 1999. During those steps the mean global surface air temperature trend remained almost constant, while the Arctic trend increased. Although the “best” CMIP6 models reproduce the increasing trend of the AA in 1980s they do not capture the sharply increasing trend of the AA after 1999 including its rapid step-like increase. We propose that the first sharp AA increase around 1986 is due to external forcing, while the second step close to 1999 is due to internal climate variability, which models cannot reproduce in the observed time.

54 ENVIRONMENTAL SCIENCES↗

Emergent Constrained Projections of Mean and Extreme Warming in China

Reliable regional temperature projections including heat extremes are essential for climate change adaptation and mitigation. Taking China as an example, simple averages from Coupled Model Intercomparison Project Phase 6 (CMIP6) models project high warming due to sampling many high climate sensitivities in the ensemble. Here, we develop an emergent constraint (EC) framework to obtain constrained mean and daily maximum temperature (TXx) warming over China by using observed global warming and local residual warming. The constrained annual mean and TXx warming over China (2.33°C [1.61–3.05°C] and 2.31°C [1.21–2.99°C]) are 0.65°C [0.04–1.76°C] and 0.63°C [–0.50–2.39°C], respectively, lower than raw projections (2.98°C [1.85–4.22°C] and 2.94°C [2.04–4.39°C]) for 2080–2099 under the intermediate-emission scenario. Approximately half model uncertainty is reduced after constraint. The land area (population) experiencing temperature extremes in our metric is 78% (85%) of the raw projections. Our results imply a lower impact of extreme heat than implied by current raw CMIP6 projections.

54 ENVIRONMENTAL SCIENCES↗