scholarly journals A New Surface Model for β-Plane Turbulence

2007 ◽  
Vol 64 (7) ◽  
pp. 2717-2725 ◽  
Author(s):  
Iordanka N. Panayotova

Abstract This paper introduces a new numerical model for studying wave–turbulence interactions in a continuously stratified rotating flow, having a uniform potential vorticity and a rigid boundary. The meridional variation in the Coriolis parameter (β effect), a channel geometry, and the first-order nonlinear terms in a small Rossby number expansion are included into the surface quasigeostrophic dynamics. The model contains important dynamical characteristics of three-dimensional flows such as advection by ageostrophic winds, and stretching and tilting of relative vorticity. Nevertheless, it has the computational economy of two-dimensional flows. Long-term direct numerical simulations are performed for decaying turbulence arising from random initial conditions. In addition to the formation of steady zonal jets, frequently reported as a possible end state under the β effect, this model exhibits several other realistic physical effects that were lacking in the previously studied β-plane models for wave–turbulence interactions. There is a significant contrast in the spatial and time scales of the formed eddies, resulting in high meridional asymmetry of the flow. Mean surface cooling and persistent large-scale blocking eddies are observed as well. The surface potential temperature variance spectrum exhibits a well-resolved k−5/3 inertial range.

Geophysics ◽  
1993 ◽  
Vol 58 (5) ◽  
pp. 692-702 ◽  
Author(s):  
Peter Hubral ◽  
Jorg Schleicher ◽  
Martin Tygel

Zero‐offset reflections resulting from point sources are often computed on a large scale in three‐dimensional (3-D) laterally inhomogeneous isotropic media with the help of ray theory. The geometrical‐spreading factor and the number of caustics that determine the shape of the reflected pulse are then generally obtained by integrating the so‐called dynamic ray‐tracing system down and up to the two‐way normal incidence ray. Assuming that this ray is already known, we show that one integration of the dynamic ray‐tracing system in a downward direction with only the initial condition of a point source at the earth’s surface is in fact sufficient to obtain both results. To establish the Fresnel zone of the zero‐offset reflection upon the reflector requires the same single downward integration. By performing a second downward integration (using the initial conditions of a plane wave at the earth’s surface) the complete Fresnel volume around the two‐way normal ray can be found. This should be known to ascertain the validity of the computed zero‐offset event. A careful analysis of the problem as performed here shows that round‐trip integrations of the dynamic ray‐tracing system following the actually propagating wavefront along the two‐way normal ray need never be considered. In fact some useful quantities related to the two‐way normal ray (e.g., the normal‐moveout velocity) require only one single integration in one specific direction only. Finally, a two‐point ray tracing for normal rays can be derived from one‐way dynamic ray tracing.


2016 ◽  
Vol 138 (7) ◽  
Author(s):  
B. M. Wilson ◽  
R. Mejia-Alvarez ◽  
K. P. Prestridge

Mach number and initial conditions effects on Richtmyer–Meshkov (RM) mixing are studied by the vertical shock tube (VST) at Los Alamos National Laboratory (LANL). At the VST, a perturbed stable light-to-heavy (air–SF6, A = 0.64) interface is impulsively accelerated with a shock wave to induce RM mixing. We investigate changes to both large and small scales of mixing caused by changing the incident Mach number (Ma = 1.3 and 1.45) and the three-dimensional (3D) perturbations on the interface. Simultaneous density (quantitative planar laser-induced fluorescence (PLIF)) and velocity (particle image velocimetry (PIV)) measurements are used to characterize preshock initial conditions and the dynamic shocked interface. Initial conditions and fluid properties are characterized before shock. Using two types of dynamic measurements, time series (N = 5 realizations at ten locations) and statistics (N = 100 realizations at a single location) of the density and velocity fields, we calculate several mixing quantities. Mix width, density-specific volume correlations, density–vorticity correlations, vorticity, enstrophy, strain, and instantaneous dissipation rate are examined at one downstream location. Results indicate that large-scale mixing, such as the mix width, is strongly dependent on Mach number, whereas small scales are strongly influenced by initial conditions. The enstrophy and strain show focused mixing activity in the spike regions.


