SEARCH · Engineering Papers

Results for “Laminarization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Enhancement of operational Flexibility of Power Plants Using IN740

Headers are crucial components within diverse industries, especially in the energy sector, as they enable the efficient transfer of fluids. The selection of materials for headers is determined by their specific applications; for example, Grade 91 and Grade 92 steels are commonly utilized. In our research, we recommend the adoption of INCONEL alloy 740 due to its exceptional robustness and heat-resistant characteristics. One critical parameter in header systems is the heat transfer coefficient, which directly influences the efficiency of heat exchange processes and can thereby impact the structural integrity of the headers. This coefficient is closely linked with factors such as the Nusselt number, which is influenced by fluid flow characteristics and thermal properties. The flow direction within headers, whether unidirectional or multidirectional—significantly affects the overall dynamics of heat transfer. This research focuses on exploring the behavior of headers, a specific pipeline system component, using ANSYS simulation software. The study aims to predict heat transfer and mechanical behavior within headers under various conditions through steady-state and transient simulations (parts 1 and 2 of the report). Key parameters such as heat transfer coefficient, velocity, and temperature are examined with the goal of optimizing header design. Part 3 of the report addresses the critical yet underexplored relationship between pressure drop and heat transfer coefficients in the transient flow regime within headers, specifically focusing on smooth horizontal circular tubes. Limited experimental work has been conducted in this area, prompting the need for comprehensive analysis. By leveraging machine learning techniques, this research aims to establish a correlation between pressure drop and heat transfer across various flow conditions, including laminar, transient, quasi-turbulent, and turbulent regimes. The data utilized for this analysis were meticulously gathered from existing literature, capturing simultaneous measurements of pressure drops and heat transfer. Part 4: In the context of power plants, header pipes are essential components that significantly influence system performance by facilitating the collection and distribution of steam. This report highlights the critical role of header pipes in enhancing reliability, efficiency, and overall power plant performance. A key aspect of this investigation is shape optimization, which aims to maximize performance while minimizing material usage. By focusing on shape optimization, this research contributes to improved efficiency and a reduced environmental footprint for power plant installations. The methodology developed in this study emphasizes optimizing header shapes to decrease reliance on expensive alloy materials and lower maintenance costs. Furthermore, a case study was conducted using a header from an operational power plant to validate the proposed optimization techniques. Finally, part 5 addresses the challenges associated with flexible operations in boiler systems, the study outlines several strategies for enhancing the durability and reliability of steam headers. Key approaches include material selection, which involves utilizing advanced materials with superior high-temperature properties and enhanced fatigue resistance to extend the lifespan of steam headers. Design modifications are also recommended, focusing on implementing changes that mitigate thermal stress and cyclic loading to reduce the likelihood of failures under varying operating conditions. Additionally, the study emphasizes the importance of regular inspection and monitoring by establishing rigorous protocols to detect early signs of damage, allowing for timely maintenance and minimizing the risks of catastrophic failures. Furthermore, operational guidelines are developed to minimize the frequency and severity of thermal transients, ensuring stable and efficient operations. Collectively, these strategies provide a robust framework for overcoming the unique challenges posed by flexible operation in boiler systems, ultimately contributing to improved reliability and performance in power plant operations. Through this comprehensive approach, the study not only enhances the understanding of pressure drop and heat transfer relationships but also promotes advancements in design and operation that will benefit the entire power generation industry.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2025

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗

DOE OSTI2024

Overview of On-Line Optical Measurements at High Pressure for Flue Gases, Particulates and Acid-dew Point of Pressurized Oxy-Combustion

Optical flow cells are critical measurement interfaces, yet sampling under harsh conditions — high pressure, high temperature, particles, moisture, or corrosive gases — makes it difficult to maintain optical quality without perturbing the measurement. To address this challenge, a new flow cell was developed using a laminar coaxial flow field that separates the purge and sample flows. A dedicated test system was built to evaluate particle size distribution (PSD) measurements using a Malvern Panalytical Insitec analyzer. Results demonstrated that the sample flow alone defines the measurement zone, while the purge flow effectively shields the optical windows from deposition, eliminating sampling bias. The flow cell enables reliable PSD measurement under high pressure and temperature in moist, corrosive environments. As a key demonstration, the instrument was successfully deployed for on-line PSD measurement of flue gas from a 100 kWth pressurized oxy-coal combustor at 15 bara.

