Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscous effects”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Size Ranges of Magnetic Domain States in Tetrataenite

Paleomagnetic studies of meteorites provide unique constraints on the evolution of magnetic fields in the early solar system. These studies rely on the identification of magnetic minerals that can retain stable magnetizations over ≳4.5 billion years (Ga). The ferromagnetic mineral tetrataenite (γ''-Fe 0.5 Ni 0.5 ) is found in iron, stony-iron and chondrite meteorite groups. Nanoscale intergrowths of tetrataenite have been shown to carry records of paleomagnetic fields, although the effect of magnetostatic interactions on their magnetic remanence acquisition remains to be fully understood. Tetrataenite can also occur as isolated, non-interacting, nanoscale grains in many meteorite groups, although the paleomagnetic potential of these grains is particularly poorly understood. Here, we aim to improve our understanding of tetrataenite magnetization to refine our knowledge of existing paleomagnetic analyses and broaden the spectrum of meteorite groups that can be used for future paleomagnetic studies. We present the results of analytical calculations and micromagnetic modeling of isolated tetrataenite grains with various geometries. We find that tetrataenite forms a stable single domain state at grain lengths between 6 and ~160 nm dependent on its elongation. It also possesses a magnetization resistant to viscous remagnetization over the lifetime of the solar system at 293 K. At larger grain sizes, tetrataenite's lowest energy state is a lamellar two-domain state, stable at Ga-scale timescales. Unlike many other magnetic minerals, tetrataenite does not form a single-vortex domain state due to its large uniaxial anisotropy. Our results show that single domain and two-domain tetrataenite grains carry an extremely stable magnetization and therefore are promising for paleomagnetic studies.

58 GEOSCIENCES↗
DOE OSTI2024

Extracting the speed of sound in quark–gluon plasma with ultrarelativistic lead–lead collisions at the LHC

Abstract Ultrarelativistic nuclear collisions create a strongly interacting state of hot and dense quark–gluon matter that exhibits a remarkable collective flow behavior with minimal viscous dissipation. To gain deeper insights into its intrinsic nature and fundamental degrees of freedom, we determine the speed of sound in an extended volume of quark–gluon plasma using lead–lead (PbPb) collisions at a center-of-mass energy per nucleon pair of 5.02 TeV. The data were recorded by the CMS experiment at the CERN LHC and correspond to an integrated luminosity of 0.607 nb −1 . The measurement is performed by studying the multiplicity dependence of the average transverse momentum of charged particles emitted in head-on PbPb collisions. Our findings reveal that the speed of sound in this matter is nearly half the speed of light, with a squared value of 0.241 ± 0.002 ( stat ) ± 0.016 ( syst ) in natural units. The effective medium temperature, estimated using the mean transverse momentum, is 219 ± 8 ( syst ) MeV . The measured squared speed of sound at this temperature aligns precisely with predictions from lattice quantum chromodynamic (QCD) calculations. This result provides a stringent constraint on the equation of state of the created medium and direct evidence for a deconfined QCD phase being attained in relativistic nuclear collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Stacked networks improve physics-informed training: Applications to neural networks and deep operator networks

Physics-informed neural networks and operator networks have shown promise for effectively solving equations modeling physical systems. However, these networks can happen to be difficult or impossible to train accurately. Here, we present a novel multifidelity framework for stacking physics-informed neural networks and operator networks that facilitates training. We successively build a chain of networks, where the output at one step can act as a low-fidelity input for training a longer chain, gradually increasing the expressivity of the learnt model. The equations imposed at each step of the iterative process can be the same or different (akin to simulated annealing). The iterative (stacking) nature of the proposed method allows us to learn progressively features of a solution which could have been hard to learn directly. Through benchmark problems including a nonlinear pendulum, the wave equation, and the viscous Burgers equation, we show how stacking can be used to improve the accuracy and reduce the required size of physics-informed neural networks and operator networks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2022

Quantifying transport and electrocatalytic reaction processes in a gastight rotating cylinder electrode reactor via integration of Computational Fluid Dynamics modeling and experiments

