Search NASA⌕ Search

SEARCH · Search NASA

Results for “critical design”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Heterogeneous hydroformylation of alkenes by Rh-based catalysts

Heterogeneously catalyzed hydroformylation reactions avert the problems of separation in homogeneous reactions, reduce the discharge of phosphorus-containing wastes, and prevent the loss of precious metals. This review summarizes the developments of catalysts for heterogeneous hydroformylation in the last decade and proposes potential modification methods for future improvements. The hydroformylation properties of different reactants are compared in order to emphasize the fact that various olefins face unique challenges. The hydroformylation of ethylene requires the enhancement of catalytic activity and chemoselectivity, whereas that of propylene needs improvements in both chemoselectivity and regioselectivity. The hydroformylation of longer-chain olefins and styrene seeks further improvement in regioselectivity. The catalytic properties of activity, selectivity, and stability can be tuned through electronic effects (bimetallic effects, phosphine ligands, and single-atom catalysts), steric hindrance (encapsulation structures, diphosphine ligands, and local functional groups), and reaction condition optimization (temperature, pressure, and solvents). Furthermore, the understanding of these effects is critical in future catalyst design for heterogeneous hydroformylation reactions.

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

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗
DOE OSTI2023

Investigation of 3D printed lightweight hybrid composites via theoretical modeling and machine learning

Hybrid composites combine two or more different fillers to achieve multifunctional or advanced material properties, such as lightweight and enhanced mechanical properties. The properties of the composites significantly depend on their microstructures, which can be tailored via advanced 3D printing processes. Understanding the process-structure-property relationships is critical to enable the design and engineering of novel hybrid composites for applications in aerospace, automotive, and protective coatings. Here, for this work, we develop 3D printable and lightweight hybrid composites and leverage the conventional design of experiments, a theoretical hybrid model, and an image-driven machine learning (ML) method to investigate their mechanical behaviors. The hybrid composites are formulated with elastomer matrix, microfillers, and thin-shell particles, enabling a significant degree of design freedom of microstructures with densities and mechanical properties varying up to 70% and 91%, respectively. Our statistical analysis indicates that the 3D printing path direction and the microfibers fraction are dominating process parameters with contribution percentages of 45.3% and 57.7% on the specific stiffness and strength, respectively. A hybrid mechanics model is developed based on a simple Weibull distribution function and classical single-filler models to effectively capture the variations in mechanical properties, however, it overestimates the values due to its statistical constraints and idealization of experimental uncertainty. The image-driven ML model leverages the microscale images directly without losing the structural details, shows more accurate predictions with experimental data, and has 48.6% lower root mean square error than the theoretical model.

3D printing↗
DOE OSTI2026

Accelerated creep profiling: a high-throughput thermal- and stress-gradient approach applied to additive and wrought stainless steel

Here, this study investigates logarithmic creep behavior at temperatures near room temperature in wrought and additively manufactured (AM) 316 L austenitic stainless steel. A novel high-throughput (HiTp) methodology employed thermal gradients to obtain stress- and temperature-dependence of creep behavior. The findings, based on 104 two-week creep tests, reveal significant differences in creep response between wrought and AM materials. Compared to the wrought material, the AM material exhibited a higher activation energy and a lower activation volume. These differences were attributed to the materials’ distinct microstructural features, e.g., dislocation cell structures and precipitate distributions, which alter the susceptibility to dislocation glide creep, particularly at stresses below the yield stress. The HiTp methodology provides both an accelerated determination of low-temperature creep mechanisms, offering critical data for the design and optimization of alloys in near-ambient, high stress applications.

316 L Stainless steel↗
DOE OSTI2023

Graded microstructure and mechanical properties of spark plasma sintered Fe-Cr alloys

