Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Minimization”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ultrasonic deterrents provide no additional benefit over curtailment in reducing bat fatalities at an Ohio wind energy facility

Wind energy is important for achieving net-zero greenhouse gas emissions but also contributes to global bat mortality. Current strategies to minimize bat mortality due to collision with wind-turbine blades fall broadly into two categories: curtailment (limiting turbine operation during high-risk periods) and deterrence (discouraging bat activity near turbines). Recently, there has been interest in combining these strategies to achieve greater reductions in bat fatalities than either strategy might achieve in isolation. To investigate the effectiveness of combining curtailment with ultrasonic deterrent minimization strategies, we deployed six ultrasonic deterrents at nacelle height on 16 experimental turbines at Avangrid Renewables’ Blue Creek Wind Energy Facility. We rotated between four conditions (normal operations, curtailment only, deterrent only, curtailment and deterrent) randomly assigned to four wind turbines each night between 15 June and 3 October 2017. We found that bat mortality at wind turbines was independent of wind speed. The effectiveness of ultrasonic acoustic deterrents varied between high-frequency-calling species (eastern red bats) and low-frequency-calling species (hoary bats, silver-haired bats, and big brown bats). When deterrents were active, mortality was twice as high for eastern red bats compared to the control. Conversely, deterrents had a weak dampening effect on bat mortality for low-frequency species. We found no additive effects on mortality reduction for turbines operating both curtailment and deterrents compared to either approach in isolation. Our findings suggest that ultrasonic acoustic deterrents may not be effective for both high and low frequency echolocating bats. The increase in fatalities of eastern red bats is alarming and underscores the importance of considering site- and species-specific effects of minimization solutions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Chapter 6: Wind Energy: Effects on Bats

Utility-scale wind energy facilities require no fuel, consume no water, and produce no greenhouse gas emissions or other pollutants during energy production. Wind power currently supplies approximately 6.4% of the electricity consumed in the United States (U.S.), with continued growth expected in the coming years. Although further expansion of wind power is anticipated to provide environmental and economic benefits, there are increasing concerns about bat fatalities occurring at wind energy facilities across North America. Recent estimates place the number of bat fatalities in the several hundreds of thousands on an annual basis, and that number is projected to rise. As a result, population-level impacts of wind turbine-caused mortality are of increasing concern. For example, modeling efforts for a widespread species, the hoary bat (Lasiurus cinereus), indicate that the population could decline by 90% within 50 years, assuming no growth in installed wind power capacity and no significant implementation of conservation measures. Given increasing demand for wind energy and increasing evidence suggesting that bats are attracted to wind turbines, the need to develop cost-effective and practical impact minimization strategies is a high priority. Current strategies include siting restrictions and operational minimization, both of which limit wind power generation. Potential solutions that do not limit power generation include broadcasting ultrasound from wind turbines or using ultraviolet (UV) light to deter bats from approaching and entering the rotor-swept area where fatalities occur. Another possible solution is the development of turbine-surface materials that reduce the relative attractiveness of wind turbine towers to bats. In this chapter, we provide a succinct summary of known impacts of wind energy on bats and present current and future research priorities. We also describe the challenges and opportunities associated with developing and implementing effective solutions to minimize wind turbine-caused bat fatality.

17 WIND ENERGY↗
DOE OSTI2022

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗
DOE OSTI2020

Quantum complexity of time evolution with chaotic Hamiltonians

We study the quantum complexity of time evolution in large-N chaotic systems, with the SYK model as our main example. This complexity is expected to increase linearly for exponential time prior to saturating at its maximum value, and is related to the length of minimal geodesics on the manifold of unitary operators that act on Hilbert space. Using the Euler-Arnold formalism, we demonstrate that there is always a geodesic between the identity and the time evolution operator e -iHt whose length grows linearly with time. This geodesic is minimal until there is an obstruction to its minimality, after which it can fail to be a minimum either locally or globally. We identify a criterion — the Eigenstate Complexity Hypothesis (ECH) — which bounds the overlap between offdiagonal energy eigenstate projectors and the k-local operators of the theory, and use it to argue that the linear geodesic will at least be a local minimum for exponential time. We show numerically that the large-N SYK model (which is chaotic) satisfies ECH and thus has no local obstructions to linear growth of complexity for exponential time, as expected from holographic duality. In contrast, we also study the case with N = 2 fermions (which is integrable) and find short-time linear complexity growth followed by oscillations. Our analysis relates complexity to familiar properties of physical theories like their spectra and the structure of energy eigenstates and has implications for the hypothesized computational complexity class separations PSPACE * BQP/poly and PSPACE * BQSUBEXP/subexp, and the “fast-forwarding” of quantum Hamiltonians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Lattice-friendly gauge completion of a composite Higgs with top partners

