Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion length”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Performance characteristics of two annular dump diffusers using suction-stabilized vortex flow control

The two diffusers employed in the investigation had the same overall area ratio but different prediffuser area ratios and suction slot geometries. Velocity profile and diffuser pressure recovery performance data were obtained at ambient pressure and temperature, with inlet Mach numbers ranging from 0.18 to 0.41 and suction rate varying from zero to 18% of total inlet mass flow rate. On the basis of the reported investigation it is concluded that suction stabilized vortex flow diffusers show promise for application in combustors because of relatively high static pressure recovery and low total pressure loss obtained in a short length. Performance obtained using a narrow angle (7 degree) prediffuser was superior to that obtained with a prediffuser having a 14 degree included angle.

Juhasz, A. J.↗
NASA NTRS1980

Ice crystal growth in a dynamic thermal diffusion chamber

Ice crystals were grown in a supersaturated environment produced by a dynamic thermal diffusion chamber, which employed two horizontal plates separated by a distance of 2.5 cm. Air was circulated between and along the 1.2 m length of the plates past ice crystals which nucleated and grew from a fiber suspended vertically between the two plates. A zoom stereo microscope with a magnification which ranged from 3X to 80X and both 35 mm still photographs and 16 mm time lapse cine films taken through the microscope were used to study the variation of the shape and linear growth rate of ice crystals as a function of the ambient temperature, the ambient supersaturation, and the forced ventilation velocity. The ambient growth conditions were varied over the range of temperature 0 to -40 C, over the range of supersaturation 4% to 50% with respect to ice, and over the range of forced ventilation velocities 0 cm/s to 20 cm/s.

Keller, V. W.↗
NASA NTRS1981

Analytical description of charged particles transport along arbitrary guiding-field configurations

A new description of focused transport was developed which is valid for arbitrary spatial dependences of both the scattering mean free path and the focusing length. In particular, it describes the supercoherent transition from the coherent modes that occurs when focusing is strong compared to scattering to the diffusive mode that occurs when focusing is weak. Because it specifies the pitch angle distribution in terms of relatively simple functions, the new formulation facilitates not only the detailed interpretation of solar particle events, but also the qualitative understanding of focused transport. In the steady state, it predicts that the angular distribution consists of a highly collimated 'strahl' which is superimposed on an anisotropic 'halo'. This picture provides new insight into observed configurations of steady-state particle fluxes.

Earl, J. A.↗
DOE OSTI2020

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1981

Application of aerial photography to the study of small scale upper ocean phenomena

The industrial waste dumped 180 n. miles south of Galveston was monitored in July 1977 by water sampling, hydrographic measurements, acoustic tracking on board two vessels, and by aerial photography. The plume of the waste diffused vertically and horizontally. Photodensitometry of aerial photos of the plume showed lateral dispersion of the plume in agreement with two other methods: acoustic tracking of the waste suspensoid and transmissometer sampling. In addition, the method showed small scale features like the lateral and longitudinal variations in the photodensity, indicating the waste concentration. This waste concentration showed periodic changes in its axial distance, with the spectral peak at about 160 m wave length. It shows a sharp increase at the windward edge of the plume as do the acoustic records. This phenomenon is explained in terms of the shearing current near the surface together with vertical diffusion. The periodic change along the axis is explained in terms of the Langmuir circulation and in terms of internal ship waves.

Ichiye, T.↗
DOE OSTI2025

Extending the capillary wave model to include the effect of bending rigidity: X-ray reflectivity and diffuse scattering

The surface roughness of a thin film at a liquid interface exhibits contributions of thermally excited fluctuations. This thermal roughness depends on temperature (𝑇), surface tension (𝛾), and elastic material properties, specifically the bending modulus (𝜅) of the film. A nonzero 𝜅 suppresses the thermal roughness at small length scales compared to an interface with zero 𝜅, as expressed by the power spectral density of the thermal roughness. The description of the x-ray scattering of the standard capillary wave model (CWM), which is valid for zero 𝜅, is extended to include the effect of 𝜅. The extended CWM (eCWM) provides a single analytical form for both the specular x-ray reflectivity (XRR) and the diffuse scattering around the specular reflection, and recovers the expression of the CWM at its zero 𝜅 limit. This theoretical approach enables the use of single-shot grazing incidence x-ray off-specular scattering (GIXOS) measurements for characterizing the structure of thin films on a liquid surface. The eCWM analysis approach decouples the thermal roughness factor from the surface scattering signal, providing direct access to the intrinsic surface-normal structure of the film and its bending modulus. Moreover, the eCWM facilitates the calculation of reflectivity at any desired resolution (pseudo-XRR approach). The transformation into pseudo-XRR provides the benefit of using widely available XRR software to perform GIXOS analysis. The extended range of the vertical scattering vector (𝑄 𝑧 ) available with the GIXOS pseudo-XRR approach allows for a higher spatial resolution than with conventional XRR. Experimental results are presented for various lipid systems, showing strong agreement between conventional specular XRR and pseudo-XRR methods. This agreement validates the proposed approach and highlights its utility for analyzing soft, thin films.

