Search NASA⌕ Search

SEARCH · Search NASA

Results for “Volume fractions”

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 109 records · Page 6

Point-particle drag, lift, and torque closure models using machine learning: Hierarchical approach and interpretability

Developing deterministic neighborhood-informed point-particle closure models using machine learning has garnered interest recently from the dispersed multiphase flow community. The robustness of neural models for this complex multibody problem is hindered by the availability of particle-resolved data. Here, the present work addresses this unavoidable limitation of data paucity by implementing two strategies: (1) by using a rotation and reflection equivariant neural network and (2) by pursuing a physics-based hierarchical machine learning approach. The resulting machine-learned models are observed to achieve a maximum accuracy of 85% and 96% in the prediction of neighbor-induced force and torque fluctuations, respectively, for a wide range of Reynolds number and volume fraction conditions considered. Furthermore, we pursue force and torque network architectures that provide universal prediction spanning a wide range of Reynolds number (0.25 ≤ Re ≤250) and particle volume fraction (0 ≤ φ ≤0.4). The hierarchical nature of the approach enables improved prediction of quantities such as streamwise torque, by going beyond binary interactions to include trinary interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of ammonia addition on suppressing soot formation in methane co-flow diffusion flames

Due to issues surrounding carbon dioxide emissions from carbon-containing fuels, there is growing interest in ammonia (NH 3 ) as an alternative combustion fuel. One attractive method of burning NH 3 is to co-fire it with hydrocarbons, such as natural gas, and in this case soot formation is possible. To begin understanding the influence of NH 3 on soot formation when co-fired with hydrocarbons, soot volume fractions and mole fractions of gas-phase species were computationally and experimentally interrogated for CH 4 flames with up to 40% NH 3 by volumetric fuel fraction. Mole fractions of gas-phase species, including C 2 H 2 and C 6 H 6 , were measured with on-line electron impact mass spectrometry, and soot volume fractions were obtained via color-ratio pyrometry. The simulations employed a detailed chemical mechanism developed for capturing nitrogen interactions with hydrocarbons during combustion. The results are compared to findings in N 2 -CH 4 flames, in order to separate thermal and dilution effects from the chemical influence of NH 3 on soot formation. Experimentally, C 2 H 2 concentrations were found to decrease slightly for the NH 3 -CH 4 flames relative to N 2 -CH 4 flames, and a stronger suppression of C 6 H 6 was found for NH 3 relative to N 2 additions. The measured results show a strong suppression of soot with the addition of NH 3 , with soot concentrations reduced by over a factor of 10 with addition of up to 20% or more NH 3 by mole fraction. The model satisfactorily captured relative differences in maximum centerline C 2 H 2 , C 6 H 6 , and soot concentrations with addition of N 2 , but was unable to match measured differences in NH 3 -CH 4 flames. Furthermore, these results highlight the need for an improved understanding of fuel-nitrogen interactions with higher hydrocarbons to enable accurate models for predicting particulate emissions from NH 3 /hydrocarbon combustion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Particle‐Size‐Dependent Lithium‐Ion Transport in PEO/LLZO Composite Electrolytes

Lithium-metal batteries with solid electrolytes can deliver higher energy density and improved safety than conventional Li-ion batteries. Among solid electrolyte candidates, polymer/ceramic composite electrolytes are attractive because they combine polymer flexibility with the high ionic conductivity of ceramics. However, whether ceramic fillers synergistically reduce polarization losses in the polymer matrix remains unclear. A central unknown is the critical polymer/ceramic interfacial resistance (Rint,crit), below which adding ceramics lowers electrolyte overpotential. Here, we present the first macroscale model framework to quantify R int,crit for composite electrolytes based on polyethylene oxide (PEO) and Ta-doped Li 7 La 3 Zr 2 O 12 (LLZO). A 1D model for DC-polarization of tri-layer cells (PEO-LiTFSI/LLZO/PEO-LiTFSI) shows that LLZO surface functionalization reduces the PEO/LLZO interfacial resistance, consistent with electrochemical impedance measurements. Extending to a 2D composite model, we show notably that Rint,crit scales linearly with LLZO particle diameter and shifts toward experimentally accessible values (e.g., 28.8 Ωcm 2 ) as particle size increases. At fixed ceramic volume fraction, larger LLZO particles reduce the number of interfacial crossings, driving more current through the ceramic phase and lowering concentration polarization. In contrast, R int,crit is largely independent of ceramic volume fraction. These results demonstrate that ceramic filler-size engineering can enable synergistic, energy-efficient transport in polymer/ceramic composite electrolytes.

