Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic materials”

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 1,351 records · Page 75

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗

In situ synchrotron tensile investigations on ultrasonic additive manufactured (UAM) zirconium

The microstructure evolution during room temperature uniaxial tensile straining of ultrasonic additive manufactured (UAM) zirconium was evaluated using the Advanced Photon Source (APS) facility. Miniature dog-bone tensile specimens of two orientations were cut from a UAM-fabricated zirconium bar for in situ synchrotron tensile tests. Wide-angle X-ray scattering (WAXS) scanning was used to unveil the changes in microstructure of the entire gauge regions throughout the straining. A series of WAXS data analysis methods were utilized to quantify both elastic and plastic deformation mechanisms within the strained specimens. Stress concentrations were identified during early stages of plastic deformation, which become candidate necking positions and eventually lead to failure. Additionally, fracture surface analysis implied that these stress concentration locations may be correlated to the fabrication defects, providing insightful guidance for future improvements of the UAM zirconium fabrication process.

36 MATERIALS SCIENCE↗

Three-dimensional strain dynamics govern the hysteresis in heterogeneous catalysis

Understanding catalysts strain dynamic behaviours is crucial for the development of cost-effective, efficient, stable and long-lasting catalysts. Here, we reveal in situ three-dimensional strain evolution of single gold nanocrystals during a catalytic CO oxidation reaction under operando conditions with coherent X-ray diffractive imaging. We report direct observation of anisotropic strain dynamics at the nanoscale, where identically crystallographically-oriented facets are qualitatively differently affected by strain leading to preferential active sites formation. Interestingly, the single nanoparticle elastic energy landscape, which we map with attojoule precision, depends on heating versus cooling cycles. The hysteresis observed at the single particle level is following the normal/inverse hysteresis loops of the catalytic performances. This approach opens a powerful avenue for studying, at the single particle level, catalytic nanomaterials and deactivation processes under operando conditions that will enable profound insights into nanoscale catalytic mechanisms.

36 MATERIALS SCIENCE↗

Crystallographic orientation and spatially resolved damage in a dispersion-hardened Al alloy

The in-situ neutron diffraction technique, in combination with both the full-field crystal elasto-viscoplastic finite element model and microstructural characterization, was used to study the deformation-induced damage anisotropy in a commercial Al alloy, subjected to uniaxial tensile and cyclic loading. The simulations capture well the crystallographic-orientation-dependent lattice strain behavior. The hard grains, e.g. those orientated with the <111> and <422> orientations parallel with the uniaxial loading direction (LD), feature large Taylor factors and seem more prone to form damage-related band structures. Their effective elastic moduli decrease with the accumulation of damage, which are different from the soft grains orientated with the <200> orientation along the LD. Correlation between the distribution of voids and that of the residual lattice strain developed after failure may exist. The maximum tensile type residual lattice strain observed after failure may be resulted from the band structure formed in the hardest <111> grains. It was revealed that the band structure triggered by the hard particles could be one of sources of damage. In addition, while the specimen was obviously damaged, a fast stress relief was evidenced after unloading from the tension, especially at the beginning of unloading. Our present investigations provide a novel method for exploring the damage mechanisms of polycrystalline materials during plastic deformation.

36 MATERIALS SCIENCE↗

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗

Domain fluctuations in a ferroelectric low-strain BaTiO 3 thin film

In this work, a ferroelectric BaTiO 3 thin film grown on a NdScO 3 substrate was studied using x-ray photon correlation spectroscopy (XPCS) to characterize thermal fluctuations near the a/b to a/c domain structure transformation present in this low-strain material, which is absent in the bulk. XPCS studies provide a direct comparison of the role of domain fluctuations in first- and second-order phase transformations. The a/b to a/c domain transformation is accompanied by a decrease in fluctuation timescales, and an increase in intensity and correlation length. Surprisingly, domain fluctuations are observed up to 25 °C above the transformation, concomitant with the growth of a/c domains and coexistence of both domain types. After a small window of stability, as the Curie temperature is approached, a/c domain fluctuations are observed, albeit slower, potentially due to the structural transformation associated with the ferroelectric to paraelectric transformation. The observed time evolution and reconfiguration of domain patterns highlight the role played by phase coexistence and elastic boundary conditions in altering fluctuation timescales in ferroelectric thin films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗

