Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “crack nucleation”

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 163 records · Page 9

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Critical heat flux on zircaloy and accident tolerant fuel cladding under prototypical conditions of pressurized and boiling water reactors

First of a kind high spatial and temporal resolution temperature measurements of Pressurized Water Reactor and Boiling Water Reactor simulated fuel pins (high heat flux cosine profile heaters) in annular rod type geometry during nucleate boiling, departure from nucleate boiling and re-wetting has been performed. Distributed fiber optic measurements of the cladding axial temperature with spatial resolutions down to 2.5 mm axially at up to three equally spaced azimuthal locations per rod have demonstrated the similarities and differences of the different accident tolerant fuel cladding materials. These experiments show the performance of accident tolerant cladding with a precise measurement of the location and evolution of dryout followed by the re-wetting phenomena (see supplementary time history of a critical heat flux event). Experimental parameters varied coolant mass fluxes from 1695 to 2712 kg/m 2 s, pressures from 10 to 20 MPa, and inlet subcooling from 10 to 55 °C. It was found that the coatings and different material have little effect on the value of the critical heat flux. Any discrepancies between the local critical heat flux verified among the cladding material are attributed to different departure from nucleate boiling phenomena captured during the experiment with the high-resolution optical fiber temperature data. Substantial differences in the performance of the different cladding material to survive a critical heat flux event were however observed by post-test X-ray and optical inspection. External evident damage and internal cracks to the simulated fuel pin having the bare zircaloy cladding were observed in regions associated with the occurrence of the critical heat flux for longer periods. Such damages were associated with zirconium oxidation, and not noticed on the Cr coated zircaloy or FeCrAl claddings. These results indicate that the use of FeCrAl or a Cr coated zircaloy as cladding in Pressurized Water Reactor and Boiling Water Reactor may increase the survivability of the cladding material for short transients.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

HRTEM Imaging and Mechanistic Insights Into Carbon Nanotube Nucleation and Growth on Fe Nanocatalysts in a Thermal Plasma

Thermal plasma decomposition of natural gas is a scalable pathway for the production of hydrogen alongside high-value carbon nanotubes (CNTs). Metals evaporate from an electrode and condense to form seed nanoparticles that nucleate and grow CNTs. However, the lack of mechanistic understanding of the CNT nucleation and growth processes in thermal plasma makes control over CNT diameter, chirality, length, and yield difficult. We debundled and separated CNTs from soots produced using iron (Fe) nanocatalysts, and distributed them on monolayer graphene for high-resolution transmission electron microscopy (HRTEM) imaging to gain mechanistic insights. Full graphene encapsulation was found for relatively small Fe nanoparticles that were molten at high temperatures. Zigzag single-wall or double-wall CNTs (SWCNTs or DWCNTs) appeared to have grown out directly from the graphene covering on the conical or cylindrical bodies of small molten Fe nanodroplets with high curvature. Also, SWCNTs likely grew out from H- or O-atom etched single-wall carbon nanocones observed on conical Fe nanoparticles. A SWCNT/DWCNT could also be generated from the cracked opening of the graphene covering on a face-centered cubic (FCC) Fe nanoparticle. A simple, plausible pathway is proposed for the growth of an open, H-passivated, zigzag SWCNT involving reaction of CH 2 and CH radicals at high temperatures.

Fe nanocatalysts↗

Simulation studies for surfaces and materials strength

Investigations were carried out in the following three areas. First, a simulation study was conducted to investigate the early stages of a nucleation and growth mechanism taking place on diamond surfaces. The energetically most favorable binding sites were calculated along with the binding energies for up to three carbon atoms deposited on the planes of a diamond crystal. Secondly, calculations were carried out to analyze energies and structural properties for restructured patterns of the surface of a diamond. Two different model functions, developed recently for carbon, were employed in the calculations. Lastly, a crack propagation process for the graphitic basal plane, was investigated using a molecular dynamics technique. Interactions among carbon atoms in the system were calculated using the Tersoff potential energy function.

Halicioglu, Timur↗

Ultrafine-grained Fe-TiB 2 high-modulus nanocomposite steel with high strength and isotropic mechanical properties by laser powder bed fusion

