Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Feedback regulation of crystal growth by buffering monomer concentration

Abstract Crystallization is a ubiquitous means of self-assembly that can organize matter over length scales orders of magnitude larger than those of the monomer units. Yet crystallization is notoriously difficult to control because it is exquisitely sensitive to monomer concentration, which changes as monomers are depleted during growth. Living cells control crystallization using chemical reaction networks that offset depletion by synthesizing or activating monomers to regulate monomer concentration, stabilizing growth conditions even as depletion rates change, and thus reliably yielding desired products. Using DNA nanotubes as a model system, here we show that coupling a generic reversible bimolecular monomer buffering reaction to a crystallization process leads to reliable growth of large, uniformly sized crystals even when crystal growth rates change over time. Buffering could be applied broadly as a simple means to regulate and sustain batch crystallization and could facilitate the self-assembly of complex, hierarchical synthetic structures.

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

Mechanistic Pathways for N 2 O Elimination from trans -R 3 Sn-O-N=N-O-SnR 3 and for Reversible Binding of CO 2 to R 3 Sn-O-SnR 3 (R = Ph, Cy)

The rate and mechanism of the elimination of N 2 O from trans-R 3 Sn-O-N=N-O-SnR 3 (R = Ph ( 1 Ph ) and R = Cy ( 1 Cy )) to form R 3 Sn-O-SnR 3 (R = Ph ( 2 Ph ) and R = Cy ( 2 Cy )) have been studied using both NMR and IR techniques to monitor the reactions in the temperature range of 39–79 °C in C 6 D 6 . Activation parameters for this reaction are ΔH ‡ = 15.8 ± 2.0 kcal·mol –1 and ΔS ‡ = –28.5 ± 5 cal·mol –1 ·K –1 for 1 Ph and ΔH ‡ = 22.7 ± 2.5 kcal·mol –1 and ΔS ‡ = –12.4 ± 6 cal·mol –1 ·K –1 for 1 Cy . Addition of O 2 , CO 2 , N 2 O, or PPh 3 to sealed tube NMR experiments did not alter in a detectable way the rate or product distribution of the reactions. Computational DFT studies of elimination of hyponitrite from trans-Me 3 Sn-O-N=N-O-SnMe 3 ( 1 Me ) yield a mechanism involving initial migration of the R 3 Sn group from O to N passing through a marginally stable intermediate product and subsequent N 2 O elimination. Reactions of 1 Ph with protic acids HX are rapid and lead to formation of R 3 SnX and trans-H 2 N 2 O 2 . Reaction of 1 Ph with the metal radical •Cr(CO) 3 C 5 Me 5 at low concentrations results in rapid evolution of N 2 O. At higher •Cr(CO) 3 C 5 Me 5 concentrations, evolution of CO 2 rather than N 2 O is observed. Addition of 1 atm or less CO 2 to benzene or toluene solutions of 2 Ph and 2 Cy resulted in very rapid reaction to form the corresponding carbonates R 3 Sn-O-C(=O)-O-SnR 3 (R = Ph ( 3 Ph ) and R = Cy ( 3 Cy )) at room temperature. Evacuation results in fast loss of bound CO 2 and regeneration of 2 Ph and 2 Cy . Variable temperature data for formation of 3 Cy yield ΔH o = –8.7 ± 0.6 kcal·mol –1 , ΔS o = –17.1 ± 2.0 cal·mol –1 ·K –1 , and ΔG o 298K = –3.6 ± 1.2 kcal·mol –1 . Furthermore, DFT studies were performed and provide additional insight into the energetics and mechanisms for the reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

Solution of mixed convection heat transfer from isothermal in-line fins

Transient and steady state combined natural and forced convective flows over two in-line finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The first method of solution is based on the 'Simple Arbitrary Lagrangian Eulerian' (SALE) technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which the velocity field is updated by the effects of all forces, and (2) an Eulerian phase that executes all advective fluxes of mass, momentum and energy. The second method of solution uses the finite element code entitled FIDAP. In the first part, comparison of the results by FIDAP, SALE, and available experimental work were done and discussed for steady state forced convection over louvered fins. Good agreements were deduced between the three sets of results especially for the flow over a single fin. In the second part and in the absence of experimental literature, the numerical predictions were extended to the transient transports and to the opposing flow where pressure drop is reversed. Results are presented and discussed for heat transfer and pressure drop in assisting and opposing mixed convection flows.

