Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Growth from a hypercooled melt near absolute stability

The stability of a solid-liquid interface in a hypercooled melt is studied, taking into account attachment kinetics, surface energy, and surface energy in the heat balance. There is a basic-state solution with the planar interface moving at constant speed. Linear-stability theory gives a long-wave absolute-stability limit. Near this point a string model is introduced in which a thermal-boundary-layer approximation is used and an evolution equation for the interface is obtained. In a limiting case this interface equation reduces to a Kuromoto-Sivashinsky equation. Comparison with experimental and numerical results are discussed, and a conceptual picture of unconstrained growth for all undercoolings is addressed.

Umantsev, A.↗
NASA NTRS2011

Wind and Temperature Spectrometry of the Upper Atmosphere in Low-Earth Orbit

Wind and Temperature Spectrometry (WATS) is a new approach to measure the full wind vector, temperature, and relative densities of major neutral species in the Earth's thermosphere. The method uses an energy-angle spectrometer moving through the tenuous upper atmosphere to measure directly the angular and energy distributions of the air stream that enters the spectrometer. The angular distribution gives the direction of the total velocity of the air entering the spectrometer, and the energy distribution gives the magnitude of the total velocity. The wind velocity vector is uniquely determined since the measured total velocity depends on the wind vector and the orbiting velocity vector. The orbiting spectrometer moves supersonically, Mach 8 or greater, through the air and must point within a few degrees of its orbital velocity vector (the ram direction). Pointing knowledge is critical; for example, pointing errors 0.1 lead to errors of about 10 m/s in the wind. The WATS method may also be applied without modification to measure the ion-drift vector, ion temperature, and relative ion densities of major ionic species in the ionosphere. In such an application it may be called IDTS: Ion-Drift Temperature Spectrometry. A spectrometer-based coordinate system with one axis instantaneously pointing along the ram direction makes it possible to transform the Maxwellian velocity distribution of the air molecules to a Maxwellian energy-angle distribution for the molecular flux entering the spectrometer. This implementation of WATS is called the gas kinetic method (GKM) because it is applied to the case of the Maxwellian distribution. The WATS method follows from the recognition that in a supersonic platform moving at 8,000 m/s, the measurement of small wind velocities in the air on the order of a few 100 m/s and less requires precise knowledge of the angle of incidence of the neutral atoms and molecules. The same is true for the case of ion-drift measurements. WATS also provides a general approach that can obtain non-equilibrium distributions as may exist in the upper regions of the thermosphere, above 500 km and into the exosphere. Finally, WATS serves as a mass spectrometer, with very low mass resolution of roughly 1 part in 3, but easily separating atomic oxygen from molecular nitrogen.

Herrero, Federico↗
DOE OSTI2023

Unifying principles for catalytic hydrotreating processes (Final Technical Report)

This project builds on the hypothesis that the hydrotreating processes for the removal of oxygen and sulfur are fundamentally similar at the atomic-scale and existing knowledge from the treatment of petroleum derived feedstock can be leveraged for the design of novel catalysts for the upgrade of bio-oil. We tested this hypothesis by comparing computed potential energy diagrams for hydrodesulfurization (HDS) of thiophene over MoS 2 with hydrodeoxygenation (HDO) of furan over MoO 3 and concluded that certain aspects, such as catalyst promotion with transition metals, are valid strategies for both reactions. On the other hand, we also noticed significant differences in the mechanism for hydrogen (H 2 ) activation, which requires sites with metallic character. While MoS 2 is known to have metallic edge states that can catalyze H 2 dissociation, this elementary step is prohibitively slow on defect-free oxides. Only in the presence of vacancies or by creating metal/oxide interfaces can efficient H 2 activation sites during HDO be formed. The need for bifunctional catalyst when it comes to efficient and selective HDO or dehydrogenation reactions was further corroborated in joint experimental and theoretical studies of the Guerbet reaction for the coupling of biomass derived oxygenates over PdCu alloys, nitrate reduction over In-promoted Pd nanoparticles, and ethylene dehydroaromatization over Ga-exchanged ZSM-5 zeolites. All of these catalytic systems have in common that catalytic sites with distinct functional requirements are needed to create a working catalyst. Detailed computational studies were carried out for HDO of m -cresol and phenol on Ru-modified TiO 2 surfaces, which allowed us to attribute catalytic activity to the metal/oxide interface. A surprising finding was that heterolytic cleavage of the H-H bond across the Ru/TiO 2 interface was critically important, despite lower barriers for homolytic H 2 activation on Ru metal. The explanation lies in the high barriers for hydrogen spillover from Ru to TiO 2 , which becomes unnecessary in the heterolytic activation pathway. Moreover, we also reported that proton transfer steps between metal and oxide sites are mediated by weakly adsorbed surface water. During attempts to develop and validate a kinetic Monte Carlo (kMC) model for HDO reactions at the Ru/TiO 2 interface, it became clear that lateral interactions are paramount to describe realistic surface chemistry and without these interactions, the reduction and hydroxylation behavior from our simulations was inconsistent with reported experiments. To assess the importance of lateral interactions in popular computational catalyst design strategies relying on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends, we compared predicted trends with those obtained from descriptor-based kMC models. We critically evaluated the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. The mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, we concluded that the MKM approach does provide an overall better entry point for computational catalyst design. Overall, this project has led to 11 peer-reviewed publications, and their scientific is impact is well illustrated by their combined 702 citations.