Developing graded microstructure and mechanical properties is critical for accelerating the design and optimization of structural materials for a wide range of applications. Here, in this study, Fe-Cr alloys were fabricated by spark plasma sintering (SPS) technique. The microstructure, microhardness, and tensile properties of the as-fabricated and thermally annealed variants were investigated. Graded grain structures were created in the as-fabricated cylinders along both axial and radial directions. Grain size was gradually reduced from the sample periphery to the center. Microhardness measurements reveal a gradual decrease of hardness towards the periphery of the samples. The as-fabricated Fe-Cr alloys show a desired combination of tensile strength and elongation, primarily due to the formation of a high density of oxides, voids, dislocations, and grain boundaries. After in-situ thermal annealing at 600 °C under 60 MPa for 2 h using SPS, the Fe-Cr alloys underwent minimal grain growth, and the graded grain structures were retained. The hardness was more uniformly distributed in the annealed variants, and the tensile strength was reduced with an increase in the total elongation, which is attributed to the dissolution of nano-sized oxide particles and the relief of the residual stress. This study demonstrates that SPS coupled with subsequent heat treatment can tailor the graded microstructure and control the mechanical properties of Fe-Cr alloys, showing potential applications in other alloy systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing surface-adsorbate interactions through active particle dynamics

Adsorbate molecules present in a reaction mixture may bind to and block catalytic sites. Measurement of the surface coverage of these molecules via adsorption isotherms is critical for modeling and design of catalytic reactions on surfaces. However, it is challenging to measure isotherms in solution in a way that is directly relevant to catalytic activity under reaction conditions, particularly since adsorbates may bind with an enormous range of surface affinity parameters. Here we used the motion of self-propelled catalytic Janus particles, which employ the decomposition of hydrogen peroxide fuel as a propulsion mechanism, to determine the effective surface coverage of thioglycerol, furfural, and ethanol on a platinum surface as a function of concentration in aqueous solution by measuring the decrease in active motion due to the blocking of active sites. For strongly adsorbing thioglycerol, this effective coverage was compared and contrasted to the total adsorbed amount measured using inductively-coupled plasma analysis. Demonstrating the broad applicability of this approach, the surface affinity of the three adsorbates spanned more than four orders of magnitude. For each species, the adsorbate-mediated attenuation of active motion occurred over a wide concentration range and was well-described by a Langmuir isotherm. The strongly interacting thioglycerol had the highest affinity towards the surface (K a = 15.5 ± 4.3 mM –1 ) and fully deactivated the active particle motion at surface saturation. Furfural had an intermediate affinity (K a = 0.42 ± 0.07 mM –1 ) but did not fully block H 2 O 2 access to the surface at apparent saturation, consistent with a maximum fractional surface coverage of θ max = 0.67. Ethanol exhibited even lower affinity (K a = 0.0025 ± 2x10 -4 mM –1 ) and its coverage saturated at only θ max = 0.38. Analysis of isotherms at elevated temperatures enabled direct extraction of the enthalpies of adsorption. The degree of surface coverage at adsorbate saturation appeared to correlate with the relative energies of adsorption for the different adsorbate species and was consistent with adsorbate saturation of one of multiple active site populations towards H 2 O 2 decomposition. Furthermore, computational investigations into solvent effects on furfural adsorption showed good quantitative agreement with the experimental results. This work leverages unique properties of active particles to explore fundamental catalysis questions and demonstrates a novel paradigm for significant and experimentally accessible multidisciplinary research.

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

Multiscale study of helium diffusion in Ni-Cr alloys: Short-range trapping versus long-range channeling

Ni-Cr alloys are widely employed as structural materials in fast nuclear reactors but are vulnerable to high-temperature helium (He) embrittlement (HTHE) under fast neutron irradiation. A comprehensive understanding of He diffusion in Ni-Cr alloys, which governs the kinetics of HTHE, is therefore essential for developing resilient materials and preventing failure. In this work, we reveal the underlying mechanisms of He diffusion in pure Ni and Ni-Cr alloys by integrating density functional theory (DFT) with atomic kinetic Monte Carlo (AKMC) simulations. Our findings uncover a non-monotonic dependence of He diffusivity on Cr concentration, contradicting the monotonic trends predicted by DFT-parameterized theories. At low Cr concentrations, He diffusion is dominated by short-range trapping, characterized by multiple trapping sites and a distinct mechanism within the first nearest neighbor of Cr, differing from that in pure Ni. At high Cr concentrations, these local traps become interconnected, forming long-range fast diffusion channels that enhance He mobility. The competition between localized trapping and extended channeling results in a diffusivity that first decreases, then increases with rising Cr content. These atomic-scale insights offer critical guidance for the design of radiation-tolerant Ni-based alloys. Moreover, the combined DFT-AKMC methodology and the concept of random walker diffusion through interconnected energy basins present a broadly applicable framework for studying transport phenomena in disordered systems.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Electrolyte-driven interphase stabilization in high-voltage sodium-ion full cells

