Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependence analysis”

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 37 records · Page 2

Solid-to-Fluid Radiative Heat Transfer Modeling for System Analysis Module

System Analysis Module (SAM) is a system-level thermal hydraulics code being developed at Argonne National Laboratory for advanced nuclear reactor analysis. In addition to a wide range of interests from the advanced reactor design community, SAM has also been adopted by the United States Nuclear Regulatory Commission's suite of codes purposed for advanced reactor licensing. Nevertheless, the code is still under active development and new capabilities are being added to address various modeling and simulation challenges for advanced reactor analysis. One such phenomenon important to the thermal behavior of some advanced reactor concepts is radiative heat transfer (radHT). Conditions, such as high temperatures and long optical paths, increase the radiative contributions from solids and coolants alike. This paper discusses the development of radHT modeling in SAM and describes the new capabilities provided for thermal analysis. Depending on the geometry and temperatures of the system at hand, as well as the coolant in question, thermal transfer due to thermal radiation will vary dramatically. Therefore, the ability to model variable radiative systems was maintained as a priority during development of SAM radHT modeling. This newly developed simulation feature provides a flexible solid-to-fluid radiative heat transfer framework necessary for SAM to perform accurate analysis for advanced reactor designs. Test cases are also presented and shown to match analytical solutions, which demonstrate the radHT model's efficacy.

Radiative heat transfer↗

Towards data-driven constitutive modelling for granular materials via micromechanics-informed deep learning

The analytical description of path-dependent elastic-plastic responses of a granular system is highly complicated because of continuously evolving microstructures and strain localisation within the system undergoing deformation. This study offers an alternative to the current analytical paradigm by developing micromechanics-informed machine-learning based constitutive modelling approaches for granular materials. A set of critical variables associated with the constitutive behaviour of granular materials are identified through an incremental stress-strain relationship analysis. Depending on the strategy to exploit the priori micromechanical knowledge, three different training strategies are explored. The first model uses only the measurable external variables to make stress predictions; the second model utilises a directed graph to link all the external strain sequences and internal microstructural evolution variables into a single prediction model comprised of a series of sub-mappings, and the third model explicitly integrates the physically important non-temporal properties with external strain paths into training through an enhanced Gated Recurrent Unit (GRU). These three models show satisfactory agreement with unseen test specimens based on multi-directional loading cases. The basic features and potential applications of each model are explained. Lastly, the key factors for constitutive training and limitations of the current work are also discussed in detail.

36 MATERIALS SCIENCE↗

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements of the time-dependent cosmic-ray Sun shadow with seven years of IceCube data: Comparison with the Solar cycle and magnetic field models

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. Here, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4σ level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of 3σ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-γ Experiment; a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

79 ASTRONOMY AND ASTROPHYSICS↗

Bimetallic Metal–Organic Framework Fe/Co-MIL-88(NH 2 ) Exhibiting High Peroxidase-like Activity and Its Application in Detection of Extracellular Vesicles

Metal–organic frameworks (MOFs) have many attractive features, including tunable composition, rigid structure, controllable pore size, and large specific surface area, and thus are highly applicable in molecular analysis. Depending on the MOF structure, a high number of un-saturated metal sites can be exposed to catalyze chemical reactions. In the present work, we report that by using both Co(II) and Fe(III) to prepare the MIL-88(NH 2 ) MOF, we can produce the bimetallic MOF that can catalyze the conversion of 3,3', 5,5"-tetramethylbenzidine (TMB) to a color product through reaction with H 2 O 2 at a higher reaction rate than the monometallic Fe-MIL-88(NH 2 ). The Michaelis constants (K m ) of the catalytic reaction for TMB and H 2 O 2 are 3-5 times smaller, and the catalytic constants (k cat ) are 5-10 times higher than those of the horse-radish peroxidase (HRP), supporting ultrahigh peroxidase-like activity. These values are also much more superior to those of the HRP-mimicking MOFs reported previously. Interestingly, the bimetallic MOF can be coupled with glucose oxidase (GOx) to trigger the cascade enzymatic reaction for highly sensitive detection of extracellular vesicles (EVs), a family of important biomarkers. Through conjugation to the aptamer that recognizes the marker protein on EV surface, the MOF can help isolate the EVs from biological matrices, which are subsequently labeled by GOx via antibody recognition. The cascade enzymatic reaction between MOF and GOx enables detection of EVs at a concentration as low as 7.8 × 10 4 particles/ml. Further, the assay can be applied to monitor EV secretion by cultured cells, and also can successfully detect the different EV quantities in the sera samples collected from cancer patients and healthy controls. Overall, we prove that the bimetallic Fe/Co-MIL-88(NH 2 ) MOF, with its high peroxidase-activity and high biocompatibility, is a valuable tool deployable in clinical assays to facility disease diagnosis and prognosis.

