scholarly journals Kicked Burgers turbulence

2000 ◽  
Vol 416 ◽  
pp. 239-267 ◽  
Author(s):  
J. BEC ◽  
U. FRISCH ◽  
K. KHANIN
Keyword(s):  
Numerical Simulations ◽  
Power Law ◽  
Large Scale ◽  
Order Term ◽  
Inviscid Limit ◽  
Order Structure ◽  
Legendre Transform ◽  
Leading Order ◽  
Space And Time ◽  
Burgers Turbulence

Burgers turbulence subject to a force f(x, t) = [sum ]jfj(x)δ(t − tj), where tj are 'kicking times' and the 'impulses' fj(x) have arbitrary space dependence, combines features of the purely decaying and the continuously forced cases. With large-scale forcing this ‘kicked’ Burgers turbulence presents many of the regimes proposed by E et al. (1997) for the case of random white-noise-in-time forcing. It is also amenable to efficient numerical simulations in the inviscid limit, using a modification of the fast Legendre transform method developed for decaying Burgers turbulence by Noullez & Vergassola (1994). For the kicked case, concepts such as ‘minimizers’ and ‘main shock’, which play crucial roles in recent developments for forced Burgers turbulence, become elementary since everything can be constructed from simple two-dimensional area-preserving Euler–Lagrange maps.The main results are for the case of identical deterministic kicks which are periodic and analytic in space and are applied periodically in time. When the space integrals of the initial velocity and of the impulses vanish, it is proved and illustrated numerically that a space- and time-periodic solution is achieved exponentially fast. In this regime, probabilities can be defined by averaging over space and time periods. The probability densities of large negative velocity gradients and of (not-too-large) negative velocity increments follow the power law with −7/2 exponent proposed by E et al. (1997) in the inviscid limit, whose existence is still controversial in the case of white-in-time forcing. This power law, which is seen very clearly in the numerical simulations, is the signature of nascent shocks (preshocks) and holds only when at least one new shock is born between successive kicks.It is shown that the third-order structure function over a spatial separation Δx is analytic in Δx although the velocity field is generally only piecewise analytic (i.e. between shocks). Structure functions of order p ≠ 3 are non-analytic at Δx = 0. For even p there is a leading-order term proportional to [mid ]Δx[mid ] and for odd p > 3 the leading-order term ∝Δx has a non-analytic correction ∝Δx[mid ]Δx[mid ] stemming from shock mergers.

Nonlinear interactions in turbulence with strong irrotational straining

1997 ◽  
Vol 337 ◽  
pp. 333-364 ◽  
Author(s):  
N.K.-R. KEVLAHAN ◽  
J.C.R. HUNT
Keyword(s):  
Large Scale ◽  
Order Term ◽  
Initial Energy ◽  
Perturbation Series ◽  
Small Scale ◽  
Bluff Bodies ◽  
Leading Order ◽  
Viscous Effects ◽  
First Case ◽  
Maximum Amplification

The rate of growth of the nonlinear terms in the vorticity equation are analysed for a turbulent flow with r.m.s. velocity u0 and integral length scale L subjected to a strong uniform irrotational plane strain S, where (u0/L)/S=ε[Lt ]1. The rapid distortion theory (RDT) solution is the zeroth-order term of the perturbation series solution in terms of ε. We use the asymptotic form of the convolution integrals for the leading-order nonlinear terms when β= exp(−St)[Lt ]1 to determine at what time t and beyond what wavenumber k (normalized on L) the perturbation series in ε fails, and hence derive the following conditions for the validity of RDT in these flows. (a) The magnitude of the nonlinear terms of order ε depends sensitively on the amplitude of eddies with large length scales in the direction x2 of negative strain. (b) If the integral of the velocity component u2 is zero the leading-order nonlinear terms increase and decrease in the same way as the linear terms, even those that decrease exponentially. In this case RDT calculations of vorticity spectra become invalid at a time tNL∼L/u0k−3 independent of ε and the power law of the initial energy spectrum, but the calculation of the r.m.s. velocity components by RDT remains accurate until t= TNL∼L/u0, when the maximum amplification of r.m.s. vorticity is ω/S∼εexp(ε−1)[Gt ]1. (c) If this special condition does not apply, the leading-order nonlinear terms increase faster than the linear terms by a factor O(β−1). RDT calculations of the vorticity spectrum then fail at a shorter time tNL∼(1/S) ln(ε−1k−3); in this case TNL∼(1/S) ln(ε−1) and the maximum amplification of r.m.s. vorticity is ω/S∼1. (d) Viscous effects dominate when t[Gt ](1/S) ln(k−1(Re/ε)1/2). In the first case RDT fails immediately in this range, while in the second case RDT usually fails before viscosity becomes important. The general analytical result (a) is confirmed by numerical evaluation of the integrals for a particular form of eddy, while (a), (b), (c) are explained physically by considering the deformation of differently oriented vortex rings. The results are compared with small-scale turbulence approaching bluff bodies where ε[Lt ]1 and β[Lt ] 1.These results also explain dynamically why the intermediate eigenvector of the strain S aligns with the vorticity vector, why the greatest increase in enstrophy production occurs in regions where S has a positive intermediate eigenvalue; and why large-scale strain S of a small-scale vorticity can amplify the small-scale strain rates to a level greater than S – one of the essential characteristics of high-Reynolds-number turbulence.