Sodium-ion batteries with a high-voltage O3-type layered oxide cathode paired with a hard carbon anode can offer high energy density; however, significant interfacial instabilities driven by electrode/electrolyte reactions limit a broader industrial adoption. Localized high concentration electrolytes (LHCEs) are a rational choice as they promote salt decomposition over solvent, forming stable, inorganic-rich electrode-electrolyte interphases (EEIs). We present here a comparison of high-voltage (4.2 V) hard carbon | NaNi 1/3 Fe 1/3 Mn 1/3 O 2 pouch cells in LHCEs and in a standard carbonate-based electrolyte by (i) examining the influence of diluent choice on the electrochemical performance of LHCEs and (ii) investigating how the electrolyte chemistry affects the composition and structure of EEIs formed. Importantly, LHCEs demonstrate superior electrochemical performance, achieving 37% higher capacity after 200 cycles (119 vs. 87 mA h g -1 ) compared to the carbonate-based electrolyte. The enhanced stabilization provided by LHCEs at the interface with high-voltage sodium layered oxide cathode is revealed by gas evolution measurements obtained through online electrochemical mass spectrometry (OEMS). Time-of-flight secondary ion mass spectrometry paired with focused ion beam and advanced statistical analyses reveal that the superior performance of LHCE stems from a robust, thin cathode electrolyte interphase formed on the sodium layered oxide cathode and a homogeneous solid electrolyte interphase formation on the hard carbon anode. Furthermore, this study highlights the critical importance of electrolyte design in interphase stabilization, which plays a key role in advancing sodium-ion batteries toward commercial viability.

25 ENERGY STORAGE↗
DOE OSTI2024

Critical role of solvation on CC13 porous organic cages for design of porous liquids

Efficient carbon capture requires the design of new materials with high CO 2 selectivity and gas adsorption capacity that can be incorporated into existing industrial processes. Porous liquids (PLs) are promising candidate materials that consist of a nanoporous host and a solvent forming a liquid with permanent porosity based on exclusion of the solvent from the interior of the nanoporous host. Stable PLs are based on solvent-nanoporous host interactions, which can be evaluated through molecular simulations. Here, time- and temperature-dependent density functional theory simulations were performed between four solvents, 2-bromophenol, 4-methylphenol, 2,4-dimethylphenol, and cyclohexanone and the CC13 porous organic cage (POC) as a prototypical PL composition. Overall, minimal reactions occurred in the PL including no changes in the POC structure. Additionally, POC-solvent coordination occurred through interactions of neighboring functional groups such as methyl/bromide and hydroxyl on the solvent molecules with the POC surface. Therefore, the location rather than the number of functional groups on the solvent molecule controls the POC-solvent interactions. Additionally, the POC pore window contracted or expanded up to 8% during solvation, which correlates with the experimental solubility and static solvent-POC binding, where solvents that caused less contraction of the POC pore window increased POC solubility. Finally, these results allow for the design of optimized POC-based PL compositions based on solvent-nanoporous host binding and variation in the pore window during solvation.

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

A modeling study of ocean thermal energy conversion resource and potential environmental effects around Kailua-Kona, Hawaii

Ocean Thermal Energy Conversion (OTEC) offers a promising renewable energy solution through a heat exchange process using the temperature difference between warm surface seawater and cold deep seawater. Because accurate resource characterization is critical for the optimal design and implementation of OTEC systems, a high-resolution numerical model is employed to better characterize the OTEC resource at Kona, Hawaii. Our model provides detailed spatial and temporal variability of the thermal gradient, which is essential for assessing the viability and efficiency of OTEC systems. The model results reveal distinct patterns and dynamics not captured by existing observations or models (e.g., lower-resolution information). These findings highlight the importance of using high-resolution models for accurate predictions of thermal gradient variability, ultimately supporting more efficient and sustainable OTEC deployment. Additionally, the study investigates the impacts of mixed water discharge from OTEC plants that can cause shock to organisms living in the surface water and potentially destabilize the water column. Understanding these effects is vital for minimizing any potential negative environmental consequences and ensuring the long-term viability of OTEC operations. Further, our model improves OTEC resource characterization, which can lead to optimal design and deployment of OTEC systems. The analysis of OTEC water discharge impacts can accelerate the development of OTEC technologies, overcoming permitting/consenting challenges. These findings contribute to the broader adoption of high-resolution modeling in ocean energy resource characterization, particularly for OTEC applications.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Benefits and challenges in deployment of low global warming potential R290 refrigerant for room air conditioning equipment in California

