Large-scale Structure and Turbulence Transport in the Young Solar Wind – Comparison of Parker Solar Probe Observations with a Global 3D Reynolds-averaged MHD Model

Author(s):  
Rohit Chhiber ◽  
Arcadi Usmanov ◽  
William Matthaeus ◽  
Melvyn Goldstein ◽  
Riddhi Bandyopadhyay

<div>Simulation results from a global <span>magnetohydrodynamic</span> model of the solar corona and the solar wind are compared with Parker Solar <span>Probe's</span> (<span>PSP</span>) observations during its first several orbits. The fully three-dimensional model (<span>Usmanov</span> <span>et</span> <span>al</span>., 2018, <span>ApJ</span>, 865, 25) is based on Reynolds-averaged mean-flow equations coupled with turbulence transport equations. The model accounts for effects of electron heat conduction, Coulomb collisions, Reynolds stresses, and heating of protons and electrons via nonlinear turbulent cascade. Turbulence transport equations for turbulence energy, cross <span>helicity</span>, and correlation length are solved concurrently with the mean-flow equations. We specify boundary conditions at the coronal base using solar synoptic <span>magnetograms</span> and calculate plasma, magnetic field, and turbulence parameters along the <span>PSP</span> trajectory. We also accumulate data from all orbits considered, to obtain the trends observed as a function of heliocentric distance. Comparison of simulation results with <span>PSP</span> data show general agreement. Finally, we generate synthetic fluctuations constrained by the local rms turbulence amplitude given by the model, and compare properties of this synthetic turbulence with PSP observations.</div>

2021 ◽  
Vol 923 (1) ◽  
pp. 89
Author(s):  
Rohit Chhiber ◽  
Arcadi V. Usmanov ◽  
William H. Matthaeus ◽  
Melvyn L. Goldstein

Abstract Simulation results from a global magnetohydrodynamic model of the solar corona and solar wind are compared with Parker Solar Probe (PSP) observations during its first five orbits. The fully three-dimensional model is based on Reynolds-averaged mean-flow equations coupled with turbulence-transport equations. The model includes the effects of electron heat conduction, Coulomb collisions, turbulent Reynolds stresses, and heating of protons and electrons via a turbulent cascade. Turbulence-transport equations for average turbulence energy, cross helicity, and correlation length are solved concurrently with the mean-flow equations. Boundary conditions at the coronal base are specified using solar synoptic magnetograms. Plasma, magnetic field, and turbulence parameters are calculated along the PSP trajectory. Data from the first five orbits are aggregated to obtain trends as a function of heliocentric distance. Comparison of simulation results with PSP data shows good agreement, especially for mean-flow parameters. Synthetic distributions of magnetic fluctuations are generated, constrained by the local rms turbulence amplitude given by the model. Properties of this computed turbulence are compared with PSP observations.


Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 368
Author(s):  
Laxman Adhikari ◽  
Gary P. Zank ◽  
Lingling Zhao

A detailed study of solar wind turbulence throughout the heliosphere in both the upwind and downwind directions is presented. We use an incompressible magnetohydrodynamic (MHD) turbulence model that includes the effects of electrons, the separation of turbulence energy into proton and electron heating, the electron heat flux, and Coulomb collisions between protons and electrons. We derive expressions for the turbulence cascade rate corresponding to the energy in forward and backward propagating modes, the fluctuating kinetic and magnetic energy, the normalized cross-helicity, and the normalized residual energy, and calculate the turbulence cascade rate from 0.17 to 75 au in the upwind and downwind directions. Finally, we use the turbulence transport models to derive cosmic ray (CR) parallel and perpendicular mean free paths (mfps) in the upwind and downwind heliocentric directions. We find that turbulence in the upwind and downwind directions is different, in part because of the asymmetric distribution of new born pickup ions in the two directions, which results in the CR mfps being different in the two directions. This is important for models that describe the modulation of cosmic rays by the solar wind.


Author(s):  
Yanfei Gao ◽  
Yangwei Liu ◽  
Luyang Zhong ◽  
Jiexuan Hou ◽  
Lipeng Lu

AbstractThe standard k-ε model (SKE) and the Reynolds stress model (RSM) are employed to predict the tip leakage flow (TLF) in a low-speed large-scale axial compressor rotor. Then, a new research method is adopted to “freeze” the turbulent kinetic energy and dissipation rate of the flow field derived from the RSM, and obtain the turbulent viscosity using the Boussinesq hypothesis. The Reynolds stresses and mean flow field computed on the basis of the frozen viscosity are compared with the results of the SKE and the RSM. The flow field in the tip region based on the frozen viscosity is more similar to the results of the RSM than those of the SKE, although certain differences can be observed. This finding indicates that the non-equilibrium turbulence transport nature plays an important role in predicting the TLF, as well as the turbulence anisotropy.


