Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamics (Mechanics)”

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 613 records · Page 34

Control of mechanical systems with rolling constraints: Application to dynamic control of mobile robots

There are many examples of mechanical systems that require rolling contacts between two or more rigid bodies. Rolling contacts engender nonholonomic constraints in an otherwise holonomic system. In this article, we develop a unified approach to the control of mechanical systems subject to both holonomic and nonholonomic constraints. We first present a state space realization of a constrained system. We then discuss the input-output linearization and zero dynamics of the system. This approach is applied to the dynamic control of mobile robots. Two types of control algorithms for mobile robots are investigated: trajectory tracking and path following. In each case, a smooth nonlinear feedback is obtained to achieve asymptotic input-output stability and Lagrange stability of the overall system. Simulation results are presented to demonstrate the effectiveness of the control algorithms and to compare the performane of trajectory-tracking and path-following algorithms.

Sarkar, Nilanjan↗

Investigation of Cycling Coal-Fired Power Plants Using High-Fidelity Models

The project delivers a well-integrated and validated simulation platform for cycling operation analysis in coal-fired power plant. Two critical mechanical components of the boiler island were analyzed through mechanical integrity assessment and economic benefit analysis. The current phase of the project focuses on the development of the integrated simulation infrastructure and prove its feasibility and effectiveness using two typical use cases. This integrated platform can help save a lot of engineering efforts for model development and simulation analysis. Through the real simulation scenarios in this document, it was demonstrated that using this platform, an analysis can be completed in approximately 2 days, while it could cost several weeks before using this platform. Going forward, the platform built in this project can be used for more boiler service applications, and it can be further enhanced with more functions/features to maximize its usage and benefits. 1) Extend component-level analysis with more use cases to cover all the major critical components of boiler island under cycling operations. A library of critical components can be developed and validated for typical pulverized coal-fired subcritical boiler units. 2) Develop predictive maintenance features based on the integrated models (Digital Twins) and engineering analysis procedures. Predictive maintenance enables each asset to be serviced based on forecast on life consumption and cost profile for replacing/welding the critical mechanical parts of the boiler. This minimizes the chance of unscheduled shutdowns and emergency services at much higher costs and prevent the fatal accidents in unit operations. 3) Develop and maintain a standard library for critical component analysis under flexible plant operations, which will include libraries of: process models, MI models for typical pressure parts, and economic models with typical plant operating data and ISO power trade data.

01 COAL, LIGNITE, AND PEAT↗

Geophysical Impacts and Spectroscopic Identification of a Hydrous Iron Sulfate on Icy Worlds

Over geologic time-scales, large volumes of exogenic sulfur ions from Io's plasma torus have been supplied to the surface of Europa and Ganymede, which, combined with recent interpretations of orbiter images, dynamical modeling, and surface-subsurface exchange, suggests further sulfur transport into the interior of the icy worlds. These observations motivate mixed-phase spectral modeling for interpreting orbiter spectroscopy data and determination of hydration states of candidate surface materials including hydrous sulfates. In this work, we present a combined experimental and theoretical study of the low temperature and high pressure vibrational spectral signature of the iron-sulfate monohydrate endmember, szomolnokite (FeSO 4 ·H 2 O). By employing synchrotron Fourier-transform infrared spectroscopy (FTIR) in the diamond anvil cell up to 23 GPa and down to 20 K, we explore the extreme range of pressure-temperature domains relevant to icy environments throughout our solar system and beyond. Combined with our density-functional theory quantum-mechanics molecular dynamics results, we demonstrate that experimentally observed infrared features in the O-H stretching region commonly associated with nH 2 O (n > 1) hydration states can be attributed to a pure monohydrate without the need for pressure-induced exsolved ice, other coexisting hydrous iron sulfates, or strong overtone and combination modes. We further discuss the possibility of lateral variations in density and shear properties on icy worlds associated with temperature variations and the high-pressure phases of kieserite group monohydrated sulfates.

Geosciences↗

