SEARCH · Engineering Papers

Results for “Shape Memory Behavior”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

464 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING↗

NASA NTRS2026

Development, Characterization, and Validation of Elevated-Temperature Constitutive Models: Deformation and Damage

This report provides a brief review of experimentally observed hereditary and nonhereditary material behavior along with background information on standard as well as advanced internal state variable constitutive modeling at elevated temperature. A description of exploratory, characterization, and validation testing is presented along with a detailed outline of what constitutes “sufficient” data content (i.e., quality and quantity) for developing or enhancing, characterizing, and validating a particular sophisticated nonlinear time- and history-dependent (hereditary) class of constitutive models known as GVIPS (generalized viscoplasticity with potential structure). The tests described are necessary to reveal a material’s behavior in the reversible (or viscoelastic) and irreversible (or viscoplastic) regimes, for the identification of both deformation and damage model parameters. Results presented are primarily for metallic materials. In all cases, both uniaxial and multiaxial tests are described, and the linkage between the specific tests and the parameters within the model that can be characterized from the results of these tests are also discussed. Discussion is also provided relative to the role information management must play relative to material data collection, analysis, maintenance, and dissemination. The need for such an information system is particularly important as both analyst and designer move toward utilizations of sophisticated, nonlinear time- and history dependent (hereditary) constitutive models. Lastly, the concept of state space and its utility in understanding and establishing constitutive models is addressed throughout. The intent behind this document is to help both the modeler and experimentalist understand each other’s specific points of view and provide guidance for both model development and characterization. Emphasis has been placed on providing the mechanician with information regarding how tests are performed and what issues to be aware of when interpreting results. It is hoped that experimentalists will take away a new perspective on the types of information that modelers and analysts are looking for from them.

Constitutive Modeling↗

NASA NTRS1982

A Survey of Methods for Measuring Thermally Induced Distortions of Test Articles Undergoing Solar Thermal Vacuum Test

Various methods of measuring the small thermally induced distortions experienced either by various points on the space vehicles or by deviations of surfaces from a known shape during solar thermal vacuum tests are examined. State-of-the-art application of both photographic and real time observation are discussed. The relative merits of each of the methods are compared and evaluated in their applications to different types of test articles and situations. Magnitudes of thermally induced distortions which may be expected to be routinely measurable by the various methods are presented and compared.

R D Rempt↗

DOE OSTI2025

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗

NASA NTRSDate not supplied

3D Woven Mid-Density Carbon Phenolic (3MDCP) Full-Scale Heatshield Development for Mars Sample Return (MSR) Earth Entry System (EES)

The Mars Sample Return (MSR) mission will be returning samples of Martian soil to Earth. The samples will be returned in the Earth Entry System (EES). The EES capsule is protected by a new thermal protection system; the 3-Dimensionally Woven, Mid-Density, Carbon Phenolic (3MDCP) material which is derived from the Heatshield for Extreme Entry Environment Technology (HEEET). The baseline 3MDCP design is a single piece thermal protection system that avoids the manufacturing and certification challenges associated with a tiled configuration. A 3MDCP heatshield begins as a flat woven preform that is formed to the final heatshield shape and then infused with phenolic resin. A NASA Ames lead team, in collaboration with Spirit Textiles (formerly TEAM Inc.), Fiber Materials Inc. (FMI, a Spirit AeroSystems Company), and Kratos SRE, has been developing and demonstrating the processes to fabricate the 3MDCP heatshield and characterize the resulting properties. The process of forming the flat woven preform into the final heatshield shape, a sphere-cone geometry with a 52.5° cone angle, involves local but substantial movement of the yarns in the weave. This results in continuous fibers across the single piece heatshield, albeit with property variations between different regions of the heatshield. This presentation will provide a high-level status of 3MDCP development for the MSR EES heatshield. This will include an overview of the manufacturing processes with an emphasis on the impact of forming on fiber orientation, material properties, and performance. Results comparing preliminary formed and flat materials will be shown. Additionally, the presentation will layout the broader plan for property testing and tie-in to aerothermal performance.

Peter Edward Marshall↗

NASA NTRS2025

Sulfur outgassing and in-gassing in lunar orange glass beads and implications for 33S “Anomaly” in the Moon

