Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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

Neutrino-Nucleus Interactions in the Precision Era

Current and future accelerator-based neutrino facilities, leveraging intense neutrino beams and advanced detectors, aim to precisely determine neutrino properties and probe signals of weakly interacting beyond the Standard Model physics. Achieving discovery-level precision and fully exploring the physics potential of these experiments critically depends on the accuracy of our understanding of fundamental underlying neutrino-nucleus interaction processes. This talk will focus on neutrino interactions spanning energies from tens of MeV to a few GeV—a complex, multi-scale and multi-process domain spanning from low-energy nuclear physics to perturbative QCD, with no unified underlying framework currently known. In this seminar, I will provide an overview of the field, discuss recent advancements, and share examples of ongoing cross-community efforts addressing these challenges.

Pandey, Vishvas [Fermilab]↗

Neutrino-Nucleus Interactions in the Precision Era

Current and future accelerator-based neutrino facilities, leveraging intense neutrino beams and advanced detectors, aim to precisely determine neutrino properties and probe signals of weakly interacting beyond the Standard Model physics. Achieving discovery-level precision and fully exploring the physics potential of these experiments critically depends on the accuracy of our understanding of fundamental underlying neutrino-nucleus interaction processes. This talk will focus on neutrino interactions spanning energies from tens of MeV to a few GeV a complex, multi-scale and multi-process domain spanning from low-energy nuclear physics to perturbative QCD, with no unified underlying framework currently known. In this talk, I will provide an overview of the field, discuss recent advancements, and share examples of ongoing cross-community efforts addressing these challenges.

