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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulations of crack extensions in arc-shaped tension specimens of uncharged and hydrogen-charged 21-6-9 austenitic stainless steels using cohesive zone modeling with varying cohesive parameters

Crack extensions in uncharged and hydrogen-charged side-grooved A(T) specimens of conventionally forged 21-6-9 austenitic stainless steels are simulated using the cohesive zone modeling approach. Two-dimensional plane strain finite element analyses with fixed cohesive parameters are first conducted to fit the experimental load-displacement curves. Similar analyses using varying cohesive parameters as functions of the crack extension are then conducted to fit the experimental load-crack extension and crack extension-displacement curves. The computational results with varying cohesive parameters can fit very well the experimental data. The computational results also indicate that the average cohesive energy for the hydrogen-charged A(T) specimen is lower than that for the uncharged A(T) specimen.

42 ENGINEERING↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Large-scale protein level comparison of Deltaproteobacteria reveals cohesive metabolic groups

Abstract Deltaproteobacteria, now proposed to be the phyla Desulfobacterota, Myxococcota, and SAR324, are ubiquitous in marine environments and play essential roles in global carbon, sulfur, and nutrient cycling. Despite their importance, our understanding of these bacteria is biased towards cultured organisms. Here we address this gap by compiling a genomic catalog of 1 792 genomes, including 402 newly reconstructed and characterized metagenome-assembled genomes (MAGs) from coastal and deep-sea sediments. Phylogenomic analyses reveal that many of these novel MAGs are uncultured representatives of Myxococcota and Desulfobacterota that are understudied. To better characterize Deltaproteobacteria diversity, metabolism, and ecology, we clustered ~1 500 genomes based on the presence/absence patterns of their protein families. Protein content analysis coupled with large-scale metabolic reconstructions separates eight genomic clusters of Deltaproteobacteria with unique metabolic profiles. While these eight clusters largely correspond to phylogeny, there are exceptions where more distantly related organisms appear to have similar ecological roles and closely related organisms have distinct protein content. Our analyses have identified previously unrecognized roles in the cycling of methylamines and denitrification among uncultured Deltaproteobacteria. This new view of Deltaproteobacteria diversity expands our understanding of these dominant bacteria and highlights metabolic abilities across diverse taxa.

Langwig, Marguerite V. (ORCID:0000000202472816)↗
NASA NTRS2011

Investigating Some Technical Issues on Cohesive Zone Modeling of Fracture

This study investigates some technical issues related to the use of cohesive zone models (CZMs) in modeling fracture processes. These issues include: why cohesive laws of different shapes can produce similar fracture predictions; under what conditions CZM predictions have a high degree of agreement with linear elastic fracture mechanics (LEFM) analysis results; when the shape of cohesive laws becomes important in the fracture predictions; and why the opening profile along the cohesive zone length needs to be accurately predicted. Two cohesive models were used in this study to address these technical issues. They are the linear softening cohesive model and the Dugdale perfectly plastic cohesive model. Each cohesive model constitutes five cohesive laws of different maximum tractions. All cohesive laws have the same cohesive work rate (CWR) which is defined by the area under the traction-separation curve. The effects of the maximum traction on the cohesive zone length and the critical remote applied stress are investigated for both models. For a CZM to predict a fracture load similar to that obtained by an LEFM analysis, the cohesive zone length needs to be much smaller than the crack length, which reflects the small scale yielding condition requirement for LEFM analysis to be valid. For large-scale cohesive zone cases, the predicted critical remote applied stresses depend on the shape of cohesive models used and can significantly deviate from LEFM results. Furthermore, this study also reveals the importance of accurately predicting the cohesive zone profile in determining the critical remote applied load.

Wang, John T.↗
NASA NTRSDate not supplied

The Interplay of Humor and Team Cohesion Over Time