Vortices in oscillating spin-up

2007 ◽  
Vol 573 ◽  
pp. 339-369 ◽  
Author(s):  
M. G. WELLS ◽  
H. J. H. CLERCX ◽  
G. J. F. VAN HEIJST
Keyword(s):  
Numerical Simulations ◽  
Power Law ◽  
Large Scale ◽  
Laboratory Experiments ◽  
Passive Scalar ◽  
Decay Rates ◽  
Small Scale ◽  
Two Dimensional ◽  
Ekman Pumping ◽  
Viscous Boundary Layer

Laboratory experiments and numerical simulations of oscillating spin-up in a square tank have been conducted to investigate the production of small-scale vorticity near the no-slip sidewalls of the container and the formation and subsequent decay of wall-generated quasi-two-dimensional vortices. The flow is made quasi-two-dimensional by a steady background rotation, and a small sinusoidal perturbation to the background rotation leads to the periodic formation of eddies in the corners of the tank by the roll-up of vorticity generated along the sidewalls. When the oscillation period is greater than the time scale required to advect a full-grown corner vortex to approximately halfway along the sidewall, dipole structures are observed to form. These dipoles migrate away from the walls, and the interior of the tank is continually filled with new vortices. The average size of these vortices appears to be largely controlled by the initial formation mechanism. Their vorticity decays from interactions with other stronger vortices that strip off filaments of vorticity, and by Ekman pumping at the bottom of the tank. Subsequent interactions between the weaker ‘old’ vortices and the ‘young’ vortices result in the straining, and finally the destruction, of older vortices. This inhibits the formation of large-scale vortices with diameters comparable to the size of the container.The laboratory experiments revealed a k−5/3 power law of the energy spectrum for small-to-intermediate wavenumbers. Measurements of the intensity spectrum of a passive scalar were consistent with the Batchelor prediction of a k−1 power law at large wavenumbers. Two-dimensional numerical simulations, under similar conditions to those in the experiments (with weak Ekman decay), were also performed and the simultaneous presence of a k−5/3 and k−3−ζ (with 0 < ζ « 1) power spectrum is observed, with the transition occurring at the wavenumber at which vorticity is injected from the viscous boundary layer into the interior. For higher Ekman decay rates, steeper spectra are obtained for the large wavenumber range, with ζ = O(1) and proportional to the Ekman decay rate. Movies are available with the online version of the paper.

scholarly journals Steady Motion of Ice Sheets

1980 ◽  
Vol 25 (92) ◽  
pp. 229-246 ◽  
Author(s):  
L. W. Morland ◽  
I. R. Johnson
Keyword(s):  
Differential Equation ◽  
Linear Differential Equation ◽  
Power Law ◽  
Perturbation Expansion ◽  
Ice Sheet ◽  
Leading Order ◽  
Plane Flow ◽  
General Law ◽  
Non Linear ◽  
Linear Differential

AbstractSteady plane flow under gravity of a symmetric ice sheet resting on a horizontal rigid bed, subject to surface accumulation and ablation, basal drainage, and basal sliding according to a shear-traction-velocity power law, is treated. The surface accumulation is taken to depend on height, and the drainage and sliding coefficient also depend on the height of overlying ice. The ice is described as a general non-linearly viscous incompressible fluid, with illustrations presented for Glen’s power law, the polynomial law of Colbeck and Evans, and a Newtonian fluid. Uniform temperature is assumed so that effects of a realistic temperature distribution on the ice response are not taken into account. In dimensionless variables a small paramter ν occurs, but the ν = 0 solution corresponds to an unbounded sheet of uniform depth. To obtain a bounded sheet, a horizontal coordinate scaling by a small factor ε(ν) is required, so that the aspect ratio ε of a steady ice sheet is determined by the ice properties, accumulation magnitude, and the magnitude of the central thickness. A perturbation expansion in ε gives simple leading-order terms for the stress and velocity components, and generates a first order non-linear differential equation for the free-surface slope, which is then integrated to determine the profile. The non-linear differential equation can be solved explicitly for a linear sliding law in the Newtonian case. For the general law it is shown that the leading-order approximation is valid both at the margin and in the central zone provided that the power and coefficient in the sliding law satisfy certain restrictions.