Pandey, Vishvas [Fermilab] (ORCID:0000000230827987↗

Corrosion of Stainless and Carbon Steel in Aqueous Piperazine for CO2 Capture

Current obstacles that prevent commercial implementation of amine-scrubbing CO2 capture are the high costs. Reducing capital costs by appropriate selection of construction materials, which requires knowledge of material corrosion performance for the process, will improve the economic feasibility of this technology. Corrosion was evaluated in three pilot plant campaigns using aqueous piperazine with the Advanced Stripper (PZAS). 316L stainless steel (SS) experienced higher corrosion than 304 SS and 2205 duplex SS, and the corrosion rate showed strong dependence on the temperature. 304 and 2205 performed well at all locations and should be good construction materials for PZAS. Degraded PZ exacerbated 316L corrosion, and removal of PZ degradation products using a carbon adsorption bed significantly reduced corrosion. Carbon steel (CS) corrosion showed a weak temperature effect because the corrosion was more dependent on the protective siderite film. The protectiveness of the films was related to fluid velocity. Ni-based alloys corroded in PZ, and the rate increased with temperature. Corrosion of C1010 CS and SS (304, 316L, 430) was measured at absorber and water wash conditions on the bench-scale. Corrosion rate decreases with increasing PZ. 8 With more than 0.003 m PZ in solution, CS has acceptable corrosion performance. Corrosion of CS increases with increasing partial pressure of CO2, suggesting loading is another dominant parameter for carbon steel corrosion. Temperature has a less significant effect than PZ concentration and loading. CS corrosion increases with increasing flow velocity at both absorber and water wash conditions. SS had little corrosion at this lower temperature. Performance of siderite (FeCO3) protective films on CS was studied at representative stripper conditions on the bench-scale. Siderite films can deposit on the surface of CS in CO2-loaded PZ solution at temperatures >100 °C and protect CS from corrosion. The protection may fail in degraded PZ. Ethylenediamine (EDA) is one of the major contributors for the loss of film protectiveness or can be the surrogate for the effect of PZ degradation on siderite film protection. A link between protectiveness and the apparent density of siderite films was discovered. The apparent density of siderite films decreases with increasing flow velocity and decreasing CO2 loading, resulting in higher corrosion of CS.

Liu, Ching-Ting↗

A FETI approach to domain decomposition for meshfree discretizations of nonlocal problems

We propose a domain decomposition method for the efficient simulation of nonlocal problems. Our approach is based on a multi-domain formulation of a nonlocal diffusion problem where the subdomains share “nonlocal” interfaces of the size of the nonlocal horizon. This system of nonlocal equations is first rewritten in terms of minimization of a nonlocal energy, then discretized with a meshfree approximation and finally solved via a Lagrange multiplier approach in a way that resembles the finite element tearing and interconnect method. Specifically, we propose a distributed projected gradient algorithm for the solution of the Lagrange multiplier system, whose unknowns determine the nonlocal interface conditions between subdomains. Several two-dimensional numerical tests on problems as large as 191 million unknowns illustrate the strong and the weak scalability of our algorithm, which outperforms the standard approach to the distributed numerical solution of the problem. Finally, this work is the first rigorous numerical study in a two-dimensional multi-domain setting for nonlocal operators with finite horizon and, as such, it is a fundamental step towards increasing the use of nonlocal models in large scale simulations.

42 ENGINEERING↗

Testing dark matter halo properties using self-similarity

ABSTRACT We use self-similarity in N-body simulations of scale-free models to test for resolution dependence in the mass function and two-point correlation functions of dark matter haloes. We use 10243 particle simulations performed with abacus , and compare results obtained with two halo finders: friends-of-friends (fof ), and rockstar . The fof mass functions show a systematic deviation from self-similarity which is explained by resolution dependence of the fof mass assignment previously reported in the literature. Weak evidence for convergence is observed only starting from haloes of several thousand particles, and mass functions are overestimated by at least as much as $20-25{{\ \rm per\ cent}}$ for haloes of 50 particles. The mass function of the default rockstar halo catalogue (with bound virial spherical overdensity mass), on the other hand, shows good convergence of the order of 50 to 100 particles per halo, with no detectable evidence at the few percent level of any systematic dependence for larger particle number. Tests show that the mass unbinding procedure in rockstar is the key factor in obtaining this much improved resolution. Applying the same analysis to the halo–halo two point correlation function, we find again strong evidence for convergence only for rockstar haloes, at separations sufficiently large so that haloes do not overlap. At these separations, we can exclude dependence on resolution at the $5-10{{\ \rm per\ cent}}$ level once haloes have of the order of 50 to 100 particles. At smaller separations results are not converged even at significantly larger particle number, and bigger simulations would be required to establish the resolution required for convergence.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Reconnection in Strongly-Magnetized, Weakly-Collisional Plasmas: Onset, Turbulence, and Energy-Partition in 3D, Plasmoid-Dominated Regimes

This project is concerned with theoretical and numerical studies of magnetic reconnection. Key areas of focus are the reconnection onset (how reconnection gets started), energy partition (how the magnetic energy gets distributed during a reconnection event), and the effect of background turbulence on reconnection. We aim to shed light on these issues in kinetic plasmas (i.e., plasmas where the collisional mean free path is long compared to the length-scales of interest), with emphasis on 3D geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Beyond-mean-field approaches for nuclear neutrinoless double beta decay in the standard mechanism

Nuclear weak decays provide important probes to fundamental symmetries in nature. A precise description of these processes in atomic nuclei requires comprehensive knowledge on both the strong and weak interactions in the nuclear medium and on the dynamics of quantum many-body systems. In particular, an observation of the hypothetical double beta decay without emission of neutrinos (0νββ) would unambiguously demonstrate the Majorana nature of neutrinos and the existence of the lepton-number-violation process. It would also provide unique information on the ordering and absolute scale of neutrino masses. The next-generation tonne-scale experiments with sensitivity up to 10 28 years after a few years of running will probably provide a definite answer to these fundamental questions based on our current knowledge on the nuclear matrix element (NME), the precise determination of which is a challenge to nuclear theory. Beyond-mean-field approaches have been frequently adapted for the study of nuclear structure and decay throughout the nuclear chart for several decades. In this review, we summarize the status of beyond-mean-field calculations of the NMEs of 0νββ decay assuming the standard mechanism of an exchange of light Majorana neutrinos. In conclusion, the challenges and prospects in the extension and application of beyond-mean-field approaches for 0νββ decay are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Feedbacks Between Surface Deformation and Permafrost Degradation in Ice Wedge Polygons, Arctic Coastal Plain, Alaska

In the past three decades, an abrupt, pan-Arctic acceleration of ice wedge melting has transformed tundra landscapes, spurring the formation of hummock-like features known as high-centered polygons (HCPs). This rapid geomorphic transition profoundly alters regional hydrology and influences surface emissions of CO 2 and CH 4 . In Arctic Alaska, most recent instances of ice wedge degradation have arrested within 15–20 years of inception, stabilizing HCP microtopography. However, feedbacks between ground surface deformation and permafrost stability are incompletely understood, limiting our capacity to predict trajectories of landscape evolution in a still warmer future. In this work, we use field data from a site near Prudhoe Bay, Alaska, to develop a modeling-based framework for assessing the strength of positive (i.e., exacerbating) feedbacks on ice wedge degradation, focusing on the importance of heterogeneity in surface drainage and microtopographic conditions. Our simulations suggest that, when troughs are narrow, positive feedbacks on ice wedge melting (associated with thermokarst pool formation) are relatively weak. Positive feedbacks are markedly stronger beneath wide troughs, such as those that form above older, larger ice wedges. Seasonal thaw abruptly accelerates once a talik begins to form beneath wide and deep thermokarst pools. Once a talik initiates, winter severity and snowpack thickness increase in importance as predictors of thaw intensity in summer. Our results indicate that meter-scale heterogeneity in polygonal microtopography potentially exerts strong, nonlinear controls on thermokarst trajectories. These findings are useful for predicting future thermokarst dynamics and for interpreting the results from coarser-resolution land surface models operating at greater spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗

Thermal effects mediating the flow induced by laser-induced optical breakdown

Here, a short focused laser pulse can generate a hot plasma, which expands then contracts and can eject a hot jet, the speed and direction of which is sensitive to the details of the plasma kernel. The coupling of thermal and chemical nonequilibrium plasma mechanisms with this hydrodynamic development is assessed with detailed simulations of a two-temperature, three-species plasma model. Time scales for electron recombination, thermal relaxation, and diffusion are compared to that of the plasma expansion to anticipate conditions in which these mechanisms might affect the vorticity generation that leads to the ultimate flow pattern. The effect of these mechanisms are analyzed through comparison with corresponding inert-gas and equilibrium models. Thermal-nonequilibrium effects are found to be weak due to rapid relaxation of the heavy-particle and electronic temperatures. In contrast, chemical equilibration occurs at a rate comparable to the expansion and thereby enhance both it and subsequent hydrodynamic mechanisms as the energy stored in ion formation during the preceding breakdown is released by electron recombination. Thermal conduction, enhanced by high-temperature free electrons, weakens the ejection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quincke Oscillations of Colloids at Planar Electrodes

Dielectric particles in weakly conducting fluids rotate spontaneously when subject to strong electric fields. Such Quincke rotation near a plane electrode leads to particle translation that enables physical models of active matter. In this Letter, we show that Quincke rollers can also exhibit oscillatory dynamics, whereby particles move back and forth about a fixed location. We explain how oscillations arise for micron-scale particles commensurate with the thickness of a field-induced boundary layer in the nonpolar electrolyte. This work enables the design of colloidal oscillators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Supersymmetric proton decay revisited

Encouraged by the advent of a new generation of underground detectors – JUNO, DUNE and Hyper-Kamiokande – that are projected to improve significantly on the present sensitivities to various baryon decay modes, we revisit baryon decay in the minimal supersymmetric SU(5) GUT. We discuss the phenomenological uncertainties associated with hadronic matrix elements and the value of the strong coupling αs – which are the most important – the weak mixing angle θ W , quark masses including one-loop renormalization effects, quark mixing and novel GUT phases that are not visible in electroweak interaction processes. We apply our analysis to a variety of CMSSM, super- and sub-GUT scenarios in which soft supersymmetry-breaking parameters are assumed to be universal at, above and below the GUT scale, respectively. In many cases, we find that the next generation of underground detectors should be able to probe models with sparticle masses that are O(10) TeV, beyond the reach of the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Alternating Direction Decomposition with Strong Bounding and Convexification (ADDSBC) for Solving Security Constrained AC Unit Commitment Problems

This project aims to develop efficient and robust computational methods for solving the security-constrained unit commitment and alternating current optimal power flow problem (SC-UC-ACOPF). The SC-UC-ACOPF problem is at the center of the short-term operation of the U.S. Power Grid. It is solved every week, every day, and every 10 minutes to plan for the optimal action of electricity generation and consumption by minimizing the generation cost and maintaining power system reliability against potential disruptions of equipment failures. In mathematical terms, SC-UC-ACOPF is a challenging large-scale mixed-integer nonlinear optimization model. This means that the decisions involve both discrete variables, e.g. the turning on and off of generators and switching of transmission lines and transformers, and continuous decisions, e.g. the amount of energy generated by each generator and the power flows in the power grid. The physics of the power flow is described by nonlinear equations involving real and reactive power and bus voltages. Another key feature is the large number of contingencies, i.e. the system needs to stay reliable in face of failure of any one equipment, such as transmission lines and generators. The U.S. power grids are extremely complicated and large scale with more than 5,000 generators, 50,000 buses, and 100,000 high-voltage transmission lines, making the SC-UC-ACOPF a very large-scale computation challenge. The research developed in this project aims to solve the SC-UC-ACOPF problems in the three timescales, i.e. weekly, daily, and every 10-min. The proposed computational methods are built on a principled algorithmic approach of decomposition and penalization. More specifically, the algorithm develops spatial and temporal decomposition by exploiting the strong temporal coupling and weak spatial coupling of the UC problem and the complementary feature, i.e. weak temporal coupling and strong spatial coupling of the ACOPF problem. The algorithm also leverages recent progresses in strong convex relaxation of ACOPF. A unique feature of the proposed approach is that it generates a valid, global upper bound on the optimal maximum profit. In this way, a global optimality gap is available to measure the quality of the solution. To further speed up computation, the research team has developed a plethora of effective heuristics to strengthen the iterative penalty-based decomposition framework. For instance, a heuristic is developed to construct inner approximations of the time coupling constraints within the time decoupled problems. Contingencies are pre-screened and low-rank matrix computation is exploited to find the almost unique solution to each contingency. A novel heuristic for line switching is proposed and tested with positive impacts on instances where line switching is beneficial. Taking a systematic approach and carefully handling every detail of the problem pays off. The TIM-GO’s performance throughout the trials and the final event was stellar. TIM-GO garnered the second highest total prize money and is ranked in the top three positions across all categories of comparison.

97 MATHEMATICS AND COMPUTING↗

Weak Form Scientific Machine Learning: Test Function Construction for System Identification

Weak form Scientific Machine Learning (WSciML) is a recently developed framework for data-driven modeling and scientific discovery. It leverages the weak form of equation error residuals to provide enhanced noise robustness in system identification via convolving model equations with test functions, reformulating the problem to avoid direct differentiation of data. The performance, however, relies on wisely choosing a set of compactly supported test functions. In this work, we mathematically motivate a novel data-driven method for constructing Single-scale-Local reference functions for creating the set of test functions. Our approach numerically approximates the integration error introduced by the quadrature and identifies the support size for which the error is minimal, without requiring access to the model parameter values. Through numerical experiments across various models, noise levels, and temporal resolutions, we demonstrate that the selected supports consistently align with regions of minimal parameter estimation error. We also compare the proposed method against the strategy for constructing Multi-scale-Global (and orthogonal) test functions introduced in our prior work, demonstrating the improved computational efficiency.

FOS: Computer and information sciences↗

CCCP and MENeaCS: (updated) weak-lensing masses for 100 galaxy clusters

ABSTRACT Large area surveys continue to increase the samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at 0.05 < z < 0.15, observed as part of the Multi Epoch Nearby Cluster Survey, and re-evaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project. Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5 per cent, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of 1 − b = 0.84 ± 0.04 (stat) ±0.05 (syst). We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 0.07. We find a gas fraction of 0.139 ± 0.014 (stat) for eight relaxed clusters in our sample, which can also be used to infer cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

Monsoon Depression Amplification by Horizontal Shear and Humidity Gradients: A Shallow Water Perspective

Transient, synoptic-scale vortices produce a large fraction of total rainfall in most monsoon regions and are often associated with extreme precipitation. However, the mechanism of their amplification remains a topic of active research. For monsoon depressions, which are the most prominent synoptic-scale vortex in the Asian–Australian monsoon, recent work has suggested that meridional gradients in zonal wind in the vortex environment may produce growth through barotropic instability, while meridional gradients in environmental humidity have also been proposed to cause amplification through coupling with precipitating convection. Here, a two-dimensional shallow water model on a sphere with parameterized precipitation is used to examine the relative role played by these two environmental gradients. By systematically varying the meridional moisture gradient and meridional wind shear for both weak, quasi-linear waves and finite-amplitude isolated vortices, we show that rotational winds in the initial vortex are amplified most strongly by meridional shear of the environmental zonal wind, while vortex precipitation rates are most sensitive to environmental moisture gradients. The growth rate in the presence of both gradients is less than the sum of growth rates in the presence of isolated gradients, as the phase relation between moisture and vorticity anomalies becomes distorted with increasing shear. In conclusion, these results suggest that background meridional gradients in both zonal wind and environmental humidity can contribute to the amplification of vortices to monsoon depression strength, but with some degree of decoupling of the dry rotational flow and the moist convection.

54 ENVIRONMENTAL SCIENCES↗

Linear simulation of kinetic electromagnetic instabilities in a tokamak plasma with weak magnetic shear

Gyrokinetic simulation and eigenvalue calculation of electromagnetic instabilities are carried out for an experimentally observed low-n mode in weak magnetic shear discharge. With different magnetic shear values, the Gyrokinetic Toroidal Code simulation of the ion temperature gradient mode is consistent with the eigenvalue code (HD7) calculation. Due to the sensitivity of the kinetic ballooning mode (KBM) to global equilibrium, the simulation of the KBM deviates from the eigenvalue results, for the ballooning representation used in HD7 satisfies its spatial scale separations. Under a flat safety factor profile, the KBM is more unstable and its mode structure tends to move with the peak of the ion temperature drive. Further simulation of the KBM in an HL-2A-like equilibrium shows that the β excitation threshold of the mode is lower than 0.2% and the dominant toroidal mode number is n = 4, which is consistent with the measured experimental spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗

Smart Inverter Stability Enhancement in Weak Grids Using Adaptive Virtual-Inductance

An adaptive virtual-inductance feedforward scheme is presented in this article to enhance the stable operation of smart inverters in weak grids. In the developed scheme, the virtual inductance is updated in real time according to an adaptation law to ensure that the inverter follows the dynamic response of a robust reference model, thereby ensuring stability without the need for any grid impedance estimation. In this article, the developed scheme is validated with experimental results obtained from testing a small-scale two-level, 208-V, 2-kW inverter.

42 ENGINEERING↗