We give an explicit example of a composite Higgs model with a pseudo-Nambu-Goldstone Higgs in which the top Yukawa coupling is generated via the partial compositeness mechanism. This mechanism requires composite top partners which are relatively light compared to the typical mass scale of the strongly coupled theory. While most studies of the phenomenology of such models have focused on a bottom-up approach with a minimal effective theory, a top-down approach suggests that that the theory should contain a limit in which an unbroken global chiral symmetry protects the mass of the top partners, and the spectrum of the partners satisfies `t Hooft matching conditions. We therefore consider a model for the UV gauge group which could provide a solution to the matching conditions, and note that the relatively light fermions and pseudo-Goldstone bosons fall into complete multiplets of a large approximate global symmetry. This implies that the spectrum of particles lighter than a few TeV is non-minimal. Our example illustrates likely features of a composite Higgs theory, and also serves as an example of a non-chiral theory with no sign problem and a possible solution to `t Hooft matching conditions. It would therefore be very interesting for a lattice exploration. We find in this example that for some low-energy parameters in the effective theory the top partners can decay into high multiplicity final states, which could be difficult for the Large Hadron Collider (LHC) to constrain. This may potentially allow for the top partners to be lighter than those in more minimal models.

Lattice Quantum Field Theory↗
DOE OSTI2021

A canonical purification for the entanglement wedge cross-section

In AdS/CFT we consider a class of bulk geometric quantities inside the entanglement wedge called reflected minimal surfaces. The areas of these surfaces are dual to the entanglement entropy associated to a canonical purification (the GNS state) that we dub the reflected entropy. From the bulk point of view, we show that half the area of the reflected minimal surface gives a reinterpretation of the notion of the entanglement wedge cross-section. We prove some general properties of the reflected entropy and introduce a novel replica trick in CFTs for studying it. The duality is established using a recently introduced approach to holographic modular flow. We also consider an explicit holographic construction of the canonical purification, introduced by Engelhardt and Wall; the reflected minimal surfaces are simply RT surfaces in this new spacetime. We contrast our results with the entanglement of purification conjecture, and finally comment on the continuum limit where we find a relation to the split property: the reflected entropy computes the von Neumann entropy of a canonical splitting type-I factor introduced by Doplicher and Longo.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Color-octet scalars in Dirac gaugino models with broken R symmetry

In this work we study the collider phenomenology of color-octet scalars ( sgluons ) in supersymmetric models with Dirac gaugino masses that feature an explicitly broken R symmetry ( R - broken models ). We construct such models by augmenting minimal R -symmetric models with a fairly general set of supersymmetric and softly supersymmetry-breaking operators that explicitly break R symmetry. We then compute the rates of all significant two-body decays and highlight new features that appear as a result of R symmetry breaking, including enhancements to extant decay rates, novel tree- and loop-level decays, and improved cross sections of single sgluon production. We demonstrate in some detail how the familiar results from minimal R -symmetric models can be obtained by restoring R symmetry. In parallel to this discussion, we explore constraints on these models from the Large Hadron Collider. We find that, in general, R symmetry breaking quantitatively affects existing limits on color-octet scalars, perhaps closing loopholes for light CP-odd (pseudoscalar) sgluons while opening one for a light CP-even (scalar) particle. Qualitatively, however, we find that — much as for minimal R -symmetric models, despite stark differences in phenomenology — scenarios with broken R symmetry and two sgluons below the TeV scale can be accommodated by existing searches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A realistic theory of E6 unification through novel intermediate symmetries