Structural deformation and mechanical response of CrS 2 , CrSe 2 and Janus CrSSe

In the framework of density functional theory (DFT), we investigate the structural deformation, and mechanical behavior of the Janus CrSSe, which has out-of-plane structural asymmetry, with conventional transition metal dichalcogenides (TMDs) CrS 2 and CrSe 2 . The Janus CrSSe could be a potential candidate for machinable optoelectronic and piezoelectric applications. Here, we predict that these compounds are chemically, mechanically, and dynamically stable with the covalent bond between the TM(Cr) and chalcogen(X=S, Se) atoms. Due to the influence of tensile strain, the Cr–X bond length of each monolayers increases and the thickness decreases. Interestingly, the in-plane stiffness, shear and layer moduli, Poisson’s ratio, ultimate bi/uni-axial stress of Janus CrSSe are in between the values of CrS 2 and CrSe 2 monolayers. Similar to TMDs, the orientation-dependent in-plane stiffness and Poisson’s ratio demonstrates the isotropic behavior in Janus CrSSe. Furthermore, it can sustain a larger value of uni/bi-axial tensile strain with the critical strain equivalent to CrX 2 monolayers. By applying higher-order strain, we have also found average elastic–plastic behavior as expected. These findings demonstrate that the Janus CrSSe monolayer is a mechanically stable and ductile compound that maintains the hybrid behavior.

36 MATERIALS SCIENCE↗

An analysis of the flow field in the region of the ASRM field joints

The flow field in the region of a solid rocket motor field joint is very important since fluid dynamic and mechanical propellant stresses can couple to cause a motor failure at a joint. Presented here is an examination of the flow field in the region of the Advanced Solid Rocket Motor (ASRM) field joints. The analyses were performed as a first step in assessing the design of the ASRM forward and aft field joints in order to assure the proper operation of the motor prior to further development of test firing. The analyses presented here were performed by employing a two-dimensional axisymmetric assumption. Fluent/BFC, a three dimensional full Navier-Stokes flow field code, was used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical algorithm. Wall functions are used to determine the character of the laminar sublayer, and a standard kappa-epsilon turbulence model is used to close the fluid dynamic equations. The analyses performed to this date verify that the ASRM field joint design operates properly. The fluid dynamic stresses at the field joints are small due to the inherent design of the field joints. A problem observed in some other solid rocket motors is that large fluid dynamic stresses are generated at the motor joint on the downstream propellant grain due to forward facing step geometries. The design of the ASRM field joints are such that this is not a problem as shown by the analyses. Also, the analyses of the inhibitor stub left protruding into the port flow from normal propellant burn back show that more information is necessary to complete these analyses. These analyses were performed as parametric analyses in relation to the height of the inhibitor stub left protruding into the motor port. A better estimate of the amount of the inhibitor stub remaining at later burn times must be determined since the height which the inhibitor stub protrudes into the port flow drastically affects the fluid dynamics induced stresses on the propellant grain at the field joints.

Dill, Richard A.↗

Effect of Mechanical Deformation on the Dielectric Electric Field in Dynamic Umbilical Cables

Wave energy converter (WEC) systems have unique umbilical cable design requirements compared to conventional marine designs. Due to their dynamic nature, WEC umbilical cables are required to handle increased loads and larger motions over longer deployment scenarios. However, the life expectancy of dynamic umbilical cables is predominantly analyzed for mechanical fatigue, with little consideration given to the effect of mechanical stress on the response of the dielectric as it has been previously assumed that the non-polymeric components will fail prior to the polymeric components. In this work, we explore the effects of mechanical bending on the concentration of the electric field in the insulation of a conceptual three-core dynamic umbilical cable to determine the limiting requirements for life expectancy. A conceptual industrial schematic of a medium voltage dynamic umbilical cable was used as the basis for the WEC system under study. Simulations on the dynamic umbilical cable were performed using COMSOL Multiphysics®. Mechanical loads up to the calculated limit (2.64 m minimum bend radius or 388 kN transverse load) were applied to the dynamic umbilical in a step-wise fashion with one end of the umbilical cable fixed. One-way coupling in three-dimensions was conducted by first determining the engineering strains, then using the engineering strains as the initial condition for determination of the Green-Lagrange strains, and lastly calculating the electric field in the deformed dielectric. Based upon the input load, deformation of the dielectric at the minimum bend radius produced an increase in the electric field of more than 12% when compared to the undeformed umbilical cable, which may lead to dielectric breakdown of the insulation prior to mechanical failure of the umbilical cable.

