Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cohesiveness”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Critical Shear Stresses for Erosion from Yield Stresses

The critical shear stress for erosion is important in many fields of engineering and science that involve suspension of particles from settled particle beds. However, determining the critical shear stress for erosion remains intricate, particularly where interparticle forces play a key role. Because these forces vary widely (Pease, et al., 2019), a more convenient means of estimating the critical shear stress for erosion would be advantageous. Here we evaluate the assertion that critical shear stresses for erosion may be derived from yield stresses for nominally cohesive particle beds across the full range of particle Reynolds numbers. The critical shear stress for erosion for non-cohesive particle beds may be estimated from the Shields diagrams or the analytical expressions from conservation of force about a particle at the interface. However, equivalent diagrams for cohesive particle beds have remained elusive. In 1992, pioneering analysis by Dade included the influence of the van der Waals force in particle forces balances to show that the critical shear stresses for erosion may be related to yield stresses subject to the limitations of small particle Reynolds numbers, flat particle beds, uniform particle diameters, and negligible electrostatic forces. Here we relax these restrictions. We find that the critical shear stress for erosion does relate to the yield stresses in a linear fashion for sufficiently large yield stresses. A criterion is developed to determine when the yield stresses may be important. Pease LF, AJ Fuher, JA Bamberger, and MJ Minette. 2024. Shields Diagrams for Cohesive Particle Beds. Manuscript. Dade WB, ARM Nowell and PA Jumars. 1992. “Predicting erosion resistance of muds.” Marine Geology, 105 (1992) 285-297.

shear stress, erosion, yield stress, slurry, settl↗
DOE OSTI2024

On plasticity-enhanced interfacial toughness in bonded joints

Here, the performance and reliability of many structures and components depend on the integrity of interfaces between dissimilar materials. Interfacial toughness Γ is the key material parameter that characterizes resistance to interfacial crack growth, and Γ is known to depend on many factors including temperature. For example, previous work showed that the toughness of an epoxy/aluminum interface decreased 40 % as the test temperature was increased from -60 °C to room temperature (RT). Interfacial integrity at elevated temperatures is of considerable practical importance. Recent measurements show that instead of continuing to decrease with increasing temperature, Γ increases when test temperature is above RT. Cohesive zone finite element calculations of an adhesively bonded, asymmetric double cantilever beam specimen of the type used to measure Γ suggest that this increase in toughness may be a result of R-curve behavior generated by plasticity-enhanced toughening during stable subcritical crack growth with interfacial toughness defined as the critical steady-state limit value. In these calculations, which used an elastic-perfectly plastic epoxy model with a temperature-dependent yield strength, the plasticity-enhanced increase in Γ above its intrinsic value Γ o depended on the ratio of interfacial strength σ* to the yield strength σyb of the bond material. There is a nonlinear relationship between Γ/Γ o and σ*/σ yb with the value Γ/Γ o increasing rapidly above a threshold value of σ*/σ yb . The predicted increase in toughness can be significant. For example, there is nearly a factor of two predicted increase in Γ/Γ o during micrometer-scale crack-growth when σ*/σ yb = 2 (a reasonable choice for σ*/σ yb ). Furthermore, contrary to other reported results, plasticity-enhanced toughening can occur prior to crack advance as the cohesive zone forms and the peak stress at the tip of the original crack tip translates to the tip of the fully formed cohesive zone. These results suggest that plasticity-enhanced toughening should be considered when modeling interfaces at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A robust interface finite element formulation for modeling brittle material failure problems

Abstract Failure of many brittle materials and structures can be modeled using interface‐oriented finite elements combined with intrinsic cohesive zone models. The discontinuous Galerkin (DG) finite element method provides an innovative framework for modeling brittle crack propagation with zero‐thickness interface elements, which can accommodate extrinsic cohesive laws to avoid the artificial compliance required in intrinsic cohesive models. However, robust formulations and implementations of DG methods are critical in alleviating the well‐known convergence issues for both crack nucleation and propagation with reduced instability. This paper presents a robust interface element formulation by modifying the incomplete interior penalty Galerkin (IIPG) method, which successfully avoids the initial element interface penetration across elements that occurs prior to crack nucleation, and thereby greatly reduces the instability issue as cracks open. We further verified and validated our implementation by using a bar tension test and a beam fracturing benchmark. The robustness of our proposed interface element method was demonstrated by a micromechanics fiber/matrix debonding problem with 64 fibers embedded in a bulk matrix.

Liu, Ruijie↗
DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗
DOE OSTI2024

A High-Performance Discrete-Element Framework for Simulating Flow and Jamming of Moisture Bearing Biomass Feedstocks

