SEARCH · Search NASA

Results for “multipliers”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗

DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗

DOE OSTI2024

Structure preserving and energy dissipative contact approaches for implicit dynamics

In this work, several structure preserving and energy dissipative contact approaches are proposed and evaluated. The time integration schemes considered are general with regard to the version of constraint type, but here the emphasis was on mortar contact. The proposed mortar contact approach conserves both linear and angular momentum for mortar contact in a novel way. The proposed time integration scheme can conserve energy or provide strict contact dissipation. In addition, the proposed scheme enforces both gap constraints and gap velocity constraints (i.e., persistency). Using a midstep time integrator often causes energy dissipation during initial impact, here this energy can be recovered in a novel way. Enforcing the gap velocity constraint mitigates the contact chatter of the contact pressure and nodes in many problems. Whereas some approaches enforce the persistency condition and gap constraints simultaneously during the solution of the equations of motion (EOM) requiring multipliers for both constraints included in the equation set, here the gap constraint is solved through the equations of motion and the persistency condition is satisfied in the time integration scheme by the velocity update after the equations of motion. It is shown that this approach is strictly dissipative in that a plastic contact condition can be achieved. Analogous to a coefficient of restitution for rigid bodies, any dissipated energy can then be returned upon release if energy conservation is desired. Structure preserving methods are good for long-time dynamics simulations and energy conserving and strictly dissipative methods can overcome stability issues associated with standard time integration algorithms such as the Newmark method.

42 ENGINEERING↗

DOE OSTI2022

Topological monopole's gauge field‐induced anomalous Hall effect in artificial honeycomb lattice

Abstract Vortex magnetic structure in artificial honeycomb lattice provides a unique platform to explore emergent properties due to the additional Berry phase curvature imparted by chiral magnetization to circulating electrons via direct interaction. We argue that while the perpendicularly aligned magnetic component leads to the quantized flux of monopole at the center of the Berry sphere, the in‐plane vortex circulation of magnetization gives rise to unexpected nontrivial topological Berry phase due to the gauge field transformation. The unprecedented effect signifies the importance of vector potential in multiply connected geometrical systems. Experimental confirmations to proposed hypotheses are obtained from Hall resistance measurements on permalloy honeycomb lattice. Investigation of the topological gauge transformation due to the in‐plane chirality reveals anomalous quasi‐oscillatory behavior in Hall resistance as a function of perpendicular field. The oscillatory nature of is owed to the fluctuation in equilibrium current as a function of Fermi wave‐vector , envisaged under the proposed new formulation in this article. Our synergistic approach suggests that artificially tunable nanostructured material provides new vista to the exploration of topological phenomena of strong fundamental importance. Key Points Demonstration of gauge field flux due to vortex magnetism. Anomalous Hall effect in magnetic honeycomb lattice. Topological monopole‐induced magnetic flux.

Guo, Jiasen↗

DOE OSTI2024

e RPCA : Robust Principal Component Analysis for Exponential Family Distributions

Abstract Robust principal component analysis (RPCA) is a widely used method for recovering low‐rank structure from data matrices corrupted by significant and sparse outliers. These corruptions may arise from occlusions, malicious tampering, or other causes for anomalies, and the joint identification of such corruptions with low‐rank background is critical for process monitoring and diagnosis. However, existing RPCA methods and their extensions largely do not account for the underlying probabilistic distribution for the data matrices, which in many applications are known and can be highly non‐Gaussian. We thus propose a new method called RPCA for exponential family distributions (), which can perform the desired decomposition into low‐rank and sparse matrices when such a distribution falls within the exponential family. We present a novel alternating direction method of multiplier optimization algorithm for efficient decomposition, under either its natural or canonical parametrization. The effectiveness of is then demonstrated in two applications: the first for steel sheet defect detection and the second for crime activity monitoring in the Atlanta metropolitan area.

Zheng, Xiaojun↗

DOE OSTI2021

Avoiding Excess Computation in Asynchronous Evolutionary Algorithms

Asynchronous evolutionary algorithms are becoming increasingly popular as a means of making full use of many processors while solving computationally expensive search and optimization problems. These algorithms excel at keeping large clusters fully utilized, but may sometimes inefficiently sample an excess of fast-evaluating solutions at the expense of higher-quality, slow-evaluating ones. We introduce a steady-state parent selection strategy, SWEET (“Selection whilE EvaluaTing”), that sometimes selects individuals that are still being evaluated and allows them to reproduce early. This gives slow-evaluating individuals that have higher fitnesses an increased ability to multiply in the population. We find that SWEET appears effective in simulated take-over time analysis, but that its benefit is confined mostly to early in the run, and our preliminary study on an autonomous vehicle controller problem that involves tuning a spiking neural network proves inconclusive.

