Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization methods”

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

Wormhole calculus, replicas, and entropies

We investigate contributions of spacetime wormholes, describing baby universe emission and absorption, to calculations of entropies and correlation functions, for example those based on the replica method. We find that the rules of the “wormhole calculus”, developed in the 1980s, together with standard quantum mechanical prescriptions for computing entropies and correlators, imply definite rules for limited patterns of connection between replica factors in simple calculations. These results stand in contrast with assumptions that all topologies connecting replicas should be summed over, and call into question the explanation for the latter. In a “free” approximation baby universes introduce probability distributions for coupling constants, and we review and extend arguments that successive experiments in a “parent” universe increasingly precisely fix such couplings, resulting in ultimately pure evolution. Once this has happened, the nontrivial question remains of how topology-changing effects can modify the standard description of black hole information loss.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Uncertain role of clouds in shaping summertime atmosphere-sea ice connections in reanalyses and CMIP6 models

Downwelling longwave radiation (DLR) driven by the atmospheric and cloud conditions in the troposphere is suggested to be a dominant factor to determine the summertime net surface energy budget over the Arctic Ocean and thus plays a key role to shape the September sea ice. We use reanalyses and the self-organizing map (SOM) method to distinguish CMIP6 model performance in replicating the observed strong atmosphere-DLR connection. We find all models can reasonably simulate the linkage between key atmosphere variables and the clear sky DLR but behave differently in replicating the atmosphere-DLR connection due to cloud forcing. In ERA5 and strongly coupled models, tropospheric high pressure is associated with decreased clouds in the mid- and high-levels and increased clouds near the surface. This out-of-phase structure indicates that DLR cloud forcing is nearly neutral, making the clear sky DLR more important to bridge JJA circulation to late-summer sea ice. In MERRA-2 and weakly coupled models, tropospheric clouds display a vertically homogeneous reduction; the cloud DLR is thus strongly reduced due to the cooling effect, which partially cancels out the clear sky DLR and makes the total DLR less efficient to translate circulation forcing to sea ice. The differences of cloud vertical distribution in CMIP6 appear to be differentiated by circulation related relative humidity. Therefore, a better understanding of the discrepancy of different reanalyses and remote sensing products is critical to comprehensively evaluate simulated interactions among circulation, clouds, sea ice and energy budget at the surface in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Uranium hydride corrosion. II Modeling spatially random growth using reaction and nucleation kinetics

A uranium hydride (UH 3 ) corrosion kinetics model has been developed, based on dimensionally restricted nucleation kinetics, which combines heretofore separate reaction and nucleation kinetics models into a single predictive construct. This theoretical framework accounts for the superficial and random pitting behaviour of UH 3 corrosion. The model can be used to generate the net corrosion surface area and volume fractions for a uranium specimen, generally accounting for its geometric form factor through the surface-to-volume ratio. Furthermore, the combined corrosion kinetics model enables the generation of quantities which may be directly compared to those obtained by non-destructive, non-contact evaluation methods.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Model predictive control for optimal dispatch of chillers and thermal energy storage tank in airports

Cost of energy consumption is one of the biggest operational cost for airports, and it is increasing from time to time as airports expand to support growing number of passengers. Various factors affect the energy consumption including efficiency of airport Heating Ventilation and Air conditioning (HVAC) systems, which in turn depends on the efficiency of individual subsystems. Here, in this paper, we present an optimal scheduling method for the central plant system at Dallas Fort Worth airport, involving chillers, pumps, and a thermal energy storage (TES) system. A model predictive control (MPC) problem is formulated to minimize both energy and demand charge costs while satisfying the cooling needs of the airport. The proposed Mixed-Integer Nonlinear Programming (MINLP) formulation includes performance curve based models for chillers and pumps and a simplified state of charge model for TES. The formulation also includes predictions of cooling load and chilled water return temperature. Simulation results for a month in summer show savings around 10% compared to the baseline. Initial recommendations based on insights from simulation results to the manual operation procedures resulted in significant savings. Field test results show a 7% chiller efficiency improvement.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Degradation of solid oxide electrolysis cells: Phenomena, mechanisms, and emerging mitigation strategies—A review

