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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Visualizing pectin polymer-polymer entanglement produced by interfacial water movement

In this report, we investigated the physical conditions for creating pectin polymer-polymer (homopolymer) entanglement. The potential role of water movement in creating pectin entanglement was investigated by placing water droplets-equivalent to the water content of two gel phase films-between two glass phase films and compressing the films at variable probe velocities. Slow probe velocity (0.5 mm/sec) demonstrated no significant debonding. Furthermore, corresponding videomicroscopy demonstrated an occasional water bridge, but no evidence of stranding or polymer entanglement. In contrast, fast probe velocity (5 mm/sec) resulted in 1) an increase in peak adhesion strength, 2) a progressive debonding curve, and 3) increased work of cohesion (p < .001). Corresponding videomicroscopy demonstrated pectin stranding and delamination between pectin films. Scanning electron microscopy images obtained during pectin debonding provided additional evidence of both stranding and delamination. We conclude that water movement can supply the motive force for the rapid chain entanglement between pectin films.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Simulation of multi-shell fullerenes using Machine-Learning Gaussian Approximation Potential

Multi-shell fullerenes ”buckyonions ” were simulated, starting from initially random configurations, using a density-functional-theory (DFT)-trained machine-learning carbon potential within the Gaussian Approximation Potential (GAP) Framework [Volker L. Deringer and Gábor Csányi, Phys. Rev. B 95, 094203 (2017)]. Fullerenes formed from seven different system sizes, ranging from 60 ~ 3774 atoms, were considered. The buckyonions are formed by clustering and layering starting from the outermost shell and proceeding inward. Inter-shell cohesion is partly due to interaction between delocalized π electrons protruding into the gallery. The energies of the models were validated ex post facto using density functional codes, VASP and SIESTA , revealing an energy difference within the range of 0.02 - 0.08 eV/atom after conjugate gradient energy convergence of the models was achieved with both methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A new approach to model geomaterials with heterogeneous properties in thermo-hydro-mechanical coupled problems

The main objective of this article is to present a new approach to model coupled thermo-hydro-mechanical problems considering geomaterials with heterogeneous properties. This approach has been implemented in the software CODE_BRIGHT and it provides the possibility of considering geomaterials with a spatially correlated heterogeneous field of porosity, following a normal distribution. This spatial correlation can be isotropic or anisotropic. An important feature of this approach is that material properties such as intrinsic permeability, thermal conductivity, diffusivity, retention curve, elastic modulus or cohesion are defined as a function of porosity and, thus, they become heterogeneous with spatial correlation and, eventually, anisotropic. A validation exercise and other basic numerical examples have been carried out to illustrate the possibilities of the proposed approach. The results, which have been compared with a homogeneous case, show that considering heterogeneous fields can be relevant in different modelling problems, especially coupled thermo-hydro-mechanical problems.

42 ENGINEERING↗
DOE OSTI2025

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗
DOE OSTI2025

Thermal and mechanical performance of gypsum composites containing fly ash cenosphere encapsulated phase change materials

Utilizing phase change materials (PCMs) within construction materials improves the thermal energy storage capabilities in the built environment. Notably, integrating PCMs into drywalls presents a promising avenue to reduce building energy use. However, incorporating PCMs into drywalls poses challenges, particularly reduced mechanical properties and workability limitations as it can adversely affect the hardening and cohesiveness of gypsum composites. To address these challenges, this study explores a novel strategy to incorporate PCMs into gypsum composites through fly ash cenosphere-encapsulated PCM, namely CenoPCM, aiming to enhance the mechanical properties and thermal conductivity over existing commercial counterparts. An experimental program was conducted to study the impact of the inclusion of CenoPCM on the mechanical and thermal performance of gypsum composites, and the results were compared to those of commercially available polymerencapsulated PCM benchmark (Micronal). The findings highlighted that PCM-charged samples exhibited up to 95.90 J/g (CenoPCM) of latent heat with 70 % CenoPCM loading. The results of mechanical tests also showed that CenoPCM outperformed its polymer-encapsulated counterpart, achieving 55 % higher mechanical strength on average. Here, this improvement is attributed to the rigid shell of fly ash cenosphere (FAC) serving as a PCM carrier, providing a skeleton-like structure to gypsum composites for higher load-bearing capacity. Additionally, hotbox tests demonstrated that the CenoPCM-charged wall panels yielded a 20 % higher thermal conductivity and a 22 % increased time lag compared to Micronal PCM-charged panels, indicating its potential to shift the peak load in buildings.