Fe-TiB 2 metal matrix composite, also called high-modulus steels (HMSs), are of great interest for applications in fuel-efficient transportation infrastructure, aerospace, and wear industries due to their high specific stiffness and yield strength. However, conventional cast Fe-TiB 2 HMSs often contain coarse and sharp-edged TiB 2 particles which easily trigger premature cracking during loading. Here, we synthesized a Fe-TiB 2 nanocomposite HMS via laser powder bed fusion (LPBF) additive manufacturing of mixed micro-sized powders of Fe, Ti, and Fe 2 B. We investigated the microstructure formation and mechanical behavior of the Fe-TiB 2 HMS. We found that in situ chemical reaction of Ti and Fe 2 B enables the formation of TiB 2 particles at nanoscale during rapid solidification of LPBF. These nanoscale TiB 2 particles can serve as heterogeneous nucleation sites and promote the formation of ultrafine and equiaxed α-Fe grains with random crystallographic textures, which differ from many other additively manufactured (AM) metal alloys characteristic of strong crystallographic textures. As such, isotropic mechanical properties were achieved in the AM Fe-TiB 2 nanocomposite HMS with a high elastic modulus of ~ 240 GPa, an exceptional yield strength of ~ 1450 MPa, and a large plasticity of ~ 20% under compression. Quantitative analysis reveals that the high yield strength primarily originates from strengthening contributions of the ultrafine grains with an average grain size of ~450 nm, the nanoscale TiB 2 reinforcing particles of 20–180 nm, and a high density of printing-induced dislocations of the order of 10 15 m –2 . In situ synchrotron high-energy X-ray diffraction unveils the load partitioning from the softer α-Fe matrix to the stiffer and stronger TiB 2 nanoparticles, contributing to the sustained strain hardening during compression. Our work not only provides a general pathway for achieving high-performance metal matrix nanocomposites by in situ chemical reaction and precipitation of ceramic nanoparticles during additive manufacturing, but also offers mechanistic insights into the deformation mechanism of nanoparticle-reinforced HMS composites.

36 MATERIALS SCIENCE↗

Pyrolysis of high-density polyethylene: Degradation behaviors, kinetics, and product characteristics

Pyrolysis is a promising technology for converting plastic waste into valuable raw materials while offering a potential solution to the global plastic pollution crisis. In this study, the thermal pyrolysis of high-density polyethylene (HDPE) is investigated in a drop tube reactor under nearly isothermal conditions. The impact of reaction temperature and gas/volatile residence time on carbon conversion and product distribution is examined across a range of 500–900°C and 3.6–32.2s, respectively. Non-condensable gas products detected by online mass spectrometry are H 2 , CH 4 , C 2 H 4 , C 2 H 6 , C 3 H 6 , and C 3 H 8 . At elevated temperatures and prolonged residence time, H 2 yield reaches as high as 8.6 wt% of the initial HDPE mass due to intensified cracking reactions of C 2 –C 3 hydrocarbons and long-chain aliphatic compounds. Consequently, pyrolysis tars consist mainly of polycyclic aromatic hydrocarbons (PAHs) with 5–7 rings, accompanied by visible coke deposition within the reactor. HDPE decomposition to volatiles is an endothermic process and it is complete at a temperature between 492°C and 525°C, depending on the heating rate employed, from non-isothermal thermogravimetric analysis and differential scanning calorimetry (TGA-DSC) measurements. The thermal degradation of HDPE pellets follows the two-dimensional nucleation growth model for conversion levels up to 0.8 with an apparent activation energy of 259–270 kJ/mol and a pre-exponential factor of 4.83 × 10 17 –1.37 × 10 19 min -1 , determined from various isoconversional methods such as Flynn-Wall-Ozawa (FWO), Kissinger-Akahira-Sunose (KAS), and Starink, along with Criado's master plots. Further, these findings provide valuable insights into optimizing process parameters and refining reactor design for pyrolysis, which can be integrated with gasification and reforming processes to enhance hydrogen production on a larger scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of B 4 C for inhibiting crystallization in silica at high temperatures

