Source author record

Rodolfo Ostilla-Mónico

Rodolfo Ostilla-Mónico appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2025arXiv

Vortex-magnetic competition and regime transitions in antiparallel flux tubes

Vortex-magnetic interactions shape magnetohydrodynamic (MHD) turbulence, influencing energy transfer in astrophysical, geophysical, and industrial systems. On the Sun, granular-scale vortex flows couple strongly with magnetic fields, channeling energy into the corona. At high Reynolds numbers, vorticity and magnetic fields are nearly frozen into the charged fluid, and MHD flows emerge from the Lorentz force mediated interactions between coherent vortex structures in matter and the field. To probe this competition in a controlled setting, we revisit the canonical problem of two antiparallel flux tubes. By varying the magnetic flux threading each tube--and thus sweeping the interaction parameter $N_i$, which gauges Lorentz-to-inertial force balance--we uncover three distinct regimes: vortex-dominated joint reconnection, instability-triggered cascade, and Lorentz-induced vortex disruption. At low $N_i$, classical vortex dynamics dominate, driving joint vortex-magnetic reconnection and amplifying magnetic energy via a dynamo effect. At moderate $N_i$, the system oscillates between vorticity-driven attraction and magnetic damping, triggering instabilities and nonlinear interactions that spawn secondary filaments and drive an energy cascade. At high $N_i$, Lorentz forces suppress vortex interactions, aligning the tubes axially while disrupting vortex cores and rapidly converting magnetic to kinetic energy. These findings reveal how the inertial-Lorentz balance governs energy transfer and coherent structure formation in MHD turbulence, offering insight into vortex-magnetic coevolution in astrophysical plasmas.

preprint2022arXiv

Transition between Boundary-Limited Scaling and Mixing-Length Scaling of Turbulent Transport in Internally Heated Convection