2007 ◽  
Vol 64 (12) ◽  
pp. 4417-4431 ◽  
Author(s):  
Chris Snyder ◽  
David J. Muraki ◽  
Riwal Plougonven ◽  
Fuqing Zhang

Abstract Vortex dipoles provide a simple representation of localized atmospheric jets. Numerical simulations of a synoptic-scale dipole in surface potential temperature are considered in a rotating, stratified fluid with approximately uniform potential vorticity. Following an initial period of adjustment, the dipole propagates along a slightly curved trajectory at a nearly steady rate and with a nearly fixed structure for more than 50 days. Downstream from the jet maximum, the flow also contains smaller-scale, upward-propagating inertia–gravity waves that are embedded within and stationary relative to the dipole. The waves form elongated bows along the leading edge of the dipole. Consistent with propagation in horizontal deformation and vertical shear, the waves’ horizontal scale shrinks and the vertical slope varies as they approach the leading stagnation point in the dipole’s flow. Because the waves persist for tens of days despite explicit dissipation in the numerical model that would otherwise damp the waves on a time scale of a few hours, they must be inherent features of the dipole itself, rather than remnants of imbalances in the initial conditions. The wave amplitude varies with the strength of the dipole, with waves becoming obvious once the maximum vertical vorticity in the dipole is roughly half the Coriolis parameter. Possible mechanisms for the wave generation are spontaneous wave emission and the instability of the underlying balanced dipole.


2021 ◽  
pp. 1-53
Author(s):  
Kedar P. Nawathe ◽  
Rui Zhu ◽  
Enci Lin ◽  
Yong Kim ◽  
Terrence Simon

Abstract The stators of the first stage of a gas turbine are exposed to severe temperatures. The coolant streams introduced to prevent the stators from thermal damage further complicate the highly three-dimensional vane passage flow. Recent results have shown that the coolant streams injected for cooling the combustor also influence the flow physics and the cooling effectiveness in the first-stage stator vanes passage. However, the effects of changing the mass flow rate of these combustor coolant streams on the passage flowfield have not been studied. As understanding the coolant transport is necessary for analyzing changes in cooling effectiveness in the vane passage, detailed aerodynamic and thermal measurements along the whole vane passage are required. This two-part paper presents such measurements taken for a variety of combustor coolant and endwall film coolant flow rates. The experiments were conducted in a low-Mach-number facility with engine-representative Reynolds numbers and large-scale high-level turbulence. The objective of the first part is to describe the flow that influences endwall and vane surface cooling effectiveness distributions, which are presented in the second part. The measurements show changes in the passage flowfield due to changes in both combustor coolant and endwall film coolant flow rates. Overall, the flow-physics remains largely unaffected by changes in coolant flow rates except in the endwall-vane surfaces region where the combustor coolant flow rate dominates changes in coolant transport. This is shown to have a high impact on endwall and vane surface cooling.


2016 ◽  
Vol 811 ◽  
pp. 701-741 ◽  
Author(s):  
W. He ◽  
R. S. Gioria ◽  
J. M. Pérez ◽  
V. Theofilis