Solid oxide electrolysis cell (SOEC) is a promising electrochemical device with high efficiency for energy storage and conversion. However, the degradation of SOEC is a significant barrier to commercial viability. In this review article, the typical degradation phenomena of SOEC are summarized, with great attention into the anodes/oxygen electrodes, including the commonly used and newly developed anode materials. Meanwhile, mechanistic investigations on the electrode/electrolyte interfaces are provided to unveil how the intrinsic factor, oxygen partial pressure , and the electrochemical operation conditions, affect the interfacial stability of SOEC. At last, this paper also presents some emerging mitigation strategies to circumvent long-term degradation, which include novel infiltration method, development of new anode materials and engineering of the microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Water treatment based on atomically engineered materials: Atomic layer deposition and beyond

Global water stress and challenges for producing sufficient supplies of fit-for-purpose water are amplifying. Atomically engineered interfaces are emerging as a powerful tool in the fabrication of advanced water treatment materials. Atomic layer deposition (ALD) and recently developed related methods, such as sequential infiltration synthesis (SIS), offer a tremendously diverse library of chemistries for interface functionalization. Thickness, stoichiometry, and physicochemical properties can be manipulated with precision. We review their fundamental physical chemistry and processing factors. ALD/SIS engineering strategies, including direct deposition, growth with intermediate layers, and secondary treatment are presented with realization of efficient water treatment. We lay out a pathway to establishing an ALD/SIS-based universal functionalization platform for water treatment, including sensitization strategies, in situ regulation, secondary reactions, and simulation/machine learning. Furthermore, we also provide a perspective on ALD/SIS-based interface engineering via synergy with other widely used interface engineering techniques to develop facile, versatile, and energy-efficient strategies for tackling increasingly complex water challenges.

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

Highly Rectifying Fluidic Diodes Based on Asymmetric Layer-by-Layer Nanofilms on Nanochannel Membranes

Nanochannel-based fluidic diodes display ion selectivity and ion-current rectification (ICR), which may prove important in energy harvesting and biosensors. This paper reports asymmetric functionalization of the outer surface of a flexible nanochannel polymer membrane to create fluidic diodes that give ICR. Layer-by-layer (LbL) adsorption with cross-linking of only the underlying part of the polyelectrolyte nanofilm leads to a porosity step across the film. The combination of a high effective surface-charge density and the porosity step in the film leads to a remarkable maximum ICR factor of ~200 with a pH gradient across the film. Incorporation of pH-sensitive polyelectrolyte components enables the ICR factor to increase an order of magnitude on going from pH 8 to pH 3. Furthermore, the coated membrane shows excellent anion-selectivity. Thus, LbL adsorption with partial cross-linking provides a simple method for creating coated nanochannel membranes that serve as pH-responsive ionic diodes for potential chemical/biosensors.

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

Dispersive and nondispersive 𝐾-matrix formalisms

The modeling of coupled-channel effects has become increasingly important due to the availability of highly precise data for a large variety of hadronic (re)scattering processes. The 𝐾-matrix is a powerful, yet comparatively simple, method to describe scattering amplitudes, including coupled-channel effects, with the aim of interpreting experimental data. Throughout the literature, a range of dispersive and nondispersive 𝐾-matrix methods are employed. Here, we compare the dispersive and nondispersive formulations in the context of the N/D method. It is shown that the methods are equivalent in the physical region under 𝐾-matrix reparametrization. Differences away from the physical region are examined. Applications to synthetic data are used to illustrate the effects of model choices concerning form factors and the application of dispersion relations, with the goal of clarifying best practices. We find no clear preference with regard to dispersive modeling. In contrast, we find that interpretational ambiguity of the bare model parameters—and even of the form of the bare model—is endemic, and recommend a thorough sampling of data and model spaces to assess conclusion robustness.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Chemical and elemental mapping of spent nuclear fuel sections by soft X-ray spectromicroscopy

Soft X-ray spectromicroscopy at the O K -edge, U N 4,5 -edges and Ce M 4,5 -edges has been performed on focused ion beam sections of spent nuclear fuel for the first time, yielding chemical information on the sub-micrometer scale. To analyze these data, a modification to non-negative matrix factorization (NMF) was developed, in which the data are no longer required to be non-negative, but the non-negativity of the spectral components and fit coefficients is largely preserved. The modified NMF method was utilized at the O K -edge to distinguish between two components, one present in the bulk of the sample similar to UO 2 and one present at the interface of the sample which is a hyperstoichiometric UO 2+ x species. The species maps are consistent with a model of a thin layer of UO 2+ x over the entire sample, which is likely explained by oxidation after focused ion beam (FIB) sectioning. In addition to the uranium oxide bulk of the sample, Ce measurements were also performed to investigate the oxidation state of that fission product, which is the subject of considerable interest. Analysis of the Ce spectra shows that Ce is in a predominantly trivalent state, with a possible contribution from tetravalent Ce. Atom probe analysis was performed to provide confirmation of the presence and localization of Ce in the spent fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