2001 ◽  
Vol 124 (1) ◽  
pp. 86-99 ◽  
Author(s):  
G. A. Gerolymos ◽  
J. Neubauer ◽  
V. C. Sharma ◽  
I. Vallet

In this paper an assessment of the improvement in the prediction of complex turbomachinery flows using a new near-wall Reynolds-stress model is attempted. The turbulence closure used is a near-wall low-turbulence-Reynolds-number Reynolds-stress model, that is independent of the distance-from-the-wall and of the normal-to-the-wall direction. The model takes into account the Coriolis redistribution effect on the Reynolds-stresses. The five mean flow equations and the seven turbulence model equations are solved using an implicit coupled OΔx3 upwind-biased solver. Results are compared with experimental data for three turbomachinery configurations: the NTUA high subsonic annular cascade, the NASA_37 rotor, and the RWTH 1 1/2 stage turbine. A detailed analysis of the flowfield is given. It is seen that the new model that takes into account the Reynolds-stress anisotropy substantially improves the agreement with experimental data, particularily for flows with large separation, while being only 30 percent more expensive than the k−ε model (thanks to an efficient implicit implementation). It is believed that further work on advanced turbulence models will substantially enhance the predictive capability of complex turbulent flows in turbomachinery.


Author(s):  
Xihai Xu ◽  
Xiaodong Li

An anisotropic component of the jet noise source model for the Reynolds-averaged Navier–Stokes equation-based jet noise prediction method is proposed. The modelling is based on Goldstein's generalized acoustic analogy, and both the fine-scale and large-scale turbulent noise sources are considered. To model the anisotropic characteristics of jet noise source, the Reynolds stress tensor is used in place of the turbulent kinetic energy. The Launder–Reece–Rodi model (LRR), combined with Menter's ω -equation for the length scale, with modified coefficients developed by the present authors, is used to calculate the mean flow velocities and Reynolds stresses accurately. Comparison between predicted results and acoustic data has been carried out to verify the accuracy of the new anisotropic source model. This article is part of the theme issue ‘Frontiers of aeroacoustics research: theory, computation and experiment’.


2012 ◽  
Vol 42 (4) ◽  
pp. 539-557 ◽  
Author(s):  
David P. Marshall ◽  
James R. Maddison ◽  
Pavel S. Berloff

Abstract A framework for parameterizing eddy potential vorticity fluxes is developed that is consistent with conservation of energy and momentum while retaining the symmetries of the original eddy flux. The framework involves rewriting the residual-mean eddy force, or equivalently the eddy potential vorticity flux, as the divergence of an eddy stress tensor. A norm of this tensor is bounded by the eddy energy, allowing the components of the stress tensor to be rewritten in terms of the eddy energy and nondimensional parameters describing the mean shape and orientation of the eddies. If a prognostic equation is solved for the eddy energy, the remaining unknowns are nondimensional and bounded in magnitude by unity. Moreover, these nondimensional geometric parameters have strong connections with classical stability theory. When applied to the Eady problem, it is shown that the new framework preserves the functional form of the Eady growth rate for linear instability. Moreover, in the limit in which Reynolds stresses are neglected, the framework reduces to a Gent and McWilliams type of eddy closure where the eddy diffusivity can be interpreted as the form proposed by Visbeck et al. Simulations of three-layer wind-driven gyres are used to diagnose the eddy shape and orientations in fully developed geostrophic turbulence. These fields are found to have large-scale structure that appears related to the structure of the mean flow. The eddy energy sets the magnitude of the eddy stress tensor and hence the eddy potential vorticity fluxes. Possible extensions of the framework to ensure potential vorticity is mixed on average are discussed.


Author(s):  
Sabine Ardey ◽  
Stefan Wolff ◽  
Leonhard Fottner