Our recent investigations have discovered inward diffusion (in-gassing) of moderately volatile elements (MVEs; e.g., Na, K and Cu) from volcanic gas into volcanic beads/droplets. In this work, we examine the distribution of sulfur in lunar orange glass beads. Our analyses reveal that sulfur exhibits a non-uniform distribution across the beads, forming “U” or “W” shaped profiles typical of in-gassing. A model developed to assess sulfur contributions from different sources (original magmatic sulfur versus atmospheric in-gassed sulfur) in the orange beads indicates that atmospheric sulfur in-gassed during eruption contributes approximately 9–24 % to the total sulfur content of an orange bead, averaging around 16 %. This in-gassed sulfur is derived from the eruption plume, where atmospheric sulfur could undergo photochemical reactions induced by UV light, leading to mass independent fractionation and a distinct sulfur isotope signature. Interestingly, a recent study discovered a small mass independent isotope fractionation of sulfur in lunar orange glass beads in drive tube 74002/1 and a lack of such mass independent isotope fractionation in black glass beads in the same lunar sample. This finding contrasts with sulfur in lunar basalts, which typically exhibit mass dependent fractionation. With our work, the observed mass independent fractionation signal in sulfur isotopes of orange beads can be attributed to the in-gassing of photolytic sulfur in the optically thin part of the eruption plume where UV light can penetrate. Using the sulfur isotope data of lunar orange beads, we estimate that the Δ33S value of atmospheric sulfur is approximately −0.18 ‰. Our study provides new insights into the complex dynamics of volatile elements in lunar volcanic processes, highlighting the role of in-gassing in shaping sulfur isotope signatures in volcanic glass beads.

Xue Su↗

DOE OSTI2026

Process control-enabled mitigation of microstructural and plastic heterogeneities in additively manufactured Grade 91 steel

Synergizing wire arc-directed energy deposition (WA-DED) additive manufacturing (AM) with particle-strengthened creep strength-enhanced ferritic (CSEF) steels enables fabrication and repair of critical power-plant components. Investigations focused on fusion-welded particle-strengthened CSEF steels, such as Grade 91 steel, have linked microstructurally heterogeneous regions—forming due to heat affected zones (HAZ)—with premature failure during elevated temperature service. Fusion-based AM, including WA-DED, likewise generates microstructurally and plastically heterogeneous regions due to spatiotemporally varying thermokinetics during deposition. However, works investigating such microstructural heterogeneities, their implications for mechanical behavior, and strategies to mitigate their formation remain scarce. This work identifies microstructurally and plastically heterogeneous regions within the WA-DED-processed Grade 91 steel. Spatial microhardness variations in the as-fabricated specimen correlate with the variation in the attributes of grain, martensitic microstructure, and precipitates across the fusion zone and HAZ. Digital image correlation-enabled tensile tests performed at 500 °C revealed pronounced deformation localization and a wave-like strain distribution, with wavelength close to the melt pool depth, indicating susceptibility of the as-fabricated components to premature creep failure. Such heterogeneity in microstructural and mechanical behavior was attributed to recurring solid-state phase transformations. Subsequently, an interlayer temperature control strategy was implemented, wherein maintaining interlayer temperature above the martensitic start temperature mitigated the heterogeneous microstructural and plastic response in the as-fabricated condition. Findings open pathways to achieving deformation-localization- and creep-resistant microstructures in WA-DED fabricated particle-strengthened CSEF steel components, reducing reliance on post-welding heat treatments—conventionally required to enhance creep resistance—and enabling on-demand, short lead-time fabrication of next-generation power-plant components.

Heat affected zones↗

NASA NTRS1998

High performance zinc anode for battery applications

An improved zinc anode for use in a high density rechargeable alkaline battery is disclosed. A process for making the zinc electrode comprises electrolytic loading of the zinc active material from a slightly acidic zinc nitrate solution into a substrate of nickel, copper or silver. The substrate comprises a sintered plaque having very fine pores, a high surface area, and 80-85 percent total initial porosity. The residual porosity after zinc loading is approximately 25-30%. The electrode of the present invention exhibits reduced zinc mobility, shape change and distortion, and demonstrates reduced dendrite buildup cycling of the battery. The disclosed battery is useful for applications requiring high energy density and multiple charge capability.

Casey, John E., Jr.↗

NASA NTRS2005

Chapter 19B - Antenna Radiation Patterns