Understanding the complexity of the multiple processes of mass, momentum, charge, and heat transport, and how these affect reaction kinetics at the electrode/electrolyte interface is one of the major challenges in the field of energy and catalysis. The rapid and rational scale-up of electrocatalytic systems to industrial scales require a detailed understanding of nonlinear transport-reaction processes, accessible only through the building of multi-physics models that capture with high fidelity the complexity of real-world devices. The gastight rotating cylinder electrode (RCE) reactor is a promising lab-scale tool that can decouple transport from intrinsic kinetics to generate data for first-principle models useful in the design of industrial, electrochemical reactors. Computational Fluid Dynamics (CFD) studies have previously been used to investigate the bulk flow in RCE reactors for simple corrosion and electroplating processes. However, the quantification of changes in local concentration within the viscous layer where catalysis takes place requires capturing the correct flow conditions inside the hydrodynamic boundary layer near the surface of the electrode. Further, this requires simulations with spatial resolution in the nm and μm scale and temporal resolutions between ms and s scales that are similar to the timescales for reactions on the electrode surface. In this study, experimental electrocatalysis is combined with CFD modeling to elucidate and parameterize the hydrodynamics in a gastight RCE reactor. CFD simulations of the electrochemical ferricyanide reduction reaction under mass transport limited conditions are used to evaluate the validity of the CFD model parameters by comparing calculated dimensionless mass transport descriptors to dimensionless correlations obtained experimentally. Justifications for assumptions and details of the simulation methods used in this study are presented to provide a detailed understanding of the effect that each model parameter has on the ability to accurately simulate electrocatalysis in RCE systems. The simulation methodology reported here is a first step towards the development of multi-scale models for the study of transport dependent electrocatalytic processes, such as the electrochemical transformation of CO 2 to fuels and chemicals.

42 ENGINEERING↗
DOE OSTI2025

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

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

Utilization of event shape in search of the chiral magnetic effect in heavy-ion collisions

The search for the chiral magnetic effect (CME) has been a subject of great interest in the field of high-energy heavy-ion collision physics, and various observables have been proposed to probe the CME. Experimental observables are often contaminated with background contributions arising from collective motions (specifically elliptic flow) of the collision system. We present a method study of event-shape engineering (ESE) that projects the CME-sensitive γ 112 correlator and its variations (γ 132 and γ 123 ) to a class of events with minimal flow. We discuss the realization of the zero-flow mode, the sensitivity on the CME signal, and the corresponding statistical significance for Au + Au, Ru + Ru, and Zr + Zr collisions at √s NN = 200 GeV with a multiphase transport model, as well as a new event generator, event-by-event anomalous-viscous fluid dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Rotational and translational diffusion of liquid n-hexane: EFP-based molecular dynamics analysis

Molecular Dynamics (MD) simulations based on the Effective Fragment Potential (EFP) method are utilized to provide a comprehensive assessment of diffusion in liquid n-hexane. We decompose translational diffusion into components along and orthogonal to the long axis of the molecule. Rotational diffusion is decomposed into tumbling and spinning motions about this axis. Our analysis yields four corresponding diffusion coefficients which are related to diagonal entries in the complete 6 × 6 diffusion tensor accounting for the three rotational and three translational degrees of freedom and for the potential coupling between them. However, coupling between different degrees of freedom is expected to be minimal for a natural choice of the molecular body-fixed axis, so then off-diagonal entries in the tensor are negligible. This expectation is supported by a hydrodynamic analysis of the diffusion tensor which treats the liquid surrounding the molecule being tracked as a viscous continuum. Thus, the EFP MD analysis provides a comprehensive characterization of diffusion and also reveals expected shortcomings of the hydrodynamic treatment, particularly for rotational diffusion, when applied to neat liquids.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Retention of Long-period Gas Giant Planets: Type II Migration Revisited

During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

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

Rates of Sea‐Level Rise Are Highly Sensitive to Ice Viscosity Parameters in Model Benchmarks

Glacier flow plays a major role in current and future rates of globally averaged sea-level rise. The viscosity of glacial ice, controlling the rate of flow, decreases as stress increases and is highly sensitive to the value of the stress exponent, $n$, in the constitutive equation for viscous flow. Glaciologists and climate modelers almost exclusively assume $n=3$ when modeling ice flow and projecting sea-level rise through forward modeling. However, recent work suggests that $n\approx 4$ better fits observations, prompting the question: How sensitive are projections of sea-level rise to the value of $n$? We use an established community ice flow model and standard benchmark experiments designed as an idealized representation of Pine Island Glacier, West Antarctica. While initializing an $n=3$ model to match observations of an $n=4$ ice sheet is possible, we find that incorrectly assuming $n=3$ when in fact $n=4$ dramatically underestimates rates of sea-level rise. The scale of this error grows nonlinearly with the magnitude of the climate forcing, acting to increase projection uncertainties. Additionally, we find that models often account for this stress-dependent rheology mismatch during model initialization in a way that masks this rheological effect in the short term while leaving model outputs vulnerable to larger biases in longer-term projections. Initializations to observations of Pine Island Glacier display similar rheology-mismatch fingerprints to our idealized example.

