Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHEAR LAYER”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Growth of organized flow coherent motions within a single-stream shear layer: 4D-PTV measurements

Abstract This study investigates the evolution of a single-stream shear layer (SSSL) originating from a wall boundary layer past a backward-facing step. Utilizing a time-resolved 3D-Particle Tracking Velocimetry (4D-PTV) technique, we track the trajectories of fluorescent particles to gain insight into the flow characteristics of the SSSL. A compact water tunnel facility ( $$\textrm{Re}_\tau =1\,240$$ Re τ = 1 240 ) is fabricated to obtain an SSSL with a perpendicular slow entrainment stream past the separation edge. A hybrid interpolation approach that combines ensemble binning and Gaussian weighting is implemented to derive minimally filtered mean and instantaneous lower- and higher-order flow field parameters. Spanwise-dominant coherent motion accompanied by finer flow scales is observed to grow due to flow entrainment through “nibbling” actions of small-scale vortices, “engulfing” by large-scale vortices, and vortex pairing events. Furthermore, the non-zero-speed stream edge grows relatively faster than the zero-speed stream edge, showing a strong asymmetry in mixing composition across a mixing layer. The SSSL reaches self-similarity at a streamwise distance of $$\approx 55\,\theta _{0}$$ ≈ 55 θ 0 , where $$\theta _0$$ θ 0 is the initial momentum thickness from the separation edge, i.e., considerably shorter than reported in previous studies. A literature comparison of growth rate parameters raises intriguing questions regarding a potential inclusive growth scaling unifying the free shear layers. A turbulent kinetic energy (TKE) budget analysis reveals a negative production region immediately downstream of the separation edge attributed to a large positive streamwise gradient of streamwise velocity. In the self-similar region, the phase-averaged flow mapping demonstrates a larger concentration of turbulence production rate around the outer edges of spanwise vortices, specifically at the intersection of braids and vortices. Furthermore, a spatial separation exists in the regions of peak production and dissipation rates within the vortex core region favoring dissipation. The braids exhibit a larger concentration of turbulence diffusion rates, indicating their function as a conduit for exchanging turbulence between neighboring coherent motions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Observation of Nonaxisymmetric Standard Magnetorotational Instability Induced by a Free-Shear Layer

The standard magnetorotational instability (SMRI) with a magnetic field component parallel to the rotation axis is widely believed to be responsible for the fast accretion in astronomical disks. In conventional base flows with a Keplerian profile or an ideal Couette profile, most studies focus on axisymmetric SMRI, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number (Rm) more than an order of magnitude larger. Here, we report that, in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI with an azimuthal mode number m = 1 can be triggered by a free-shear layer in the base flow at Rm ≳1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally [Wang et al., Nat. Commun. 13, 4679 (2022)] is the nonaxisymmetric m =1 SMRI. Further, our finding has astronomical implications as free-shear layers are ubiquitous in celestial systems, such as the disk-star boundary layer, the solar tachocline, and the edge of planet-opened gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Hypersonic wake velocity measurements using acetone molecular tagging velocimetry

Molecular tagging velocimetry, a minimally intrusive laser diagnostic, was used to find quantitative off-body velocity measurements in the wake of both strut-mounted and free-flight spheres in hypersonic flow. Acetone gas seeded in the flow was excited using the fourth harmonic of a pulse-burst Nd:YAG laser. The experimental tests were validated and compared to simulations using a continuous Galerkin flow solver and an adaptive mesh refinement process. Near the centerline in the wake of the sphere, velocities between −105 and 65 m/s were observed due to reverse flow in the viscous shear layer. Similarly, the simulated data showed negative velocity values in the shear layer region. Outside of the shear layer, average velocities of 730 m/s were observed in both the experimental and simulated results. The overall average uncertainty for the strut-mounted and free-flight cases was estimated to have an accuracy of ±4% (±35.5m/s or 4.9%). In conclusion, this non-intrusive technique provides a velocity map of the wake behind a sphere under hypersonic conditions.

42 ENGINEERING↗
DOE OSTI2026

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗
DOE OSTI2025

Oceanic mixing and waves in the presence of a suspended canopy

