Search NASA⌕ Search

SEARCH · Search NASA

Results for “materials and molecular modeling”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Characterizing the Self-Assembly Properties of Monoolein Lipid Isosteres

Living cells feature lipid compartments which exhibit a variety of shapes and structures that assist essential cellular processes. Many natural cell compartments frequently adopt convoluted nonlamellar lipid architectures that facilitate specific biological reactions. Improved methods for controlling the structural organization of artificial model membranes would facilitate investigations into how membrane morphology affects biological functions. Monoolein (MO) is a single-chain amphiphile which forms nonlamellar lipid phases in aqueous solution and has wide applications in nanomaterial development, the food industry, drug delivery, and protein crystallization. However, even if MO has been extensively studied, simple isosteres of MO, while readily accessible, have seen limited characterization. An improved understanding of how relatively minor changes in lipid chemical structure affect self assembly and membrane topology could instruct the construction of artificial cells and organelles for modeling biological structures and facilitate nanomaterial-based applications. Here, we investigate the differences in self-assembly and large-scale organization between MO and two MO lipid isosteres. We show that replacing the ester linkage between the hydrophilic headgroup and hydrophobic hydrocarbon chain with a thioesther or amide functional group results in the assembly of lipid structures with different phases not resembling those formed by MO. Using light and cryo-electron microscopy, small-angle X-ray scattering, and infrared spectroscopy, we demonstrate differences in the molecular ordering and large-scale architectures of the self-assembled structures made from MO and its isosteric analogues. These results improve our understanding of the molecular underpinnings of lipid mesophase assembly and may facilitate the development of MO-based materials for biomedicine and as model lipid compartments.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2013

Raster-Based Approach to Solar Pressure Modeling

An algorithm has been developed to take advantage of the graphics processing hardware in modern computers to efficiently compute high-fidelity solar pressure forces and torques on spacecraft, taking into account the possibility of self-shading due to the articulation of spacecraft components such as solar arrays. The process is easily extended to compute other results that depend on three-dimensional attitude analysis, such as solar array power generation or free molecular flow drag. The impact of photons upon a spacecraft introduces small forces and moments. The magnitude and direction of the forces depend on the material properties of the spacecraft components being illuminated. The parts of the components being lit depends on the orientation of the craft with respect to the Sun, as well as the gimbal angles for any significant moving external parts (solar arrays, typically). Some components may shield others from the Sun. The purpose of this innovation is to enable high-fidelity computation of solar pressure and power generation effects of illuminated portions of spacecraft, taking self-shading from spacecraft attitude and movable components into account. The key idea in this innovation is to compute results dependent upon complicated geometry by using an image to break the problem into thousands or millions of sub-problems with simple geometry, and then the results from the simpler problems are combined to give high-fidelity results for the full geometry. This process is performed by constructing a 3D model of a spacecraft using an appropriate computer language (OpenGL), and running that model on a modern computer's 3D accelerated video processor. This quickly and accurately generates a view of the model (as shown on a computer screen) that takes rotation and articulation of spacecraft components into account. When this view is interpreted as the spacecraft as seen by the Sun, then only the portions of the craft visible in the view are illuminated. The view as shown on the computer screen is composed of up to millions of pixels. Each of those pixels is associated with a small illuminated area of the spacecraft. For each pixel, it is possible to compute its position, angle (surface normal) from the view direction, and the spacecraft material (and therefore, optical coefficients) associated with that area. With this information, the area associated with each pixel can be modeled as a simple flat plate for calculating solar pressure. The vector sum of these individual flat plate models is a high-fidelity approximation of the solar pressure forces and torques on the whole vehicle. In addition to using optical coefficients associated with each spacecraft material to calculate solar pressure, a power generation coefficient is added for computing solar array power generation from the sum of the illuminated areas. Similarly, other area-based calculations, such as free molecular flow drag, are also enabled. Because the model rendering is separated from other calculations, it is relatively easy to add a new model to explore a new vehicle or mission configuration. Adding a new model is performed by adding OpenGL code, but a future version might read a mesh file exported from a computer-aided design (CAD) system to enable very rapid turnaround for new designs

Wright, Theodore W. II↗
DOE OSTI2025

Correlating processing variables to material properties in recycled polypropylene: A data‐driven approach

