Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VOID COEFFICIENT”

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 73 records · Page 4

High-Sensitivity Low-Energy Ion Spectroscopy with Sub-Nanometer Depth Resolution Reveals Oxidation Resistance of MoS 2 Increases with Film Density and Shear-Induced Nanostructural Modifications of the Surface

For decades, density has been attributed as a critical aspect of the structure of sputter-deposited nanocrystalline molybdenum disulfide (MoS 2 ) coatings impacting oxidation resistance and wear resistance. Despite its importance, there are few examples in the literature that explicitly investigate the relationship between the density and oxidation behaviors of MoS 2 coatings. Aging and oxidation are primary considerations for the use of MoS 2 coatings in aerospace applications as they inevitably experience prolonged storage in water and oxygen-rich environments prior to use. Oxidation that is either limited to the first few nanometers or through the bulk of the coating can result in seizure due to high initial coefficients of friction or component failure from excessive wear. High-sensitivity low-energy ion spectroscopy (HS-LEIS) and Rutherford backscattering spectrometry (RBS) are both used to understand the extent of oxidation throughout the first ~10 nanometers of the surface of pure sputtered nanocrystalline MoS 2 coatings after high-temperature aging and how it is impacted by the density of coatings as measured by RBS. Results show that low-density coatings (ρ = 3.55 g/cm 3 ) exhibit a more columnar microstructure and voiding, which act as pathways for oxidative species to penetrate and interact with edge sites, causing severe surface and subsurface oxidation. Furthermore, HS-LEIS of surfaces sheared prior to oxidation reveals that the oxidation resistance of low-density MoS 2 coatings can be significantly improved by shear-induced reorientation of the surface microstructure to a basal orientation and elimination of pathways for oxygen into the bulk through compaction of surface and subsurface voids.

36 MATERIALS SCIENCE↗

Combining Spectroscopic Ellipsometry with Transmission Spectroscopy to Derive Accurate Optical Constants for Organic Materials

A method for deriving the optical constants of organic powdered materials in the mid-infrared spectral range is introduced using variable angle spectroscopic ellipsometry and transmission spectroscopy. The approach uses pressed pellets of the powder and is applied to organic solids, which have both strong and weak infrared absorption features. Many powders have significant voids and do not press into smooth, homogenous pellets. To account for pellet non-idealities and to accurately measure both n and k, three different forms of pellets were pressed and measured: A pure analyte pellet, a mixed analyte/KBr pellet with a large analyte percentage, and a KBr transmission pellet with only a small analyte percentage. Using all three pellets in a multi-sample analysis involving both ellipsometric and transmission data, the complex refractive index (n/k) values can be derived for many organic compounds. This method is illustrated to calculate the optical constants for anhydrous lactose from 6000-400 cm-1. The transmission measurements improve the spectral fitting of weak absorption features, and the multi-sample analysis enables a better determination of the significant void space that is present in the pure pellet, leading to lower values for both n and k if not properly accounted for in the multi-oscillator model used to fit the ellipsometric data.

Optical constants, refractive index, extinction co↗

A diffusion–deformation model with damage for polymer undergoing rapid decompression failure

Green hydrogen is emerging as one of the most promising alternatives to replace fossil fuels. While hydrogen gas has a good energy density by weight, its poor energy density by volume requires it to be stored under high pressure for commercial use. The hydrogen infrastructure developed to handle this high pressure hydrogen contains multiple components made of polymers as sealing agents. Although polymers do not react chemically with hydrogen gas, they undergo a mechanical failure when high pressure hydrogen gas is suddenly depressurized. This phenomenon known as 'rapid decompression failure' occurs due to the diffusion of hydrogen through polymer and getting trapped inside pre-existing cavities or voids. Here, in this paper, a continuum mechanics-based coupled diffusion-deformation-damage model was developed to predict the hydrogen distribution, stress distribution, and damage propagation inside the polymer while it undergoes rapid decompression failure. The polymer was modeled as a hyperelastic material because it represents the nonlinear material response observed in uniaxial tensile tests perfectly. The effects of hydrogen diffusivity, pre-existing cavity size, cavity location, applied hydrogen pressure, and depressurization rate on damage initiation were studied. It was found that the coefficient of diffusion plays an important role in damage initiation and damage was mostly concentrated in the inside areas rather than near the surface. Experiments were conducted with EPDM polymer which agreed well with the predicted trends using the given model. The effect of adding carbon black and silica filler particles and plasticizer to the pure EPDM polymer was also studied. It was found that damage during RDF decreases with the addition of fillers, but increases with the addition of the plasticizer. Finally, the damage evolution in the presence of two cavities was also studied, and was found that the interaction of stress fields around the cavities alters the damage occurring during RDF.