36 MATERIALS SCIENCE↗

In Situ Diffraction and Ex Situ Transmission X‐Ray Microscopy Studies of Solid‐State Upcycling for NMC Cathodes

Upcycling of recycled LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) cathodes offers an economical route to produce cathode materials with increased energy density (i.e., LiNi 0.8 Mn 0.1 Co 0.1 O 2 , NMC811) that meet the performance needs of present-day electric vehicles. In this work, solid-state upcycling of NMC622 via calcination with Ni(OH) 2 and LiOH was monitored using in situ synchrotron powder X-ray diffraction measurements. Sequential Rietveld refinements indicate that the calcination proceeds by initially converting Ni(OH) 2 to a rocksalt NiO phase followed by lithiation of NiO to form LiNiO 2 (LNO), with both NMC and LNO phases present in nearly equal proportions at the calcination endpoint. Variable-energy transmission X-ray microscopy tomograms of upcycled samples reveal that the NMC and LNO domains are intermixed at sub-micron length scales. Depth-dependent analysis of multi-elemental fitting maps matches the expected NMC811 composition at the secondary particle level and indicates that transition metal diffusion is not limited by the secondary particle size.

cathode upcycling↗

Machine learning elastic constants of multi-component alloys

The present manuscript explores application of machine learning methods for determining elastic constants and other derived mechanical properties of multi-component alloys. Here, a number of machine learning models, including linear regression, neural network and random forest based models, are trained and tested on a dataset of binary alloys generated using density functional theory (DFT) calculations and spanning over a large number of elemental species in the periodic table. Starting with a wide range of simple and easily accessible compositionally-averaged elemental features, a correlation-based feature selection strategy was used to systematically down-select a set of most relevant features towards the prediction of the elasticity tensor components. The true predictive performance and the associated uncertainties of the models were established by testing on unseen data and bootstrapping, respectively. A single and pair-wise feature partial dependence analysis was performed to visualize the average property trends in the multi-dimensional feature space in order to further understand the achieved predictive performance. The utility of the trained model is further demonstrated by obtaining sufficiently accurate yet highly efficient approximations for bulk modulus, Young’s modulus, shear modulus and Poisson’s ratio for alloys beyond the binary space (i.e., two-component alloys) on which the model was originally trained. More importantly, we test and validate the predictive performance of the developed model directly against the experimentally measured elastic constants of technologically relevant multi-component alloys (such as, Ni- and Ti-based alloys). Finally, utility of such a data-enabled route is demonstrated by predicting the possible range of various elastic properties for vast composition space available within the five component Ni-Cr-Fe-Mo-W alloy system in a high-throughput manner.

36 MATERIALS SCIENCE↗

Spectra-Orthogonal Optical Anisotropy in Wafer-Scale Molecular Crystal Monolayers

