Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow Structure Interaction”

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

Scaling approach to nuclear structure in high-energy heavy-ion collisions

In high-energy heavy-ion collisions, the initial condition of the produced quark-gluon plasma (QGP) and its evolution are sensitive to collective nuclear structure parameters describing the shape and radial profiles of the nuclei. Here, we find a general scaling relation between these parameters and many experimental observables such as elliptic flow, triangular flow, and particle multiplicity distribution. In particular, the ratios of observables between two isobar systems depend only on the differences of these parameters, but not on the details of the final state interactions, hence offering a new way to constrain the QGP initial condition. Using this scaling relation, we show how the structure parameters of $^{96}_{44}$Ru and $^{96}_{40}$Zr conspire to produce the rich centrality dependences of these ratios, as measured by the STAR Collaboration. Our scaling approach demonstrates that isobar collisions are a precision tool to probe the initial condition of heavy-ion collisions, as well as the collective nuclear structures, including the neutron skin, of the atomic nuclei across energy scales.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Amino Acid Residues Controlling Domain Interaction and Interdomain Electron Transfer in Cellobiose Dehydrogenase

The function of cellobiose dehydrogenase (CDH) in biosensors, biofuel cells, and as a physiological redox partner of lytic polysaccharide monooxygenase (LPMO) is based on its role as an electron donor. Before donating electrons to LPMO or electrodes, an interdomain electron transfer from the catalytic FAD-containing dehydrogenase domain to the electron shuttling cytochrome domain of CDH is required. This study investigates the role of two crucial amino acids located at the dehydrogenase domain on domain interaction and interdomain electron transfer by structure-based engineering. The electron transfer kinetics of wild-type Myriococcum thermophilum CDH and its variants M309A, R698S, and M309A/R698S were analyzed by stopped-flow spectrophotometry and structural effects were studied by small-angle X-ray scattering. The data show that R698 is essential to pull the cytochrome domain close to the dehydrogenase domain and orient the heme propionate group towards the FAD, while M309 is an integral part of the electron transfer pathway – its mutation reducing the interdomain electron transfer 10-fold. Structural models and molecular dynamics simulations pinpoint the action of these two residues on the domain interaction and interdomain electron transfer.

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

Failure investigation and mitigation after experimental research reactor fuel plate deformation in an irradiation device

Experimental research reactor fuel testing is conducted in the Belgian Reactor 2 (BR2) of the Belgian Nuclear Research Centre (SCK CEN) in dedicated irradiation vehicles or rigs. One such vehicle allows flat full-size fuel plates to be irradiated by inserting them into slotted baskets that captures a narrow portion of the longitudinal edges of the plates. The motion of the fuel plates within the baskets is possible within the narrow slots and thus, the plate is considered to be unattached. The design intentionally omits fixing mechanisms of the fuel plates to the baskets to facilitate the inspection and repositioning of the plates between the irradiation cycles and the accommodation of thermal expansion of the plates in the lateral direction. However, loosely inserted fuel plates have weak structural boundary conditions allowing for larger out-of-plane deflections caused by hydrodynamic loads exerted by the flowing coolant, as compared to those of fixed plates. Unexpected large deformations of plates occurred in several irradiation cycles that further resulted in a loss of cladding integrity. These deformations could not be attributed to a single source. This triggered a series of thermal hydraulic, structural, and fluid-structure interaction analyses aiming at understanding the observed phenomenon. The analyses revealed that, for a certain combination of unfavorable manufacturing and assembly tolerances, fuel plate edges could escape out of the slots in the irradiation basket due to the hydrodynamic load. Subsequently, the plate could become wedged inside the basket coolant channel opening. This resulted in reduced coolant flow and accelerated temperature increase and thermal expansion of the plate while under irradiation. This unfavorable feedback loop could then lead to excessive plate surface temperatures, deformed plates and cladding failure, as was observed in the experiments. These analyses not only provided a probable cause of the fuel plate failures, but also resulted in a new and improved design of the irradiation basket to avoid these issues in the future. In conclusion, a series of recent successful irradiations confirm that the sources of failures were identified correctly, and the implemented mitigations were adequate.

BR2↗
DOE OSTI2026

Spatiotemporal structure of edge harmonic oscillation and its role in ELM-free QH-mode at KSTAR

