Convection Model Research Articles

The impact of Pangean subducted oceans on mantle dynamics: Passive piles and the positioning of deep mantle plumes

Seismic imaging of the Earth’s interior reveals plumes originating from relatively hot regions of the lowermost mantle, surrounded by cooler material thought to be remnants of ancient subducted oceans. Currently, there is no clear consensus on the internal composition of the hot regions, with end-member conditions being that they are thermo-chemical in nature or purely thermal plume clusters. Previous modelling studies have shown a range of scenarios where deep chemical heterogenities or purely thermal anomalies are essential in developing appropriate present-day mantle dynamics. Here, we add to this discus- sion by quantifying the location of rising mantle plumes using numerical 3-D global mantle convection models constrained by 410 million years of palaeo-ocean evolution (encompassing the formation and breakup of supercontinent Pangea). Our study compares numerical simulations with purely thermal convection to those where a deep thermo-chemical anomaly is laterally mobile. The results show that models both with and without large-scale chemical heterogeneities can generate appropriate present-day plume dynamics, which illustrate the power of sinking ocean plates to stir mantle ow and control the thermal evolution of the mantle. Our models add to the discussion on bottom-up and top-down mantle dynamics, indicating the difficulty in unravelling the processes using numerical models alone.

Numerical Simulation and Neural Network Model for Hydromagnetic Nanofluid Convection in a Porous Wavy Channel with Thermal Non-Equilibrium Model

Abstract The present study aims to analyze the nonfluid MHD convective heat transfer in a porous wavy channel with a local thermal non-equilibrium (LTNE) model. Such a model finds applications related to enhancement in thermal performance, increasing the heat transfer coefficient in the compact design of heat exchangers for the aerospace and automotive industries and elevation in the efficiency of the solar collector. A sinusoidal porous wavy LTNE channel containing nanofluid and subjected to the induced and applied magnetic fields is considered. A uniform magnetic field is applied orthogonal to the channel and the induced magnetic field effects are considered due to the large magnetic Reynolds number. The momentum, continuity, energy, and nanoparticle volume fraction equations constitute the coupled nonlinear system of differential equations and are solved using the Galerkin finite element method. The reliability of the technique is assessed by comparing the proposed procedure with the results from earlier sources. A detailed analysis is presented to determine the effects of different physical parameters arising in the system on temperature, nanoparticle concentration, and velocity profiles. As an illustration, the findings exhibit that increasing the modified diffusivity ratio increases the values of the nanoparticle volume fraction whereas, reducing the modified diffusivity ratio enhances the temperature distribution. A higher value of thermal Rayleigh number presents a significant involvement of buoyancy forces, potentially resulting in the development of convective currents. A higher Nield number indicates more effective heat transport from the solid surface to the nanofluid. Consequently, there is a minimal thermal difference between the solid surface and the bulk nanofluid. Effective heat transmission enhances the nanofluid ability to absorb heat and generates a more consistent dispersion of temperature inside the fluid. The performance of the designed algorithms of the artificial neural network, namely, the Levenberg – Marquardt algorithm, in the problem under consideration is evaluated and the methodology is found reasonably precise with the matching of order around 6 to 7 decimal places of accuracy.

Toward a Comprehensive Grid of Cepheid Models with MESA. I. Uncertainties of the Evolutionary Tracks of Intermediate-mass Stars

Helium-burning stars, in particular Cepheids, are especially difficult to model, as the choice of free parameters can greatly impact the shape of the blue loops—the part of the evolutionary track at which the instability strip is crossed. Contemporary one-dimensional stellar evolution codes, like Modules for Experiments in Stellar Astrophysics (MESA), come with a large number of free parameters that allow us to model the physical processes in stellar interiors under many assumptions. The uncertainties that arise from this freedom are rarely discussed in the literature despite their impact on the evolution of the model. We calculate a grid of evolutionary models with MESA, varying several controls, like solar mixture of heavy elements, mixing-length theory prescription, nuclear reaction rates, the scheme to determine convective boundaries, atmosphere model, and temporal and spatial resolution, and quantify their impact on age and location of the evolutionary track on the H-R diagram from the main sequence until the end of core helium burning. Our investigation was conducted for a full range of masses and metallicities expected for classical Cepheids. The uncertainties are significant, especially during core helium burning, reaching or exceeding the observational uncertainties of logTeff and logL for detached eclipsing binary systems. For ≥9 M ⊙ models, thin convective shells develop and evolve erratically, not allowing the models to converge. A careful inspection of Kippenhahn diagrams and convergence study is advised for a given mass and metallicity, to assess how severe this problem is and to what extent it may affect the evolution.