Khalilollahi, Amir↗
DOE OSTI2011

Analytical and computational study of the ideal full two-fluid plasma model and asymptotic approximations for Hall-magnetohydrodynamics

The 5-moment two-fluid plasma model uses Euler equations to describe the ion and electron fluids and Maxwell's equations to describe the electric and magnetic fields. Two-fluid physics becomes significant when the characteristic spatial scales are on the order of the ion skin depth and characteristic time scales are on the order of the ion cyclotron period. The full two-fluid plasma model has disparate characteristic speeds ranging from the ion and electron speeds of sound to the speed of light. Two asymptotic approximations are applied to the full two-fluid plasma to arrive at the Hall-MHD model, namely negligible electron inertia and infinite speed of light. The full two-fluid plasma model and the Hall-MHD model are studied for applications to an electromagnetic plasma shock, geospace environmental modeling (GEM challenge) magnetic reconnection, an axisymmetric Z-pinch, and an axisymmetric field reversed configuration (FRC).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dimetalloylene (M‐E‐M) Complexes of Heavier Main Group Elements Ge, Sn, Pb, Bi via Cleavage of E‐X Bonds (X=N(SiMe 3 ) 2 , O t Bu) with an Iridium Hydride

Abstract Reactions of the Ir V hydride [ Me BDI Dipp ]IrH 4 {BDI=(Dipp)NC(Me)CH(Me)CN(Dipp); Dipp=2,6‐ i Pr 2 C 6 H 3 } with E[N(SiMe 3 ) 2 ] 2 (E=Sn, Pb) afforded the unusual dimeric dimetallotetrylenes ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐E) 2 in good yields. Moreover, ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐Ge) 2 was formed in situ from thermal decomposition of [ Me BDI Dipp ]Ir(H) 2 Ge[N(SiMe 3 ) 2 ] 2 . These reactions are accompanied by liberation of HN(SiMe 3 ) 2 and H 2 through the apparent cleavage of an E−N(SiMe 3 ) 2 bond by Ir−H. In a reversal of this process, ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐E) 2 reacted with excess H 2 to regenerate [ Me BDI Dipp ]IrH 4 . Varying the concentrations of reactants led to formation of the trimeric ([ Me BDI Dipp ]IrH 2 ) 3 ( μ 2 ‐E) 3 . The further scope of this synthetic route was investigated with group 15 amides, and ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐Bi) 2 was prepared by the reaction of [ Me BDI Dipp ]IrH 4 with Bi(NMe 2 ) 3 or Bi(O t Bu) 3 to afford the first example of a “naked” two‐coordinate Bi atom bound exclusively to transition metals. A viable mechanism that accounts for the formation of these products is proposed. Computational investigations of the Ir 2 E 2 (E=Sn, Pb) compounds characterized them as open‐shell singlets with confined nonbonding lone pairs at the E centers. In contrast, Ir 2 Bi 2 is characterized as having a closed‐shell singlet ground state.

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

Best of both worlds: Enforcing detailed balance in machine learning models of transition rates

