Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic fracture”

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 577 records · Page 32

Molecular dynamics simulations of ultrafast radiation induced melting at metal–semiconductor interfaces

Metal–semiconductor contacts in silicon carbide (SiC) diodes endure damages at the interface when exposed to harsh radiation environments. Due to the rapid rise in temperature and ultrafast cooling that follows the radiation impact, the structural properties of the materials can be altered through melting, recrystallization, and amorphization. A detailed understanding of the material failure modes at the interface is lacking, specifically at the nanoscale. We use molecular simulations to investigate the ultrafast melting at tungsten (W)–SiC interfaces following radiation damage and apply deep learning techniques to track the transient evolution of the local molecular structures. We show that W near the radiation track undergoes melting and, eventually, most of it recrystallizes with a noticeable degree of undercooling, while SiC is rendered permanently amorphous. The observation of local undercooling in W films is important as it can affect the device performance even before the bulk melting temperature of the material is reached. We also show that at high temperatures, the interface undergoes a fracture-like failure. The results presented here are significant in understating the different failure modes of SiC diode materials.

Ashwin Ravichandran↗

A multiscale cohesive law for carbon fiber networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. In this study, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

36 MATERIALS SCIENCE↗

Analyzing Mechanical integrity of Microwave resonant Cavity flowmeter

The microwave resonant cavity flow sensor has been proposed to address the difficulties of measuring high-temperature molten salt flow in reactor vessels. The flow sensor consists of a metallic cylindrical cavity coupled to a microwave waveguide. The cylindrical cavity has one of the flat surfaces thin enough to undergo microscopic deflection due to dynamic pressure of the exerted by the flowing fluid. Deflection of the membrane causes a shift in the microwave resonant frequency. Fluid velocity can be estimated from the calibration curve. The objective of this report is to investigate mechanical integrity of the flowmeter’s membrane. The membrane is a thin stainless-steel circular plate with a clamped edge, subjected to uniform load. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. If the maximum stress is greater than the yield strength, the plate will undergo plastic deformation, and the sensor will be disabled. A stress greater than the ultimate tensile strength will cause the plate to fracture. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Stresses in both the analytic model and COMSOL numerical model were calculated for stainless steel 316 thin circular plates under room temperature conditions, where the fluid is water. The diameter of the circular plate is approximately one inch, and thickness is 10mil. Fluid flow velocity is in the range from 0.5m/s to 1m/s. The results from the two models are similar. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Nanomechanical investigation of the interplay between pore morphology and crack orientation of amorphous silica

Porous amorphous silica (a-SiO 2 ) is of both fundamental and practical interests, as they exhibit a large specific surface area and tunable porous network. However, the brittle nature of a-SiO 2 and the presence of pre-existing cracks at both micro- and nano-scales lead to complex mechanical behavior. Here, we systematically investigate the effects of pre-existing crack and its orientation on the mechanical properties of a-SiO 2 with varying pore shapes using reactive molecular dynamics simulations. We demonstrate that pore shape will primarily influence the Young’s modulus (E) and critical energy release rate (G IC ). We further investigate the impact of pore shape and crack orientation by local characterization of the structural parameters. By defining the high stress and inter-mediate regions, the overall mechanical properties are found to be greatly influenced by the pore shape which can be reflected through the spatial distribution of von Mises stress. Overall, G IC is found to increase with the increase of ligament length (also known as pore wall thickness). Meanwhile, the effect of the pre-existing crack on the crack propagation process is confirmed by analyzing the density distribution evolution. These results highlight the interplay between pore morphology and crack orientation in controlling the fracture behaviors in brittle porous materials.

42 ENGINEERING↗

In situ characterization of foreign object damage (FOD) in environmental-barrier-coated silicon carbide (SiC) ceramic

Environmental barrier coatings (EBCs) protect advanced ceramics and ceramic matrix composites (CMCs) from oxidation and corrosion in gas turbine engine environments. Foreign object damage (FOD), where debris impact the protective coatings, is a critical hazard which limits the turbine durability. Despite previous efforts to understand FOD in EBCs, a detailed understanding of the fundamental transient damage mechanisms is still lacking. In the current work, the real-time FOD behavior of a Mullite / Silicon EBC was visualized via a dynamic synchrotron X-ray source in phase contrast imaging (PCI) configuration. Prior to the in situ FOD experiments, the microstructure and composition of the coating were respectively characterized using a scanning electron microscope (SEM) and X-ray diffraction (XRD). The variation in the properties of the debris was modeled by ~ 1.5 mm diameter Partially Stabilized Zirconia (PSZ) and Silicon Nitride (Si 3 N 4 ) spheres. A modified light-gas gun setup, synchronized with the X-ray beam, was used to propel the projectiles at velocities ranging between 300 and 355 m·s -1 . Coated samples were impacted under a fully backed support configuration and at normal incidence. Coating penetration and delamination, as well as projectile deformation at the bond coat resulted for FOD by PSZ spheres. Comparatively, projectile fracture, with subsequent rebound of fragments, as well as complete coating penetration and delamination at the substrate interface occurred for FOD by Si 3 N 4 spheres. It should be noted that only cone cracking was observed for FOD by PSZ spheres while back surface cracking was present for both projectile types. Qualitative understanding of the driving forces for the observed damage mechanisms were provided.