Abstract We propose a non-supersymmetric E 6 GUT with the scalar sector consisting of650⨁351′ ⨁27. Making use of the first representation for the initial symmetry breaking to an intermediate stage, and the latter two representations for second-stage breaking to the Standard Model and a realistic Yukawa sector, this theory represents the minimal E 6 GUT that proceeds through one of the intermediate stages that are novel compared to SU(5) or SO(10) GUT: trinification SU(3) C × SU(3) L × SU(3) R , SU(6) × SU(2) and flipped SO(10) × U(1). We analyze these possibilities under the choice of vacuum that preserves a ℤ 2 “spinorial parity”, which disentangles the chiral and vector-like fermions of E 6 and provides a dark matter candidate in the form of a (scalar) inert doublet. Three cases are shown to consistently unify under the extended survival hypothesis (with minimal fine-tuning): trinification symmetry SU(3) C × SU(3) L × SU(3) R with eitherLRorCRparity, and SU(6) CR × SU(2) L . Although the successful cases give a large range for proton lifetime estimates, all of them include regions consistent with current experimental bounds and within reach of forthcoming experiments. The scenario investigated in this paper essentially represents the unique (potentially) viable choice in the class of E 6 GUTs proceeding through a novel-symmetry intermediate stage, since non-minimal alternatives seem to be intrinsically non-perturbative.

Physics↗
DOE OSTI2023

Search for leptoquarks decaying into the bτ final state in pp collisions at $\sqrt{\textrm{s}}$ = 13 TeV with the ATLAS detector

A search for leptoquarks decaying into the $bτ$ final state is performed using Run 2 proton-proton collision data from the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb —1 at $\sqrt{\textrm{s}}$ = 13 TeV recorded by the ATLAS detector. The benchmark models considered in this search are vector leptoquarks with electric charge of 2/3e and scalar leptoquarks with an electric charge of 4/3e. No significant excess above the Standard Model prediction is observed, and 95% confidence level upper limits are set on the cross-section times branching fraction of leptoquarks decaying into bτ. For the vector leptoquark production two models are considered: the Yang-Mills and Minimal coupling models. In the Yang-Mills (Minimal coupling) scenario, vector leptoquarks with a mass below 1.58 (1.35) TeV are excluded for a gauge coupling of 1.0 and below 2.05 (1.99) TeV for a gauge coupling of 2.5. In the case of scalar leptoquarks, masses below 1.28 (1.53) TeV are excluded for a Yukawa coupling of 1.0 (2.5). Finally, an interpretation of the results with minimal model dependence is performed for each of the signal region categories, and limits on the visible cross-section for beyond the Standard Model processes are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Axion-gauge coupling quantization with a twist

The possible couplings of an axion to gauge fields depend on the global structure of the gauge group. If the Standard Model gauge group is minimal, or equivalently if fractionally charged color-singlet particles are forbidden, then the QCD axion’s Chern-Simons couplings to photons and gluons obey correlated quantization conditions. Specifically, the photon coupling can have a fractional part which is a multiple of 1/3, but which is determined by the gluon coupling. A consequence of this result is that, among all theories with a minimal gauge group and minimal axion coupling to gluons, the smallest possible axion-photon amplitude |g aγγ | arises for E/N = 8/3. This provides a new motivation for experiments targeting this axion-photon coupling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Reflected entropy in random tensor networks. Part III. Triway cuts

For general random tensor network states at large bond dimension, we prove that the integer Rényi reflected entropies (away from phase transitions) are determined by minimal triway cuts through the network. This generalizes the minimal cut description of bipartite entanglement for these states. A natural extrapolation away from integer Rényi parameters, suggested by the triway cut problem, implies the holographic conjecture S R = 2EW, where S R is the reflected entropy and EW is the entanglement wedge cross-section. Minimal triway cuts can be formulated as integer programs which cannot be relaxed to find a dual maximal flow/bit-thread description. This sheds light on the gap between the existence of tripartite entanglement in holographic states and the bipartite entanglement structure motivated by bit-threads. In particular, we prove that the Markov gap that measures tripartite entanglement is lower bounded by the integrality gap of the integer program that computes the triway cut.

AdS-CFT correspondence↗
DOE OSTI2022

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Closure of Internal Porosity in Continuous Casting Bloom During Heavy Reduction Process