The slow microstructural evolution of materials often plays a key role in determining material properties. When the unit steps of the evolution process are slow, direct simulation approaches such as molecular dynamics become prohibitive and Kinetic Monte-Carlo (kMC) algorithms, where the state-to-state evolution of the system is represented in terms of a continuous-time Markov chain, are instead frequently relied upon to efficiently predict long-time evolution. The accuracy of kMC simulations however relies on the complete and accurate knowledge of reaction pathways and corresponding kinetics. This requirement becomes extremely stringent in complex systems such as concentrated alloys where the astronomical number of local atomic configurations makes the a priori tabulation of all possible transitions impractical. Machine learning models of transition kinetics have been used to mitigate this problem by enabling the efficient on-the-fly prediction of kinetic parameters. While conventional KMC methods based on transition state theory naturally yield reversible dynamics that exactly obey the detailed balance criterion, providing strong guarantees on the properties of the stationary distribution, many recently-proposed ML-based approaches to barrier predictions provide no such guarantees. In this study, we derive conditions under which physics-informed ML architectures exactly enforce the detailed balance condition by construction, even when relying on non-extensive descriptions of states in terms of local environments around mobile defects. In conclusion, using the diffusion of a vacancy in a concentrated alloy as an example, we show that such ML architectures also exhibit superior performance in terms of prediction accuracy, demonstrating that the imposition of physical constraints can facilitate the accurate learning of barriers at no increase in computational cost.

36 MATERIALS SCIENCE↗
DOE OSTI2015

HybriSol Hybrid nanostructures for high-energy-density solar thermal fuels (Final Report)

In this project, we used theory and experiment to expand the concept of solar thermal fuels, in which a material absorbs sunlight and stores the energy in chemical form for later release as heat, like a solar cell plus a battery plus a resistor, but all in a single material that can undergo many reversible cycles. This technology can help solve the problem that solar energy is only available intermittently, namely when it is daytime and sunny, which is one of the problems holding back more widespread deployment of solar electricity generation. It also has many other promising applications in which solar energy is used to provide a clean and controllable source of heat on demand. De-icing of windshields and solar cookers for off-grid use are two particular applications we studied in this work. Our work advanced solar thermal fuels by theoretically identifying new candidate materials and approaches for higher performance, making new materials and validating predictions of their improved properties, and furthering the fundamental understanding of how these materials work via computer simulations. The prospects are good for solar thermal fuels since their basic functionality (photochromism) is a phenomenon that has long been known in chemistry, and our work has advanced the materials closer to the properties that will enable practical applications.

14 SOLAR ENERGY↗
DOE OSTI2012

Engineering of a novel Ca{sup 2+}-regulated kinesin molecular motor using a calmodulin dimer linker

Highlights: Black-Right-Pointing-Pointer Engineered kinesin-M13 and calmodulin involving single cysteine were prepared. Black-Right-Pointing-Pointer CaM mutant was cross-linked to dimer by bifunctional thiol reactive reagent. Black-Right-Pointing-Pointer Kinesin-M13 was dimerized via CaM dimer in the presence of calcium. Black-Right-Pointing-Pointer Function of the engineered kinesin was regulated by a Ca{sup 2+}-calmodulin dimer linker. -- Abstract: The kinesin-microtubule system holds great promise as a molecular shuttle device within biochips. However, one current barrier is that such shuttles do not have 'on-off' control of their movement. Here we report the development of a novel molecular motor powered by an accelerator and brake system, using a kinesin monomer and a calmodulin (CaM) dimer. The kinesin monomer, K355, was fused with a CaM target peptide (M13 peptide) at the C-terminal part of the neck region (K355-M13). We also prepared CaM dimers using CaM mutants (Q3C), (R86C), or (A147C) and crosslinkers that react with cysteine residues. Following induction of K355-M13 dimerization with CaM dimers, we measured K355-M13 motility and found that it can be reversibly regulated in a Ca{sup 2+}-dependent manner. We also found that velocities of K355-M13 varied depending on the type and crosslink position of the CaM dimer used; crosslink length also had a moderate effect on motility. These results suggest Ca{sup 2+}-dependent dimerization of K355-M13 could be used as a novel molecular shuttle, equipped with an accelerator and brake system, for biochip applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Signal amplification in a solid-state sensor through asymmetric many-body echo