Scott, Eric↗

DOE OSTI2020

Search for lepton flavour violating decays of a neutral heavy Higgs boson to $\mu\tau$ and e$\tau$ in proton-proton collisions at $\sqrt{s}=$ 13 TeV

A search for lepton flavour violating decays of a neutral non-standard-model Higgs boson in the μτ and eτ decay modes is presented. The search is based on proton-proton collisions at a center of mass energy $ \sqrt{s} $ = 13 TeV collected with the CMS detector in 2016, corresponding to an integrated luminosity of 35.9 fb$^{−1}$. The τ leptons are reconstructed in the leptonic and hadronic decay modes. No signal is observed in the mass range 200–900 GeV. At 95% confidence level, the observed (expected) upper limits on the production cross section multiplied by the branching fraction vary from 51.9 (57.4) fb to 1.6 (2.1) fb for the μτ and from 94.1 (91.6) fb to 2.3 (2.3) fb for the eτ decay modes.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Universality (beyond leading log) of soft radiative corrections to $\hat q$ in $p_\perp$ broadening and energy loss

It has been known for many years that soft radiation can give potentially large double-logarithm corrections to $p_\perp$ broadening of a high-energy particle traveling through QCD matter, but that this soft radiation correction can be absorbed into an effective value $\hat q_{eff}$ for the medium $p_\perp$-broadening parameter $\hat q$. Here ``soft'' means high energy compared to medium scales but soft compared to the original high-energy particle traveling through the medium. A similar situation arises in the case of soft corrections to hard splitting of a high-energy particle, such as hard g → gg, where double logarithms can also be absorbed using the same effective $\hat q_{eff}$. In this paper, I study whether the same holds true for potentially large, subleading, {\it single}-logarthim corrections. The correspondence is more indirect for single logarithms, but I show (in the large-$N_c$ limit) that single logarithms from soft radiation in the case of $p_\perp$ broadening also determine the single logarithms from soft radiative corrections to hard g → gg splitting. Along the way, there is an interesting variation of the original BDMPS-Z calculation of splitting rates in the q(hat) approximation. I also discuss how, for soft-radiative corrections to hard splitting processes, there are two different types of ''$\hat q_{eff}$'' that come into play, which differ by ''$i\pi$'' terms that multiply single logarithms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Differential $t\overline{t}$ cross-section measurements using boosted top quarks in the all-hadronic final state with 139 fb -1 of ATLAS data

Measurements of single-, double-, and triple-differential cross-sections are presented for boosted top-quark pair-production in 13 TeV proton–proton collisions recorded by the ATLAS detector at the LHC. The top quarks are observed through their hadronic decay and reconstructed as large-radius jets with the leading jet having transverse momentum (pT) greater than 500 GeV. The observed data are unfolded to remove detector effects. The particle-level cross-section, multiplied by the $t\overline{t}$ $\rightarrow$ $WWb$$\overline{b}$ branching fraction and measured in a fiducial phase space defined by requiring the leading and second-leading jets to have p T > 500 GeV and p T > 350 GeV, respectively, is 331 ± 3(stat.) ± 39(syst.) fb. This is approximately 20% lower than the prediction of ${398}^{+48}_{-49}$ fb by POWHEG+PYTHIA 8 with next-to-leading-order (NLO) accuracy but consistent within the theoretical uncertainties. Results are also presented at the parton level, where the effects of top-quark decay, parton showering, and hadronization are removed such that they can be compared with fixed-order next-to-next-to-leading-order (NNLO) calculations. The parton-level cross-section, measured in a fiducial phase space similar to that at particle level, is 1.94 ± 0.02(stat.) ± 0.25(syst.) pb. This agrees with the NNLO prediction of ${1.96}^{+0.02}_{-0.17}$ pb. Reasonable agreement with the differential cross-sections is found for most NLO models, while the NNLO calculations are generally in better agreement with the data. The differential cross-sections are interpreted using a Standard Model effective field-theory formalism and limits are set on Wilson coefficients of several four-fermion operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Simultaneous CTEQ-TEA extraction of PDFs and SMEFT parameters from jet and $t\bar{t}$ data