TULIP: An RNA-seq-based Primary Tumor Type Prediction Tool Using Convolutional Neural Networks

Background: With cancer as one of the leading causes of death worldwide, accurate primary tumor type prediction is critical in identifying genetic factors that can inhibit or slow tumor progression. There have been efforts to categorize primary tumor types with gene expression data using machine learning, and more recently with deep learning, in the last several years. Methods In this paper, we developed four 1-dimensional (1D) Convolutional Neural Network (CNN) models to classify RNA-seq count data as one of 17 highly represented primary tumor types or 32 primary tumor types regardless of imbalanced representation. Additionally, we adapted the models to take as input either all Ensembl genes (60,483) or protein coding genes only (19,758). Unlike previous work, we avoided selection bias by not filtering genes based on expression values. RNA-seq count data expressed as FPKM-UQ of 9,025 and 10,940 samples from The Cancer Genome Atlas (TCGA) were downloaded from the Genomic Data Commons (GDC) corresponding to 17 and 32 primary tumor types respectively for training and validating the models. Results: All 4 1D-CNN models had an overall accuracy of 94.7% to 97.6% on the test dataset. Further evaluation indicates that the models with protein coding genes only as features performed with better accuracy compared to the models with all Ensembl genes for both 17 and 32 primary tumor types. For all models, the accuracy by primary tumor type was above 80% for most primary tumor types. Conclusions: We packaged all 4 models as a Python-based deep learning classification tool called TULIP (TUmor CLassIfication Predictor) for performing quality control on primary tumor samples and characterizing cancer samples of unknown tumor type. Further optimization of the models is needed to improve the accuracy of certain primary tumor types.

Jones, Sara↗
DOE OSTI2020

Outgrowth of erlotinib-resistant subpopulations recapitulated in patient-derived lung tumor spheroids and organoids

A model that recapitulates development of acquired therapeutic resistance is needed to improve oncology drug development and patient outcomes. To achieve this end, we established methods for the preparation and growth of spheroids from primary human lung adenocarcinomas, including methods to culture, passage, monitor growth, and evaluate changes in mutational profile over time. Primary lung tumor spheroids were cultured in Matrigel® with varying concentrations of erlotinib, a small molecule kinase inhibitor of epidermal growth factor receptor (EGFR) that is ineffective against KRAS mutant cells. Subtle changes in spheroid size and number were observed within the first two weeks of culture. Spheroids were cultured for up to 24 weeks, during which time interactions between different cell types, movement, and assembly into heterogeneous organoid structures were documented. Allele-specific competitive blocker PCR (ACB-PCR) was used to quantify low frequency BRAF V600E, KRAS G12D, KRAS G12V, and PIK3CA H1047R mutant subpopulations in tumor tissue residue (TR) samples and cultured spheroids. Mutant subpopulations, including multiple mutant subpopulations, were quite prevalent. Twelve examples of mutant enrichment were found in eight of the 14 tumors analyzed, based on the criteria that a statistically-significant increase in mutant fraction was observed relative to both the TR and the no-erlotinib control. Of the mutants quantified in erlotinib-treated cultures, PIK3CA H1047 mutant subpopulations increased most often (5/14 tumors), which is consistent with clinical observations. Thus, this ex vivo lung tumor spheroid model replicates the cellular and mutational tumor heterogeneity of human lung adenocarcinomas and can be used to assess the outgrowth of mutant subpopulations. Spheroid cultures with characterized mutant subpopulations could be used to investigate the efficacy of lung cancer combination therapies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Utilizing Sensitivity and Correlation Coefficients from MCNP and Whisper to Guide Microreactor Experiment Design

When designing experiments for full-scale reactor systems, MCNP®* and Whisper can be used to create neutronic models and compare the similarity of two nuclear systems via correlation coefficients for κ eff , effective multiplication factor. This thesis applies this framework to a conceptual heat-pipe, yttrium-hydride moderated microreactor system and experiments. The framework is intended as a supplement to other neutronics/thermal/multiphysics analyses and provides a concrete method to measure the neutronic similarity of two systems. By analyzing the shared nuclear data uncertainty, as well as sensitivity to nuclear data over all neutron energies, highly informative experiments can be designed to aid in the development of microreactor and other advanced reactor technologies and systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Barriers and Opportunities for Energy Technology Adoption in Juneau, Alaska