Team bonding and interpersonal relationships are increasingly important as we continue understanding the effects of the pandemic on team dynamics. Specifically, team cohesion, defined as having a sense of belonging and commitment towards the group, has been found to be important for key team outcomes such as performance, relationships, and productivity (Beal, Cohen, Burke, & McLendon, 2003; Bell & Marentette, 2011). Research has stressed the benefits of positive humor on cohesion, with humor creating positive emotions and decreasing social distance (e.g., Robert & Wilbanks, 2012; Robinson & Smith-Lovin, 2001; Romero & Pescosolido, 2008). However, research around the impact of negative humor on cohesion is scarce. This study aims to examine the impact of both positive and negative humor on cohesion for teams over time. The longitudinal design allows us to examine the interplay of humor and cohesion throughout a team’s lifecycle to guide relationship development and management (Romero & Pescosolido, 2008). Humor is defined as remarks recognized as amusing or funny (Martin, 2007). Positive humor can be categorized as good-natured aimed to create mutual agreement, while negative humor is defined as aggressive, hostile, and unreciprocated (Coan & Gottman, 2007). Research has concluded the many benefits of positive humor, showing it can increase open communication, social support, and create a pleasant environment (Martineau, 1972; Robinson & Smith-Lovin, 2001; Romero & Pescosolido, 2008). The Humor Wheel model suggests that positive humor is contagious due to the positive affect produced and mimicked (Lehmann-Willenbrock & Allen 2014). Research has suggested positive humor would increase over time based on the norm of reciprocity; those treated positively are likely to respond similarly (Gouldner, 1960). Multiple parties’ involvement in humor exchanges highlights the interpersonal nature of humor. Thus, long-term consequences on relationships should be considered, such as cohesion (Robert & Wilbanks, 2012). It would be instrumental to understand how humor interactions could weaken or strengthen group cohesion through relational and emotional bonds over time (Martineau, 1972; Romero & Pescosolido, 2008). Positive humor should fundamentally increase team cohesion (Mesmer-Magnus, Glew, & Viswesvaran, 2012; Walter & Baruch, 2008) by promoting group harmony, team member attractiveness, and ingroup identity (Holmes, 2006). Thus, it is expected there is a positive relationship between positive humor and team cohesion over time (hypothesis 1). Humor is not always beneficial, especially when it elicits negative emotions (Dikkers, Doosje, & de Lange, 2011). Negative humor can be used to humiliate, tease, or insult (Coan & Gottman, 2007; Collinson, 1988; Romero & Cruthirds, 2006). It tends to cause distress and damage relationships, and in turn, decrease organizational performance and production (Avolio, Howell, & Sosik, 1999; Janes & Olson, 2000; Yerkes, 2001). Research concludes negative humor may create subgroups and decrease cooperation (Romero & Arendt, 2011). Research has also indicated positive implications for negative humor under specific circumstances. Negative humor could help to solidify ingroup identity and strengthen bonds when directed outside the team (Terrion & Ashforth, 2002; Thomae & Pina, 2015). For example, sport teams may use negative humor towards their opponent before a game, which emphasizes their differences and lowers the status and power of their opponent, ultimately promoting solidarity and cohesion inside the group. For this study, it is anticipated that negative humor toward other teammates will decrease team cohesion over time due to the damage to relationships and negative affect produced (Duncan, Smeltzer, & Leap, 1990; Maples, Dupey, Torres‐Rivera, Phan, Vereen, & Garrett, 2001). We also predicted that negative humor focused on something or someone outside the team will increase cohesion by enhancing ingroup ties and creating a common identity (Thomae & Pina, 2015). Thus, the relationship between negative humor and cohesion will be moderated by the target of the humor. The relationship will be negative when the target is within the team and positive when the target is outside (hypothesis 2).

Melissa M. Vazquez↗
DOE OSTI2024

Intrinsic and environmental drivers of pairwise cohesion in wild Canis social groups