Ceno PCM↗
DOE OSTI2022

Leveraging SETS resilience capabilities for safe-to-fail infrastructure under climate change

As the rehabilitation of infrastructure is outpaced by changes in the profile, frequency, and intensity of extreme weather events, infrastructure’s service disruptions and failures become increasingly likely. Safe-to-fail approaches for infrastructure planning and design improve the capacity of cities to adapt for uncertain climate futures by identifying social, ecological, and technological systems (SETS) capabilities to prepare for potential failure scenarios. In this work, we argue for transforming infrastructure planning and design to effectively utilize safe-to-fail approaches by navigating the opportunities and trade-offs of SETS resilience capabilities. From a technological vantage point, traditional infrastructure planning approaches account for social and ecological domains as external design conditions rather than embedded system characteristics. Safe-to-fail approaches directly challenge the isolation of the technological domain by necessitating a recognition that SETS domains are interconnected and interdependent in infrastructure systems, as such risks and system capabilities for resilience must be managed cohesively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

LAVA 1.0: A general-purpose python toolkit for calculation of material properties with LAMMPS and VASP

Here, we introduce LAVA, a general-purpose python toolkit to provide user-friendly and high throughput calculations for material properties using both LAMMPS and VASP. It contains a set of classes as well as pre-processing and post-processing functions to prepare, execute and extract information from LAMMPS/VASP simulations. An overview of the program structure is provided. The current version contains modules to calculate elastic and mechanical properties such as lattice constant, cohesive energy, cold curve, elastic constants, bulk and shear modulus, volume conserving/non-conserving deformation path, vacancy/interstitial formation energy, surface energy, stacking fault energy, melting point, radial distribution function, and thermal expansion. These functionalities are demonstrated for different interatomic potentials and DFT calculations in Al.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ensemble Monte Carlo calculations with five novel moves

We introduce five novel types of Monte Carlo (MC) moves that brings the number of moves of ensemble MC calculations from three to eight. So far such calculations have relied on affine invariant stretch moves that were originally introduced by Christen (2007), walk moves by Goodman and Weare (2010) and quadratic moves by Militzer (2023). Ensemble MC methods have been very popular because they harness information about the fitness landscape from a population of walkers rather than relying on expert knowledge. Here we modified the affine method and employed a simplex of points to set the stretch direction. We adopt the simplex concept to quadratic moves. We also generalize quadratic moves to arbitrary order. Finally, we introduce directed moves that employ the values of the probability density while all other types of moves rely solely on the location of the walkers. We apply all algorithms to the Rosenbrock density in 2 and 20 dimensions and to the ring potential in 12 and 24 dimensions. We evaluate their efficiency by comparing error bars, autocorrelation time, travel time, and the level of cohesion that measures whether any walkers were left behind. Our code is open source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A renewable lignin-based thermoplastic adhesive for steel joining

Adhesive bonding of metals has become increasingly relevant in recent years due to the demand for reducing weight and improving performance in structural applications such as automobiles and aerospace. We developed renewable thermoplastic adhesives from technical organosolv lignin isolated from hardwood biomass and acrylonitrile butadiene co-polymer rubber (NBR) for joining steel substrates. NBR33, NBR41 and NBR51 with acrylonitrile molar ratios of 33, 41 and 51%, respectively, were blended with lignin to form two-phase thermoplastic adhesives, and their adhesion, viscoelastic and surface properties were measured. Lignin content in the compositions were varied, ranging from 40% to 80% (w/w), to alter toughness, stiffness, and surface energy characteristics of the material. Better interaction or reactivity between the lignin and NBR phases was observed with greater nitrile content in NBR, leading to greater modulus and stiffness of the adhesive. Simultaneously, increasing the proportion of lignin reduced toughness and improved stiffness, with the highest adhesive strength of 13.1 MPa measured in a 60% lignin loading ratio with NBR51. Surface energy measurements revealed that total surface energy (sum of polar and dispersive surface energy) raised with lignin loading, suggesting that both surface energy and matrix strength play a critical role in the adhesive properties of the synthesized materials. A finite element-based cohesive zone model (CZM) was developed and implemented to study the failure strength of the adhesively bonded joint. Here, this study demonstrates the viability of lignin as a valuable building block for adhesives, not only due to its inherent chemical structure and rigidity, but also for its surface energy characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Design of biphenylene-derived tunable dirac materials