08 HYDROGEN↗
DOE OSTI2022

Modeling Heat Transfer and Reaction Kinetics of Biomass in Pyrolysis Feeding Systems

Biomass feedstocks can be used to produce low carbon-intensity chemicals and liquid fuels by displacing fossil-fuel sources. Pyrolysis is the thermal decomposition of biomass in the absence of oxygen and can be used to produce crude biofuels. A typical feeding mechanism for pyrolysis reactors is by lock hopper followed by a horizontal auger feeder. Particle agglomeration and plugging of the auger screw are regularly occurring phenomena that is detrimental to the process since it results in the need for periodic inspection and cleaning. The underlying mechanisms leading to agglomeration and plugging are unclear, but a working hypothesis is that heat from the reactor raises the temperature of the biomass in the feeder to the point at which preliminary decomposition reactions occur and produce "sticky" products. Using principles of heat and mass transfer, we derived one-dimensional differential equations for heat flow from the pyrolysis reactor through the auger feeder. Temperature profiles were solved from the system of one-dimensional differential equations using the shooting method. Subsequently, we used the temperature profiles as a one-way coupling to a reaction model for the biomass in the feeder. Kinetic models were adopted from the literature and used to predict the formation of intermediates, products, char, extractives, and metaplastics as well as phase changes along the auger feeder. Species were identified that correlated with softening of the biomass, causing it to become sticky and plug the feeder. We will discuss the details of this work, our findings, and the direction of future work that includes fully coupled transport and reaction-kinetics simulations.

BIOMASS FUELS↗
DOE OSTI2025

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstractions by $CH_3\dot{O}$ radical

H-atom abstraction by methoxy radical ($CH_3\dot{O}$) plays an important role in capturing the kinetics of reactions between gasoline components and alcohols. This study focuses on determining the reaction rates and thermodynamic properties of methoxy radical reactions with five gasoline fuel components: n-heptane, iso-octane, 1-hexene, cyclopentane and toluene. Electronic structure calculations were per-formed for all the stationary points with M06-2X/6 -311 ++ g(d,p) method. G3 composite method with atomization method is used for determining Δ f H 0 of all the closed shell and radical species, using which the necessary thermodynamic data of all the species was determined. Coupled cluster theory QCISD(T)/cc-pVXZ (where X = D and T) and Møller-Plesset perturbation theory MP2/cc-pVXZ (where X = D, T and Q) were used to calculate single point energies. Subsequently, rate constants for all hydrogen atom ab-straction channels have been performed using conventional transition state theory with unsymmetric tunneling corrections. A systematic comparison of rates for abstraction from different sites within the same species and same site from different species is done in order to get insights into this reaction class. Here, the computed thermodynamic properties and rate constants were incorporated into a recent gasoline mechanism to investigate the impact of the calculations performed in this work. A shift in predicted NTC (negative temperature coefficient) behavior and significant reduction in model reactivity is observed upon incorporating the rates calculated herein.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Irradiation damage reduces alloy corrosion rate via oxide space charge compensation effects

Radiation effects in materials often compound and accelerate other detrimental phenomena such as embrittlement, oxidation and creep. However, irradiation can also decrease the oxidation rate, for instance with ZrNb alloys nuclear fuel cladding. In this study, we rationalize this observation on Zr-0.5Nb alloy by introducing a mechanism based on oxide space charge modification, resulting from irradiation enhanced Nb clustering. Here, this mechanism is investigated using a multiscale approach: from the macroscale, to determine post-irradiation oxidation kinetics, to the atomic scale, using in-situ atom probe tomography sample oxidation, to observe elemental solute redistribution across the oxide/metal interface. The mechanism is further supported by high resolution transmission electron microscopy characterization and density functional theory calculations. A point defect model is proposed to account for oxide space charge effects and their changes under irradiation. This integrated, multiscale experimental and modeling approach challenges the current paradigm on irradiation effects and how they can potentially improve materials performance in extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