Abstract Polypropylene (PP) is one of the most widely used plastics, yet its recycling remains limited, with less than 1% of solid waste PP being reprocessed. Mechanical recycling through extrusion is the most practical method, but inconsistent reprocessing conditions introduce variability in material properties. While temperature, screw speed, and residence time influence the thermomechanical stress applied during reprocessing, there are no standardized guidelines for optimizing these parameters. This study examines how these factors shape the properties of recycled PP, using conditions designed to mimic post‐industrial recycled (PIR) scrap. Residence time was measured using colorimetric tracking and correlated with molecular weight, viscosity, and mechanical properties over multiple extrusion cycles. Data‐driven modeling, including response surface methodology, support vector machines, and artificial neural networks, identified processing temperature as the dominant factor in material degradation, followed by residence time. Mechanical properties remained stable, while viscosity decreased predictably with increasing residence time. By linking reprocessing conditions to property evolution, this study provides a method to optimize processing parameters and reduce variability in recycled PP. These findings help manufacturers improve process control, making recycled PP more predictable for reuse in manufacturing. Highlights Study of PIR‐quality PP without additives or compatibilizers. Residence time analysis shows processing temperature drives PP property changes. Mark‐Houwink enables quick molecular weight checks for quality control. Models predict mechanical and rheological shifts in reprocessing. Optimized processing parameters minimize property degradation in recycling.