Amorphous silica has numerous high temperature applications due to its inherent thermal shock resistance and low coefficient of thermal expansion (CTE). However, at the high temperatures required for processing (>1200 °C), the metastable β-cristobalite phase preferentially forms and is accompanied by a volume change, and potential cracking, upon conversion to the low temperature α-cristobalite phase. The CTE of the final crystalline phase is an order of magnitude higher than its amorphous counterpart. Here experimental results demonstrate that small additions of B 4 C (3.5 wt%) effectively inhibit silica crystallization in powder and sintered gel-cast forms up to 22h at temperatures as high as 1500°C as confirmed via x-ray diffraction. Further, the oxidation of B 4 C to B 2 O 3 and its subsequent melt and evaporation disrupts the nucleation and growth of the cristobalite phase. The mechanism for crystallization inhibition is further explored through optical microscopy to probe changes in surface morphology.

36 MATERIALS SCIENCE↗

A first-principles study of hydrogen surface coverage on δ-Pu (100), (111), and (110) surfaces

Hydriding corrosion of plutonium leads to surface cracking, pitting, and ultimately structural failure. Laboratory experiments demonstrate that hydriding begins on the surface or near the subsurface of plutonium. However, there has not yet been a systematic evaluation of hydrogen surface coverage on plutonium. In this work, we compute the surface energies of the low facet surfaces of face-centered cubic δ-Pu. The adsorption free energies of expected hydrogen structures at low and high coverage are presented along with the likely progression for filling sites as the H 2 partial pressure increases. Implications for near-equilibrium pressure hydride nucleation and non-equilibrium millibar pressure hydriding are discussed.

36 MATERIALS SCIENCE↗

Aging effects on the fracture toughness of SiC whisker reinforced 2XXX aluminum alloys

The effect of aging (at 150 C) time on the fracture toughness behavior of a 2XXX alloy (Al-3.55Cu-1.29Mg-0.01Fe-trace Mn) reinforced with 5 vol pct F-8 SiC whiskers was investigated by measuring hardness and electrical conductivity followed by fracture toughness tests on center-cracked specimens. The ageing time-hardening response plots showed that, independent of whisker orientation, the initial rapid increase in hardness was followed by a more gradual increase, with a broad hardness peak between 32 and 128 hrs of aging. Coincident with the hardness changes, the electrical conductivity initially decreased, reached a minimum, and then increased at aging times beyond 32 hrs. Examination by SEM indicated that the initial increase in hardness and decrease in conductivity was due to the GPB zone formation, while the subsequent increase in electrical conductivity and decrease in hardness (overaging) was due to S nucleation and growth.

Ratnaparkhi, P. L.↗

The evolution and growth of Al2O3 scales on beta-NiAl

The formation and growth of Al2O3 scales on (beta)-NiAl were studied using electron microscopy and other analytical techniques to gain an understanding of the oxidation properties of (beta)-NiAl and of alumina-forming alloys, in general. The transient and mature stages of oxidation were studied as well as the transformation stage during which the oxide scale transforms from metastable Al2O3 phases to the thermodynamically stable alpha-Al2O3 phase. The transient oxidation stages were studied at 800 deg C and for short times at 1100C. At 800C, the scales consist predominantly of delta-Al2O3 which forms by cation vacancy ordering in the defective spinel lattice of gamma-Al2O3. At 1100C, a fast-growth morphology of theta-Al2O3 forms as a surface layer over delta-Al2O3. For both oxidation temperatures, the scales are often epitaxially oriented with respect to the metal. The transient scales grow by outward cation diffusion as evidenced by surface growth morphologies. The transformation to alpha-Al2O3 occurs within 1 hour at 1100C by a nucleation and radial growth process. The large volume decrease associated with the transformation results in a highly strained alpha-Al2O3 microstructure. A change in scale growth mechanism from outward cation to inward anion diffusion allows transient surface morphologies to be smoothed by surface diffusion. The mature stage of oxidation involves the growth of an alpha-Al2O3 scale having the lacey morphology formed as a result of the gamma yields alpha transformation. Growth of the scale occurs by counterdiffusion along grain boundaries resulting in ridges formed by impingement of alpha-Al2O3 nuclei during the transformation stage. Also scale growth occurs by inward oxygen diffusion through healed cracks; the cracks result from transformation stresses. The measured growth rates of scales having the lacey morphology are an order of magnitude less than fine-grained alpha-Al2O3 scales. Metal orientations were found to have a large effect on oxide morpholoies during all stages of oxidation.