Elastic Screening of Pseudogauge Fields in Graphene

Lattice deformations in graphene couple to the low-energy electronic degrees of freedom as effective scalar and gauge fields. Using molecular dynamics simulations, we show that the optical component of the displacement field, i.e., the relative motion of different sublattices, contributes at equal order as the acoustic component and effectively screens the pseudogauge fields. In particular, we consider twisted bilayer graphene and corrugated monolayer graphene. In both cases, optical lattice displacements significantly reduce the overall magnitude of the pseudomagnetic fields. For corrugated graphene, optical contributions also reshape the pseudomagnetic field and significantly modify the electronic bands near charge neutrality. Previous studies based on continuum elasticity, which ignores this effect, have therefore systematically overestimated the strength of the strain-induced pseudomagnetic field. Furthermore, our results have important consequences for the interpretation of experiments and design of straintronic applications.

36 MATERIALS SCIENCE↗

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗

Modeling Sorption of Water and Vanadium Cations by Ion-Exchange Membranes

We develop a macroscopic model for sorption of ions and solvent by ion-exchange membranes and use it to investigate how different membranes behave when exposed to electrolytes for vanadium flow batteries. The model combines a classical expression for the Gibbs free energy of a moderately concentrated electrolyte containing several ions with a term that describes elastic swelling of a polymer membrane. Parameters describing pairwise interactions between mobile ions and fixed membrane ions were fit to published experimental data and discussed in the context of the competition for sites in ion-exchange membranes. The model provides insights into why sorption of vanadium cations by Nafion is lower than predicted by Donnan theory and why switching to an anion-exchange membrane does not dramatically reduce vanadium sorption.

Electrochemistry↗

Meso-scale modeling and damage analysis of carbon/epoxy woven fabric composite under in-plane tension and compression loadings

The mechanical properties and damage behaviors of carbon/epoxy woven fabric composite under in-plane tension and compression are studied at the meso-scale level through experiment and simulation. An efficient representative volume element (RVE) modeling method with consistent mesh, high yarn volume fraction and realistic geometry is proposed. The material constitutive laws with plasticity, tension-compression asymmetry and damage evolution are established for the three components - yarn, matrix and interface, respectively. Significantly different mechanical properties and damage evolutions are observed depending on loading conditions and initial geometry characteristics. It shows a non-linear stress-strain curve with clear transition region and intensive damage in tension, while a quasi-linear behavior up to facture is observed in compression with little damage prior to final fracture. Furthermore, compared to the constant Poisson's ratio with straining in compression, a dramatic increase in Poisson's ratio appears in tension. Simulation shows damage mechanisms including transverse damage, matrix damage and delamination, which all play critical roles in the property evolution. In particular, the rapid damage accumulation after elastic deformation destroys the strong bonds and causes the easy deformation of transverse yarns which results in the transition region and large Poisson's ratio in tension. All the mechanical behaviors and damage evolutions are well captured and explained with the current RVE model.

42 ENGINEERING↗

Technical Performance of Refractory Liners for Molten Chloride Salt Thermal Energy Storage Systems