Animals within social groups respond to costs and benefits of sociality by adjusting the proportion of time they spend in close proximity to other individuals in the group (cohesion). Variation in cohesion between individuals, in turn, shapes important group-level processes such as subgroup formation and fission–fusion dynamics. Although critical to animal sociality, a comprehensive understanding of the factors influencing cohesion remains a gap in our knowledge of cooperative behavior in animals. We tracked 574 individuals from six species within the genus Canis in 15 countries on four continents with GPS telemetry to estimate the time that pairs of individuals within social groups spent in close proximity and test hypotheses regarding drivers of cohesion. Pairs of social canids (Canis spp.) varied widely in the proportion of time they spent together (5%–100%) during seasonal monitoring periods relative to both intrinsic characteristics and environmental conditions. The majority of our data came from three species of wolves (gray wolves, eastern wolves, and red wolves) and coyotes. For these species, cohesion within social groups was greatest between breeding pairs and varied seasonally as the nature of cooperative activities changed relative to annual life history patterns. Across species, wolves were more cohesive than coyotes. For wolves, pairs were less cohesive in larger groups, and when suitable, small prey was present reflecting the constraints of food resources and intragroup competition on social associations. Pair cohesion in wolves declined with increased anthropogenic modification of the landscape and greater climatic variability, underscoring challenges for conserving social top predators in a changing world. We show that pairwise cohesion in social groups varies strongly both within and across Canis species, as individuals respond to changing ecological context defined by resources, competition, and anthropogenic disturbance. Our work highlights that cohesion is a highly plastic component of animal sociality that holds significant promise for elucidating ecological and evolutionary mechanisms underlying cooperative behavior.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗
DOE OSTI2022

Calibration of cohesive parameters for fracture test simulations of sub-sized bend specimens of AM steels

In this work, crack extensions in sub-sized side-grooved bend specimens of AM steels are examined with changing cohesive parameters. Two-dimensional FE analyses are first conducted to simulate crack extensions by nodal release. The simulation results indicate the near-tip maximum opening stress and separation work rate increase and then slightly decrease for the bend specimens. The hydrogen-charged bend specimens have lower near-tip maximum opening stresses and average separation work rates than those of the uncharged bend specimens. The changing near-tip maximum opening stresses and the separation work rates with increasing crack extension are determined as references for the changing cohesive strengths and energies, respectively. Two-dimensional FE analyses with different changing cohesive parameters are then performed. The simulation results with calibrated cohesive parameters can match well with the test results. The hydrogen-charged specimens have the lower changing cohesive strength vs crack growth curves compared with those of the uncharged specimens. The hydrogen-charged specimens have either lower or similar changing cohesive energy vs crack growth curves compared with those of the uncharged specimens. Finally, a three-dimensional FE analysis is conducted to simulate crack extension by nodal release in one bend specimen to examine the stress triaxialities in front of the growing crack fronts at different locations from the center plane to the planes near the side groove front with increasing crack extension to identify possible physical reasons for the changing cohesive parameters.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Relating Cohesive Zone Model to Linear Elastic Fracture Mechanics

The conditions required for a cohesive zone model (CZM) to predict a failure load of a cracked structure similar to that obtained by a linear elastic fracture mechanics (LEFM) analysis are investigated in this paper. This study clarifies why many different phenomenological cohesive laws can produce similar fracture predictions. Analytical results for five cohesive zone models are obtained, using five different cohesive laws that have the same cohesive work rate (CWR-area under the traction-separation curve) but different maximum tractions. The effect of the maximum traction on the predicted cohesive zone length and the remote applied load at fracture is presented. Similar to the small scale yielding condition for an LEFM analysis to be valid. the cohesive zone length also needs to be much smaller than the crack length. This is a necessary condition for a CZM to obtain a fracture prediction equivalent to an LEFM result.

Wang, John T.↗
DOE OSTI2024

A Modeling framework for flocculated cohesive sediment transport in the current bottom boundary layer

Cohesive sediment transport, where its settling velocity is controlled by the flocculation process, is a crucial component in determining biochemical cycles, fate of pollutants, and morphodynamics in many aquatic ecosystems. In this study, a modeling framework is presented to investigate how flocculation influences cohesive sediment transport in the current bottom boundary layer in dilute conditions, consistent with the calibration range of the flocculation model. From a local analysis of floc dynamics in homogenous turbulence, we identify that the floc size distribution is mainly controlled by floc cohesion and yield strength. The uncertainty in fractal dimension plays a minor role for the floc size but it influences the resulting floc density and settling velocity. The transport analysis in the current boundary layer shows that the flocculation process alters the vertical distribution of the settling velocity and hence the sediment concentration with a strong dependence on cohesion, floc yield strength, and floc structure. When the flocs are more susceptible to breaking, a well-mixed concentration profile is obtained. In contrast, for flocs with higher cohesion or yield strength, higher concentration with a sharp gradient is observed close to the bed. Overall, the settling velocity exhibits a low vertical variability within 20 % of the depth-averaged value except near the bed. Further, this suggests that using a depth-averaged settling velocity yields acceptable predictions of the sediment concentration profiles, especially for flocs with lower cohesion.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

