SEARCH · Search NASA

Results for “coupling coefficient”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Mesoscopic modeling of transport resistances in a polymer-electrolyte fuel-cell catalyst layer: Analysis of hydrogen limiting currents

Understanding transport resistances in a polymer-electrolyte fuel cell (PEFC) catalyst layer (CL) is essential to mitigate the unexpected voltage loss when using low loadings of precious metals. In this paper, we explore through mesoscopic modeling the quantification analyses of the transport resistances in CL as derived using hydrogen-pump limiting current. Numerical treatments on the conjugated interfacial conditions at interfaces of ionomer/pore and Pt/ionomer are proposed to describe the mesoscopic transport processes of hydrogen and proton. Characterizations of the reconstructed microstructure of CL are performed. Parameter analyses on the influences of the critical transport properties such as the permeation coefficient and the dissolution and adsorption reaction rates at the surfaces of ionomer/pore and Pt/ionomer on the local transport resistance are presented. It is found that the local transport resistance is mainly originated from the diffusion resistance of the ionomer thin-film, which is more resistive than its bulk analogue with its permeation coefficient fitted to be 5.9% of the bulk one. The interfacial transport resistances and the diffusion resistance are coupled. The local transport resistance increases with I/C ratio due to thicker ionomer coated on the particles. Higher Pt/C ratio and bare carbon fraction lead to higher local transport resistance since the ionomer loading relative to Pt roughness factor decreases. The local transport resistance decreases with the porosity. Furthermore, the contribution of pores to the CL resistance, which decreases with the porosity, is comparatively small at low loadings.

25 ENERGY STORAGE↗

DOE OSTI2015

Effects of radiative heat transfer on the turbulence structure in inert and reacting mixing layers

We use large-eddy simulation to study the interaction between turbulence and radiative heat transfer in low-speed inert and reacting plane temporal mixing layers. An explicit filtering scheme based on approximate deconvolution is applied to treat the closure problem arising from quadratic nonlinearities of the filtered transport equations. In the reacting case, the working fluid is a mixture of ideal gases where the low-speed stream consists of hydrogen and nitrogen and the high-speed stream consists of oxygen and nitrogen. Both streams are premixed in a way that the free-stream densities are the same and the stoichiometric mixture fraction is 0.3. The filtered heat release term is modelled using equilibrium chemistry. In the inert case, the low-speed stream consists of nitrogen at a temperature of 1000 K and the highspeed stream is pure water vapour of 2000 K, when radiation is turned off. Simulations assuming the gas mixtures as gray gases with artificially increased Planck mean absorption coefficients are performed in which the large-eddy simulation code and the radiation code PRISSMA are fully coupled. In both cases, radiative heat transfer is found to clearly affect fluctuations of thermodynamic variables, Reynolds stresses, and Reynolds stress budget terms like pressure-strain correlations. Source terms in the transport equation for the variance of temperature are used to explain the decrease of this variance in the reacting case and its increase in the inert case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Experimental observation of electron bounce resonance through electron energy distribution measurement in a finite size inductively coupled plasma

The electron bounce resonance was experimentally investigated in a low pressure planar inductively coupled plasma. The electron energy probability functions (EEPFs) were measured at different chamber heights and the energy diffusion coefficients were calculated by the kinetic model. It is found that the EEPFs begin to flatten at the first electron bounce resonance condition, and the plateau shifts to a higher electron energy as the chamber height increases. The plateau which indicates strong electron heating corresponds not only to the electron bounce resonance condition but also to the peaks of the first component of the energy diffusion coefficients. As a result, the plateau formation in the EEPFs is mainly due to the electron bounce resonance in a finite inductive discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Theoretical formulation of phase space microbunching instability in the presence of intrabeam scattering for single-pass or recirculation accelerators