Large-eddy simulations are analysed to determine the influence of suspended canopies, such as those formed in macroalgal farms, on ocean mixed layer (OML) deepening and internal wave generation. In the absence of a canopy, we show that Langmuir turbulence, when compared with wind-driven shear turbulence, results in a deeper OML and more pronounced internal waves beneath the OML. Subsequently, we examine simulations with suspended canopies of varying densities located in the OML, in the presence of a background geostrophic current. Intensified turbulence occurs in the shear layer at the canopy’s bottom edge, arising from the interaction between the geostrophic current and canopy drag. Structures resembling Kelvin–Helmholtz (KH) instability emerge as the canopy shear layer interacts with the underlying stratification, radiating internal waves beneath the OML. Both intensified turbulence and lower-frequency motions associated with KH-type structures are critical factors in enhancing mixing. Consequently, the OML depth increases by up to a factor of two compared with cases without a canopy. Denser canopies and stronger geostrophic currents lead to more pronounced KH-type structures and internal waves, stronger turbulence and greater OML deepening. Additionally, vertical nutrient transport is enhanced as the OML deepens due to the presence of the canopy. Considering that the canopy density investigated in this study closely represents offshore macroalgal farms, these findings suggest a mechanism for passive nutrient entrainment conducive to sustainable farming. Overall, this study reveals the intricate interactions between the suspended canopy, turbulent mixing and stratification, underscoring their importance in reshaping OML characteristics.

Bo, Tong (ORCID:0000000260300561)↗
DOE OSTI2025

Environmental Influences on Deep Convective Upscale Growth Rate in Central Argentina From a Convection‐Permitting Simulation

This study uses a convection‐permitting model simulation to describe the environmental conditions under which convective upscale growth occurs in central Argentina, particularly examining environmental parameters when deep convection initially forms that could differentiate the rate of initial upscale growth. Simulated mesoscale convective systems (MCSs) are separated into slow and rapid growth by the rate of spatial growth from convection initiation until reaching the MCS scale. A low‐level jet (LLJ) is found more frequently near the deep convection that experiences rapid growth to an MCS, but its presence alone is not predictive of rapid growth. Using spatially‐averaged parameters, we find that rapid growth to MCSs also occurs in environments that are significantly more thermodynamically favorable with greater low‐level moisture and instability. Fewer significant differences are found in the kinematic environment with only the 0–2 km vertical wind shear magnitude being significantly larger for rapid growth MCSs compared to slow growth MCSs, potentially related to LLJs often peaking near this height. When focusing only on MCSs with the slowest and fastest growth rates, elevated‐layer shear is significantly smaller for very rapid growth MCSs, suggesting elevated‐layer shear may help discriminate between the upper and lower bounds of growth rate. Finally, when upscale growth occurs near the Sierras de Córdoba (SDC) with a LLJ present, rapid growth is also supported by favorable wind shear orientation. However, this does not hold for upscale growth occurring away from the SDC, highlighting the importance of interpreting shear direction relative to the orientation of features initiating deep convection.