Electronic spins of nitrogen–vacancy centres in diamond constitute a promising system for micro- and nanoscale magnetic sensing, because of their operation under ambient conditions, ease of placement in close proximity to sensing targets and biological compatibility. At high densities, the electronic spins interact through dipolar coupling, which typically limits but can also potentially enhance sensing performance. Here we report the experimental demonstration of many-body signal amplification in a solid-state, room-temperature quantum sensor. Our approach uses time-reversed two-axis-twisting interactions, engineered through dynamical control of the quantization axis and Floquet engineering in a two-dimensional ensemble of nitrogen–vacancy centres. We observe that optimal amplification occurs when the backward evolution time equals twice the forward evolution time, in sharp contrast to the conventional Loschmidt echo. These observations can be understood as resulting from an underlying time-reversed mirror symmetry of the microscopic dynamics, providing key insights into signal amplification and opportunities for practical entanglement-enhanced quantum sensing.

quantum information↗
DOE OSTI2016

Real-time microstructure imaging by Laue microdiffraction: A sample application in laser 3D printed Ni-based superalloys

Synchrotron-based Laue microdiffraction has been widely applied to characterize the local crystal structure, orientation, and defects of inhomogeneous polycrystalline solids by raster scanning them under a micro/nano focused polychromatic X-ray probe. In a typical experiment, a large number of Laue diffraction patterns are collected, requiring novel data reduction and analysis approaches, especially for researchers who do not have access to fast parallel computing capabilities. In this article, a novel approach is developed by plotting the distributions of the average recorded intensity and the average filtered intensity of the Laue patterns. Visualization of the characteristic microstructural features is realized in real time during data collection. As an example, this method is applied to image key features such as microcracks, carbides, heat affected zone, and dendrites in a laser assisted 3D printed Ni-based superalloy, at a speed much faster than data collection. Such analytical approach remains valid for a wide range of crystalline solids, and therefore extends the application range of the Laue microdiffraction technique to problems where real-time decision-making during experiment is crucial (for instance time-resolved non-reversible experiments).

36 MATERIALS SCIENCE↗
DOE OSTI2014

Influence of longitudinal isotope substitution on the thermal conductivity of carbon nanotubes: Results of nonequilibrium molecular dynamics and local density functional calculations

We report reverse nonequilibrium molecular dynamics calculations of the thermal conductivity of isotope substituted (10,10) carbon nanotubes (CNTs) at 300 K. {sup 12}C and {sup 14}C isotopes both at 50% content were arranged either randomly, in bands running parallel to the main axis of the CNTs or in bands perpendicular to this axis. It is found that the systems with randomly distributed isotopes yield significantly reduced thermal conductivity. In contrast, the systems where the isotopes are organized in patterns parallel to the CNTs axis feature no reduction in thermal conductivity when compared with the pure {sup 14}C system. Moreover, a reduction of approximately 30% is observed in the system with the bands of isotopes running perpendicular to the CNT axis. The computation of phonon dispersion curves in the local density approximation and classical densities of vibrational states reveal that the phonon structure of carbon nanotubes is conserved in the isotope substituted systems with the ordered patterns, yielding high thermal conductivities in spite of the mass heterogeneity. In order to complement our conclusions on the {sup 12}C-{sup 14}C mixtures, we computed the thermal conductivity of systems where the {sup 14}C isotope was turned into pseudo-atoms of 20 and 40 atomic mass units.

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

Reaction Temperature Manipulation as a Process Intensification Approach for CO 2 Absorption

Reactor temperature manipulation to increase product yields of chemical reactions is a known technique used in many industrial processes. In the case of exothermic chemical reactions, the well-known Le Chatelier’s principle predicts that a decrease in temperature will displace the chemical reaction toward the formation of products by increasing the value of the equilibrium constant. The reverse is true for endothermic reactions. Reactor temperature manipulation in an industrial system, however, affects the values of many variables, including physical properties, transport parameters, reaction kinetic parameters, etc. In the case of reactive absorption, some variables change with increasing temperatures due to solute absorption, while others change in such a way that the solute absorption rate decreases. For example, temperature drop increases product formation for exothermic reactions but reduces the value of transport parameters, leading to decreasing interfacial concentrations and absorption rates. Therefore, temperature manipulation strategies must be designed carefully to achieve the process goals. In this work, we theoretically study the use of temperature as a tool to increase CO 2 absorption by solvents in a semi-batch reactor. A computer code has been developed and validated using reported experimental data. Calculated results demonstrate an increase in absorbed CO 2 of more than 28% with respect to the highest temperature used. Despite high agitation and high gas flow rate, the system is mass transfer controlled at short times, becoming kinetically controlled as time increases. An operating strategy to decrease cooling energy costs is also proposed. This study reveals that reactor temperature manipulation can be an effective process to improve CO 2 absorption by solvents in two-phase semi-batch reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Configuration and Intraseasonal Duration of Interannual Anomalies of the Great Plains Low-Level Jet

