Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weight function”

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

Radiolytic degradation of formic acid and formate in aqueous solution: modeling the final stages of organic mineralization under advanced oxidation process conditions

The successful use of advanced oxidation processes to treat aqueous solutions containing undesirable organic species requires the degradation of these species to lower molecular weight, lower hazard compounds. Safe application of this technology requires a thorough understanding of the mechanisms of degradation. These oxidative transformations are mainly initiated by the reactions of reactive oxygen species, particularly hydroxyl radicals. These react with organic molecules to generate carbon-centered radicals. In the presence of dissolved oxygen, the carbon-centered radicals are next converted to peroxyl radicals, which then decay to lower molecular weight species by multiple mechanistic pathways. Formic acid and its conjugate base formate are the last stable chemical species produced immediately before the complete mineralization of any organic molecule undergoing oxidative degradation in aqueous solution. Once understood, the radical-induced chemistry of formic acid/formate under these conditions has wide applicability in all advanced oxidation technologies. To develop this quantitative knowledge, here we have performed a series of 60 Co gamma irradiation studies on aqueous formic acid/formate over different pH and solution conditions. The measured species concentration changes, as a function of applied dose, are compared with the predictions of a kinetic computer model constructed from literature reactions and reported rate coefficients. The excellent agreement found between the results and modeling gives confidence in the mechanism presented here and provide the first complete computer model for the radiolytic degradation of formic acid in water.

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

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗
DOE OSTI2022

A Novel Process for Converting Coal to High-Value Polyurethane Products

Battelle has demonstrated a patented process for making high-value, polyurethane (PU) foam from coal, based on preparing liquefied coal via direct liquefaction, converting it to polyols as an intermediate via ozonation, and then making PU foams from these polyols. This process represents a breakthrough in innovative utilization of U.S. coals, and is applicable to bituminous as well as sub-bituminous coals. The resulting PU foam products are projected to have an extremely high value (i.e., over $\$5,000$ /ton), with nearly 100% of carbon utilization from coalderived liquid feedstock, and 31.5% to 43.5% of the carbon in the PU foam polyol product being bio-based. The targeted products represent an extremely large (i.e., over $80 billion/year), existing PU foam market, which could expand into making coatings and adhesives. The process can further help reduce petroleum imports, while improving the economics of PU foam production. This work was completed with funding from the National Energy Technology Laboratory (NETL), with cost share from the State of Ohio’s Ohio Development Services Agency (ODSA) and others, and has advanced the process to 10 kg/day continuous scale and thus to Technology Readiness Level (TRL) 5. A total of 48 coal-based polyols were prepared and evaluated. The initial 28 polyols focused on range finding for ideal conditions. The later 20 polyols were produced as part of process optimizations. These optimizations were targeted around a continuous ozonolysis process to evaluate extended time reactions and to create the necessary intermediate for production of 1-gallon samples of polyol. The most unique attribute of Battelle’s polyol is in the utilization of coal’s aromaticity to gain final foam rigidity. Typically, polyols depend on the isocyanate fraction and cross-linking to gain rigidity. By utilizing coal, we were able to maintain rigidity while reducing the overall hydroxyl value of the polyol. This is important as lower hydroxyl value leads to greater percent weight of the coal-based polyol because less isocyanate is required for foaming. This leads to greater foam cost savings. This report provides the details, process, and process cost models of the conversion of coal to polyols and further to PU foams. Battelle’s process begins with coal liquids. These liquids can be obtained by two processes: coal coking or pyrolysis to produce coal tar, and Battelle’s biobased coal-to-liquids (CTL) process to produce heavy syncrude after liquifying >85% coal. After liquification, Battelle utilizes ozonolysis to create functionalization on the polyaromatic coal structure. The functionalization is then converted to the final polyester polyol through transesterification, or to hydroxyamide polyol through amidification. Equivalent or better standard properties have been obtained for 2 lb/ft 3 density rigid, water and freon-alternative blown foams, including compressive strength, density, R-value, and dimensional stability. Target applications for these foams are insulation, packaging, and energy-absorbing foams. Some exploratory testing also showed promise for adhesives applications. A detailed economic analysis showed that Battelle’s polyol process is economical, at a 140 metric tons per day (MTD) polyol production scale. An attractive return on investment (ROI) at competitive pricing validates the process is ready for a pilot-plant demonstration. A scale-up plan is provided.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Structure of Lipopolysaccharide from Liberibacter crescens Is Low Molecular Weight and Offers Insight into Candidatus Liberibacter Biology