Sasaki, Clayton R. S. [Univ. of Washington, Seattl↗
DOE OSTI2024

Understanding the L-H isotope effect at the DIII-D tokamak and advancements in synthetic turbulence diagnostics

Abstract It is determined that while heat flux differences between hydrogen and deuterium isotope experiments result from natural differences in carbon impurity content at DIII-D, it is not the origin of the low to high confinement mode (L-H) transition isotope effect. More specifically, a two times larger edge radial electric field in hydrogen compared to deuterium is uncovered and believed to play an important role. The origin of this radial electric field difference is determined to have two possible origins: differences in poloidal rotation and turbulent Reynolds stress in the closed field line region, and increased outer strike point temperatures and space potentials on open field lines. Experimental observations from both profile and turbulence diagnostics are supported by nonlinear gyrokinetic simulations using the code CGYRO. Simulations illustrated heat transport isotope effects in the plasma edge and shear layer resulting from differences in impurity content, electron non-adiabaticity, and main ion mass dependent E × B shear stabilization. Turbulence prediction comparisons from flux-matched CGYRO simulations to experimental measurements including electron temperature, density and velocity fluctuations are found to be in good agreement with available data. A dedicated DIII-D experiment in hydrogen was performed to seed more carbon than naturally occurring, to match deuterium experiments, and possibly reduce the L-H power threshold based on gyro-kinetic predictions. To our surprise, while ion temperature gradient (ITG) turbulence was stabilized, nodiscernible change in L-H power threshold were observed in these special hydrogen experiments. In particular, it is noticed that the edge radial electric field and Reynolds stress were observed as nearly unchanging in the presence of ITG stabilization. These experimental data have enabled a more comprehensive picture of the multitude of isotope effects at play in fusion experiments, and the important potential connection between the confined and unconfined plasma regions in regulating L-H transition dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comparative studies of cross-phase dynamics in turbulent momentum flux and particle flux at the tokamak edge

Turbulent transport events, including turbulent transport flux of momentum (i.e., turbulent momentum flux or Reynolds stress) and turbulent transport flux of particle (i.e., turbulent particle flux), have important effects on the confinement performance of magnetic confinement fusion devices. Poloidal Reynolds stress is the ensemble average of the product of radial velocity fluctuations and poloidal velocity fluctuations, i.e., $\langle {\widetilde{v}}_{r}{\widetilde{v}}_{\theta }\rangle$. Turbulent particle flux is the ensemble average of the product of radial velocity fluctuations and density fluctuations, i.e., $\langle \widetilde{n}{\widetilde{v}}_{r}\rangle$. Changes in either amplitude of fluctuations or cross phase between fluctuations can cause changes in turbulent transport. In this paper, cross-phase dynamics in the Reynolds stress and turbulent particle flux at the tokamak edge are studied in detail. Reynolds stress and turbulent particle flux are, respectively, written as the product of fluctuation amplitudes and an average cross-phase factor. The mathematical expressions of the average cross-phase factors are derived. The average cross-phase factors and the power spectra of cross phase are obtained using experimental measurement data. It is found that the cross-phase dynamics in Reynolds stress and particle flux are very different. Reynolds stress is found to be more sensitive to cross phase than particle flux is. In the strong $E\times B$ shear layer, spatial slips of cross phase lead to the obvious radial gradient of Reynolds stress. In the no/weak $E\times B$ shear region, the cross phase in Reynolds stress tends to lock. Here, phase locking refers to that the power spectra of phase tend to distribute around a fixed phase which does not change with radial position, while phase slip means that the power spectra of cross phase tend to distribute around a phase that varies with radial position. Phase slip or locking mainly describes the central phase weighted by the power spectra, while the phase scattering mainly describes the dispersion of the power spectrum distribution of the phase. The increased scattering of cross phase, which indicates the power spectra distribution of the phase is more dispersed, contributes to the decreased Reynolds stress for higher collisionality. The cross phase in particle flux tends to lock in both strong and no/weak shear regions. The degree of scattering of cross phase in the particle flux does not change obviously as collisionality increases. For higher collisionality, it is the increased density fluctuation amplitude rather than cross-phase dynamics that leads to the increased particle flux. The underlying physical mechanism that causes Reynolds stress and particle flux to exhibit different phase dynamics is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗
DOE OSTI2026

Impact of atmospheric turbulence on performance and loads of wind turbines: knowledge gaps and research challenges

Wind energy harvesting from the atmosphere takes place in the atmospheric boundary layer. The boundary layer shear and buoyancy create three-dimensional turbulent eddies spanning a range of scales that form a continuous forward cascade of kinetic energy to the smallest scales of motion where energy is dissipated. Large-scale atmospheric circulations modulate the boundary layer turbulence, characterized by coherence and intermittency. As wind turbines grow in size and the integrated control of both turbines and wind farms spans greater distances, the relationship between the scales of atmospheric turbulence and the design and operation of wind energy facilities has entered new territory. The boundary layer turbulence impacts both wind turbine power production and turbine loads. Optimizing wind turbine and wind farm performance requires an understanding of how turbulence affects both wind turbine efficiency and reliability. While the characteristics of atmospheric boundary layer turbulence have been observed and studied in detail over the last few decades, there are still significant gaps in our understanding of the impact of turbulence on wind power resources and wind farm operations. This paper outlines the current state of turbulence research relevant to wind energy applications and points to gaps in our knowledge that need to be addressed to effectively utilize wind resources.

Kosović, Branko [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2025

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗
DOE OSTI2025

Stable low diffusion flux splitting schemes on unstructured meshes

Shock instabilities are shown to manifest in modern low-diffusion flux-vector splitting (FVS) schemes when used on unstructured meshes, or situations where shocks do not align with the mesh lines. These instabilities occur irrespective of the Mach number of the shock. Three types of dissipative mechanisms that suppress these instabilities are presented. These mechanisms are carefully designed in order to affect only problematic regions of the flux-splittings. The AUSM + and LDFSS schemes are stabilized using the proposed modifications. It is shown that the added dissipation improves the shock behavior of AUSM and LDFSS on unstructured meshes. It is also shown that the AUSM + -up scheme is prone to the “carbuncle” instability, a specific type of shock instability, when used on unstructured meshes. The modifications proposed in this work do not lead to carbuncle instabilities for the problems considered here. Furthermore, the modified schemes are shown to satisfy certain properties that are crucial for accurate shear layer computations, such as stationary contact preservation. Using benchmark problems, it is demonstrated that despite the diffusion added for stabilization, these schemes are not overly diffusive. Furthermore, due to these advantages, the modified FVS schemes presented here are promising candidates for high-speed compressible flow computations on unstructured meshes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗
DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2026

Suppression of ion temperature gradient modes by Alfvén activity above a drive threshold in DIII-D

A recent study demonstrates that the suppression of ion temperature gradient (ITG) modes can occur during the nonlinear evolution of toroidicity-induced Alfvén eigenmodes (TAEs) in the DIII-D tokamak [Du et al., Phys. Rev. Lett. 135, 265101 (2025)]. In that work, ITG was suppressed by the formation of a narrow, TAE-induced shear flow layer, whose shearing rate exceeded the ITG decorrelation rate. Here, the shear flow arises from an imbalance between Reynolds and Maxwell stress forces, as the TAE departs from the conventional shear Alfvén wave polarization. This follow-up paper systematically identifies the plasma conditions required for robust ITG suppression through a series of comparative experiments. The results show that TAEs routinely suppress ITG turbulence in plasmas with higher local safety factor (q), elevated fast-ion beta, and larger populations of fast ions on passing orbits, i.e., conditions consistent with the substantial TAE drive. Database analysis further reveals the existence of a threshold in TAE drive for ITG suppression. That is, once the drive exceeds a critical value within a favorable q window, the system undergoes the nonlinear bifurcation process, characterized by a sharp increase in TAE saturated amplitudes at a nearly fixed fast ion drive, suppression of ITG turbulence, and the formation of an internal thermal transport barrier at the localized radii.

Du, X. D. [General Atomics, San Diego, CA (United ↗
DOE OSTI2024

Overview of the recent experimental research on the J-TEXT tokamak

The J-TEXT capability is enhanced compared to two years ago with several upgrades of its diagnostics and the increase of electron cyclotron resonance heating (ECRH) power to 1 MW. With the application of electron cyclotron wave (ECW), the ECW assisted plasma startup is achieved; the tearing mode is suppressed; the toroidal injection of 300 kW ECW drives around 24 kA current; fast electrons are generated with toroidal injected ECW and the runaway current conversion efficiency increases with ECRH power. The mode coupling between 2/1 and 3/1 modes are extensively studied. The coupled 2/1 and 3/1 modes usually lead to major disruption. Their coupling can be either suppressed or avoided by external resonant magnetic perturbation fields and hence avoids the major disruption. It is also found that the 2/1 threshold of external field is significantly reduced by a pre-excited 3/1 mode, which can be either a locked island or an external kink mode. The disruption control is studied by developing prediction methods capable of cross tokamak application and by new mitigation methods, such as the biased electrode or electromagnetic pellet injector. The high-density operation and related disruptions are studied from various aspects. Approaching the density limit, the collapse of the edge shear layer is observed and such collapse can be prevented by applying edge biasing, leading to an increased density limit. The density limit is also observed to increase, if the plasma is operated in the poloidal divertor configuration or the plasma purity is increased by increasing the pre-filled gas pressure or ECRH power during the start-up phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Visualizing the mass transfer flow in direct-impact accretion

We use a variety of visualization techniques to display the interior and surface flows in a double white dwarf binary undergoing direct-impact mass transfer and evolving dynamically to a merger. The structure of the flow can be interpreted in terms of standard dynamical, cyclostrophic, and geostrophic arguments. We describe and showcase some visualization and analysis techniques of potential interest for astrophysical hydrodynamics. In the context of R Coronae Borealis stars, we find that mixing of accretor material with donor material at the shear layer between the fast accretion belt and the slower rotating accretor body will always result in some dredge-up. We also discuss briefly some potential applications to other types of binaries.

79 ASTRONOMY AND ASTROPHYSICS↗