We developed and verified a high-performance open-source discrete element method (DEM) solver with simultaneously-supported feedstock-specific interaction models, including bonded-sphere, liquid bridge, cohesion, and non-linear contact models. Our solver uses parallel data structures on hybrid central and graphics processing unit (CPU/GPU) architectures, with favorable strong scaling performance observed for large problem sizes comprised of (100 M particles), and 4X single-node GPU speedup. The particles for corn stover feedstock were conceptualized and calibrated based on experimental measurements and results. Sensitivity analyses demonstrate that the mass flow rate from a wedge hopper is governed primarily by moisture content, friction coefficient, and cohesion energy density. The model is used to reproduce experimentally observed hopper jamming results, highlighting that the experimental no-flow trends can only be achieved by using non-spherical particles, liquid bridge and cohesion models, highlighting the importance of using concurrent feedstock specialized models for the effective representation of biomass material handling problems.

bioenergy↗
DOE OSTI2020

Dueling Backbones: Comparing Peptoid and Peptide Analogues of a Mussel Adhesive Protein

Ensembles of amino acid side chains often dominate the interfacial interactions of intrinsically disordered proteins; however, backbone contributions are far from negligible. Using a combination of nanoscale force measurements and molecular dynamics simulations, we demonstrated with analogous mussel-mimetic adhesive peptides and peptoids both 34 residues long that highly divergent adhesive/cohesive outcomes can be achieved on mica surfaces by altering backbone chemistry only. The Phe, Tyr, and Dopa containing peptoid variants used in this study deposited as dehydrated and incompressible films that facilitated analysis of peptoid side chain contributions to adhesion and cohesion. For example, whereas Phe and Dopa peptoids exhibited similar cohesion, Dopa peptoids were ~3 times more adhesive than Phe peptoids on mica. Compared with the peptides, Phe peptoid achieved only ~20% of Phe containing peptide adhesion, but the Dopa peptoids were >2-fold more adhesive than the Dopa peptides. Cation-π interactions accounted for some but not all of the cohesive interactions. Finally, our results were corroborated by molecular dynamics simulations and highlight the importance of backbone chemistry and the potential of peptoids or peptoid/peptide hybrids as wet adhesives and primers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Density functional theory calculation of helium segregation and decohesion effect in W110/Ni111 interphase boundary

Density functional theory (DFT) is employed to study the preferential distribution and decohesion effect of He in a W-NiFe composite consisting of W particles embedded in an Ni-based solid solution matrix. A slab containing {110}<100>W//{111}<110>Ni interface is employed as a surrogate model for the W-NiFe system. Firstly, the fracture energy of the W/Ni interphase boundary (IB) (4.37 J/m 2 ) is higher than the cleavage energy of Ni{111} (3.82 J/m^2) and lower than the cleavage energy of W{110} (6.60 J/m 2 ). The comparison shows that the cohesion of the IB is stronger than the Ni{111} planes that are away from the IB. However, the cohesion between the Ni{111} planes adjacent to the IB is found to be the weakest in this system, with a cleavage energy of 3.11 J/m 2 . Subsequently, formation energy of He is calculated in the Ni slab, W slab, and various interstitial sites in the IB. The DFT calculations show that He is significantly more stable in Ni than in W, by about 1.75 eV. Interestingly, He does not prefer to segregate at the IB as compared to bulk Ni. Nevertheless, it prefers to segregate to the region between the Ni{111} planes adjacent to the IB, and decreases the cohesion of the already weakest region. Based on an estimated amount of He gas production in 5 years under first wall neutron irradiation (neutron flux of 1.04x10 15 n/cm 2 /s), the He will decrease the cleavage energy of the weakest region by 21.2%, 15.4%, and 12.2% at 800, 1000, and 1200 C, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Erosion and cathodic arc plasma of Nb–Al cathodes: composite versus intermetallic

Many properties of cathodic arcs from single-element cathodes show a correlation to the cohesive energy of the cathode material. For example, the burning voltage, the erosion rate, or, to a lesser extent, plasma properties like electron temperatures or average ion energy and charge states. For multi-element cathodes, various phases with different cohesive energies can initially be present in the cathode, or form due to arc exposure, complicating the evaluation of such correlations. To test the influence of morphology and phase composition of multi-element cathodes on cathodic arc properties, a Nb–Al cathode model system was used that includes: pure Nb and Al cathodes; intermetallic Nb 3 Al, Nb 2 Al and NbAl 3 cathodes; and three composite Nb–Al cathodes with atomic ratios corresponding to the stoichiometric ratios of the intermetallic phases. Pulsed cathodic arc plasmas from these cathodes were examined using a mass-per-charge and energy-per-charge analyzer, showing that charge-state-resolved ion energy distributions of plasmas from the intermetallic and corresponding composite cathodes are nearly identical. An examination of converted layers of eroded cathodes using x-ray diffraction and scanning electron microscopy indicates the formation of a surface layer with similar phase composition for intermetallic and their corresponding composite cathode types. The average arc voltages do not follow the trend of cohesive energies of Nb, Al and intermetallic Nb–Al phases, which have been calculated using density functional theory. Possible reasons for this effect are discussed based on the current knowledge of multi-element arc cathodes and their arc plasma available in literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reduction of Ultrafine Hematite for a Spouting Bed Chemical Looping Process - Presentation