36 MATERIALS SCIENCE↗
NASA NTRS1993

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗
DOE OSTI2024

Uncovering fast solid-acid proton conductors based on dynamics of polyanion groups and proton bonding strength

Achieving high proton conductivity in inorganic solids is key for advancing many electrochemical technologies, including low-energy nano-electronics and energy-efficient fuel cells and electrolyzers. A quantitative understanding of the physical traits of a material that regulate proton diffusion is necessary for accelerating the discovery of fast proton conductors. In this work, we have mapped the structural, chemical and dynamic properties of solid acids to the elementary steps of the Grotthuss mechanism of proton diffusion. Our approach combines ab initio molecular dynamics simulations, analysis of phonon spectra and atomic structure calculations. We have identified the donor–hydrogen bond lengths and the acidity of polyanion groups as key descriptors of local proton transfer and the vibrational frequencies of the cation framework as the key descriptor of lattice flexibility. The latter facilitates rotations of polyanion groups and long-range proton migration in solid acid proton conductors. The calculated lattice flexibility also correlates with the experimentally reported superprotonic transition temperatures. Using these descriptors, we have screened the Materials Project database and identified potential solid acid proton conductors with monovalent, divalent and trivalent cations, including Ag + , Sr 2+ , Ba 2+ and Er 3+ cations, which go beyond the traditionally considered monovalent alkali cations (Cs + , Rb + , K + , and NH 4 + ) in solid acids.

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

Effect of boattail geometry on the acoustics of parallel baffles in ducts

Sound attenuation and total pressure drop of parallel duct baffles incorporating certain boattail geometries were measured in the NASA Ames Research Center 7- by 10-Foot Wind Tunnel. The baseline baffles were 1.56 m long and 20 cm thick, on 45-cm center-to-center spacings, and spanned the test section from floor to ceiling. Four different boattails were evaluated: a short, smooth (nonacoustic) boattail; a longer, smooth boattail; and two boattails with perforated surfaces and sound-absorbent filler. Acoustic measurements showed the acoustic boattails improved the sound attenuation of the baffles at approximately half the rate to be expected from constant-thickness sections of the same length; that is, 1.5 dB/n, where n is the ratio of acoustic treatment length to duct passage width between baffles. The aerodynamic total pressure loss was somewhat sensitive to tail geometry. Lengthening the tails to reduce the diffusion half-angle from 11 to 5 degrees reduced the total pressure loss approximately 9%. Perforating the boattails, which increased the surface roughness, did not have a large effect on the total pressure loss. Aerodynamic results are compared with a published empirical method for predicting baffle total pressure drop.

Soderman, P. T.↗
NASA NTRS1984

The effect of convection upon charged particle transport in random magnetic fields

In a coordinate system moving with the plasma and random magnetic fields of a wind that blows with constant velocity in the direction of the guiding field, transport of energetic particles is described by a Boltzmann equation which is similar to the one that describes unconvected transport. Although this formulation is mathematically identical to that developed by Luhmann, which refers to the system where the guiding field is static, there are both practical and fundamental reasons to adopt the new approach. It leads to first-order approximate transport equations which are similar to those that apply in the absence of convection. However, these equations are more general than Parker's description of diffusion and convection, for they describe the coherent modes of transport that appear when the mean free path is large compared to the scale length for spatial variations of the guiding field, and they are valid for arbitrary wind velocity. The latter characteristic opens up new possibilities for analyzing particle transport in relativistic flows seen in some astronomical objects.

Earl, J. A.↗
NASA NTRS1992

Hit rates and radiation doses to nuclei of bone lining cells from alpha-particle-emitting radionuclides

Factors relating the local concentration of a bone-seeking alpha-particle emitter to the mean hit rate have been determined for nuclei of bone lining cells using a Monte Carlo procedure. Cell nuclei were approximated by oblate spheroids with dimensions and location taken from a previous histomorphometric study. The Monte Carlo simulation is applicable for planar and diffuse labels at plane or cylindrical bone surfaces. Additionally, the mean nuclear dose per hit, the dose mean per hit, the mean track segment length and its second moment, the percentage of stoppers, and the frequency distribution of the dose have been determined. Some basic features of the hit statistics for bone lining cells have been outlined, and the consequences of existing standards of radiation protection with regard to the hit frequency to cell nuclei are discussed.

NASA Discipline Musculoskeletal↗
NASA NTRS2004

[Using Molecular Simulations to Understand Complex Nanoscale Dynamic Phenomena in Polymer Solutions]