The exploration of carbon allotropes has unveiled a series of two-dimensional (2D) materials with unique electronic and mechanical properties, yet the need for stable structures with tailored electronic properties persists. Here, in this study, we introduce a new class of 2D carbon allotropes derived from the biphenylene network (BPN), incorporating acetylenic linkages to tune their structural and electronic characteristics. Through density functional theory calculations, we identified ten novel BPN-derived structures that exhibit both energetic and dynamic stability, confirmed by cohesive energy and phonon spectrum analyses. Among them, BPN-02 and BPN04 are metallic, featuring critically-tilted type-III Dirac cones under ~ 5 % biaxial strain, while BPN-22 is a semiconductor with a band gap of 0.95 eV and exhibits highly anisotropic carrier mobility. Additionally, these structures demonstrate significant anisotropy in their elastic properties, further distinguishing them from other 2D carbon materials like graphene. Our findings suggest that these novel BPN-based structures have strong potential for next-generation electronic and optoelectronic applications, providing new avenues for the design and synthesis of advanced carbon materials.

2D carbon allotropes↗
DOE OSTI2026

Third-body stabilization of supercritical CO 2 in CO oxidation: development and application of a ReaxFF force field for the CO/O/CO 2 system

Supercritical CO 2 (scCO 2 ) plays a crucial role as a solvent in separation processes, advanced power cycles, and materials processing. Nonetheless, the atomistic comprehension of how the dense scCO 2 matrix influences the fundamental reaction of carbon monoxide (CO) is still insufficiently explored. Experimental studies and molecular dynamics (MD) simulations frequently fail to detect the highly reactive, transient intermediates, such as atomic oxygen (O), that drive these reactions. Here, to address this issue, we have developed a novel ReaxFF reactive force field for the CO 2 /CO/O system. The force field parameters were calibrated using density functional theory and second-order Møller-Plesset calculations to model CO 2 crystal properties, intermolecular interactions, bond dissociation curves, and reaction energy barriers. The force field reproduces the cohesive energy of the CO 2 crystal, the pressure characteristics of bulk scCO 2 , the equation-of-state behavior over a wide pressure–density range, the pressure dependence of the C–O bond length under compression, and the structural properties of liquid and scCO 2 , as documented by experiments, ab-initio MD, and prominent non-reactive models. The force field was subsequently applied to study the CO + O → CO 2 reaction. In a dilute environment, the reaction is inefficient as the newly formed CO 2 rapidly dissociates due to excess kinetic and potential energy acquired from the exothermic reaction. Conversely, in a dense scCO 2 environment, the surrounding matrix acts as an efficient third body, stabilizing the emerging CO2 product via molecular collisions. Statistical analysis confirms an average excess energy dissipation of 133.9 ± 3.6 kcal/mol over 112.4 ± 17.9 ps. Kinetic energy decomposition reveals that ∼ 92% of the excess kinetic energy is stored in internal (rotational and vibrational) degrees of freedom. This ReaxFF force field establishes a mechanistic foundation for third-body stabilization in dense reactive environments.

