High-resolution Numerical Simulations Research Articles

High-resolution observations and numerical simulations provide insight into solar atmospheric heating

High-resolution observations and numerical simulations provide insight into solar atmospheric heating

On the Sensitivity of Large-Eddy Simulations of the Atmospheric Boundary Layer Coupled with Realistic Large-Scale Dynamics

Abstract We present a new ensemble of 36 numerical experiments aimed at comprehensively gauging the sensitivity of nested large-eddy simulations (LES) driven by large-scale dynamics. Specifically, we explore 36 multiscale configurations of the Weather Research and Forecasting (WRF) Model to simulate the boundary layer flow over the complex topography at the Perdigão field site, with five nested domains discretized at horizontal resolutions ranging from 11.25 km to 30 m. Each ensemble member has a unique combination of the following input factors: (i) large-scale initial and boundary conditions, (ii) subgrid turbulence modeling in the gray zone of turbulence, (iii) subgrid-scale (SGS) models in LES, and (iv) topography and land-cover datasets. We probe their relative importance for LES calculations of velocity, temperature, and moisture fields. Variance decomposition analysis unravels large sensitivities to topography and land-use datasets and very weak sensitivity to the LES SGS model. Discrepancies within ensemble members can be as large as 2.5 m s−1 for the time-averaged near-surface wind speed on the ridge and as large as 10 m s−1 without time averaging. At specific time points, a large fraction of this sensitivity can be explained by the different turbulence models in the gray zone domains. We implement a horizontal momentum and moisture budget routine in WRF to further elucidate the mechanisms behind the observed sensitivity, paving the way for an increased understanding of the tangible effects of the gray zone of turbulence problem. Significance Statement Several science and engineering applications, including wind turbine siting and operations, weather prediction, and downscaling of climate projections, call for high-resolution numerical simulations of the lowest part of the atmosphere. Recent studies have highlighted that such high-resolution simulations, coupled with large-scale models, are challenging and require several important assumptions. With a new set of numerical experiments, we evaluate and compare the significance of different assumptions and outstanding challenges in multiscale modeling (i.e., coupling large-scale models and high-resolution atmospheric simulations). The ultimate goal of this analysis is to put each individual assumption into the wider perspective of a realistic problem and quantify its relative importance compared to other important modeling choices.

A Numerical Study of Clear-Air Turbulence over North China on 6 June 2017

On 6 June 2017, four severe clear-air turbulence (CAT) events were observed over northern China within 3 h. These events mainly occurred at altitudes between 8.1 and 9.5 km. The characteristics and possible mechanisms of the CAT events in the different regions are investigated here using the weather research and forecasting (WRF) model. The simulated wind and temperature fields in a 27 km coarse domain were found to be in good agreement with those of the ERA5 (European Centre for Medium-Range Weather Forecasts Reanalysis v5) and the observed soundings of operational radiosondes over northern China. In terms of synoptic features, the region where the turbulence occurred is characterized by a southwest–northeast upper-level jet stream. The upper-level jet stream observed at an altitude of 10.4 km consistently moved eastwards, with a maximum wind speed of 61.7 m/s. Simultaneously, the upper-level front–jet system on the cyclonic shear side of the upper-level jet stream also exhibited an eastward motion. The developed upper-level front–jet system induced significant vertical wind shear (VWS) and tropopause folding in the vicinity of these CAT events. Despite the high stability resulting from tropopause folding, the presence of strong VWS (1.90 × 10−2 s−1–2.55 × 10−2 s−1) led to a low Richardson number (Ri) (0.24–0.88) and caused Kelvin–Helmholtz instability (KHI), which ultimately induced CAT. Although a standard numerical weather forecast resolution of tens of kilometers is adequate to capture turbulence for most CAT events, it is still necessary to use high-resolution numerical simulations (such as 3 km) to calculate more accurate CAT indices (such as Ri) for CAT prediction in some specific cases.

Trade Off between Hydrodynamic and Thermodynamic Forces at the Liquid-Liquid Interface.