36 MATERIALS SCIENCE↗

Effect of Cyclic Thermo-Mechanical Loads on Fatigue Reliability in Polymer Matrix Composites

A methodology to compute probabilistic fatigue life of polymer matrix laminated composites has been developed and demonstrated. Matrix degradation effects caused by long term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress dependent multi-factor interaction relationship developed at NASA Lewis Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability- based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/- 45/90)(sub s) graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical cyclic loads and low thermal cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical cyclic loads and high thermal cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Shah, A. R.↗

Probabilistic Simulation for Combined Cycle Fatigue in Composites

A methodology to compute probabilistic fatigue life of polymer matrix laminated composites has been developed and demonstrated. Matrix degradation effects caused by long term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress dependent multifactor interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/- 45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical cyclic loads and low thermal cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical cyclic loads and high thermal cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

Probabilistic Simulation of Combined Thermo-Mechanical Cyclic Fatigue in Composites

A methodology to compute probabilistically-combined thermo-mechanical fatigue life of polymer matrix laminated composites has been developed and is demonstrated. Matrix degradation effects caused by long-term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress-dependent multifactor-interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability-integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/-45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical-cyclic loads and low thermal-cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical-cyclic loads and high thermal-cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

Probabilistic Simulation of Combined Thermo-Mechanical Cyclic Fatigue in Composites

A methodology to compute probabilistically-combined thermo-mechanical fatigue life of polymer matrix laminated composites has been developed and is demonstrated. Matrix degradation effects caused by long-term environmental exposure and mechanical/thermal cyclic loads are accounted for in the simulation process. A unified time-temperature-stress-dependent multifactor-interaction relationship developed at NASA Glenn Research Center has been used to model the degradation/aging of material properties due to cyclic loads. The fast probability-integration method is used to compute probabilistic distribution of response. Sensitivities of fatigue life reliability to uncertainties in the primitive random variables (e.g., constituent properties, fiber volume ratio, void volume ratio, ply thickness, etc.) computed and their significance in the reliability-based design for maximum life is discussed. The effect of variation in the thermal cyclic loads on the fatigue reliability for a (0/+/-45/90)s graphite/epoxy laminate with a ply thickness of 0.127 mm, with respect to impending failure modes has been studied. The results show that, at low mechanical-cyclic loads and low thermal-cyclic amplitudes, fatigue life for 0.999 reliability is most sensitive to matrix compressive strength, matrix modulus, thermal expansion coefficient, and ply thickness. Whereas at high mechanical-cyclic loads and high thermal-cyclic amplitudes, fatigue life at 0.999 reliability is more sensitive to the shear strength of matrix, longitudinal fiber modulus, matrix modulus, and ply thickness.

Chamis, Christos C.↗

Near-surface temperature gradients and their effects on thermal-infrared emission spectra of particulate planetary surfaces

The infrared energy emitted from a planetary surface is generated within a finite depth determined by the material's absorption skin depth. This parameter varies significantly with wavelength in the infrared but has an average value of around 50 microns for most geologic materials. In solid rock, heat transfer is efficient enough so that this 50 micron zone of the near surface from which the radiation emanates will be more or less isothermal. In particulate materials, however, heat transfer is more complicated and occurs via a combination of mechanisms, including solid conduction within grains and across grain contacts, conduction through the interstitial gas, and thermal radiation within individual particles and across the void spaces in between grains. On planets with substantial atmospheres, the gas component dominates the heat transfer and tends to mitigate near-surface thermal gradients. However, on airless bodies, the gas component is absent and heat transfer occurs via solid conductions and radiation. If the particles are small relative to the average absorption skin depth, then the top 50-100 microns or so of the surface will be cooled by radiation to space allowing the creation of significant near-surface thermal gradients. In those regions of the spectrum where the absorption coefficient is low, the emission will come from the deeper, warmer parts of the medium, whereas in regions of high absorption, the emission will emanate from shallower, cooler parts of the medium. The resulting emission spectrum will show non-compositional features as a result of the thermal structure in the material. We have modeled the heat transfer in a particulate medium in order to determine the magnitude of near-surface thermal gradients for surfaces on airless bodies and on Mars. We use the calculated thermal structure to determine the effects it has on the infrared emission spectrum of the surface.