25 ENERGY STORAGE↗

The Application of Differential Scanning Calorimetry to Investigate Precipitation Behavior in Nickel-Base Superalloys Under Continuous Cooling and Heating Conditions

In this study, a suite of experimental tools and fast-acting, numerical-simulation techniques was used to quantify the precipitation behavior of three nickel-base superalloys: IN-100, LSHR, and 718. Experimental methods comprised differential scanning calorimetry (DSC) to establish the specific heat as a function of temperature and selected direct-resistance heating trials (using a Gleeble ® machine) to obtain samples for microstructural analysis. For the DSC experiments, each alloy was cooled at a prescribed constant rate (between 5 and 20 K/min) after an initial soak/equilibration in the high-temperature, single-phase (supersolvus) temperature regime. On-heating DSC trials beginning at ambient temperature were also performed on alloy 718 in three different starting conditions: super- δ -solvus solution treated and water quenched (denoted as ST), solution treated and aged (STA), and solution treated and overaged (STOA). DSC results, revealing the thermal signatures associated with the kinetics of precipitation of γ ' (IN-100, LSHR) or γ ' and γ " (718), were interpreted using a previously-developed fast-acting routine that treats concurrent nucleation, growth, coarsening, and dissolution. For these simulations, special attention was paid to various thermo-kinetic input parameters including equilibrium solvus-approach curves, bulk free energies of transformation, matrix-precipitate interface energies, and effective diffusivities. For the γ - γ ' superalloys (IN-100 and LSHR), estimates of precipitate volume fraction as a function of temperature from the specific-heat data revealed semi-quantitative agreement with simulation predictions. For the γ - γ '- γ " superalloy (718), simulation predictions of precipitate volume fractions were converted to specific heat as a function of temperature and showed semi-quantitative agreement with the direct measurements.

36 MATERIALS SCIENCE↗

High performance Nb/TiNi nanocomposites produced by packaged accumulative roll bonding

Taking full advantage of intrinsic high strength of nano-reinforcements has been proven difficult in a composite until the concept of strain matching was experimental verified in a nanowire Nb/TiNi in-situ composite, which was obtained via solidifying eutectic Nb/TiNi and subsequent severe wire drawing. However, the volume fraction of the nano-reinforcement was dictated by eutectic reaction, which severely limited the development of high-performance composites. Here, the martensitic transformation nanocomposites were successfully acquired with increased volume fraction and changeable morphology of the nano-reinforcement through hot packaged ARB (PARB) and subsequent wire drawing. The nano-reinforcement with lamellar configuration in the composite could exhibit a large lattice strain of similar to 3.8%, which demonstrated the applicability and effectiveness of the strain matching principle in nano-lamellar composite. The obtained composites show a yield strength over 2.0 GPa and a quasi-linear elastic strain as larger as similar to 6.5% with adjustable apparent elastic modulus from 30 to 50 GPa.

42 ENGINEERING↗

Tunable anomalous Hall conductivity through volume-wise magnetic competition in a topological kagome magnet

Magnetic topological phases of quantum matter are an emerging frontier in physics and material science. Along these lines, several kagome magnets have appeared as the most promising platforms. Here, we explore magnetic correlations in the kagome magnet Co 3 Sn 2 S 2 . Using muon spin-rotation, we present evidence for competing magnetic orders in the kagome lattice of this compound. Our results show that while the sample exhibits an out-of-plane ferromagnetic ground state, an in-plane antiferromagnetic state appears at temperatures above 90 K, eventually attaining a volume fraction of 80% around 170 K, before reaching a non-magnetic state. Strikingly, the reduction of the anomalous Hall conductivity (AHC) above 90 K linearly follows the disappearance of the volume fraction of the ferromagnetic state. We further show that the competition of these magnetic phases is tunable through applying either an external magnetic field or hydrostatic pressure. Our results taken together suggest the thermal and quantum tuning of Berry curvature induced AHC via external tuning of magnetic order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Additively Manufactured, Lightweight, Low-Cost Composite Vessels for Compressed Natural Gas Fuel Storage