Viscous fingering (VF) instability has been investigated in the case of a partially miscible binary system by nonlinear numerical simulations. Partially miscible fluid systems offer the possibility of phase separation coupled with VF instability. The thermodynamics of such systems are governed by the Margules parameter (interaction parameter) as well as the fluid concentrations. Kinetics of the decomposition is also influenced by dynamical parameters such as the viscosity of the fluid, which incidentally also affects the hydrodynamic forces. Here, we explore the effects of concentration and Margules parameter in order to ascertain the trade-offs incurred between hydrodynamic and thermodynamic effects at the interface as well as the thermodynamics of the bulk. Based on the Gibb's free energy versus concentration curve, we select concentrations (i) outside spinodal and binodal regions, (ii) within binodal but outside the spinodal, and (iii) within the spinodal curve. We solve the modified Cahn-Hilliard-Hele-Shaw equation employing the COMSOL Multiphysics software. Applying high-resolution numerical simulations, we show a strong dependence of the thermodynamic forces on the concentration of the mixtures. Rapid phase separation and hence a faster rate of droplet formation have been found when the concentration lies inside the spinodal region. Further, we have investigated the correlation between the fractal dimension and dynamics of the system. The spatiotemporal studies presented in this work clearly illustrate the competition between hydrodynamic and thermodynamic forces and provide insights on the kinetics of decomposition and growth of interfacial instabilities.

Amplification of Supersonic Microjets by Resonant Inertial Cavitation-Bubble Pair.

We reveal for the first time by experiments that within a narrow parameter regime, two cavitation bubbles with identical energy generated in antiphase develop a supersonic jet. High-resolution numerical simulation shows a mechanism for jet amplification based on toroidal shock wave and bubble necking interaction. The microjet reaches velocities in excess of 1000 m s^{-1}. We demonstrate that potential flow theory established for Worthington jets accurately predicts the evolution of the bubble gas-liquid interfaces unifying compressible and incompressible jet amplification.

Effects of Orography on the High-Temperature Event on 22 June 2023 in North China

An extreme high-temperature event occurred in North China on 22 June 2023, with the maximum temperature reaching 41.8 °C. The high-temperature centers preferentially occurred at the foothills of the Taihang and Yanshan Mountains, indicating an important role of the underlying orography. In the present work, we study the orographic effects of this extreme high-temperature event according to high-resolution numerical simulations using the Weather Research and Forecasting model. The results show that the presence of the mountains in North China contributed notably to the high-temperature event, which can enhance the 2 m air temperature by up to 3 °C. In the daytime, the enhancement of temperature is primarily due to the diabatic heating of sensible heat flux at the terrain surface caused by solar shortwave radiation, whereas the well-known foehn effect has little contribution. Indeed, the dynamically forced downslope flow of foehn is totally suppressed by the upslope flow of the thermally driven mountain-plain circulation. In the nighttime, the sensible heat flux at the terrain surface changes to weakly negative, given the cooling of land surface longwave radiation. As a result, the enhancement of near-surface temperature at the terrain foothill is dominated by the adiabatic warming of downslope flow. Yet, the near-surface temperature far away from the mountain is enhanced by the subsidence warming of a synoptic anomalous anti-cyclone, which is induced by the diabatic heating over the mountains in the daytime. These findings help improve the understanding of the thermal and dynamical effects of orography on the occurrence of high-temperature events.

DESIGN OF THE POSITION OF THE BLOCKING MATERIAL IN THE WATERFLOODED HIGH-PERMEABLE INTERLAYER OF OIL RESERVOIR FOR FIVE-SPOT FLOODING SCHEME

The comparative analysis of the effectiveness of various arrangements of blocking material in a thin, high-permeable water-cut layer of an oil reservoir is presented in order to reduce unproductive injection and in- crease oil recovery. Efficiency calculation was performed using high-resolution numerical simulation in the vertical section of the typical stream tube for a five-spot waterflooding scheme. Three types of location of isolation intervals are considered for two representative ratios of the viscosity of the water and oil phases.