36 MATERIALS SCIENCE↗

The delineation and interpretation of the Earth's gravity field

The geoid and topographic fields of the central Pacific were delineated and shown to correlate closely at intermediate wavelengths (500 to 2500 km). The associated admittance shows that anomalies having wavelengths less than about 1000 km are probably supported by the elastic strength of the lithosphere. Larger wavelength anomalies are due to dynamic effects in the sublithosphere. Direct modeling of small scale convection in the asthenosphere shows that the amplitudes of observed geoid and topographic anomalies can be independently matched, but that the observed admittance cannot. Only by imposing an initial regional variation in the thermal regime is it possible to match the admittance. It is proposed that this variation may be due to differences in the onset time of convection beneath the lithosphere of different ages. That is, convection beneath thickening lithosphere is strongly dependent on the rate of thickening (V) relative to the rise time for convection. The critical Rayleigh number contains the length scale K/V, where K is thermal diffusivity. Young, fast growing lithosphere stabilizes the underlying asthenosphere unless it has an unusually low viscosity. Lithosphere of different age, separated by fracture zones, will go unstable at different times, producing regional horizontal temperature gradient that may strongly influence convection. Laboratory and numerical experiments are proposed to study this form of convection and its influence on the geoid.

Marsh, Bruce D.↗

Large-scale wind turbine structures

The purpose of this presentation is to show how structural technology was applied in the design of modern wind turbines, which were recently brought to an advanced stage of development as sources of renewable power. Wind turbine structures present many difficult problems because they are relatively slender and flexible; subject to vibration and aeroelastic instabilities; acted upon by loads which are often nondeterministic; operated continuously with little maintenance in all weather; and dominated by life-cycle cost considerations. Progress in horizontal-axis wind turbines (HAWT) development was paced by progress in the understanding of structural loads, modeling of structural dynamic response, and designing of innovative structural response. During the past 15 years a series of large HAWTs was developed. This has culminated in the recent completion of the world's largest operating wind turbine, the 3.2 MW Mod-5B power plane installed on the island of Oahu, Hawaii. Some of the applications of structures technology to wind turbine will be illustrated by referring to the Mod-5B design. First, a video overview will be presented to provide familiarization with the Mod-5B project and the important components of the wind turbine system. Next, the structural requirements for large-scale wind turbines will be discussed, emphasizing the difficult fatigue-life requirements. Finally, the procedures used to design the structure will be presented, including the use of the fracture mechanics approach for determining allowable fatigue stresses.

Spera, David A.↗

Osteoporosis Treatments Affect Bone Matrix Maturation in a Rat Model of Induced Cortical Remodeling

ABSTRACT To test how osteoporosis drugs affect bone matrix maturation during cortical bone remodeling, 72 pregnant rats were switched from a 0.4% to a 0.01% calcium diet at parturition for a 23‐day lactation period. At weaning, eight dams were sacrificed to establish baseline values, while the remaining dams were returned to 0.4% calcium and treated with vehicle (saline), sodium fluoride (NaF), zoledronic acid (ZA), or sclerostin antibody (Scl‐Ab) for either 7 or 28 days (eight animals per group per time point). Femora were examined by μCT, dynamic histomorphometry, Fourier transform infrared imaging, and three‐point bending of notched specimens. Cortical porosity decreased in all groups from baseline to day 28. Intracortical mineralizing surface (MS/BS) and mineral apposition rate (MAR), as well as the mineral‐to‐matrix ratio were unaffected by treatment, but intracortical crystallinity was increased in the ZA group at day 10 compared with vehicle. Cortical area increased in all groups over 28 days mainly because of an addition of bone at the endocortical surface. Endocortical MS/BS did not vary among the groups, but endocortical MAR was suppressed in the NaF group at day 2 and elevated in the Scl‐Ab group at day 4 compared with vehicle. Endocortical mineral‐to‐matrix ratio was increased at days 5 and 10 following NaF treatment and endocortical crystallinity was increased at day 5 following ZA treatment compared with vehicle. Fracture toughness did not differ among the groups. Thus, the treatments affected matrix maturation more strongly at the endocortical then intracortical envelope. In this model of induced remodeling, the bone formation phase is synchronized at multiple sites, facilitating study of the effects of drugs or other bone‐targeting agents on matrix maturation independent of their effects on the initiation of remodeling. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