Huanglongbing (HLB) disease, also known as citrus greening disease, was first reported in the US in 2005. Since then, the disease has decimated the citrus industry in Florida, resulting in billions of dollars in crop losses and the destruction of thousands of acres of citrus groves. The causative agent of citrus greening disease is the phloem limited pathogen Candidatus Liberibacter asiaticus. As it has not been cultured, very little is known about the structural biology of the organism. Liberibacter are part of the Rhizobiaceae family, which includes nitrogen-fixing symbionts of legumes as well as the Agrobacterium plant pathogens. To better understand the Liberibacter genus, a closely related culturable bacterium (Liberibacter crescens or Lcr) has attracted attention as a model organism for structural and functional genomics of Liberibacters. Given that the structure of lipopolysaccharides (LPS) from Gram-negative bacteria plays a crucial role in mediating host-pathogen interactions, we sought to characterize the LPS from Lcr. We found that the major lipid A component of the LPS consisted of a pentaacylated molecule with a β-6-GlcN disaccharide backbone lacking phosphate. The polysaccharide portion of the LPS was unusual compared to previously described members of the Rhizobiaceae family in that it contained ribofuranosyl residues. The LPS structure presented here allows us to extrapolate known LPS structure/function relationships to members of the Liberibacter genus which cannot yet be cultured. It also offers insights into the biology of the organism and how they manage to effectively attack citrus trees.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Bullying among children with heart conditions, National Survey of Children’s Health, 2018–2020

Abstract Children with chronic illnesses report being bullied by peers, yet little is known about bullying among children with heart conditions. Using 2018–2020 National Survey of Children’s Health data, the prevalence and frequency of being bullied in the past year (never; annually or monthly; weekly or daily) were compared between children aged 6–17 years with and without heart conditions. Among children with heart conditions, associations between demographic and health characteristics and being bullied, and prevalence of diagnosed anxiety or depression by bullying status were examined. Differences were assessed with chi-square tests and multivariable logistic regression using predicted marginals to produce adjusted prevalence ratios and 95% confidence intervals. Weights yielded national estimates. Of 69,428 children, 2.2% had heart conditions. Children with heart conditions, compared to those without, were more likely to be bullied (56.3% and 43.3% respectively; adjusted prevalence ratio [95% confidence interval] = 1.3 [1.2, 1.4]) and bullied more frequently (weekly or daily = 11.2% and 5.3%; p < 0.001). Among children with heart conditions, characteristics associated with greater odds of weekly or daily bullying included ages 9–11 years compared to 15–17 years (3.4 [2.0, 5.7]), other genetic or inherited condition (1.7 [1.0, 3.0]), ever overweight (1.7 [1.0, 2.8]), and a functional limitation (4.8 [2.7, 8.5]). Children with heart conditions who were bullied, compared to never, more commonly had anxiety (40.1%, 25.9%, and 12.8%, respectively) and depression (18.0%, 9.3%, and 4.7%; p < 0.01 for both). Findings highlight the social and psychological needs of children with heart conditions.

Cardiovascular System & Cardiology↗
DOE OSTI2021

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗
DOE OSTI2024

Copacabana: a probabilistic membership assignment method for galaxy clusters

Cosmological analyses using galaxy clusters in optical/near-infrared photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for all galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89 per cent membership accuracy with a mild dependence on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass–cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities (⁠μ * ⁠) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation at low redshift for a Legacy Survey of Space and Time type survey. These results indicate the potential of Copacabana and μ * to be used in cosmological analyses of optically selected clusters in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Phenomenology of TMD parton distributions in Drell-Yan and 𝑍 0 boson production in a hadron structure oriented approach