Despite the fact that the low-level jet of the southern Great Plains (the GPLLJ) of the U.S. is primarily a nocturnal phenomenon that virtually vanishes during the daylight hours, it is one of the most persistent and stable climatological features of the low-level continental flow during the warm-season months, May through August. We have used significant-level data to validate the skill of the GEOS-1 Data Assimilation System (DAS) in realistically detecting this jet and inferring its structure and evolution. We have then carried out a 15-year reanalysis with the GEOS-1 DAS to determine its climatology and mean diurnal cycle and to study its interannual variability. Interannual anomalies of the meridional flow associated with the GPLLJ are much smaller than the mean diurnal fluctuations, than random intraseasonal anomalies, and than the mean wind itself. There are three maxima of low-level meridional flow variance over the Great Plains and the Gulf of Mexico: a 1.2 m2 s-2 peak over the southeast Texas, to the east and south of the mean velocity peak, a 1.0 m2 s-2 peak over the western Gulf of Mexico, and a .8 m2 s-2 peak over the upper Great Plains (UGP), near the Nebraska/South Dakota border. Each of the three variance maxima corresponds to a spatially coherent, jet-like pattern of low-level flow interannual variability. There are also three dominant modes of interannual variability corresponding to the three variance maxima, but not in a simple one-to-one relationship. Cross-sectional profiles of mean southerly wind over Texas remain relatively stable and recognizable from year to year with only its eastward flank showing significant variability. This variability, however, exhibits a distinct, biennial oscillation during the first six to seven years of the reanalysis period and only then. This intermittent biennial oscillation (IBO, one of the three modes discussed in the previous paragraph) in the lowlevel flow is restricted to the region surrounding eastern Texas and is also evident in the NCEP/NCAR reanalysis data set from about 1978 to 1985 or 1986 and again from 1995 to 2000. It is evident as well in surface pressure in both the GEOS-1 and NCEP/NCAR sets. The interannual anomalies do not necessarily persist uniformly throughout an entire season, but can fluctuate from one part of the season to the next. To estimate the characteristic sub-seasonal time scales for coherence of these fluctuations, we have taken the weekly anomaly of low-level wind at each point of the domain from the climatological average for that given point and that given week of the season and computed the covariance of its fluctuations over all weeks and over all years with the weekly climatological anomaly of the meridional wind at each of the three reference points discussed above. The typical duration of a coherent interannual anomaly within a given warm season increases with decreasing latitude from 2 to 3 weeks over the UGP, to 6 to 7 weeks over eastern Texas. Coherence over the western Gulf of Mexico is intermediate between the two with a typical duration of 4 to 5 weeks. There appears to be evidence that the interannual anomalies over Texas the Gulf propagate to the UGP after a week and those over the Gulf propagate there after 2 to 3 weeks. There also appears to be some reverse propagation of interannual anomalies over the UGP to Texas and to the Gulf after a period of about one week. The interannual anomalies in southerly flow over eastern Texas seem to correlate well with interannual anomalies of surface temperature and (negative) ground wetness and over western Texas.

Helfand, H. M.↗
DOE OSTI2020

Arctic Tipping Points Triggering Global Change (LDRD Final Report)