59 BASIC BIOLOGICAL SCIENCES↗

Mineralized sclerites in the gorgonian coral Leptogorgia chilensis as a natural jamming system

The soft corals (Cnidaria, Octocorallia), a diverse group of colonial marine invertebrates, can reversibly tune their body stiffness in response to external stimuli. This capability is attributed to their dynamic skeletal systems, which consist of thousands of mineralized skeletal elements, called sclerites, embedded within a gel-like matrix that swells/deswells and unjams/jams the sclerites, thus modulating skeletal stiffness. While sclerite morphology is widely used for species identification, its role in the mechanical performance of a soft coral’s skeletal system is largely unknown. Here, we investigated structure-jamming relationships in sclerite-based skeletal architectures using the red gorgonian octocoral Leptogorgia chilensis as a model system. The sclerites of L. chilensis exhibit a shaft-like geometry with two axial branches and two sets of triradiate side branches, which are aligned with the crystallographic symmetry of the constituent magnesium-containing calcite. By combining multiscale three-dimensional (3D) structural characterization, parametric geometrical modeling, 3D printing, mechanical testing, and discrete element simulations, we demonstrate how sclerite geometry achieves a balanced jamming performance in terms of stiffness, weight, strength, and fracture resistance in comparison to alternative geometries parametrically modified from the native sclerites (e.g., changes in the length and number of side branches). Here, we also found that these performance metrics are achieved through the effective interlocking among side and axial branches, which is further enhanced by the fractal-like microscopic spikes on the branch tips. The findings in this natural jamming system offer insights for designing synthetic mechanotunable material architectures for a wide range of applications, from soft robotics to mechanical dampeners.

36 MATERIALS SCIENCE↗

Exemption of lattice collapse in Ni–MnO 2 birnessite regulated by the structural water mobility

Lattice collapse and associated mechanical fracture frequently occur in Li-intercalated metal oxide cathodes at the deep charge state upon Li-ion removal, governed by chemical compositions and the resulting electron density of the oxygen atoms. However, similar lattice collapse for metal oxide electrodes in aqueous storage and its mitigation have not been well studied. In this work, we reported the lattice collapse of MnO 2 layered birnessite during the aqueous de-sodiation process at high voltage due to the structural water motion, as evidenced by in situ XRD. Unlike non-aqueous Li-intercalated electrodes, Ni-dopants mitigated the lattice collapse of birnessite at deep charge states. Moreover, density functional theory (DFT) calculations showed that Ni doping induces charge depletion of lattice oxygen due to its higher electronegativity than Mn. This charge reduction, in turn, yields a significant decrease in electrostatic repulsion between the oxygens belonging to lattice and structural water. Classical molecular dynamics simulations based on atomic charges obtained from DFT elucidate that Ni-doping facilitates immobilization of structural water in (Ni)MnO 2 and prevents lattice collapse upon Na-ion removal. (Ni)MnO 2 exempted from lattice collapse shows an improved storage capacity relative to MnO 2 while maintaining similar cycling stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microstructure dependent shear localization in 316L SS

Shear localization is characterized by the concentration of plastic deformation within narrow bands of material and constitutes the precursor to failure in materials subjected to high strain rates. These bands have been observed experimentally, but classical continuum material models have not been able to fully capture their development. The aim of this research is to investigate the role of microstructural differences in the dynamic shear response of wrought and additively manufactured (AM) 316L SS. Compact Forced Simple Shear specimens for Digital Image Correlation (CFSS-DIC) of wrought and additively manufactured 316L SS were tested on a Split Hopkinson Pressure Bar (SHPB) in conjunction with infrared thermometry measurements. Results indicate that the microstructures of wrought and AM 316L SS exhibit differences that affect the shear response of this materials. In conclusion, CFSS-DIC samples tested on the SHPB did not exhibit recrystallization inside the shear region, fact confirmed by thermometry measurements.

36 MATERIALS SCIENCE↗

Calcium silicate hydrate colloid at different humidities: Microstructure, deformation mechanism, and mechanical properties