Estela‐García, John E. [Polymer Engineering Center↗
NASA NTRS2014

Computational Design of Materials: Planetary Entry to Electric Aircraft and Beyond

NASA's projects and missions push the bounds of what is possible. To support the agency's work, materials development must stay on the cutting edge in order to keep pace. Today, researchers at NASA Ames Research Center perform multiscale modeling to aid the development of new materials and provide insight into existing ones. Multiscale modeling enables researchers to determine micro- and macroscale properties by connecting computational methods ranging from the atomic level (density functional theory, molecular dynamics) to the macroscale (finite element method). The output of one level is passed on as input to the next level, creating a powerful predictive model.

Materials Design↗
DOE OSTI2023

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non‐Covalent Interactions in a Uranyl Tetrahalide Model System

Abstract Together with the synthesis and experimental characterization of 14 hybrid materials containing [UO 2 X 4 ] 2− (X=Cl − and Br − ) and organic cations, we report on novel methods for determining correlation trends in their formation enthalpy (ΔH f ) and observed vibrational signatures. Δ H f values were analyzed through isothermal acid calorimetry and a Density Functional Theory+Thermodynamics (DFT+T) approach with results showing good agreement between theory and experiment. Three factors (packing efficiency, cation protonation enthalpy, and hydrogen bonding energy [ ]) were assessed as descriptors for trends in Δ H f . Results demonstrated a strong correlation between and Δ H f , highlighting the importance of hydrogen bonding networks in determining the relative stability of solid‐state hybrid materials. Lastly, we investigate how hydrogen bonding networks affect the vibrational characteristics of uranyl solid‐state materials using experimental Raman and IR spectroscopy and theoretical bond orders and find that hydrogen bonding can red‐shift U≡O stretching modes. Overall, the tightly integrated experimental and theoretical studies presented here bridge the trends in macroscopic thermodynamic energies and spectroscopic features with molecular‐level details of the geometry and electronic structure. This modeling framework forms a basis for exploring 3D hydrogen bonding as a tunable design feature in the pursuit of supramolecular materials by rational design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non‐Covalent Interactions in a Uranyl Tetrahalide Model System

Abstract Together with the synthesis and experimental characterization of 14 hybrid materials containing [UO 2 X 4 ] 2− (X=Cl − and Br − ) and organic cations, we report on novel methods for determining correlation trends in their formation enthalpy (ΔH f ) and observed vibrational signatures. Δ H f values were analyzed through isothermal acid calorimetry and a Density Functional Theory+Thermodynamics (DFT+T) approach with results showing good agreement between theory and experiment. Three factors (packing efficiency, cation protonation enthalpy, and hydrogen bonding energy [ ]) were assessed as descriptors for trends in Δ H f . Results demonstrated a strong correlation between and Δ H f , highlighting the importance of hydrogen bonding networks in determining the relative stability of solid‐state hybrid materials. Lastly, we investigate how hydrogen bonding networks affect the vibrational characteristics of uranyl solid‐state materials using experimental Raman and IR spectroscopy and theoretical bond orders and find that hydrogen bonding can red‐shift U≡O stretching modes. Overall, the tightly integrated experimental and theoretical studies presented here bridge the trends in macroscopic thermodynamic energies and spectroscopic features with molecular‐level details of the geometry and electronic structure. This modeling framework forms a basis for exploring 3D hydrogen bonding as a tunable design feature in the pursuit of supramolecular materials by rational design.

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

The Development of Catalysts for Upgrading of Pyrolysis Vapor for Refinery Feedstocks and Intermediates (CRADA Final Report)

Catalytic fast pyrolysis (CFP) is a versatile technology platform to convert biomass into fungible hydrocarbon transportation fuels and chemical co-products. Key technical barriers to reaching this goal include increasing the product yields and achieving the desired fuel properties for gasoline, diesel, and jet range fuels or blendstocks that would be suitable for introduction into existing refinery unit operations. Overcoming these barriers will require durable catalysts that are effective at upgrading and stabilizing biomass pyrolysis vapors. Towards these goals, this CRADA leveraged NREL experience as a leader in biomass pyrolysis research and Johnson Matthey's (JM) experience as a leader in the production of advanced catalytic materials. The scope spanned CFP catalyst development, characterization, multi-scale reaction testing, and computational modeling. CRADA benefits to DOE, Participant, and U.S. Taxpayer: Assists laboratory in achieving programmatic scope, Uses the laboratory’s core competencies. The purpose of this CRADA was to develop and deploy catalysts for biomass CFP to help achieve cost-competitive biofuels and bio-based products. This was accomplished through a close collaboration between biomass conversion researchers at NREL and catalyst development researchers at JM. Summary of Research Results: Focus Area 1. Foundational research on catalytic conversion and deactivation: Key interactions between pyrolysis vapors and heterogeneous catalysts were probed through catalyst characterization, model compound reaction testing, and atomistic-scale computational modeling. Catalyst development focused on multifunctional materials, which include zeolites, oxides, carbides, and nitrides. Computational modeling identified reaction mechanisms and elucidated surface chemistry to test hypotheses regarding mechanisms of deoxygenation, coupling, cracking, dehydration, coke formation, hydrogen transfer, and aromatic ring reactions. This information was used to design multifunctional catalysts to increase product yields, control product selectivity, and reduce deactivation during CFP and downstream processing steps. The results served to increase fundamental understanding of key catalyst attributes and durability features for the upgrading of biomass pyrolysis vapors. Model compound experiments confirmed the importance of metal-acid bifunctionality for the deoxygenation of lignin-derived phenolic species under hydrodeoxygenation conditions. This insight led to the development of catalysts such as Pt/TiO2 and Mo2C, which were confirmed as high-performing materials during subsequent bench-scale experiments using biomass-derived pyrolysis vapors. This focus area also led to the identification of important catalyst deactivation mechanisms associated with the deposition of inorganic contaminants such as potassium. The molecular-level insight from model compound experiments and computational modeling, shown in Figure 1, informed the development of regeneration procedures that have been shown to be effective for restoration of > 90% of initial catalyst activity. This understanding has subsequently been translated to other catalyst systems, including zeolite materials that can be operated without requirements for co-fed hydrogen.

09 BIOMASS FUELS↗
DOE OSTI2020

Methodology for Generating Covariance Data of Thermal Neutron Scattering Cross Sections

This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Learning Molecular Mixture Property Using Chemistry-Aware Graph Neural Network

Recent advances in machine learning (ML) are expediting materials discovery and design. One significant challenge facing ML for materials is the expansive combinatorial space of potential materials formed by diverse constituents and their flexible configurations. This complexity is particularly evident in molecular mixtures, a frequently explored space for materials, such as battery electrolytes. Owing to the complex structures of molecules and the sequence-independent nature of mixtures, conventional ML methods have difficulties in modeling such systems. Here, we present MolSets, a specialized ML model for molecular mixtures, to overcome the difficulties. Representing individual molecules as graphs and their mixture as a set, MolSets leverages a graph neural network and the deep sets architecture to extract information at the molecular level and aggregate it at the mixture level, thus addressing local complexity while retaining global flexibility. We demonstrate the efficacy of MolSets in predicting the conductivity of lithium battery electrolytes and highlight its benefits in the virtual screening of the combinatorial chemical space. Published by the American Physical Society 2024

Zhang, Hengrui (ORCID:0000000231831654)↗
DOE OSTI2024

A simple method for floating graphene oxide films facilitates nanoscale investigations of ion and water adsorption

Graphene oxide (GO) is a promising material for separations. Nanoscale GO thin films at the air/water interface are excellent experimental models to understand molecular-scale interactions of ions and water with GO. However, the characteristics of GO, such as functional groups and flake size, also affect the thin film properties making it difficult to make systematic studies with GO thin films. This paper reports a simple, reliable, and quick method of preparing ultra-thin GO films, irrespective of their origin, and demonstrates the new opportunities possible with the utilization of this method. The total amount of GO used to form the thin film is significantly less compared to previous examples in the literature, minimizing the dissolved GO in the subphase. X-ray reflectivity (XR) studies show that the majority of the GO film has 1.5 nm thickness over a macroscopic area (~100 cm 2 ) with very small roughness. Sum frequency generation (SFG) spectroscopy measurements show that H 2 O and D 2 O interact differently with GO films, a property that was not observed before. SFG data show that functional groups vary significantly between different commercially available GO samples. The differences are also characterized with XR at high resolution. X-ray fluorescence near total reflection (XFNTR) measurements show that these differences strongly affect ion adsorption and interfacial water behavior near GO, which are vital properties in separation applications. The results pave the way for future studies to elucidate the complex separation mechanisms with GO.

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

Electronic Structure Theory and Novel Materials

This grant supported research on electronic structure and materials theory, with focus on three main issues: (i) novel techniques to deal with correlation in the electronic ground-state, (ii) topological materials, (iii) the phase diagram of lattice spin models. Regarding (i), we applied to the homogeneous electron liquid an approach that we previously developed in the context of molecular systems. In this scheme the electronic occupation probabilities and the natural spin orbitals are used to construct an approximate two-body density matrix for the electronic ground-state. Regarding (ii) we used standard electronic structure methods based on density functional theory to model topological materials and interpret experimental observations. Finally, regarding (iii) we further developed a numerical approach to compute the renormalized couplings within real space renormalization group theory in the context of lattice spin models. The main findings were the following. (i) We found that with our approximate two-body density matrix, which works well for small molecules, is not sufficiently accurate for condensed phase systems. Missing a systematic way of improving on the adopted approximations, we decided not to pursue this approach. (ii) We performed two studies. In one, we investigated the influence of Te defects on the topological properties of a WTe2 monolayer, finding that while Te vacancies, even in modest concentration, destroy the topological character, Te adatoms do not, consistent with a recent experiment. In another study, we predicted Weyl semimetal character and strong anomalous Hall effect in the Heusler compensated ferrimagnet Ti2MnAl. (iii) We developed a new Monte Carlo method to do real space renormalization group calculations for lattice spin models. We subsequently extended the scheme to deal with lattice spin models in presence of quenched disorder, finding that the approach can distinguish systems with finite and strong disorder. In the finite disorder case, the method allows one to find with good approximation the critical coupling distribution and the critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS 2 channel can be significantly altered. Further, we found that larger cations with weaker hydration energy (K + , Rb + , and Cs + ) exhibited a distinctive CN decrease under confinement, while smaller cations (Li + and Na + ) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K + , Rb + , and Cs + ) in MoS 2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS 2 channel.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Antiviral discovery using sparse datasets by integrating experiments, molecular simulations, and machine learning