The Arctic is warming and feedbacks in the coupled Earth system may be driving the Arctic to tipping events that could have critical downstream impacts for the rest of the globe. In this project we have focused on analyzing sea ice variability and loss in the coupled Earth system Summer sea ice loss is happening rapidly and although the loss may be smooth and reversible, it has significant consequences for other Arctic systems as well as geopolitical and economic implications. Accurate seasonal predictions of sea ice minimum extent and long-term estimates of timing for a seasonally ice-free Arctic depend on a better understanding of the factors influencing sea ice dynamics and variation in this strongly coupled system. Under this project we have investigated the most influential factors in accurate predictions of September Arctic sea ice extent using machine learning models trained separately on observational data and on simulation data from five E3SM historical ensembles. Monthly averaged data from June, July, and August for a selection of ice, ocean, and atmosphere variables were used to train a random forest regression model. Gini importance measures were computed for each input feature with the testing data. We found that sea ice volume is most important earlier in the season (June) and sea ice extent became a more important predictor closer to September. Results from this study provide insight into how feature importance changes with forecast length and illustrates differences between observational data and simulated Earth system data. We have additionally performed a global sensitivity analysis (GSA) using a fully coupled ultra- low resolution configuration E3SM. To our knowledge, this is the first global sensitivity analysis involving the fully-coupled E3SM Earth system model. We have found that parameter variations show significant impact on the Arctic climate state and atmospheric parameters related to cloud parameterizations are the most significant. We also find significant interactions between parameters from different components of E3SM. The results of this study provide invaluable insight into the relative importance of various parameters from the sea ice, atmosphere and ocean components of the E3SM (including cross-component parameter interactions) on various Arctic-focused quantities of interest (QOIs).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Analysis of plasma measurements for the Geotail mission

The CPI plasma measurements from the Geotail spacecraft are currently used by a number of scientists in support of varied research projects. Measurements from the CPI hot plasma analyzer have been processed to compute one-minute averages of plasma densities, temperatures, and velocities for the period 1 Oct. 1992 - 31 July 1993. These parameters are used in a preliminary survey of the magnetotail for distances from earth of 10 to 210 earth radii. The cold plasmas within the magnetotail drift towards the midplane and are thought to be a principal source for the hot plasma sheet. A remarkable result from Geotail is the observation of cold ion beams coexisting as distinct components in the presence of hot plasma-sheet plasmas. The development and evolution of plasmoids is a topic of considerable interest in studies of the magnetotail and magnetospheric substorms. A number of possible plasmoids have been identified in the Geotail data set based upon reversals of the Z component of the magnetic field as observed with the MGF instrumentation. An intensive study of the electron and ion velocity distributions and the plasma parameters for several of these events reveals unexpected features. Collaborative work with several Geotail instrument teams and other researchers is also ongoing.

Frank, Louis A.↗
DOE OSTI2022

Reverse-mode differentiation in arbitrary tensor network format: with application to supervised learning.

This paper describes an efficient reverse-mode differentiation algorithm for contraction operations of tensor networks that may have arbitrary and unconventional network topologies. The approach leverages the tensor contraction tree of Evenbly and Pfeifer (2014), which provides an instruction set for the contraction sequence of a network. We show that this tree can be efficiently leveraged for differentiation of a full tensor network contraction using a recursive scheme that exploits (1) the bilinear property of contraction and (2) the property that trees have single path from root to leaves. While differentiation of tensor-tensor contraction is already possible in most automatic differentiation packages, we show that exploiting these two additional properties in the specific context of contraction sequences can improve efficiency. Following a description of the algorithm and computational complexity analysis, we investigate its utility for gradient-based supervised learning for low-rank function recovery and for fitting real-world unstructured datasets. We demonstrate improved performance over alternating least-squares optimization approaches and the capability to handle heterogeneous and arbitrary tensor network formats. When compared to alternating minimization algorithms, we find that the gradient-based approach requires a smaller oversampling ratio (number of samples compared to number model parameters) for recovery. This increased efficiency extends to fitting unstructured data of varying dimensionality and when employing a variety of tensor network formats. Here, we show improved learning using the hierarchical Tucker method over the tensor-train in high-dimensional settings on a number of benchmark problems.

97 MATHEMATICS AND COMPUTING↗