This report presents findings from a qualitative study examining barriers and opportunities for air source heat pump (ASHP) and electric vehicle (EV) adoption in Juneau, Alaska, with a particular focus on manufactured and multifamily housing. The analysis draws on community insights from end users and middle actors to better understand how technology adoption unfolds in contexts with distinct logistical, infrastructural, and housing constraints. The report is organized according to key barriers and opportunities identified through stakeholder input, providing a structured understanding of adoption dynamics across technologies and housing types. These insights are intended to inform program design and support more effective electrification strategies tailored to local conditions. The study team employed qualitative methods to capture both in-depth and high-level perspectives on technology adoption. Data collection included: 1) two 2-hour focus groups with a total of five end users and seven middle actors, enabling detailed and structured discussion and 2) ten semistructured interviews with manufactured home owners, multifamily landlords, and one tenant, providing complementary insights across housing contexts. Focus groups captured accounts of shared challenges and opportunities while interviews offered more concise reflections on individual experiences. Together, these methods enabled a more comprehensive understanding of both systemic barriers and lived experiences with ASHPs and EVs. The findings reveal that adoption of electrification technologies is shaped by a combination of economic, logistical, and informational factors that vary across housing types, technology characteristics, user groups, and other demographic factors. Addressing these factors requires tailored strategies that reflect local conditions and user experiences. The insights in this report can provide a foundation for organizations such as AEL&P to refine program design, support more effective outreach, and anticipate shifts in energy demand associated with increased electrification. More broadly, the study highlights the importance of incorporating community perspectives when developing electrification initiatives to ensure they are both practical and responsive to real-world constraints.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Unbalanced Nested Random Effects Estimation of Variance Components

To investigate the contributing factors of variance in the measurement of an iodine 127 (127I) sample, we implement a nested random effects analysis of variance (ANOVA). Historically, the reported uncertainty on a measurement of 127I has been obtained by methods of forward uncertainty propagation because there is typically only one replicate of a given sample for which to estimate the uncertainty. When samples are processed there are several types of quality control (QC) standards analyzed along-side the unknown samples with two to five replicates for each. Assuming the variance observed in these replicate QC standards is representative of that of the unknown samples, we use these data in a nested random effects ANOVA to estimate the total uncertainty of a measurement. We demonstrate this approach with two sets of measurements from Idaho National Laboratory and compare the results with the forward uncertainty propagation approach. Variance component estimates for the coarsest level of the nesting structure were most imprecise because replicates were most limited at these levels. We find that the results largely agree between forward propagation and ANOVA, and the greatest contributors of variance are due to instrument variation and chemical processing, with human processing being among the smallest contributors. This analysis provides reassurance that the reported uncertainties using forward propagation are reasonable and the process is well controlled. We propose a future designed experiment to increase replicates at the coarsest level of the hierarchy to improve estimates of these variance components.

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

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR LiPON

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu12-xMxSb4S13, where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m-1 K-1) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive.In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu10.5NiZn0.5Sb4S13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu14Sb4S13 is revealed by Cu self-diffusivity, nuclear density map and "nudged elastic band" (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu10Zn2Sb4S13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu12Sb4S13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li3PO4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability

25 ENERGY STORAGE↗
DOE OSTI2022

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗
DOE OSTI2023

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗
DOE OSTI2020

Risk-Informed Applications and Science – Past, Present, & Future

This is a presentation to be made at the University of Tokyo. The Idaho National Laboratory is creating the next-generation reliability- and risk-assessment methods and tools that support risk-informed decision-making by combining physics-based models with probabilistic quantification approaches. Integrating these two worlds of physics and probability using a simulation framework leads us to predictions based upon an approach called “computational risk assessment” which will serve as the technical basis for the future of reliability and risk approaches. The driving factors for this new approach includes: temporal (timing issues), spatial (location issues), mechanistic (physics issues), and topological (complexity issues). By combining phenomenology directly with stochastic quantification, we are able to perform advanced uncertainty analysis directly on both parameters and models. In addition, new tools allow for the creation of reduced order models that mimic high-fidelity engineering models while still permitting the realization of thousands of simulation iterations in short run-times on current workstation-class computers. While these advanced methods and tools can provide increased realism in our engineering safety and risk approaches, their greater benefit is to provide a risk-informed engineering framework for design and operation. This talk will briefly review the history of risk assessment and risk-informed applications, describe some of the current research and development found at the Idaho National Laboratory, and discuss potential future applications and approaches for advanced methods and tools.

97 MATHEMATICS AND COMPUTING↗