Henderson, B. G.↗

Enhanced flow boiling heat transfer on chromium coated zircaloy-4 using cold spray technique for accident tolerant fuel (ATF) materials

Flow boiling heat transfer tests were conducted to evaluate the Critical Heat Flux (CHF) and Heat Transfer Coefficient (HTC) of conventional and accident tolerant fuel (ATF) cladding materials, i.e., bare Zircaloy-4 (Zr4), and Zircaloy-4 coated with Chromium using physical vapor deposition (PVD) (Zr4-Cr-PVD) and cold spray process (Zr4-Cr-CS), respectively. The tests were performed on a single heater rod with a uniform heat flux profile, at the atmospheric pressure, inlet temperature (24 °C) and mass flow rate (750 kg/m 2 s). HTC’s were increased by 5.2% and 3.3% on Zr4-Cr-CS compared to bare Zr4 and Zr4-Cr-PVD. Improved HTCwas attributed to the increased void fraction and high roughness (Ra = 532 nm) of the Zr4-Cr-CS. In detail, the micro-cavities, which act as bubble seeds, were entirely distributed on the Zr4-Cr-CS, and these micro-cavities generated smaller and faster bubbles, thus the void fraction increased by 10.9% compared to bare Zr4. HTC is improved by these smaller and faster bubbles which can remove heat from the heater rod surface more efficiently. To identify the exact location of CHF, the surface temperature profile was measured using advanced fiber-optic sensors which have high temporal/spatial resolution (distance between point to point: ~2.5 mm, frequency ~100 Hz). CHF occurred at 80–95% along the heated length and showed an 11.6% reduction on Zr4-Cr-CS compared to bare Zr4. We hypothesize the CHF reduction was caused by the lower wettability of the Zr4-Cr-CS which reduces the liquid supply to the surface, and vigorous bubble accumulation owing to high void fraction near the outlet, which results in early formation of the vapor film. During the post-CHF quenching phase, the rougher Zr4-Cr-CS showed larger cooling rates compared to the bare Zr4 and Zr4-Cr-PVD, preventing of oxidation by chromium layer. In conclusion, we demonstrate that the high roughness on the Zr4-Cr-CS increased the number of micro-cavities on the surface, resulting in 5.2% increase in HTC and improvement in the quenching heat transfer performance whereas CHF was 11.6% reduced compared to bare Zr4. Overall statistically there is little effect in the coating with regards to CHF and slight improvement in HTC.

42 ENGINEERING↗

High-Strain Rate Spall Strength Measurement for CoCrFeMnNi High-Entropy Alloy

In this study, we experimentally investigate the high stain rate and spall behavior of Cantor high-entropy alloy (HEA), CoCrFeMnNi. First, the Hugoniot equations of state (EOS) for the samples are determined using laser-driven CoCrFeMnNi flyers launched into known Lithium Fluoride (LiF) windows. Photon Doppler Velocimetry (PDV) recordings of the velocity profiles find the EOS coefficients using an impedance mismatch technique. Following this set of measurements, laser-driven aluminum flyer plates are accelerated to velocities of 0.5–1.0 km/s using a high-energy pulse laser. Upon impact with CoCrFeMnNi samples, the shock response is found through PDV measurements of the free surface velocities. From this second set of measurements, the spall strength of the alloy is found for pressures up to 5 GPa and strain rates in excess of 106 s−1. Further analysis of the failure mechanisms behind the spallation is conducted using fractography revealing the occurrence of ductile fracture at voids presumed to be caused by chromium oxide deposits created during the manufacturing process.

36 MATERIALS SCIENCE↗

Hygrothermal Aging and Thermomechanical Characterization of As-Manufactured Tidal Turbine Blade Composites