Controlling the spectral and polarization response of two-dimensional (2D) crystals is vital for developing ultrathin platforms for compact optoelectronic devices. However, independently tuning optical anisotropy and spectral response remains challenging in conventional semiconductors due to the intertwined nature of their lattice and electronic structures. Here, we report spectra-orthogonal optical anisotropy─where polarization anisotropy is tuned independently of spectral response─in wafer-scale, one-atom-thick 2D molecular crystal (2DMC) monolayers synthesized on monolayer transition-metal dichalcogenide (TMD) crystals. Utilizing the concomitant spectral consistency and structural tunability of perylene derivatives, we demonstrate tunable optical polarization anisotropy in 2DMCs with similar spectral profiles, as confirmed by room-temperature scanning tunneling microscopy and cross-polarized reflectance microscopy. Additional angle-dependent analysis of the single-crystal and polycrystalline molecular domains reveals an epitaxial relationship between the 2DMC and TMD. In conclusion, our results establish a scalable, molecule-based 2D crystalline platform for unique and tunable functionalities unattainable in covalent 2D solids.

2D materials↗

Light induced ion migration studies in perovskite solar cell using nonlinear impedance spectroscopy

Complex interactions between mobile ions and charge carriers in perovskite solar cells (PSCs) make it challenging to fully understand their dynamic interplay. Exposure to light further complicates these interactions, altering the system’s dynamics and inducing nonlinear effects that lead to changes in the J−V curve. Understanding these effects is crucial for improving the operational stability of PSCs. Impedance spectroscopy (IS) is a powerful technique for evaluating relaxation processes in the frequency domain; however, it is limited in capturing nonlinear contributions. Here, in this work, nonlinear impedance spectroscopy (NLIS) is employed to analyze the higher harmonic response to AC perturbation, both in the dark and after short-term light exposure. A shift in the low-frequency (LF) higher harmonic peak is observed after open-circuit light exposure, attributed to an altered electric field suggesting ion re-distribution, whereas closed-circuit exposure shows no LF shift, indicating minimal ion movement. Additionally, light exposure reduces higher-order admittance, more notably in open-circuit conditions, suggesting decreased recombination. Temperature-dependent analysis was conducted to characterize the activation energy of migrating species, identifying iodide as the dominant migrating ion.

14 SOLAR ENERGY↗

Autonomous Nanoparticle Synthesis Guided by In Situ Multiscale Structural Characterization

Autonomous synthesis platforms promise rapid exploration of vast parameter spaces; yet, integrating in situ structural characterization in closed-loop synthesis optimization remains challenging. We demonstrate a realization of such a closed-loop platform coupled with a droplet-flow microreactor, in situ X-ray scattering methods (SAXS/WAXS), and Gaussian process optimization to synthesize citrate-reduced Au nanoparticles with targeted characteristics. The system efficiently explored ∼19,000 synthesis recipes through 365 experiments, achieving precise control over size (4–60 nm) and polydispersity (σ < 0.11) across large citrate/gold ratios, exceeding traditional synthesis boundaries (1–10). Beyond confirming classical Turkevich–Frens trends, partial-dependence analysis revealed strong nonlinear coupling among precursor, citrate, and pH effects. Combining quantitative SAXS/WAXS analysis with electron microscopy characterization, we uncovered that crystallite size (d c ) and particle size (d) follow d c = 0.18d + β, where synthesis chemistry controls the intercept β while maintaining a universal slope. This parallel-band structure enables independent tuning of crystallite domain size at fixed particle diameter through a combination of chloride, gold precursor, citrate, and pH contributions (cross-validated Spearman ρ = 0.7 ± 0.1). High-resolution electron microscopy shows multiple lattice-fringe orientations within single particles, directly confirming polycrystalline domains and the ability to tune d c at the fixed d. The platform’s validation includes indistinguishable static versus flowing measurements, stable droplet transport at 100 °C, and <5% run-to-run variation, establishing a robust framework for mapping and controlling multiscale nanoparticle structure across expansive chemical spaces. In conclusion, the developed closed-loop platform can be applied to a borad range of nanosyntheis processes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

The curious case of the structural phase transition in SnSe insights from neutron total scattering