Atmospheric-river-induced precipitation in California as simulated by the regionally refined Simple Convective Resolving E3SM Atmosphere Model (SCREAM) Version 0

Abstract. Using the regionally refined mesh (RRM) configuration of the US Department of Energy's Simple Cloud-Resolving Energy Exascale Earth System Model (E3SM) Atmosphere Model (SCREAM), we simulate and evaluate four meteorologically distinct atmospheric river events over California. We test five different RRM configurations, each differing in terms of the areal extent of the refined mesh and the resolution (ranging from 800 m to 3.25 km). We find that SCREAM RRM generally has a good representation of the AR-generated precipitation in CA, even for the control simulation which has a very small 3 km refined patch, and is able to capture the fine-scale regional distributions that are controlled largely by the fine-scale topography of the state. It is found that SCREAM generally has a wet bias over topography, most prominently over the Sierra Nevada mountain range, with a corresponding dry bias on the lee side. We find that refining the resolution beyond 3 km (specifically 1.6 km and 800 m) has virtually no benefit towards reducing systematic precipitation biases but that improvements can be found when increasing the areal extent of the upstream refined mesh. However, these improvements are relatively modest and only realized if the size of the refined mesh is expanded to the scale where employing RRM no longer achieves the substantial cost benefit it was intended for.

Marginal stability analyses for thermochemical convection and its implications for the dynamics of continental lithosphere and core–mantle boundary regions

SUMMARY Significant compositional differences may exist in the lithospheric mantle and above the core–mantle boundary (CMB) relative to the ambient mantle. The intrinsic density differences may affect the development of thermal boundary layer (TBL) instabilities associated with lithospheric delamination and formation of thermochemical plumes. In this study, we explored the instability of two-layer thermochemical fluid using two different techniques: marginal stability analysis with a propagator-matrix method and finite element modelling. We investigated both the instabilities in lithospheric mantle (i.e. lithospheric instability) and the mantle above the CMB (i.e. plume-forming instability) using a background temperature Tbg(z) with the TBL. For lithospheric instability, we found that two-layer fluid with free-slip boundary conditions mainly undergoes the same three different convective modes (i.e. two oscillatory convection modes and one layered convection regime) as that with no-slip boundary conditions reported in Jaupart et al. However, with free-slip boundary conditions, the transitions between these convection modes occur at larger values of buoyancy number B. Free-slip boundary conditions lead to smaller critical Rayleigh number Rac, but larger convective wavelength and oscillation frequency ωc, compared with those with no-slip boundary conditions. Our numerical modelling results demonstrate that Rac and ωc predicted from the classical marginal stability analyses using Tbg(z) with TBL temperature may have significant errors when the oscillatory period is comparable with or larger than the timescale of lithospheric thermal diffusion that causes Tbg(z) to vary with time significantly. In this case, using a more gently sloped background temperature profile ignoring the TBL temperature, the stability analysis predicts more accurate stability conditions, thus presenting an effective remedy to the stability analysis. For plume-forming instability, because of the reduced viscosity in the hot and compositionally dense bottom layer, the transition to the layered convection occurs at significantly smaller B values, and in the oscillatory convection regime, Rac is larger but ωc is smaller, compared with those for lithospheric instability. Finally, our study provides a successful benchmark of numerical models of thermochemical convection by comparing Rac and ωc from numerical models with those from the marginal stability analysis.

Quiet Please: Detrending Radial Velocity Variations from Stellar Activity with a Physically Motivated Spot Model