Modern civil and military aircraft are equipped with a variety of communication devices, radio navigation equipment, and air traffic control systems. For all of these devices appropriate antennas must be available to transmit and receive the signals. As a result, as many as 30 antennas, and sometimes even more, are mounted around today's aircraft. For this reason, it is necessary to know the capabilities of the receiving/transmitting equipment on board. These capabilities are driven by the antenna characteristics. Therefore, coverage, shading, and beam pointing pattern data are necessary for optimum electronic coverage. The aircraft antenna has to convert the available power density of the electromagnetic field to an electric voltage at its connector. This voltage level must be sufficiently high to operate the connected equipment. The reciprocity principle states that it makes no difference if the antenna is receiving or transmitting. For the certification of the aircraft antennas it must be proven that the antennas are at least generating the minimum receiver input voltage which is required for each radio service. The prime condition is, of course, that the specified field power densities of the different radio navigation and communication services are available. An important property of a radio frequency link is the electro-magnetic field intensity at every point in space for a given output power of the antenna. As the propagation of radio frequency waves in free space is well known, the spatial distribution of the field intensity needs only be measured at one spatial sphere around the antenna. The information is usually given as distributions along the circumference of flat sections through this sphere: each of these is called an ARP. Several ARPs are usually required to describe the complete spatial antenna pattern of an antenna. In addition to mathematical modeling, measurements on sub-scale models, and static measurements on full size models or aircraft on the ground, dynamic measurements on aircraft in flight play the most important role in the aircraft antenna testing and certification. [19B-1] The shape of an ARP, for one given frequency, is determined by the shape of the antenna and the shape and material of the surface it is mounted on. As the directly transmitted waves interfere with waves reflected by the aircraft skin with its complex geometry, and the surface material parameters are only roughly known, it is not possible to predict the ARP with the required accuracy. In ground measurements the earth's surface also acts as a reflector thereby causing the ground ARP to be different from the in-flight ARP. As the in-flight ARP is the ARP we are actually interested in it becomes clear that it is necessary to conduct in-flight ARP measurements. The Flight Test Engineer must be aware of the needs of the specialists who are establishing the test program for measuring antenna patterns and radar cross sections. These tests will require special test equipment and dedicated flights to obtain the data that they require. This Section provides an introduction to the principles of determining antenna patterns and the flight techniques for determining both antenna patterns and radar cross section. Reference 19B-1 provides detailed information.

Helmut Bothe↗

NASA NTRS1994

Control of an Axisymmetric Jet Using Vortex Generators

The results of an experimental investigation on the effect of vortex generators, in the form of small tabs at the nozzle exit on the evolution of a jet, are reported in this paper. Primarily tabs of triangular shape are considered, and the effect is studied up to a jet Mach number of 1.8. Each tab is found to produce a dominant pair of counter-rotating streamwise vortices having a sense of rotation opposite to that expected from the wrapping of the boundary layer. This results in an inward indentation of the mixing layer into the core of the jet. A triangular-shaped tab with its apex leaning downstream, referred to as a delta tab, is found to be the most effective in producing such vortices, with a consequential large influence on the overall jet evolution. Two delta tabs, spaced 180 deg apart, completely bifurcate the jet. Four delta tabs stretch the mixing layer into four 'fingers,' resulting in a significant increase in the jet mixing downstream. For six delta tabs the mixing layer distortion settles back to a three finger configuration through an interaction of the streamwise vortices. The tabs are found to be equally effective in jets with turbulent or laminar initial boundary layers. Two sources of streamwise vorticity are postulated for the flow under consideration. One is the upstream 'pressure hill,' generated by the tab, which constitutes the main contributor of vorticity to the dominant pair. Another is due to vortex filaments shed from the sides of the tab and reoriented downstream by the mean shear of the mixing layer. Depending on the orientation of the tab, the latter source can produce a vortex pair having a sense of rotation opposite to that of the dominant pair. In the case of the delta tab, vorticity from the two sources add, explaining the strong effect in that configuration.

Streamwise Vorticity↗

DOE OSTI2024

Adaptive Scalpel Scanning Probe Microscopy for Enhanced Volumetric Sensing in Tomographic Analysis