For a better understanding of the turbulence structures attached to film cooling jets, mean flow velocities and turbulent fluctuations were measured by means of 3D hot wire anemometry in the injection zone of a linear, large scale, high pressure turbine cascade with leading edge film cooling. Near the stagnation point, the blades are equipped with one row of film cooling holes each on the suction and pressure side. Two basically different coolant jet situations are compared: On the suction side the jet features the ordinary kidney vortex. On the pressure side, the jet separates completely from the blade surface since it is located above a large recirculation zone created by a locally adverse pressure gradient and a flow separation near the pressure side injection. The measured Reynolds stresses were analyzed with regard to turbulence production and diffusion. The Bousinesque Hypothesis was tested and could not be confirmed. It was found that the turbulence is highly anisotropic. In addition to the brief description of the experimental set up and the acquired data, given in this paper, the complete information are published as a test case (Ardey and Fottner, 1998) that is directly accessible via internet at: http://www.unibw-muenchen.de/campus/LRT12/welcome.html


2010 ◽  
Vol 644 ◽  
pp. 193-216 ◽  
Author(s):  
IN MEI SOU ◽  
EDWIN A. COWEN ◽  
PHILIP L.-F. LIU

The velocity field and turbulence structure within the surf and swash zones forced by a laboratory-generated plunging breaking wave were investigated using a particle image velocimetry measurement technique. Two-dimensional velocity fields in the vertical plane from 200 consecutive monochromatic waves were measured at four cross-shore locations, shoreward of the breaker line. The phase-averaged mean flow fields indicate that a shear layer occurs when the uprush of the bore front interacts with the downwash flow. The turbulence characteristics were examined via spectral analysis. The larger-scale turbulence structure is closely associated with the breaking-wave- and the bore-generated turbulence in the surf zone; then, the large-scale turbulence energy cascades to smaller scales, as the turbulent kinetic energy (TKE) evolves from the outer surf zone to the swash zone. Smaller-scale energy injection during the latter stage of the downwash phase is associated with the bed-generated turbulence, yielding a −1 slope in the upper inertial range in the spatial spectra. Depth-integrated TKE budget components indicate that a local TKE equilibrium exists during the bore-front phases and the latter stage of the downwash phases in the outer surf zone. The TKE decay resembles the decay of grid turbulence during the latter stage of the uprush and the early stage of the downwash, as the production rate is small because of the absence of strong mean shear during this stage of the wave cycle as well as the relatively short time available for the growth of the bed boundary layer.


1989 ◽  
Vol 111 (1) ◽  
pp. 97-102 ◽  
Author(s):  
D. R. Ballal ◽  
T. H. Chen ◽  
W. J. Schmoll

Turbulence measurements were performed on a 45 deg conical flame stabilizer with a 31 percent blockage ratio, mounted coaxially at the mouth of a circular pipe and supplied with a turbulent premixed methane-air mixture at a Reynolds number of 2.85 × 104. A two-component LDA system was used in the measurement of mean velocities, turbulence intensities, Reynolds stresses, skewness, and kurtosis. It was found that combustion accelerates mean-flow velocities but damps turbulence intensity via the processes of turbulent dilatation and viscous dissipation due to heat release. Measurements in the axial direction showed that the length of the recirculation zone was nearly doubled as a result of combustion. Also, the region around the downstream stagnation point where streamlines meet and velocities change direction was found to be highly turbulent. Skewness and kurtosis data indicated that large-scale eddies carrying fresh combustible mixture are entrained into the high-shear region surrounding the recirculation zone. Finally, a discussion of turbulence-combustion interaction is presented to explain these experimental results.


2022 ◽  
Vol 924 (2) ◽  
pp. 41
Author(s):  
Bernard J. Vasquez ◽  
Sergei A. Markovskii ◽  
Charles W. Smith

Abstract Three-dimensional hybrid kinetic simulations are conducted with particle protons and warm fluid electrons. Alfvénic fluctuations initialized at large scales and with wavevectors that are highly oblique with respect to the background magnetic field evolve into a turbulent energy cascade that dissipates at proton kinetic scales. Accompanying the proton scales is a spectral magnetic helicity signature with a peak in magnitude. A series of simulation runs are made with different large-scale cross helicity and different initial fluctuation phases and wavevector configurations. From the simulations a so-called total magnetic helicity peak is evaluated by summing contributions at a wavenumber perpendicular to the background magnetic field. The total is then compared with the reduced magnetic helicity calculated along spacecraft-like trajectories through the simulation box. The reduced combines the helicity from different perpendicular wavenumbers and depends on the sampling direction. The total is then the better physical quantity to characterize the turbulence. On average the ratio of reduced to total is 0.45. The total magnetic helicity and the reduced magnetic helicity show intrinsic variability based on initial fluctuation conditions. This variability can contribute to the scatter found in the observed distribution of solar wind reduced magnetic helicity as a function of cross helicity.


Sign in / Sign up

Export Citation Format

Share Document