We present a first practical implementation of a recently proposed hadron structure oriented (HSO) approach to transverse momentum dependent (TMD) phenomenology applied to Drell-Yan-like processes, including lepton pair production at moderate 𝑄 2 and 𝑍 0 boson production. We compare and contrast general features of our methodology with other common practices and emphasize aspects that we view as improvements. Separately, we discuss the importance of treatments that preserve various theoretical properties of the nonperturbative parts of TMD parametrizations for future applications that have the extraction of hadron structure as a primary goal. These include the HSO’s preservation of a basic TMD parton-model-like framework even while accounting for full TMD factorization and evolution, explicit preservation of the integral relationship between TMD and collinear parton density functions, and the ability to meaningfully compare different theoretical models of nonperturbative TMD parton distributions. In our examples, we show that there is significant sensitivity at moderate 𝑄 2 to both the form of the nonperturbative transverse momentum dependence and the parametrization of collinear parton densities. However, we also find that evolving to 𝑄 2 =𝑀$^2_𝑍$, without fitting, results in a satisfactory postdiction of existing data for 𝑍0 production, nearly independently of the modeling of nonperturbative transverse momentum behavior. We argue that this demonstrates that moderate 𝑄 measurements should be given greater weight than high 𝑄 measurements in extractions of nonperturbative transverse momentum dependence. We also obtain new extractions of the nonperturbative Collins-Soper kernel within the HSO approach. We discuss its features and compare with some earlier extractions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Characterization of Radiation Induced Defects Across Scales

The Whole Powder Pattern Modelling (WPPM) approach of diffraction line profile analysis (LPA) was applied for the first time on PuO 2 in order to understand the evolution of the microstructural properties as a function of process conditions. Crystallite size distributions and defect concentrations were determined by WPPM on three different PuO 2 samples - including one with a bimodal distribution and compared to the result of more common LPA analyses used in PuO 2 literature studies (e.g. Scherrer equation, double Voigt method). The conclusion was that a lognormal distribution of spherical crystallites with dislocations most closely represent the microstructure of PuO 2 produced by heating Pu 2 (C 2 O 4 ) 3 (H 2 O) 6 •3H 2 O between 650°C and 850°C. The volume-weighted mean column length (Lvol) parameter obtained from the double Voigt method for the monomodal samples were the same as the Lvol obtained from the WPPM within uncertainty. The two LPA methods did not agree on the Lvol parameters of the bimodal distribution sample with only the WPPM getting values that were consistent with the individual components. The WPPM can potentially provide more insight into future studies of the influences on the microstructure properties of PuO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ab Initio Polariton Spectra of ZnTPP Molecules Collectively Coupled Inside an Optical Cavity

Exciton-polaritons are quasi-particles formed by the quantum mechanical hybridization of electronic and photonic excitations. Despite extensive investigations, a fundamental understanding of molecular polariton spectra and the polariton delocalization from an ab initio theoretical perspective remains elusive. We simulate experimentally measured linear transmission spectroscopy of many Zinc(II) tetraphenylporphyrin (ZnTPP) molecules collectively coupled to a cavity from first principles. Our theoretical approach incorporates many low-lying electronic excitations in ZnTPP molecules, as well as collective light-matter couplings between ZnTPP and the quantized radiation modes, both of which are shown to be the key to accurately recovering the experimental spectra. We further analyzed to what extent the polariton and dark states are delocalized over many molecules, for the first time, using fully ab initio descriptions of the molecules. We finally investigate the line width as a function of detuning, providing new theoretical insights into the experimentally observed motional narrowing behavior. Our work presents first-ofits- kind theoretical studies on molecular polariton spectra, offering a new perspective on molecular polariton formation in realistic ab initio molecular systems whose rich, many-state nature provides spectral features enabled by the high density of electronic states beyond simple quantum optics models.

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

Quaternary i-MAX Phases (Mo 2/3 RE 1/3 ) 2 AlC (RE: Dy, Tb, Er): Experimental Characterization and First-Principles Insights into their Fundamental Properties