CatTestHub: A benchmarking database of experimental heterogeneous catalysis for evaluating advanced materials

The ability to quantitatively compare newly evolving catalytic materials and technologies is hindered by the widespread availability of catalytic data collected in a consistent manner. While certain catalytic chemistries have been widely studied across decades of scientific research, quantitative comparisons based on literature information is hindered by variability in reaction conditions, types of reported data, and reporting procedures. Here, we present CatTestHub, an open-access database dedicated to benchmarking experimental heterogeneous catalysis data. Combining systematically reported catalytic activity data for selected probe chemistries, with relevant material characterization and reactor configuration information, the database provides a collection of catalytic benchmarks for distinct classes of active site functionality. Through key choices in data access, availability, and traceability, CatTestHub seeks to balance the fundamental information needs of chemical catalysis and the FAIR data design principles. Details of the database architecture and the means through which to navigate it are presented, highlighting examples of catalytic insights readily drawn from the available benchmarking data. In its current iteration, CatTestHub spans over 250 unique experimental data points, collected over 24 solid catalysts, that facilitated the turnover of 3 distinct catalytic chemistries. Here, a roadmap is presented through which to expand the open-access platform that serves as a community wide benchmark, primarily through continuous addition of kinetic information on select catalytic systems by members of the heterogeneous catalysis community at large.

Benchmark↗
DOE OSTI2021

Reversible electrowetting transitions on superhydrophobic surfaces

An electric field applied across the interface has been shown to enable transitions from the Cassie to the Wenzel state on superhydrophobic surfaces with miniature corrugations. Molecular dynamics (MD) simulations manifest the possibility of reversible cycling between the two states when narrow surface wells support spontaneous expulsion of water in the absence of the field. With approximately 1 nm sized wells between the surface asperities, the response times to changes in the electric field are of O(0.1) ns, allowing up to GHz frequency of the cycle. Because of the orientation preferences of interfacial water in contact with the solid, the phenomenon depends on the polarity of the field normal to the interface. The threshold field strength for the Cassie-to-Wenzel transition is significantly lower for the field pointing from the aqueous phase to the surface; however, once in the Wenzel state, the opposite field direction secures tighter filling of the wells. Considerable hysteresis revealed by the delayed water retraction at decreasing field strength indicates the presence of moderate kinetic barriers to expulsion. Known to scale approximately with the square of the length scale of the corrugations, these barriers preclude the use of increased corrugation sizes while the reduction of the well diameter necessitates stronger electric fields. Field-controlled Cassie-to-Wenzel transitions are therefore optimized by using superhydrophobic surfaces with nanosized corrugations. Abrupt changes indicate a high degree of cooperativity reflecting the correlations between the wetting states of interconnected wells on the textured surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Characterization-curing-property studies of HBRF 55A resin formulations

Characterization curing property investigations on HBRF 55A resin formulations are reported. The initial studies on as received cured samples cut from a full-size FWC are reviewed. Inadequacies of as-received and aged samples are pointed out and additional electron microscopic evidence is offered. Characterization of as-received ingredients of HBRF 55A formulation is described. Specifically, Epon 826, Epon 828, EpiRez 5022, RD-2 and various amines, including Tonox and Tonox 60.40, were characterized. Cure kinetics of various formulations are investigated. Changes in physical/thermal properties (viscosity, specific heat, thermal conductivity and density) during cure are described.

Pearce, E. M.↗
DOE OSTI2021

Bulk Diffusion of Cl through O Vacancies in α-Cr 2 O 3 : A Density Functional Theory Study

Many metals form passive oxides which are broken or weakened by anionic-induced degradation leading to material failure and anionic diffusion is an important step of this degradation process. The diffusion of anionic species through oxides involves a combination of diffusion along grain boundaries, cracks, and channels, as well as bulk diffusion via point defects which is the focus of this study. Using density functional theory, we study bulk diffusion of Cl through O vacancies in α -Cr 2 O 3 as a model system for passive metal oxides. Little is known about Cl diffusion and bonding characteristics, so we benchmark our work through comparison to numerous studies on O diffusion in α -Cr 2 O 3 to analyze similarities and differences between the O and Cl diffusion in the passivation and degradation process of α -Cr 2 O 3 respectively. Unlike O diffusion, the lowest diffusion barrier for Cl is cross-plane diffusion between two (0001) planes through a vacant cation site but the much shorter in-plane diffusion path within the same coordination polyhedron has 35% higher barrier. This work provides the basis for considering the contributions of Cl bulk-diffusion in the overall diffusion kinetics of Cl through α -Cr 2 O 3 .