Recasting phenomenological Lagrangians in terms of SM effective field theory (SMEFT) provides a valuable means of connecting potential BSM physics at momenta well above the electroweak scale to experimental signatures at lower energies. In this work we jointly fit the Wilson coefficients of SMEFT operators as well as the PDFs in an extension of the CT18 global analysis framework, obtaining self-consistent constraints to possible BSM physics effects. Global fits are boosted with machine-learning techniques in the form of neural networks to ensure efficient scans of the full PDF+SMEFT parameter space. We focus on several operators relevant for top-quark pair and jet production at hadron colliders and obtain constraints on the Wilson coefficients with Lagrange Multiplier scans. We find mild correlations between the extracted Wilson coefficients, PDFs, and other QCD parameters, and see indications that these correlations may become more prominent in future analyses based on data of higher precision. This work serves as a new platform for joint analyses of SM and BSM physics based on the CTEQ-TEA framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Measurements of the branching fractions of ${\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}$, ${\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+}$, and ${\Xi }_{c}^{+}\to {\Xi }^{0}{K}+$ at Belle and Belle II

Using 983.0 fb −1 and 427.9 fb −1 data samples collected with the Belle and Belle II detectors at the KEKB and SuperKEKB asymmetric energy e + e − colliders, respectively, we present studies of the Cabibbo-favored ${\Xi }_{c}^{+}$ decays ${\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}$ and ${\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+}$, and the singly Cabibbo-suppressed decay ${\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+}$. The ratios of branching fractions of ${\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}$ and ${\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+}$ relative to that of ${\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}$ are measured for the first time, while the ratio $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+})/\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+})$ is also determined and improved by an order of magnitude in precision. The measured branching fraction ratios are $\begin{array}{c}\frac{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0}\right)}{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}\right)}=0.067\pm 0.007\pm 0.003,\\ \frac{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+}\right)}{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}\right)}=0.251\pm 0.005\pm 0.010,\\ \frac{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+}\right)}{\mathcal{B}\left({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+}\right)}=0.017\pm 0.003\pm 0.001.\end{array}$ Additionally, the ratio $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+})/\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+})$ is measured to be 0.068 ± 0.010 ± 0.004. Here, the first and second uncertainties are statistical and systematic, respectively. Multiplying the ratios by the branching fraction of the normalization mode, $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+})=(2.9\pm 1.3)%$, we obtain the following absolute branching fractions $\begin{array}{c}\mathcal{B}({\Xi }_{c}^{+}\to {\Sigma }^{+}{K}_{S}^{0})=(0.194\pm 0.021\pm 0.009\pm 0.087)\text{%},\\ \mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{\pi }^{+})=(0.728\pm 0.014\pm 0.027\pm 0.326)\text{%},\\ \mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{0}{K}^{+})=(0.049\pm 0.007\pm 0.003\pm 0.022)\text{%},\end{array}$ where the third uncertainties are from $\mathcal{B}({\Xi }_{c}^{+}\to {\Xi }^{-}{\pi }^{+}{\pi }^{+})$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Drell-Yan angular lepton distributions at small x from TMD factorization.

The Drell-Yan process is studied in the framework of TMD factorization in the Sudakov region s » Q 2 » $ {q}_{\perp}^2 $ corresponding to recent LHC experiments with Q 2 of order of mass of Z-boson and transverse momentum of DY pair ~ few tens GeV. The DY hadronic tensors are expressed in terms of quark and quark-gluon TMDs with $ \frac{1}{Q^2} $ and $ \frac{1}{N_c^2} $ accuracy. It is demonstrated that in the leading order in N c the higher-twist quark-quark-gluon TMDs reduce to leading-twist TMDs due to QCD equation of motion. The resulting hadronic tensors depend on two leading-twist TMDs: f 1 responsible for total DY cross section, and Boer-Mulders function $ {h}_1^{\perp } $ . The corresponding qualitative and semi-quantitative predictions seem to agree with LHC data on five angular coefficients A 0 – A 4 of DY pair production. The remaining three coefficients A 5 – A 7 are determined by quark-quark-gluon TMDs multiplied by extra $ \frac{1}{N_c} $ so they appear to be relatively small in accordance with LHC results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Normalization of ZZ instanton amplitudes in type 0B minimal superstring theory

Abstract We study ZZ instanton corrections in the (2,4k)$$ \mathcal{N} $$ N = 1 minimal superstring theory with the type 0B GSO projection, which becomes the type 0B$$ \mathcal{N} $$ N = 1 super-JT gravity in thek→ ∞ limit. Each member of the (2,4k) family of theories has two phases distinguished by the sign of the Liouville bulk cosmological constant. The worldsheet method for computing the one-loop normalization constant multiplying the instanton corrections gives an ill-defined answer in both phases. We fix these divergences using insights from string field theory and find finite, unambiguous results. Each member of the (2,4k) family of theories is dual to a double-scaled one-matrix integral, where the double-scaling limit can be obtained starting either from a unitary matrix integral with a leading one-cut saddle point, or from a hermitian matrix integral with a leading two-cut saddle point. The matrix integral exhibits a gap-closing transition, which is the same as the double-scaled Gross-Witten-Wadia transition whenk= 1. We also compute instanton corrections in the double-scaled matrix integral for allkand in both phases, and find perfect agreement with the string theory results.