Spencer, Mychal P.↗

Microscopic ergodicity breaking governs the emergence and evolution of elasticity in glass-forming nanoclay suspensions

We report a study combining x-ray photon correlation spectroscopy (XPCS) with in situ rheology to investigate the microscopic dynamics and mechanical properties of aqueous suspensions of the synthetic hectorite clay Laponite, which is composed of charged, nanometer-scale, disk-shaped particles. The suspensions, with particle concentrations ranging from 3.25 to 3.75 wt%, evolve over time from a fluid to a soft glass that displays aging behavior. The XPCS measurements characterize the localization of the particles during the formation and aging of the soft-glass state. The fraction of localized particles f 0 increases rapidly during the early formation stage and grows more slowly during subsequent aging, while the characteristic localization length r loc steadily de- creases. Despite the strongly varying rates of aging at different concentrations, both f 0 and r loc scale with the elastic shear modulus G' in a manner independent of concentration. During the later aging stage, the scaling between r loc and G' agrees quantitatively with a prediction of naive mode coupling theory. Breakdown of agreement with the theory during the early formation stage indi- cates the prevalence of dynamic heterogeneity, suggesting the soft solid forms through precursors of dynamically localized clusters.

36 MATERIALS SCIENCE↗

Investigating Bacterial-Fungal Interactions using Fungal Highway Columns in Diverse Environments and Substrates

Bacterial-fungal interactions (BFIs) play an integral role in shaping microbial community composition, biogeochemical functions, spatial dynamics, and microbial dispersal. Mycelial networks created by filamentous fungi or other filamentous microorganisms (e.g., Oomycetes) act as 'fungal highways' that can be utilized by bacteria for transport throughout heterogeneous environments, greatly facilitating their mobility and granting them access to regions that may be challenging or impossible to reach on their own (e.g., due to air pockets within the soil). Several devices and experimental protocols have been created to study these fungal highways, including fungal highway columns. The fungal highway column designed by our group can be used for a variety of in situ or in vitro applications, as well as with diverse environmental and host-associated sample types. Herein, we describe the methods for performing experiments with these columns, including designing, printing, sterilizing, and preparing the devices. The options for analyzing data obtained from the use of these devices are also discussed here, and troubleshooting advice regarding potential pitfalls associated with experiments using fungal highway columns is offered. These devices can be used to gain a more comprehensive understanding of the diversity, mechanisms, and dynamics of fungal highway BFIs to provide valuable insights into the structural and functional dynamics within complex environments (e.g., soils) and across diverse habitats in which bacteria and fungi co-exist.

59 BASIC BIOLOGICAL SCIENCES↗

Ultrafast kinetics of the antiferromagnetic-ferromagnetic phase transition in FeRh

Understanding how fast short-range interactions build up long-range order is one of the most intriguing topics in condensed matter physics. FeRh is a test specimen for studying this problem in magnetism, where the microscopic spin-spin exchange interaction is ultimately responsible for either ferro- or antiferromagnetic macroscopic order. Femtosecond laser excitation can induce ferromagnetism in antiferromagnetic FeRh, but the mechanism and dynamics of this transition are topics of intense debates. Employing double-pump THz emission spectroscopy has enabled us to dramatically increase the temporal detection window of THz emission probes of transient states without sacrificing any loss of resolution or sensitivity. It allows us to study the kinetics of emergent ferromagnetism from the femtosecond up to the nanosecond timescales in FeRh/Pt bilayers. Our results strongly suggest a latency period between the initial pump-excitation and the emission of THz radiation by ferromagnetic nuclei.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics of the Trans-Neptune Region: Apsidal Waves in the Kuiper Belt