Numerical simulation of droplet dispersion within meso-porous membranes

Analysis of membrane processes in fluid processing, and their main influencing operating conditions are relevant in a variety of industrial applications. Increasing regulatory scrutiny and environmental considerations are forcing industries across all sectors, from food and pharma to oil and gas, to further understand and optimise the handling and formulation of liquid systems for efficient process design. In a generic setup for emulsification and liquid formulation the flow and dispersion behaviour of a liquid oil droplet on its way through a porous water filled membrane is analysed. A set of high-resolution numerical simulations of a single oil droplet dispersed in water through a porous membrane structure with varying contact angles is performed. In this work cluster analysis of volume-of-fluid simulation results to obtain statistical droplet size distributions is conducted and further analysed to highlight the effect of the contact angle as well as pressure drop on the dynamics of the system. It is observed that based on the membrane surface activity the droplet behaviour changes from filtration with coalescence when the membrane is lipophilic to emulsification with droplet break-up when the membrane is lipophobic. Furthermore, the pressure drop is identified as a key factor for the dynamics of the droplet process and the frame in which it occurs. These results highlight that the membrane wettability is a determining factor for the emulsification or filtration effectiveness of a membrane for various applications.

Stellar feedback in the star formation–gas density relation: Comparison between simulations and observations

Context. The impact of stellar feedback on the Kennicutt–Schmidt (KS) law, which relates the star formation rate (SFR) to the surface gas density, is a topic of ongoing debate. The interpretation of high-resolution observations of individual clouds is challenging due to the various processes at play simultaneously and inherent biases. Therefore, a numerical investigation is necessary to understand the role of stellar feedback and identify observable signatures. Aims. In this study we investigate the impact of stellar feedback on the KS law, aiming to identify distinct signatures that can be observed and analysed. By employing magnetohydrodynamic simulations of an isolated cloud, we specifically isolate the effects of high-mass star radiation feedback and protostellar jets. High-resolution numerical simulations are a valuable tool for isolating the impact of stellar feedback on the star formation process, while also allowing us to assess how observational biases may affect the derived relation. Methods. We used high-resolution (<0.01 pc) magnetohydrodynamic numerical simulations of a 104 M⊙ cloud and followed its evolution under different feedback prescriptions. The set of simulations contained four types of feedback: one with only protostellar jets, one with ionising radiation from massive stars (>8 M⊙), one with the combination of the two, and one without any stellar feedback. In order to compare these simulations with the existing observational results, we analysed their evolution by adopting the same techniques applied in the observational studies. Then, we simulated how the same analyses would change if the data were affected by typical observational biases: counting young stellar objects (YSO) to estimate the SFR, the limited resolution for the column density maps, and a sensitivity threshold for detecting faint embedded YSOs. Results. Our analysis reveals that the presence of stellar feedback strongly influences the shape of the KS relation and the star formation efficiency per free-fall time (ϵff). The impact of feedback on the relation is primarily governed by its influence on the cloud’s structure. Additionally, the evolution of ϵff throughout the star formation event suggests that variations in this quantity can mask the impact of feedback in observational studies that do not account for the evolutionary stage of the clouds. Although the ϵff measured in our clouds is higher than what is usually observed in real clouds, upon applying prescriptions to mimic observational biases we recover a good agreement with the expected values. From that, we can infer that observations tend to underestimate the total SFR. Moreover, this likely indicates that the physics included in our simulations is sufficient to reproduce the basic mechanisms that contribute to setting ϵff. Conclusions. We demonstrate the interest of employing numerical simulations to address the impact of early feedback on star formation laws and to correctly interpret observational data. This study will be extended to other types of molecular clouds and ionising stars, sampling different feedback strengths, to fully characterise the impact of H II regions on star formation.

Generalizability of transformer-based deep learning for multidimensional turbulent flow data