Two- and three-dimensional modal and non-modal instability mechanisms of steady spanwise-homogeneous laminar separated flow over airfoil profiles, placed at large angles of attack against the oncoming flow, have been investigated using global linear stability theory. Three NACA profiles of distinct thickness and camber were considered in order to assess geometry effects on the laminar–turbulent transition paths discussed. At the conditions investigated, large-scale steady separation occurs, such that Tollmien–Schlichting and cross-flow mechanisms have not been considered. It has been found that the leading modal instability on all three airfoils is that associated with the Kelvin–Helmholtz mechanism, taking the form of the eigenmodes known from analysis of generic bluff bodies. The three-dimensional stationary eigenmode of the two-dimensional laminar separation bubble, associated in earlier analyses with the formation on the airfoil surface of large-scale separation patterns akin to stall cells, is shown to be more strongly damped than the Kelvin–Helmholtz mode at all conditions examined. Non-modal instability analysis reveals the potential of the flows considered to sustain transient growth which becomes stronger with increasing angle of attack and Reynolds number. Optimal initial conditions have been computed and found to be analogous to those on a cascade of low pressure turbine blades. By changing the time horizon of the analysis, these linear optimal initial conditions have been found to evolve into the Kelvin–Helmholtz mode. The time-periodic base flows ensuing linear amplification of the Kelvin–Helmholtz mode have been analysed via temporal Floquet theory. Two amplified modes have been discovered, having characteristic spanwise wavelengths of approximately 0.6 and 2 chord lengths, respectively. Unlike secondary instabilities on the circular cylinder, three-dimensional short-wavelength perturbations are the first to become linearly unstable on all airfoils. Long-wavelength perturbations are quasi-periodic, standing or travelling-wave perturbations that also become unstable as the Reynolds number is further increased. The dominant short-wavelength instability gives rise to spanwise periodic wall-shear patterns, akin to the separation cells encountered on airfoils at low angles of attack and the stall cells found in flight at conditions close to stall. Thickness and camber have quantitative but not qualitative effect on the secondary instability analysis results obtained.


2020 ◽  
Vol 500 (1) ◽  
pp. 1018-1028
Author(s):  
Shogo Masaki ◽  
Takahiro Nishimichi ◽  
Masahiro Takada

ABSTRACT To generate initial conditions for cosmological N-body simulations, one needs to prepare a uniform distribution of simulation particles, the so-called pre-initial condition (pre-IC). The standard method to construct the pre-IC is to place the particles on the lattice grids evenly spaced in the three-dimensional spatial coordinates. However, even after the initial displacement of each particle according to cosmological perturbations, the particle distribution remains to display an artificial anisotropy. Such an artefact causes systematic effects in simulations at later time until the evolved particle distribution sufficiently erases the initial anisotropy. In this paper, we study the impacts of the pre-IC on the anisotropic separate universe simulation, where the effect of large-scale tidal field on structure formation is taken into account using the anisotropic expansion in a local background (simulation volume). To quantify the impacts, we compare the simulations employing the standard grid pre-IC and the glass one, where the latter is supposed to suppress the initial anisotropy. We show that the artificial features in the grid pre-IC simulations are seen until z ∼ 9, while the glass pre-IC simulations appear to be stable and accurate over the range of scales we study. From these results we find that a coupling of the large-scale tidal field with matter clustering is enhanced compared to the leading-order prediction of perturbation theory in the quasi-non-linear regime in the redshift range 5 ≲ z ≲ 15, indicating the importance of tidal field on structure formation at such high redshifts, e.g. during the epoch of reionization.


Author(s):  
Josef Hasslberger ◽  
Peter Katzy ◽  
Thomas Sattelmayer ◽  
Lorenz R. Boeck

For the purpose of nuclear safety analysis, a reactive flow solver has been developed to determine the hazard potential of large-scale hydrogen explosions. Without using empirical transition criteria, the whole combustion process (including DDT) is computed within a single solver framework. In this paper, we present massively parallelized three-dimensional explosion simulations in a full-scale pressurized water reactor of the Konvoi type. Several generic DDT scenarios in globally lean hydrogen-air mixtures are examined to assess the importance of different input parameters. It is demonstrated that the explosion process is highly sensitive to mixture composition, ignition location and thermodynamic initial conditions. Pressure loads on the confining structure show a profoundly dynamic behavior depending on the position in the containment.


2007 ◽  
Vol 585 ◽  
pp. 395-420 ◽  
Author(s):  
P. LAVOIE ◽  
L. DJENIDI ◽  
R. A. ANTONIA