This project will develop a process to combine AM via direct ink writing (DIW) technology for CFC printing and to use design optimization tools pioneered at LLNL, with advances in resin/composite formulation enabled by chemical and nano-material modification to produce lightweight low-cost CNG tanks. Our approach will yield sub-scale prototype composite pressure tanks equivalent to Type-5 CNG vessel designs that demonstrate a potential cost-benefit advantage. Central to our vision is using agile AM and design based on computationally informed DIW of both short and continuous CF, further coupled with high-performance thermoset polymer matrixes modified by emergent nanomaterials. Our single-stage, multi-material AM technology, combined with a decreased volume fraction of CF and an increased proportion of economically advantaged short fiber, all together drive the reduction in manufacturing time and overall cost. Importantly, reductions in continuous fiber and overall fiber volume fraction will be achieved without detriment to the mechanical strength of the composite vessel. This will be achieved by employing a single process using multi materials grading involving a thermoset resin “ink” modified with aligned nanoplatelets to leverage the efficient tortuous-path gas barrier effect, printed as an inner flexible gas barrier as the initial stage in our manufacturing process before compositionally grading the AM feedstock in real-time to transition to a rigid, structural CF-filled resin. The proposed hybrid construction is projected to achieve pressure ratings at a service range of 2,900–3,600 psi with a 3× burst safety factor comparable to conventional filament-wound composite tanks with an estimated 30–50% reduction in total manufacturing cost.

03 NATURAL GAS↗

Phase Identification in Synchrotron X-ray Diffraction Patterns of Ti–6Al–4V Using Computer Vision and Deep Learning

X-ray diffraction patterns contain information about the atomistic structure and microstructure (defect population) of materials, extracting detailed information from diffraction patterns is complex, demanding and relies on prior knowledge. Here, we hypothesize that deep-learning techniques can help to perform an effective and accurate analysis with high throughput rates. To demonstrate this concept, we applied a novel deep learning framework to determine the evolution of the β-phase volume fraction in a Ti–6Al–4V alloy during heat-treatment from video sequences of 2D diffraction patterns recorded in transmission and with highly monochromatic radiation in a synchrotron beamline. In particular, we studied the impact of network design on prediction reliability and computational performance. Networks of different architectures were trained using 3008 experimental 2D patterns. A well-tuned model was found to reproduce the phase fractions of another experimental data set, consisting of 1100 diffraction patterns, with a mean-square error as small as 2.6 x 10 -4 . The average prediction error of β-phase volume fraction was within 1.6 x 10 -2 (in each diffraction pattern) of the values obtained by conventional methods. Our work demonstrates that convolutional neural networks can evaluate high energy X-ray diffraction patterns with a remarkable level of reliability. Furthermore, it demonstrates the significance of network design on the reliability of predictions and computational performance. The most complex models do not necessarily result in highest accuracy and may even fail to learn from the data.

36 MATERIALS SCIENCE↗

Design Principles to Govern Electrode Fabrication for the Lithium Trivanadate Cathode

A full depth of discharge mathematical model for the lithium trivanadate cathode, considering lithiation of the layered α-phase, phase change, and lithiation of the rock-salt like β-phase at lower potentials, is developed. The coupled electrode-scale and crystal-scale model is fit to electrochemical data, and additionally validated with operando EDXRD. There is good agreement between the simulated and measured spatial variation of the volume fraction of the β-phase. This mathematical model is used to guide electrode fabrication, accounting for both ionic and electronic transport effects. Values of three design parameters—electrode thickness, porosity, and volume fraction of conductor—are identified, and the sensitivity of the energy density to these design parameters is quantified. The model is also used to investigate electrode design to create electrodes that deliver the maximum achievable energy density under the constraint that the α to β-phase transition is avoided, since phase change has been demonstrated to reduce cycle life. The energy density sacrificed to avoid phase change decreases at higher discharge rates, but the target values for electrode fabrication remain the same as those when optimizing the electrode for the full depth of discharge.

25 ENERGY STORAGE↗

Review—Meeting Fuel Cell Catalyst Requirements for Heavy-Duty Vehicle Applications

Catalyst requirements for proton exchange membrane (PEM) fuel cells differ by applications. Commercial heavy-duty vehicle (HDV) applications consume more H 2 fuel and demand higher durability than many others and the total cost of ownership (TCO) of the vehicle is largely related to the performance and durability of catalysts. This article is written to bridge the gap between the industrial requirements and academic activity for advanced cathode catalysts with an emphasis on durability. From a materials perspective, the underlying nature of the carbon support, Pt-alloy crystal structure, stability of the alloying element, cathode ionomer volume fraction, and catalyst-ionomer interface play a critical role in improving performance and durability. We provide our perspective on four major approaches, namely, mesoporous carbon supports, ordered PtCo intermetallic alloys, thrifting ionomer volume fraction, and shell-protection strategies that are currently being pursued. While each approach has its merits and demerits, their key developmental needs for future are highlighted.

Ramaswamy, Nagappan (ORCID:0000000234302758)↗

The effects of Fe and Al on the phase transformations and mechanical behavior of β-Ti alloy Ti-11at.%Cr