Deep learning has been going through rapid advancement and becoming useful in scientific computation, with many opportunities to be applied to various fields, including but not limited to fluid flows and fluid–structure interactions. High-resolution numerical simulations are computationally expensive, while experiments are equally demanding and encompass instrumentation constraints for obtaining flow, acoustics and structural data, particularly at high flow speeds. This paper presents a Transformer-based deep learning method for turbulent flow time series data. Turbulent signals across spatiotemporal and geometrical variations are investigated. The pressure signals are coarsely-grained, and the Transformer creates a fine-grained pressure signal. The training includes data across spatial locations of compliant panels with static deformations arising from the aeroelastic effects of shock-boundary layer interaction. Different training approaches using the Transformer were investigated. Evaluations were carried out using the predicted pressure signal and their power spectra. The Transformer's predicted signals show promising performance. The proposed method is not limited to pressure fluctuations and can be extended to other turbulent or turbulent-like signals.

Low-Order Moments of Velocity Gradient Tensors in Two-Dimensional Isotropic Turbulence

In isotropic turbulence, symmetry of different directions can reduce the number of independent components for velocity gradient tensors. In three-dimensional isotropic turbulence, the independent components under either incompressible or compressible conditions have already been analyzed in the literature. However, for two-dimensional isotropic turbulence, they are still unclear. We derive rigorously the independent components for velocity gradient tensors of two-dimensional isotropic turbulence and give physical explanations. These theoretical results are validated using high-resolution direct numerical simulations (DNSs) of two-dimensional compressible turbulence. Results show that the present DNS setup is still not sufficient to capture the isotropy of third-order moments, suggesting that more investigations on determining the smallest scale and improving the numerical schemes for two-dimensional compressible turbulence are required.

Planetary Boundary Layer Flow over Complex Terrain during a Cold Surge Event: A Case Study

Planetary boundary layer (PBL) flow over complex terrain during a cold surge event was investigated using 3-hourly radiosonde measurements in upwind, near ridge, and downwind locations of mountains in the northeastern part of Republic of Korea and using a high-resolution (333-m) numerical simulation. A cold surge occurred on 23 January 2018 and lasted for 4 days. We analyzed onset day of the cold surge when air temperature dropped rapidly. Analysis of the radiosonde data shows that the PBL was characterized by an adiabatic layer with strong capping inversion in early morning and evening as well as during daytime in the upwind and near-ridge sites. The PBL flow at the near-ridge site was strongest among three sites except at 0600 local standard time (LST) when the PBL flow in the lee was strongest. We performed high-resolution (333-m) numerical simulations using the Weather Research and Forecasting (WRF) model. The adiabatic PBL in the upwind site at 0600 LST was simulated, although its depth was underestimated. The model reproduced the strong low-level wind at 0600 LST and large wind shear during the daytime in the lee, but it did not capture the exact timing of the large wind shear. The model showed overall good performance in simulating the vertical profile of the virtual potential temperature and wind below 2 km above ground level at the three sites, with a high index of agreement (IOA) except for the wind at 1200 and 1500 LST in the lee. To examine the cause for the different behavior of PBL flow in the lee of mountains between 0600 LST and the daytime, we calculated the Froude number for PBL flow using radiosonde measurements based on reduced gravity shallow water (RGSW) theory. At 0600 LST, the upwind Froude number F0 was close to 1, while during the daytime, it was much lower than 1. The observed lee flow behavior was consistent with the flow regime change of a single layer over an obstacle with changing F0; the flow with a propagating lee jump changes into that with a stationary lee jump with decreasing F0. Numerical simulation shows that the steepening of streamlines of lee-wave field leads to a jump-like structure in the lee of mountains during the daytime.

Simulations of a tornadic supercell event in the south of Spain: Sensitivity to initial and boundary conditions and microphysics parameterizations

