MFIX-Exa: CFD-DEM Simulations of Thermodynamics and Chemical Reactions in Multiphase Flows
For presentation at the AIChE Annual Meeting, Orlando, FL, November 5-10, 2023.
Porcu, Roberto↗
SEARCH · Search NASA
Results for “scale analysis”
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.
Discussion on molecular dynamics (MD) simulations of the asphalt materials
The application of asphalt materials in pavement engineering has been increasingly widespread and sophisticated over the past several decades. Variations in the properties of asphalt binder during mixing, transportation, and paving can affect the performance of asphalt pavement. However, the asphalt material is a non-homogeneous and complex organic substance, consisting of various molecules with widely various molecular weights, elemental compositions, and structures. This complexity leads to difficulties for researchers to clearly and immediately understand the properties of asphalt materials and their variations. The multi-scale research approach combines macroscopic experimental data and microscopic simulation results from a practical engineering perspective. It helps to improve the understanding of asphalt materials. The molecular dynamics (MD) simulation proposes a corresponding molecular model of asphalt material based on experimental data, and the simulation algorithm is able to derive properties similar to those of real asphalt. Here, this paper provides a comprehensive review of the current studies on MD simulation of asphalt materials, including modeling, properties, and multi-scale analysis. As a key part of the computational simulation, this paper discusses the typical asphalt binder and asphalt-aggregate interface models constructed by different groups, and also presents their differences from real samples and their feasibility based on fundamental properties. After the introduction of molecular models, the extensive work made by researchers based on molecular models is categorically reviewed and discussed. The strengths and weaknesses of MD simulation methods in the study of asphalt materials are also summarized in order to provide the reader with a more comprehensive understanding of the relevant contents and to guide subsequent research.
36 MATERIALS SCIENCE↗
Uniform secondary current distribution at disk electrodes under Tafel kinetics enabled by concentric current-shielding rings
Non-uniform secondary current distribution at rotating disk electrodes (RDE) is a common problem when using resistive electrolyte media or large applied currents. In a recent publication, we have shown that auxiliary electrodes such as the ring of a rotating ring disk electrode (RRDE) can help suppress current non-uniformities at the disk enabling reliable electroanalytical measurements. However, this previous work considered linear kinetics where current distribution non-uniformities were moderate. In the present contribution, we consider current distribution non-uniformities encountered under Tafel kinetics. We show, for the case of 2,5-dihydroxy-1,4-benzoquinone (DHBQ) reduction, that optimally chosen ring conditions serve to provide effective shielding at the disk edge rendering the overall disk current distribution to be uniform. Numerical modeling and scaling analysis (using the Wagner number) are presented to aid a user in determining the optimal ring current density for achieving uniform disk current distribution under Tafel kinetics. This approach is especially useful when studying soluble-soluble redox transitions for which, unlike deposit distribution in electrodeposition, the current distribution non-uniformity is not visually apparent.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Gas hydrate deposit formation in transient flowloop tests and mitigation with a surface treatment
Thermodynamic and thermal mitigation strategies are traditionally utilized to prevent gas hydrate formation but may be rendered ineffective for certain unplanned, extended shut-in scenarios. Unfortunately, the current knowledge regarding transient gas hydrate formation has been extremely limited to only a few studies and unable to address this complex issue adequately. This work, conducted in a lab-scale flowloop, aims to add to the understanding of this area by delving into the implications of transient conditions on gas hydrate formation mechanisms, while also exploring a possible mitigation strategy of a surface treatment. Specifically, the effects of phase separation on gas hydrate formation locations, transportability, and plugging are presented for an oil-dominated system. An omniphobic surface treatment is shown to affect initial gas hydrate deposition and its wall adherence strength, as well as the overall plugging outcome. Overall, this work combines observational information with process data measurements to provide a better understanding of phase separated gas hydrate systems and passive deposition mitigation strategies. This can further aid advanced gas hydrate modeling efforts, which currently do not adequately address transient flow effects on phase separation/dispersion. Here, it offers representative gas hydrate testing for coatings, which have exclusively been limited to non-flowing and/or bench-scale analysis, to give confidence in their potential application as advanced deposition prevention tools.
03 NATURAL GAS↗
Unified descriptor for enhanced critical heat flux during pool boiling of hemi-wicking surfaces
Boiling heat transfer is dictated by interfacial phenomena at the three-phase contact line where vapor bubbles form on the surface. Structured surfaces have shown significant enhancement in critical heat flux (CHF) during pool boiling by tailoring interfacial phenomena. This CHF enhancement has been primarily explained by two structural effects: roughness, which extends the contact line length at the bubble base, and wickability, the ability to imbibe liquid through surface structures by capillary pumping. Here, in this work, we show that CHF enhancement on structured surfaces cannot be described by roughness or wickability alone. This result was confirmed using systematically designed micropillar surfaces with controlled roughness and wickability. Further, we performed a scaling analysis and derived a unified descriptor, which represents the combined effects of thin film density and volumetric wicking rate. This unified descriptor shows a reasonable correlation with CHF values with our experiments and literature data. This work provides important insights in understanding the role of surface structures on CHF enhancement, thereby providing guidelines for the systematic design of surface structures for enhanced pool boiling heat transfer.
42 ENGINEERING↗
Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling
The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Gas-driven short disconnection mitigates thermal runaway in Li-ion batteries under mechanical abuse
Mechanical abuse poses a critical safety risk to Li-ion batteries by inducing internal short circuits that initiate thermal runaway. Remarkably, voltage recovery frequently emerges during thermal runaway initiation, a phenomenon that conventional theories fail to explain. Here, our study develops a new multiphysics mechanism to explain voltage recovery, termed gas-driven short disconnection, whereby internal gas pressure mechanically disengages short-circuit contacts and causes the cell voltage to rebound. This mechanism incorporates gas generation and its structural impact on the short circuit. Real-time optical and thermal imaging and X-ray computed tomography reveal fluid-structure interaction between internal gas flow and adjacent shorting contacts. We establish a mechanistic framework linking gas-driven short disconnection to cell-level voltage and temperature responses, elucidating the extension–truncation pattern of voltage recovery. Furthermore, thermal regime maps show that a voltage recovery duration exceeding 5 s correlates with limited temperature rise below 150 °C, indicating that sustained short-circuit disconnection suppresses Joule heating. Additionally, a dimensionless criterion is deduced from scaling analysis for physical plausibility of gas-driven short disconnection in mechanically abused cells. This finding inspires smart venting control, which regulates gas release to maintain the internal pressure while dissipating gas enthalpy, thereby providing a device-level strategy for thermal runaway mitigation.
Gas-driven short disconnection↗
The role of geographical spreaders in infectious pattern formation and front propagation speeds
The pattern formation and spatial spread of infectious populations are investigated using a kernel-based Susceptible–Infectious–Recovered (SIR) model applicable across a wide range of basic reproduction numbers R O . The goal is to examine the role of geographical spreaders on transient spatial pattern formation of infectious populations and the associated maximum invasive front speeds c max . In the simulations conducted here, geographical spreaders are defined as a portion of the infected population Φ experiencing high mobility between identical communities. The spatial organization of the infected population and c max are determined when the infections are randomly initiated in space within multiple communities. For small but finite , scaling analysis and numerical simulations in 1-dimension suggest that when the spreading kernel is Gaussian-shaped, where is the inverse of the infectious duration. This finding for agrees with a diffusion-based representation of mobility in 1-D. Numerical simulations in 2-D across wide-ranging suggest that , the variance of the spatial kernel describing mobility of long-distance geographical spreaders across communities, determines the spatial organization of infections across communities. When (long-distance mobility, where is the minimum spatial extent defining adjacent communities), the infectious population will experience a transient but spatially coherent pattern with a wavelength that can be derived from the spreading kernel properties. Moreover, the 2-D simulations for the bounded kernel suggest that attainment of is also dictated by but the magnitude is not sensitive to unlike diffusion-based models.
60 APPLIED LIFE SCIENCES↗
Integrated Systems-to-Atoms (S2A) Framework for Designing Resilient and Efficient Hydrogen Infrastructure Solutions
The success of a future clean hydrogen infrastructure will depend on technology performance, operating conditions, and system configuration, which will be integrated to meet specific end-use requirements. We introduce a Systems-to-Atoms framework to cosimulate material, device, and system design for clean hydrogen storage and transport. Here we present a demonstration scenario in which hydrogen is distributed using a liquid organic hydrogen carrier to 700 bar refueling stations. Cross-scale analysis on hydrogen release at the refueling station reveals that while a high reaction pressure (up to 100 bar) may not enhance the catalyst performance, it can lead to system-level cost savings (up to 5% of total refueling station cost) due to savings in downstream gas compression. Further, a typically low-performing dehydrogenation catalyst material (e.g., copper) may be preferred over a high-performing catalyst (e.g., palladium) at different operating conditions when considering overall system affordability and supply chain resilience.
08 HYDROGEN↗
High Magnetic Anisotropy and Magnetocaloric Effects in Single-Crystal Cr 2 Te 3
Here, we report a systematic investigation of anisotropic magnetocaloric effects in single-crystal Cr 2 Te 3 . Single-crystal samples are synthesized by chemical vapor transport and characterized by X-ray and Laue diffraction methods. The maximum magnetic entropy change –ΔS M max is 4.50 J kg –1 K –1 for the easy c-axis (3.36 J kg –1 K –1 for the hard axis along ab-plane), and the relative cooling power (RCP) is 296.7 J kg –1 for the easy c-axis (183.84 J kg –1 for the hard axis along ab-plane) for a magnetic field change of 9 T near the Curie temperature. The magneto-crystalline anisotropy constant K u is estimated to be 580.12 kJ m –3 at 140 K, decreasing to 148.60 kJ m –3 at 168 K. Meanwhile, the maximum of the rotational magnetic entropy change –ΔS M R (T, H) between the c-axis and the ab-plane is about 1.14 J kg –1 K –1 for magnetic-field change of 9 T. The critical exponents are estimated by analyzing magnetocaloric effects, which indicate a 2D-Ising type magnetic system. The accuracy of estimated critical exponents is verified by scaling analysis. The maximum magnetic entropy change –ΔS M max ≈ 5.25 J kg –1 K –1 (along the c-axis) and the corresponding adiabatic temperature change ΔT ad ≈ 3.31 K (along the c-axis) are estimated by analyzing heat capacity measurements with a magnetic field up to 9 T.
36 MATERIALS SCIENCE↗
Screening Cu-Zeolites for Methane Activation Using Curriculum-Based Training
Machine learning (ML), when used synergistically with atomistic simulations, has recently emerged as a powerful tool for accelerated catalyst discovery. However, the application of these techniques has been limited by the lack of interpretable and transferable ML models. In this work, we propose a curriculum-based training (CBT) philosophy to systematically develop reactive machine learning potentials (rMLPs) for high-throughput screening of zeolite catalysts. Our CBT approach combines several different types of calculations to gradually teach the ML model about the relevant regions of the reactive potential energy surface. The resulting rMLPs are accurate, transferable, and interpretable. We further demonstrate the effectiveness of this approach by exhaustively screening thousands of [CuOCu] 2+ sites across hundreds of Cu-zeolites for the industrially relevant methane activation reaction. Specifically, this large-scale analysis of the entire International Zeolite Association (IZA) database identifies a set of previously unexplored zeolites (i.e., MEI, ATN, EWO, and CAS) that show the highest ensemble-averaged rates for [CuOCu] 2+ -catalyzed methane activation. We believe that this CBT philosophy can be generally applied to other zeolite-catalyzed reactions and, subsequently, to other types of heterogeneous catalysts. Thus, this represents an important step toward overcoming the long-standing barriers within the computational heterogeneous catalysis community.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The breakdown of both strange metal and superconducting states at a pressure-induced quantum critical point in iron-pnictide superconductors
Here we report the first observation of the concurrent breakdown of the strange metal (SM) normal state and superconductivity at a pressure-induced quantum critical point in Ca 10 (Pt 4 As 8 )((Fe 0.97 Pt 0.03 ) 2 As 2 ) 5 superconductor. We find that, upon suppressing the superconducting state, the power exponent (α) changes from 1 to 2, and the slope of the temperature-linear resistivity per FeAs layer (A $\Box$ ) gradually diminishes. At a critical pressure, A $\Box$ and superconducting transition temperature (T c ) go to zero concurrently, where a quantum phase transition from a superconducting state with a SM normal state to a non-superconducting Fermi liquid state occurs. Scaling analysis reveals that the change of A $\Box$ with T c obeys the relation of T c ~ (A $\Box$ ) 0.5 , similar to what is seen in other chemically doped unconventional superconductors. These results suggest that there is a simple but powerful organizational principle of connecting the SM normal state with the high-T c superconductivity.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Surface triggered stabilization of metastable charge-ordered phase in SrTiO 3
Charge ordering (CO), characterized by a periodic modulation of electron density and lattice distortion, has been a fundamental topic in condensed matter physics, serving as a potential platform for inducing novel functional properties. The charge-ordered phase is known to occur in a doped system with high d-electron occupancy, rather than low occupancy. Here, we report the realization of the charge-ordered phase in electron-doped (100) SrTiO 3 epitaxial thin films that have the lowest d-electron occupancy i.e., d 1 -d 0 . Theoretical calculation predicts the presence of a metastable CO state in the bulk state of electron-doped SrTiO 3 . Atomic scale analysis reveals that (100) surface distortion favors electron-lattice coupling for the charge-ordered state, and triggering the stabilization of the CO phase from a correlated metal state. This stabilization extends up to six unit cells from the top surface to the interior. Our approach offers an insight into the means of stabilizing a new phase of matter, extending CO phase to the lowest electron occupancy and encompassing a wide range of 3d transition metal oxides.
36 MATERIALS SCIENCE↗
Closed-form theoretical model of the secondary drop size in partial coalescence—Capturing pertinent timescales and viscous forces
One of the most important outcomes of partial drop coalescence is the ratio of the secondary drop radius to the primary drop radius, known as the drop ratio, ri. This ratio is thought to be approximately constant and independent of physical parameters of the fluids involved. However, here we reveal that ambient fluid viscosity can alter the size of the secondary drop and the drop ratio consequently. Using scaling analysis, we derive a model that predicts the behavior of the drop ratio as a function of the Ohnesorge number, a dimensionless ratio of viscous to inertial forces. In addition, we present our experimental results of coalescing drops on a planar interface under the influence of surface tension gradients. A high-speed digital camera is used to observe the evolution of drops as they coalesce with a bulk liquid. We show that this process is influenced by the surface tension gradient between the drop and the bulk liquid. The ratio of the secondary drop to the primary drop in partial coalescence is smaller than the reported values for coalescence without a surface tension gradient. The analytical model derived through this study is based on a new modified Ohnesorge number that includes surface tension gradients. Our analytical model is compared against other models and the results illustrate good agreement with our experimental findings and experimental data in the literature.
42 ENGINEERING↗
Scattering evidence of positional charge correlations in polyelectrolyte complexes
Polyelectrolyte complexation plays an important role in materials science and biology. The internal structure of the resultant polyelectrolyte complex (PEC) phase dictates properties such as physical state, response to external stimuli, and dynamics. Small-angle scattering experiments with X-rays and neutrons have revealed structural similarities between PECs and semidilute solutions of neutral polymers, where the total scattering function exhibits an Ornstein–Zernike form. In spite of consensus among different theoretical predictions, the existence of positional correlations between polyanion and polycation charges has not been confirmed experimentally. Here, we present small-angle neutron scattering profiles where the polycation scattering length density is matched to that of the solvent to extract positional correlations among anionic monomers. The polyanion scattering functions exhibit a peak at the inverse polymer screening radius of Coulomb interactions, q* ≈ 0.2 Å –1 . This peak, attributed to Coulomb repulsions between the fragments of polyanions and their attractions to polycations, is even more pronounced in the calculated charge scattering function that quantifies positional correlations of all polymer charges within the PEC. Screening of electrostatic interactions by adding salt leads to the gradual disappearance of this correlation peak, and the scattering functions regain an Ornstein–Zernike form. Experimental scattering results are consistent with those calculated from the random phase approximation, a scaling analysis, and molecular simulations.
74 ATOMIC AND MOLECULAR PHYSICS↗
An implementation of neural simulation-based inference for parameter estimation in ATLAS
Neural simulation-based inference (NSBI) is a powerful class of machine-learning-based methods for statistical inference that naturally handles high-dimensional parameter estimation without the need to bin data into low-dimensional summary histograms. Such methods are promising for a range of measurements, including at the Large Hadron Collider, where no single observable may be optimal to scan over the entire theoretical phase space under consideration, or where binning data into histograms could result in a loss of sensitivity. This work develops a NSBI framework for statistical inference, using neural networks to estimate probability density ratios, which enables the application to a full-scale analysis. It incorporates a large number of systematic uncertainties, quantifies the uncertainty due to the finite number of events in training samples, develops a method to construct confidence intervals, and demonstrates a series of intermediate diagnostic checks that can be performed to validate the robustness of the method. As an example, the power and feasibility of the method are assessed on simulated data for a simplified version of an off-shell Higgs boson couplings measurement in the four-lepton final states. This approach represents an extension to the standard statistical methodology used by the experiments at the Large Hadron Collider, and can benefit many physics analyses.
frequentist statistics↗
Manipulating high-temperature superconductivity by oxygen doping in Bi2Sr2CaCu2O8+δ thin flakes
Abstract Harnessing the fascinating properties of correlated oxides requires precise control of their carrier density. Compared to other methods, oxygen doping provides an effective and more direct way to tune the electronic properties of correlated oxides. Although several approaches, such as thermal annealing and oxygen migration, have been introduced to change the oxygen content, a continuous and reversible solution that can be integrated with modern electronic technology is much in demand. Here, we report a novel ionic field-effect transistor using solid Gd-doped CeO2 as the gate dielectric, which shows a remarkable carrier-density-tuning ability via electric-field-controlled oxygen concentration at room temperature. In Bi2Sr2CaCu2O8+δ (Bi-2212) thin flakes, we achieve a reversible superconductor–insulator transition by driving oxygen ions in and out of the samples with electric fields, and map out the phase diagram all the way from the insulating regime to the over-doped superconducting regime by continuously changing the oxygen doping level. Scaling analysis indicates that the reversible superconductor–insulator transition for the Bi-2212 thin flakes follows the theoretical description of a two-dimensional quantum phase transition. Our work provides a route for realizing electric-field control of phase transition in correlated oxides. Moreover, the configuration of this type of transistor makes heterostructure/interface engineering possible, thus having the potential to serve as the next-generation all-solid-state field-effect transistor.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