The effects of initial conditions on grid turbulence are investigated for low to moderate Reynolds numbers. Four grid geometries are used to yield variations in initial conditions and a secondary contraction is introduced to improve the isotropy of the turbulence. The hot-wire measurements, believed to be the most detailed to date for this flow, indicate that initial conditions have a persistent impact on the large-scale organization of the flow over the length of the tunnel. The power-law coefficients, determined via an improved method, also depend on the initial conditions. For example, the power-law exponent m is affected by the various levels of large-scale organization and anisotropy generated by the different grids and the shape of the energy spectrum at low wavenumbers. However, the results show that these effects are primarily related to deviations between the turbulence produced in the wind tunnel and true decaying homogenous isotropic turbulence (HIT). Indeed, when isotropy is improved and the intensity of the large-scale periodicity, which is primarily associated with round-rod grids, is decreased, the importance of initial conditions on both the character of the turbulence and m is diminished. However, even in the case where the turbulence is nearly perfectly isotropic, m is not equal to −1, nor does it show an asymptotic trend in x towards this value, as suggested by recent analysis. Furthermore, the evolution of the second- and third-order velocity structure functions satisfies equilibrium similarity only approximately.


2015 ◽  
Vol 27 (4) ◽  
pp. 561-572 ◽  
Author(s):  
Lanming Jin ◽  
Gaoming Jiang

Purpose – Multilayer weft knitted fabrics possess many advantages, such as strongly stereoscopic patterns, soft handling and adjustable thickness of apparel and home textiles use. However, it is difficult to predict the final visual effects before the productive process because of the three-dimensional (3D) effect caused by the connecting yarn of the fabric. The purpose of this paper is to realize a realistic simulation of the fabric. Design/methodology/approach – The authors applied to the curve and surface model to simulate the knitted fabric, instead of previous single loop model by NURBS. Macro simulation is more suitable for the fabric with the 3D effect because of the quick, real and convenient simulation. This research includes experiments on the structural parameters concerning the regular sag of multilayer weft knitted fabrics, and analysis of parameter data and the simulation process with the aim of realizing a computer simulation of the fabric, especially with a sense of reality. The Digital Elevation Model was also applied to build a simulated 3D model. Findings – To obtain the values for the change rules, different samples were used and the outputs of the model were found to be close to the experimental results. The thickest and thinnest lengths and the changing curves between them were established. Patterned simple multilayer weft knitted fabric could be simulated through the results of the research. It is possible to simulate different real fabrics using their materials and expected effects. The authors are going to improve the model to simulate the complicate large-scale jacquard fabrics in further research. Practical implications – The results will be useful for establishing a computer surface simulation system for stereo perception of fabrics. Originality/value – The authors put forward the concept of surface warpage degree (R). It is an important factor affecting the fabric stereo feeling. The larger the value of R, the stronger the stereo sense of the fabric. It could be applied to most 3D fabric. A thickness difference testing method was proposed to characterize the stereo perception of fabrics. It is possible to simulate different real fabrics quickly without the model of the woven loop.


Author(s):  
Tom Dörffel ◽  
Ariane Papke ◽  
Rupert Klein ◽  
Natalia Ernst ◽  
Piotr K. Smolarkiewicz

AbstractPäschke et al. (J Fluid Mech, 2012) studied the nonlinear dynamics of strongly tilted vortices subject to asymmetric diabatic heating by asymptotic methods. They found, inter alia, that an azimuthal Fourier mode 1 heating pattern can intensify or attenuate such a vortex depending on the relative orientation of the tilt and the heating asymmetries. The theory originally addressed the gradient wind regime which, asymptotically speaking, corresponds to vortex Rossby numbers of order unity in the limit. Formally, this restricts the applicability of the theory to rather weak vortices. It is shown below that said theory is, in contrast, uniformly valid for vanishing Coriolis parameter and thus applicable to vortices up to low hurricane strengths. An extended discussion of the asymptotics as regards their physical interpretation and their implications for the overall vortex dynamics is also provided in this context. The paper’s second contribution is a series of three-dimensional numerical simulations examining the effect of different orientations of dipolar diabatic heating on idealized tropical cyclones. Comparisons with numerical solutions of the asymptotic equations yield evidence that supports the original theoretical predictions of Päschke et al. In addition, the influence of asymmetric diabatic heating on the time evolution of the vortex centerline is further analyzed, and a steering mechanism that depends on the orientation of the heating dipole is revealed. Finally, the steering mechanism is traced back to the correlation of dipolar perturbations of potential temperature, induced by the vortex tilt, and vertical velocity, for which diabatic heating not necessarily needs to be responsible, but which may have other origins.


Sign in / Sign up

Export Citation Format

Share Document