A chloride-based molten-salt system that uses a ternary blend of MgCl2/KCl/NaCl is investigated to provide higher temperature thermal energy storage capability. Despite higher thermal stability, molten chlorides present several unique challenges, including the design of internal refractory-ceramic liners to prevent the corrosion and thermal stress of alloy tank shells. This work discusses issues and potential solutions related to containment of molten chloride salt, specifically the optimization of the refractory material at the molten salt interface (hot face). The down-selected hot face candidate limits permeation of salt through the material and forms a highly stable secondary surface phase in equilibrium with the molten salt. A mortar is created using the corrosion resistant hot face brick. Brick and mortar composite are subjected to mechanical stress/strain analysis, in order to calculate composite material properties and better inform thermomechanical models. The U.S. Department of Energy Generation 3 (DOE Gen3) program seeks to develop higher efficiency CSP plants that can provide cost-competitive, flexible power in the U.S. electric grid. The proposed Gen3 Liquid Pathway CSP plant closely resembles the configuration of current nitrate salt power towers with two-tank storage (Gen2). The differences between Gen2 and Gen3 include the types of compatible materials used in salt storage tank construction. Stainless steel loses strength at Gen3 temperatures, and although nickel superalloys would be capable of withstanding sustained high temperatures, these materials are prohibitively expensive at scale. Uninsulated tank shells also pose a significant risk as common steels are highly susceptible to chemical attack from molten chloride salt. To address these concerns, refractory-ceramic based containment materials are proposed to line the inside of the hot and cold storage tanks. In doing so, stainless or carbon steel shells may be used in construction depending on the level of insulation provided. The composition of the internal liner requires careful consideration to maximize the efficacy of multiple parameters including corrosion resistance, strength at operating temperature, durability, and cost. This is particularly true for the material at the interface with the salt, known as the "hot face", which is responsible for protecting the insulating layers between the tank shell and the hot face brick layer. The molten salt in this system is superheated over 300 °C above its freeze temperature. Therefore, unlike other industrial processes which use refractory-lined vessels, it is not expected that a freeze plane will develop in the hot face. Therefore, the hot face must be designed to withstand chemical corrosion and inhibit permeation of molten salt into the insulating layers. A down selection was performed to identify a hot face candidate best equipped to maintain thermal, mechanical, and chemical integrity when exposed to molten salt over extended periods of time. Long-duration chemical capability experiments were conducted with the down selected hot face refractory fully immersed in molten chloride salt for up to 3000 hours. The average salt penetration does not exceed 100 microns. When extrapolated to 20 and 30 years of continuous exposure, the expected salt penetration depth is approximately 2.0 mm and 2.9 mm, respectively. A magnesium-rich secondary phases develops at the salt/refractory interface. X-ray diffraction identifies the material as forsterite (Mg2SiO4), which is reported to form synthetically in molten chloride salt solutions. These results suggest the selected hot face will adequately inhibit salt permeation. While there is optimism that the hot face brick will inhibit salt permeation, mortar joints are typically the weakest point of a refractory brick liner. From a thermochemical perspective, differing thermal expansion coefficients may result in the mortar and brick to grow independent of each other, creating gaps through which molten salt can penetrate. To address this issue, NREL has developed an in-house mortar composed of the down selected hot face brick that has been shown to be compatible with the salt. Compressive stress/strain analyses have been performed on the brick/ mortar composites to generate stress/strain curves. Modulus of elasticity and Poisson's ratio of the composite may be calculated from the stress/strain curves, in order to provide more representative data to finite element mechanical models for accurate approximation of stress on the tank shell and the amount of thermal expansion expected within the tank liner.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Impact of anisotropy on TRISO fuel performance

Manufacturing of tristructural isotropic (TRISO) particles involves the deposition of pyrolytic carbon (PyC) and silicon carbide (SiC) layers using the fluidized bed chemical vapor deposition (CVD) process. The CVD process is known to generate polycrystalline layers with crystallographic textures, which imparts anisotropic thermophysical properties to the layers. Past studies have shown the risk for particle failure increases with an increase in anisotropy. The limit beyond which the anisotropy of PyC layers becomes unacceptable due to failure risk has been identified as a high-priority knowledge gap. This work presents a first systematic study on the effects of anisotropic thermal and mechanical properties on TRISO fuel performance. This computational study, performed using the fuel performance code BISON, investigates how the anisotropy in elasticity and thermal properties affect the stresses, temperature, and failure of a TRISO particle. The influence of other factors, such as operating temperature and particle geometry on the anisotropy effects, also has been analyzed. The studies utilize the recently published anisotropic elasticity and thermal behavior models for TRISO PyC and SiC layers implemented using tensors with full anisotropic capability. The spherical TRISO particles with anisotropic properties were found to have greater maximum tensile stress and significantly higher failure probability than the spherical particles with isotropic properties. In conclusion, the fuel performance predicted using these recently developed models was found to be comparable with the performance obtained using the historical models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Piezosuperconductivity: Novel effects in noncentrosymmetric superconductors