This study investigates the hygrothermal aging behavior and thermomechanical properties of as-manufactured glass fiber-reinforced epoxy and thermoplastic composite tidal turbine blades. The blades were previously deployed in a marine environment and subsequently analyzed through a comprehensive suite of material characterization techniques, including hygrothermal aging, dynamic mechanical analysis (DMA), tensile testing and X-ray computed tomography (XCT). Hygrothermal aging experiments revealed that while thermoplastic composites exhibited lower overall water absorption (0.78% vs. 0.47%), they had significantly higher diffusion coefficients than epoxy (2.1 vs. 12.1 × 10 −13 m 2 s −1 ), suggesting faster saturation in operational environments. DMA results demonstrated that water ingress caused plasticization in epoxy matrices, reducing the glass transition temperature and increasing damping (112 °C to 104 °C), while thermoplastic composites showed more stable thermal behavior (87 °C glass transition temperature). Tensile testing revealed substantial reductions in ultimate strength (>40%) for both materials after prolonged water exposure, with minimal change in elastic modulus, highlighting the role of matrix degradation over fiber reinforcement. XCT image analysis showed that both composites were manufactured with high quality: no large voids or cracks were present, and the degree of misalignment was low. These findings inform future marine renewable energy composite designs by emphasizing the critical influence of moisture on long-term structural integrity and the need for optimized material systems in harsh marine environments. This work provides a rare real-world comparison of epoxy and recyclable thermoplastic tidal turbine blades, showing how laboratory aging tests and advanced imaging reveal the influence of material and manufacturing choices on long-term marine durability.

16 TIDAL AND WAVE POWER↗

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗

An electromagnetic noncontacting sensor for thickness measurement in a dispersive medium

This paper describes a general purpose imaging technology developed by the U.S. Bureau of Mines (USBM) that, when fully implemented, will solve the general problem of 'seeing into the earth.' A first-generation radar coal thickness sensor, the RCTS-1, has been developed and field-tested in both underground and highwall mines. The noncontacting electromagnetic technique uses spatial modulation created by moving a simple sensor antenna in a direction along each axis to be measured while the complex reflection coefficient is measured at multiple frequencies over a two-to-one bandwidth. The antenna motion imparts spatial modulation to the data that enables signal processing to solve the problems of media, target, and antenna dispersion. Knowledge of the dielectric constant of the media is not necessary because the electrical properties of the media are determined automatically along with the distance to the target and thickness of each layer of the target. The sensor was developed as a navigation guidance sensor to accurately detect the coal/noncoal interface required for the USBM computer-assisted mining machine program. Other mining applications include the location of rock fractures, water-filled voids, and abandoned gas wells. These hazards can be detected in advance of the mining operation. This initiating technology is being expanded into a full three-dimensional (3-D) imaging system that will have applications in both the underground and surface environment.

Chufo, Robert L.↗

Probabilistic Micromechanics and Macromechanics for Ceramic Matrix Composites