For solar-type stars, spots and their associated magnetic regions induce radial velocity perturbations through the Doppler rotation signal and the suppression of convective blueshift, collectively known as rotation modulation. We developed the Rotation–Convection (RC) model: a method of detrending and characterizing rotation modulation using only cross–correlation functions or one-dimensional spectra without the need for continuous high-cadence measurements. The RC method uses a simple model for the anomalous radial velocity induced by an active region and has two inputs: stellar flux (or a flux proxy) and the relative radial velocity between strongly and weakly absorbed wavelengths (analogous to the bisector–inverse slope). On NEID solar data (3 month baseline), the RC model lowers the amplitude of rotationally modulated stellar activity to below the meter–per–second level. For the standard star HD 26965, the RC model detrends the activity signal to the meter–per–second level for HARPS, EXPRES, and NEID observations, even though the temporal density and time span of the observations differ by an order of magnitude between the three data sets. In addition to detrending, the RC model also characterizes the rotation–modulation signal. From comparison with the Solar Dynamics Observatory, we confirmed that the model accurately recovers and separates the rotation and convection radial velocity components. We also mapped the amplitude of the rotation and convection perturbations as a function of height within the stellar atmosphere. Probing stellar atmospheres with our revised spot model will fuel future innovations in stellar activity mitigation, enabling robust exoplanet detection.

Revisiting the natural convection effects at ultra-low redox concentration solutions: The influence of viscosity and diameter of wire electrode

Natural convection could arise even at ultra-low redox concentration solutions (1–10 mM). Models such as convection–diffusion layer model and spontaneous convection model have been established to describe this phenomenon. However, the driving forces as well as the parameters that influence this natural convection effects are still not clear. Herein we investigated the effects of viscosity on natural convection by introducing sodium alginate (SA), which enhanced viscosity without changing the diffusion coefficient of the redox couple. Resultantly, it allowed us to obtain the relationship between microscopic flow of solutions and the thickness of natural convection layer. Moreover, wire electrodes with various diameters were also tested to reveal the natural convection effects on mass transfer. An empirical equation was established to describing the influences of solution viscosity and diameter of wire electrode on the thickness of natural convection layer in ultra-low redox concentration solutions.

Thermal convection in a Brinkman–Darcy and Kelvin–Voigt fluid of order one with coupled stresses effect

AbstractOur paper thoroughly examines a model for the thermal convection of specific viscoelastic fluids within a porous medium characterized by Brinkman–Darcy properties. These fluids demonstrate first‐order Kelvin–Voigt behavior as they permeate the medium. This investigation considers the influence of couple stresses on the system. We calculate the critical Rayleigh number required for linear instability theory, focusing on both stationary and oscillatory types of convection. Additionally, we conduct the nonlinear stability theory. In the latter cases, it is demonstrated that oscillatory convection can occur. A detailed investigation is conducted into the characteristics of both stationary and oscillatory convection branches, including the points at which the transition occurs.

A Mixed Convection Model for Estimating the Critical Velocity to Prevent Smoke Backlayering in Tunnels

AbstractA novel mathematical model for the critical ventilation velocity to prevent smoke backlayering in tunnels is presented, addressing limitations of prior approaches. The basis of the model is a rigorous characterisation of the physical processes by the characteristic quantities. Empirical parameters within the new model are determined, to align with results from both full-size and small-scale tunnel experiments. Data from numerical simulations (CFD, Computational Fluid Dynamics), validated by known test data, are then used to estimate the effects of tunnel slope and other parameters on the critical velocity. The model is seen to approximate the critical velocity well, following all trends identified by test data and CFD parameter studies. The empirically calibrated equation permits prediction of the critical velocity beyond the narrow range of tunnel geometries where known results already give an answer. The resulting equation has practical application for tunnel design.

The modelling error in multi-dimensional time-dependent solute transport models

Starting from full-dimensional models of solute transport, we derive and analyze multi-dimensional models of time-dependent convection, diffusion, and exchange in and around pulsating vascular and perivascular networks. These models are widely applicable for modelling transport in vascularized tissue, brain perivascular spaces, vascular plants and similar environments. We show the existence and uniqueness of solutions to both the full- and the multi-dimensional equations under suitable assumptions on the domain velocity. Moreover, we quantify the associated modelling errors by establishing a-priori estimates in evolving Bochner spaces. In particular, we show that the modelling error decreases with the characteristic vessel diameter and thus vanishes for infinitely slender vessels. Numerical tests in idealized geometries corroborate and extend upon our theoretical findings.