climate sensitivity↗
DOE OSTI2025

Evolution of intermittent filaments in the scrape-off layer of NSTX

Filamentary structures naturally arise from background turbulence in the scrape-off layer (SOL) of plasmas, leading to significant particle and heat transport that can degrade overall plasma confinement. This enhanced transport can contribute to unacceptably high heat loads on plasma-facing components. As such, understanding the physics of SOL plasma filaments is critical for predicting and mitigating their effects in future fusion devices. On the National Spherical Torus Experiment (NSTX), plasma filaments—commonly referred to as “blobs”—were investigated using the gas-puff imaging (GPI) diagnostic in the edge and SOL regions. The analysis involved identifying, segmenting, and tracking the characteristic contours of the blobs in each frame of the GPI video sequences. Their evolution was characterized through shape descriptors, velocity, and angular velocity derived from their contour coordinates. The results indicate that as the blob area increases, their shapes become more concave and less circular, suggesting reduced structural stability in larger blobs. This result aligns with previous theoretical results where it was shown that larger blobs are more susceptible to instabilities [Krasheninnikov et al., J. Plasma Phys. 74, 679–717 (2008) and D'Ippolito et al., Phys. Plasmas 18, 060501 (2011)]. A positive correlation was observed between radial velocity and radial position, suggesting radially outward acceleration of the filaments, potentially driven by decreasing viscous drag toward the far SOL. Interestingly, blobs in background SOL turbulence exhibited minimal spinning in contrast to filaments originating from edge localized modes, which show substantial rotation during their paths [Lampert et al., Phys. Plasmas 29, 102502 (2022)]. Statistical analysis of the solidity and total curvature shape descriptors, along with their temporal evolution, revealed relatively broad, near-Gaussian distributions. This suggests that blob morphology is strongly influenced by stochastic turbulent processes in the surrounding plasma environment. Blob parameters were also compared with bulk plasma and radial profile measurements. Notable trends were found between blob rotation and poloidal velocity with collisionality and line-integrated density. These findings contribute to a deeper understanding of blob dynamics and provide valuable insights for refining SOL turbulence models.

Covariance and correlation↗
DOE OSTI2021

Effect of CMAS viscosity on the infiltration depth in thermal barrier coatings of different microstructures

Calcium magnesium aluminum silicate (CMAS) is one of the leading concerns for the gas turbine industry. The effects of CMAS viscosity and the coating microstructure on CMAS infiltration depth were explored by conducting a time dependent interaction study. Three CMAS compositions were used from literature, and their viscosities predicted through FactSage viscosity module were drastically different. The interaction was carried out on three different TBCs synthesized using the solution precursor plasma spray process (SPPS): two of the TBCs were made of yttrium aluminum garnet (YAG) having different microstructures that promote different modes of CMAS infiltration, and one TBC was made of gadolinium zirconate (GZO). All samples had stress relieving vertical cracks and different intensities of horizontally banded porosity, (inter pass boundaries IPBs). A concentration of 100 mg/cm 2 of CMAS was applied on the TBCs which were then subjected to a 5-minute interaction at 1300 °C. Samples were analyzed using scanning electron microscopy (SEM), electron dispersive X-ray spectroscopy (EDXS), and transmission electron microscopy (TEM). Low viscosity CMAS readily penetrated the TBCs while more viscous CMAS showed less penetration. The depth of CMAS infiltration depended on the coating microstructure. In the YAG with IPBs, the CMAS spread horizontally in the IPBs before infiltrating deeper, resulting in reduced infiltration depth compared to other samples in spite of having wider vertical cracks. TEM and EDXS analysis were performed to investigate the phases present in the CMAS-TBC interaction region in YAG. Two regions were chosen, the top TBC surface in direct contact with the sea of CMAS, and the region at the CMAS penetration ended within the coating. Finally, the results showed that no secondary phases like apatite were observed in YAG, thus it can be concluded that the arrest of CMAS happened solely because of CMAS viscosity and the short infiltration time.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Twin-screw granulation: Mechanistic understanding of the effect of material properties on key granule quality attributes through the analysis of mixing dynamics and granulation rate mechanisms

This study focuses on understanding the effect of material properties on granule quality attributes through the analysis of mixing dynamics and granulation rate mechanisms. Powder wettability, binder viscosity, and liquid-to-solid (L/S) ratio were the factors that were investigated in this study. Further, the mixing occurring inside the twin-screw granulator (TSG) was quantitatively assessed by obtaining the axial dispersion coefficient from the experimentally measured residence time distribution (RTD) curves. It was observed that the quality of the nuclei fed to the kneading zone significantly affected the mixing dynamics. The quality of nuclei was governed by nucleation kinetics, which in turn was principally affected by the liquid saturation of the nuclei and the ratio of drop penetration time and encounter time, which in turn were affected by the L/S ratio and binder viscosity respectively. The hydrophobicity of the blend mainly affected the extent of nucleation. The type of nuclei entering the kneading zone and mixing dynamics in the TSG also determined whether the granulation growth mechanism was “layering-dominant” or “viscous-dominant”. It was also shown that the resultant granule quality attributes were a reflection of the growth mechanisms. Ultimately, a mechanistic link between material properties, mixing dynamics, granulation rate mechanisms, and granule quality attributes was established.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Viscoelastic necking dynamics between attractive microgels

Hypothesis: Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer–solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. Here we hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales. Experiments: We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics. Findings: We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

42 ENGINEERING↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2023

Study on Electrostatic Separation of Quinoline Insolubles from Coal Tar Pitch

The feasibility of electrical separation in the removal of quinoline insoluble (QI) particles from coal tar pitch (CTP) was experimentally investigated. QI particle involvement prohibits the effective fabrication of high-quality value-added products, such as carbon fibers, from CTP. A substantial and sustainable CTP market exists around the world; therefore, a strong incentive exists to develop a viable technical approach to effectively and economically remove QI particles from CTP. The electrical separation method shows promise to achieve this technical goal (Cao et al., 2012). Even with the given setup (wire- cylinder adapted), critical issues remain to be addressed for this method to be applicable to the CTP: (1) identifying the wash oil used in the original QI separation, (2) understanding the mechanism of QI separation, (3) characterizing the deposit, and (4) identifying the QI. This study integrated a set of experimental and analysis techniques into the electrical separation tests to address the aforementioned issues. Key findings are as follows: 1. Electric current responses: • The electric current level of the CTP–wash oil solution reported by Cao et al. (2012) can be attained by using a mixture of 25% quinoline and 75% toluene for the CTPs examined in this study under the same electrical load condition. • The electric current tends to decrease during the test period because of the decreasing number of charged particles. • The reversed field corresponded to the configuration of electrostatic precipitation for positive corona discharge. The high electric field can result in dielectric breakdown and lead to an abrupt current surge. 2 Deposit response and solvent candidates: • Deposit of particles in CTP mixture can be effectively implemented in a wire-cylinder configuration as proposed and examined in this study. • Deposit depends on the solvents. Among the solvents tested, two- and three-part solvents that included quinoline and ethanol (i.e., 25% quinoline and 75% toluene; 50% quinoline and 50% toluene; and 33% wash oil, 33% BTX, and 34% ethanol) produced the highest deposit weight. • The deposit process examined in this study is derived from the charged particles. An appreciable relation exists between deposit weight and electric charge. 3. QI removal efficiency (RE): • The RE of the electrostatic separation can reach as high as 76.2% for Carbores, and 61.6% for Koppers. • The RE can be further enhanced if the field level increases from 0.16 kV/mm that was used in the current study to 0.23 kV/mm, according to the relation established between the RE and the electric field. 4. Testing of the modeling system: • The motion of particles is originally driven by the charge-based electric force in the cases tested. The solid particle separation mechanism is similar to that of QI deposition in CTP. • The mechanical movement of solid particles can be strongly affected by the gravitational and viscous forces in the electric field. 5. EDS analysis and chemical compositions: • The main chemical composition of deposit QI matches that of as-prepared CTP QI. The deposit QI is derived from the same group of the CTP QI. • In addition to carbon, the QI contains oxygen, sodium, aluminum, silicon, iron, and sulfur. The work for the near future is also discussed.

01 COAL, LIGNITE, AND PEAT↗