Microbunching instability (MBI) has been one of the most challenging issues in designs of high-brightness beam transport lines for single-pass or recirculating accelerators. Although the intrabeam scattering (IBS) has long been studied in lepton or hadron storage rings as a slow diffusion process or in high-intensity proton linear accelerators as one mechanism for the beam halo, the effects of IBS on single-pass or recirculating electron accelerators have drawn attention only in the recent two decades due to emergence of linac-based or energy-recovery-linac-based fourth-generation light sources, which require high-quality electron beams during the beam transport. In this paper we develop a theoretical formulation of microbunching instability in the presence of IBS for single-pass or recirculation accelerators. To quantify MBI with inclusion of IBS, we start from the Vlasov-Fokker-Planck (VFP) equation, combining both collective interactions and incoherent IBS effects. The linearized VFP equation and the corresponding friction and diffusion coefficients are derived. The evolutions of the resultant density and energy modulations are formulated as a set of coupled integral equations. The theoretical formulation is then applied to a recirculating beamline design. The results from the semianalytical calculation are compared and show good agreement with massive particle tracking simulations.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Quarkonium polarization in medium from open quantum systems and chromomagnetic correlators

Here, by considering the Markovian condition and applying the Wigner transformation upon the diagonal spin components of the quarkonium density matrix with the semiclassical expansion, we systematically derive the Boltzmann transport equation for quarkonia with polarization dependence in the quantum optical limit. Unlike the spin-independent collision terms governed by certain chromoelectric field correlators, new gauge invariant correlators of chromomagnetic fields determine the recombination and dissociation terms with polarization dependence at the order we are working. We also derive a Lindblad equation describing the in-medium transitions between spin-singlet and spin-triplet heavy quark-antiquark pairs in the quantum Brownian motion limit. The Lindblad equation is governed by new transport coefficients defined in terms of the chromomagnetic field correlators. Our formalism is generic and valid for both weakly coupled and strongly coupled quark gluon plasmas. It can be further applied to study spin alignment of vector quarkonia in heavy ion collisions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Comparative Evaluation of U.S. Brand and Generic Intravenous Sodium Ferric Gluconate Complex in Sucrose Injection: Physicochemical Characterization

The objective of this study was to evaluate physicochemical equivalence between brand (i.e., Ferrlecit) and generic sodium ferric gluconate (SFG) in sucrose injection by conducting a series of comparative in vitro characterizations using advanced analytical techniques. The elemental iron and carbon content, thermal properties, viscosity, particle size, zeta potential, sedimentation coefficient, and molecular weight were determined. There was no noticeable difference between brand and generic SFG in sucrose injection for the above physical parameters evaluated, except for the sedimentation coefficient determined by sedimentation velocity analytical ultracentrifugation (SV-AUC) and molecular weight by asymmetric field flow fractionation-multi-angle light scattering (AFFF-MALS). In addition, brand and generic SFG complex products showed comparable molecular weight distributions when determined by gel permeation chromatography (GPC). The observed minor differences between brand and generic SFG, such as sedimentation coefficient, do not impact their biological activities in separate studies of in vitro cellular uptake and rat biodistribution. Coupled with the ongoing clinical study comparing the labile iron level in healthy volunteers, the FDA-funded post-market studies intended to illustrate comprehensive surveillance efforts ensuring safety and efficacy profiles of generic SFG complex in sucrose injection, and also to shed new light on the approval standards on generic parenteral iron colloidal products.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

An efficient algorithm for incompressible N-phase flows

We present an efficient algorithm within the phase field framework for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids, with possibly very different physical properties such as densities, viscosities, and pairwise surface tensions. The algorithm employs a physical formulation for the N-phase system that honors the conservations of mass and momentum and the second law of thermodynamics. We present a method for uniquely determining the mixing energy density coefficients involved in the N-phase model based on the pairwise surface tensions among the N fluids. Our numerical algorithm has several attractive properties that make it computationally very efficient: (i) it has completely de-coupled the computations for different flow variables, and has also completely de-coupled the computations for the (N−1) phase field functions; (ii) the algorithm only requires the solution of linear algebraic systems after discretization, and no nonlinear algebraic solve is needed; (iii) for each flow variable the linear algebraic system involves only constant and time-independent coefficient matrices, which can be pre-computed during pre-processing, despite the variable density and variable viscosity of the N-phase mixture; (iv) within a time step the semi-discretized system involves only individual de-coupled Helmholtz-type (including Poisson) equations, despite the strongly-coupled phase–field system of fourth spatial order at the continuum level; (v) the algorithm is suitable for large density contrasts and large viscosity contrasts among the N fluids. Extensive numerical experiments have been presented for several problems involving multiple fluid phases, large density contrasts and large viscosity contrasts. In particular, we compare our simulations with the de Gennes theory, and demonstrate that our method produces physically accurate results for multiple fluid phases. We also demonstrate the significant and sometimes dramatic effects of the gravity, density ratios, pairwise surface tensions, and drop sizes on the N-phase configurations and dynamics. The numerical results show that the method developed herein is capable of dealing with N-phase systems with large density ratios, large viscosity ratios, and pairwise surface tensions, and that it can be a powerful tool for studying the interactions among multiple types of fluid interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Strong enhancement of surface diffusion by nonlinear surface acoustic waves