In Quiescent H-mode (QH-mode), edge-localized modes (ELMs) are naturally replaced by a low-n edge harmonic oscillation (EHO), yet the self-regulating transport mechanism driven by the EHO remains insufficiently understood. Using high-spatiotemporal-resolution imaging diagnostics on KSTAR—electron cyclotron emission imaging and broadband electron cyclotron emission—we resolve the eigenmode structure of the EHO and elucidate its regulatory role in edge transport. The EHO is localized within the pedestal near the maximum pressure gradient, propagates in the ion-diamagnetic direction, and its radial envelope expands with increasing shear, suggesting that rotational shear is closely associated with the structural evolution of the EHO. Information-theoretic Transfer Entropy analysis identifies a distinct ‘dual-stabilization’ regulatory interaction pattern: the EHO is associated with enhanced outward energy transport to limit the pressure gradient, while the background shear flow is associated with preferential regulatory influence on the EHO saturation amplitude and energy flux. These results demonstrate that the QH-mode pedestal is sustained by a shear-associated regulation mechanism, where the coupling between rotational shear and mode structure appears closely linked to the sustained ELM-free state.

ECEI↗
DOE OSTI2020

ExaWind: A multifidelity modeling and simulation environment for wind energy

We introduce the open-source ExaWind modeling and simulation environment for wind energy. The primary physics codes of ExaWind are Nalu-Wind and OpenFAST. Nalu-Wind is a wind-focused computational fluid dynamics (CFD) code that is coupled to the whole-turbine simulation code OpenFAST. The ExaWind environment was created under U.S. Department of Energy funding to achieve the highest-fidelity simulations of wind turbines and wind farms to date, with the goal of enabling disruptive changes to turbine and plant design and operation. Innovation will be gleaned through better understanding of the complex flow dynamics in wind farms, including wake evolution and the impact of wakes on downstream turbines and turbulent flow from complex terrain. High-fidelity predictive simulations employ hybrid turbulence models, geometry/boundary-layer-resolving CFD meshes, atmospheric turbulence, nonlinear structural dynamics, and fluid-structure interaction. While there is an emphasis on very high-fidelity simulations (e.g., blade resolved with full fluid-structure coupling), the ExaWind environment supports lower-fidelity modeling capabilities including actuator-line and -disk methods. Important in the development of ExaWind codes is that the codes scale well on today's largest petascale supercomputers and on the next-generation platforms that will enable exascale computing.

17 WIND ENERGY↗
DOE OSTI2025

Bridging microscopic dynamics and rheology in the yielding of charged colloidal suspensions

The yielding of soft materials is critical to many natural and industrial processes, yet experimental insights into microscopic aspects of yielding are limited. This study combines angle X-ray scattering, X-ray photon correlation spectroscopy, and in situ rheology (Rheo-SAXS-XPCS) with fast lubrication dynamics simulations to examine how interparticle interactions influence yielding in charged colloidal suspensions. By tuning attraction through salt addition, we compare repulsive and attractive systems under deformation. Repulsive suspensions yield uniformly with Andrade-like creep and minimal structural change. In contrast, attractive suspensions show complex behaviors, including shear banding, delayed yielding, and resolidification, governed by transient dynamics at shear band interfaces. These results directly link microscopic particle dynamics to macroscopic flow and demonstrate how interaction potentials control rheological behavior. This work offers a framework for designing soft materials with tailored properties for applications in coatings, food processing, drug delivery, and other technologies requiring precise mechanical control.

Molecular Dynamics Simulation↗
DOE OSTI2021

Directing the far-from-equilibrium assembly of nanoparticles in confined liquid crystals by hydrodynamic fields

The assembly of nematic colloids relies on long-range elastic interactions that can be manipulated through external stimuli. Confinement and the presence of a hydrodynamic field alter the defect structures and the energetic interactions between the particles. Here, the assembly landscape of nanoparticles embedded in a nematic liquid crystal confined in a nanochannel under a pressure-driven flow is determined. The dynamics of the liquid crystal tensor alignment field is determined through a Poisson-Bracket framework, namely the Stark–Lubensky equations, coupled with the zero-Reynolds momentum equations and the liquid crystal Landau-de Gennes free energy functional. A second order semi-implicit time integration and a three-dimensional Galerkin finite element method are used to resolve flow and nematic fields under several conditions. In general, the zero Reynolds flow displaces the defects around the particles in the upstream direction and renders the surface anchoring ineffective when the flow strength dominates over the nematic elasticity. More importantly, the potential of mean force for particle assembly is non-monotonic independent of surface anchoring. Our results show that the confinement length scale determines the repulsion/attraction transition between colloids, while the flow strength modifies the static defect structure surrounding the particles and determines the magnitude of the energetic barrier for successful assembly. In the attractive regime, the particles move at different rates through the nematic until one particle eventually catches up with the other. This process occurs against or along the direction of flow depending on the flow strength. Ultimately, these results provide a template for engineering and controlling the transport and assembly of nanoparticles under far-from equilibrium conditions in anisotropic media.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanical and Thermal Forcing for Upslope Flows and Cumulus Convection over the Sierras de Córdoba

