scholarly journals Temperature gradient driven heat flux closure in fluid simulations of collisionless reconnection

2018 ◽  
Vol 84 (3) ◽  
Author(s):  
F. Allmann-Rahn ◽  
T. Trost ◽  
R. Grauer

Recent efforts to include kinetic effects in fluid simulations of plasmas have been very promising. Concerning collisionless magnetic reconnection, it has been found before that damping of the pressure tensor to isotropy leads to good agreement with kinetic runs in certain scenarios. An accurate representation of kinetic effects in reconnection was achieved in a study by Wang et al. (Phys. Plasmas, vol. 22, 2015, 012108) with a closure derived from earlier work by Hammett and Perkins (Phys. Rev. Lett., vol. 64, 1990, 3019). Here, their approach is analysed on the basis of heat flux data from a Vlasov simulation. As a result, we propose a new local closure in which heat flux is driven by temperature gradients. This way, a more realistic approximation of Landau damping in the collisionless regime is achieved. Previous issues are addressed and the agreement with kinetic simulations in different reconnection set-ups is improved significantly. To the authors’ knowledge, the new fluid model is the first to perform well in simulations of the coalescence of large magnetic islands.

2008 ◽  
Vol 130 (12) ◽  
Author(s):  
David T. W. Lin ◽  
Hung Yi Li ◽  
Wei Mon Yan

An inverse solution scheme based on the conjugate gradient method with the minimization of the object function is presented for estimating the unknown wall heat flux of conjugated forced convection flows between two corotating disks from temperature measurements acquired within the flow field. The validity of the proposed approach is demonstrated via the estimation of three time- and space-dependent heat flux profiles. A good agreement is observed between the estimated results and the exact solution in every case. In general, the accuracy of the estimated results is found to improve as the temperature sensors are moved closer to the unknown boundary surface and the error in the measured temperature data is reduced.


2000 ◽  
Author(s):  
Christian Proulx ◽  
Daniel R. Rousse ◽  
Rodolphe Vaillon ◽  
Jean-François Sacadura

Abstract This article presents selected results of a study comparing two procedures for the treatment of collimated irradiation impinging on one boundary of a participating one-dimensional plane-parallel medium. These procedures are implemented in a CVFEM used to calculate the radiative heat flux and source. Both isotropically and anisotropically scattering media are considered. The results presented show that both procedures provide results in good agreement with those obtained using a Monte Carlo method, when the collimated beam impinges normally.


Author(s):  
Y. Bouaichaoui ◽  
R. Kibboua ◽  
M. Matkovič

The knowledge of the onset of subcooled boiling in forced convective flow at high liquid velocity and subcooling is of importance in thermal hydraulic studies. Measurements were performed under various conditions of mass flux, heat flux, and inlet subcooling, which enabled to study the influence of different boundary conditions on the development of local flow parameters. Also, some measurements have been compared to the predictions by the three-dimensional two-fluid model of subcooled boiling flow carried out with the computer code ANSYS-CFX-13. A computational method based on theoretical studies of steady state two phase forced convection along a test section loop was released. The calculation model covers a wide range of two phase flow conditions. It predicts the heat transfer rates and transitions points such as the Onset of Critical Heat Flux.


1978 ◽  
Vol 33 (7) ◽  
pp. 749-760 ◽  
Author(s):  
G. E. J. Eggermont ◽  
P. W. Hermans ◽  
L. J. F. Hermans ◽  
H. F. P. Knaap ◽  
J. J. M. Beenakker

In a rarefied polyatomic gas streaming through a rectangular channel, an external magnetic field produces a heat flux perpendicular to the flow direction. Experiments on this “viscom agnetic heat flux” have been performed for CO, N2, CH4 and HD at room temperature, with different orientations of the magnetic field. Such measurements enable one to separate the boundary layer contribution from the purely bulk contribution by means of the theory recently developed by Vestner. Very good agreement is found between the experimentally determined bulk contribution and the theoretical Burnett value for CO, N2 and CH4 , yet the behavior of HD is found to be anomalous.


2020 ◽  
Vol 2020 ◽  
pp. 1-14
Author(s):  
Keltoum Chahour ◽  
Rajae Aboulaich ◽  
Abderrahmane Habbal ◽  
Nejib Zemzemi ◽  
Chérif Abdelkhirane

Fractional flow reserve (FFR) has proved its efficiency in improving patient diagnosis. In this paper, we consider a 2D reconstructed left coronary tree with two artificial lesions of different degrees. We use a generalized fluid model with a Carreau law and use a coupled multidomain method to implement Windkessel boundary conditions at the outlets. We introduce our methodology to quantify the virtual FFR and conduct several numerical experiments. We compare FFR results from the Navier–Stokes model versus generalized flow model and for Windkessel versus traction-free outlet boundary conditions or mixed outlet boundary conditions. We also investigate some sources of uncertainty that the FFR index might encounter during the invasive procedure, in particular, the arbitrary position of the distal sensor. The computational FFR results show that the degree of stenosis is not enough to classify a lesion, while there is a good agreement between the Navier–Stokes model and the non-Newtonian flow model adopted in classifying coronary lesions. Furthermore, we highlight that the lack of standardization while making FFR measurement might be misleading regarding the significance of stenosis.