Physics↗

DOE OSTI2021

Search for the doubly charmed baryon $ {\varXi}_{cc}^{+} $ in the $ {\varXi}_c^{+}{\pi}^{-}{\pi}^{+} $ final state

A search for the doubly charmed baryon $ {\varXi}_{cc}^{+} $ is performed in the $ {\varXi}_c^{+}{\pi}^{-}{\pi}^{+} $ invariant-mass spectrum, where the $ {\varXi}_c^{+} $ baryon is reconstructed in the pK - π + final state. The study uses proton-proton collision data collected with the LHCb detector at a centre- of-mass energy of 13 TeV, corresponding to a total integrated luminosity of 5.4 fb -1 . No significant signal is observed in the invariant-mass range of 3.4–3.8 GeV/ c 2 . Upper limits are set on the ratio of branching fractions multiplied by the production cross-section with respect to the $ {\varXi}_{cc}^{++} $ →( $ {\varXi}_c^{+} $ → pK - π + ) π + decay for different $ {\varXi}_{cc}^{+} $ mass and lifetime hypotheses in the rapidity range from 2.0 to 4.5 and the transverse momentum range from 2.5 to 25 GeV/ c . The results from this search are combined with a previously published search for the $ {\varXi}_{cc}^{+} $ → $ {\varLambda}_c^{+} $ K - π + decay mode, yielding a maximum local significance of 4.0 standard deviations around the mass of 3620 MeV/ c 2 , including systematic uncertainties. Taking into account the look-elsewhere effect in the 3.5–3.7 GeV/ c 2 mass window, the combined global significance is 2.9 standard deviations including systematic uncertainties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Spectra for Toeplitz Operators Associated with a Constrained Subalgebra

A two-point algebra is a set of bounded analytic functions on the unit disk that agree at two distinct points $a, b \in \mathbb{D}$. This algebra serves as a multiplier algebra for the family of Hardy Hilbert spaces $H^2_t := \{ f \in H^2 : tf(b)\}$ where $t \in \mathbb{C} \cup \{\infty\}$. We show that various spectra of certain Toeplitz operators acting on these spaces are connected.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Broadband Fourier-Transform-Detected EPR at W-Band

This manuscript reports development of wideband Fourier transform (FT) electron paramagnetic resonance (EPR) detection capability at the uniquely high frequency of 94 GHz. Here, this technique is based on the HiPER spectrometer developed at the University of St. Andrews, into which we have integrated an arbitrary waveform generator (AWG) that is used to modulate the output from a solid-state multiplier chain prior to amplification, generating up to 1 kW microwave power with 1 GHz bandwidth. Benchmark experiments are presented for a standard TEMPOL radical, which comprises a 500 MHz broad EPR spectrum at 94 GHz. Using a single frequency-swept microwave pulse, efficient inversion of this spectrum is demonstrated, enabling frequency-dependent studies of the magnetization relaxation in the time-domain, again via broadband echo detection. From these measurements, an anisotropy in the spin–lattice relaxation time, T 1 , can be determined with ease for TEMPOL. In addition, we implement the FT detection scheme for multi-dimensional (electron–electron double resonance, or ELDOR) experiments, demonstrating the full capabilities of the HiPER spectrometer. As an example, we present FT-detected version of the ELDOR-detected nuclear magnetic resonance technique for the TEMPOL radical. We also describe potential applications of the new instrumentation for studies of gadolinium qubits.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

DOE OSTI2021

Timed to save: the added value of accounting for hourly incidence of electricity savings from residential space-conditioning measures

Abstract Previous research has recognized that the value of measures that reduce electricity usage depends upon the timing of the savings generated, but the lack of hourly savings shapes has limited the demonstration of this concept. We develop empirical hourly savings shapes for residential space-conditioning measures from nearly 18,000 efficiency projects in California and show how they combine with the diurnal and seasonal variation in electricity system costs. We find that these measures (cooling replacements; windows, doors, and skylights; and other envelope measures) tend to save electricity when system costs are highest and that the hourly savings account for 1.4–1.5 times as much value as non-time-sensitive estimates of efficiency would predict. We present these impact multipliers for each measure to quantify the additional value revealed by the time-sensitive approach. We show that this additional value is similar in an evolving electricity grid with storage, rather than natural gas generation, as the marginal resource.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