The sensitivity of lowermost mantle anisotropy to past mantle convection

It is widely believed that seismic anisotropy in the lowermost mantle is caused by the flow-induced alignment of anisotropic crystals such as post-perovskite. What is unclear, however, is whether the anisotropy observations in the lowermost mantle hold information about past mantle flow, or if they only inform us about the present-day flow field. To investigate this, we compare the general and seismic anisotropy calculated using Earth-like mantle convection models where one has a time-varying flow, and another where the present-day flow is constant throughout time. To do this, we track a post-perovskite polycrystal through the flow fields and calculate texture development using the sampled strain rate and the visco-plastic self-consistent approach. We assume dominant slip on (001) and test the effect of the relative importance of this glide plane over others by using three different plasticity models with different efficiencies at developing texture. We compare the radial anisotropy parameters and the anisotropic components of the elastic tensors produced by the flow field test cases at the same location. We find, under all ease-of-texturing cases, the radial anisotropy is very similar (difference <2%) in the majority of locations and in some regions, the difference can be very large (>10%). The same is true when comparing the elastic tensors directly. Varying the ease-of-texture development in the crystal aggregate suggests that easier-to-texture material may hold a stronger signal from past flow than harder-to-texture material. Our results imply that broad-scale observations of seismic anisotropy such as those from seismic tomography, 1-D estimates and normal mode observations, will be mainly sensitive to present-day flow. Shear-wave splitting measurements, however, could hold information about past mantle flow. In general, mantle memory expressed in anisotropy may be dependent on path length in the post-perovskite stability field. Our work implies that, as knowledge of the exact causative mechanism of lowermost mantle anisotropy develops, we may be able to constrain both present-day and past mantle convection.

Global tropical and extra-tropical tropospheric ozone trends and radiative forcing deduced from satellite and ozonesonde measurements for the period 2005–2020

Tropospheric ozone (TPO) is considered as a "near-term climate forcer”, whose impact on climate depends on its radiative forcing (RF), which is a change in the Earth's energy flux. Here, we use the ground-based and satellite measurements during the period 2005–2020 to deduce the trends of TPO, which is significantly positive in the tropical and extra-tropical northern hemisphere (0.2–0.5 DU/yr) and southern extra-tropics (0.1–0.2 DU/yr). Furthermore, the trends derived using a multiple linear regression model (MLR) also confirm these estimates, which are about 0.05–0.1 DU/yr and the regions with higher trends (>0.06 DU/yr) are statistically significant. We also use a standalone Rapid Radiative Transfer Model coupled with a convective model (Radiative-Convective Model; RCM) to assess the climate forcing of ozone using its vertical profiles from the Modern Era Retrospective Reanalysis (MERRA)-2 reanalysis. The estimated temperature rise due to the radiative forcing of ozone in the tropical troposphere (1000–100 hPa) is about 0.2–0.3 °C for the study period. In brief, there is a positive trend in the tropospheric ozone in the tropics and extra-tropics, which is a great concern for regional warming, public health and ecosystem dynamics.

Mixed convection model for predicting the Nusselt number of oscillating vertical plates

In this study, we developed a mixed convection model for predicting the Nusselt number of periodically oscillating vertical plates. When a hot or cold vertical plate oscillates, a complex mixed convection flow occurs, a combination of natural and forced convection flow caused by temperature difference and oscillation, respectively. Here, to simplify the analysis, the natural convection flow is treated as a forced convection flow with a certain stream velocity. Then, by superimposing this flow on the forced convection, the mixed convection flow is approximated to a single effective forced convection flow. Based on this approximation, two closed-form correlations of the Nusselt number are suggested, each applicable to vertically and horizontally oscillating plates, respectively. To check the validity of the suggested correlations, we constructed an experimental apparatus to measure the heat transfer coefficient of oscillating plates. We validated the correlation applicable to horizontally oscillating plates using experimental data obtained from this apparatus, while the correlation developed for vertically oscillating plates was validated using experimental data obtained by previous researchers. Our model and experimental data are in good agreement within 8 % for horizontally oscillating plates and 13 % for vertically oscillating plates. To understand the influence of the Grashof and Reynolds numbers on the Nusselt number of oscillating plates, contour maps of the Nusselt number were presented with respect to the Grashof and Reynolds numbers, and then analyzed. Moreover, we suggest convection regime maps with respect to the Richardson and Prandtl numbers, which can be used to calculate the convection regime of oscillating plates. With the absence of empirical constants, the current model is applicable to all the Prandtl numbers in the range of Re<6×104 and 104<Gr×Pr<109. For this reason, our developed model can be widely used for heat transfer analysis of oscillating vertical plates under various conditions.