Large Eddy Simulations of Staggered Parallel-Plate Flows

2005 ◽  
Author(s):  
M. V. Pham ◽  
F. Plourde ◽  
S. K. Doan
Keyword(s):  
Heat Transfer ◽  
Numerical Simulations ◽  
Heat Transfer Enhancement ◽  
Parallel Plate ◽  
Large Scale ◽  
Primary Objective ◽  
Transfer Enhancement ◽  
Turbulent Structures ◽  
Wide Range ◽  
Large Eddy

Heat transfer enhancement is a subject of major concern in numerous fields of industry and research. Having received undivided attention over the years, it is still studied worldwide. Given the exponential growth of computing power, large-scale numerical simulations are growing steadily more realistic, and it is now possible to obtain accurate time-dependent solutions with far fewer preliminary assumptions about the problems. As a result, an increasingly wide range of physics is now open for exploration. More specifically, it is time to take full advantage of large eddy simulation technique so as to describe heat transfer in staggered parallel-plate flows. In fact, from simple theory through experimental results, it has been demonstrated that surface interruption enhances heat transfer. Staggered parallel-plate geometries are of great potential interest, and yet many numerical works dedicated to them have been tarnished by excessively simple assumptions. That is to say, numerical simulations have generally hypothesized lengthwise periodicity, even though flows are not periodic; moreover, the LES technique has not been employed with sufficient frequency. Actually, our primary objective is to analyze turbulent influence with regard to heat transfers in staggered parallel-plate fin geometries. In order to do so, we have developed a LES code, and numerical results are compared with regard to several grid mesh resolutions. We have focused mainly upon identification of turbulent structures and their role in heat transfer enhancement. Another key point involves the distinct roles of boundary restart and the vortex shedding mechanism on heat transfer and friction factor.

Numerical simulations of Rayleigh–Bénard convection for Prandtl numbers between 10−1 and 104 and Rayleigh numbers between 105 and 109

2010 ◽  
Vol 662 ◽  
pp. 409-446 ◽  
Author(s):  
G. SILANO ◽  
K. R. SREENIVASAN ◽  
R. VERZICCO
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Nusselt Number ◽  
Numerical Simulations ◽  
Multiple Solutions ◽  
Large Scale ◽  
Prandtl Numbers ◽  
Rayleigh Benard Convection ◽  
Rayleigh Numbers ◽  
Rayleigh Benard

We summarize the results of an extensive campaign of direct numerical simulations of Rayleigh–Bénard convection at moderate and high Prandtl numbers (10−1 ≤ Pr ≤ 104) and moderate Rayleigh numbers (105 ≤ Ra ≤ 109). The computational domain is a cylindrical cell of aspect ratio Γ = 1/2, with the no-slip condition imposed on all boundaries. By scaling the numerical results, we find that the free-fall velocity should be multiplied by $1/\sqrt{{\it Pr}}$ in order to obtain a more appropriate representation of the large-scale velocity at high Pr. We investigate the Nusselt and the Reynolds number dependences on Ra and Pr, comparing the outcome with previous numerical and experimental results. Depending on Pr, we obtain different power laws of the Nusselt number with respect to Ra, ranging from Ra2/7 for Pr = 1 up to Ra0.31 for Pr = 103. The Nusselt number is independent of Pr. The Reynolds number scales as ${\it Re}\,{\sim}\,\sqrt{{\it Ra}}/{\it Pr}$, neglecting logarithmic corrections. We analyse the global and local features of viscous and thermal boundary layers and their scaling behaviours with respect to Ra and Pr, and with respect to the Reynolds and Péclet numbers. We find that the flow approaches a saturation state when Reynolds number decreases below the critical value, Res ≃ 40. The thermal-boundary-layer thickness increases slightly (instead of decreasing) when the Péclet number increases, because of the moderating influence of the viscous boundary layer. The simulated ranges of Ra and Pr contain steady, periodic and turbulent solutions. A rough estimate of the transition from the steady to the unsteady state is obtained by monitoring the time evolution of the system until it reaches stationary solutions. We find multiple solutions as long-term phenomena at Ra = 108 and Pr = 103, which, however, do not result in significantly different Nusselt numbers. One of these multiple solutions, even if stable over a long time interval, shows a break in the mid-plane symmetry of the temperature profile. We analyse the flow structures through the transitional phases by direct visualizations of the temperature and velocity fields. A wide variety of large-scale circulation and plume structures has been found. The single-roll circulation is characteristic only of the steady and periodic solutions. For other regimes at lower Pr, the mean flow generally consists of two opposite toroidal structures; at higher Pr, the flow is organized in the form of multi-jet structures, extending mostly in the vertical direction. At high Pr, plumes mainly detach from sheet-like structures. The signatures of different large-scale structures are generally well reflected in the data trends with respect to Ra, less in those with respect to Pr.