Abstract At elevated temperatures SnSe is reported to undergo a structural transition from the low symmetry orthorhombic GeS-type to a higher symmetry orthorhombic TlI-type. Although increasing symmetry should likewise increase lattice thermal conductivity, many experiments on single crystals and polycrystalline materials indicate that this is not the case. Here we present temperature dependent analysis of time-of-flight (TOF) neutron total scattering data in combination with theoretical modeling to probe the local to long-range evolution of the structure. We report that while SnSe is well characterized on average within the high symmetry space group above the transition, over length scales of a few unit cells SnSe remains better characterized in the low symmetry GeS-type space group. Our finding from robust modeling provides further insight into the curious case of a dynamic order-disorder phase transition in SnSe, a model consistent with the soft-phonon picture of the high thermoelectric power above the phase transition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interpretable machine learning-guided design of Fe-based soft magnetic alloys

Here, we present a machine learning (ML) guided approach to predict saturation magnetization (𝑀 S ) and coercivity (𝐻 C ) in Fe-rich soft magnetic alloys, particularly Fe-Si-B systems. ML models trained on experimental data reveal that increasing Si and B content reduces 𝑀 S from 1.81 T (DFT ≈ 2.04 T) to ≈1.54 T (DFT ≈ 1.56T) in Fe-Si-B, which is attributed to decreased magnetic density and structural modifications. Experimental validation of ML predicted magnetic saturation on Fe-1Si-1B (2.09 T), Fe-5Si-5B (2.01 T), and Fe-10Si-10B (1.54 T) alloy compositions further supports our findings. These trends are consistent with density functional theory predictions, which link increased electronic disorder and band broadening to lower 𝑀 S values. Experimental validation on selected alloys confirms the predictive accuracy of the ML model, with good agreement across compositions. Beyond predictive accuracy, detailed uncertainty quantification and model interpretability including through feature importance and partial dependence analysis reveal that 𝑀 S is governed by a nonlinear interplay between Fe content and early transition metal ratios, while 𝐻 C is more sensitive to processing conditions such as ribbon thickness and thermal treatment windows. The ML framework was further applied to Fe-Si-B/Cr/Cu/Zr/Nb alloys in a pseudoquaternary compositional space, which shows comparable magnetic properties to NANOMET (Fe 84.8 ⁢Si 0.5 ⁢B 9.4 ⁢Cu 0.8⁢ P 3.5 ⁢C 1 ), FINEMET (Fe 73.5 ⁢Si 13.5 ⁢B 9 Cu 1 ⁢Nb 3 ), NANOPERM (Fe 88 ⁢Zr 7⁢ B 4 ⁢Cu 1 ), and HITPERM (Fe 44 ⁢Co 44 ⁢Zr 7⁢ B 4 ⁢Cu 1 . Our findings demonstrate the potential of the ML framework for accelerated search of high-performance soft magnetic materials.

density functional theory↗

Benchmarking Variables for Checkpointing in HPC Applications

Checkpoint/Restart (C/R) is a widely used fault tolerance mechanism in converged systems of cloud, edge, and HPC. However, users often rely on their experience to determine which variables to checkpoint, as there is currently no benchmark that can provide a reference. This can result in checkpointing redundant or even incorrect variables. To address this issue, we propose a benchmark suite that includes critical variables for checkpointing, which have been manually identified, and a method for identifying those critical variables, with 20 representative HPC applications. Our method involves analyzing data dependency between variables to identify critical variables analytically. We verify the identified variables' correctness with a widely used C/R library FTI by an ablation study. With our benchmark suite and data dependency analysis, HPC practitioners now have a reference for identifying checkpointing variables and better knowledge of what kind of variables to checkpoint.

Fu, Xiang↗

Ambipolar Transport in Polycrystalline GeSn Transistors for Complementary Metal-Oxide-Semiconductor Applications

Group-IV alloy GeSn is a promising material for electronic and optoelectronic applications due to its compatibility with both Si substrates and established Si fabrication processes. This study focuses on polycrystalline GeSn (10% Sn), which offers a cost-effective, large-area, and versatile alternative to epitaxial GeSn. We demonstrate ambipolar transport behavior in polycrystalline GeSn thin film transistors, achieving electron and hole field-effect mobilities reaching up to 0.05 cm 2 /Vs and 2.05 cm 2 /Vs, respectively. Through temperature-dependent analysis, we elucidate the underlying mechanism of this phenomenon, which we attribute to quantum tunneling between the Schottky barrier contact and the channel, as well as potential barriers between the grain boundaries of this polycrystalline film, thereby advancing the understanding of polycrystalline GeSn's electrical properties. Furthermore, this work highlights the potential of ambipolar transport as a technique to employ towards the development of GeSn complementary metal-oxide-semiconductor field-effect transistors, promising to simplify and reduce the cost of GeSn manufacturing processes for edge computing and sensing applications.

42 ENGINEERING↗

Measurement of the Higgs boson production rate in association with top quarks in final states with electrons, muons, and hadronically decaying tau leptons at $\sqrt{s} =$ 13 TeV

The rate for Higgs (${\mathrm{H}} $) bosons production in association with either one (${\mathrm{t}} {\mathrm{H}} $) or two (${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $) top quarks is measured in final states containing multiple electrons, muons, or tau leptons decaying to hadrons and a neutrino, using proton–proton collisions recorded at a center-of-mass energy of $13\,\text {Te}\text {V} $ by the CMS experiment. The analyzed data correspond to an integrated luminosity of 137$\,\text {fb}^{-1}$. The analysis is aimed at events that contain ${\mathrm{H}} \rightarrow {\mathrm{W}} {\mathrm{W}} $, ${\mathrm{H}} \rightarrow {\tau } {\tau } $, or ${\mathrm{H}} \rightarrow {\mathrm{Z}} {\mathrm{Z}} $ decays and each of the top quark(s) decays either to lepton+jets or all-jet channels. Sensitivity to signal is maximized by including ten signatures in the analysis, depending on the lepton multiplicity. The separation among ${\mathrm{t}} {\mathrm{H}} $, ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and the backgrounds is enhanced through machine-learning techniques and matrix-element methods. The measured production rates for the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ and ${\mathrm{t}} {\mathrm{H}} $ signals correspond to $0.92 \pm 0.19\,\text {(stat)} ^{+0.17}_{-0.13}\,\text {(syst)} $ and $5.7 \pm 2.7\,\text {(stat)} \pm 3.0\,\text {(syst)} $ of their respective standard model (SM) expectations. The corresponding observed (expected) significance amounts to 4.7 (5.2) standard deviations for ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and to 1.4 (0.3) for ${\mathrm{t}} {\mathrm{H}} $ production. Assuming that the Higgs boson coupling to the tau lepton is equal in strength to its expectation in the SM, the coupling $y_{{\mathrm{t}}}$ of the Higgs boson to the top quark divided by its SM expectation, $\kappa _{{\mathrm{t}}}=y_{{\mathrm{t}}}/y_{{\mathrm{t}}}^{\mathrm {SM}}$, is constrained to be within $-0.9< \kappa _{{\mathrm{t}}}< -0.7$ or $0.7< \kappa _{{\mathrm{t}}}< 1.1$, at 95% confidence level. This result is the most sensitive measurement of the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ production rate to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Forward sensitivity analysis and mode dependent control for closure modeling of Galerkin systems

Model reduction by projection-based approaches is often associated with losing some of the important features that contribute towards the dynamics of the retained scales. As a result, a mismatch occurs between the predicted trajectories of the original system and the truncated one. We put forth a framework to apply a continuous time control signal in the latent space of the reduced order model (ROM) to account for the effect of truncation. We set the control input using parameterized models by following energy transfer principles. Our methodology relies on observing the system behavior in the physical space and using the projection operator to restrict the feedback signal into the latent space. Then, we leverage the forward sensitivity method (FSM) to derive relationships between the feedback and the desired mode-dependent control. We test the performance of the proposed approach using two test cases, corresponding to viscous Burgers and vortex merger problems at high Reynolds number. Results show that the ROM trajectory with the applied FSM control closely matches its target values in both the data-dense and data-sparse regimes.

97 MATHEMATICS AND COMPUTING↗