Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “laser directionality”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Lasergate: A windowless gas target for enhanced laser preheat in magnetized liner inertial fusion

At the Z Facility at Sandia National Laboratories, the magnetized liner inertial fusion (MagLIF) program aims to study the inertial confinement fusion in deuterium-filled gas cells by implementing a three-step process on the fuel: premagnetization, laser preheat, and Z-pinch compression. In the laser preheat stage, the Z-Beamlet laser focuses through a thin polyimide window to enter the gas cell and heat the fusion fuel. However, it is known that the presence of the few μm thick window reduces the amount of laser energy that enters the gas and causes window material to mix into the fuel. These effects are detrimental to achieving fusion; therefore, a windowless target is desired. The Lasergate concept is designed to accomplish this by “cutting” the window and allowing the interior gas pressure to push the window material out of the beam path just before the heating laser arrives. In this work, we present the proof-of-principle experiments to evaluate a laser-cutting approach to Lasergate and explore the subsequent window and gas dynamics. Further, an experimental comparison of gas preheat with and without Lasergate gives clear indications of an energy deposition advantage using the Lasergate concept, as well as other observed and hypothesized benefits. While Lasergate was conceived with MagLIF in mind, the method is applicable to any laser or diagnostic application requiring direct line of sight to the interior of gas cell targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Double cylinder target design for study of hydrodynamic instabilities in multi-shell ICF

Cylindrical implosions are used to study hydrodynamic instability growth for inertial confinement fusion (ICF) applications, as the cylindrical geometry allows for easier diagnostic access while retaining convergence effects. In this work, we use the established cylindrical implosion platform [Palaniyappan et al., Phys. Plasmas 27, 042708 (2020)] to inform the double shell ICF campaign [Montgomery et al., Phys. Plasmas 25, 092706 (2018)]. We present a design for a double cylindrical target as an analogue to the double shell ICF capsule in order to study hydrodynamic instability growth on the high-Z inner shell. Our design work is done with two-dimensional (2D) Eulerian radiation-hydrodynamics simulations, considering the axial uniformity of the implosion and feasibility of measuring the instability growth of pre-seeded single mode sinusoidal perturbations. We discuss in depth the design for a target to be directly driven at the OMEGA laser facility [Boehly et al., Opt. Commun. 133, 495 (1997)]. We evaluate the design for axial implosion symmetry and visibility of instability growth using synthetic radiographs constructed from the simulations, as the instability growth on the inner cylinder is experimentally measured using x-ray radiography of the implosion. We find that the seeded perturbation growth on the inner cylinder should be visible in an experiment, even with axial implosion asymmetry and preheat. We compare our 2D simulations with linear theory predictions for perturbation growth and show that a cylinder with lower azimuthal mode number (mode-20) perturbations compares more favorably with linear theory, while a cylinder with higher azimuthal mode number (mode-40) perturbations at the same starting amplitude saturates and is over-predicted by linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Spatiotemporally resolved measurements of electric field around a piezoelectric transformer using electric-field induced second harmonic (E-FISH) generation

When a piezoelectric transformer (PT) is actuated at its second harmonic frequency by a low input voltage, the generated electric field at the distal end can be sufficient to breakdown the surrounding gas, making them attractive power sources for non-equilibrium plasma generation. Understanding the potential and electric fields produced in the surrounding medium by the PT is important for effectively designing and using PT plasma devices. In this work, the spatiotemporally resolved characteristics of the electric field generated by a PT operating in open air have been investigated using the femtosecond electric field-induced second harmonic generation (E-FISH) method. Here, electric field components were determined by simultaneously conducting E-FISH measurements with the incident laser polarized in two orthogonal directions relative to the PT crystal. Results of this work demonstrate the spatial distribution of electric field around the PT’s output distal end and how it evolves as a function of time. Notably, the strongest electric field appears on the face of the PT’s distal surface, near the top and bottom edges and decreases by approximately 70% over 3 mm. The time delay between the PT’s input voltage and measured electric field indicates that there is an about 0.45 phase difference between the PT’s input voltage and output signal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Femtosecond pulse parameter estimation from photoelectron momenta using machine learning