Calcium silicate hydrate (C-S-H) is the main hydration product of Portland cement, while fundamental questions persist regarding its colloidal structure and water content. In this work, we developed a full atomistic C-S-H colloid model and investigated the physical and mechanical properties at different gel water concentrations. Water confined in nanometer gel pores degrades the mechanical properties of C-S-H because of its hydrolysis effect. A “ductile to brittle” transition was observed in the C-S-H colloid model with increasing gel water contents, evolving from transgranular failure to intergranular fracture. The gel water confined in nanometer-sized pores displays glassy characteristics, and it acts as a lubricant and facilitates structure distortions during the deformation. We elucidated the topology effect of the hydrogen bond network on the evolution of the viscoelasticity of C-S-H at different humidities. In conclusion, this work sheds light on the role of water in C-S-H colloids from nanostructure perspectives.

36 MATERIALS SCIENCE↗

Seismoball: A Small Europa Orbiter Drop-Off Probe for Early Exploration of the Europan Surface

Recent magnetometry data received from Galileo indicate that the most likely explanation for the magnetic signature there is indeed a global conducting layer below the surface. This conducting layer is well matched by a salty, mineral rich strata beneath the Europan ice crust or a salt water ocean. Galileo imaging results show a variety of terrain types thought to contain young material; for example, lineaments, chaotic terrain, and eruption features. Additionally, Galileo images have shown indications of areas of up-welling where subsurface material periodically gets pushed to the surface due to the forces of fracturing, butting, and refreezing of the ice sheet. While Europa Orbiter will provide close-flyby high resolution images, as well as magnetometry, spectroscopy and other remote sensing data of the surface, it will not be able to provide essential engineering data like surface hardness and surface ice structure needed to support eventual landed missions. Additionally, ice chemical composition at microscopic scales can only be studied in detail through in situ instrumentation. Seismoball is a small probe designed to be injected into a surface intersect orbit around Europa. Using small reverse thrusters, the probe will be capable of nulling the high horizontal injection velocity as it approaches the 2 km surface injection altitude, thus allowing it to fall to the surface at an impact velocity of < 100m/sec (much less than the DS-2 impact velocities). The external breakaway thruster structure and crushable exterior shell absorb the impact energy while allowing the science instrument suite to remain intact. JPL has already started analyzing the entry dynamics and designing/building a small, low mass probe which will withstand the impact g-forces and fit as a 'carry-on' on board the Europa Orbiter. The probe will carry a suite of 5-6 micro-instruments for imaging the surface (both microscopic and far-field), surface and shallow subsurface ice temperatures, surface hardness, crustal dynamics and periodicity, and compositional chemistry. If selected, this flight development activity will provide a unique science opportunity and adjunct to the primary Orbiter science mission. The final flight system will be designed to accommodate orbiter mass, volume, and power interface constraints, as well as entry dynamics, g-load mitigation, and arbitrary landing orientation.

Tamppari, L.↗

Insights into controlling bacterial cellulose nanofiber film properties through balancing thermodynamic interactions and colloidal dynamics

In recent years, nanocellulose has emerged as a sustainable and environmentally friendly alternative to traditional petroleum-derived structural polymers. Sourced either from plants, algae, or bacteria, nanocellulose can be processed into colloid, gel, film and fiber forms. However, the required fundamental understanding of process parameters that govern the morphology and structure–property relationships of nanocellulose systems, from colloidal suspensions to bulk materials, has not been developed and generalized for all forms of cellulose. This further hinders the more widespread adoption of this biopolymer in applications. Our study investigates the dispersion of cellulose nanofibers (CNFs) produced by a bacterial–yeast co-culture, in solvents, highlighting the role of thermodynamic interactions in influencing their colloidal behavior. By adjusting Hansen solubility parameters, we controlled the thermodynamic relationship between CNFs and solvents across various concentrations, studying the dilute to semi-dilute regimes. Rheological measurements revealed that the threshold at which a concentration-based regime transition occurs is distinctly solvent-dependent. Complementing rheological analysis with small angle X-ray scattering and zeta potential measurements, our findings reveal that enhancing CNF–solvent interactions increases excluded volume in the dilute regime, emphasizing the importance of the balance between fiber–fiber and fiber–solvent interactions. Moreover, we investigated the transition from colloidal to solid state by creating films from dispersions with varying interaction parameters in semi-dilute regimes. Through mechanical testing and scanning electron microscopy imaging of the fracture surfaces, we highlight the significance of electrokinetic effects in such transitions, as dispersions with higher electrokinetic stabilization gave rise to stronger and tougher films despite having less favorable thermodynamic interaction parameters. Finally, our work provides insights into the thermodynamic and electrokinetic interplay that governs bacterial CNF dispersion, offering a foundation for future application and a deeper understanding of nanocellulose's colloidal and structure-property relationships.