Alloy composition is important for developing desired microstructures in beta-titanium (β-Ti) alloys during thermomechanical processing. In this work, the relatively low-cost alloying elements iron (Fe) and aluminum (Al) were added to a base Ti-11at%chromium (Cr) alloy and in-situ and ex-situ examination of the microstructural evolution and mechanical properties during tension, hardness, elevated temperature X-ray diffraction (XRD) and resonance ultrasound spectroscopy (RUS) were performed. The 400 °C XRD revealed that Ti-11at%Cr underwent β-to-ω and β-to-α transformations. Adding 0.85at% Fe reduced the volume fraction of the ω- and α-phases, and adding 5.3at% Al inhibited the β-to-ω transformation. The 400 °C RUS showed that the alloys containing the ω phase exhibited an increase in shear modulus (G) of ~140%, while the ω-free alloys exhibited an increase of only ~102%. The hardness and strength values of the ω-containing Ti-11Cr and Ti-11Cr-0.85Fe increased with increased ω-phase volume fraction when the ω/α ratio was greater than 0.25. The relationship between ω/α volume ratios and mechanical properties is discussed.

36 MATERIALS SCIENCE↗

Energy absorption of architectured PrintCast interpenetrating composites in tension

Additively manufactured (AM) metal-metal composites consisting of PrintCasted 316 L austenitic stainless-steel lattice structures infiltrated with A356 casting alloy, have recently been developed for use in high energy absorption systems with potential applications ranging from static load bearing to dynamic blast containment structures. This system has a unique mechanical behavior as the volume fraction of lattice increases showing a transition from localized to de-localized failure and dramatic increase in energy absorption capability. In this work, PrintCast A356/316 L composite tensile specimens were produced with lattice volume fractions ranging from 20 % to 50 % to capture the range of this behavior. Finite element simulations support neutron diffraction measurements of stress state. Results illustrate that in tension, the reinforcement material is in tension while the matrix support material is in compression, information offering significant insight into the transition to de-localized failure. Moreover, the simulation results provide further insight into how interfacial bonding (or lack of bonding) affects the energy absorption capabilities of the PrintCast composites.

36 MATERIALS SCIENCE↗

An automated procedure built on MTEX for reconstructing deformation twin hierarchies from electron backscattered diffraction datasets of heavily twinned microstructures

Here we present a set of algorithms built on the MTEX and MATLAB graph toolboxes for automatic reconstruction of deformation twin hierarchies from Electron Backscatter Diffraction (EBSD) datasets with a focus on developing methods for heavily twinned microstructures (twin fractions >0.5). The algorithms address key issues arising at large strains, mainly: missing twin relationships, grouping of heavily deformed grain fragments into families of similar orientation originating from a single initial grain, identification of parent fragments for large twin volume fractions, and classification of families having twin relationships with multiple families. To facilitate the development of these algorithms, large-grained ultra-high purity α-Ti deformed in compression along two directions is investigated. Graphs are utilized to handle non-local geometric merging and to represent relationships throughout the reconstruction process. When determining if a grain fragment is from the undeformed microstructure, the combined metrics of the fragment's orientation volume fraction in the initial texture and the directed graph centrality measure of out-closeness (the number of nodes reached in a graph from a given node) are essential. To address automation in reconstructing the sequence of twinning and relating fragments originating from a single grain in the initial microstructure, the twin family tree is formulated as a minimum spanning tree emanating from the initial grain family. A scheme constructing the distances associated with twin relationship comprising the spanning tree is developed, and a novel quasi-directional Prim spanning tree algorithm is used to determine the twin family tree. The procedure is demonstrated to significantly improve the level of automation in reconstructing twin hierarchies in heavily twinned microstructure compared to other methodologies in literature. The procedure can readily be applied to analyses of twinning in metals, as well as provide an approach for routinely extracting twin statistics at larger deformation levels than previously possible. Significantly, the procedure is demonstrated to be capable of identifying third generation twinning in α-Ti microstructures.

36 MATERIALS SCIENCE↗

Granular packings with sliding, rolling, and twisting friction

Intuition tells us that a rolling or spinning sphere will eventually stop due to the presence of friction and other dissipative interactions. The resistance to rolling and spinning or twisting torque that stops a sphere also changes the microstructure of a granular packing of frictional spheres by increasing the number of constraints on the degrees of freedom of motion. We perform discrete element modeling simulations to construct sphere packings implementing a range of frictional constraints under a pressure-controlled protocol. Mechanically stable packings are achievable at volume fractions and average coordination numbers as low as 0.53 and 2.5, respectively, when the particles experience high resistance to sliding, rolling, and twisting. Only when the particle model includes rolling and twisting friction were experimental volume fractions reproduced.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A numerical investigation of coal particle modeling for the inlet section of a pressurized oxy-combustion burner