Rare earth (RE)-based materials have unique electronic, magnetic, and optical properties, leading to the recent discovery of atomically layered solids with the chemical formula (M' 2/3 RE 1/3 ) 2 AlC, which have since garnered significant attention in the scientific community. This study aims to synthesize, characterize, and investigate the structural and thermal stability of the RE i-MAX phases. We prepared i-MAX phases using molybdenum (Mo) as M′ and RE elements as Dy, Tb, and Er, namely (Mo 2/3 Dy 1/3 ) 2 AlC, (Mo 2/3 Tb 1/3 ) 2 AlC, and (Mo 2/3 Er 1/3 ) 2 AlC. Structural characterization through x-ray diffraction (XRD) and Raman spectroscopy confirms the formation of the RE-based i-MAX phase, along with the presence of minor impurity phases in the alloys. Thermogravimetric analysis (TGA) conducted up to 1000°C under ambient conditions reveals that the i-MAX phases remain thermally stable up to approximately 450°C, beyond which oxidation leads to a noticeable weight gain in all samples. Differential scanning calorimetry (DSC) measurements during heating and cooling cycles show endothermic and exothermic peaks for (Mo 2/3 Dy 1/3 ) 2 AlC i-MAX in the 410–420°C range, indicating a temperature-induced minor atomic arrangement. In contrast, these peaks are absent in the Tb- and Er-based i-MAX phases. These findings offer valuable insights into the thermal behavior and stability of these i-MAX phases under thermal stress, contributing to a deeper understanding of their unique properties. Furthermore, first-principles density functional theory (DFT) calculations were performed to investigate the electronic and optical properties of the i-MAX phases. The results reveal their metallic nature, with pronounced contributions from Mo and RE elements near the Fermi level and within the conduction band.

Rare earth↗
DOE OSTI2021

Experimental and Theoretical Study of the Electronic Structures of Lanthanide Indium Perovskites LnInO 3

Ternary lanthanide indium oxides LnInO 3 (Ln = La, Pr, Nd, Sm) were synthesized by high-temperature solid-state reaction and characterized by Xray powder diffraction. Rietveld refinement of the powder patterns showed the LnInO 3 materials to be orthorhombic perovskites belonging to the space group Pnma, based on almost-regular InO 6 octahedra and highly distorted LnO 12 polyhedra. Experimental structural data were compared with results from density functional theory (DFT) calculations employing a hybrid Hamiltonian. Valence region X-ray photoelectron and K-shell X-ray emission and absorption spectra of the LnInO 3 compounds were simulated with the aid of the DFT calculations. Photoionization of lanthanide 4f orbitals gives rise to a complex final-state multiplet structure in the valence region for the 4f n compounds PrInO 3 , NdInO 3 , and SmInO 3 , and the overall photoemission spectral profiles were shown to be a superposition of final-state 4f n-1 terms onto the cross-section weighted partial densities of states from the other orbitals. The occupied 4f states are stabilized in moving across the series Pr-Nd-Sm. Band gaps were measured using diffuse reflectance spectroscopy. These results demonstrated that the band gap of LaInO 3 is 4.32 eV, in agreement with DFT calculations. This is significantly larger than a band gap of 2.2 eV first proposed in 1967 and based on the idea that In 4d states lie above the top of the O 2p valence band. However, both DFT and X-ray spectroscopy show that In 4d is a shallow core level located well below the bottom of the valence band. Band gaps greater than 4 eV were observed for NdInO 3 and SmInO 3 , but a lower gap of 3.6 eV for PrInO 3 was shown to arise from the occupied Pr 4f states lying above the main O 2p valence band.

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

Precisely segmented PEEK–ionene + ionic liquid composite membranes for CO 2 separation

Ionenes, polymers with ionic groups incorporated directly within the backbone are a highly versatile class of materials, although they have received much less attention than polyelectrolytes which have ionic groups pendant from the polymer backbone. By designing ionenes that incorporate robust properties of poly(ether ether ketone) (PEEK), we have achieved new imidazolium-containing PEEK–ionene architectures that create opportunities for enhanced CO 2 separation membranes. To achieve these materials, an new imidazole-functionalized PEEK oligomer (ImK(EEK) 2 KIm) was synthesized through facile and straightforward routes. This molecule was then polymerized via condensation reactions with two different aromatic (α,α'-dibromo-p-xylene) and aliphatic (dibromohexyl containing-bis(imidazolium)hexane dibromide salt) comonomers and exchanged with bistriflimide ([Tf 2 N]) anion to obtain two unique PEEK–ionenes containing distinct PEEK and ionic segments: p([ImK(EEK) 2 Im-p-xyl][Tf 2 N] 2 ) and p([ImK(EEK) 2 (ImC 6 ) 3 ][Tf 2 N] 4 ). While the neat PEEK–ionenes exhibited high molecular weight but were not able to form high-quality films, adding stoichiometric amounts of “free” IL (1-methyl-3-butylimidazolium bistriflimide, [C 4 mIm][Tf 2 N]), greatly improved the flexibility and processability of the resultant membranes. Further, the structure–property relationships of bulk PEEK–ionenes and corresponding composites were extensively characterized by different analytical techniques (thermogravimetric analysis, differential scanning calorimetry, X-ray diffraction, and solid-state NMR). In conclusion, the gas separation properties were investigated, with the PEEK–ionenes + IL composites exhibiting CO 2 permeabilities of 14–94 barrer and good CO 2 /N 2 permselectivities of 26–35, indicating that new designs of segmented ionenes and composites with ILs a promising material design strategy for developing gas separation membranes.

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

High Impedance Fault Detection Through Quasi-Static State Estimation: A Parameter Error Modeling Approach

This paper presents a model for detecting high impedance faults using parameter error modeling and a two step per-phase weighted-least squares state estimation process. The proposed scheme leverages the use of Phasor Measurement Units and synthetic measurements to identify per-phase power flow and injection measurements which indicate a parameter error through ?2 Hypothesis Testing applied to the composed measurement error. Although current and voltage waveforms are commonly analyzed for high-impedance fault detection, wide area power flow and injection measurements, which are already inherent to the state estimation process, also show promise for real-world high-impedance fault detection applications. The error distributions after detection share the measurement function error spread observed in proven parameter error diagnostics and can be applied to high-impedance fault identification. Further, this error spread across measurement functions related to the fault will be clearly discerned from measurement error. Case studies are performed on the IEEE 33-Bus Distribution System along with the proposed model in Simulink.

Cooper, Austin↗
DOE OSTI2023

Linear Viscoelastic Properties of Putative Cyclic Polymers Synthesized by Reversible Radical Recombination Polymerization (R3P)

Linear viscoelastic properties in both melt and solution states are reported for a series of poly(3,6-dioxa-1,8-octanedithiol) (polyDODT) made by reversible radical recombination polymerization (R 3 P) under conditions designed to produce linear (LDODT), cyclic (RDODT) and linear-cyclic mixtures (LRDODT). PolyDODT is amorphous (T g < -50 °C) and highly flexible (entanglement molecular weight M e,lin ≈ 1850 g/mol for LDODT). PolyDODT’s low T g and low M e,lin enable characterization over a wide dynamic range and a wide range of dimensionless weight-average molecular weight Z w = M w /M e,lin . Measurements at temperatures from -57 °C to 100 °C provide up to 18 deades of reduced frequency, which is necessary to characterize RDODT melts with Z w from 23 to 300. The two highest molecular weight polymers in the present RDODT series have such high M w (406k and 556k g/mol) that mass spectrometry, NMR spectroscopy, and even chemical assays for chain ends are unable to rule out up to 2 mol% linear contaminant. By studying the samples in solution (using dilution to reduce Z w ) we could compare their dynamics with those of previously established high-purity polystyrene (PS) rings (limited to Z w ≤ 13.6). RDODT solutions with Z w < 15 (concentrations <5wt% for RDODT-406k and 556k) have dynamic moduli G^* that accord with LCCC-purified PS rings in terms of the frequency dependence (including the absence of a plateau), the progression of shapes of G* as a function of Z w , and the linear scaling of their zero-shear viscosity η0 with M w . The shape of G* as a function of Z w for solutions of RDODT-406k and -556k also accord with lower M w RDODT melts (which have ≤ 1.3mol% linear contaminant). Thus, measurement of the linear viscoelastic properties of appropriate concentrations of high M w (>200k g/mol) putative cyclic polymers, in which linear chains evade spectroscopic detection, may provide an alternative means (though not fully proven) of validation of sample purity. When Z w > 15 (including all seven RDODT melts and eight of their solutions), G* has a rubbery plateau. This suggests that the onset of entanglement-like behavior in rings requires 4-5-fold greater Z w than is required for linear chains. Further, the plateau moduli of RDODT samples are indistinguishable from G N o of the corresponding LDODT (melt or matched-concentration solutions). In entangled linear polymers, the observation that G N o is independent of Z w follows from limitations on lateral fluctuations due to neighboring chains becoming independent of position along a given chain. The present results for RDODT suggest that this holds for sufficiently long endless chains, too. While the RDODT have the same G N o as entangled LDODT, when Z w > 60 the terminal relaxation, if reached at all, of RDODT extends to orders of magnitude lower frequency than an entangled linear polymer of the same Z w . Consequently, the viscosity of RDODT with Z w > 60 increases with Z w much more strongly than the 3.4-power observed for entangled linear polymers. Lastly, these novel polymers, with disulfide-linked backbone and broad relaxation time distribution may prove important in relation to biodegradable elastomers and materials with exceptional low-frequency dissipation, extending at least 12 decades below the onset of the rubbery plateau.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2020

Viscoelastic Response of Dispersed Entangled Polymer Melts

Any polymer synthesis routes result in molecular weight dispersity DAI, depending on the polymerization mechanism. The lowest dispersity polymers are those made by anionic and atom-transfer radical polymerization, which exhibit relatively narrow distributions DM=g,„1A172 — 1.02-1.04. This small divergence from monodispersity results in significant number of molecules that differ in their molecular weight from the average and their impact on the viscoelastic response remains an open question. Here the effects of low dispersity on stress relaxation and shear viscosity of entangled melts are studied using a coarse-grained model for polyethylene. Polymer melts with chain length dispersity set to follow a Schulz-Zimm distribution of Dm = 1.0 — 1.16 for an entangled polyethylene melt using molecular dynamic simulations. The systems were studied to times up to 800 ,us which is beyond the terminal time. These systems are compared to those with binary and ternary distributions. The stress relaxation functions were extracted from the Green-Kubo relation and from stress-strain relaxation following a uniaxial extension. We find on the entanglement time scale, both the mean squared displacement and the stress relaxation are independent of D m . The tube diameter and the rubbery plateau for entangled dynamics do not depend on D m . At longer times, the terminal relaxation time decreases as the faster motion of the shorter chains results in constraint release for the longer chains. Probing this low dispersity regime opens the way to evaluate the degree of dispersity that affects mechanical properties.

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

A coarse-grain force field based on quantum mechanics (CGq FF) for molecular dynamics simulation of poly(ethylene glycol)- block -poly(ε-caprolactone) (PEG- b -PCL) micelles

In order to provide the means to predict from molecular dynamics (MD) simulations the structures of copolymer-based micelles in solution, we developed coarse grain force field (CGq FF) parameters for poly(ethylene glycol) (PEG) and for poly(ε-caprolactone) (PCL). A key advance here is the use of quantum mechanics to train the parameters describing the non-bonded (NB) interactions between the CG beads. The functional forms are the same as the MARTINI CG FF so standard MD codes can be used. Our CGq FF describes well the experimentally observed properties for the polymer–air and polymer–water interfaces, indicating the accuracy of the NB interactions. Additionally, the structural properties (density, radius of gyration ( R g ), and end-to-end distance ( h )) match both experiment and all atom (AA) simulations. We illustrate the application of this CGq FF by following the formation of a spherical micelle from 250 chains of PEG 23 - b -PCL 9 diblock copolymer, each block with molecular weight of 1000 Daltons (10 500 beads, corresponding to 123 250 atoms), in a water box with 119 139 water beads (426 553 water molecules).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