From S-N to the Paris Law with a New Mixed-Mode Cohesive Fatigue Model

The relationship between fatigue life and fatigue crack propagation rate is explored with a new cohesive damage model. The parameters of the model are obtained from idealizations of S-N diagrams used in engineering design. The model is based on the hypothesis that both stable tearing damage and damage due to cyclic loading are representations of a density of microcracks and, therefore, a single damage variable can describe the state of damage. This assumption implies that the quasi-static cohesive law that describes tearing is also the envelope of the fatigue damage. Fatigue damage within the cohesive envelope is assumed to accumulate at a rate that depends on the displacement jumps. The fatigue model was implemented as a UMAT subroutine for Abaqus cohesive elements by adding fatigue damage accumulation to a cohesive model based on the Turon mixed-mode cohesive laws. The analyses were conducted using a simplified cyclic loading procedure in which the maximum applied load is kept constant and the computational expense of cycling the load is avoided. The predicted propagation rates in mode I and mixed mode were compared to experimental results for IM7/8552 graphite/epoxy tape. Several aspects of the results were investigated, including the effect of R-curves, the stress ratio R, and the difference between displacement control and force control.

Davila, Carlos G.↗
NASA NTRS2018

A Cohesive Fatigue Model Based on the S-N Diagram

The relationship between fatigue life and fatigue crack propagation rate is explored with a new cohesive damage model. The parameters of the model are obtained from idealizations of S-N diagrams used in engineering design. The model is based on the hypothesis that both stable tearing damage and damage due to cyclic loading are representations of a density of microcracks and, therefore, a single damage variable can describe the state of damage. This assumption implies that the quasi-static cohesive law that describes tearing is also the envelope of the fatigue damage. Fatigue damage within the cohesive envelope is assumed to accumulate at a rate that depends on the displacement jumps. The fatigue model was implemented as a UMAT subroutine for Abaqus cohesive elements by adding fatigue damage accumulation to a cohesive model based on the Turon quasi-static cohesive law. The analyses were conducted using a simplified cyclic loading procedure in which the maximum applied load is kept constant and the computational expense of cycling the load is avoided. The predicted propagation rates for a double cantilever beam (DCB) specimen were compared to experimental results for IM7/8552 graphite/epoxy tape.

Dávila, Carlos G.↗
DOE OSTI2020

A Custom Machine to Convey Radiologically Impacted Cohesive Soils to the Orion ScanSortSM System - 20382

An established method to reduce the volume of material required to be disposed of as radioactive waste is conveyor-based sorting of potentially radiologically impacted soils. With the Orion ScanSortSM system, potentially impacted soils are conveyed beneath radiation detectors and sorted into above and below criteria bins based on the detectors' response. Confident measurements can only be efficiently achieved when the soil column being conveyed has a reasonably consistent geometry. Highly cohesive soils present a very significant challenge as they tend to clump and adhere to the conveyor hopper, strike-off bar, and belt skirting. This causes voids and valleys in the soil column that affect the measurement geometry and reduce the confidence in the measurement. The reduction in confidence necessitates a longer residence time (i.e. reduction in conveyor belt speed and processing rates) or that the soils in highly unfavorable geometries be dispositioned as impacted or segregated for resurvey (should site logistics support that option). These outcomes may increase the volume of material required to be dispositioned as radiological waste and have a negative impact on project cost and schedule. To minimize impact to sorting operations, Wood designed and built a customized Extruder, to optimize the soil column prior to assay. The design was developed and refined to minimize the volume of highly cohesive soils presented in highly unfavorable geometries in the soil column in order to confidently assess and disposition such soils in a high-production environment. The Extruder is a customized 90 cm wide flat conveyor that has an oversized hopper with a pair of motor-driven rollers that force the soils in the hopper through an opening of adjustable height. The soils are extruded through the opening to produce a soil column that is ∼75 cm wide and ∼8-18 cm deep with a design rate of 225 metric tons/hr. The Extruder was recently deployed with the Orion ScanSort{sup SM} System to a site in northeast Ohio. The site has highly cohesive soils that result from the weathering of glacial sediments and consist of clay and clay loam, resulting in poor drainage and high moisture content. The project infrastructure only supported two bins of material (above and below criteria), which necessitated that soils with unfavorable geometries be discharged into the above criteria bin, along with material determined to have been contaminated. The Extruder was extremely successful in producing a stable soil column with favorable geometries for radiological assay. During production, less than 0.5% (by mass) of the soils processed were presented in unsatisfactory geometries requiring disposition as radioactive waste. The Extruder did not get clogged due to the cohesive soils (as is typical with conventional conveyors) and required very little maintenance, again minimizing impacts to sorting operations. The conveyor belt speed was operated at 14 cm/s with a typical belt loading of ∼109 kg/m resulting in a mean process rate of ∼49 metric tons/hr. The process rate was constrained by other project logistics, rather than by the capabilities of the Extruder, which was operated at less than 20% of the maximum design speed (75 cm/s). It is concluded that the Extruder is highly likely to produce a stable highly cohesive soil column suitable for efficient and confident radiological assay in production environments of 250 short tons/hr or more. The Extruder has the capability to drastically reduce the duration and costs of projects with radiologically impacted highly cohesive soils. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
DOE OSTI2024