Doychak, J. K.↗

Dynamic properties of FeCrMnNi, a High Entropy Alloy

The goal of this paper was to assess the dynamic properties of a new class of materials, High entropy alloys (HEAs). Specifically, gas gun experiments coupled with recovery were performed on a specific HEA, FeCrMnNi, to measure its equation of state (EOS) and its spall strength, which is a measure of the stress required to nucleate voids under shock loading conditions. While there has been a plethora of work performed to assess the mechanical properties of these HEAs under uniaxial stress conditions as a function of strain rate, investigation of its properties in dynamic extremes remains rare. The current work fills this gap in knowledge for this novel class of materials. Our results show that the measured EOS for this material at one velocity was in reasonable agreement with an estimated Hugoniot. Furthermore, while the spall strength for this material was found to be ~1.9 GPa, with some variation based on sample location, the main failure mode wasbrittle. This brittle failure leading to formation of cracks in the material is different than the ductile failure observed in pure iron and its alloys. DOE/NV/03624--1298

36 MATERIALS SCIENCE↗

Dynamic properties of FeCrMnNi, a high entropy alloy

The goal of this paper was to assess the dynamic properties of a new class of materials, High entropy alloys (HEAs). Specifically, gas gun experiments coupled with recovery were performed on a specific HEA, FeCrMnNi, to measure its equation of state (EOS) and its spall strength, which is a measure of the stress required to nucleate voids under shock loading conditions. While there has been a plethora of work performed to assess the mechanical properties of these HEAs under uniaxial stress conditions as a function of strain rate, investigation of its properties in dynamic extremes remains rare. The current work fills this gap in knowledge for this novel class of materials. Our results show that the measured EOS for this material at one velocity was in reasonable agreement with an estimated Hugoniot. Furthermore, while the spall strength for this material was found to be ~1.9 GPa, with some variation based on sample location, the main failure mode was brittle. This brittle failure leading to formation of cracks in the material is different than the ductile failure observed in pure iron and its alloys.

36 MATERIALS SCIENCE↗

On relationship between microfissuring and microstructure in the HAZ of Inconel 718

Inconel 718, as well as many other metals, is susceptible to intergranular hot cracking in the weld heat-affected-zone (HAZ). These cracks form near the solidus temperature of the metal during the welding process. This problem is of particular concern to NASA/MSFC because the SSME is primarily constructed of welded Inconel 718. The present program studied microfissuring in the weld HAZ of Inconel 718 by simulating HAZ thermal cycles with the Gleeble machine. Previous researchers have studied the instantaneous mechanical properties of the HAZ using the Gleeble. The present study examines the instantaneous microstructure of the HAZ. This approach showed that second phase structures, high in niobium, melt and wet the grain boundaries in the HAZ during welding. It is postulated that the resultant HAZ grain boundaries, enriched in niobium, act as preferred sites for microfissure nucleation as the weld zone cools.

Thompson, R. G.↗

Toward a Theory of the Evolution of Drop Morphology and Splintering by Freezing

Abstract The drop freezing process is described by a phase-field model. Two cases are considered: when the freezing is triggered by central nucleation and when nucleation occurs on the drop surface. Depending on the environmental temperature and drop size, different morphological structures develop. Detailed dendritic growth was simulated at the first stage of drop freezing. Independent of the nucleation location, a decrease in temperature within the range from ~−5° to −25°C led to an increase in the number of dendrites and a decrease in their width and the interdendritic space. At temperatures lower than about −25°C, a planar front developed following surface nucleation, while dendrites formed a granular-like structure with small interdendritic distances following bulk nucleation. An ice shell grew in at the same time (but slower) as dendrites following surface nucleation, while it started forming once the dendrites have reached the drop surface in the case of central nucleation. The formed ice morphology at the first freezing stage predefined the splintering probability. We assume that stresses needed to break the ice shell arose from freezing of the water in the interdendritic spaces. Under this assumption, the number of possible splinters/fragments was proportional to the number of dendrites, and the maximum rate of splintering/fragmentation occurred within a temperature range of about −10° to −20°C, is in agreement with available laboratory and in situ measurements. At temperatures < −25°C, freezing did not lead to the formation of significant stresses, making splintering unlikely. The number of dendrites increased with drop size, causing a corresponding increase in the number of splinters. Examples of morphology that favors drop cracking are presented, and the duration of the freezing stages is evaluated. Sensitivity of the freezing process to the surface fluxes is discussed.