Oware Sarfo, Kofi↗
DOE OSTI2023

Ultrafast Optical Nanoscopy of Carrier Dynamics in Silicon Nanowires

Carrier distribution and dynamics in semiconductor materials often govern their physical properties that are critical to functionalities and performance in industrial applications. The continued miniaturization of electronic and photonic devices calls for tools to probe carrier behavior in semiconductors simultaneously at the picosecond time and nanometer length scales. Here, in this paper, we report pump–probe optical nanoscopy in the visible-near-infrared spectral region to characterize the carrier dynamics in silicon nanostructures. By coupling experiments with the point-dipole model, we resolve the size-dependent photoexcited carrier lifetime in individual silicon nanowires. We further demonstrate local carrier decay time mapping in silicon nanostructures with a sub-50 nm spatial resolution. Our study enables the nanoimaging of ultrafast carrier kinetics, which will find promising applications in the future design of a broad range of electronic, photonic, and optoelectronic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Understanding copper diffusion in CuInSe 2 with first-principles based atomistic and continuum models

We investigate the diffusion of copper in CuInSe 2 using thermodynamic and kinetic models based on density functional theory calculations, attempting to reconcile large differences in reported experimental diffusivities. We find that observations of rapid chemical diffusion can be explained by large thermodynamic factors, which we calculate using a compositionally constrained model of intrinsic point defect formation. We further characterize how copper diffusion coefficients depend on material synthesis conditions and exhibit their variation across the CuInSe 2 secondary phase diagram. In doing so, we identify stable off-stoichiometries that are dominated by either vacancy- or interstitial-mediated diffusion mechanisms. These results are employed in the development of a continuum reaction–diffusion model, which we use to simulate experimental depth profiles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A North Sea in Situ Evaluation of the Fitch Wind Farm Parameterization Within the Mellor-Yamada-Nakanishi-Niino and 3D Planetary Boundary Layer Schemes