59 BASIC BIOLOGICAL SCIENCES↗

Low-temperature failure mechanism of [001] niobium micropillars under uniaxial tension

The plasticity of body-centered cubic (bcc) metals is dependent of temperature as well as sample dimension at the micrometer scale, but the effects of cryogenic temperature on the plasticity and the related failure process in micron-sized bcc metals have not been studied under uniaxial tension. In this work, we utilized in situ cryogenic micro-tensile tests, transmission electron microscopy, and dislocation dynamic simulations to examine the plasticity and failure processes of [001]-oriented bcc niobium micropillars. Our study reveals that a strong suppression of cross-slip at low temperatures prevents dislocation multiplication and leads to a dislocation starvation state, at which no mobile dislocation exists due to the rapid annihilation of dislocations at free surfaces. New dislocations are then nucleated until stress concentration at a slip step creates a micro-crack, the propagation of which leads to catastrophic failure. As a result, this unique failure process results from the combined effects of sample dimension and temperature.

36 MATERIALS SCIENCE↗

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [111] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.

Yamakov, V.↗

Enhanced Interfacial Bonding of Graft Copolymers

To understand how thermoplastic welding strength can be tuned through chemical modifications and macromolecular topology, we combined coarse-grained molecular dynamics (MD) simulations with experimental validation. Our simulations examined the diffusion dynamics of both linear and graft polymers across representative interfaces, revealing that diffusion-controlled interdigitation follows a power law, with the exponent decreasing from 0.34 to 0.11 as grafting density increases from 7.5 to 196% (with side chains grafted to both sides of a monomer unit). The addition of side chains enhances welding efficiency, as dense bottlebrush polymers with high grafting density reach maximum rupture strength faster than linear polymers. However, their saturated rupture strength is lower. This observation is subsequently corroborated by experimental lap-shear tests comparing linear polyethylene with octene grafted polyethylene elastomers. Our MD simulations show that unlike linear polymers, where backbone entanglements dominate, the grafted side chains introduce mechanisms in addition to entanglement dilution. The rapid interdigitation of side chains creates a dense mesh of entropic van der Waals contacts, which can also enhance the film welding. Furthermore, our MD simulations reveal a brittle rupture behavior in linear and comb-like (mildly grafted) polymers, while bottlebrush (densely grafted) polymers display elastomeric behavior with a pronounced stress plateau prior to fracture. Our simulations deconvolute the influence of polymer topology on deformation behavior. The rate of polymer deformation becomes lower than the applied strain rate prior to rupture, and the onset of this deviation is progressively delayed from linear to bottlebrush polymers. This trend highlights the critical role of molecular architecture in governing the mechanical response. In conclusion, these results provide deeper insight into the underlying welding mechanisms of topological polymers and present a potential approach for mitigating the interface anisotropy that is inherent in advanced manufacturing techniques such as fused filament fabrication.

graft copolymers↗

Effects of oxygen on the tensile properties and fracture behavior of the Ta-10 W refractory alloy at quasi-static and impact strain rates

This study investigates effects of pre-charged oxygen at levels up to ∼1750 ppm on the tensile properties of the refractory alloy Ta-10 W. Tests were conducted at temperatures to 1100 °C and strain rates from quasi-static (∼10−3 s−1) to impact (∼103 s−1). At all temperatures, ductility decreased with increasing oxygen concentration and could be visually identified as a ductile-to-brittle transition (DBT). Additionally, for a given oxygen concentration, ductility was lower at impact strain rate than at quasi-static. The DBT was accompanied by a change from ductile transgranular to more brittle modes, intergranular and transgranular cleavage. The oxygen concentrations associated with the DBT depend on test temperature and strain rate. However, the transitions themselves are not sharp and cannot be defined as occurring at a specific oxygen concentration. At impact strain rates, ductility is relatively insensitive to test temperatures in the range 300–1100 °C with similar embrittling effects of oxygen; however, oxygen has a more severe embrittling effect at room temperature. This behavior is reversed at quasi-static strain rates where the embrittling effect of oxygen is more severe at elevated temperatures than at room temperature. It appears this difference is related to the kinetics of oxygen diffusion to grain boundaries and crack tips. Specimens that were not pre-charged with oxygen but tensile tested in various partial pressures of oxygen underwent dynamic embrittlement during the test, with the degree of embrittlement increasing as the partial pressure of oxygen in the environment increased. Consistent with such a mechanism, ductility increased with increasing strain rate for a given oxygen partial pressure, which would be the reverse if only intrinsic dislocation mobilities were controlling.

Miller, Roger [ORNL] (ORCID:0000000172155851)↗