Abstract Deep learning models have provided huge interpretation power for image-like data. Specifically, convolutional neural networks (CNNs) have demonstrated incredible acuity for tasks such as feature extraction or parameter estimation. Here we test CNNs on strong-field ionization photoelectron spectra, training on theoretical data sets to ‘invert’ experimental data. Pulse characterization is used as a ‘testing ground’, specifically we retrieve the laser intensity, where ‘traditional’ measurements typically lead to 20% uncertainty. We report on crucial data augmentation techniques required to successfully train on theoretical data and return consistent results from experiments, including accounting for detector saturation. The same procedure can be repeated to apply CNNs in a range of scenarios for strong-field ionization. Using a predictive uncertainty estimation, reliable laser intensity uncertainties of a few percent can be extracted, which are consistently lower than those given by traditional techniques. Using interpretability methods can reveal parts of the distribution that are most sensitive to laser intensity, which can be directly associated with holographic interferences. The CNNs employed provide an accurate and convenient ways to extract parameters, and represent a novel interpretational tool for strong-field ionization spectra.

Szołdra, Tomasz (ORCID:0000000228970506)↗
DOE OSTI2020

Pathways to excitation of atoms with bicircular laser pulses

We study the excitation of the hydrogen atom by bichromatic circularly polarized laser pulses using numerical solutions of the time-dependent Schrödinger equation. The results are in agreement with the selection rules for multiphoton processes in such fields, namely, excited states are populated in which orbital angular momentum and magnetic quantum numbers are either both odd or both even, independent of the relative helicity, peak intensity, and pulse duration of the pulses. For co-rotating pulses the results show that excitation predominantly proceeds to states with magnetic quantum number of the same helicity as the laser pulses. Besides pathways via direct photon absorption from the ground state our results indicate that a transfer of population among the Rydberg states occurs via Λ-type transitions. In the case of counter-rotating pulses the largest excitation probability is found for Rydberg states that differ in magnetic quantum number by Δm=±3. This pattern allows us to estimate how many photons from each of the two bichromatic fields have been absorbed. Finally, we confirm that a population in Rydberg states beyond a maximum orbital angular quantum number is unlikely.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Development and experimental study of an automated laser-foil-printing additive manufacturing system

Purpose This paper aims to present the development and experimental study of a fully automated system using a novel laser additive manufacturing technology called laser foil printing (LFP), to fabricate metal parts layer by layer. The mechanical properties of parts fabricated with this novel system are compared with those of comparable methodologies to emphasize the suitability of this process. Design/methodology/approach Test specimens and parts with different geometries were fabricated from 304L stainless steel foil using an automated LFP system. The dimensions of the fabricated parts were measured, and the mechanical properties of the test specimens were characterized in terms of mechanical strength and elongation. Findings The properties of parts fabricated with the automated LFP system were compared with those of parts fabricated with the powder bed fusion additive manufacturing methods. The mechanical strength is higher than those of parts fabricated by the laser powder bed fusion and directed energy deposition technologies. Originality/value To the best knowledge of authors, this is the first time a fully automated LFP system has been developed and the properties of its fabricated parts were compared with other additive manufacturing methods for evaluation.

Engineering↗
DOE OSTI2025

Convergent Manufacturing of Large-Scale Components for Nuclear Applications, via Additive Manufacturing and Powder Metallurgy Hot Isostatic Pressing