The phenomenon of acoustic activation of surface diffusion is investigated in a combined computational and experimental study. The ability of pulsed laser-generated surface acoustic waves (SAWs) to enhance the mobility of small atomic clusters is demonstrated by directly tracking, with fluorescence microscopy, individual Au 8 clusters moving on a (111) silicon substrate. In this work, a 19-fold increase in the effective diffusion coefficient is measured in room temperature experiments in the presence of SAWs generated by nanosecond pulse laser irradiation at a 100 Hz repetition rate. A strong enhancement of cluster mobility by SAWs is also observed in large-scale molecular dynamics simulations of surface diffusion of small atomic clusters. The analysis of the computational results demonstrates that the nonlinear sharpening of SAWs and the corresponding enrichment of the SAW spectra by high frequency harmonics which are capable of dynamic coupling to the cluster vibrations are responsible for the efficient acoustic activation of surface mobility in the simulations. The increase in the effective diffusion coefficient is proportional to the number of the SAW pulses passing through the diffusion region per unit time and a dramatic 4500-fold diffusion enhancement (corresponds to an equivalent temperature increase by ~10 3 K ) is predicted in the simulations for 15 GHz SAWs. The ability of SAWs to affect atomic-level surface processes has far-reaching implications for the design of new techniques where the acoustic energy serves as an effective substitution for thermal activation in applications where heating must be avoided or rapid switching of surface conditions is required.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2017

Ideas and perspectives: how coupled is the vegetation to the boundary layer?

Understanding the sensitivity of transpiration to stomatal conductance is critical to simulating the water cycle. This sensitivity is a function of the degree of coupling between the vegetation and the atmosphere and is commonly expressed by the decoupling factor. The degree of coupling assumed by models varies considerably and has previously been shown to be a major cause of model disagreement when simulating changes in transpiration in response to elevated CO 2 . The degree of coupling also offers us insight into how different vegetation types control transpiration fluxes, which is fundamental to our understanding of land–atmosphere interactions. To explore this issue, we combined an extensive literature summary from 41 studies with estimates of the decoupling coefficient estimated from FLUXNET data. We found some notable departures from the values previously reported in single-site studies. There was large variability in estimated decoupling coefficients (range 0.05–0.51) for evergreen needleleaf forests. This is a result that was broadly supported by our literature review but contrasts with the early literature which suggests that evergreen needleleaf forests are generally well coupled. Estimates from FLUXNET indicated that evergreen broadleaved forests were the most tightly coupled, differing from our literature review and instead suggesting that it was evergreen needleleaf forests. We also found that the assumption that grasses would be strongly decoupled (due to vegetation stature) was only true for high precipitation sites. These results were robust to assumptions about aerodynamic conductance and, to a lesser extent, energy balance closure. Hence, these data form a benchmarking metric against which to test model assumptions about coupling. Our results identify a clear need to improve the quantification of the processes involved in scaling from the leaf to the whole ecosystem. Progress could be made with targeted measurement campaigns at flux sites and greater site characteristic information across the FLUXNET network.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Phenomenological Constitutive Modeling of High-Temperature Flow Behavior Incorporating Individual and Coupled Effects of Processing Parameters in Super-austenitic Stainless Steel