Cheng, Mao↗

DOE OSTI2023

Hypersonic Fluid–Structure Interaction on a Cone–Slice–Ramp Geometry

In this study, fluid–structure interactions were measured between a representative control surface and the hypersonic flow deflected by it. The control surface is simplified as a spanwise finite ramp placed on a longitudinal slice of a cone. The front surface of the ramp contains a thin panel designed to respond to the unsteady fluid loading arising from the shock-wave/boundary-layer interactions. Experiments were conducted at Mach 5 and Mach 8 with ramps of different angles. High-speed schlieren captured the unsteady flow dynamics and accelerometers behind the thin panel measured its structural response. Panel vibrations were dominated by natural modes that were excited by the broadband aerodynamic fluctuations arising in the flowfield. However, increased structural response was observed in two distinct flow regimes: 1) attached or small separation interactions, where the transitional regime induced the strongest panel fluctuations. This was in agreement with the observation of increased convective undulations or bulges in the separation shock generated by the passage of turbulent spots, and 2) large separated interactions, where shear layer flapping in the laminar regime produced strong panel response at the flapping frequency. In addition, panel heating during the experiment caused a downward shift in its natural mode frequencies.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Assessment of the Amplification Factor Transport Transition Model for High-Mach Number Flows

In its present state, the Amplification Factor Transport (AFT) model provides an efficient alternative for high-fidelity stability methods for complex geometries subject to low-speed conditions considered in the industrial-production CFD environment. Its performance has not yet been documented for the typical transition scenarios in hypersonic conditions. The goal of the present study is to apply the AFT model to a yawed- and flared-cone configuration, for which the laminar flows respectively support the crossflow and second-mode instability mechanisms, and to compare the amplification results against high-fidelity results computed with Linear Parabolized Stability Equations (LPSE). High-resolution viscous-flow solutions are provided with OVERFLOW and COFFE. For the yawed-cone case, the AFT model yields the largest amplification in a region of the flow where the expected crossflow instability is not amplified (in the leeward symmetry plane) according to LPSE. For the flared-cone case with a wall-temperature slightly lower than the adiabatic value, the AFT model yields extremely large amplification factors (in excess of 100), exceeding the computed LPSE amplification of second-mode instability results by a factor 14. Upon considering a cold-wall condition, the model shows much smaller amplification factors, which is the opposite with respect to the expected behavior of the second mode. The AFT model does not adequately capture the crossflow and second-mode instability mechanisms in its present state.

Groot, Koen J.↗

DOE OSTI2024

1D Simulation of Avalanche to Streamer to Spark Transition of Plasma Discharge in Ammonia-Air Combustion

The backward problem of plasma assisted combustion emphasizes evaluating the effect of the evolving thermochemical state on the plasma discharge. This paper investigates the dependence of avalanche to streamer to spark formation dynamics and kinetics on the gas composition and temperature at different points in an ammonia-air premixed laminar flame using a self-consistent multigrid-based 1D plasma solver. Different values of ..alpha.., the coefficient for effective ionization events per unit length, have been reported for electron avalanches in air and stoichiometric NH3-air mixtures. The streamer inception has been shown to obey the Meek's criterion. An exponential reduction in streamer and spark formation time has been observed from plasma simulations at different points in the unburnt, pre-heat zone, reaction zone and the fully burnt regions of the premixed flame. While the enhancement of the reduced electric field with increasing temperature affects effective ionization, there exists a minimum breakdown field for streamer formation, which does not vary proportionally with the changing number density of the gas. The change in the mixture from reactants (NH3, O2, N2) to products of complete combustion of ammonia in air (N2, H2O) has also been shown to affect the streamer and spark formation. Finally, the major pathways during the streamer and spark phases which are responsible for producing important radicals used in combustion of NH3 are also discussed.

ADVANCED PROPULSION SYSTEMS,INORGANIC, ORGANIC, PH↗

DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗

DOE OSTI2025

Mesoscale fractal whey protein particles derived from microscale linear-shaped protein assemblies (Part 1): Manufacturing method and particle characteristics

Whey protein isolates (WPI) are widely used in processed foods for their versatile functional properties. Modifying the structural properties of proteins by assembling them into mesoscale or microscale particles may improve their functionality and broaden their applications. This study aims to manufacture and characterize mesoscale whey protein particles (WPP) derived from WPI. Two types of WPP, WPP1 (0.05 mL/min) and WPP2 (0.25 mL/min), were prepared through a multistep approach involving liquid antisolvent (LAS) precipitation, heat treatment, and microfluidization. Liquid antisolvent precipitation was performed by injecting a 20% (wt/vol) WPI dispersion (pH 7) into an ethanol-glycerol mixture (75:25, vol/vol) under laminar flow, followed by heat treatment at 80°C for 20 min as a particle hardening step. This process produced stable fiber- and ribbon-shaped whey protein assemblies (WPA), which served as precursors to WPP. Subsequent microfluidization (150 MPa, 6 passages) reduced the size of WPA, yielding mesoscale WPP with irregular morphologies and a more uniform size distribution, as revealed by microscopy and dynamic light scattering. ζ-Potential and fluorescence labeling indicated higher surface charge and surface hydrophobicity of WPP compared with untreated WPI. The WPP showed internal mass fractal and surface fractal structures at larger length scales, analyzed using small-angle X-ray scattering. Fourier transform infrared spectroscopy demonstrated an increased fraction of intermolecular β-sheets in WPP, suggesting that hydrogen bonding contributed to their formation. Gel electrophoresis confirmed that disulfide bonds served as the primary cross-links stabilizing the WPP structure. Furthermore, turbidity measurements showed that WPP exhibited superior colloidal phase stability compared with untreated WPI and maintained high colloidal stability under both acidic and neutral pH conditions.

Antisolvent precipitation↗

DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗

DOE OSTI2022

Numerical and Experimental Study of a Novel Additively Manufactured Metal-Polymer Composite Heat-Exchanger for Liquid Cooling Electronics

In order to meet increasing power-dissipation requirements of the electronics industry, compact, low-cost, and lightweight heat exchangers (HXs) are desired. With proper design, materials, and manufacture, polymer composite heat exchangers could meet these requirements. This paper presents a novel crossflow air-to-water, low-cost, and lightweight metal-polymer composite HX. This HX, which is entirely additively manufactured, utilizes a novel cross-media approach that provides direct heat exchange between air and liquid sides by using connecting fins. A robust numerical model was developed, which includes the dimensional effects of additive manufacturing. The study consists of a simplified 3D CFD model based on ellipsoidal-shaped staggered tube banks for the laminar range. It then uses an analytical approach to compute entire HX performance. The model is validated experimentally within 8% for thermal performance, 12% for air-side impedance, and 18% for water-side impedance. Finally, HX is compared with a conventional CPU radiator and performs within 10% of the conventional unit for reasonable flow rates and pressure-drop ranges. Moreover, HX also provides added design and cost advantages over the conventional unit, which makes the HX a potential candidate for electronic cooling applications.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2025

A Novel Optical Instrument for On-Line Measurement of Particle Size Distribution—Application to Clean Coal Technologies

A flow cell is a critical measurement interface for many optical instruments. However, the flows are often sampled under harsh conditions, such as under high pressure and/or high temperature, in the presence of particles, moisture, vapors with high dew points or corrosive gases. Therefore, obtaining a high-optical-quality flow cell that does not perturb the measurement is a significant challenge. To address this challenge, we proposed a new flow cell that employs a unique laminar coaxial flow field (for the purge and sample flows). A test system was built to conduct particle size distribution (PSD) measurements with no sampling bias using a state-of-the-art analyzer (Malvern Panalytical Insitec). The results revealed that the measurement zone is well defined solely by the sample flow, and the optical windows are well protected by the purge flow, with minimal risk of any depositions from the sample flow. Using this flow cell, the Insitec can successfully measure PSD under high pressure and temperature under moist, corrosive conditions without generating any sampling bias. Importantly, we successfully applied this flow cell for on-line PSD measurement for the flue gas of a 100 kWth pressurized oxy-coal combustor operating at 15 bara.