Powder metallurgy (PM)–hot isostatic pressing (PM-HIP) has long been recognized as a powerful route for producing fully dense, near net shape metallic components. By consolidating powders under high temperature and pressure, HIP provides isotropic properties, uniform microstructures, and scalability to complex geometries that are vital for sectors such as aerospace, energy, and nuclear power. Yet despite these advantages, the technology has remained constrained by costly trial and error canister fabrication, limitations of conventional forging, and incomplete knowledge about how the canister design influences final part properties. Additive manufacturing (AM), by contrast, thrives on design freedom and geometric flexibility but struggles with speed, scalability, and cost when applied to very large structures. The research presented in this report investigated how a convergent manufacturing approach, combining AM with PM-HIP, can merge the strengths of both technologies, leveraging AM’s flexibility for canister design and HIP’s consolidation capability to deliver reliable, large, and complex parts. The work progressed through three case studies that built on one another in scale and complexity. Small cylindrical canisters fabricated by conventional methods, laser powder bed fusion, and directed energy deposition were filled with stainless steel powders and subjected to HIP. The resulting parts demonstrated near-full density and mechanical properties on par with wrought stainless steel, showing for the first time that AM canisters can be a direct substitute for conventional ones without sacrificing quality. The next step involved a medium-scale, noncentrosymmetric T-valve, which is an enclosed, multibranch geometry that tested the limits of AM + PM-HIP integration. The T-valve achieved predictable shrinkage and uniform densification, confirming feasibility for enclosed designs. However, this study also revealed oxide inclusions and interfacial challenges at the AM + HIP boundary, underscoring the critical importance of controlling interface chemistry and employing robust, in situ strategies, such as melt pool monitoring and thermal monitoring, coupled with nondestructive evaluation techniques such as x-ray computed tomography. Finally, the effort culminated in fabricating a large-scale impeller weighing nearly 2000 lb and spanning 5 ft in diameter. Produced via multirobot wire arc AM and hot isostatic pressed to near-full density, the impeller validated industrial-scale feasibility. Predictive models closely matched experimental shrinkage, tensile properties were spatially uniform across the component, and the AM + PM-HIP interface proved mechanically sound despite the presence of oxide-decorated prior particle boundaries. This large-scale demonstration is a major milestone, showing that hybrid AM + PM‑HIP can reliably deliver components at reactor-relevant scales. Collectively, these studies charted a logical pathway: small-scale work built scientific confidence, medium-scale work highlighted opportunities and challenges, and large-scale work proved industrial impact. The overarching conclusion of this report is that AM + PM-HIP should not be seen as a replacement for forging but as a complementary pathway that provides the US with flexibility, resilience, and new options for manufacturing nuclear-grade components. Looking ahead, several directions emerge as critical to sustaining progress. Predictive modeling must become faster, more accessible, and more accurate, with digital twins and machine learning reducing reliance on trial and error. Powders and alloys must be optimized for HIP, with improved cleanliness, reduced oxides, and tailored chemistries that enhance creep, fatigue, and irradiation resistance. Interfaces between AM and HIP regions must be better engineered through coatings, machining strategies, and surface treatments to mitigate oxide formation and ensure reliable bonding to explore opportunities for HIP of targeted compositional parts, as well as multimaterial HIP cladding applications. Monitoring and nondestructive evaluation need to expand, incorporating multimodal sensors, x-ray computed tomography, and real-time data integration through platforms such as Pelican. At the same time, the pathway to industrial adoption requires techno-economic analysis, machinability studies, and qualification frameworks aligned with industry and regulatory standards. Finally, workforce and academic engagement must be strengthened. Programs that train technicians and engineers for US Navy and US Department of Energy manufacturing challenges should be paired with academic partnerships to support fundamental research, with open sharing of non-export-controlled data to accelerate innovation and build the next generation of experts. In conclusion, this report demonstrates that hybrid AM + PM-HIP is scientifically viable and strategically important. By combining the design agility of AM with the consolidation strength of HIP and embedding modeling, monitoring, and workforce development, this approach provided a transformative new capability for US manufacturing. The path forward is clear: hybrid AM + PM-HIP is not just a promising research direction but is also potentially an industrially relevant pathway that can reshape how nuclear-grade components are designed, qualified, and deployed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Compact absorptivity measurement system for additive manufacturing