Abstract The upslope flow processes affecting the vertical extent of orographic cumulus convection are examined using observations from the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. Specifically, clear air returns from the U.S. Department of Energy (DOE) second-generation C-band scanning Atmospheric Radiation Measurement (ARM) precipitation radar (CSAPR2) are used to characterize the structure and variability of the ridge-normal (i.e., up/downslope) flow components, which transport mass to the crest of Argentina’s Sierras de Córdoba and contribute to convective initiation. Data are compiled for the entire CACTI period (October–April), including days with clear skies, shallow cumuli, cumulus congestus, and deep convection. To examine shared variability among >70 000 radar scans, we use (i) a principal component analysis (PCA) to isolate modes of variability in the upslope flow and (ii) composite analysis based on convective outcomes, determined from GOES-16 satellite observations. These data are contextualized with observed surface sensible heat fluxes, thermodynamic profiles, and synoptic-scale analysis. Results indicate distinct thermally and mechanically forced upslope flow modes, modulated by diurnal heating and synoptic-scale variations, respectively. In some instances, there is a superposition of thermal and mechanical forcing, yielding either deeper or shallower upslope flow. The composite analyses based on satellite data show that successively deeper convective outcomes are associated with successively deeper upslope flow layers that more readily transport mass to the ridge crest in conjunction with lower lifting condensation levels, facilitating convective initiation. These results help to isolate the forcing mechanisms for orographic convection and thus provide a foundation for parameterizing orographic convective processes in coarse resolution models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Energy-dependent implementation of secondary electron emission models in continuum kinetic sheath simulations

The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Highly Accurate Experimental Heave Decay Tests with a Floating Sphere: A Public Benchmark Dataset for Model Validation of Fluid–Structure Interaction

Highly accurate and precise heave decay tests on a sphere with a diameter of 300 mm were completed in a meticulously designed test setup in the wave basin in the Ocean and Coastal Engineering Laboratory at Aalborg University, Denmark. The tests were dedicated to providing a rigorous benchmark dataset for numerical model validation. The sphere was ballasted to half submergence, thereby floating with the waterline at the equator when at rest in calm water. Heave decay tests were conducted, wherein the sphere was held stationary and dropped from three drop heights: a small drop height, which can be considered a linear case, a moderately nonlinear case, and a highly nonlinear case with a drop height from a position where the whole sphere was initially above the water. The precision of the heave decay time series was calculated from random and systematic standard uncertainties. At a 95% confidence level, uncertainties were found to be very low—on average only about 0.3% of the respective drop heights. Physical parameters of the test setup and associated uncertainties were quantified. A test case was formulated that closely represents the physical tests, enabling the reader to do his/her own numerical tests. The paper includes a comparison of the physical test results to the results from several independent numerical models based on linear potential flow, fully nonlinear potential flow, and the Reynolds-averaged Navier–Stokes (RANS) equations. A high correlation between physical and numerical test results is shown. The physical test results are very suitable for numerical model validation and are public as a benchmark dataset.

42 ENGINEERING↗
DOE OSTI2021

Shear-free mixing to achieve accurate temporospatial nanoscale kinetics through scanning-SAXS: ion-induced phase transition of dispersed cellulose nanocrystals