Chemical looping combustion research has been a topic of great interest since the mid-1990s, although the concept can be traced much further back. It is a combustion technology for fuel utilization that uses a two-reactor system to produce a carbon free effluent stream from the air reactor and a carbon dioxide rich stream from the fuel reactor with oxygen shuttled between the reactors with an oxygen carrier. Chemical looping combustion has been studied extensively because of this inherent separation of the carbon dioxide combustion product from the atmospheric nitrogen over that past 25 years with hundreds of research papers and tens of books. The oxygen carrier is, in general, a solid such as naturally occurring hematite or an engineered carrier like a copper and iron oxide supported on alumina. System technical and economic research activities conducted at the National Energy Technology Laboratory have shown that the costs associated with the carrier, initial and replacement, are a significant factor in the economic viability of the technology. To overcome this economic challenge, the work reported on in this paper investigates the use of a waste iron oxide as the chemical looping oxygen carrier. This waste material is quite unique as it has a particle size of 0.7 um. Given the particle size, the reactions are quite fast. However, the development of a contactor reactor to mix the fuel with the carrier is a fluidization challenge as the material is a Geldart group C material and highly cohesive. Normally, Geldart group C materials are fluidized with either a vibrating bed. This technique is acceptable for low temperature pharmaceuticals production but does not lend itself to much larger reactors used for power production or operation at temperatures in the order of 1100 to 1200 C. Therefore, a hydrodynamic solution was needed to break up the cohesive structures to replace the vibration. In this instance, it was hypothesized that a spouting bed of large Geldart group B or Geldart group D particles could provide the necessary energy to break up the cohesive structures of the Geldart group C material. This paper presents the reduction performance for this novel carrier material in a laboratory spouting bed reactor.

Breault, Ronald↗
DOE OSTI2022

Transient creep‐fatigue crack growth in creep‐brittle materials: Application to Alloy 718

Abstract This study presents incremental finite element computations of creep‐fatigue crack growth in Alloy 718 at in air. Alloy 718 is representative of creep‐brittle materials, in which viscoplastic deformation is restricted near the tip of a growing crack. The computations predict crack growth using a unique combination of an irreversible cohesive zone formulation and a strain gradient viscoplastic material model based on the Kocks–Mecking formulation. Cohesive zone damage parameters are estimated using sustained loading and constant‐amplitude cyclic loading experiments. Computations of crack extension under three different waveforms containing overloads all predict post‐overload retardation. The amount of retardation depends strongly on the overload ratio, consistent with experiments in the literature using similar waveforms. Analysis of the crack‐tip fields demonstrates retardation is associated with unloading in the highly deformed material near the advancing crack tip. Dynamic recovery and geometrically necessary dislocations are shown to significantly influence post‐overload crack extension.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Geomechanical and Hydrogeological Evaluation of a Shallow Hydraulic Fracture at the Devine Fracture Pilot Site, Medina County, Texas

UT-Austin’s Devine Fracture Pilot Site, 50 miles southwest of San Antonio, Texas, has been targeted for a comprehensive, multidisciplinary development of fracture diagnostic techniques that are cross-validated by ground-truth data acquisition near a recently created, 175-ft-deep, horizontal hydraulic fracture (Ahmadian et al. 2018 Demonstration of proof of concept of electromagnetic geophysical methods for high resolution illumination of induced fracture networks. In Proceedings of the SPE Hydraulic Fracturing Technology Conference and Exhibition, The Woodlands, Texas, USA, 23–25 January 2018. SPE-189858-MS.). To evaluate the fracture diagnostic methods at this site, we conducted injection tests with a predefined volumetric flow-rate profile, resembling a diagnostic fracture injection test on September 2020. Subsequently, we developed hydrogeological and geomechanical models based on flow-rate and bottomhole-pressure measurements. History-matching efforts using a simplified layer-cake hydrogeological model resulted in the field-scale formation permeability of 9.87 × 10 –15 m 2 (10 mD) and Darcy-scale fracture permeability. The analysis of the bottomhole pressure and injection-rate history showed that (1) the newly created horizontal fracture was closed adjacent to the injection well pre-injection and (2) the initial pump-pressure increase at a nominal volumetric injection rate led to near-well fracture reopening, fluid conductivity increase, and abrupt injection-rate increase. To overcome hydrogeological-model limitations of predicting fracture reopening throughout injection, we extended the modeling to a finite-element, poroelastic analysis of horizontal-fracture growth using a cohesive-zone model. Using this fracture-reopening model, we improved the history match of the transient-pressure response during the experiment by adjusting the hydromechanical properties. Furthermore, post-injection pressure transient analyses helped reduce uncertainty in the overburden-stress gradient, and the initial hydraulic-fracturing simulation verified the plausibility of achieving the surveyed propped fracture area.