High global warming potential gases (“high GWP”) are the fastest growing sector of greenhouse gas emissions in the world and in California and are primarily used as refrigerant gases in refrigeration and cooling equipment. Hydrofluorocarbons (HFCs) refrigerants are the dominant type of high GWP gases with GWP values thousands of times larger than CO 2 on a 100-year timescale. Refrigerant-grade propane (“R290”) has a very low GWP (GWP = 3.3) with good thermodynamic properties and good cooling equipment performance but the flammability of any leaked refrigerant makes equipment design, handling, and maintenance critical factors to manage. This paper focuses on the potential climate benefits and costs of transitioning to R290 refrigerant in small room air conditioning (AC) units, specifically window AC, packaged terminal AC/heat pumps (PTAC/PTHP), and mini-split heat pumps. Overall climate impact for a transition to all three types of air conditioning units in the 2022–2051 timeframe is found to be from 15 to 64 million metric tons of greenhouse gas (GHG) savings in California with a cost of saved CO 2 eq that ranges from $14.50 per ton of CO 2 eq saved to -$50.30 per ton of CO 2 eq saved (net savings) depending on whether the baseline refrigerant is R32 or R410A and depending on the relative energy efficiency for R290 units compared to baseline units.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Computing the shape gradient of stellarator coil complexity with respect to the plasma boundary

Coil complexity is a critical consideration in stellarator design. The traditional two-step optimization approach, in which the plasma boundary is optimized for physics properties and the coils are subsequently optimized to be consistent with this boundary, can result in plasma shapes which cannot be produced with sufficiently simple coils. To address this challenge, we propose a method to incorporate considerations of coil complexity in the optimization of the plasma boundary. Coil complexity metrics are computed from the current potential solution obtained with the REGCOIL code (Landreman, Nucl. Fusion , vol. 57, 2017, 046003). While such metrics have previously been included in derivative-free fixed-boundary optimization (Drevlak et al. , Nucl. Fusion , vol. 59, 2018, 016010), we compute the local sensitivity of these metrics with respect to perturbations of the plasma boundary using the shape gradient (Landreman & Paul, Nucl. Fusion , vol. 58, 2018, 076023). We extend REGCOIL to compute derivatives of these metrics with respect to parameters describing the plasma boundary. In keeping with previous research on winding surface optimization (Paul et al. , Nucl. Fusion , vol. 58, 2018, 076015), the shape derivatives are computed with a discrete adjoint method. In contrast with the previous work, derivatives are computed with respect to the plasma surface parameters rather than the winding surface parameters. To further reduce the resolution required to compute the shape gradient, we present a more efficient representation of the plasma surface which uses a single Fourier series to describe the radial distance from a coordinate axis and a spectrally condensed poloidal angle. This representation is advantageous over the standard cylindrical representation used in the VMEC code (Hirshman & Whitson, Phys. Fluids , vol. 26, 1983, pp. 3553–3568), as it provides a uniquely defined poloidal angle, eliminating a null space in the optimization of the plasma surface. In comparison with previous spectral condensation methods (Hirshman & Breslau, Phys. Plasmas , vol. 5, 1998, p. 2664), the modified poloidal angle is obtained algebraically rather than through the solution of a nonlinear optimization problem. The resulting shape gradient highlights features of the plasma boundary that are consistent with simple coils and can be used to couple coil and fixed-boundary optimization.

Physics↗
DOE OSTI2023

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cobalt(II) Phthalocyanine Substituents Tune the Electrocatalytic CO 2 Conversion to Methanol