The present disclosure relates to a powdered bed fusion additive manufacturing (PBFAM) apparatus. The apparatus uses a container for holding a quantity of powdered material. The container has a bottom wall for supporting the powdered material, wherein the bottom wall is made from the same material as the powdered material. A temperature sensing subsystem is coupled to a portion of the container for detecting a temperature of the container. A laser generates an optical beam directed at the powdered material held by the container for melting the powdered material. A controller receives temperature information from the temperature sensing subsystem and determines an absorptivity of the powdered material based on the temperature information.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unzipping hBN with ultrashort mid-infrared pulses

Manipulating the nanostructure of materials is critical for numerous applications in electronics, magnetics, and photonics. However, conventional methods such as lithography and laser writing require cleanroom facilities or leave residue. We describe an approach to creating atomically sharp line defects in hexagonal boron nitride (hBN) at room temperature by direct optical phonon excitation with a mid-infrared pulsed laser from free space. We term this phenomenon “unzipping” to describe the rapid formation and growth of a crack tens of nanometers wide from a point within the laser-driven region. Formation of these features is attributed to the large atomic displacement and high local bond strain produced by strongly driving the crystal at a natural resonance. This process occurs only via coherent phonon excitation and is highly sensitive to the relative orientation of the crystal axes and the laser polarization. Its cleanliness, directionality, and sharpness enable applications such as polariton cavities, phonon-wave coupling, and in situ flake cleaving.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hot-electron preheat and mitigation in polar-direct-drive experiments at the National Ignition Facility

Target preheat by superthermal electrons from laser–plasma instabilities is a major obstacle to achieving thermonuclear ignition via direct-drive inertial confinement fusion at the National Ignition Facility (NIF). Polar-direct-drive surrogate plastic implosion experiments were performed on the NIF to quantify preheat levels at ignition-relevant scale and develop mitigation strategies. Here, the experiments were used to infer the hot-electron temperature, energy fraction, divergence, and to directly measure the spatial hot-electron energy deposition profile inside the imploding shell. Silicon layers buried in the ablator are shown to mitigate the growth of laser–plasma instabilities and reduce preheat, providing a promising path forward for ignition designs at an on-target intensity of about 10 15 W/cm 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Optimization Methodology of Polar Direct-Drive Illumination for the National Ignition Facility

Improved laser illumination uniformity drives shocks and implosions to create more extreme high energy density environments. Predominantly, the geometry of experiments that can be performed is dictated by the layout of beams at laser facilities, limiting inter-facility and multiscale investigations. This Letter presents the first automated, algorithmic approach for generating illumination configurations for high energy density experiments. The method is demonstrated in comparison to a polar direct drive solid target experiment at the National Ignition Facility. The new illumination configuration is simulated to create greater than ×3 higher peak pressure and almost ×2 higher density by maintaining better shock uniformity. Furthermore, the optimization process is performed with reduced computational expense and isotropic plasma profiles, while accounting for the impact of cross-beam energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structural basis of transcription: RNA polymerase II substrate binding and metal coordination using a free-electron laser

Catalysis and translocation of multisubunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near-atomic resolution and precise arrangement of key active site components have been elusive. Here, we present the free-electron laser (FEL) structures of a matched ATP-bound Pol II and the hyperactive Rpb1 T834P bridge helix (BH) mutant at the highest resolution to date. The radiation-damage-free FEL structures reveal the full active site interaction network, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg 2+ , and, more importantly, a putative third (site C) Mg 2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structures indicate that the third Mg 2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/BH interactions induce conformational changes that could allow translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the T834P mutant reveal rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Science & Technology - Other Topics↗
DOE OSTI2023

The importance of laser wavelength for driving inertial confinement fusion targets. II. Target design