54 ENVIRONMENTAL SCIENCES↗

Characterization of the tip field of a discrete dislocation pileup for the development of physically based micromechanics

It is shown, on the basis of calculations by Eshelby et al. (1951), Armstrong et al. (1966), and Chou and Li (1969), that a single parameter, such as the force on the leading dislocation (F), the crack extension force, or the stress intensity factor, is capable of characterizing uniquely the entire tip field of a discrete dislocation pileup, including the positions of mobile dislocations behind the locked leading dislocation at the tip. Conversely, the position of the i-th mobile dislocation X(i) is related to the value of F and is capable of characterizing the entire stress, strain, and displacement fields at the tip of a discrete dislocation pileup. If the interactions between dislocations are linear elastic, the measured positions of the mobile dislocations can be used to determine the value of F, which can then be used as a quantitative measure of the strength of a dislocation barrier resisting the propagation of a microslip or the nucleation of a microfracture.

Gao, Q.↗

An assessment of spall failure modes in laser powder bed fusion fabricated stainless steel 316L with low-volume intentional porosity

Here, this paper reports on spall failure and damage modes in Laser Powder Bed Fusion fabricated Stainless Steel 316L (SS316L) with intentional levels of low-volume (1–5 vol. %) porosity and pore sizes of 200, 350, and 500 μm. The fabricated specimens were subjected to uniaxial-strain plate-impact loading at ~4.5 GPa, to initiate incipient spall failure. Analysis of velocimetry profiles measured using multi-probe photon-Doppler velocimetry coupled with post-mortem analysis of soft-recovered samples reveals local suppression of spall failure (termed as spall-dominated) as a function of porosity, as the failure mechanism transitions from spall-centered tensile stress dominated to a pore-centered microstructure-dominated damage mode involving void/crack nucleation and growth at pre-existing pores. The critical porosity level where the suppression of spall failure is first observed, as well as the spall location, is dependent on both the volume fraction and the size of the initially fabricated pores. In samples of 500 μm pore size, the suppression of spall failure is observed with as little as 1 vol. % porosity, while samples with smaller pores (200 μm) still experience spall-centered tensile stress dominated failure with higher levels (5 vol. %) of porosity. In the case of pore-centered microstructure-dominated failure, spall damage can occur but the spall plane is shifted toward the rear free surface, or more generally in areas further away from the region with pores. Highly heterogeneous deformation twinning, shear banding, grain rotation, and cracking are observed in the vicinity of pre-existing pores and expected spall failure sites.

36 MATERIALS SCIENCE↗

Coupled phase field damage and crystal plasticity analysis of intragranular fracture: The role of crystallographic orientation and voids

Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.

Aluminum alloy↗

The effects of zinc addition on the environmental stability of Al-Li alloys

It was found that relatively small addition of Zn can improve the stress corrosion cracking (SCC) resistance of Al-Li alloys. However, the mechanism by which this is accomplished is unclear. The role that Zn plays in altering the behavior of Alloy 8090 is investigated. Early results suggest that Zn additions increase the volume fraction of delta(Al3Li) precipitation and differential scanning calorimetry (DSC) on these alloys confirms this. The four alloys studied each had initial compositions lying in the 8090 window and had varying amounts of Zn added to them. Alloy 8090, like other Al-Li alloys, displays a delta' precipitate free zone (PFZ) upon artificial aging along the grain and subgrain boundaries. However Zn additions greatly decreased or eliminated a delta' PFZ after 100 hours at 160 C. This implies that the subgrain boundary precipitation kinetics are being altered and suppressed. Furthermore, there appears to be a window of Zn concentration above which a delta ' PFZ can reappear with the nucleation and growth of a currently unidentified precipitate on the boundaries. Polarization experiments were performed and the results presented. The experiments were performed in deaerated 3.5 w/o NaCl in both as received (T3) condition and at peak aging of 100 hours at 160 C. The aging profile was determined via Vickers Hardness tests.

Kilmer, Raymond J.↗