To investigate the closure behavior of internal porosity (also referred to as internal void) in continuous casting bloom during heavy reduction (HR) and thus provide theoretical guidance for minimizing this kind of internal defect more effectively with HR, a three-dimensional (3D) mechanical model was developed based on the predicted temperature field by a 2D heat transfer model. With this 3D mechanical model, closure behaviors of internal porosity in continuous casting bloom during HR at and after the strand solidification end under different process conditions were numerically studied. It was found that the void axis length decreased significantly along the bloom thickness direction and increased slightly along the casting and bloom width directions after HR, and the influence of the initial void size on the void closure was not obvious. With a decrease of temperature difference between the bloom surface and center, HR efficiency for minimizing internal void decreased, while the required reduction force significantly increased. Compared with blooms with a uniform temperature distribution of 1100 °C, the void closure index after HR implemented at the strand solidification end was increased by ~ 25 pct. Compared with a conventional flat roll, the application of a convex roll during HR could contribute to minimizing the internal porosity more effectively and significantly enhance the reduction capacity of the withdrawal and straightening units. The void closure index of η{sub s} and η{sub v} (where η{sub s} and η{sub v} were defined based on the variation of the void aspect ratio and the void volume, respectively) was closely related to the equivalent strain (ε{sub eq}) and the hydrostatic integration parameter (Q), respectively, and two mathematical equations were derived to quantitatively describe the relationship of η{sub s} − ε{sub eq} and η{sub v} − Q.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Using Digital Waveforms to Mitigate Solvent Clustering During Mass Filter Analysis of Proteins

With advances in the precision of digital electronics, waveform generation technology has progressed to a state that enables the creation of m/z filters that are purely digitally driven. These advances present new methods of performing mass analyses that provide information from a chemical system that are inherently difficult to achieve by other means. One notable characteristic of digitally driven mass filters is the capacity to transmit ions at m/z ratios that vastly exceed the capabilities of traditional resonant systems. However, the capacity to probe ion m/z ratios that span multiple orders of magnitudes across multiple orders of magnitude presents a new set of issues requiring a solution. In the present work, when probing multiply charged protein species beyond m/z 2000 using a gentle atmospheric pressure interface, the presence of solvent adducts and poorly resolved multimers can severely degrade spectral fidelity. Increasing energy imparted into a target ion population is one approach minimizing these clusters; however, the use of digital waveform technology provides an alternative that maximizes ion transport efficiency and simultaneously minimizes solvent clustering. In addition to the frequency of the applied waveform, digital manipulation also provides control over the duty cycle of the target waveform. This work examines the conditions and approach leading to optimal digital waveform operation to minimize solvent clustering. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Investigation of rotating mode behavior in BWR out-of-phase limit cycle oscillations – Part 2: TRACE/PARCS model and physical explanation

Previous neutronic/thermal-hydraulic (TH) coupled numerical simulations using full-core TRACE/PARCS and SIMULATE-3K boiling water reactor (BWR) models have shown evidence of a specific “rotating mode” behavior (steady rotation of the symmetry line, i.e. constant phase shift of approximately 90° between the first two azimuthal modes) in BWR out-of-phase limit cycle oscillations, regardless of initial conditions and even if the first two azimuthal modes have different natural frequencies. This suggests a nonlinear coupling between these modes; otherwise, the phase shift between these modes would change at a constant rate during the limit cycle. The previous paper (“Part 1”) presented a series of results to examine this rotating behavior with a reduced-order model. The goal of the present study is to provide additional analyses of the predicted rotating mode behavior using higher-fidelity numerical modeling, as well as a physical explanation for why this mode is favored over side-to-side or other oscillatory behaviors from a TH perspective. Results are presented using TRACE and TRACE/PARCS for a small number of parallel channels, which confirmed that the conclusions developed from the reduced-order model remain applicable when applying a full two-fluid, six-equation, finite-volume modeling approach. From these results, a physical explanation has been put forth to explain why the rotating symmetry line behavior is preferred from a TH standpoint, demonstrating that predominantly out-of-phase unstable systems are most unstable when the variation in the total inlet flow rate is minimized (which minimizes the effective single-phase to two-phase pressure drop ratio) and that the rotating mode is the most successful in minimizing this total flow rate variation as compared with the side-to-side case or any other oscillation pattern. The conclusion is that the rotating mode will be favored for any out-of-phase unstable system of parallel channels with no neutronic feedback or relatively weak neutronic feedback. Here, previous analyses have indicated that systems with sufficiently strong neutronic coupling may favor the side-to-side oscillation mode over the rotating mode; this topic is left as a subject of future investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Geometry and speed of evolution for a spin-s system with long-range zz-type Ising interaction