We describe details of radiation-hydrodynamics simulations of directly driven targets for inertial confinement fusion using laser drivers with different laser wavelengths. Of particular interest here are comparisons of frequency-tripled glass (laser wavelength 351 nm) lasers with the argon fluoride (193 nm) and krypton fluoride (248 nm) excimer lasers and the effects that these laser wavelengths have on the target designs. We explore the effect these drivers have on the compromise involved between lowering laser plasma instabilities (LPIs) or hydrodynamic instabilities while providing high gains and seek to quantify this trade-off. Short-wavelength drivers have significant advantages, primarily in using less power and energy to drive targets. Additionally, they expand the allowed operating regime that is constrained by LPI avoidance and the production of higher pressures needed for more hydrodynamically stable targets. Potential disadvantages to shorter drive wavelengths, such as increased imprint, are examined and found to be unimportant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Spinning Hall Probe magnetic compass

In a large range of physics and space experiments, the direction of the magnetic field needs to be determined with an accuracy on the level of one milli radian or better. We have proposed a new type of magnetic compass - a rotation-based one, which provides an alternating signal from the Hall probe proportional to the value of the magnetic field transverse to the axis of rotation. The spinning Hall probe device has been realized. Alignment of the spinning axis for a minimum (zero) value of the signal allows us to find the direction of the magnetic field. The device does not require calibration and is free of any drift problems. The measurement of the rotation axis direction was accomplished by means of a laser and a flat mirror attached to the rotor. The constructed prototype achieved an accuracy for the magnetic field direction in the experiment with the polarized He-3 target on the level of one milli radian.

Wojtsekhowski, Bogdan↗
DOE OSTI2020

Direct observation of pore formation mechanisms during LPBF additive manufacturing process and high energy density laser welding

Laser powder bed fusion (LPBF) is a 3D printing technology that can print parts with complex geometries that are unachievable by conventional manufacturing technologies. However, pores formed during the printing process impair the mechanical performance of the printed parts, severely hindering their widespread application. Here, we report six pore formation mechanisms that were observed during the LPBF process. Our results in this study reconfirm three pore formation mechanisms - keyhole induced pores, pore formation from feedstock powder and pore formation along the melting boundary during laser melting from vaporization of a volatile substance or an expansion of a tiny trapped gas. We also observe three new pore formation mechanisms: (1) pore trapped by surface fluctuation, (2) pore formation due to depression zone fluctuation when the depression zone is shallow and (3) pore formation from a crack. The results presented here provide direct evidence and insight into pore formation mechanisms during the LPBF process, which may guide the development of pore elimination/mitigation approaches. Since certain laser processing conditions studied here are similar to the situations in high energy density laser welding, the results presented here also have implications for laser welding.

42 ENGINEERING↗
DOE OSTI2025

Depth-Resolved Lithiated Gradients in Pristine and Laser-Ablated Anodes During Fast Charging

Laser ablating 3D electrode microstructures is a technique to improve Li-ion battery fast-charge performance. This technique has been theoretically proposed and electrochemically validated previously in the literature. The fundamental principle underlying laser ablation is that the ablated features reduce Li-ion transport pathways, improving access to the electrode active material near the current collector. This, in turn, promotes more homogeneous electrode utilization. The present study seeks to directly affirm the physics attributed to laser ablation using operando high-speed synchrotron X-ray diffraction. In this study, depth-resolved graphite lithiation gradients are measured operando during high-rate (15 min) charging. The depth-resolved lithiation dynamics of both ablated and non-ablated anodes are compared. The results highlight that the laser-ablated graphite electrode has notably more homogeneous utilization as compared to the non-ablated electrode. Additionally, the ablated electrode has a significant delay in reaching the maximum graphite lithiation at the separator, indicating less propensity for lithium plating. During low rate delithiation/discharge (2 hr), the two cells' lithiation gradients converge. Notably, a calibrated physics-based electrochemical model accurately reflects experimental findings, suggesting the potential to use pseudo-4D models not only to optimize laser ablation parameters in fast-charge capable electrodes but also to guide fast charging protocols that avoid lithium plating.

25 ENERGY STORAGE↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