Wind resource assessments and wind power forecasts that account for wind farm wakes are sensitive to the choice of planetary boundary layer (PBL) scheme. This work compares the one-dimensional Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme with a three-dimensional PBL (3DPBL) scheme, evaluating predictions made with both schemes against two sets of North Sea in situ observations of wind farm wakes. The optimal PBL scheme varies based on the observations (FINO1 tower vs. aircraft), the quantity of interest (wind speed vs. turbulence kinetic energy [TKE]), and the error metric (bias, centered root mean square error [cRMSE], R2, and earth mover's distance [EMD]). Whereas 3DPBL wind speeds outperform MYNN wind speeds with respect to the cRMSE at the FINO1 site located at a single point within the turbine rotor layer, 3DPBL TKE bias is larger than MYNN TKE bias when compared to aircraft observations taken 100 m above a wind farm. Wind speeds in the aircraft region are ambiguous with regard to which PBL scheme is optimal. Aircraft MYNN wind speeds outperform 3DPBL wind speeds with respect to R2 and cRMSE but underperform with respect to bias and EMD. Future evaluations across broader temporal and spatial scales may offer further insight into model differences.

17 WIND ENERGY↗
NASA NTRS1981

The 1980 rocket coronagraph measurements of the solar wind acceleration region

Spectroscopic measurements of temperatures, densities and flow velocities in the solar wind acceleration region provide critical empirical constraints on solar stellar wind theory. Preliminary results of an analysis of H I Lyman-alpha and white light measurements made on 16 February 1980 in a polar coronal region are reported. The hydrogen kinetic temperatures in the observed region were found to be nearly constant with T sub HI approximately equals 10 degrees K at heliocentric distances between 1.5 and 4 solar radii. The outflow velocities were found to be subsonic indicating that the critical point in the observed region was located at r approximately greater than 4 solar radii.

Withbroe, G. L.↗
NASA NTRS1988

Variable response load limiting device

An energy absorbing device used as a load limiting member in a structure to control its response to applied loads is described. It functions by utilizing a spool assembly having flanged ends and an interior cavity of sufficiently large diameter to cause it to deform plastically at a prescribed load. In application, the spool is utilized as a pivot point for the legs of an airplane seat. When properly designed and integrated into the seat arrangement the spool will twist about its axis, deforming plastically when the impact load exceeds the spool yield value. Through this deformation, the spool absorbs the kinetic energy of the movement of the seat at a substantially constant rate, thereby controlling the level of loads transmitted to the seat occupant. By proper sizing and collection of materials, it is possible to control load response in a predictable manner.

Mcsmith, Dwight D.↗
NASA NTRS2004

Montmorillonite, oligonucleotides, RNA and origin of life

Na-montmorillonite prepared from Volclay by the titration method facilitates the self-condensation of ImpA, the 5'-phosphorimidazolide derivative of adenosine. As was shown by AE-HPLC analysis and selective enzymatic hydrolysis of products, oligo(A)s formed in this reaction are 10 monomer units long and contain 67% 3',5'-phosphodiester bonds (Ferris and Ertem, 1992a). Under the same reaction conditions, 5'-phosphorimidazolide derivatives of cytidine, uridine and guanosine also undergo self-condensation producing oligomers containing up to 12-14 monomer units for oligo(C)s to 6 monomer units for oligo(G)s. In oligo(C)s and oligo(U)s, 75-80% of the monomers are linked by 2',5'-phosphodiester bonds. Hexamer and higher oligomers isolated from synthetic oligo(C)s formed by montmorillonite catalysis, which contain both 3',5'- and 2',5'-linkages, serve as catalysts for the non-enzymatic template directed synthesis of oligo(G)s from activated monomer 2-MeImpG, guanosine 5'-phospho-2-methylimidazolide (Ertem and Ferris, 1996). Pentamer and higher oligomers containing exclusively 2',5'-linkages, which were isolated from the synthetic oligo(C)s, also serve as templates and produce oligo(G)s with both 2',5'- and 3',5'-phosphodiester bonds. Kinetic studies on montmorillonite catalyzed elongation rates of oligomers using the computer program SIMFIT demonstrated that the rate constants for the formation of oligo(A)s increased in the order of 2-mer < 3-mer < 4-mer ... < 7-mer (Kawamura and Ferris, 1994). A decameric primer, dA(pdA)8pA bound to montmorillonite was elongated to contain up to 50 monomer units by daily addition of activated monomer ImpA to the reaction mixture (Ferris, Hill and Orgel, 1996). Analysis of dimer fractions formed in the montmorillonite catalyzed reaction of binary and quaternary mixtures of ImpA, ImpC, 2-MeImpG and ImpU suggested that only a limited number of oligomers could have formed on the primitive Earth rather than equal amounts of all possible isomers (Ertem and Ferris, 2000). Formation of phosphodiester bonds between mononucleotides by montmorillonite catalysis is a fascinating discovery, and a significant step forward in efforts to find out how the first RNA-like oligomers might have formed in the course of chemical evolution. However, as has been pointed out in several publications, these systems should be regarded as models rather than a literal representation of prebiotic chemistry (Orgel, 1998; Joyce and Orgel, 1999; Schwartz, 1999).

RNA/chemical synthesis↗
NASA NTRS1986

The velocity dispersion of the giant molecule clouds: A viscous origin

The observations of interstellar cloud motion show that the cloud velocity dispersion is nearly constant, to within a factor of 2, for clouds covering at least three orders of magnitude in mass. For example, the Giant Molecular Cloud (GMC) with typical masses of approx 5 x 10 to the 5th power solar masses have a one-dimensional planar, cloud-cloud, random velocity dispersion of approx 3-7 km/s. The HI clouds, of approx. 400 solar masses each, on the other hand, have a typical one-dimensional velocity dispersion of approx 6 km/s. Clearly, the clouds are not in kinetic energy equipartion. The GMC spatial distribution in the galactic disk is not that of an isolated, 3-D system; rather, the GMCs exhibit a very thin disk (approx nearly a monolayer) distribution; with the ratio of the diameter of a typical GMC to the vertical scale-height of the GMC distribution being approx 50 pc/150 pc = 0.3. The supernova shocks, which can accelerate the low mass clouds, are extremely ineffective in accelerating the GMCs because of the much larger mass/area ratio for the GMCs. The above points suggest that the GMCs do not constitute an isolated, 3-D system - rather, they indicate that the dynamics of the GMCs is mainly determined by the fact that they are located in a differentially rotating galactic disk, and that, as for particles in planetary rings, viscosity is the primary energy input. Specifically, it is proposed that gravitational scattering of the massive clouds off each other in the differentially rotating galactic disk constitutes an effective gravitational viscosity, which causes an increase in the random kinetic energy of the GMCs at the expense of their ordered, rotational kinetic energy.

Jog, C. J.↗