02 PETROLEUM↗
DOE OSTI2022

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deformation and Delamination of Dynamically Bulging Bilayered Films

We employ finite element simulations to study the dynamic bulging of films under shock tube conditions and compare these responses with quasi-static bulging. Both scenarios subject the film to time-varying pressure. Two film types are considered: uniform material films and bilayers with a relatively soft interior layer and a hard exterior layer. The interface between the materials in the bilayer cases is modeled as a cohesive zone that follows a bilinear normal traction-separation cohesive law. Additionally, bilayer films with preexisting cracks in the interior layer are analyzed. Parametric studies are also conducted to explore the effects of loading rate and specimen size. The mechanical response of all materials is assumed to be strain-rate independent, focusing solely on the inertial effects in the response of the films under dynamic loading, and the absence of such effects under quasi-static loading. The results indicate that quasi-static bulging exhibits a predictable plate bending like deformation. In contrast, depending on the loading rate and specimen size, dynamic bulging can give rise to elastic waves, and a sequence of deformation processes including initial uniform acceleration, followed by bending-like deformation, and finally through-thickness biaxial stretching. Dynamic loading also leads to larger bulge formation and greater deformation than quasi-static loading. A key finding is the role of preexisting cracks in driving delamination under dynamic loading, which is absent in quasi-static loading. This work enhances our understanding of dynamic bulging and motivates further research on characterizing delamination in layered films using dynamic bulge tests.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modeling the interfacial failure and resulting mechanical properties of z-pinned additively manufactured composites

The mechanical properties of Fused Filament Fabrication (FFF) parts are limited in the build direction, albeit advantages such as design flexibility and in-house customization. Here, in this paper, a novel z-pinning approach that enhances mechanical properties in the build direction by depositing material across multiple layers within the part was investigated through numerical simulations and validated experimentally. A finite element model for z-pinned composite structures was developed by assigning fiber-orientation dependent material properties obtained using a micromechanics approach to beads and pins. The bead-to-bead and pin-to-pin adhesions within the z-pinned structures were modeled using a cohesive traction separation law. The properties of cohesive elements for carbon fiber-reinforced polylactic acid (CF-PLA) z-pinned composites were calibrated using tensile experiments. The elastic modulus and tensile strength of CF-PLA z-pinned composites in the build direction were predicted with the developed numerical model. The numerical investigation on various geometrical parameters revealed that the largest pin volume increases the stiffness and tensile strength by 40% and thus, has the greatest influence on the mechanical properties of the z-pinned composites. The effect of z-pin geometrical parameters on the mechanical properties was summarized to aid in the design of z-pinned additively manufactured composite structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗
DOE OSTI2020

Perspectives on intrinsic toughening strategies and passivation of perovskite films with organic additives

We report on the inclusion of 5-aminovaleric acid (5-AVA)—a bulky, long-chained organic molecule—as an additive to enhance the mechanical integrity of hybrid perovskite films. We studied a range of organic cation additives and found that 5-AVA was the most effective at improving cohesion energy, a key metric of mechanical reliability. MAPbI 3 films reinforced with 5-AVA added in 5% concentration increased cohesion energy 12-fold from 0.53 J/m 2 to 6.04 J/m 2 , an effect which is attributed to increased plasticity and crack deflection around grain boundaries of the additive-containing perovskite. The addition of 5-AVA also improves V oc in perovskite solar cells and carrier lifetimes with a minimal decrease in PCE, attributed to passivation reducing defect and trap densities. A Tauc plot analysis of the bandgap shows that 5-AVA increases the band gap of the perovskite, correlating with reduced film stress compared to MAPbI 3 . Finally, the usage of 5-AVA improves the intrinsic thermomechanical reliability, but this improvement comes at the penalty of slightly reduced device efficiency due to reduced charge extraction from the presence of the bulky, insulating organic additive.

14 SOLAR ENERGY↗