A novel constitutive model has been developed for predicting flow responses of super-austenitic stainless steel over a wide range of strains (0.05-0.6), temperatures (1173-1423 K) and strain rates (0.001-1 s{sup −1}). Further, the predictability of this new model has been compared with the existing Johnson–Cook (JC) and modified Zerilli–Armstrong (M-ZA) model. The JC model is not befitted for flow prediction as it is found to be exhibiting very high (~ 36%) average absolute error (δ) and low (~ 0.92) correlation coefficient (R). On the contrary, the M-ZA model has demonstrated relatively lower δ (~ 13%) and higher R (~ 0.96) for flow prediction. The incorporation of couplings of processing parameters in M-ZA model has led to exhibit better prediction than JC model. However, the flow analyses of the studied alloy have revealed the additional synergistic influences of strain and strain rate as well as strain, temperature, and strain rate apart from those considered in M-ZA model. Hence, the new phenomenological model has been formulated incorporating all the individual and synergistic effects of processing parameters and a ‘strain-shifting’ parameter. The proposed model predicted the flow behavior of the alloy with much better correlation and generalization than M-ZA model as substantiated by its lower δ (~ 7.9%) and higher R (~ 0.99) of prediction.

36 MATERIALS SCIENCE↗

DOE OSTI2014

Experimental characterization and modeling of non-linear coupling of the LHCD power on Tore Supra

To achieve steady state operation on future tokamaks, in particular on ITER, the unique capability of a LHCD system to efficiently drive off-axis non-inductive current is needed. In this context, it is of prime importance to study and master the coupling of LH wave to the core plasma at high power density (tens of MW/m{sup 2}). In some specific conditions, deleterious effects on the LHCD coupling are sometimes observed on Tore Supra. At high power the waves may modify the edge parameters that change the wave coupling properties in a non-linear manner. In this way, dedicated LHCD experiments have been performed using the LHCD system of Tore Supra, composed of two different conceptual designs of launcher: the Fully Active Multijunction (FAM) and the new Passive Active Multijunction (PAM) antennas. A nonlinear interaction between the electron density and the electric field has been characterized in a thin plasma layer in front of the two LHCD antennas. The resulting dependence of the power reflection coefficient with the LHCD power, leading occasionally to trips in the output power, is not predicted by the standard linear theory of the LH wave coupling. Therefore, it is important to investigate and understand the possible origin of such non-linear effects in order to avoid their possible deleterious consequences. The PICCOLO-2D code, which self-consistently treats the wave propagation in the antenna vicinity and its interaction with the local edge plasma density, is used to simulate Tore Supra discharges. The simulation reproduces very well the occurrence of a non-linear behavior in the coupling observed in the LHCD experiments. The important differences and trends between the FAM and the PAM antennas, especially a larger increase in RC for the FAM, are also reproduced by the PICCOLO-2D simulation. The working hypothesis of the contribution of the ponderomotive effect in the non-linear observations of LHCD coupling is therefore validated through this comprehensive modeling for the first time on the FAM and PAM antennas on Tore Supra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Precision test of gauge/gravity duality in D0-brane matrix model at low temperature

We test the gauge/gravity duality between the matrix model and type IIA string theory at low temperatures with unprecedented accuracy. To this end, we perform lattice Monte Carlo simulations of the Berenstein-Maldacena-Nastase (BMN) matrix model, which is the one-parameter deformation of the Banks-Fischler-Shenker-Susskind (BFSS) matrix model, taking both the large N and continuum limits. We leverage the fact that sufficiently small flux parameters in the BMN matrix model have a negligible impact on the energy of the system while stabilizing the flat directions so that simulations at smaller N than in the BFSS matrix model are possible. Hence, we can perform a precision measurement of the large N continuum energy at the lowest temperatures to date. The energy is in perfect agreement with supergravity predictions including estimations of α'-corrections from previous simulations. At the lowest temperature where we can simulate efficiently (T = 0.25λ 1/3 , where λ is the ’t Hooft coupling), the difference in energy to the pure supergravity prediction is less than 10%. Furthermore, we can extract the coefficient of the 1/N 4 corrections at a fixed temperature with good accuracy, which was previously unknown.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Nonrad: Computing nonradiative capture coefficients from first principles