Controlling nanoscale tip‐induced material removal is crucial for achieving atomic‐level precision in tomographic sensing with atomic force microscopy (AFM). While advances have enabled volumetric probing of conductive features with nanometer accuracy in solid‐state devices, materials, and photovoltaics, limitations in spatial resolution and volumetric sensitivity persist. This work identifies and addresses in‐plane and vertical tip‐sample junction leakage as sources of parasitic contrast in tomographic AFM, hindering real‐space 3D reconstructions. Novel strategies are proposed to overcome these limitations. First, the contrast mechanisms analyzing nanosized conductive features are explored when confining current collection purely to in‐plane transport, thus allowing reconstruction with a reduction in the overestimation of the lateral dimensions. Furthermore, an adaptive tip‐sample biasing scheme is demonstrated for the mitigation of a class of artefacts induced by the high electric field inside the thin oxide when volumetrically reduced. This significantly enhances vertical sensitivity by approaching the intrinsic limits set by quantum tunneling processes, allowing detailed depth analysis in thin dielectrics. The effectiveness of these methods is showcased in tomographic reconstructions of conductive filaments in valence change memory, highlighting the potential for application in nanoelectronics devices and bulk materials and unlocking new limits for tomographic AFM.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Passive Confirmation of the Presence Of High-Explosive Material Via Neutron Transmission Spectroscopy

The goal of this project was to explore a novel approach for identifying presence and potential types of high explosives (HE) in treaty-controlled items with passive neutron sources by using passive neutron transmission spectroscopy. We predominantly look to measure relative elemental abundance of Carbon, Hydrogen, Oxygen, and Nitrogen (CHON) since most relevant materials are certain mix of these elements, as shown in Table 1. Previously, most studies exploring potential identification of CHON elemental content of targets in the vicinity of passive neutron source were focused solely on gamma spectroscopy using high resolution gamma detectors. Using neutron transmission spectroscopy with pulse-shape discrimination (PSD) capable organic scintillators is not only a novel idea in of itself, but arguably necessary for accurate identification of the CHON elemental content of an interrogated target. While high resolution gamma spectroscopy in principle can determine a presence of hydrogen and nitrogen, it is effectively blind to a quantitatively measuring areal densities of C, O. The initially proposed approach, described in detail in the project proposal, was to leverage previous Monte Carlo study on neutron transmission spectroscopy using slab shaped target interrogated with well-collimated (“pencil”) neutron beam. In this study, a simulated test CO 2 target was interrogated by a pencil neutron beam and transmitted neutron spectra were measured using PSD capable organic scintillator using MLEM based unfolding technique. To establish relative elemental content of carbon and oxygen, the unfolded neutron spectrum was fit with a parametrized combination of their respective elemental spectral templates. The individual spectral template was calculated by convoluting ENDF neutron crosssection with the known resolution of the PSD detector used in the study. This approach in the studied configuration successfully quantitatively established relative elemental composition of carbon and oxygen.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Development of large-scale, 3D Printed high temperature ceramic material

Through this collaborative effort, a new 3D printing platform for ceramics called laser-induced slip casting (LIS) was explored and developed for improving the processing and manufacture of silicon carbide (SiC), a high temperature ceramic material. This method prints layers of ceramic slips/slurries with subsequent selective-laser heating to dry each layer to build a 3D structure. The completed and dried part is then sintered. This manufacturing process is inherently lower cost than alternative ceramic printing methods such as binder jet technology or stereolithography when it comes to the feedstock material but is more expensive than robocasting or direct ink writing. However, it has the potential to make more controlled parts with less defects compared to robocasting. The largest cost is the heating source for the printer. With this technique, there is potential to make large ceramic parts in a near-net shape. Further, there is no known commercial manufacturing of 3D printed ceramics that creates a large range of different high-density ceramics at large scale. The goals and outcomes for the development of the new printing technology were: 1) assessing the technology with alumina by characterizing coupons, 2) printing and sintering of silicon carbide (SiC), and 3) characterization and properties testing of materials printed and a scale-up of SiC part(s). Through this collaborative project, it was sought to provide the best solution to achieving highly dense and near-net shaped ceramic parts at large scale and reduced cost. The goal was to produce high density sintered materials for applications in defense, energy generating systems, armor, and wear parts using 3D printing methods. Variables including powder particle sizes, dispersant molecular weights (MWs), binders, printing parameters, and post processing were explored as well as the characterization of the physical properties (add specifics here).

36 MATERIALS SCIENCE↗

DOE OSTI2026

Interfacial magnetization enhancement and the perpendicularly magnetized component in Fe16N2 thin films