Heat transport in turbulent thermal convection increases with the thermal forcing, but in almost all studies the rate of this increase is slower than it would be if transport became independent of the molecular diffusivities -- the heat transport scaling is slower than the mixing-length (or `ultimate') scaling. In configurations driven by either thermal boundary conditions or internal heating, thermal boundary layers instead lead to a boundary-limited (or `classical') scaling. With net-zero internal heating and cooling in different regions, mixing-length scaling can occur because heat need not cross a boundary. We report numerical simulations in which heating and cooling are unequal, as in various natural systems, at a Prandtl number of unity. As heating and cooling rates are made closer, the scaling exponent of heat transport varies from its boundary-limited value to its mixing-length value.

preprint2020arXiv

Double maxima of angular momentum transport in $η=0.91$ TC turbulence

We use experiments and direct numerical simulations to probe the phase-space of low-curvature Taylor--Couette (TC) flow in the vicinity of the ultimate regime. The cylinder radius ratio is fixed at $η=r_i/r_o=0.91$. Non-dimensional shear drivings (Taylor numbers $\text{Ta}$) in the range $10^7\leq\text{Ta}\leq10^{11}$ are explored for both co- and counter-rotating configurations. In the $\text{Ta}$ range $10^8\leq\text{Ta}\leq10^{10}$, we observe two local maxima of the angular momentum transport as a function of the cylinder rotation ratio, which can be described as either as "co-" and "counter-rotating" due to their location or as "broad" and "narrow" due to their shape. We confirm that the broad peak is accompanied by the strengthening of the large-scale structures, and that the narrow peak appears once the driving (Ta) is strong enough. As first evidenced in numerical simulations by Brauckmann \emph{et al.}~(2016), the broad peak is produced by centrifugal instabilities and that the narrow peak is a consequence of shear instabilities. We describe how the peaks change with $\text{Ta}$ as the flow becomes more turbulent. Close to the transition to the ultimate regime when the boundary layers (BLs) become turbulent, the usual structure of counter-rotating Taylor vortex pairs breaks down and stable unpaired rolls appear locally. We attribute this state to changes in the underlying roll characteristics during the transition to the ultimate regime. Further changes in the flow structure around $\text{Ta}\approx10^{10}$ cause the broad peak to disappear completely and the narrow peak to move. This second transition is caused when the regions inside the BLs which are locally smooth regions disappear and the whole boundary layer becomes active.

preprint2019arXiv

Large-scale structures in high Reynolds number rotating Waleffe flow

We perform direct numerical simulations of rotating turbulent Waleffe flow, the flow between two parallel plates with a sinusoidal streamwise shear driving force, to study the formation of large-scale structures and the mechanisms for momentum transport. We simulate different cyclonic and anti-cyclonic rotations in the range of dimensionless rotation numbers (inverse Rossby numbers) $R_Ω$ $\in$ $[-0.16, 2.21]$, and fix the Reynolds number to $Re=3.16\times 10^3$, large enough such that the shear transport is almost entirely due to Reynolds stresses and viscous transport is negligible. We find an optimum rotation in anti-cyclonic regime at $R_Ω=0.63$, where a given streamwise momentum transport in the wall-normal direction is achieved with minimum mean energy of the streamwise flow. We link this optimal transport to the strength of large scale structures, as was done in plane Couette by Brauckmann \& Eckhardt (J. Fluid Mech., 815, 2017). Furthermore, we explore the large-scale structures and their behaviour under spanwise rotation, and find disorganized large structures at $R_Ω=0$ but highly organized structures in the anti-cyclonic regime, similar to the rolls in rotating plane Couette and turbulent Taylor Couette flow. We compare the large scale structures of plane Couette flow and Waleffe flow, and observe that the streamwise vorticity is localized inside the cores of the rolls. We show that the rolls take energy from the mean flow at long time-scales, and relate these structures to eigenvalues of the streamfunction.

preprint2015arXiv

A multiple--resolution strategy for Direct Numerical Simulation of scalar turbulence

In this paper a numerical procedure to simulate low diffusivity scalar turbulence is presented. The method consists of using a grid for the advected scalar with a higher spatial resolutions than that of the momentum. The latter usually requires a less refined mesh and integrating both fields on a single grid tailored to the most demanding variable, produces an unnecessary computational overhead. A multiple resolution approach is used also in the time integration in order to maintain the stability of the scalars on the finer grid. The method is the more advantageous the less diffusive the scalar is with respect to momentum, therefore it is particularly well suited for large Prandtl or Schmidt number flows. However, even in the case of equal diffusivities the present procedure gives CPU time and memory occupation savings. The reason is that the absence of the pressure term in the scalar equation leads to much steeper gradients in the scalar field as compared to the velocity field.

preprint2015arXiv

Inertial waves and mean velocity profiles in a rotating pipe and a circular annulus with axial flow

In this paper we solve the inviscid inertial wave solutions in a circular pipe or annulus rotating constantly about its axis with moderate angular speed. The solutions are constructed by the so-called helical wave functions. We reveal that the mean velocity profiles must satisfy certain conditions to accommodate the inertial waves at the bulk region away from boundary. These conditions require the axial and azimuthal components of the mean velocity take the shapes of the zeroth and first order Bessel functions of the first kind, respectively. The theory is then verified by data obtained from direct numerical simulations for both rotating pipe and circular annulus, and excellent agreement is found between theory and numerical results. Large scale vortex clusters are found in the bulk region where the mean velocity profiles match the theoretical predictions. The success of the theory in rotating pipe, circular annulus, and streamwise rotating channel suggests that such inertial waves are quite common in wall bounded flow with background rotation.

preprint2015arXiv

Logarithmic mean temperature profiles and their connection to plume emissions in turbulent Rayleigh-Bénard convection

Two-dimensional simulations of Rayleigh-Bénard convection at $Ra = 5\times10^{10}$ show that vertical logarithmic mean temperature profiles can be observed in regions of the boundary layer where thermal plumes are emitted. The profile is logarithmic only in these regions and not in the rest of the boundary layer where it is sheared by the large scale wind and impacted by plumes. In addition, the logarithmic behavior is not visible in the horizontal average. The findings reveal that the temperature profiles are strongly connected to thermal plume emission and support a perception that parts of the boundary layer can be turbulent, while others are not. The transition to the ultimate regime, in which the boundary layers are considered to be fully turbulent, can therefore be understood as a gradual increases in fraction of the plume-emitting ('turbulent') regions of the boundary layer.

preprint2015arXiv

Salinity transfer in bounded double diffusive convection

The double diffusive convection between two parallel plates is numerically studied for a series of parameters. The flow is driven by the salinity difference and stabilized by the thermal field. Our simulations are directly compared to experiments by Hage and Tilgner (\emph{Phys. Fluids} 22, 076603 (2010)) for several sets of parameters and reasonable agreement is found. This in particular holds for the salinity flux and its dependence on the salinity Rayleigh number. Salt fingers are present in all simulations and extend through the entire height. The thermal Rayleigh number seems to have minor influence on salinity flux but affects the Reynolds number and the morphology of the flow. Next to the numerical calculation, we apply the Grossmann-Lohse theory for Rayleigh-Bénard flow to the current problem without introducing any new coefficients. The theory successfully predicts the salinity flux both with respect to the scaling and even with respect to the absolute value for the numerical and experimental results.