We study the evolution of a spin-s system described by the long-range zz-type Ising interaction. The Fubini–Study metric of the quantum state manifold defined by this evolution is obtained. We explore the topology of this manifold and show that it corresponds to a sphere. Exploration of the Riemannian curvature allows us to determine the manifold geometry. Also we calculate the speed of evolution of the system and represent the curvature by means of this speed. This is important for an experimental measurement of the curvature. The conditions for achieving the minimal and maximal values of the speed of evolution are obtained. Also we examine the geometry of state manifold and speed of evolution of spin system in the thermodynamic limit. We propose the physical system of methane molecule for application of our considerations. Finally, we study the influence of an external magnetic field on the metric of state manifold and on the speed of evolution. In this case we obtain the conditions for achieving the minimal possible speed of evolution. For some predefined initial states the orientations of magnetic fields to reach the minimal and maximal values of the speed are found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Improved sampling technique to collect natural gas from hydrate-bearing pressure cores

High quality gas compositional data are an important factor in interpreting the genetic source of natural gas hosted in hydrate-bearing sediments and other subsurface systems. In order to accurately characterize the composition of gas samples degassed from hydrate-bearing pressure cores, one must use a reproducible sampling technique that minimizes artifacts of the sampling process. Herein, we review sediment core degassing techniques and compare data obtained from a commonly used degassing approach, which we term the standard quantitative degassing (SQD) technique, to our newly developed modified quantitative degassing (MQD) method designed to minimize atmospheric contamination and gas-water interactions. The SQD method allows sample gas to interact with water in a bubbling chamber, which we hypothesize could alter the gas composition following mixing with water or dissolved gases in the bubbling chamber. Whereas, the MQD method allows for the collection of sample gas prior to the bubbling chamber. To compare the SQD and MQD methods, we performed a side-by-side comparison of noble (He, Ne, Ar, Kr, and Xe), major (H 2 , N 2 , O 2 , and CO 2 ), and hydrocarbon (CH 4 , C 2 H 6 , C 3 H 8 , i-C 4 H 10 , C 4 H 10 , i-C 5 H 12 , C 5 H 12 ) gas concentrations and select isotopic compositions obtained using both sample collection techniques. Gas samples were collected from hydrate-bearing pressure cores recovered and maintained under hydrate stable conditions from the northern Gulf of Mexico during the UT-GOM 2 -1 Expedition. The MQD method displayed significantly lower concentrations of atmospheric gases, higher proportions of hydrocarbon gases, lower ratios of C 1 /C 2 + , and heavier stable carbon and hydrogen isotopes of methane than the SQD method. These results demonstrate that the MQD method reduced air contamination and minimized alteration of the hydrocarbon gases. Finally, we conclude this method may be important for future work that seeks to determine the composition of natural gas from pressure cores using quantitative degassing experiments, especially those seeking to measure major (e.g., N 2 ) and noble gases.

58 GEOSCIENCES↗
DOE OSTI2018

Generation and application of multivariate polynomial quadrature rules

The search for multivariate quadrature rules of minimal size with a specified polynomial accuracy has been the topic of many years of research. Finding such a rule allows accurate integration of moments, which play a central role in many aspects of scientific computing with complex models. The contribution of this paper is twofold. First, we provide novel mathematical analysis of the polynomial quadrature problem that provides a lower bound for the minimal possible number of nodes in a polynomial rule with specified accuracy. We give concrete but simplistic multivariate examples where a minimal quadrature rule can be designed that achieves this lower bound, along with situations that showcase when it is not possible to achieve this lower bound. Our second contribution is the formulation of an algorithm that is able to efficiently generate multivariate quadrature rules with positive weights on non-tensorial domains. Our tests show success of this procedure in up to 20 dimensions. We test our method on applications to dimension reduction and chemical kinetics problems, including comparisons against popular alternatives such as sparse grids, Monte Carlo and quasi Monte Carlo sequences, and Stroud rules. The quadrature rules computed in this paper outperform these alternatives in almost all scenarios.

97 MATHEMATICS AND COMPUTING↗