Here we study effects in noncentrosymmetric superconductors arising from their unique coupling of Cooper-pair condensate and elasticity. We show that although the much discussed Lifshitz coupling is not observable in a uniform bulk state, it strikingly endows dislocations with a fractional magnetic flux. We also predict a generation of voltage-free strain by a DC current in a P- and T-breaking Josephson junction. Viewing superconductors through the lens of higher-form symmetries we identify the Lifshitz coupling as a chemical potential for the approximately conserved winding number, drawing an analogy with pyroelectric insulators.

36 MATERIALS SCIENCE↗

Reactive extrusion of frontally polymerizing continuous carbon fiber reinforced polymer composites

The manufacturing of carbon fiber-reinforced polymer (CFRP) composites demands rapid and energy-efficient strategies. Frontal polymerization (FP) enables the manufacturing of CFRP using dicyclopentadiene (DCPD) thermoset polymer which meets these requirements. In this work, we introduce reactive extrusion of CFRP (RE-CFRP), where two rollers provide localized heat and pressure to sustain the curing reaction and the consolidation of a continuous carbon fiber tow pre-impregnated with DCPD. We study the effect of the extrusion speed, temperature, and compaction force on the properties of the produced CFRP. Mechanical testing confirms that the resulting fiber volume fraction and the elastic modulus are similar to bulk cured tows. A homogenized thermo-chemical model is developed to study the effect of the process parameters on the polymerization reaction. The process produces hollow woven composite tubes directly via extrusion and in situ curing. Overall, this process offers advantages in curing, tooling, speed, and energy.

36 MATERIALS SCIENCE↗

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sustainable Shape Memory Elastomers with Reduced Melt Viscosity and Enhanced Stiffness

Melt reactive processing of lignin with nitrile rubber is a promising approach to synthesizing shape memory materials. The strong intramolecular interactions in lignin macromolecular structures, caused by π–π stacking in aromatic rings and hydrogen bonding, often result in large phase separation or low miscibility with rubbers. In this study, we investigated the chemical and molecular characteristics, as well as the stiffness and complex viscosity, of modified kraft lignin melt-reacted with an acrylonitrile/butadiene copolymer containing 41% acrylonitrile (NBR41). To enhance the macromolecular compatibility of kraft lignin with NBR41, kraft lignin was cross-linked with poly(propylene glycol) diglycidyl ether (PPDE) and trimethylolpropane triglycidyl ether (TTE), both rich in epoxy reactive groups capable of forming chemical bonds with hydroxyl and carboxyl groups. Here, our findings demonstrate that the modification of kraft lignin with PPDE and TTE resulted in significantly increased stiffness of the composites. The elastic modulus of NBR41-Kraft lignin-PPDE and NBR41-Kraft lignin-TTE increased by 82 and 162%, respectively. Both the yield strength and Young’s modulus of these two samples showed dramatic improvements. Specifically, the yield strength and Young’s modulus of NBR41-Kraft lignin-TTE increased nearly 4 and 3-fold, respectively, compared to the control sample. Interestingly, despite significant improvements in mechanical properties, the viscosity of NBR41-Kraft lignin-PPDE was substantially lower than that of the control sample. At 210 °C and an angular frequency of 1 rad/s, the complex viscosity of NBR41-Kraft lignin was approximately 100.25 ± 4.77 kPa·s, while that of NBR41-Kraft lignin-PPDE was significantly lower at 56 ± 0.93 kPa·s. These findings were validated through Fourier transform infrared spectroscopy, scanning electron microscopy, dynamic mechanical analysis, thermal characterization, rheological tests, and quasi-elastic neutron scattering techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular insights into the complex mechanics of plant epidermal cell walls

Computational analysis of cell walls Layers of intertwined fibers make up plant cell walls. The various types of fibers respond differently to deformation. Cellulose microfibrils, for example, can stretch or curve, changing their end-to-end length, and can also slide past each other, reorient relative directions, and bundle with neighboring microfibrils. Zhang et al. developed a computational model based on observations of onion skin epidermis that describes how these complex changes in space govern cell wall mechanics. The results inform ways to engineer multifunctional fibrous materials. Science , this issue p. 706

Science & Technology - Other Topics↗