Processes ◽  
2020 ◽  
Vol 8 (7) ◽  
pp. 810
Author(s):  
Jade Gesare Abuga ◽  
Tiri Chinyoka

The flow of viscoelastic fluids may, under certain conditions, exhibit shear-banding characteristics that result from their susceptibility to unusual flow instabilities. In this work, we explore both the existing shear banding mechanisms in the literature, namely; constitutive instabilities and flow-induced inhomogeneities. Shear banding due to constitutive instabilities is modelled via either the Johnson–Segalman or the Giesekus constitutive models. Shear banding due to flow-induced inhomogeneities is modelled via the Rolie–Poly constitutive model. The Rolie–Poly constitutive equation is especially chosen because it expresses, precisely, the shear rheometry of polymer solutions for a large number of strain rates. For the Rolie–Poly approach, we use the two-fluid model wherein the stress dynamics are coupled with concentration equations. We follow a computational analysis approach via an efficient and versatile numerical algorithm. The numerical algorithm is based on the Finite Volume Method (FVM) and it is implemented in the open-source software package, OpenFOAM. The efficiency of our numerical algorithms is enhanced via two possible stabilization techniques, namely; the Log-Conformation Reformulation (LCR) and the Discrete Elastic Viscous Stress Splitting (DEVSS) methodologies. We demonstrate that our stabilized numerical algorithms accurately simulate these complex (shear banded) flows of complex (viscoelastic) fluids. Verification of the shear-banding results via both the Giesekus and Johnson-Segalman models show good agreement with existing literature using the DEVSS technique. A comparison of the Rolie–Poly two-fluid model results with existing literature for the concentration and velocity profiles is also in good agreement.


2019 ◽  
Vol 864 ◽  
pp. 519-553 ◽  
Author(s):  
Grace A. Cox ◽  
Christopher J. Davies ◽  
Philip W. Livermore ◽  
James Singleton

Motivated by the dynamics within terrestrial bodies, we consider a rotating, strongly thermally stratified fluid within a spherical shell subject to a prescribed laterally inhomogeneous heat-flux condition at the outer boundary. Using a numerical model, we explore a broad range of three key dimensionless numbers: a thermal stratification parameter (the relative size of boundary temperature gradients to imposed vertical temperature gradients), $10^{-3}\leqslant S\leqslant 10^{4}$, a buoyancy parameter (the strength of applied boundary heat-flux anomalies), $10^{-2}\leqslant B\leqslant 10^{6}$, and the Ekman number (ratio of viscous to Coriolis forces), $10^{-6}\leqslant E\leqslant 10^{-4}$. We find both steady and time-dependent solutions and delineate the regime boundaries. We focus on steady-state solutions, for which a clear transition is found between a low $S$ regime, in which buoyancy dominates the dynamics, and a high $S$ regime, in which stratification dominates. For the low-$S$ regime, we find that the characteristic flow speed scales as $B^{2/3}$, whereas for high-$S$, the radial and horizontal velocities scale respectively as $u_{r}\sim S^{-1}$, $u_{h}\sim S^{-3/4}B^{1/4}$ and are confined within a thin layer of depth $(SB)^{-1/4}$ at the outer edge of the domain. For the Earth, if lower mantle heterogeneous structure is due principally to chemical anomalies, we estimate that the core is in the high-$S$ regime and steady flows arising from strong outer boundary thermal anomalies cannot penetrate the stable layer. However, if the mantle heterogeneities are due to thermal anomalies and the heat-flux variation is large, the core will be in a low-$S$ regime in which the stable layer is likely penetrated by boundary-driven flows.


2007 ◽  
Vol 14 (4) ◽  
pp. 373-383 ◽  
Author(s):  
D. Borgogno ◽  
T. Passot ◽  
P. L. Sulem

Abstract. Non-propagating magnetic hole solutions in anisotropic plasmas near the mirror instability threshold are investigated in numerical simulations of a fluid model that incorporates linear Landau damping and finite Larmor radius corrections calculated in the gyrokinetic approximation. This FLR-Landau fluid model reproduces the subcritical mirror bifurcation recently identified on the Vlasov-Maxwell system, both by theory and numerics. Stable magnetic hole solutions that display a polarization different from that of Hall-MHD solitons are indeed obtained slighlty below threshold, while magnetic patterns and spatio-temporal chaos emerge when the system is maintained in a mirror unstable regime.


Sign in / Sign up

Export Citation Format

Share Document