Time-resolved in situ characterization of well-defined mixing processes using small-angle X-ray scattering (SAXS) is usually challenging, especially if the process involves changes of material viscoelasticity. In specific, it can be difficult to create a continuous mixing experiment without shearing the material of interest; a desirable situation since shear flow both affects nanoscale structures and flow stability as well as resulting in unreliable time-resolved data. Here, we demonstrate a flow-focusing mixing device for in situ nanostructural characterization using scanning-SAXS. Given the interfacial tension and viscosity ratio between core and sheath fluids, the core material confined by sheath flows is completely detached from the walls and forms a zero-shear plug flow at the channel center, allowing for a trivial conversion of spatial coordinates to mixing times. With this technique, the time-resolved gel formation of dispersed cellulose nanocrystals (CNCs) was studied by mixing with a sodium chloride solution. It is observed how locally ordered regions, so called tactoids, are disrupted when the added monovalent ions affect the electrostatic interactions, which in turn leads to a loss of CNC alignment through enhanced rotary diffusion. The demonstrated flow-focusing scanning-SAXS technique can be used to unveil important kinetics during structural formation of nanocellulosic materials. However, the same technique is also applicable in many soft matter systems to provide new insights into the nanoscale dynamics during mixing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Heterogeneous Supramolecular Catalysis through Immobilization of Anionic M 4 L 6 Assemblies on Cationic Polymers

Although most of the currently developed supramolecular catalysts that emulate enzymatic reactivity with unique selectivity and activity through specific host-guest interactions work under homogeneous conditions, enzymes in nature can operate under heterogeneous conditions as membrane-bound enzymes. In order to develop such a heterogeneous system, an immobilized chiral supramolecular cluster Ga 4 1 6 ( 2 ) was introduced into cross-linked polymers with cationic functionalities. These heterogeneous supramolecular catalysts were used in aza-Prins and aza-Cope reactions and successfully applied to continuous-flow reactions. Additionally, they showed high durability and maintained high turnovers for long periods of time. In sharp contrast to the majority of examples of heterogenized homogeneous catalysts, the newly developed catalysts showed enhanced activity and robustness compared to those exhibited by the corresponding soluble cluster catalyst. An enantioenriched cluster was also immobilized to enable asymmetric catalysis, and activity and enantioselectivity of the supported chiral catalyst were maintained during recovery and reuse experiments and under a continuous-flow process. Significantly, the structure of the ammonium cations in the polymers affected stability, reactivity, and enantioselectivity, which is consistent with the hypothesis that the cationic moieties in the polymer support interact with cluster as an exohedral protecting shell, thereby influencing their catalytic performance.

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

Quantitative analysis of zonal flow influence on turbulent plasmas driven by trapped electrons

The role of zonal flow (ZF) in the turbulence saturation of trapped electron modes (TEM) in magnetically confined plasmas is revisited. Here, in this study, we examine ZF excitation and saturation mechanisms in TEM turbulence using detailed free-energy transfer diagnostics from nonlinear gyrokinetic simulations. When the ion channel becomes subdominant, electrons shift to transferring energy to the zonal component, which results in zonal flows always playing a significant though subdominant role, in the case of the temperature gradient driven TEM. The velocity-space structure of energy transfer reveals that trapped electrons directly couple to ZF. Regarding the saturation physics, the zonal flow advection, drift-wave–drift-wave interactions, and stable modes at the pump waves are compared. When ZF saturation is weak, drift-wave–drift-wave interactions play a larger role in maintaining turbulence saturation. The findings elucidate the roles of nonzonal to zonal flow coupling, stable modes, and nonzonal wave-wave interactions, and provide a basis for improving reduced models of turbulent transport.