A tornadic supercell near the village of Campillos (Malaga province) occurred on 26 August 2019 with several associated damages. In this study, we analyze the sensitivity of the Weather Research and Forecasting (WRF) model to initial and boundary conditions and microphysics parameterizations when reproducing this particular event. A total of 20 cases were simulated with four initial and boundary conditions (ERA5, GDAS/FNL, GFS, HRES-IFS) and five microphysics parameterizations (WDM6, Goddard, Morrison, Thompson, WSM6). Forecasts are quantitatively compared with 10-min observations at 153 weather stations for the following variables: temperature at 2 m (T2), relative humidity at 2 m (RH2), wind speed at 10 m (U10) and accumulated precipitation in 3 h (PCP3h). Similarly, the spatial verification method MODE (Method for Object-Based Diagnostic Evaluation) has been used to evaluate the reproduction of convective storms by high-resolution numerical simulations. Our results indicate that microphysics plays a key role in this event, being more relevant than initial and boundary conditions. In particular, Goddard configuration combined with ERA5 or IFS initial and boundary conditions has shown the best overall performance. Simulation based on Goddard and ERA5 captured features to diagnose the development of severe convection. Therefore, WRF model, with the correct parameter sets, was able to reproduce the environmental conditions for supercell formation. The spatial resolution used in the simulations was appropriate to explicitly generate a supercellular storm, although it was not sufficiently high to resolve the formation of tornadoes. Nevertheless, the model showed signs of an environment favourable to tornadogenesis.

Numerical investigation of mixed-phase turbulence induced by a plunging jet

In nature and engineering applications, water jet plunging acts as a key process causing interface breaking and generating mixed-phase turbulence. In this paper, high-resolution numerical simulations of the plunging of a water jet into a quiescent pool were performed to investigate the statistical properties of mixed-phase turbulence, with a special focus on the closure problem of the Reynolds-averaged equation. We conducted phase-resolved simulations, with the air–water interface captured using a coupled level-set and volume-of-fluid method. Various cases were performed to analyse the effects of the Froude number and Reynolds number. The simulation results showed that the turbulence statistics are insensitive to the Reynolds number under investigation, while the Froude number influences the flow properties significantly. To investigate the closure problem of the mean momentum equation, the turbulent kinetic energy (TKE) and turbulent mass flux (TMF) and their transport equations were analysed further. It was discovered that the balance relationship of the TKE budget terms remained similar to many single-phase turbulent flows. The TMF is an additional unclosed term in mixed-phase turbulence over the single-phase turbulence. Our simulation results showed that the production term in its transport equation was highly correlated to TKE. Based on this finding, a closure model for the production term of TMF was further proposed.

Numerical investigation of detonation initiation in a modeled rotating detonation engine

In experimental studies, single-wave mode and two counter-rotating wave mode are often observed in rotating detonation combustors. To investigate the mechanism behind different propagation modes, high-resolution numerical simulations of two-dimensional detonation in hydrogen/air mixtures are conducted by solving the reactive Navier–Stokes equations with a detailed chemical mechanism. The numerical results show that the occurrence of the dual-wave detonation propagation mode is positively influenced by an increase in both the channel width and the initial pressure. The dual-wave modes are observed when increasing the channel width, and it is found that the dual-wave modes are caused by increasing the residual premixed gas height near the inner wall. When increasing the initial pressure, the initial peak detonation heat release increases, which leads to the increase in the hot spot intensity formed, and it is found that the dual-wave modes are mainly caused by the interactions between the initial detonation wave and the inner wall. However, the initial equivalent ratio appears to have a relatively minor impact on the detonation propagation mode due to a relatively narrow range variation of physical properties. The peak heat release rate exerts a greater influence on the change of the propagation mode than the induction time does through a wider range test on rotating detonation engines' working condition. Moreover, the velocities and the cell sizes of detonation waves propagating in different directions with different channel widths are also analyzed, revealing that the characteristics of the detonation waves propagating in different directions are nearly the same.

Three-dimensional receptivity of a high-speed blunt cone to different types of freestream disturbances