Network- and CFD/CAA-modelling of the high frequency flame response in multi-jet combustors

Low order networks are widely used for linear stability analysis of combustors in the low frequency limit. High frequency stability analysis, however, is limited to cost-intensive numerical or experimental methods, since derivation of analytical solutions is either cumbersome or impossible. The article at hand provides a quasi-two-port network model for the effective modal acoustic pressure and axial velocity normalised with the transverse acoustic field for cylindrical combustors. This network modelling approach includes transfer matrices of acoustic area jumps, ducts for longitudinal, standing and spinning transverse and mixed mode wave propagation. The purely acoustic transfer matrices are validated with a generic non-reactive experiment. On the basis of phase-locked [Formula: see text] images of an engine-similar multi-jet combustor with a forced T1 mode, a locally distributed flame response model is derived, which is reduced to a global flame transfer matrix. A locally resolved convective flame response model is implemented in a numerical model in order to verify the provided theory by the comparison of the analytical and numerical flame transfer matrix for the high-frequency regime.

RCEMIP-II: mock-Walker simulations as phase II of the radiative–convective equilibrium model intercomparison project

Abstract. The radiative–convective equilibrium (RCE) model intercomparison project (RCEMIP) leveraged the simplicity of RCE to focus attention on moist convective processes and their interactions with radiation and circulation across a wide range of model types including cloud-resolving models (CRMs), general circulation models (GCMs), single-column models, global cloud-resolving models, and large-eddy simulations. While several robust results emerged across the spectrum of models that participated in the first phase of RCEMIP (RCEMIP-I), two points that stand out are (1) the strikingly large diversity in simulated climate states and (2) the strong imprint of convective self-aggregation on the climate state. However, the lack of consensus in the structure of self-aggregation and its response to warming is a barrier to understanding. Gaining a deeper understanding of convective aggregation and tropical climate will require reducing the degrees of freedom with which convection can vary. Therefore, we propose phase II of RCEMIP (RCEMIP-II) that utilizes a prescribed sinusoidal sea surface temperature (SST) pattern to provide a constraint on the structure of convection and move one critical step up the model hierarchy. This so-called “mock-Walker” configuration generates features that resemble observed tropical circulations. The specification of the mock-Walker protocol for RCEMIP-II is described, along with example results from one CRM and one GCM. RCEMIP-II will consist of five required simulations: three simulations with the same three mean SSTs as in RCEMIP-I but with an SST gradient and two additional simulations at one of the mean SSTs with different values of the SST gradients. We also test the sensitivity to the imposed SST gradient and the domain size. Under weak SST gradients, unforced self-aggregation emerges across the entire domain, similar to what was found in RCEMIP. As the SST gradient increases, the convective region narrows and is more confined to the warmest SSTs. At warmer mean SSTs and stronger SST gradients, low-frequency variability in the convective aggregation emerges, suggesting that simulations of at least 200 d may be needed to achieve robust equilibrium statistics in this configuration. Simulations with different domain sizes generally have similar mean statistics and convective structures, depending on the value of the SST gradient. The prescribed SST boundary condition is the only difference in the set-up between RCEMIP-II and RCEMIP-I, which enables comparison between the two; however, we also welcome participation in RCEMIP-II from models that did not participate in RCEMIP-I.

Near‐surface wind profiles from numerical model predictions. Part II: Verifications against Australia‐wide surface wind observations