Cheng, Mao (ORCID:0000000213273900)↗

DOE OSTI2020

Cryogenic Liquid Jets for High Repetition Rate Discovery Science

In this paper, we present a detailed procedure for the operation of continuous, micron-sized cryogenic cylindrical and planar liquid jets. When operated as described here, the jet exhibits high laminarity and stability for centimeters. Successful operation of a cryogenic liquid jet in the Rayleigh regime requires a basic understanding of fluid dynamics and thermodynamics at cryogenic temperatures. Theoretical calculations and typical empirical values are provided as a guide to design a comparable system. This report identifies the importance of both cleanliness during cryogenic source assembly and stability of the cryogenic source temperature once liquefied. The system can be used for high repetition rate laser-driven proton acceleration, with an envisioned application in proton therapy. Other applications include laboratory astrophysics, materials science, and next-generation particle accelerators.

42 ENGINEERING↗

DOE OSTI2022

Formation and Reconnection of Electron Scale Current Layers in the Turbulent Outflows of a Primary Reconnection Site

We simulate with 3D particle in cell, the spontaneous formation of turbulent outflows in an initially laminar 3D reconnecting current layer. We observe the formation of many secondary current layers and reconnection sites in the outflow. The approach we follow is to study each individual feature within the turbulent outflow. To identify all clusters of current in the outflow we use a clustering technique widely used in unsupervised machine learning: density-based spatial clustering of applications with noise. Once the clusters are identified we measure their size and compute reconnection indicators to establish which are undergoing reconnection. With this analysis we establish that the size of the current clusters reaches all the way from its initial system scale down to subelectron skin depth scale. We observe that the smaller current clusters are more prone to reconnecting and to releasing energy. We then find the process of reconnection of the smaller current cluster to be of the recently observed electron-only type that leaves the ions essentially unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Winds and Disk Turbulence Exert Equal Torques on Thick Magnetically Arrested Disks

Abstract The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds,that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration (1.2 × 10 5 r g /c) simulation of a rapidly rotating (a= 0.9) black hole feeding from a thick (H/r∼ 0.3), adiabatic, magnetically arrested disk (MAD), whose dynamically important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (i.e., winds and jets, respectively). By carefully disentangling the various angular momentum transport processes within the system, we demonstrate the novel result that disk winds and disk turbulence both extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates of L ̇ ∝ r 0.9 and M ̇ ∝ r 0.4 , respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma-ray bursts.

Astronomy & Astrophysics↗

DOE OSTI2025

Transition to Petschek Reconnection in Subrelativistic Pair Plasmas: Implications for Particle Acceleration

While relativistic magnetic reconnection in pair plasmas has emerged in recent years as a candidate for the origin of radiation from extreme astrophysical environments, the corresponding subrelativistic pair-plasma regime has remained less explored, leaving open the question of how relativistic physics affects reconnection. In this paper, we investigate the differences between these regimes by contrasting two-dimensional particle-in-cell simulations of reconnection in pair plasmas with relativistic magnetization (σ ≫ 1) and subrelativistic magnetization (σ < 1). By utilizing unprecedentedly large domain sizes and outflow boundary conditions, we demonstrate that lowering the magnetization results in a change in the reconnection geometry from a plasmoid chain to a Petschek geometry, where laminar exhausts bounded by slow-mode shocks emanate from a single diffusion region. We attribute this change to the reduced plasmoid production rate in the low-σ case: When the secondary tearing rate is sufficiently low, plasmoids are too few in number to prevent the system from relaxing into a stable Petschek configuration. This geometric change also affects particle energization: We show that while high-σ plasmoid chains generate power-law energy spectra, low-σ Petschek exhausts merely heat incoming plasma and yield negligible nonthermal acceleration. These results have implications for predicting the global current sheet geometry and the resulting energy spectra in a variety of systems.

Plasma astrophysics↗

DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