Song, Jiheon [Hanyang University, Seoul (Korea, Re↗
DOE OSTI2023

Improved Platform for the Benchtop Study of Complex Radiological Environments that Affect Fallout Formation: DHS Feasbility Study (Final Report)

Fallout analysis is a crucial field of study as it provides insight into the structural interactions that occur between radiological components and their surroundings following a nuclear event. Due to the complexity of actual fallout samples, setups like the plasma flow reactor (PFR) at LLNL have been developed to synthesize fallout in a controlled environment. While the PFR allows for control of parameters such as initial temperature, cooling rate, and concentration of radiological species, controlling the concentration of water has remained a challenge as samples are introduced into the setup via aqueous solution. To that end, our goal with this project was to develop methodologies to (1) monitor and (2) control the presence of water in the PFR. We found that a commercially available FTIR and MCT detector was sufficient to monitor the presence of both condensed and vapor-phase water, while a desolvating nebulizer could be used to control the amount of water injected into the system. In total, the results of this feasibility study have expanded the ability of the flow reactor to investigate key sensitivities that influence fallout formation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Analysis of the Aeroelastic Dynamics of Lightweight Flexible Variations of the SNL-NRT Turbine

Current trends show that wind turbines are growing in size to meet a rising demand for renewable energy generation, and their upscaled rotors have inherently become more flexible to maintain a proportionally lighter design. This is because larger rotors must be less massive relative to their diameter to minimize the levelized cost of energy (LCOE), which means that blades that are notably less stiff are produced as a result. These structural changes to blades are often reflected in their compromised aeroelastic stability and amplified deformation during operation, which has the potential to decrease the blade’s expected lifetime and the performance of the machine overall. Variations in blade flexibility are also known to influence vortex-wake structures downstream of the turbine, causing patterns of velocity deficit to evolve in ways that affect the performance of other turbines in the farm. This research explores how the increased flexibility of modern utility-scale wind turbine blades influences rotor aeroelastic behavior and interactions with farm flow. High-fidelity simulations of Sandia National Laboratories’ (SNL) National Rotor Testbed (NRT) wind turbine are presented. Flexible variations of the NRT baseline blade are simulated in a variety of realistic operational conditions typically expected at the SNL’s SWiFT facility in Lubbock, Texas. Solutions are then compared to investigate how specific changes to the structural properties of the NRT baseline blade’s design and construction can influence its aeroelastic response at the rotor and the evolution of the turbine’s wake.

Farrell, Alayna (ORCID:000000023555720X)↗
DOE OSTI2024

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

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

bioenergy↗
DOE OSTI2024

Impacts of Irradiation Structural Behavior on Thermal Hydraulics Safety Analysis to Support MURR LEU Conversion

The University of Missouri Research Reactor (MURR) located in Columbia, Missouri is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is actively collaborating with U. S. Department of Energy (DOE) National Nuclear Security Administration (NNSA) Material Management and Minimization (M3) Office of Reactor Conversion and Uranium Supply to convert from the use of highly enriched uranium (HEU; ≥20 wt% U-235) to low-enriched uranium (LEU; <20 wt% U-235) fuel. A new type of very high-density LEU fuel based on an alloy of uranium and 10 wt% molybdenum (U-10Mo) is expected to allow the conversion to LEU of MURR, as well as four other USHPRR. MURR has been working with the Reactor Conversion Pillar at Argonne to perform fuel element design and fuel cycle performance analyses, steady-state thermal hydraulics safety analysis, and accident safety analyses in preparation for the conversion of MURR and to support a preliminary safety analysis report for conversion to LEU fuels. Subsequent analyses have also been performed, including transition cycles where all-fresh LEU fuel elements are introduced upon conversion and progressing through reactor operations the core is brought to equilibrium. Thermal hydraulics safety analyses performed as part of the above have employed an assumption on channel gap reduction due to burnup-related phenomena including fuel swelling, irradiation creep, and oxide layer buildup. Recently, a series of structural analyses have been performed on the MURR LEU fuel plates and an element due to significant differences between the plate and element designs of the MURR HEU and LEU fuels. In addition, NUREG-1537 indicates that structural phenomena are to be evaluated. Two separate types of structural analyses were performed for the MURR LEU fuel element: fluid-structure interaction (FSI) and irradiation thermo-mechanical. The FSI analysis evaluated the effects of hydraulic forces on the MURR LEU fuel element to quantify the flow-induced plate deflection, and a minimal impact to the channel gap thickness was predicted under prototypic and bounding conditions. The irradiation thermo-mechanical analysis evaluated the effects of fuel swelling, irradiation creep, and thermal expansion for the MURR LEU plates and the element for prototypic thermal and irradiation conditions based on a high-fidelity approach multiphysics approach. Overall, this thermo-mechanical analysis predicts smaller gap thickness changes in previously limiting regions. Larger changes are predicted in the middle of channels, and for end channels where power density is not typically a maximum. An additional thermo-mechanical analysis was performed for the outermost HEU fuel plate, which showed a similar magnitude of deflection as the outermost LEU plate. Due to substantial differences between the channel gap reductions assumed for the previous safety analyses and those predicted by the irradiation thermo-mechanical analysis, a need to evaluate their impact on the thermal hydraulics safety analyses arose. This report presents the results from the steady-state safety analyses for normal operation as well as the accident analyses for the two most limiting accident scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