AbstractThe new scheme for deriving the near‐surface wind profiles discussed in Ma (in review) is applied to an Australian Bureau of Meteorology operational convective scale model over various domains. Both the new and conventional schemes' diagnostic 10‐m winds are then verified against Australia‐wide automatic weather station observations. Analyses of bulk statistics reveal that the new scheme's 10‐m wind forecasts have generally better accuracy than the current conventional scheme with a consistent reduction of biases over all domains. A widely recognised diurnal bias pattern of surface wind speed over the land is substantially reduced, and the inclusion of Ekman spiral effect on the 10‐m wind marginally improves statistics of the wind direction during the nighttime. The new scheme introduces no systemic bias, given the histogram of a bulk mean bias is analogiased to a Gaussian distribution, and moves the distribution of diagnostic wind speed closer to that observed.

Research on temperature distribution rule and influencing factors of integrated mechanized stope in potassium salt mine

Managing heat damage presents a significant challenge in deep mechanized stopes. Investigating and analyzing temperature distribution patterns and influential factors within the stope environment holds profound importance for enhancing thermal conditions. This study establishes a coupling model of convection and heat transfer, encompassing air flow, heat transfer, species transport and water evaporation. Additionally, a physical model is developed using the Integrated Mechanized Stope at the Laos Kaiyuan Mining potash salt mine as a representative case study. The accuracy of the numerical model is confirmed through comparison with measured data and numerical simulations. Results from numerical simulations indicate that heat dissipation from electromechanical equipment is the primary cause of stope temperature rise. In the stope stabilized airflow region, temperatures near the bottom plate are lower compared to those near the top plate, with an average temperature difference of 1.09 °C. Variation in key parameters reveals that stope temperature decreases with reduced initial airflow temperatures and surrounding rock temperatures, while it increases with diminished initial air volume. Variance analysis indicates the order of influence of key parameters on stope temperature, with initial airflow temperature exerting the greatest impact, followed by surrounding rock temperature and initial air volume. The optimal parameter scheme for controlling heat damage in the stope includes an initial airflow temperature of 22 °C, a surrounding rock temperature of 35 °C, and an initial air volume of 20 kg/s. It can reduce the average temperature of the stope's main working area to 26.21 °C. This study offers a theoretical foundation and practical reference for managing heat damage in underground mining Integrated Mechanized Stopes.

Development of a Reaction–Diffusion–Convection Model and Its Application to NO Oxidation in Reactors

Development of a Reaction–Diffusion–Convection Model and Its Application to NO Oxidation in Reactors

Ancient Stratified Thermochemical Piles Due To High Intrinsic Viscosity

AbstractThe two Large Low Velocity Provinces (LLVPs) in the lowermost Earth mantle are thought to affect large‐scale heat and material transport, governing mantle evolution. LLVPs have been interpreted as thermochemical piles of recycled oceanic crust (ROC) and/or other dense rock types. However, the role of ROC intrinsic viscosity in pile formation and related effects on mantle evolution remain poorly understood. Using mantle convection models, we show that, while ROC intrinsic density controls pile formation, intrinsic viscosity determines whether piles are internally convecting or stratified. Only high‐viscosity, stratified piles can preserve material over several billions of years. Pile stratification is therefore required to reconcile geochemical evidence for the survival of ancient reservoirs. Compositionally layered piles are also consistent with geophysical observations that point to vertical gradients in LLVP properties. As mineral physics constraints point to low‐viscosity ROC, our results suggest that LLVPs may be partly formed by early basal‐magma‐ocean cumulates.

Experimental and analytical investigation of intensive cooling of local distributed heat sources by means of forced flow of liquid

In this work a two-phase forced flow in channels of small geometric dimensions is analyzed. The Navier-Stokes equation was used to create a theoretical model describing the phenomenon of a evaporation of a flowing vapor and the Maxwell equation was used to describe the flow of a thin liquid layer. Both equations are coupled and overall create an original theoretical model describing the phenomenon of the evaporation of a liquid layer. This model of a liquid forced convection cooling the flat element with distributed internal heat sources is theoretically and experimentally analyzed. A flow of the cooling and evaporating liquid is forced by a suction of the vacuum pump. The cooling process is investigated theoretically for two specific cases with both, a constant and a cyclic heat emission. The impact of the liquid layer thickness on the cooling intensity is analyzed for both cases.