Point defects in semiconductor crystals provide a means for carriers to recombine nonradiatively. This recombination process impacts the performance of devices. We present the Nonrad code that implements the first-principles approach of Alkauskas et al. (2014) [8] for the evaluation of nonradiative capture coefficients based on a quantum-mechanical description of the capture process. An approach for evaluating electron-phonon coupling within the projector augmented wave formalism is presented. We also show that the common procedure of replacing Dirac delta functions with Gaussians can introduce errors into the resulting capture rate, and implement an alternative scheme to properly account for vibrational broadening. Finally, we assess the accuracy of using an analytic approximation to the Sommerfeld parameter by comparing with direct numerical evaluation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Emergency core cooling system performance criteria for Multi-Layered Silicon Carbide nuclear fuel cladding

Multi-layered Silicon Carbide (SiC) cladding has been proposed to replace the current Zircaloy (Zr) cladding in existing nuclear power plants. Due to its ceramic property, SiC may have a different failure mechanism from that of metallic Zr clad. This characteristic should be considered in designing the plant’s safety system performance, such that the fuel’s integrity is reasonably assured during accident scenarios. This paper aims at optimizing the Emergency Core Cooling System (ECCS) performance criteria in mitigating LBLOCA accidents at a Pressurized Water Reactor (PWR) using the uranium dioxide fuel and duplex or triplex SiC cladding structures. Clad failures by brittle fractures and thermal degradation were modeled in a fuel analysis code. This fuel model was coupled with the RELAP5 plant model by adapting the fuel’s mesh configuration and formulating the effective heat transfer coefficients. A model to quantify the likelihood of clad functional failure as a combination of individual layer fractures at various axial positions was proposed. Updated SiC material properties were gathered from recent reference studies. Clad reliability estimates using the proposed model were compared with the estimates obtained using other models. The results were used to relax the criterion of required actuation timing of ECCS, such that the risk of the new plant design is less or equal to that of the initial plant configuration. Triplex clad was more reliable than duplex clad and Zr, reducing the conditional core damage probability by three orders of magnitude in the reference LBLOCA scenario. Compressive stresses due to irradiation swelling played the determining role in the vulnerabilities of each clad design. The duplex clad was vulnerable to pressure-loading stresses while triplex clad was prone to fail due to thermal shocks experienced during the LOCA reflood phase. Additionally, it was concluded that the most vulnerable axial fuel rod location in the reference LBLOCA scenario was located at the highest axial power peaking location. A sensitivity study on the ECCS timing density was conducted. The results showed that ECCS actuation timing could be extended from the initial technical specification, without increasing the fuel damage likelihood. The time extension was distributed in the order of several hundreds of seconds, which was comparable with time extensions offered by FeCrAl and Cr-coated-Zr clad designs in LBLOCA scenarios. This time extension may potentially reduce the failure probabilities of EDGs and pumps, hence improving the plant safety further.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Thermomechanical analysis and modeling of a high-temperature light trapping planar cavity solar receiver

The development of durable particle-based high-temperature solar receivers is critical for advancing concentrating solar-thermal (CST) technologies to enable high-efficiency power generation and industrial process heat. Here, this study presents a computational framework to evaluate the thermomechanical performance of a proposed enclosed light-trapping planar cavity receiver designed for particle-based thermal energy systems. The receiver incorporates absorptive cavities and fluidized particle-bed channels to enhance heat capture and reduce thermal losses. Finite element analysis (FEA) is employed to assess stress, strain, and creep-fatigue behavior under concentrated solar flux using realistic thermal boundary conditions derived from coupled system models and experimental assembly parameters. The analysis investigates the influence of particle-to-wall heat transfer coefficients (HTC) ranging from 800 to 1800 W/m 2 .K on the thermomechanical response of six candidate high-temperature alloys: Alloy 740H, Alloy 282, Alloy 617, 316H, Alloy 230, and 800H. Results show that increasing HTC reduces thermal gradients, leading to lower stresses and strains and extended minimum predicted creep life. While all materials satisfy fatigue life requirements under the investigated conditions, significant differences in creep resistance are observed. Alloy 740H consistently exhibits the longest minimum predicted creep life and the most favorable durability margins, followed by Alloy 282, with the remaining materials showing reduced creep resistance under identical loading. The reported creep lives are conservative lower-bound estimates intended for comparative material evaluation. This framework highlights the critical roles of material selection and geometry optimization in improving mechanical durability and reliability of solar-thermal receivers, forming a foundation for future experimental validation and design optimization.