The first half of the project concentrated on molecular simulation studies of the translocation of model molecules for single-stranded DNA through a nanosized pore. This has resulted in the publication, Translocation of a polymer chain across a nanopore: A Brownian dynamics simulation study, by Pu Tian and Grant D. Smith, JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 21 1 DECEMBER 2003, which is attached to this report. In this work we carried out Brownian dynamics simulation studies of the translocation of single polymer chains across a nanosized pore under the driving of an applied field (chemical potential gradient) designed to mimic an electrostatic field. The translocation process can be either dominated by the entropic barrier resulted from restricted motion of flexible polymer chains or by applied forces (or chemical gradient). We focused on the latter case in our studies. Calculation of radius of gyration of the translocating chain at the two opposite sides of the wall shows that the polymer chains are not in equilibrium during the translocation process. Despite this fact, our results show that the one-dimensional diffusion and the nucleation model provide an excellent description of the dependence of average translocation time on the chemical potential gradients, the polymer chain length and the solvent viscosity. In good agreement with experimental results and theoretical predictions, the translocation time distribution of our simple model shows strong non-Gaussian characteristics. It is observed that even for this simple tube-like pore geometry, more than one peak of translocation time distribution can be generated for proper pore diameter and applied field strengths. Both repulsive Weeks-Chandler-Anderson and attractive Lennard-Jones polymer-nanopore interaction were studied. Attraction facilitates the translocation process by shortening the total translocation time and dramatically improve the capturing of polymer chain. The width of the translocation time distribution was found to decrease with increasing temperature, increasing field strength, and decreasing pore diameter.

Smith, Grant↗
NASA NTRS2003

Translocation of a Polymer Chain across a Nanopore: A Brownian Dynamics Simulation Study

We carried out Brownian dynamics simulation studies of the translocation of single polymer chains across a nanosized pore under the driving of an applied field (chemical potential gradient). The translocation process can be either dominated by the entropic barrier resulted from restricted motion of flexible polymer chains or by applied forces (or chemical gradient across the wall), we focused on the latter case in our studies. Calculation of radius of gyrations at the two opposite sides of the wall shows that the polymer chains are not in equilibrium during the translocation process. Despite this fact, our results show that the one-dimensional diffusion and the nucleation model provide an excellent description of the dependence of average translocation time on the chemical potential gradients, the polymer chain length and the solvent viscosity. In good agreement with experimental results and theoretical predictions, the translocation time distribution of our simple model shows strong non-Gaussian characteristics. It is observed that even for this simple tubelike pore geometry, more than one peak of translocation time distribution can be generated for proper pore diameter and applied field strengths. Both repulsive Weeks-Chandler-Anderson and attractive Lennard-Jones polymer-nanopore interaction were studied, attraction facilitates the translocation process by shortening the total translocation time and dramatically improve the capturing of polymer chain. The width of the translocation time distribution was found to decrease with increasing temperature, increasing field strength, and decreasing pore diameter.

Tian, Pu↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2020

Competition between Hückel’s Rule and Jahn–Teller Distortion in Small Carbon Rings: A Quantum Monte Carlo Study

Structural properties and energetics of carbon rings are studied with the diffusion Monte Carlo (DMC) method. Our DMC-based geometry optimization reveals that both polyynic C 4n and cumulenic C 4n+2 rings exhibit bond length alternations for n >= 3, which is understood to be due to Jahn-Teller distortions. The bond length alternation even in a cumulenic (4n + 2) carbon ring was experimentally observed in a recently synthesized C 18 molecule. From a comparison of the DMC cohesive energies of C 4n with those of C 4n + 2 , we present a comprehensive picture of the competition between Huckel's rule and Jahn-Teller distortion in small carbon rings; the former is more dominant than the latter for n < 5 where C 4n + 2 rings are more stable than C 4n , while C 4n rings are as stable as C 4n + 2 for n < 5 where dimerization effects due to Jahn-Teller distortion are more important.

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

Accelerating network layouts using graph neural networks

Graph layout algorithms used in network visualization represent the first and the most widely used tool to unveil the inner structure and the behavior of complex networks. Current network visualization software relies on the force-directed layout (FDL) algorithm, whose high computational complexity makes the visualization of large real networks computationally prohibitive and traps large graphs into high energy configurations, resulting in hard-to-interpret “hairball” layouts. Here we use Graph Neural Networks (GNN) to accelerate FDL, showing that deep learning can address both limitations of FDL: it offers a 10 to 100 fold improvement in speed while also yielding layouts which are more informative. We analytically derive the speedup offered by GNN, relating it to the number of outliers in the eigenspectrum of the adjacency matrix, predicting that GNNs are particularly effective for networks with communities and local regularities. Finally, we use GNN to generate a three-dimensional layout of the Internet, and introduce additional measures to assess the layout quality and its interpretability, exploring the algorithm’s ability to separate communities and the link-length distribution. The novel use of deep neural networks can help accelerate other network-based optimization problems as well, with applications from reaction-diffusion systems to epidemics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