A detailed study of depth-resolved magnetization and microstructure of Fe16N2 thin films on MgO (001) and MgAl2O4 (001) substrates and Fe and Cr seed layers is presented. Two aspects of the magnetic properties of Fe16N2 thin films are discussed. First, magnetization enhancement at the interface is observed. Strain and nitrogen deficiency are discussed as possible interfacial mechanisms contributing to this enhancement. Second, the perpendicularly magnetized component (PMC) is identified in Fe16N2 thin films. Correlation with microstructural observations suggests that the PMC is associated with V-shaped grains that are not fully confined within the continuous Fe16N2 layer.

Hang, Xudong [University of Minnesota]↗

NASA NTRS1997

Radiation Temperature and Extinction of Transient Gaseous Diffusion Flames in Microgravity

The absence of buoyancy-induced flows in micro-g and the resulting increase in the reactant residence time significantly alters the fundamentals of many combustion processes. Substantial differences between 1-g and micro-g flames have been reported in experiments on candle flames, flame spread over solids, droplet combustion and others. These differences are more basic than just in the visible flame shape. Longer residence times and higher concentration of combustion products in the flame zone create a thermochemical environment which changes the flame chemistry and the heat and mass transfer processes. Processes such as flame radiation (and its interaction with flame chemistry), that are often ignored under normal gravity, become very important and sometimes even controlling. This is particularly true for conditions at extinction of a micro-g diffusion flame.

Arvind Atreya↗

NASA NTRS2026

Hydrology Copilot: A Cloud-Native Ai System for Hydrological Data Analysis

The emergence of AI-driven Earth observation systems promises to broaden access to petabyte-scale geospatial data beyond domain specialists. However, translating this vision into operational scientific infrastructure requires addressing fundamental challenges in data virtualization, code transparency, and domain-specific reasoning. We present Hydrology Copilot, a cloud-native AI framework for natural-language-driven analysis of Earth observation data. To demonstrate operational capabilities at scale, we implement the system using NASA's North American Land Data Assimilation System version 3 (NLDAS-3), which provides surface meteorological forcing and land-surface model output across North and Central America at 1-km resolution, from which drought diagnostics are derived. The system integrates five core contributions: (1) scalable data virtualization using Kerchunk-based cloud optimized access, achieving a 1.5 to 4.6 times improvement in I/O latency across benchmark queries spanning regional single-day extractions (4.6 times speedup) to continental monthly aggregations (1.5 times speedup); (2) transparent code generation through Microsoft Azure AI Foundry agents that expose executable Python workflows for scientific verification; (3) persistent conversational memory enabling multi-turn analytical discourse across sessions; (4) intelligent query validation that enforces dataset boundaries and resolves ambiguous requests before execution; and (5) a multi-agent architecture coordinating query parsing, code generation, and visualization. We evaluate the system through drought-monitoring workflows, demonstrating reliable code generation, accurate results validated against reference computations and the operational U.S. Drought Monitor, and efficient operation across increasingly complex tasks. By bridging natural-language interfaces with rigorous hydrological analysis, Hydrology Copilot advances beyond proof-of-concept demonstrations to provide a deployable framework for operational Earth science applications.

Data virtualization↗

DOE OSTI2024

Molecules to Masterpieces: Bridging Materials Science and the Arts

Art and materials innovation have always been intertwined, dating back to the earliest human creations. In modern times, however, the increasing specialization of materials science often restricts artists' access to cutting-edge materials. Here, the materials science aspects of an art-science collaboration between artist Kimsooja and the Wiesner Lab at Cornell University, are detailed. The project involves the development of a custom-made iridescent block copolymer coating by means of self-assembly, originally applied to transparent window panels of a façade for the ≈14 m tall art installation: A Needle Woman: Galaxy Is a Memory, Earth is a Souvenir by artist Kimsooja. After several exhibitions in the US and Europe, the installation is now part of the permanent museum collection at Yorkshire Sculpture Park in Wakefield, UK. Full characterization of the solution blade-cast coatings show shear aligned, standing up lamellar morphologies that behave as volume-phase gratings with periodicities between 300 and 400 nm. Coatings are also applied to foldable (origami) paper and converted into iridescent porous ceramic materials. Furthermore, it is hoped this work inspires and informs communities across materials science, the arts, and architecture.

Architectural nanomaterials↗

DOE OSTI2025

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