14 SOLAR ENERGY↗

DOE OSTI2017

Molecular Chemistry of Atmospheric Brown Carbon Inferred from a Nationwide Biomass Burning Event

Lag Ba'Omer, a nationwide bonfire festival in Israel, was chosen as a case study to investigate the influence of a major biomass burning event on the light absorption properties of atmospheric brown carbon (BrC). The chemical composition and optical properties of BrC chromophores were investigated using a high performance liquid chromatography (HPLC) platform coupled to photo diode array (PDA) and high resolution mass spectrometry (HRMS) detectors. Substantial increase of BrC light absorption coefficient was observed during the night-long biomass burning event. Most chromophores observed during the event were attributed to nitroaromatic compounds, comprising 28 elemental formulas of at least 63 structural isomers. The NAC, in combination, accounted for 50-80% of the total visible light absorption (> 400 nm) by solvent extractable BrC. The results highlight that NAC, particular nitrophenols, are important light absorption contributors of biomass burning organic aerosol (BBOA), suggesting that night time chemistry of ·NO 3 and N 2 O 5 with particles may play a significant role in atmospheric transformations of BrC. Nitrophenols and related compounds were especially important chromophores of BBOA. The absorption spectra of the BrC chromophores are influenced by the extraction solvent and solution pH, implying that the aerosol acidity is an important factor controlling the light absorption properties of BrC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Superconductivity in CaBi 2

We observed superconductivity with critical temperature T c = 2.0 K in self-flux-grown single crystals of CaBi 2 . This material adopts the ZrSi 2 structure type with lattice parameters a = 4.696(1) Å, b = 17.081(2) Å and c = 4.611(1) Å. The crystals of CaBi 2 were studied by means of magnetic susceptibility, specific heat and electrical resistivity measurements. The heat capacity jump at T c is ΔC/γT c = 1.41, confirming bulk superconductivity; the Sommerfeld coefficient γ = 4.1 mJ mol -1 K -2 and the Debye temperature Θ D = 157 K. The electron–phonon coupling strength is λ el–ph = 0.59, and the thermodynamic critical field H c is low, between 111 and 124 Oe CaBi 2 is a moderate coupling type-I superconductor. Our results of electronic structure calculations are reported and charge densities, electronic bands, densities of states and Fermi surfaces are discussed, focusing on the effects of spin–orbit coupling and electronic property anisotropy. Furthermore, we find a mixed quasi-2D + 3D character in the electronic structure, which reflects the layered crystal structure of the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2014

Viscosity, relaxation time, and dynamics within a model asphalt of larger molecules

The dynamics properties of a new “next generation” model asphalt system that represents SHRP AAA-1 asphalt using larger molecules than past models is studied using molecular simulation. The system contains 72 molecules distributed over 12 molecule types that range from nonpolar branched alkanes to polar resins and asphaltenes. Molecular weights range from 290 to 890 g/mol. All-atom molecular dynamics simulations conducted at six temperatures from 298.15 to 533.15 K provide a wealth of correlation data. The modified Kohlrausch-Williams-Watts equation was regressed to reorientation time correlation functions and extrapolated to calculate average rotational relaxation times for individual molecules. The rotational relaxation rate of molecules decreased significantly with increasing size and decreasing temperature. Translational self-diffusion coefficients followed an Arrhenius dependence. Similar activation energies of ∼42 kJ/mol were found for all 12 molecules in the model system, while diffusion prefactors spanned an order of magnitude. Viscosities calculated directly at 533.15 K and estimated at lower temperatures using the Debye-Stokes-Einstein relationship were consistent with experimental data for asphalts. The product of diffusion coefficient and rotational relaxation time showed only small changes with temperature above 358.15 K, indicating rotation and translation that couple self-consistently with viscosity. At lower temperatures, rotation slowed more than diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