The role of apsidal density waves propagating in a primordial trans-Neptune disk (i.e., Kuiper belt) is investigated. It is shown that Neptune launches apsidal waves at its secular resonance near 40 AU that propagate radially outward, deeper into the particle disk. The wavelength of apsidal waves is considerably longer than waves that might be launched at Lindblad resonances, because the pattern speed, g(sub s), resulting from the apsis precession of Neptune is much slower than its mean motion, Omega(sub s). If the early Kuiper belt had a sufficient surface density, sigma, the disk's wave response to Neptune's secular perturbation would have spread the disturbing torque radially over a collective scale lambda(sub *) approx. = r(2(mu)(sub d)Omega/ absolute value of r dg/dr)(sup 1/2), where mu(sub d)equivalent pi(sigma)r(exp 2)/(1 solar mass) and Omega(r) and g(r) are respectively the mean motion and precession frequency of the disk particles. This results in considerably smaller eccentricities at resonance than had the disk particles been treated as noninteracting test particles. Consequently, particles are less apt to be excited into planet-crossing orbits, implying that the erosion timescales reported by earlier test-particle simulations of the Kuiper belt may be underestimated. It is also shown that the torque the disk exerts upon the planet (due to its gravitational attraction for the disk's spiral wave pattern) damps the planet's eccentricity and further inhibits the planet's ability to erode the disk. Key words: celestial mechanics, stellar dynamics - comets: general minor planets, asteroids

Ward, William R.↗

On Practical Aspects of Variational Consistency in Contact Dynamics

Usage of contact mechanics methodologies is a pervasive modeling requirement in dynamic simulations. While for some trivial problems, solutions taken from analytical geometry are available, use of a finite element framework is common to achieve formulation generality. This work explores two dynamic contact formulations: one based on the traditional node-to-segment (NTS) approach, and a variationally consistent segment-to-segment (STS) mortar formulation. The NTS formulation employed here enforces the constraints kinematically (i.e., the interpenetration is enforced to the solver tolerance), whereas the mortar approach uses Lagrange multipliers to enforce the contact constraints. Both approaches are implemented in the open-source finite element framework Multiphysics Object-Oriented Simulation Environment (MOOSE). The results highlight two relevant contact-interface-related dynamic phenomena in finite element simulations. First, stabilization of contact constraints is discussed, taking into account the evolution of the total energy in a benchmark problem. Second, the influence of finite element discretization on both of the aforementioned contact formulations is analyzed by exercising a large-deformation example with continuous relative sliding. Variationally consistent contact approaches such as the mortar formulation lead to improved energy preservation and avoid spurious excitation of the system's frequencies. This is especially relevant in settings where inertia and vibrations are of importance.

42 ENGINEERING↗

Understanding the dynamics of vacuole trafficking and secretion for enhanced plant storage reserves

Vacuoles are essential in plants and serve roles in both protein and organelle degradation for reutilization and storage of high energy molecules such as storage proteins which are critical for seed germination under normal and stress conditions and impact biomass production and nutrition. In addition, plants secrete many different molecules such as cell wall components, lipids, waxes and others that are important for plant protection, structure and biomass. Vacuole accumulation and secretion are closely linked and highly regulated pathways within the endomembrane trafficking network. We have discovered many drug-like synthetic small molecules that can suppress vacuole trafficking and enhance secretion; other chemicals can enhance vacuole trafficking and suppress secretion. The ability of our compounds to modulate these pathways indicates that these pathways are tightly coupled and highly regulated. Thus, to achieve the goals of enhanced energy reserves or enhanced secretion it is essential to understand how this dynamic is controlled. Bioactive chemicals have provided us unique insight into natural mechanisms controlling the dynamics of vacuole trafficking and secretion that was not possible previously. One such compound, ES2, negatively modulates secretion while enhancing vacuole trafficking. This is the first compound discovered that targets an essential secretory complex known as the exocyst and will reveal important regulatory domains controlling secretion and vacuole trafficking. During the current funding period, we determined the structure of the EXO70A1 in plants and identifying its binding site. To understand the dynamics of vacuole trafficking it is key to have compounds that inhibit vacuole trafficking as well as enhancing it (ES2). We have used sophisticated image analysis to identify another novel compound ES17 that enhances secretion while negatively modulating vacuole trafficking (opposite of ES2). We have already made rapid progress. Our preliminary data indicates that ES17 very specifically targets a late stage of vesicle fusion with the tonoplast which is highly regulated. Thus, ES17 is an ideal compound to dissect the regulation of vacuole targeting dynamics. In addition, we have now synthesized a biotinylated version of ES17 that is fully functional, permitting the use of chemical proteomics to identify the cognate target with a high probability of success. We propose to investigate the ES17 target which can negatively modulate vacuole trafficking and combine this with knowledge of the ES2 target which can promote vacuole trafficking. Dissecting the mechanisms controlling the dynamics of secretion vs vacuole accumulation and the use of synthetic modulators will permit us to understand and regulate trafficking to enhance secretion to the plasma membrane and cell wall or the accumulation of biological energy reserves for increased biomass, nutritional value, stress resistance, and the secretion of specialized plant products.

59 BASIC BIOLOGICAL SCIENCES↗

A High-Speed Rotational Diamond Anvil Cell for In Situ Analysis of Hierarchical Microstructural Evolution of Metallic Alloys during Extreme Shear Deformation

High speed shear deformation is ubiquitous in engineering applications, ranging from material processing methods such as friction stir processing/extrusion and in tribological contacts. However, analyzing the microstructural evolution of materials while they are undergoing high speed shear deformation have been a long-standing challenge. This led to predominant reliance on ex situ microscopy before and after shear deformation. But ex situ microscopy lacks the ability to analyze dynamic and transient hierarchical microstructural evolution mechanisms that could occur during shear deformation of materials. Therefore, to better understand the dynamic mechanisms of mass and energy transfer in materials under shear deformation, we developed a first of its kind high-speed rotational diamond anvil cell (HS-RDAC) for synchrotron-based in situ high-energy x-ray diffraction (XRD). We studied the time resolved lattice strain evolution, XRD peak broadening and changes in spatial variation of shear deformation induced alloying in pure metal and metal alloy sheets and powder mixture using the HS-RDAC. These in situ results were combined with detailed ex situ microstructural characterization before and after the shear deformation using transmission electron microscopy and atom probe tomography, which revealed the different stages of evolution of a shear deformation induced hierarchical nanostructure. Multiscale computational simulations including computational fluid dynamics, crystal plasticity, molecular dynamic simulation and density functional theory uncovered the mechanisms behind morphological changes, evolution of defect structures and changes in driving force for shear deformation induced intermixing. In conclusion, this in situ HS-RDAC capability, in combination with ex situ microstructural characterization and computational simulations, can provide new insights into the hierarchical microstructural evolution pathway during shear deformation.

36 MATERIALS SCIENCE↗

Coordination and thermophysical properties of select trivalent lanthanides in LiCl–KCl

The coordination chemistry of various fission and decay products, such as actinides and lanthanides, are crucial to the commercial deployment of molten salt reactors as they can affect the thermophysical properties. In this report we examined the structure, coordination environment, and physical properties such as the density and the vibrational density of states for three lanthanide species, namely Ce, Eu, and Sm in the LiCl–KCl eutectic system using a combination of quantum mechanics simulations and spectroscopic experiments. Quantum mechanics molecular dynamics (QM-MD) modelling was employed to determine the physical properties of each system resulting in accurate local coordination of each species. Then, the vibrational density of states (DOS) was determined using a two-phase thermodynamic modelling which was then compared to the experimentally obtained Raman spectra of the species in molten LiCl–KCl having the eutectic composition. We find that Ce 3+ , Eu 3+ and Sm 3+ all adopt octahedral local coordination environments in the eutectic salt composition in good agreement with experimental results. Ce3+ is found to fluctuate between an octahedral six-coordinated and a seven-coordinated structure due to the increased local proximity of Cl in the eutectic salt, resulting in a lower fluidicity/diffusivity than the other trivalent lanthanides studied. The thermophysical properties of the eutectic composition with trivalent lanthanides were not significantly different from the pure eutectic salt composition, but several changes were noted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synergistic Density Functional Theory and Molecular Dynamics Approach to Elucidate PNIPAM–Water Interaction Mechanisms

This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(N-isopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.

08 HYDROGEN↗

Experimental and modeling studies of IPDI-based polyurea elastomers – The role of hard segment fraction

Segmented polyureas (PUa) are industrially important class of polymers widely used in coatings, sealant, and adhesive applications. Here, in this paper, we report synthesis, characterization, and modeling of Isophorone Diisocyanate-Diethyl-Toluene-Diamine-Polyether amine (IPDI-DETDA-PO PUa) with varied hard segment contents of 20, 30, and 40 weight percent. For each of the three materials, we study its structure and phase behavior using FTIR, DSC, and TEM, and clearly show the presence of microphase separation between the hard and soft nanodomains. We then measure the linear viscoelastic response of the PUa-s using DMA (frequency sweeps at multiple temperatures). The DMA data are shown to obey the time-temperature superposition. Finally, we develop a new micromechanical model describing the DMA results; the model describes a phase-separated PUa as two “Fractional-order Maxwell gels” branches, connected in parallel, with the first FMG branch representing the “percolated hard phase" and the second one modeling the “filled soft phase". In agreement with the earlier thermodynamic theories, the volume-fraction of the percolated hard phase is related to the hard segment weight-fraction (HSWF), defined as the combined mass of IPDI and DETDA normalized to the total mass of the polymer. The data and model are found to be in a good qualitative and quantitative agreement.

36 MATERIALS SCIENCE↗

Histidine methyltransferase SETD3 methylates structurally diverse histidine mimics in actin

Actin histidine N τ -methylation by histidine methyltransferase SETD3 plays an important role in human biology and diseases. Here, we report integrated synthetic, biocatalytic, biostructural, and computational analyses on human SETD3-catalyzed methylation of actin peptides possessing histidine and its structurally and chemically diverse mimics. Our enzyme assays supported by biostructural analyses demonstrate that SETD3 has a broader substrate scope beyond histidine, including N-nucleophiles on the aromatic and aliphatic side chains. Quantum mechanical/molecular mechanical molecular dynamics and free-energy simulations provide insight into binding geometries and the free energy barrier for the enzymatic methyl transfer to histidine mimics, further supporting experimental data that histidine is the superior SETD3 substrate over its analogs. This work demonstrates that human SETD3 has a potential to catalyze efficient methylation of several histidine mimics, overall providing mechanistic, biocatalytic, and functional insight into actin histidine methylation by SETD3.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-Step Crystallization of Self-Organized Spiral Eutectics

A method for the solidification of metallic alloys involving spiral self-organization is presented as a new strategy for producing large-area chiral patterns with emergent structural and optical properties, with attention to the underlying mechanism and dynamics. This study reports the discovery of a new growth mode for metastable, two-phase spiral patterns from a liquid metal. Crystallization proceeds via a non-classical, two-step pathway consisting of the initial formation of a polytetrahedral seed crystal, followed by ordering of two solid phases that nucleate heterogeneously on the seed and grow in a strongly coupled fashion. Crystallographic defects within the seed provide a template for spiral self-organization. These observations demonstrate the ubiquity of defect-mediated growth in multi-phase materials and establish a pathway toward bottom-up synthesis of chiral materials with an inter-phase spacing comparable to the wavelength of infrared light. Given that liquids often possess polytetrahedral short-range order, our results are applicable to many systems undergoing multi-step crystallization.

36 MATERIALS SCIENCE↗