The properties of ceramic matrix composites (CMC's) are known to display a considerable amount of scatter due to variations in fiber/matrix properties, interphase properties, interphase bonding, amount of matrix voids, and many geometry- or fabrication-related parameters, such as ply thickness and ply orientation. This paper summarizes preliminary studies in which formal probabilistic descriptions of the material-behavior- and fabrication-related parameters were incorporated into micromechanics and macromechanics for CMC'S. In this process two existing methodologies, namely CMC micromechanics and macromechanics analysis and a fast probability integration (FPI) technique are synergistically coupled to obtain the probabilistic composite behavior or response. Preliminary results in the form of cumulative probability distributions and information on the probability sensitivities of the response to primitive variables for a unidirectional silicon carbide/reaction-bonded silicon nitride (SiC/RBSN) CMC are presented. The cumulative distribution functions are computed for composite moduli, thermal expansion coefficients, thermal conductivities, and longitudinal tensile strength at room temperature. The variations in the constituent properties that directly affect these composite properties are accounted for via assumed probabilistic distributions. Collectively, the results show that the present technique provides valuable information about the composite properties and sensitivity factors, which is useful to design or test engineers. Furthermore, the present methodology is computationally more efficient than a standard Monte-Carlo simulation technique; and the agreement between the two solutions is excellent, as shown via select examples.

Murthy, Pappu L. N.↗

Understanding radiation effects in friction stir welded MA956 using ion irradiation and a rate theory model

We report an outstanding challenge in the manufacturing and joining of oxide dispersion strengthened steels is retaining the nanofeatures in the alloy throughout the fabrication and welding process. MA956 was friction stir welded with two different sets of welding parameters, resulting in a medium and high heat input. After welding, 5 MeV Fe ++ ion irradiations were performed at doses ranging from 50 to 200 dpa in the temperature range of 400 to 500°C. Post-irradiation characterization was performed with scanning transmission electron microscopy and energy-dispersive x-ray spectroscopy to investigate the Y-Al-O dispersoids, voids, and dislocations. After welding, the dispersoid microstructure coarsened, resulting in fewer and larger dispersoids regardless of heat input. After irradiation, the dispersoid behavior in the welded material was sensitive to temperature, exhibiting growth behavior attributed to Ostwald coarsening at 500 °C but a mixture of nucleation and more muted growth at 400 and 450 °C, attributed to competing mechanisms of radiation-enhanced diffusion and Ostwald coarsening. Void swelling correlated to heat input; being more prevalent in the welded conditions occurring at lower doses and in higher values relative to the base material. The low values of swelling despite microstructure coarsening caused by welding demonstrate the excellent swelling resistance of MA956, even after welding with the highest swelling values of 0.5% noted in the stir zone high heat input condition at 450 °C, 200 dpa. The dislocation behavior was inconsistent: the strongest trend was that network density was higher for welded versus base material, and an increase in loop diameter with temperature was observed. A rate theory model based on the observed microstructure suggests at high temperature interstitial loss to sinks was more likely to be dominant compared to mutual annihilation via point defect recombination, because of an increase of the radiation diffusion coefficient with temperature regardless of initial welded microstructure.

36 MATERIALS SCIENCE↗

Computer Simulations of Small Molecules in Membranes: Insights from Computer Simulations into the Interactions of Small Molecules with Lipid Bilayers

Two of Ernest Overton's lasting contributions to biology are the Meyer-Overton relationship between the potency of an anesthetic and its solubility in oil, and the Overton rule which relates the permeability of a membrane to the oil-water partition coefficient of the permeating molecule. A growing body of experimental evidence, however, cannot be reconciled with these theories. In particular, the molecular nature of membranes, unknown to Overton, needs to be included in any description of these phenomena. Computer simulations are ideally suited for providing atomic-level information about the behavior of small molecules in membranes. The authors discuss simulation studies relevant to Overton's ideas. Through simulations it was found that anesthetics tend to concentrate at interfaces and their anesthetic potency correlates better with solubility at the water-membrane interface than with solubility in oil. Simulation studies of membrane permeation revealed the anisotropic nature of the membranes, as evidenced, for example, by the highly nonuniform distribution of free volume in the bilayer. This, in turn, influences the diffusion rates of solutes, which increase with the depth in the membrane. Small solutes tend to move by hopping between voids in the bilayer, and this hopping motion may be responsible for the deviation from the Overton rule of the permeation rates of these molecules.

Pohorille, Andrew↗

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crack Turning and Arrest Mechanisms for Integral Structure

In the course of several years of research efforts to predict crack turning and flapping in aircraft fuselage structures and other problems related to crack turning, the 2nd order maximum tangential stress theory has been identified as the theory most capable of predicting the observed test results. This theory requires knowledge of a material specific characteristic length, and also a computation of the stress intensity factors and the T-stress, or second order term in the asymptotic stress field in the vicinity of the crack tip. A characteristic length, r(sub c), is proposed for ductile materials pertaining to the onset of plastic instability, as opposed to the void spacing theories espoused by previous investigators. For the plane stress case, an approximate estimate of r(sub c), is obtained from the asymptotic field for strain hardening materials given by Hutchinson, Rice and Rosengren (HRR). A previous study using of high order finite element methods to calculate T-stresses by contour integrals resulted in extremely high accuracy values obtained for selected test specimen geometries, and a theoretical error estimation parameter was defined. In the present study, it is shown that a large portion of the error in finite element computations of both K and T are systematic, and can be corrected after the initial solution if the finite element implementation utilizes a similar crack tip discretization scheme for all problems. This scheme is applied for two-dimensional problems to a both a p-version finite element code, showing that sufficiently accurate values of both K(sub I) and T can be obtained with fairly low order elements if correction is used. T-stress correction coefficients are also developed for the singular crack tip rosette utilized in the adaptive mesh finite element code FRANC2D, and shown to reduce the error in the computed T-stress significantly. Stress intensity factor correction was not attempted for FRANC2D because it employs a highly accurate quarter-point scheme to obtain stress intensity factors.

Pettit, Richard↗