Receptivity to freestream disturbances is the initial stage of the boundary-layer transition process, which can determine the final path of boundary-layer disturbance triggering transition. At present, researches on the receptivity of two-dimensional boundary layer to disturbances with zero incident angles, are relatively sufficient. In fact, the freestream disturbances often propagate into the boundary layer in the form of non-zero incident angles, resulting in a component of spatial disturbances in the circumferential direction of rotating bodies (such as a cone). It is a receptivity problem with a distinct three-dimensional feature. However, research on this three-dimensional receptivity problem is relatively scarce. The preliminary work only studied the three-dimensional receptivity to low-frequency incident slow acoustic waves. There has not been a systematic study on the three-dimensional receptivity to different types of freestream disturbances. The paper conducts a study on the three-dimensional receptivity of a blunt cone to different freestream disturbances. Firstly, a high-resolution numerical simulation method is used to compute the three-dimensional receptivity process by introduce freestream disturbances with 15 degree incident angles. The freestream disturbances include fast acoustic wave, slow acoustic wave, entropy wave, and vortex wave. Their frequencies are selected as dimensionless 1.1 and 5, corresponding to the frequencies of the first mode and second mode, respectively. Then, the phase velocity and shape function of the boundary-layer disturbances at each circumferential position for the numerical results are obtained by Fourier transform. To interpret the receptivity mechanisms, the corresponding results by linear stability analysis are obtained for comparisons. It is found that, the first and second modes of the boundary layer can be effectively excited by the incident slow acoustic waves; It is difficult for the incident fast acoustic waves to excite unstable modes in the boundary layer; The incident entropy wave and vortex wave are difficult to excite the first mode at low frequency, but can excite the second mode at high frequency. Furthermore, the incident angle of the freestream disturbances can cause the differences in the receptivity at different circumferential positions of the cone, which can be reflected in two ways. One is the difference in the dominant disturbance form at different circumferential positions; The second is the difference in the amplitude of boundary-layer disturbances. Under different disturbance types and frequencies, these differences between different circumferential positions exhibit different results. The strongest receptivity may occur on the incident front, the incident back, and the incident side. These phenomena may be the result of the combined action of the upstream head disturbances and the disturbances on the incident front.

The geodynamic origin of Los Humeros volcanic field in Mexico: insights from numerical simulations

Compared to normal arc-related volcanic eruptions, the formation of a volcanic caldera is a relatively atypical event. During caldera formation a series of large volumes of magma are erupted, reducing the structural support for the rock above the magma chamber and creating a large depression at the surface called caldera. Los Humeros volcanic field (LHVF) represents one of the largest volcanic calderas in Mexico. It is located some 400 km from the trench at the eastern edge of the Trans Mexican Volcanic Belt where the depth to the Cocos slab is more than 300 km. In this study we employ high-resolution two-dimensional thermomechanical numerical simulations of magma intrusions and a horizontal tectonic strain rate to better understand the influence of crustal deformation for the formation of Los Humeros caldera. A minimum number of three thermal anomaly pulses of hydrated mantle material (with diameter of 15 km or more) and a regional strain rate of 7.927 × 10–16 s−1 are required for magma to reach the surface. Modeling results show that regional extension coupled with deep thermal anomalies (with a temperature excess of ΔT ≥ 100 °C) that come in a specific chain-type sequence produce surface deformation patterns similar to LHVF. We propose an asthenospheric sub-slab deep source (> 300 km depth) for the thermal anomalies where previous studies showed the existence of a gap or tear in the Cocos slab.

Vertical Exchange Induced by Mixed Layer Instabilities