Concerns over climate change have led to numerous efforts to develop low-carbon technologies, and pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in the power industry. Designing an effective burner for such technologies is vital for the development of these and other coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the pressurized oxy-combustor burner is close to or even higher than the maximum limit for commercial CFD codes (e.g., ANSYS FLUENT). At these high particle volume fractions, the interactions among particles, fluid flow and wall need to be re-evaluated. The present computational work is a first step in a systematic analysis of the influence of various parameters, like method of particle release, release location, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In the POC process, pulverized coal is burned under elevated pressure in an O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing Ansys FLUENT commercial sodtware using the Reynolds-Averaged Navier-Stokes (RANS) approach. It is revealed that for this pilot-scale, pressurized burner, velocity profiles in the near-wall region exhibit some anomalies. In order to investigate the influence of particle loading in the near-wall region, particle release location was investigated. The numerical simulations incorporate the coupling between the turbulent flow and the particles. A sensitivity investigation of particle release location found that by tuning the release location, the velocity profile can be consistent with the pure gas flow velocity profile. More importantly, the particle releasing location also affects flame stability. Particle size was also found to have a significant impact on particle trajectory, flame stability and temperature. Finally, Large Eddy Simulations (LES) were performed and compared with the results from two-dimensional RANS.

Li, Lei↗

Effect of Cooling Rate on Nano-Eutectic Formation in Laser Surface Remelted and Rare Earth Modified Hypereutectic Al-20Si Alloys

Laser Surface Remelting (LSR) was applied to arc-melted Al-20Si-0.2Sr, Al-20Si-0.2Ce, and Al-20Si hypereutectic alloys to refine microstructures. Experiments revealed that microstructures in the melt pool varied from fully eutectic to a mixture of Al dendrites and inter-dendritic eutectic. We calculated cooling rates using the Eagar-Tsai model and correlated cooling rates with characteristic microstructures, revealing that a cooling rate on the order of 104 K/s could lead to maximized fully eutectic microstructure morphology. Due to rapid solidification, the Si composition in the LSR eutectic was measured at 18.2 wt.%, higher than the equilibrium eutectic composition of 12.6 wt.%Si. Compared to Al-20Si, Ce addition had no significant effect on the volume fraction of the fully eutectic structure but refined Si fibers to approximately 30 nm in diameter. Sr addition did not further refine the diameter of eutectic Si fibers compared to Al-20Si but increased the volume fraction of the fully eutectic microstructure morphology. The refinement ratio (φ) of the Si fiber diameter from the bottom of the melt pool to the surface for the three alloys was similar, at around 28%. The established correlation between the cooling rate and the size and morphology of the microstructure within the melt pool will enable tailoring of the microstructure in laser-processed as well as deposited alloys for high strength and plasticity.

36 MATERIALS SCIENCE↗

Development of a TRISO assessment case for BISON based on transient experiments in the NSRR

The tristructural isotropic (TRISO) fuel assessment and validation database currently available in BISON is based mainly on steady-state irradiation and high-temperature furnace testing. Transient assessment cases are needed to support ongoing U.S. industry efforts to design and deploy commercial reactors that utilize TRISO fuels. Most available historical transient tests involving TRISO fuels were characterized by power densities, temperatures, and energy depositions that were highly conservative with respect to typical high-temperature gas-cooled reactor (HTGR) accident conditions. Nonetheless, assessment of BISON's predictive capabilities, adaptation of material properties toward high-temperature and high-particle-power regimes, and development of a systematic approach to validating BISON for TRISO transient applications now motivate the modeling of historical transient tests. This work describes the development of 1-D thermal models of transient experiments conducted at the Nuclear Safety Research Reactor (NSRR) in Japan, with BISON predictions of energy deposition, UO2 melting onset, and molten volume fractions being compared against experimental measurements. Melting is accounted for by defining an effective specific heat capacity for UO2—one that leverages the material's heat of melting. BISON's predictions of energy deposition and molten volume fraction are in good agreement with the experimental data from the low-energy-deposition tests; however, BISON tends to overpredict melting at higher energy depositions. Though overly conservative compared to the operating conditions expected for near-term TRISO-fueled reactor applications, these simulations effectively exercise BISON TRISO models over a wider range of conditions than those encompassed by the existing assessment database. Potential contributors to the observed discrepancies are noted and additional future developments proposed.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