Computational methods have demonstrated success in identifying virucidal agents, effectively contributing to the discovery of novel virucidal molecules. In this study, we developed a machine learning (ML) model, trained on a small dataset, to predict inhibitors of human enterovirus 71 (EV71), a pathological agent that causes severe disease in children and immunocompromised adults. Despite the dataset’s limitation, comprising of only 36 compounds tested, our ML framework demonstrated significant predictive capability. Notably, experimental validation revealed that five out of the eight compounds predicted by our model from the Chinese cosmetic material list exhibited virucidal activity. The inhibitor effects displayed by the main active compounds were further confirmed by molecular dynamics simulation. This underscores the potential of our AI-driven approach to bypass data constraints in identifying active molecules against viral pathogens.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Direct and Inverse Photoelectron Spectroscopy Evidence for a Revised Picture of Electronic States of Negative Polarons in n‐Doped C 60

Abstract Determining the electronic levels associated with polarons, the fundamental charge carriers in organic semiconductors, is key to understanding the charge transport properties of these materials. Recent findings challenge the traditional view of these electronic levels by highlighting the importance of intra‐molecular Coulomb interactions in polarons. Experimental evidence was previously presented for a revised model of the negative polaron in the case of the polymer semiconductor poly(NDI2OD‐T 2 ); there, the addition of an excess electron was seen to lead to the emergence of a singly occupied state within the energy gap of the undoped material and an unoccupied state above the edge of the conduction states . Here, focus is on a small‐molecule semiconductor, C 60 , and spectral evidence is provided of a similar picture for the new states appearing upon polaron formation. Specifically, direct and inverse photoemission spectroscopy is used to investigate the density of states in C 60 films n‐doped with two dimeric dopants. The Coulomb interaction energy (Hubbard U ) of the C 60 anion is experimentally determined to be ≈1.1 eV, a value that aligns closely with theoretical predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