Abstract Submesoscale turbulence in the upper ocean consists of fronts, filaments, and vortices that have horizontal scales on the order of 100 m to 10 km. High-resolution numerical simulations have suggested that submesoscale turbulence is associated with strong vertical motion that could substantially enhance the vertical exchange between the thermocline and mixed layer, which may have an impact on marine ecosystems and climate. Theoretical, numerical, and observational work indicates that submesoscale turbulence is energized primarily by baroclinic instability in the mixed layer, which is most vigorous in winter. This study demonstrates how such mixed layer baroclinic instabilities induce vertical exchange by drawing filaments of thermocline water into the mixed layer. A scaling law is proposed for the dependence of the exchange on environmental parameters. Linear stability analysis and nonlinear simulations indicate that the exchange, quantified by how much thermocline water is entrained into the mixed layer, is proportional to the mixed layer depth, is inversely proportional to the Richardson number of the thermocline, and increases with increasing Richardson number of the mixed layer. The results imply that the tracer exchange between the thermocline and mixed layer is more efficient when the mixed layer is thicker, when the mixed layer stratification is stronger, when the lateral buoyancy gradient is stronger, and when the thermocline stratification is weaker. The scaling suggests vigorous exchange between the permanent thermocline and deep mixed layers in winter, especially in mode water formation regions. Significance Statement This study examines how instabilities in the surface layer of the ocean bring interior water up from below. This interior–surface exchange can be important for dissolved gases such as carbon dioxide and oxygen as well as nutrients fueling biological growth in the surface ocean. A scaling law is proposed for the dependence of the exchange on environmental parameters. The results of this study imply that the exchange is particularly strong if the well-mixed surface layer is thick, lateral density gradients are strong (such as at fronts), and the stratification below the surface layer is weak. These theoretical findings can be implemented in boundary layer parameterization schemes in global ocean models and improve our understanding of the marine ecosystem and how the ocean mediates climate change.

Stars Bisected by Relativistic Blades

We consider the dynamics of an equatorial explosion powered by a millisecond magnetar formed from the core collapse of a massive star. We study whether these outflows—generated by a priori magneto-centrifugally driven, relativistic magnetar winds—might be powerful enough to produce an ultrarelativistic blade (“lamina”) that successfully carves its way through the dense stellar interior. We present high-resolution numerical special-relativistic hydrodynamic simulations of axisymmetric centrifugally driven explosions inside a star and follow the blast wave propagation just after breakout. We estimate the engine requirements to produce ultrarelativistic lamina jets and comment on the physicality of the parameters considered. We find that sufficiently collimated—half-opening angle θ r ≤ 0.°2—laminas successfully break out of a compact progenitor at ultrarelativistic velocities () and extreme isotropic energies (E k,iso ∼ 5 × 1052 erg) within a few percent of the typical spin-down period for a millisecond magnetar. The various phases of these ultrathin outflows, such as collimation shocks, Kelvin–Helmholtz instabilities, and lifetime, are discussed, and we speculate on the observational signatures echoed by this outflow geometry.

The birth and early evolution of a low-mass protostar

Context. Understanding the collapse of dense molecular cloud cores to stellar densities and the subsequent evolution of the protostar is of importance to model the feedback effects such an object has on its surrounding environment, as well as describing the conditions with which it enters the stellar evolutionary track. This process is fundamentally multi-scale, both in density and in spatial extent, and requires the inclusion of complex physical processes such as self-gravity, turbulence, radiative transfer, and magnetic fields. As such, it necessitates the use of robust numerical simulations. Aims. We aim to model the birth and early evolution of a low-mass protostar. We also seek to describe the interior structure of the protostar and the radiative behavior of its accretion shock front. Methods. We carried out a high resolution numerical simulation of the collapse of a gravitationally unstable 1 M⊙ dense molecular cloud core to stellar densities using 3D radiation hydrodynamics under the gray flux-limited diffusion approximation. We followed the initial isothermal phase, the first adiabatic contraction, the second gravitational collapse triggered by the dissociation of H2 molecules, and ≈247 days of the subsequent main accretion phase. Results. We find that the subcritical radiative behavior of the protostar’s shock front causes it to swell as it accretes matter. We also find that the protostar is turbulent from the moment of its inception despite its radiative stability. This turbulence causes significant entropy mixing inside the protostar, which regulates the swelling. Furthermore, we find that the protostar is not fully ionized at birth, but the relative amount of ionized material within it increases as it accretes matter from its surroundings. Finally, we report in the appendix the results of the first 3D calculations involving a frequency-dependent treatment of radiative transfer, which has not produced any major differences with its gray counterpart.