Cobalt phthalocyanine (Co(II)Pc) and its derivatives are promising molecular electrocatalysts for the electrochemical reduction of CO 2 to CO and methanol (CH 3 OH). Despite increasing interest, a detailed mechanistic understanding of how ligand substituents influence catalytic activity, selectivity, and efficiency remains limited. In this study, we employ density functional theory (DFT) to systematically investigate the influence of electron-donating groups (EDGs) and electron-withdrawing groups (EWGs) on the electronic structure and redox properties of the Co(II)Pc electrocatalyst, and to elucidate CO 2 RR mechanistic pathways. Our results reveal that EWGs cause a positive shift in the reduction potentials, favor CO 2 binding over protonation of the Co metal center and promote downstream methanol formation at mild potentials. EDGs show opposite trends including favorable protonation steps, promoting a negative shift in the reduction potential, and facilitating the hydrogen evolution reaction (HER), which competes with the desired CO 2 RR pathway. Notably, CO dissociation is thermodynamically and kinetically unfavorable across all systems, positioning the redox potential versus CO dissociation energy as a key factor for methanol selectivity. Furthermore, these insights provide a predictive framework for rational catalyst design and underscore the critical role of electronic tuning in advancing molecular electrocatalysts for sustainable CO 2 conversion.

Alcohols↗
DOE OSTI2024

NCAP: Noncanonical Amino Acid Parameterization Software for CHARMM Potentials

Noncanonical Amino Acids (NCAAs) provide numerous avenues for introduction of novel functionality to peptides and proteins. NCAAs can be incorporated through solid phase synthesis or genetic code expansion in conjugation with heterologous expression of the encoded protein modification. Due to the difficulty of synthesis, wide chemical space and lack of empirically resolved structures modeling the effects of NCAA mutation is critical for rational protein design. To evaluate the structural and functional perturbations NCAAs introduce we utilize molecular potentials that describe the forces in protein structure. Most potentials such as CHARMM are designed to model canonical residues but can be parameterized in include novel NCAAs. Here, in this work, we introduce NCAP a software package to generate CHARMM compatible parameters from quantum chemical calculation. Unlike currently available tools NCAP is designed to recognize NCAA structure and automatically bridge the gap between DFT calculations and potential parameters. For our software we discuss workflow, validation against canonical parameter sets and comparison to published NCAA-protein structures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evaluation of Force Fields for Molecular Dynamics Simulations of Platinum in Bulk and Nanoparticle Forms

Here, understanding the size- and shape-dependent properties of platinum nanoparticles is critical for enabling the design of nanoparticle-based applications with optimal and potentially tunable functionality. Toward this goal, we evaluated nine different empirical potentials with the purpose of accurately modeling faceted platinum nanoparticles using molecular dynamics simulation. First, the potentials were evaluated by computing bulk and surface properties - surface energy, lattice constant, stiffness constants, and the equation of state - and comparing these to prior experimental measurements and quantum mechanics calculations. Then, the potentials were assessed in terms of the stability of cubic and icosahedral nanoparticles with faces in the {100} and {111} planes, respectively. Although none of the force fields predicts all the evaluated properties with perfect accuracy, one potential - the embedded atom method formalism with a specific parameter set - was identified as best able to model platinum in both bulk and nanoparticle forms.

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

Coarse-Grained Density Functional Theory Predictions via Deep Kernel Learning

Scalable electronic predictions are critical for soft materials design. Recently, the Electronic Coarse-Graining (ECG) method was introduced to renormalize all-atom quantum chemical (QC) predictions to coarse-grained (CG) resolutions using deep neural networks (DNNs). While DNNs can learn complex representations that prove challenging for kernel-based methods, they are susceptible to overfitting and the overconfidence of uncertainty estimations. Here, we develop ECG within a GPU-accelerated Deep Kernel Learning (DKL) framework to enable CG QC predictions using range-separated hybrid density functional theory (DFT), obtaining a 107 speedup relative to naive all-atom QC. By treating the predicted electronic properties as random Gaussian Processes, DKL incorporates CG mapping degeneracy by learning the distribution of electronic energies as a function of CG configuration. DKL-ECG accurately reproduces molecular orbital energies from range-separated DFT while facilitating efficient training via active learning using the uncertainties provided by DKL. Further, we show that while active learning algorithms enable efficient sampling of a more diverse configurational space relative to random sampling, all explored query methods exhibit comparable performance for the examined system. We attribute this result to the significant overlap of the feature space and output property distributions across multiple temperatures.

97 MATHEMATICS AND COMPUTING↗