Chowdhury, Emdadul Haque [Pennsylvania State Univ.↗
DOE OSTI2026

Geomechanical properties of the Meletta sandstone - the high-temperature heat storage reservoir rock of DeepStor

The DeepStor project aims at storing excess heat at temperatures up to 140 °C in the depleted Leopoldshafen oil field at a depth of about 1300 m. In order to gain knowledge on the target horizons, the different layers of the Meletta sandstone, samples cored in a block retrieved in a quarry near Nussloch where the Meletta sandstone outcrops were studied. Several petrophysical properties were investigated including mineralogy, porosity, permeability, thermal conductivity, P- and S-wave velocities and quality factor. A mechanical study focused on compressive strength under uniaxial and triaxial stress conditions, tensile strength, critical pressure and stress-dependence of physical properties. Our results show that the Meletta sandstone is heterogeneous, anisotropic, mechanically weak, stress-sensitive and prone to water weakening. Petrophysical measurements on few Meletta sandstone cores retrieved in boreholes at about 1250 m depth showed that the outcrop samples are significantly less cohesive and weaker. Therefore, one should consider the Nussloch outcrop as a good reservoir analog only with great caution. Several models were applied to account for our experimental results on anisotropy, mechanical behavior and stress-dependence of the outcrop samples. Our study provides a valuable data set which can feed numerical models to simulate the behavior of the host formation during heat exchange operations in the DeepStor project.

Geomechanics↗
DOE OSTI2020

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Crushing behaviors of unidirectional carbon fiber reinforced plastic composites under dynamic bending and axial crushing loading

The dynamic crushing characteristics of unidirectional carbon fiber reinforced plastic composites under two loading types, dynamic three-point bending and axial crushing, are investigated by experiment and finite element simulation. Hat section samples with two different layup orientations are tested at various impact velocities to investigate its effects. The experiment results show that delamination plays a critical role in dynamic bending deformation. Different layup orientations lead to similar crushing bending behaviors, but remarkable variations in delamination pattern evolutions. In contrast, layup orientation exerts more significant effects on failure modes and energy absorption in axial crushing. A finite element model with multiple layers of thick shell element and cohesive element is established, and simulations are performed for both dynamic three-point bending and axial crushing. The model can capture the crushing behaviors in the two loading types reasonably well. The damage evolution and energy absorption through various mechanisms at different loading conditions are discussed.

42 ENGINEERING↗
DOE OSTI2022

Switching from batch to continuous granulation: A case study of metoprolol succinate ER tablets

Continuous manufacturing (CM) has been used to produce several immediate release drug products. No extended-release (ER) product manufactured employing CM technology has been approved yet. Herein this study investigated the critical aspects of switching from the batch mode of high shear granulation to the continuous operation of twin-screw granulation for extended-release tablets. Metoprolol succinate ER tablets was used as a model ER formulation for this purpose. A central composite design (CCD) was employed to determine the effects of high shear granulator (HSG) parameters, namely impeller speed, granulation time, and binder liquid feeding rate, on the critical granulation characteristics important for product performance. These critical granulation characteristics served as a guide for switching from the batch processing to the continuous operation for achieving the same breaking strength and dissolution for this ER metoprolol tablets. The granulation time was the most critical factor affecting the bulk properties of granules which contributed to tablet dissolution. The higher density and lower compressibility of granules were attained at the longest granulation time of 5.4 min with the fastest liquid feeding rate of 75 g/min. The granules’ density was the primary factor negatively affecting the dissolution of metoprolol tablets. However, the breaking strength of tablets confounded the effect of granules density on metoprolol dissolution. Switching the processing parameters of high shear granulation to twin-screw granulation achieved similar dissolution profiles (F2 greater than 50). The screw speed was not found to affect bulk properties of granules. The root cause of granulation failures in twin-screw granulation, such as premature consolidation, excessive swelling, poor cohesion, inconsistent shearing effects, and formation of deformed agglomerates, were identified. In conclusion, the use of critical granulation characteristics through a performance-based approach of ER tablets facilitated the switching of manufacturing of an ER formulation form batch to continuous operation.

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

Optimization of screw design for continuous wet granulation: A case study of metoprolol succinate ER tablets

This study aimed at understanding the effect of screw design on the critical characteristics of granules and tablets of an extended-release (ER) formulation for twin screw granulation process. The screw design parameters assessed included number of kneading elements (KEs) per kneading zone, distance separating kneading zones, staggering angle (SA) of kneading elements and number of sizing elements (SEs). These input variables were varied using a design of experiment (DoE) approach to manufacture granules. Particle size distribution (PSD), flow and bulk properties of the granules, breaking strength and dissolution of tablets manufactured using these granules were characterized. The results of least square fitting showed that KEs, SA, and SEs of the screws significantly (p -values < 0.05) affected the PSD, cohesion, compressibility (CPS), conditioned bulk density (CBD) and permeability of the granules. The KEs and SEs significantly (p -value < 0.05) affected the dissolution, which was attributed to their effects on CPS and CBD of the granules. The distance between kneading zones had no significant effect on granules and tablet characteristics. Finally, these results may be used to further study the interaction of the identified critical screw design parameters with other processing parameters for continuous manufacturing of this ER matrix-based tablet formulation.

60 APPLIED LIFE SCIENCES↗