Unified cohesive zone model (UCZM) for fracturing and fragmenting solids

Here, a Unified Cohesive Zone Model (UCZM), which inherits most of the advantages while overcoming the shortcomings of existing Cohesive Zone Models (CZMs), is proposed. Similar to the traditional extrinsic CZM approach, UCZM dynamically inserts the cohesive elements into the system based on local material states (e.g., stress, strain). However, the transition from continua to discontinua is smoothly achieved, thereby eliminating the “time-discontinuous” issue seen in the extrinsic CZM. Moreover, within the novel UCZM framework, the point of transition from continua to discontinua is controllable through the introduction of crack initialization criteria. As a result, the UCZM allows any material models (e.g., elastic, plastic, damage models) for continuum solids and for discrete fracture behavior to work together. In essence, both an enhanced extrinsic cohesive zone model and an intrinsic cohesive zone model can be represented by the proposed unified model. The proposed UCZM has been verified through different numerical examples. The work demonstrates that the UCZM is a highly effective approach for modeling fracture and fragmentation processes in solids.

42 ENGINEERING↗
DOE OSTI2020

Multiphysics for nuclear energy applications using a cohesive computational framework

With the recent development of advanced numerical algorithms, software design, and low-cost high-performance computer hardware, reliance on coupled multiphysics to predict the behavior of complex physical systems is beginning to become standard practice. This is especially true in nuclear energy applications where strong nonlinear interdependencies exist between reactor physics, radiation transport, multi-scale nuclear fuels performance, thermal fluids, etc. Resolving these nonlinear dependencies requires choices in multiphysics software approaches. Two main multiphysics modeling and simulation approaches have emerged. The first is based upon "code coupling" where disparate physics codes of different software design, code languages, and spatial and temporal integration schemes are coupled together with relatively complex data passing interfaces. The second multiphysics software approach is to employ a "cohesive" framework where all physics applications are developed with a common software design, i.e., data structures, syntax, input format, integrated spatial and temporal discretization schemes, etc. In this paper we present the Multiphysics Object-Oriented Simulation Environment (MOOSE) development and runtime framework and describe the framework's cohesive modeling and simulation multiphysics approach. Then, a "cohesive-like" extension of the MOOSE framework is presented where MOOSE-based physics software applications are efficiently coupled to non-MOOSE (external) physics codes to form multiphysics applications using MOOSE's unique interface capabilities. Finally, several examples of MOOSE's cohesive and cohesive-like multiphysics applications will be demonstrated. These multiphysics demonstrations will incorporate both MOOSE-based applications and external codes, including Nek5000, RELAP-7, TRACE, BISON, and Pronghorn.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2006

Molecular-dynamics Simulation-based Cohesive Zone Representation of Intergranular Fracture Processes in Aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics simulations. A molecular-dynamics model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat 99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element an atomistic analog to a continuum cohesive zone model element - the results from the molecular-dynamics simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion. Keywords: Crack-tip plasticity; Cohesive zone model; Grain boundary decohesion; Intergranular fracture; Molecular-dynamics simulation

Yamakov, Vesselin I.↗