Cold welding behavior of metallic glass nanowires: Insights from large-scale numerical simulations

2021 ◽  
Author(s):  
Yuhang Zhang ◽  
Jiejie Li ◽  
Hongjian Zhou ◽  
Yiqun Hu ◽  
Suhang Ding ◽  
...  
Keyword(s):  
Metallic Glass ◽  
Numerical Simulations ◽  
Large Scale ◽  
Cold Welding

scholarly journals Enhancing ecosystem restoration efficiency through spatial and temporal coordination

2015 ◽  
Vol 112 (19) ◽  
pp. 6236-6241 ◽  
Author(s):  
Thomas M. Neeson ◽  
Michael C. Ferris ◽  
Matthew W. Diebel ◽  
Patrick J. Doran ◽  
Jesse R. O’Hanley ◽  
...  
Keyword(s):  
Great Lakes ◽  
Large Scale ◽  
Expert Knowledge ◽  
Spatial Scales ◽  
Small Scale ◽  
Space And Time ◽  
Ecosystem Connectivity ◽  
Breeding Grounds ◽  
Local Expert ◽  
Conservation Projects

In many large ecosystems, conservation projects are selected by a diverse set of actors operating independently at spatial scales ranging from local to international. Although small-scale decision making can leverage local expert knowledge, it also may be an inefficient means of achieving large-scale objectives if piecemeal efforts are poorly coordinated. Here, we assess the value of coordinating efforts in both space and time to maximize the restoration of aquatic ecosystem connectivity. Habitat fragmentation is a leading driver of declining biodiversity and ecosystem services in rivers worldwide, and we simultaneously evaluate optimal barrier removal strategies for 661 tributary rivers of the Laurentian Great Lakes, which are fragmented by at least 6,692 dams and 232,068 road crossings. We find that coordinating barrier removals across the entire basin is nine times more efficient at reconnecting fish to headwater breeding grounds than optimizing independently for each watershed. Similarly, a one-time pulse of restoration investment is up to 10 times more efficient than annual allocations totaling the same amount. Despite widespread emphasis on dams as key barriers in river networks, improving road culvert passability is also essential for efficiently restoring connectivity to the Great Lakes. Our results highlight the dramatic economic and ecological advantages of coordinating efforts in both space and time during restoration of large ecosystems.

scholarly journals Synoptic Control over Orographic Precipitation Distributions during the Olympics Mountains Experiment (OLYMPEX)

2018 ◽  
Vol 146 (4) ◽  
pp. 1023-1044 ◽  
Author(s):  
David J. Purnell ◽  
Daniel J. Kirshbaum
Keyword(s):  
Numerical Simulations ◽  
Physical Interpretation ◽  
Large Scale ◽  
Local Development ◽  
Moisture Flux ◽  
Heavy Precipitation ◽  
Static Stability ◽  
Orographic Precipitation ◽  
Areal Extent ◽  
Precipitation Enhancement

The synoptic controls on orographic precipitation during the Olympics Mountains Experiment (OLYMPEX) are investigated using observations and numerical simulations. Observational precipitation retrievals for six warm-frontal (WF), six warm-sector (WS), and six postfrontal (PF) periods indicate that heavy precipitation occurred in both WF and WS periods, but the latter saw larger orographic enhancements. Such enhancements extended well upstream of the terrain in WF periods but were focused over the windward slopes in both PF and WS periods. Quasi-idealized simulations, constrained by OLYMPEX data, reproduce the key synoptic sensitivities of the OLYMPEX precipitation distributions and thus facilitate physical interpretation. These sensitivities are largely explained by three upstream parameters: the large-scale precipitation rate [Formula: see text], the impinging horizontal moisture flux I, and the low-level static stability. Both WF and WS events exhibit large [Formula: see text] and I, and thus, heavy orographic precipitation, which is greatly enhanced in amplitude and areal extent by the seeder–feeder process. However, the stronger stability of the WF periods, particularly within the frontal inversion (even when it lies above crest level), causes their precipitation enhancement to weaken and shift upstream. In contrast, the small [Formula: see text] and I, larger static stability, and absence of stratiform feeder clouds in the nominally unsaturated and convective PF events yield much lighter time- and area-averaged precipitation. Modest enhancements still occur over the windward slopes due to the local development and invigoration of shallow convective showers.

Sign in / Sign up

Export Citation Format

Share Document