Large eddy simulation of turbine wakes using higher-order methods

NASA Astrophysics Data System (ADS)

Deskos, Georgios; Laizet, Sylvain; Piggott, Matthew D.; Sherwin, Spencer

2017-11-01

Large eddy simulations (LES) of a horizontal-axis turbine wake are presented using the well-known actuator line (AL) model. The fluid flow is resolved by employing higher-order numerical schemes on a 3D Cartesian mesh combined with a 2D Domain Decomposition strategy for an efficient use of supercomputers. In order to simulate flows at relatively high Reynolds numbers for a reasonable computational cost, a novel strategy is used to introduce controlled numerical dissipation to a selected range of small scales. The idea is to mimic the contribution of the unresolved small-scales by imposing a targeted numerical dissipation at small scales when evaluating the viscous term of the Navier-Stokes equations. The numerical technique is shown to behave similarly to the traditional eddy viscosity sub-filter scale models such as the classic or the dynamic Smagorinsky models. The results from the simulations are compared to experimental data for a Reynolds number scaled by the diameter equal to ReD =1,000,000 and both the time-averaged stream wise velocity and turbulent kinetic energy (TKE) are showing a good overall agreement. At the end, suggestions for the amount of numerical dissipation required by our approach are made for the particular case of horizontal-axis turbine wakes.

Impact of Mantle Wind on Subducting Plate Geometry and Interplate Pressure: Insights From Physical Modelling.

NASA Astrophysics Data System (ADS)

Boutelier, D.; Cruden, A. R.

2005-12-01

New physical models of subduction investigate the impact of large-scale mantle flow on the structure of the subducted slab and deformation of the downgoing and overriding plates. The experiments comprise two lithospheric plates made of highly filled silicone polymer resting on a model asthenosphere of low viscosity transparent silicone polymer. Subduction is driven by a piston that pushes the subducting plate at constant rate, a slab-pull force due to the relative density of the slab, and a basal drag force exerted by flow in the model asthenosphere. Large-scale mantle flow is imposed by a second piston moving at constant rate in a tunnel at the bottom of the experiment tank. Passive markers in the mantle track the evolution of flow during the experiment. Slab structure is recorded by side pictures of the experiment while horizontal deformation is studied via passive marker grids on top of both plates. The initial mantle flow direction beneath the overriding plate can be sub-horizontal or sub-vertical. In both cases, as the slab penetrates the mantle, the mantle flow pattern changes to accommodate the subducting high viscosity lithosphere. As the slab continues to descend, the imposed flow produces either over- or under-pressure on the lower surface of the slab depending on the initial mantle flow pattern (sub-horizontal or sub-vertical respectively). Over-pressure imposed on the slab lower surface promotes shallow dip subduction while under-pressure tends to steepen the slab. These effects resemble those observed in previous experiments when the overriding plate moves horizontally with respect to a static asthenosphere. Our experiments also demonstrate that a strong vertical drag force (due to relatively fast downward mantle flow) exerted on the slab results in a decrease in strain rate in both the downgoing and overriding plates, suggesting a decrease in interplate pressure. Furthermore, with an increase in drag force deformation in the downgoing plate can switch from compression to extension. The density contrast between the downgoing plate and asthenosphere is varied from 0% to ~2% in order to investigate the relative contributions of mantle flow and slab pull force on the geometry of the slab and tectonic regime (compressional or extensional).

Internal Gravity Waves Forced by an Isolated Mountain

NASA Astrophysics Data System (ADS)

Nikitina, L.; Campbell, L.

2009-12-01

Density-stratified fluid flow over topography such as mountains, hills and ridges may give rise to internal gravity waves which transport and distribute energy away from their source and have profound effects on the general circulation of the atmosphere and ocean. Much of our knowledge of internal gravity wave dynamics has been acquired from theoretical studies involving mathematical analyses of simplified forms of the governing equations, as well as numerical simulations at varying levels of approximation. In this study, both analytical and numerical methods are used to examine the nonlinear dynamics of gravity waves forced by an isolated mountain. The topography is represented by a lower boundary condition on a two-dimensional rectangular domain and the waves are represented as a perturbation to the background shear flow, thus allowing the use of weakly-nonlinear and multiple-scale asymptotic analyzes. The waves take the form of a packet, localized in the horizontal direction and comprising a continuous spectrum of horizontal wavenumbers centered at zero. For horizontally-localized wave packets, such as those forced by a mountain range with multiple peaks, there are generally two horizontal scales, the fast (short) scale which is defined by the oscillations within the packet and the slow (large) scale which is defined by the horizontal extent of the packet. In the case of an isolated mountain that we examine here, the multiple-scaling procedure is simplified by the absence of a fast spatial scale. The problem is governed by two small parameters that define the height and width of the mountain and approximate solutions are derived in terms of these parameters. Numerical solutions are also carried out to simulate nonlinear critical-level interactions such as the transfer of energy to the background flow by the wave packet, wave reflection and static instability and, eventually, wave breaking leading to turbulence. It is found that for waves forced by an isolated mountain the time frame within which these nonlinear effects become significant depends on both the mountain height and width and that they begin to occur at least an order of magnitude later and the configuration thus remains stable longer than in the case of waves forced by a mountain range of equivalent height.

Numerical studies of dispersion due to tidal flow through Moskstraumen, northern Norway

NASA Astrophysics Data System (ADS)

Lynge, Birgit Kjoss; Berntsen, Jarle; Gjevik, Bjørn

2010-08-01

The effect of horizontal grid resolution on the horizontal relative dispersion of particle pairs has been investigated on a short time scale, i.e. one tidal M 2 cycle. Of particular interest is the tidal effect on dispersion and transports in coastal waters where small-scale flow features are important. A three-dimensional ocean model has been applied to simulate the tidal flow through the Moskstraumen Maelstrom outside Lofoten in northern Norway, well known for its strong current and whirlpools (Gjevik et al., Nature 388(6645):837-838, 1997; Moe et al., Cont Shelf Res 22(3):485-504, 2002). Simulations with spatial resolution down to 50 m have been carried out. Lagrangian tracers were passively advected with the flow, and Lyapunov exponents and power law exponents have been calculated to analyse the separation statistics. It is found that the relative dispersion of particles on a short time scale (12-24 h) is very sensitive to the grid size and that the spatial variability is also very large, ranging from 0 to 100 km2 over a distance of 100 m. This means that models for prediction of transport and dispersion of oil spills, fish eggs, sea lice etc. using a single diffusion coefficient will be of limited value, unless the models actually resolves the small-scale eddies of the tidal current.

Applicability of a diffusion model to lateral transport in the terrestrial and lunar exospheres.

NASA Technical Reports Server (NTRS)

Hodges, R. R., Jr.

1972-01-01

Kinetic theory is used to determine a series expansion of the vertical flux of particles in an exosphere in terms of time and space derivatives of particle concentration, exobase velocity, and temperature. For sufficiently large scale variations of these parameters in time and space, the series can be truncated to a form that is similar to a diffusion equation. Owing to this analogy, it is possible to unite the mathematical description of molecular diffusion, which governs thermospheric flow, and the corresponding exospheric equation by using effective transport coefficients which change smoothly with altitude through the transition from thermosphere to exosphere. A new definition of the exobase for lateral flow emerges from the analogy of exospheric and thermospheric diffusion, as the altitude where the horizontal mean free path length equals the mean horizontal extent of ballistic trajectories of the transported gas, as opposed to the scale height of the dominant gas which determines the exobase for escape. It is shown that the approximation of exospheric lateral flow as a diffusion process is applicable to global scale problems concerning terrestrial helium and heavier gases, and lunar gases heavier than helium.

Momentum flux measurements: Techniques and needs, part 4.5A

NASA Technical Reports Server (NTRS)

Fritts, D. C.

1984-01-01

The vertical flux of horizontal momentum by internal gravity waves is now recognized to play a significant role in the large-scale circulation and thermal structure of the middle atmosphere. This is because a divergence of momentum flux due to wave dissipation results in an acceleration of the local mean flow towards the phase speed of the gravity wave. Such mean flow acceleration are required to offset the large zonal accelerations driven by Coriolis torques acting on the diabatic meridional circulation. Techniques and observations regarding the momentum flux distribution in the middle atmosphere are discussed.

Collection Efficiency and Ice Accretion Characteristics of Two Full Scale and One 1/4 Scale Business Jet Horizontal Tails

NASA Technical Reports Server (NTRS)

Bidwell, Colin S.; Papadakis, Michael

2005-01-01

Collection efficiency and ice accretion calculations have been made for a series of business jet horizontal tail configurations using a three-dimensional panel code, an adaptive grid code, and the NASA Glenn LEWICE3D grid based ice accretion code. The horizontal tail models included two full scale wing tips and a 25 percent scale model. Flow solutions for the horizontal tails were generated using the PMARC panel code. Grids used in the ice accretion calculations were generated using the adaptive grid code ICEGRID. The LEWICE3D grid based ice accretion program was used to calculate impingement efficiency and ice shapes. Ice shapes typifying rime and mixed icing conditions were generated for a 30 minute hold condition. All calculations were performed on an SGI Octane computer. The results have been compared to experimental flow and impingement data. In general, the calculated flow and collection efficiencies compared well with experiment, and the ice shapes appeared representative of the rime and mixed icing conditions for which they were calculated.

Relative dispersion of clustered drifters in a small micro-tidal estuary

NASA Astrophysics Data System (ADS)

Suara, Kabir; Chanson, Hubert; Borgas, Michael; Brown, Richard J.

2017-07-01

Small tide-dominated estuaries are affected by large scale flow structures which combine with the underlying bed generated smaller scale turbulence to significantly increase the magnitude of horizontal diffusivity. Field estimates of horizontal diffusivity and its associated scales are however rare due to limitations in instrumentation. Data from multiple deployments of low and high resolution clusters of GPS-drifters are used to examine the dynamics of a surface flow in a small micro-tidal estuary through relative dispersion analyses. During the field study, cluster diffusivity, which combines both large- and small-scale processes ranged between, 0.01 and 3.01 m2/s for spreading clusters and, -0.06 and -4.2 m2/s for contracting clusters. Pair-particle dispersion, Dp2, was scale dependent and grew as Dp2 ∼ t1.83 in streamwise and Dp2 ∼ t0.8 in cross-stream directions. At small separation scale, pair-particle (d < 0.5 m) relative diffusivity followed the Richardson's 4/3 power law and became weaker as separation scale increases. Pair-particle diffusivity was described as Kp ∼ d1.01 and Kp ∼ d0.85 in the streamwise and cross-stream directions, respectively for separation scales ranging from 0.1 to 10 m. Two methods were used to identify the mechanism responsible for dispersion within the channel. The results clearly revealed the importance of strain fields (stretching and shearing) in the spreading of particles within a small micro-tidal channel. The work provided input for modelling dispersion of passive particle in shallow micro-tidal estuaries where these were not previously experimentally studied.

Size Effect of the 2-D Bodies on the Geothermal Gradient and Q-A Plot

NASA Astrophysics Data System (ADS)

Thakur, M.; Blackwell, D. D.

2009-12-01

Using numerical models we have investigated some of the criticisms on the Q-A plot of related to the effect of size of the body on the slope and reduced heat flow. The effects of horizontal conduction depend on the relative difference of radioactivity between the body and the country rock (assuming constant thermal conductivity). Horizontal heat transfer due to different 2-D bodies was numerically studied in order to quantify resulting temperature differences at the Moho and errors on the predication of Qr (reduced heat flow). Using the two end member distributions of radioactivity, the step model (thickness 10km) and exponential model, different 2-D models of horizontal scale (width) ranging from 10 -500 km were investigated. Increasing the horizontal size of the body tends to move observations closer towards the 1-D solution. A temperature difference of 50 oC is produced (for the step model) at Moho between models of width 10 km versus 500 km. In other words the 1-D solution effectively provides large scale averaging in terms of heat flow and temperature field in the lithosphere. For bodies’ ≤ 100 km wide the geotherms at shallower levels are affected, but at depth they converge and are 50 oC lower than that of the infinite plate model temperature. In case of 2-D bodies surface heat flow is decreased due to horizontal transfer of heat, which will shift the Q-A point vertically downward on the Q-A plot. The smaller the size of the body, the more will be the deviation from the 1-D solution and the more will be the movement of Q-A point downwards on a Q-A plot. On the Q-A plot, a limited points of bodies of different sizes with different radioactivity contrast (for the step and exponential model), exactly reproduce the reduced heat flow Qr. Thus the size of the body can affect the slope on a Q-A plot but Qr is not changed. Therefore, Qr ~ 32 mWm-2 obtained from the global terrain average Q-A plot represents the best estimate of stable continental mantle heat flow.

Revisiting Gill's Circulation. Dynamic Response to Diabatic Heating of Different Horizontal Extents

NASA Astrophysics Data System (ADS)

Reboredo, B.; Bellon, G.

2017-12-01

The horizontal extent of diabatic heating associated with the MJO is thought to be crucial to its development, and the inability of GCMs to simulate the spatial, horizontal organization of clouds is considered a leading hypothesis to explain their limited capacity to simulate MJO events. This prevents the MJO large-circulation response from developing and feeding back on the development of clouds. We apply mid-tropospheric heating of different size in simple linear and non-linear models of the tropical atmosphere following Gill's seminal work on heat-induced tropical circulations. Results show that there is a scale for which the characteristic circulation {Γ c} for the vertical advection of moisture to produce the latent heat mean {Q} gives a rough estimate of the real world MJO scale. Overturning circulation flow rates above {Γ c} account for a circulation that transports more moisture than necessary to be maintained, and below {Γ c}, circulation would not transport enough moisture to maintain circulation. This dynamic scale might constrain the size of the spatially-organised convection necessary to the development of an MJO event. However, other effects are expected to modulate this scale, such as vertical advection of moisture anomalies, horizontal advection, evaporation, radiative heating, and sensible heat fluxes.

Inferring Upper Ocean Dynamics from Horizontal Wavenumber Spectra in the Southern California Current System

NASA Astrophysics Data System (ADS)

Chereskin, T. K.; Gille, S. T.; Rocha, C. B.; Menemenlis, D.

2016-02-01

At the largest horizontal scales (> 100 km), the surface kinetic energy of the ocean appears dominated by a regime of balanced geostrophic motions. At the smallest scales, it transitions to a regime where unbalanced motions (such as internal waves, mixed-layer instabilities, etc.) dominate the surface kinetic energy. The length scale at which the transition occurs depends on the relative energies of balanced and unbalanced motions, which in turn display significant geographic variability. Wavenumber spectra in the upper ocean have been hypothesized to have slopes consistent with either quasi-geostrophic (QG) or surface quasi-geostrophic (SQG) theory. In previous analyses of repeat-track shipboard acoustic Doppler Current profiler (ADCP) velocity observations in the Gulf Stream and the Antarctic Circumpolar Current, spectral slopes were more consistent with QG than SQG theory for length scales between 40 km and 200 km. For scales less than 40 km, the spectra deviated from both QG and SQG theory, and this was attributed in part to internal wave effects. A spectral Helmholtz decomposition was used to split the kinetic energy spectra into rotational and divergent components, identified with balanced and ageostrophic motions, respectively. The California Current System (CCS) provides a contrasting environment characterized by a weak mean flow and an energetic meso- and submeso- scale. It is a nonlinear regime where the amplitude of eddies can be as large as the total steric height increase across the California Current, and hence southward flow in the CCS can, and often is, disrupted by its eddies. This study uses 10 years of shipboard ADCP observations collected on the quarterly cruises of the California Cooperative Oceanic Fisheries Investigations. Horizontal wavenumber spectra from 36 cruises along 6 repeated tracks in the southern CCS that extend from the coast to the subtropical gyre are used to diagnose the dominant governing dynamics at meso- to submeso- scales (10-200 km), with particular attention to the partition into balanced and ageostrophic flows.

Emergence of magnetic flux generated in a solar convective dynamo

NASA Astrophysics Data System (ADS)

Chen, Feng; Rempel, Feng, Matthias; Fan, Yuhong

2016-10-01

We present a realistic numerical model of sunspot and active region formation through the emergence of flux tubes generated in a solar convective dynamo. The magnetic and velocity fields in a horizontal layer near the top boundary of the solar convective dynamo simulation are used as a time-dependent bottom boundary to drive the radiation magnetohydrodynamic simulations of the emergence of the flux tubes through the upper most layer of the convection zone to the photosphere. The emerging flux tubes interact with the convection and break into small scale magnetic elements that further rise to the photosphere. At the photosphere, several bipolar pairs of sunspots are formed through the coalescence of the small scale magnetic elements. The sunspot pairs in the simulation successfully reproduce the fundamental observed properties of solar active regions, including the more coherent leading spots with a stronger field strength, and the correct tilts of the bipolar pairs. These asymmetries originate from the intrinsic asymmetries in the emerging fields imposed at the bottom boundary, where the horizontal fields are already tilted. The leading sides of the emerging flux tubes are up against the downdraft lanes of the giant cells and strongly sheared downward. This leads to the stronger field strength of the leading polarity fields. We find a prograde flow in the emerging flux tube, which is naturally inherited from the solar convective dynamo simulation. The prograde flow gradually becomes a diverging flow as the flux tube rises. The emerging speed is similar to upflow speed of convective motions. The azimuthal average of the flows around a (leading) sunspot reveals a predominant down flow inside the sunspots and a large-scale horizontal inflow at the depth of about 10 Mm. The inflow pattern becomes an outflow in upper most convection zone in the vicinity of the sunspot, which could be considered as moat flows.

Multi-Scale Observation and Modelling of Energy and Matter Exchange in the Atmospheric Boundary-Layer (ScaleX Campaigns)

NASA Astrophysics Data System (ADS)

Zeeman, M. J.; Wolz, K.; Adler, B.; Brenner, C.; De Roo, F.; Emeis, S.; Kalthoff, N.; Mauder, M.; Schäfer, K.; Wohlfahrt, G.; Zhao, P.

2016-12-01

We investigated biosphere-atmosphere exchange processes in relation to the atmospheric boundary-layer (ABL) flow in a shallow valley. Land-use heterogeneity and topography can force local atmospheric flow patterns, including local circulations. Such flow patterns can impair current techniques for the quantification and source attribution of surface-exchange fluxes due to flux-divergence, advection and decoupling. Wind field, temperature and humidity structures in the ABL were observed in high resolution with spatially distributed observations in a 1 km3 experimental domain. Remote-sensing observations of wind, temperature and particles in the ABL (Raman-lidar; RASS; ceilometer; microwave radiometer; 3D Doppler-lidar) were combined with a high-resolution network of in-situ observations that included vertical and horizontal profiles of wind, temperature, carbon dioxide, methane and water vapor concentrations. The experiments were co-located with the long-term eddy covariance (EC) observatory Fendt (DE-Fen; ICOS, TERENO) and were part of international cooperative efforts in 2015 and 2016 (the ScaleX campaigns). The gathered experimental data offers a scale-transcending insight in local flow patterns in mountainous terrain and their influence on surface-exchange fluxes of energy and matter as observed by EC and flux-gradient methodology. In addition, the data is used for validation of Large-Eddy Simulations in complex terrain using PALM-LES. Within this modelling framework, virtual measurements are conducted to further assess the importance of three-dimensional advective and horizontal turbulent transport terms.

Nonlinear anelastic modal theory for solar convection

NASA Technical Reports Server (NTRS)

Latour, J.; Toomre, J.; Zahn, J.-P.

1983-01-01

Solar envelope models are developed using single-mode anelastic equations as a description of turbulent convection which provide estimates for the variation with depth of the largest convective cellular flows, with horizontal sizes comparable to the total depth of the convection zone. These models can be used to describe compressible motions occurring over many density scale heights. Single-mode anelastic solutions are obtained for a solar envelope whose mean stratification is nearly adiabatic over most of its vertical extent because of the enthalpy flux explicitly carried by the big cell, while a subgrid scale representation of turbulent heat transport is incorporated into the treatment near the surface. It is shown that the single-mode equations allow two solutions for the same horizontal wavelength which are distinguished by the sense of the vertical velocity at the center of the three-dimensional cell. It is found that the upward directed flow experiences large pressure effects which can modify the density fluctuations so that the sense of the buoyancy force is changed, with buoyancy braking actually achieved near the top of the convection zone. It is suggested that such dynamical processes may explain why the amplitudes of flows related to the largest scales of convection are so weak in the solar atmosphere.

Characterizing the Severe Turbulence Environments Associated With Commercial Aviation Accidents. Part 1; 44 Case Study Synoptic Observational Analyses

NASA Technical Reports Server (NTRS)

Kaplan, Michael L.; Huffman, Allan W.; Lux, Kevin M.; Charney, Joseph J.; Riordan, Allan J.; Lin, Yuh-Lang; Proctor, Fred H. (Technical Monitor)

2002-01-01

A 44 case study analysis of the large-scale atmospheric structure associated with development of accident-producing aircraft turbulence is described. Categorization is a function of the accident location, altitude, time of year, time of day, and the turbulence category, which classifies disturbances. National Centers for Environmental Prediction Reanalyses data sets and satellite imagery are employed to diagnose synoptic scale predictor fields associated with the large-scale environment preceding severe turbulence. These analyses indicate a predominance of severe accident-producing turbulence within the entrance region of a jet stream at the synoptic scale. Typically, a flow curvature region is just upstream within the jet entrance region, convection is within 100 km of the accident, vertical motion is upward, absolute vorticity is low, vertical wind shear is increasing, and horizontal cold advection is substantial. The most consistent predictor is upstream flow curvature and nearby convection is the second most frequent predictor.

Understanding thermal circulations and near-surface turbulence processes in a small mountain valley

NASA Astrophysics Data System (ADS)

Pardyjak, E.; Dupuy, F.; Durand, P.; Gunawardena, N.; Thierry, H.; Roubin, P.

2017-12-01

The interaction of turbulence and thermal circulations in complex terrain can be significantly different from idealized flat terrain. In particular, near-surface horizontal spatial and temporal variability of winds and thermodynamic variables can be significant event over very small spatial scales. The KASCADE (KAtabatic winds and Stability over CAdarache for Dispersion of Effluents) 2017 conducted from January through March 2017 was designed to address these issues and to ultimately improve prediction of dispersion in complex terrain, particularly during stable atmospheric conditions. We have used a relatively large number of sensors to improve our understanding of the spatial and temporal development, evolution and breakdown of topographically driven flows. KASCADE 2017 consisted of continuous observations and fourteen Intensive Observation Periods (IOPs) conducted in the Cadarache Valley located in southeastern France. The Cadarache Valley is a relatively small valley (5 km x 1 km) with modest slopes and relatively small elevation differences between the valley floor and nearby hilltops ( 100 m). During winter, winds in the valley are light and stably stratified at night leading to thermal circulations as well as complex near-surface atmospheric layering. In this presentation we present results quantifying spatial variability of thermodynamic and turbulence variables as a function of different large -scale forcing conditions (e.g., quiescent conditions, strong westerly flow, and Mistral flow). In addition, we attempt to characterize highly-regular nocturnal horizontal wind meandering and associated turbulence statistics.

Dual-Doppler lidar observation of horizontal convective rolls and near-surface streaks

NASA Astrophysics Data System (ADS)

Iwai, Hironori; Ishii, Shoken; Tsunematsu, Nobumitsu; Mizutani, Kohei; Murayama, Yasuhiro; Itabe, Toshikazu; Yamada, Izumi; Matayoshi, Naoki; Matsushima, Dai; Weiming, Sha; Yamazaki, Takeshi; Iwasaki, Toshiki

2008-07-01

Dual-Doppler lidar and heliborne sensors were used to investigate the three-dimensional (3D) structure of the wind field over Sendai Airport in June 2007. The 3D structures of several-hundred-meter-scale horizontal convective rolls (HCRs) in the sea-breeze layer were observed by the dual-Doppler lidar. The scale of the HCRs determined by the heliborne sensors roughly agreed with that determined by the dual-Doppler lidar. Analysis of the dual-Doppler lidar data showed that the region of upward flow in the HCRs originated in near-surface low-speed streaks. This structure is consistent with the results of large-eddy simulations of the atmospheric boundary layer. The aspect ratios of the HCRs were close to those predicted by linear theories.

Proper horizontal photospheric flows below an eruptive filament

NASA Astrophysics Data System (ADS)

Schmieder, Brigitte; Mein, Pierre; Mein, Nicole; Roudier, Thierry; Chandra, Ramseh

An analysis of the proper motions using SDO/HMI continuum images with the new version of the coherent structure tracking (CST) algorithm developed to track the granules as well as the large scale photospheric flows, was perfomed during three hours in a region containing a large filament channel on September 17, 2010. Supergranules were idenfied in the filament channel. Diverging flows inside the supergranules are similar in and out the filament channel. Using corks, we derived the passive scalar points and produced maps of cork distribution. The anchorage structures with the photosphere (feet) of the filament are located in the areas of converging flows with accumulations of corks. Averaging the velocity vectors for each latitude we defined a profile of the differential rotation. We conclude that the coupling between the convection and magnetic field in the photosphere is relatively strong. The filament experienced the convection motions through its feet. On a large scale point-of-view the differential rotation induced a shear of 0.1 km/s in the filament. On a small scale point-of-view convection motions favored the interaction/cancellation of the parasitic polarities at the base of the feet with the surrounding network explaining the brightenings,/jets and the eruption that were observed in the EUV filament.

Partitioning dynamics of unsaturated flows in fractured porous media: Laboratory studies and three-dimensional multi-scale smoothed particle hydrodynamics simulations of gravity-driven flow in fractures

NASA Astrophysics Data System (ADS)

Kordilla, J.; Bresinsky, L. T.; Shigorina, E.; Noffz, T.; Dentz, M.; Sauter, M.; Tartakovsky, A. M.

2017-12-01

Preferential flow dynamics in unsaturated fractures remain a challenging topic on various scales. On pore- and fracture-scales the highly erratic gravity-driven flow dynamics often provoke a strong deviation from classical volume-effective approaches. Against the common notion that flow in fractures (or macropores) can only occur under equilibrium conditions, i.e., if the surrounding porous matrix is fully saturated and capillary pressures are high enough to allow filling of the fracture void space, arrival times suggest the existence of rapid preferential flow along fractures, fracture networks, and fault zones, even if the matrix is not fully saturated. Modeling such flows requires efficient numerical techniques to cover various flow-relevant physics, such as surface tension, static and dynamic contact angles, free-surface (multi-phase) interface dynamics, and formation of singularities. Here we demonstrate the importance of such flow modes on the partitioning dynamics at simple fracture intersections, with a combination of laboratory experiments, analytical solutions and numerical simulations using our newly developed massively parallel smoothed particle hydrodynamics (SPH) code. Flow modes heavily influence the "bypass" behavior of water flowing along a fracture junction. Flows favoring the formation of droplets exhibit a much stronger bypass capacity compared to rivulet flows, where nearly the whole fluid mass is initially stored within the horizontal fracture. This behavior is demonstrated for a multi-inlet laboratory setup where the inlet-specific flow rate is chosen so that either a droplet or rivulet flow persists. The effect of fluid buffering within the horizontal fracture is presented in terms of dimensionless fracture inflow so that characteristic scaling regimes can be recovered. For both cases (rivulets and droplets), flow within the horizontal fracture transitions into a Washburn regime until a critical threshold is reached and the bypass efficiency increases. For rivulet flows, the initial filling of the horizontal fracture is described by classical plug flow. Meanwhile, for droplet flows, a size-dependent partitioning behavior is observed, and the filling of the fracture takes longer.

Large-scale dynamos in rapidly rotating plane layer convection

NASA Astrophysics Data System (ADS)

Bushby, P. J.; Käpylä, P. J.; Masada, Y.; Brandenburg, A.; Favier, B.; Guervilly, C.; Käpylä, M. J.

2018-05-01

Context. Convectively driven flows play a crucial role in the dynamo processes that are responsible for producing magnetic activity in stars and planets. It is still not fully understood why many astrophysical magnetic fields have a significant large-scale component. Aims: Our aim is to investigate the dynamo properties of compressible convection in a rapidly rotating Cartesian domain, focusing upon a parameter regime in which the underlying hydrodynamic flow is known to be unstable to a large-scale vortex instability. Methods: The governing equations of three-dimensional non-linear magnetohydrodynamics (MHD) are solved numerically. Different numerical schemes are compared and we propose a possible benchmark case for other similar codes. Results: In keeping with previous related studies, we find that convection in this parameter regime can drive a large-scale dynamo. The components of the mean horizontal magnetic field oscillate, leading to a continuous overall rotation of the mean field. Whilst the large-scale vortex instability dominates the early evolution of the system, the large-scale vortex is suppressed by the magnetic field and makes a negligible contribution to the mean electromotive force that is responsible for driving the large-scale dynamo. The cycle period of the dynamo is comparable to the ohmic decay time, with longer cycles for dynamos in convective systems that are closer to onset. In these particular simulations, large-scale dynamo action is found only when vertical magnetic field boundary conditions are adopted at the upper and lower boundaries. Strongly modulated large-scale dynamos are found at higher Rayleigh numbers, with periods of reduced activity (grand minima-like events) occurring during transient phases in which the large-scale vortex temporarily re-establishes itself, before being suppressed again by the magnetic field.

The structure of turbulent flow around vertical plates containing holes and attached to a channel bed

NASA Astrophysics Data System (ADS)

Basnet, K.; Constantinescu, G.

2017-11-01

High-resolution, 3-D large eddy simulations are conducted to study the physics of flow past 2-D solid and porous vertical plates of height H mounted on a horizontal surface (no bottom gap) with a fully developed, turbulent incoming flow. The porous plate consists of an array of spanwise-oriented, identical solid cylinders of rectangular cross section. The height of the solid cylinders and the spacing between the solid cylinders, corresponding to the plate's "holes," are kept constant for any given configuration, as the present study considers only plates of uniform porosity. The paper discusses how the mean flow and turbulence structure around the vertical plate, the unsteady forces acting on the plate, the dynamics of the large-scale turbulent eddies, the spectral content of the wake, and the distribution of the bed friction velocity on the horizontal channel bed vary as a function of the plate porosity (0% < P < 36%), the relative spacing between the solid elements of the porous plate (d/H), and the roughness of the channel bed surface. Simulation results are used to explain how the bleeding flow affects the dynamics on the larger billow eddies advected in the separated shear layer (SSL) forming at the top of the plate and the wake structure. It is found that the main recirculation eddy in the wake remains attached to the plate for P < 30%. For larger porosities, the main recirculation eddy forms away from the porous plate. The energy of the billows advected in the SSL decays monotonically with increasing plate porosity. For cases when the recirculation eddy remains attached to the plate, the larger billows advected in the downstream part of the SSL are partially reinjected inside the main recirculation eddy as a result of their interaction with the channel bed. This creates a feedback mechanism that induces large-scale disturbances of the spanwise-oriented vortex tubes advected inside the upstream part of the SSL. Results also show that the mean drag coefficient and the root-mean-square of the drag coefficient fluctuations increase mildly with increasing d/H. Meanwhile, varying d/H has a negligible effect on the position and size of the main recirculation eddy. The presence of large-scale roughness elements (2-D ribs) at the bed results in the decrease of the mean drag coefficient of the plate and, in the case of a solid plate, in a large decrease of the frequency of the large-scale eddies advected in the SSL.

Tidal bending of ice shelves as a mechanism for large-scale temporal variations in ice flow

NASA Astrophysics Data System (ADS)

Rosier, Sebastian H. R.; Hilmar Gudmundsson, G.

2018-05-01

GPS measurements reveal strong modulation of horizontal ice shelf and ice stream flow at a variety of tidal frequencies, most notably a fortnightly (Msf) frequency not present in the vertical tides themselves. Current theories largely fail to explain the strength and prevalence of this signal over floating ice shelves. We show how well-known non-linear aspects of ice rheology can give rise to widespread, long-periodic tidal modulation in ice shelf flow, generated within ice shelves themselves through tidal flexure acting at diurnal and semidiurnal frequencies. Using full-Stokes viscoelastic modelling, we show that inclusion of tidal bending within the model accounts for much of the observed tidal modulation of ice shelf flow. Furthermore, our model shows that, in the absence of vertical tidal forcing, the mean flow of the ice shelf is reduced by almost 30 % for the geometry that we consider.

Detection and reconstruction of large scale flow structures in a river by means of empirical mode decomposition combined with Hilbert transform

NASA Astrophysics Data System (ADS)

Franca, Mário J.; Lemmin, Ulrich

2014-05-01

The occurrence of large scale flow structures (LSFS) coherently organized throughout the flow depth has been reported in field and laboratory experiments of flows over gravel beds, especially under low relative submergence conditions. In these, the instantaneous velocity is synchronized over the whole vertical profile oscillating at a low frequency above or below the time-averaged value. The detection of large scale coherently organized regions in the flow field is often difficult since it requires detailed simultaneous observations of the flow velocities at several levels. The present research avoids the detection problem by using an Acoustic Doppler Velocity Profiler (ADVP), which permits measuring three-dimensional velocities quasi-simultaneously over the full water column. Empirical mode decomposition (EMD) combined with the application of the Hilbert transform is then applied to the instantaneous velocity data to detect and isolate LSFS. The present research was carried out in a Swiss river with low relative submergence of 2.9, herein defined as h/D50, (where h is the mean flow depth and D50 the bed grain size diameter for which 50% of the grains have smaller diameters). 3D ADVP instantaneous velocity measurements were made on a 3x5 rectangular horizontal grid (x-y). Fifteen velocity profiles were equally spaced in the spanwise direction with a distance of 10 cm, and in the streamwise direction with a distance of 15 cm. The vertical resolution of the measurements is roughly 0.5 cm. A measuring grid covering a 3D control volume was defined. The instantaneous velocity profiles were measured for 3.5 min with a sampling frequency of 26 Hz. Oscillating LSFS are detected and isolated in the instantaneous velocity signal of the 15 measured profiles. Their 3D cycle geometry is reconstructed and investigated through phase averaging based on the identification of the instantaneous signal phase (related to the Hilbert transform) applied to the original raw signal. Results for all the profiles are consistent and indicate clearly the presence of LSFS throughout the flow depth with impact on the three components of the velocity profile and on the bed friction velocity. A high correlation of the movement is found throughout the flow depth, thus corroborating the hypothesis of large-scale coherent motion evolving over the whole water depth. These latter are characterized in terms of period, horizontal scale and geometry. The high spatial and temporal resolution of our ADVP was crucial for obtaining comprehensive results on coherent structures dynamics. EMD combined with the Hilbert transform have previously been successfully applied to geophysical flow studies. Here we show that this method can also be used for the analysis of river dynamics. In particular, we demonstrate that a clean, well-behaved intrinsic mode function can be obtained from a noisy velocity time series that allowed a precise determination of the vertical structure of the coherent structures. The phase unwrapping of the UMR and the identification of the phase related velocity components brings new insight into the flow dynamics Research supported by the Swiss National Science Foundation (2000-063818). KEY WORDS: large scale flow structures (LSFS); gravel-bed rivers; empirical mode decomposition; Hilbert transform

Three-Dimensional Structure of the Circulation Induced by a Shoaling Topographic Wave

NASA Astrophysics Data System (ADS)

Mizuta, G.; Hogg, N. G.

2003-12-01

Rectification of Rossby wave energy has been proposed as a mechanism for the maintenance of the recirculation cell of the Gulf Stream (Hogg 1988; Rizzoli et al. 1995). We investigated the three-dimensional structure of potential-vorticity flux and a mean flow induced by a topographic wave incident over a bottom slope analytically and numerically, focusing on the limit that bottom friction is the dominant dissipation process. In this limit it is shown that the topographic wave cannot be a steady source of the potential vorticity outside the bottom Ekman layer. Instead, the distribution of potential vorticity is determined from the initial transient of the topographic wave. This potential vorticity and the heat flux by the topographic wave at the bottom determine the mean flow, and give a relation between the horizontal and vertical scales of the mean flow. When the horizontal scale of the mean flow is larger than the internal deformation radius, the mean flow is almost constant with depth independent of whether or not the topographic wave is trapped near the bottom. Then the mean flow at the bottom is proportional to the divergence of vertically integrated Reynolds stress ∫ -D0 /line{u'v'} dz. This divergence, which is caused by bottom friction, is large when the group velocity, cg and the vertical scale, μ -1 of the wave motion are small. Thus the mean flow tends to be large where cg and μ -1 become small, and decreases as the topographic wave is dissipated by bottom friction. Since bottom friction also dissipates the mean flow, the mean flow asymptotes to a constant value as the friction becomes zero. These features of the potential-vorticity flux and the mean flow are reproduced in numerical experiments. It is also shown from the numerical experiment that the distribution of the mean flow depends on the amplitude of the wave because of the Doppler shift of the wave by the mean flow. These feature of the mean flow are preserved when we used stratification and bottom topography resembling to those over the continental slope near the Gulf Stream. The transport of the mean flow is about 20 Sv when the wave amplitude is about 2 cm/s. These numbers are similiar to those observed in the Gulf Stream region.

Effects of Time-Dependent Inflow Perturbations on Turbulent Flow in a Street Canyon

NASA Astrophysics Data System (ADS)

Duan, G.; Ngan, K.

2017-12-01

Urban flow and turbulence are driven by atmospheric flows with larger horizontal scales. Since building-resolving computational fluid dynamics models typically employ steady Dirichlet boundary conditions or forcing, the accuracy of numerical simulations may be limited by the neglect of perturbations. We investigate the sensitivity of flow within a unit-aspect-ratio street canyon to time-dependent perturbations near the inflow boundary. Using large-eddy simulation, time-periodic perturbations to the streamwise velocity component are incorporated via the nudging technique. Spatial averages of pointwise differences between unperturbed and perturbed velocity fields (i.e., the error kinetic energy) show a clear dependence on the perturbation period, though spatial structures are largely insensitive to the time-dependent forcing. The response of the error kinetic energy is maximized for perturbation periods comparable to the time scale of the mean canyon circulation. Frequency spectra indicate that this behaviour arises from a resonance between the inflow forcing and the mean motion around closed streamlines. The robustness of the results is confirmed using perturbations derived from measurements of roof-level wind speed.

Tests of dynamic Lagrangian eddy viscosity models in Large Eddy Simulations of flow over three-dimensional bluff bodies

NASA Astrophysics Data System (ADS)

Tseng, Yu-Heng; Meneveau, Charles; Parlange, Marc B.

2004-11-01

Large Eddy Simulations (LES) of atmospheric boundary-layer air movement in urban environments are especially challenging due to complex ground topography. Typically in such applications, fairly coarse grids must be used where the subgrid-scale (SGS) model is expected to play a crucial role. A LES code using pseudo-spectral discretization in horizontal planes and second-order differencing in the vertical is implemented in conjunction with the immersed boundary method to incorporate complex ground topography, with the classic equilibrium log-law boundary condition in the new-wall region, and with several versions of the eddy-viscosity model: (1) the constant-coefficient Smagorinsky model, (2) the dynamic, scale-invariant Lagrangian model, and (3) the dynamic, scale-dependent Lagrangian model. Other planar-averaged type dynamic models are not suitable because spatial averaging is not possible without directions of statistical homogeneity. These SGS models are tested in LES of flow around a square cylinder and of flow over surface-mounted cubes. Effects on the mean flow are documented and found not to be major. Dynamic Lagrangian models give a physically more realistic SGS viscosity field, and in general, the scale-dependent Lagrangian model produces larger Smagorinsky coefficient than the scale-invariant one, leading to reduced distributions of resolved rms velocities especially in the boundary layers near the bluff bodies.

Large eddy simulation of incompressible turbulent channel flow

NASA Technical Reports Server (NTRS)

Moin, P.; Reynolds, W. C.; Ferziger, J. H.

1978-01-01

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Large-Eddy Simulations of Atmospheric Flows Over Complex Terrain Using the Immersed-Boundary Method in the Weather Research and Forecasting Model

NASA Astrophysics Data System (ADS)

Ma, Yulong; Liu, Heping

2017-12-01

Atmospheric flow over complex terrain, particularly recirculation flows, greatly influences wind-turbine siting, forest-fire behaviour, and trace-gas and pollutant dispersion. However, there is a large uncertainty in the simulation of flow over complex topography, which is attributable to the type of turbulence model, the subgrid-scale (SGS) turbulence parametrization, terrain-following coordinates, and numerical errors in finite-difference methods. Here, we upgrade the large-eddy simulation module within the Weather Research and Forecasting model by incorporating the immersed-boundary method into the module to improve simulations of the flow and recirculation over complex terrain. Simulations over the Bolund Hill indicate improved mean absolute speed-up errors with respect to previous studies, as well an improved simulation of the recirculation zone behind the escarpment of the hill. With regard to the SGS parametrization, the Lagrangian-averaged scale-dependent Smagorinsky model performs better than the classic Smagorinsky model in reproducing both velocity and turbulent kinetic energy. A finer grid resolution also improves the strength of the recirculation in flow simulations, with a higher horizontal grid resolution improving simulations just behind the escarpment, and a higher vertical grid resolution improving results on the lee side of the hill. Our modelling approach has broad applications for the simulation of atmospheric flows over complex topography.

Granular slumping on a horizontal surface

NASA Astrophysics Data System (ADS)

Lajeunesse, E.; Monnier, J. B.; Homsy, G. M.

2005-10-01

We report the results of an experimental investigation of the flow induced by the collapse of a column of granular material (glass beads of diameter d) over a horizontal surface. Two different setups are used, namely, a rectangular channel and a semicircular tube, allowing us to compare two-dimensional and axisymmetric flows, with particular focus on the internal flow structure. In both geometries the flow dynamics and the deposit morphologies are observed to depend primarily on the initial aspect ratio of the granular column a =Hi/Li, where Hi is the height of the initial granular column and Li its length along the flow direction. Two distinct regimes are observed depending on a: an avalanche of the column flanks producing truncated deposits for small a and a column free fall leading to conical deposits for large a. In both geometries the characteristic time scale is the free fall of the granular column τc=√Hi/g . The flow initiated by Coulomb-like failure never involves the whole granular heap but remains localized in a surface layer whose size and shape depend on a and vary in both space and time. Except in the vicinity of the pile foot where the flow is pluglike, velocity profiles measured at the side wall are identical to those commonly observed in steady granular surface flows: the velocity varies linearly with depth in the flowing layer and decreases exponentially with depth in the static layer. Moreover, the shear rate is constant, γ˙=0.3√g /d , independent of the initial aspect ratio, the flow geometry, position along the heap, or time. Despite the rather complex flow dynamics, the scaled deposit height Hf/Li and runout distance ΔL /Li both exhibit simple power laws whose exponents depend on a and on the flow geometry. We show that the physical origin of these power laws can be understood on the basis of a dynamic balance between acceleration, pressure gradient, and friction forces at the foot of the granular pile. Two asymptotic behaviors can be distinguished: the flow is dominated by friction forces at small a and by pressure forces at large a. The effect of the flow geometry is determined primarily by mass conservation and becomes important only for large a.

Characterization of meter-scale spatial variability of riverbed hydraulic conductivity in a lowland river (Aa River, Belgium)

NASA Astrophysics Data System (ADS)

Ghysels, Gert; Benoit, Sien; Awol, Henock; Jensen, Evan Patrick; Debele Tolche, Abebe; Anibas, Christian; Huysmans, Marijke

2018-04-01

An improved general understanding of riverbed heterogeneity is of importance for all groundwater modeling studies that include river-aquifer interaction processes. Riverbed hydraulic conductivity (K) is one of the main factors controlling river-aquifer exchange fluxes. However, the meter-scale spatial variability of riverbed K has not been adequately mapped as of yet. This study aims to fill this void by combining an extensive field measurement campaign focusing on both horizontal and vertical riverbed K with a detailed geostatistical analysis of the meter-scale spatial variability of riverbed K . In total, 220 slug tests and 45 standpipe tests were performed at two test sites along the Belgian Aa River. Omnidirectional and directional variograms (along and across the river) were calculated. Both horizontal and vertical riverbed K vary over several orders of magnitude and show significant meter-scale spatial variation. Horizontal K shows a bimodal distribution. Elongated zones of high horizontal K along the river course are observed at both sections, indicating a link between riverbed structures, depositional environment and flow regime. Vertical K is lognormally distributed and its spatial variability is mainly governed by the presence and thickness of a low permeable organic layer at the top of the riverbed. The absence of this layer in the center of the river leads to high vertical K and is related to scouring of the riverbed by high discharge events. Variograms of both horizontal and vertical K show a clear directional anisotropy with ranges along the river being twice as large as those across the river.

Agitation, Mixing, and Transfers Induced by Bubbles

NASA Astrophysics Data System (ADS)

Risso, Frédéric

2018-01-01

Bubbly flows involve bubbles randomly distributed within a liquid. At large Reynolds number, they experience an agitation that can combine shear-induced turbulence (SIT), large-scale buoyancy-driven flows, and bubble-induced agitation (BIA). The properties of BIA strongly differ from those of SIT. They have been determined from studies of homogeneous swarms of rising bubbles. Regarding the bubbles, agitation is mainly caused by the wake-induced path instability. Regarding the liquid, two contributions must be distinguished. The first one corresponds to the anisotropic flow disturbances generated near the bubbles, principally in the vertical direction. The second one is the almost isotropic turbulence induced by the flow instability through a population of bubbles, which turns out to be the main cause of horizontal fluctuations. Both contributions generate a k-3 spectral subrange and exponential probability density functions. The subsequent issue will be to understand how BIA interacts with SIT.

Migration of the Cratering Flow-Field Center with Implications for Scaling Oblique Impacts

NASA Technical Reports Server (NTRS)

Anderson, J. L. B.; Schultz, P. H.; Heineck, J. T.

2004-01-01

Crater-scaling relationships are used to predict many cratering phenomena such as final crater diameter and ejection speeds. Such nondimensional relationships are commonly determined from experimental impact and explosion data. Almost without exception, these crater-scaling relationships have used data from vertical impacts (90 deg. to the horizontal). The majority of impact craters, however, form by impacts at angles near 45 deg. to the horizontal. While even low impact angles result in relatively circular craters in sand targets, the effects of impact angle have been shown to extend well into the excavation stage of crater growth. Thus, the scaling of oblique impacts needs to be investigated more thoroughly in order to quantify fully how impact angle affects ejection speed and angle. In this study, ejection parameters from vertical (90 deg.) and 30 deg. oblique impacts are measured using three-dimensional particle image velocimetry (3D PIV) at the NASA Ames Vertical Gun Range (AVGR). The primary goal is to determine the horizontal migration of the cratering flow-field center (FFC). The location of the FFC at the time of ejection controls the scaling of oblique impacts. For vertical impacts the FFC coincides with the impact point (IP) and the crater center (CC). Oblique impacts reflect a more complex, horizontally migrating flow-field. A single, stationary point-source model cannot be used accurately to describe the evolution of the ejection angles from oblique impacts. The ejection speeds for oblique impacts also do not follow standard scaling relationships. The migration of the FFC needs to be understood and incorporated into any revised scaling relationships.

Effect of Unsaturated Flow Modes on Partitioning Dynamics of Gravity-Driven Flow at a Simple Fracture Intersection: Laboratory Study and Three-Dimensional Smoothed Particle Hydrodynamics Simulations

NASA Astrophysics Data System (ADS)

Kordilla, Jannes; Noffz, Torsten; Dentz, Marco; Geyer, Tobias; Tartakovsky, Alexandre M.

2017-11-01

In this work, we study gravity-driven flow of water in the presence of air on a synthetic surface intersected by a horizontal fracture and investigate the importance of droplet and rivulet flow modes on the partitioning behavior at the fracture intersection. We present laboratory experiments, three-dimensional smoothed particle hydrodynamics (SPH) simulations using a heavily parallelized code, and a theoretical analysis. The flow-rate-dependent mode switching from droplets to rivulets is observed in experiments and reproduced by the SPH model, and the transition ranges agree in SPH simulations and laboratory experiments. We show that flow modes heavily influence the "bypass" behavior of water flowing along a fracture junction. Flows favoring the formation of droplets exhibit a much stronger bypass capacity compared to rivulet flows, where nearly the whole fluid mass is initially stored within the horizontal fracture. The effect of fluid buffering within the horizontal fracture is presented in terms of dimensionless fracture inflow so that characteristic scaling regimes can be recovered. For both cases (rivulets and droplets), the flow within the horizontal fracture transitions into a Washburn regime until a critical threshold is reached and the bypass efficiency increases. For rivulet flows, the initial filling of the horizontal fracture is described by classical plug flow. Meanwhile, for droplet flows, a size-dependent partitioning behavior is observed, and the filling of the fracture takes longer. For the case of rivulet flow, we provide an analytical solution that demonstrates the existence of classical Washburn flow within the horizontal fracture.

Flow Environment Study Near the Empennage of a 15-Percent Scale Helicopter Model

NASA Technical Reports Server (NTRS)

Gorton, Susan Althoff; Berry, John D.; Hodges, W. Todd; Reis, Deane G.

2000-01-01

Development of advanced rotorcraft configurations has highlighted a need for high-quality experimental data to support the development of flexible and accurate analytical design tools. To provide this type of data, a test program was conducted in the Langley 14- by 22-Foot Subsonic Tunnel to measure the flow near the empennage of a 15-percent scale powered helicopter model with an operating tail fan. Three-component velocity profiles were measured with laser velocimetry (LV) one chord forward of the horizontal tail for four advance ratios to evaluate the effect of the rotor wake impingement on the horizontal tail angle of attack. These velocity data indicate the horizontal tail can experience unsteady angle of attack variations of over 30 degrees due to the rotor wake influence. The horizontal tail is most affected by the rotor wake above advance ratios of 0.10. Velocity measurements of the flow on the inlet side of the tail fan were made for a low-speed flight condition using conventional LV techniques. The velocity data show an accelerated flow near the tail fan duct, and vorticity calculations track the passage of main rotor wake vortices through the measurement plane.

Numerical Modeling and Experimental Analysis of Scale Horizontal Axis Marine Hydrokinetic (MHK) Turbines

NASA Astrophysics Data System (ADS)

Javaherchi, Teymour; Stelzenmuller, Nick; Seydel, Joseph; Aliseda, Alberto

2013-11-01

We investigate, through a combination of scale model experiments and numerical simulations, the evolution of the flow field around the rotor and in the wake of Marine Hydrokinetic (MHK) turbines. Understanding the dynamics of this flow field is the key to optimizing the energy conversion of single devices and the arrangement of turbines in commercially viable arrays. This work presents a comparison between numerical and experimental results from two different case studies of scaled horizontal axis MHK turbines (45:1 scale). In the first case study, we investigate the effect of Reynolds number (Re = 40,000 to 100,000) and Tip Speed Ratio (TSR = 5 to 12) variation on the performance and wake structure of a single turbine. In the second case, we study the effect of the turbine downstream spacing (5d to 14d) on the performance and wake development in a coaxial configuration of two turbines. These results provide insights into the dynamics of Horizontal Axis Hydrokinetic Turbines, and by extension to Horizontal Axis Wind Turbines in close proximity to each other, and highlight the capabilities and limitations of the numerical models. Once validated at laboratory scale, the numerical model can be used to address other aspects of MHK turbines at full scale. Supported by DOE through the National Northwest Marine Renewable Energy Center.

Stratified turbulence diagnostics for high-Reynolds-number momentum wakes

NASA Astrophysics Data System (ADS)

Diamessis, Peter; Zhou, Qi

2017-11-01

We analyze a large-eddy simulation (LES) dataset of the turbulent wake behind a sphere of diameter D translating at speed U in a linearly stratified Boussinesq fluid with buoyancy frequency N. These simulations are performed at Reynolds numbers Re ≡ UD / ν ∈ { 5 ×103 , 105 , 4 ×105 } and various Froude numbers Fr ≡ 2 U /(ND) . The recently obtained data at Re = 4 ×105 , the highest Re attained so far in either simulation or laboratory, and Fr ∈ { 4 , 16 } enable us to systematically investigate the effects of Reynolds number on this prototypical localized stratified turbulent shear flow. Our analysis focuses on the time evolution of various diagnostics of stratified turbulence, such as the horizontal and vertical integral length scales, turbulent kinetic energy and its dissipation rate ɛ, and the local rate of shear between the spontaneously formed layers of vorticity within the larger-scale quasi-horizontal flow structures. This leads to a discussion of the transitions between distinct stratified flow regimes (Brethouwer et al. 2007) in the appropriately defined phase diagram, and we highlight the dynamical role of the Gibson number Gi = ɛ /(νN2) , and its dependence on the body-based Reynolds number Re . ONR Grants N00014-13-1-0665 and N00014-15-1-2513.

Measuring mixing efficiency in experiments of strongly stratified turbulence

NASA Astrophysics Data System (ADS)

Augier, P.; Campagne, A.; Valran, T.; Calpe Linares, M.; Mohanan, A. V.; Micard, D.; Viboud, S.; Segalini, A.; Mordant, N.; Sommeria, J.; Lindborg, E.

2017-12-01

Oceanic and atmospheric models need better parameterization of the mixing efficiency. Therefore, we need to measure this quantity for flows representative of geophysical flows, both in terms of types of flows (with vortices and/or waves) and of dynamical regimes. In order to reach sufficiently large Reynolds number for strongly stratified flows, experiments for which salt is used to produce the stratification have to be carried out in a large rotating platform of at least 10-meter diameter.We present new experiments done in summer 2017 to study experimentally strongly stratified turbulence and mixing efficiency in the Coriolis platform. The flow is forced by a slow periodic movement of an array of large vertical or horizontal cylinders. The velocity field is measured by 3D-2C scanned horizontal particles image velocimetry (PIV) and 2D vertical PIV. Six density-temperature probes are used to measure vertical and horizontal profiles and signals at fixed positions.We will show how we rely heavily on open-science methods for this study. Our new results on the mixing efficiency will be presented and discussed in terms of mixing parameterization.

Kinematic dynamo action in square and hexagonal patterns.

PubMed

Favier, B; Proctor, M R E

2013-11-01

We consider kinematic dynamo action in rapidly rotating Boussinesq convection just above onset. The velocity is constrained to have either a square or a hexagonal pattern. For the square pattern, large-scale dynamo action is observed at onset, with most of the magnetic energy being contained in the horizontally averaged component. As the magnetic Reynolds number increases, small-scale dynamo action becomes possible, reducing the overall growth rate of the dynamo. For the hexagonal pattern, the breaking of symmetry between up and down flows results in an effective pumping velocity. For intermediate rotation rates, this additional effect can prevent the growth of any mean-field dynamo, so that only a small-scale dynamo is eventually possible at large enough magnetic Reynolds number. For very large rotation rates, this pumping term becomes negligible, and the dynamo properties of square and hexagonal patterns are qualitatively similar. These results hold for both perfectly conducting and infinite magnetic permeability boundary conditions.

Analysis of Doppler Lidar Data Acquired During the Pentagon Shield Field Campaign

DOE Office of Scientific and Technical Information (OSTI.GOV)

Newsom, Rob K.

2011-04-14

Observations from two coherent Doppler lidars deployed during the Pentagon Shield field campaign are analyzed in conjunction with other sensors to characterize the overall boundary-layer structure, and identify the dominant flow characteristics during the entire two-week field campaign. Convective boundary layer (CBL) heights and cloud base heights (CBH) are estimated from an analysis of the lidar signal-to-noise-ratio (SNR), and mean wind profiles are computed using a modified velocity-azimuth-display (VAD) algorithm. Three-dimensional wind field retrievals are computed from coordinated overlapping volume scans, and the results are analyzed by visualizing the flow in horizontal and vertical cross sections. The VAD winds showmore » that southerly flows dominate during the two-week field campaign. Low-level jets (LLJ) were evident on all but two of the nights during the field campaign. The LLJs tended to form a couple hours after sunset and reach maximum strength between 03 and 07 UTC. The surface friction velocities show distinct local maxima during four nights when strong LLJs formed. Estimates of the convective boundary layer height and residual layer height are obtained through an analysis of the vertical gradient of the lidar signal-to-noise-ratio (SNR). Strong minimum in the SNR gradient often develops just above the surface after sunrise. This minimum is associated with the developing CBL, and increases rapidly during the early portion of the daytime period. On several days, this minimum continues to increase until about sunset. Secondary minima in the SNR gradient were also observed at higher altitudes, and are believed to be remnants of the CBL height from previous days, i.e. the residual layer height. The dual-Doppler analysis technique used in this study makes use of hourly averaged radial velocity data to produce three-dimensional grids of the horizontal velocity components, and the horizontal velocity variance. Visualization of horizontal and vertical cross sections of the dual-Doppler wind retrievals often indicated a jet-like flow feature over the Potomac River under southerly flow conditions. This linear flow feature is roughly aligned with the Potomac River corridor to the south of the confluence with the Anatostia River, and is most apparent at low levels (i.e. below ~150 m MSL). It is believed that this flow arises due to reduced drag over the water surface and when the large scale flow aligns with the Potomac River corridor. A so-called area-constrained VAD analysis generally confirmed the observations from the dual-Doppler analysis. When the large scale flow is southerly, wind speeds over the Potomac River are consistently larger than the at a site just to the west of the river for altitudes less than 100 m MSL. Above this level, the trend is somewhat less obvious. The data suggest that the depth of the wind speed maximum may be reduced by strong directional shear aloft.« less

Deep Convection, Magnetism and Solar Supergranulation

NASA Astrophysics Data System (ADS)

Lord, J. W.

We examine the effect of deep convection and magnetic fields on solar supergranulation. While supergranulation was originally identified as a convective flow from relatively great depth below the solar surface, recent work suggests that supergranules may originate near the surface. We use the MURaM code to simulate solar-like surface convection with a realistic photosphere and domain size up to 197 x 197 x 49 Mm3. This yields nearly five orders of magnitude of density contrast between the bottom of the domain and the photosphere which is the most stratified solar-like convection simulations that we are aware of. Magnetic fields were thought to be a passive tracer in the photosphere, but recent work suggests that magnetism could provide a mechanism that enhances the supergranular scale flows at the surface. In particular, the enhanced radiative losses through long lived magnetic network elements may increase the lifetime of photospheric downflows and help organize low wavenumber flows. Since our simulation does not have sufficient resolution to resolve increased cooling by magnetic bright points, we artificially increase the radiative cooling in elements with strong magnetic flux. These simulations increase the cooling by 10% for magnetic field strength greater than 100 G. We find no statistically significant difference in the velocity or magnetic field spectrum by enhancing the radiative cooling. We also find no differences in the time scale of the flows or the length scales of the magnetic energy spectrum. This suggests that the magnetic field is determined by the flows and is largely a passive tracer. We use these simulations to construct a two-component model of the flows: for scales smaller than the driving (integral) scale (which is four times the local density scale height) the flows follow a Kolmogorov (k-5/3) spectrum, while larger scale modes decay with height from their driving depth (i.e. the depth where the wavelength of the mode is equal to the driving (integral) scale). This model reproduces the MURaM results well and suggests that the low wavenumber power in the photosphere imprints from below. In particular, the amplitude of the driving (integral) scale mode at each depth determines how much power imprints on the surface flows. This is validated by MURaM simulations of varying depth that show that increasing depths contribute power at a particular scale (or range of scales) that is always at lower wavenumbers than shallower flows. The mechanism for this imprinting remains unclear but, given the importance of the balances in the continuity equation to determining the spectrum of the flows, we suggest that pressure perturbations in the convective upflows are the imprinting mechanism. By comparing the MURaM simulations to SDO/HMI observations (using the coherent structure tracking code to compute the inferred horizontal velocities on both data sets), we find that the simulations have significant excess power for scales larger than supergranulation. The only way to match observations is by using an artificial energy flux to transport the solar luminosity for all depths greater than 10 Mm below the photosphere (down to the bottom of the domain at 49 Mm depth). While magnetic fields from small-scale dynamo simulations help reduce the rms velocity required to transport the solar luminosity below the surface, this provides only a small reduction in low wavenumber power in the photosphere. The convective energy transport in the Sun is constrained by theoretical models and the solar radiative luminosity. The amplitude or scale of the convective flows that transport the energy, however, are not constrained. The strong low wavenumber flows found in these local simulations are also present in current generation global simulations. While local or global dynamo magnetic fields may help suppress these large-scale flows, the magnetic fields must be substantially stronger throughout the convection domains for these simulations to match observations. The significant decrease in low wavenumber flow amplitude in the artificial energy flux simulation that matches the observed photospheric horizontal velocity spectrum suggests that convection in the Sun transports the solar luminosity with much weaker large-scale flows. This suggests that we do not understand how convective transport works in the Sun for depths greater than 10 Mm below the photosphere.

Controls on groundwater flow in the Bengal Basin of India and Bangladesh: Regional modeling analysis

USGS Publications Warehouse

Michael, H.A.; Voss, C.I.

2009-01-01

Groundwater for domestic and irrigation purposes is produced primarily from shallow parts of the Bengal Basin aquifer system (India and Bangladesh), which contains high concentrations of dissolved arsenic (exceeding worldwide drinking water standards), though deeper groundwater is generally low in arsenic. An essential first step for determining sustainable management of the deep groundwater resource is identification of hydrogeologic controls on flow and quantification of basin-scale groundwater flow patterns. Results from groundwater modeling, in which the Bengal Basin aquifer system is represented as a single aquifer with higher horizontal than vertical hydraulic conductivity, indicate that this anisotropy is the primary hydrogeologic control on the natural flowpath lengths. Despite extremely low hydraulic gradients due to minimal topographic relief, anisotropy implies large-scale (tens to hundreds of kilometers) flow at depth. Other hydrogeologic factors, including lateral and vertical changes in hydraulic conductivity, have minor effects on overall flow patterns. However, because natural hydraulic gradients are low, the impact of pumping on groundwater flow is overwhelming; modeling indicates that pumping has substantially changed the shallow groundwater budget and flowpaths from predevelopment conditions. ?? Springer-Verlag 2009.

Controls on groundwater flow in the Bengal Basin of India and Bangladesh: regional modeling analysis

NASA Astrophysics Data System (ADS)

Michael, Holly A.; Voss, Clifford I.

2009-11-01

Groundwater for domestic and irrigation purposes is produced primarily from shallow parts of the Bengal Basin aquifer system (India and Bangladesh), which contains high concentrations of dissolved arsenic (exceeding worldwide drinking water standards), though deeper groundwater is generally low in arsenic. An essential first step for determining sustainable management of the deep groundwater resource is identification of hydrogeologic controls on flow and quantification of basin-scale groundwater flow patterns. Results from groundwater modeling, in which the Bengal Basin aquifer system is represented as a single aquifer with higher horizontal than vertical hydraulic conductivity, indicate that this anisotropy is the primary hydrogeologic control on the natural flowpath lengths. Despite extremely low hydraulic gradients due to minimal topographic relief, anisotropy implies large-scale (tens to hundreds of kilometers) flow at depth. Other hydrogeologic factors, including lateral and vertical changes in hydraulic conductivity, have minor effects on overall flow patterns. However, because natural hydraulic gradients are low, the impact of pumping on groundwater flow is overwhelming; modeling indicates that pumping has substantially changed the shallow groundwater budget and flowpaths from predevelopment conditions.

Reorientations of the large-scale flow in turbulent convection in a cube

NASA Astrophysics Data System (ADS)

Foroozani, N.; Niemela, J. J.; Armenio, V.; Sreenivasan, K. R.

2017-03-01

Large-eddy simulations of turbulent Rayleigh-Bénard convection were conducted for a fluid of Prandtl number Pr=0.7 confined in a cube, for Rayleigh numbers of 106 and 108. The model solves the unsteady Navier-Stokes equations under the Boussinesq approximation, using a dynamic Smagorinsky model with a Lagrangian averaging technique for the subgrid terms. Under fully developed conditions the flow topology is characterized by a large-scale circulation (LSC) developing in a plane containing one of the diagonals of the cell, while two counter-rotating vortices consequently develop in the other diagonal plane, resulting in a strong inflow at the horizontal midplane. This flow structure is not static, with the LSC undergoing nonperiodic reorientations, or switching, between the two diagonal planes; hence, we supplement the observations of the three-dimensional time-averaged flow structures with single point measurements (time series) to shed light on the dynamics of the reorientations. For all observations, this switching results from a lateral rotation of the LSC in which some finite time spent in a transient state where the large-scale circulation is parallel to one set of side walls; there are, importantly, no observations consistent with so-called cessations of the LSC, in which it decays and then reforms in another plane without such a rotation. The average switching rate for the LSC is in excellent agreement with the results of Bai et al. [K. Bai, D. Ji, and E. Brown, Phys. Rev. E 93, 023117 (2016), 10.1103/PhysRevE.93.023117].

Reorientations of the large-scale flow in turbulent convection in a cube.

PubMed

Foroozani, N; Niemela, J J; Armenio, V; Sreenivasan, K R

2017-03-01

Large-eddy simulations of turbulent Rayleigh-Bénard convection were conducted for a fluid of Prandtl number Pr=0.7 confined in a cube, for Rayleigh numbers of 10^{6} and 10^{8}. The model solves the unsteady Navier-Stokes equations under the Boussinesq approximation, using a dynamic Smagorinsky model with a Lagrangian averaging technique for the subgrid terms. Under fully developed conditions the flow topology is characterized by a large-scale circulation (LSC) developing in a plane containing one of the diagonals of the cell, while two counter-rotating vortices consequently develop in the other diagonal plane, resulting in a strong inflow at the horizontal midplane. This flow structure is not static, with the LSC undergoing nonperiodic reorientations, or switching, between the two diagonal planes; hence, we supplement the observations of the three-dimensional time-averaged flow structures with single point measurements (time series) to shed light on the dynamics of the reorientations. For all observations, this switching results from a lateral rotation of the LSC in which some finite time spent in a transient state where the large-scale circulation is parallel to one set of side walls; there are, importantly, no observations consistent with so-called cessations of the LSC, in which it decays and then reforms in another plane without such a rotation. The average switching rate for the LSC is in excellent agreement with the results of Bai et al. [K. Bai, D. Ji, and E. Brown, Phys. Rev. E 93, 023117 (2016)PLEEE81539-375510.1103/PhysRevE.93.023117].

Sheared bioconvection in a horizontal tube

NASA Astrophysics Data System (ADS)

Croze, O. A.; Ashraf, E. E.; Bees, M. A.

2010-12-01

The recent interest in using microorganisms for biofuels is motivation enough to study bioconvection and cell dispersion in tubes subject to imposed flow. To optimize light and nutrient uptake, many microorganisms swim in directions biased by environmental cues (e.g. phototaxis in algae and chemotaxis in bacteria). Such taxes inevitably lead to accumulations of cells, which, as many microorganisms have a density different to the fluid, can induce hydrodynamic instabilites. The large-scale fluid flow and spectacular patterns that arise are termed bioconvection. However, the extent to which bioconvection is affected or suppressed by an imposed fluid flow and how bioconvection influences the mean flow profile and cell transport are open questions. This experimental study is the first to address these issues by quantifying the patterns due to suspensions of the gravitactic and gyrotactic green biflagellate alga Chlamydomonas in horizontal tubes subject to an imposed flow. With no flow, the dependence of the dominant pattern wavelength at pattern onset on cell concentration is established for three different tube diameters. For small imposed flows, the vertical plumes of cells are observed merely to bow in the direction of flow. For sufficiently high flow rates, the plumes progressively fragment into piecewise linear diagonal plumes, unexpectedly inclined at constant angles and translating at fixed speeds. The pattern wavelength generally grows with flow rate, with transitions at critical rates that depend on concentration. Even at high imposed flow rates, bioconvection is not wholly suppressed and perturbs the flow field.

Thermal convection of liquid metal in a long inclined cylinder

NASA Astrophysics Data System (ADS)

Teimurazov, Andrei; Frick, Peter

2017-11-01

The turbulent convection of low-Prandtl-number fluids (Pr=0.0083 ) in a long cylindrical cell, heated at one end face and cooled at the other, inclined to the vertical at angle β , 0 ≤β ≤π /2 with step π /20 , is studied numerically by solving the Oberbeck-Boussinesq equations with the large-eddy-simulation approach for small-scale turbulence. The cylinder length is L =5 D , where D is the diameter. The Rayleigh number, determined by the cylinder diameter, is of the order of 5 ×106 . We show that the structure of the flow strongly depends on the inclination angle. A stable large-scale circulation (LSC) slightly disturbed by small-scale turbulence exists in the horizontal cylinder. The deviation from a horizontal position provides strong amplification of both LSC and small-scale turbulence. The energy of turbulent pulsations increases monotonically with decreasing inclination angle β , matching the energy of the LSC at β ≈π /5 . The intensity of the LSC has a wide, almost flat, maximum for an inclined cylinder and slumps approaching the vertical position, in which the LSC vanishes. The dependence of the Nusselt number on the inclination angle has a maximum at β ≈7 π /20 and generally follows the dependence of the intensity of LSC on the inclination. This indicates that the total heat transport is highly determined by LSC. We examine the applicability of idealized thermal boundary conditions (BCs) for modeling a real experiment with liquid sodium flows. Therefore, the simulations are done with two types of temperature BCs: fixed face temperature and fixed heat flux. The intensity of the LSC is slightly higher in the latter case and leads to a corresponding increase of the Nusselt number and enhancement of temperature pulsations.

A Composite Medium Approximation for Moisture Tension-Dependent Anisotropy in Unsaturated Layered Sediments

NASA Astrophysics Data System (ADS)

Pruess, K.

2001-12-01

Sedimentary formations often have a layered structure in which hydrogeologic properties have substantially larger correlation length in the bedding plane than perpendicular to it. Laboratory and field experiments and observations have shown that even small-scale layering, down to millimeter-size laminations, can substantially alter and impede the downward migration of infiltrating liquids, while enhancing lateral flow. The fundamental mechanism is that of a capillary barrier: at increasingly negative moisture tension (capillary suction pressure), coarse-grained layers with large pores desaturate more quickly than finer-grained media. This strongly reduces the hydraulic conductivity of the coarser (higher saturated hydraulic conductivity) layers, which then act as barriers to downward flow, forcing water to accumulate and spread near the bottom of the overlying finer-grained material. We present a "composite medium approximation" (COMA) for anisotropic flow behavior on a typical grid block scale (0.1 - 1 m or larger) in finite-difference models. On this scale the medium is conceptualized as consisting of homogeneous horizontal layers with uniform thickness, and capillary equilibrium is assumed to prevail locally. Directionally-dependent relative permeabilities are obtained by considering horizontal flow to proceed via "conductors in parallel," while vertical flow involves "resistors in series." The model is formulated for the general case of N layers, and implementation of a simplified two-layer (fine-coarse) approximation in the multiphase flow simulator TOUGH2 is described. The accuracy of COMA is evaluated by comparing numerical simulations of plume migration in 1-D and 2-D unsaturated flow with results of fine-grid simulations in which all layers are discretized explicitly. Applications to water seepage and solute transport at the Hanford site are also described. This work was supported by the U.S. Department of Energy under Contract No. DE-AC03-76SF00098 through Memorandum Purchase Order 248861-A-B2 between Pacific Northwest National Laboratory and Lawrence Berkeley National Laboratory.

Large-scale dynamo action precedes turbulence in shearing box simulations of the magnetorotational instability

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bhat, Pallavi; Ebrahimi, Fatima; Blackman, Eric G.

Here, we study the dynamo generation (exponential growth) of large-scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal (x–y) averaging, we also demonstrate the presence of large-scale fields when vertical (y–z) averaging is employed instead. By computing space–time planar averaged fields and power spectra, we find large-scale dynamo action in the early MRI growth phase – a previously unidentified feature. Non-axisymmetric linear MRI modes with low horizontal wavenumbers and vertical wavenumbers near that of expected maximal growth, amplify the large-scale fields exponentially before turbulence and high wavenumbermore » fluctuations arise. Thus the large-scale dynamo requires only linear fluctuations but not non-linear turbulence (as defined by mode–mode coupling). Vertical averaging also allows for monitoring the evolution of the large-scale vertical field and we find that a feedback from horizontal low wavenumber MRI modes provides a clue as to why the large-scale vertical field sustains against turbulent diffusion in the non-linear saturation regime. We compute the terms in the mean field equations to identify the individual contributions to large-scale field growth for both types of averaging. The large-scale fields obtained from vertical averaging are found to compare well with global simulations and quasi-linear analytical analysis from a previous study by Ebrahimi & Blackman. We discuss the potential implications of these new results for understanding the large-scale MRI dynamo saturation and turbulence.« less

Large-scale dynamo action precedes turbulence in shearing box simulations of the magnetorotational instability

DOE PAGES

Bhat, Pallavi; Ebrahimi, Fatima; Blackman, Eric G.

2016-07-06

Here, we study the dynamo generation (exponential growth) of large-scale (planar averaged) fields in unstratified shearing box simulations of the magnetorotational instability (MRI). In contrast to previous studies restricted to horizontal (x–y) averaging, we also demonstrate the presence of large-scale fields when vertical (y–z) averaging is employed instead. By computing space–time planar averaged fields and power spectra, we find large-scale dynamo action in the early MRI growth phase – a previously unidentified feature. Non-axisymmetric linear MRI modes with low horizontal wavenumbers and vertical wavenumbers near that of expected maximal growth, amplify the large-scale fields exponentially before turbulence and high wavenumbermore » fluctuations arise. Thus the large-scale dynamo requires only linear fluctuations but not non-linear turbulence (as defined by mode–mode coupling). Vertical averaging also allows for monitoring the evolution of the large-scale vertical field and we find that a feedback from horizontal low wavenumber MRI modes provides a clue as to why the large-scale vertical field sustains against turbulent diffusion in the non-linear saturation regime. We compute the terms in the mean field equations to identify the individual contributions to large-scale field growth for both types of averaging. The large-scale fields obtained from vertical averaging are found to compare well with global simulations and quasi-linear analytical analysis from a previous study by Ebrahimi & Blackman. We discuss the potential implications of these new results for understanding the large-scale MRI dynamo saturation and turbulence.« less

Ground-based LiDAR Measurements of Actively Inflating Pahoehoe Flows, Kilauea Volcano, Hawaii: Implications for Emplacement of Basaltic Units on Mars

NASA Astrophysics Data System (ADS)

Byrnes, J. M.; Finnegan, D. C.; Nicoll, K.; Anderson, S. W.

2007-05-01

Remote sensing datasets enable planetary volcanologists to extract information regarding eruption processes. Long-lived effusive eruptions at sites such as Kilauea Volcano (HI) provide opportunities to collect rich observational data sets, including detailed measurements of topography and extrusion rates, that allow comparisons between lava flow surface morphologies and emplacement conditions for use in interpreting similar morphological features associated with planetary lava flows. On Mars, the emplacement of basaltic lava flows is a volumetrically and spatially important process, creating both large-scale and small-scale surface morphologies. On Earth, low effusion rate eruptions on relatively horizontal slopes tend to create inflated lava flows that display hummocky topography. To better understand the processes involved in creating observed surface characteristics, we repeatedly measured the surface topography of an actively flowing and inflating basaltic unit within the Pu'u O'o flow field over a 5-day period. We used a ground-based laser-scanner (LiDAR) system that provided vertical and horizontal accuracies of 4 mm. Comparing DEMs from repeated laser scans yielded the magnitudes and styles of constructional processes, allowing us to quantify the relationship between pre- and post-emplacement surface topography. Our study site (roughly 200 m x 200 m) experienced about 5 m of vertical inflation over a 3 day period and created a new hummocky surface containing several tumuli. The temporal and spatial patterns of inflation were complex and showed no obvious relationship with underlying topography. High-precision morphometric measurements acquired using ground-based LiDAR affords us the opportunity to capture the essential boundary conditions necessary for evaluating and comparing high-resolution planetary data sets, such as those acquired by the MOC, HRSC, and HiRISE instruments.

Large perceptual distortions of locomotor action space occur in ground-based coordinates: Angular expansion and the large-scale horizontal-vertical illusion.

PubMed

Klein, Brennan J; Li, Zhi; Durgin, Frank H

2016-04-01

What is the natural reference frame for seeing large-scale spatial scenes in locomotor action space? Prior studies indicate an asymmetric angular expansion in perceived direction in large-scale environments: Angular elevation relative to the horizon is perceptually exaggerated by a factor of 1.5, whereas azimuthal direction is exaggerated by a factor of about 1.25. Here participants made angular and spatial judgments when upright or on their sides to dissociate egocentric from allocentric reference frames. In Experiment 1, it was found that body orientation did not affect the magnitude of the up-down exaggeration of direction, suggesting that the relevant orientation reference frame for this directional bias is allocentric rather than egocentric. In Experiment 2, the comparison of large-scale horizontal and vertical extents was somewhat affected by viewer orientation, but only to the extent necessitated by the classic (5%) horizontal-vertical illusion (HVI) that is known to be retinotopic. Large-scale vertical extents continued to appear much larger than horizontal ground extents when observers lay sideways. When the visual world was reoriented in Experiment 3, the bias remained tied to the ground-based allocentric reference frame. The allocentric HVI is quantitatively consistent with differential angular exaggerations previously measured for elevation and azimuth in locomotor space. (PsycINFO Database Record (c) 2016 APA, all rights reserved).

Large perceptual distortions of locomotor action space occur in ground-based coordinates: Angular expansion and the large-scale horizontal-vertical illusion

PubMed Central

Klein, Brennan J.; Li, Zhi; Durgin, Frank H.

2015-01-01

What is the natural reference frame for seeing large-scale spatial scenes in locomotor action space? Prior studies indicate an asymmetric angular expansion in perceived direction in large-scale environments: Angular elevation relative to the horizon is perceptually exaggerated by a factor of 1.5, whereas azimuthal direction is exaggerated by a factor of about 1.25. Here participants made angular and spatial judgments when upright or on their sides in order to dissociate egocentric from allocentric reference frames. In Experiment 1 it was found that body orientation did not affect the magnitude of the up-down exaggeration of direction, suggesting that the relevant orientation reference frame for this directional bias is allocentric rather than egocentric. In Experiment 2, the comparison of large-scale horizontal and vertical extents was somewhat affected by viewer orientation, but only to the extent necessitated by the classic (5%) horizontal-vertical illusion (HVI) that is known to be retinotopic. Large-scale vertical extents continued to appear much larger than horizontal ground extents when observers lay sideways. When the visual world was reoriented in Experiment 3, the bias remained tied to the ground-based allocentric reference frame. The allocentric HVI is quantitatively consistent with differential angular exaggerations previously measured for elevation and azimuth in locomotor space. PMID:26594884

Magnetic flux concentrations from turbulent stratified convection

NASA Astrophysics Data System (ADS)

Käpylä, P. J.; Brandenburg, A.; Kleeorin, N.; Käpylä, M. J.; Rogachevskii, I.

2016-04-01

Context. The formation of magnetic flux concentrations within the solar convection zone leading to sunspot formation is unexplained. Aims: We study the self-organization of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains representing the uppermost 8.5-24 Mm of the solar convection zone with the horizontal size of the domain varying between 34 and 96 Mm. The density contrast in the 24 Mm deep models is more than 3 × 103 or eight density scale heights, corresponding to a little over 12 pressure scale heights. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow in set-ups where no small-scale dynamos are present. In the most highly stratified cases we employ the reduced sound speed method to relax the time step constraint arising from the high sound speed in the deep layers. We model radiation via the diffusion approximation and neglect detailed radiative transfer in order to concentrate on purely magnetohydrodynamic effects. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface in cases with moderate and high density stratification, corresponding to domain depths of 12.5 and 24 Mm. The size of the concentrations increases as the box size increases and the largest structures (20 Mm horizontally near the surface) are obtained in the models that are 24 Mm deep. The field strength in the concentrations is in the range of 3-5 kG, almost independent of the magnitude of the imposed field. The amplitude of the concentrations grows approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in most of the domain, which is unfavourable for the operation of the negative effective magnetic pressure instability (NEMPI). Simulations in which a passive vector field is evolved do not show a noticeable difference from magnetohydrodynamic runs in terms of the growth of the structures. Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is a tangling of the large-scale field rather than an instability. One plausible mechanism that can explain both the linear growth and the concentration of the flux in the regions of converging flow pattern is flux expulsion. A possible reason for the absence of NEMPI is that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign. Furthermore, there may not be sufficient scale separation, which is required for NEMPI to work. Movies associated to Figs. 4 and 5 are available in electronic form at http://www.aanda.org

Subsurface Supergranular Vertical Flows as Measured Using Large Distance Separations in Time-Distance Helioseismology

NASA Technical Reports Server (NTRS)

Duvall, T. L., Jr.; Hanasoge, S. M.

2012-01-01

As large-distance rays (say, 10-24 deg) approach the solar surface approximately vertically, travel times measured from surface pairs for these large separations are mostly sensitive to vertical flows, at least for shallow flows within a few Mm of the solar surface. All previous analyses of supergranulation have used smaller separations and have been hampered by the difficulty of separating the horizontal and vertical flow components. We find that the large separation travel times associated with upergranulation cannot be studied using the standard phase-speed filters of time-distance helioseismology. These filters, whose use is based upon a refractive model of the perturbations,reduce the resultant travel time signal by at least an order of magnitude at some distances. More effective filters are derived. Modeling suggests that the center-annulus travel time difference in the separation range 10-24 deg is insensitive to the horizontally diverging flow from the centers of the supergranules and should lead to a constant signal from the vertical flow. Our measurement of this quantity for the average supergranule, 5.1 s, is constant over the distance range. This magnitude of signal cannot be caused by the level of upflow at cell centers seen at the photosphere of 10 m/s extended in depth. It requires the vertical flow to increase with depth. A simple Gaussian model of the increase with depth implies a peak upward flow of 240 m/s at a depth of 2.3 Mm and a peak horizontal flow of 700 m/s at a depth of 1.6 Mm.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kitada, Toshihiro; Okamura, Kiyoshi; Tanaka, Setsu

Japanese large cities, mostly located in coastal region, have rapidly expanded in the last three decades. People in the region now experience warmer and unpleasant thermal environment in summer season, supposedly because of the extensive urbanization. Especially, under fine weather with light synoptic-scale gradient wind, the highest temperature often appears in rather inland area. An explanation to this phenomenon is horizontal transport of heat by the sea breeze which is strongly heated during its passage over the coastal urban area. This study investigated how and how much urbanization in Nohbi Plain of central Japan influences local wind and temperature distributionmore » over the plain and its surrounding region, Nohbi Plain which faces to the Pacific Ocean and is surrounded by the Japanese Alps, and thus is in a typical topographic situation in Japan. For the study we performed numerical simulations using our meso-scale meteorological model which includes k-c model for turbulence. One of the major results obtained is to have clarified a hierarchy in natural and artificial topography of various scales on their contributions to formation of characteristic diurnal pattern of wind and temperature distributions in the plain area. The Japanese Alps, the largest topographic feature in central Japan and often called the roof of Japan, gave the most important influence on the wind, although the mountains are located quite far, around 200 km, from Nohbi Plain. The way for the influence of the high mountains was to warm air mass in upper layer but below 2 km in altitude, over the Nohbi Plain through heat transport due to a large scale circulation consisting of the {open_quotes}flows from plain to plateau in lower layer and plateau to plain in upper layer{close_quotes}. Urbanization in the coastal Nohbi Plain showed little effect on the flow pattern, but caused large temperature rises at surface level in the inland area because of horizontal heat transfer by the sea breeze.« less

Comparison of Large eddy dynamo simulation using dynamic sub-grid scale (SGS) model with a fully resolved direct simulation in a rotating spherical shell

NASA Astrophysics Data System (ADS)

Matsui, H.; Buffett, B. A.

2017-12-01

The flow in the Earth's outer core is expected to have vast length scale from the geometry of the outer core to the thickness of the boundary layer. Because of the limitation of the spatial resolution in the numerical simulations, sub-grid scale (SGS) modeling is required to model the effects of the unresolved field on the large-scale fields. We model the effects of sub-grid scale flow and magnetic field using a dynamic scale similarity model. Four terms are introduced for the momentum flux, heat flux, Lorentz force and magnetic induction. The model was previously used in the convection-driven dynamo in a rotating plane layer and spherical shell using the Finite Element Methods. In the present study, we perform large eddy simulations (LES) using the dynamic scale similarity model. The scale similarity model is implement in Calypso, which is a numerical dynamo model using spherical harmonics expansion. To obtain the SGS terms, the spatial filtering in the horizontal directions is done by taking the convolution of a Gaussian filter expressed in terms of a spherical harmonic expansion, following Jekeli (1981). A Gaussian field is also applied in the radial direction. To verify the present model, we perform a fully resolved direct numerical simulation (DNS) with the truncation of the spherical harmonics L = 255 as a reference. And, we perform unresolved DNS and LES with SGS model on coarser resolution (L= 127, 84, and 63) using the same control parameter as the resolved DNS. We will discuss the verification results by comparison among these simulations and role of small scale fields to large scale fields through the role of the SGS terms in LES.

Material transport in a wind and buoyancy forced mixed layer

NASA Astrophysics Data System (ADS)

Mensa, J. A.; Özgökmen, T.; Poje, A. C.; Imberger, J.

2016-02-01

Flows in the upper ocean mixed layer are responsible for the transport and dispersion of biogeochemical tracers, phytoplankton and buoyant pollutants, such as hydrocarbons from an oil spill. Material dispersion in mixed layer flows subject to diurnal buoyancy forcing and weak winds (|u10|=5ms-1) are investigated using a non-hydrostatic model. Both purely buoyancy-forced and combined wind- and buoyancy-forced flows are sampled using passive tracers, as well as 2D and 3D particles to explore characteristics of horizontal and vertical dispersion. It is found that the surface tracer patterns are determined by the convergence zones created by convection cells within a time scale of just a few hours. For pure convection, the results displayed the classic signature of Rayleigh-Benard cells. When combined with a wind stress, the convective cells become anisotropic in that the along-wind length scale gets much larger than the cross-wind scale. Horizontal relative dispersion computed by sampling the flow fields using both 2D and 3D passive particles is found to be consistent with the Richardson regime. Relative dispersion is an order of magnitude higher and 2D surface releases transition to Richardson regime faster in the wind-forced case. We also show that the buoyancy-forced case results in significantly lower amplitudes of scale-dependent horizontal relative diffusivity, kD(l), than those reported by Okubo (1970), while the wind- and buoyancy forced case shows a good agreement with Okubo's diffusivity amplitude, and scaling consistent with Richardson's 4/3rd law, kD(l) l4/3. The modelling results provide a framework for measuring material dispersion by mixed layer flow in future observational programs.

Material transport in a convective surface mixed layer under weak wind forcing

NASA Astrophysics Data System (ADS)

Mensa, Jean A.; Özgökmen, Tamay M.; Poje, Andrew C.; Imberger, Jörg

2015-12-01

Flows in the upper ocean mixed layer are responsible for the transport and dispersion of biogeochemical tracers, phytoplankton and buoyant pollutants, such as hydrocarbons from an oil spill. Material dispersion in mixed layer flows subject to diurnal buoyancy forcing and weak winds (| u10 | = 5m s-1) are investigated using a non-hydrostatic model. Both purely buoyancy-forced and combined wind- and buoyancy-forced flows are sampled using passive tracers, as well as 2D and 3D particles to explore characteristics of horizontal and vertical dispersion. It is found that the surface tracer patterns are determined by the convergence zones created by convection cells within a time scale of just a few hours. For pure convection, the results displayed the classic signature of Rayleigh-Benard cells. When combined with a wind stress, the convective cells become anisotropic in that the along-wind length scale gets much larger than the cross-wind scale. Horizontal relative dispersion computed by sampling the flow fields using both 2D and 3D passive particles is found to be consistent with the Richardson regime. Relative dispersion is an order of magnitude higher and 2D surface releases transition to Richardson regime faster in the wind-forced case. We also show that the buoyancy-forced case results in significantly lower amplitudes of scale-dependent horizontal relative diffusivity, kD(ℓ), than those reported by Okubo (1970), while the wind- and buoyancy-forced case shows a good agreement with Okubo's diffusivity amplitude, and the scaling is consistent with Richardson's 4/3rd law, kD ∼ ℓ4/3. These modeling results provide a framework for measuring material dispersion by mixed layer flows in future observational programs.

Dispersive stresses in wind farms

NASA Astrophysics Data System (ADS)

Segalini, Antonio; Braunbehrens, Robert; Hyvarinen, Ann

2017-11-01

One of the most famous models of wind farms is provided by the assumption that the farm can be approximated as a horizontally-homogeneous forest canopy with vertically-varying force intensity. By means of this approximation, the flow-motion equations become drastically simpler, as many of the three-dimensional effects are gone. However, the application of the horizontal average operator to the RANS equations leads to the appearance of new transport terms (called dispersive stresses) originating from the horizontal (small-scale) variation of the mean velocity field. Since these terms are related to the individual turbine signature, they are expected to vanish outside the roughness sublayer, providing a definition for the latter. In the present work, an assessment of the dispersive stresses is performed by means of a wake-model approach and through the linearised code ORFEUS developed at KTH. Both approaches are very fast and enable the characterization of a large number of wind-farm layouts. The dispersive stress tensor and its effect on the turbulence closure models are investigated, providing guidelines for those simulations where it is impossible to resolve the farm at a turbine scale due to grid requirements (as, for instance, mesoscale simulations).

Kinetic energy budgets in areas of convection

NASA Technical Reports Server (NTRS)

Fuelberg, H. E.

1979-01-01

Synoptic scale budgets of kinetic energy are computed using 3 and 6 h data from three of NASA's Atmospheric Variability Experiments (AVE's). Numerous areas of intense convection occurred during the three experiments. Large kinetic energy variability, with periods as short as 6 h, is observed in budgets computed over each entire experiment area and over limited volumes that barely enclose the convection and move with it. Kinetic energy generation and transport processes in the smaller volumes are often a maximum when the enclosed storms are near peak intensity, but the nature of the various energy processes differs between storm cases and seems closely related to the synoptic conditions. A commonly observed energy budget for peak storm intensity indicates that generation of kinetic energy by cross-contour flow is the major energy source while dissipation to subgrid scales is the major sink. Synoptic scale vertical motion transports kinetic energy from lower to upper levels of the atmosphere while low-level horizontal flux convergence and upper-level horizontal divergence also occur. Spatial fields of the energy budget terms show that the storm environment is a major center of energy activity for the entire area.

Pseudo-incompressible, finite-amplitude gravity waves: wave trains and stability

NASA Astrophysics Data System (ADS)

Schlutow, Mark; Klein, Rupert

2017-04-01

Based on weak asymptotic WKB-like solutions for two-dimensional atmospheric gravity waves (GWs) traveling wave solutions (wave trains) are derived and analyzed with respect to stability. A systematic multiple-scale analysis using the ratio of the dominant wavelength and the scale height as a scale separation parameter is applied on the fully compressible Euler equations. A distinguished limit favorable for GWs close to static instability, reveals that pseudo-incompressible rather than Boussinesq theory applies. A spectral expansion including a mean flow, combined with the additional WKB assumption of slowly varying phases and amplitudes, is used to find general weak asymptotic solutions. This ansatz allows for arbitrarily strong, non-uniform stratification and holds even for finite-amplitude waves. It is deduced that wave trains as leading order solutions can only exist if either some non-uniform background stratification is given but the wave train propagates only horizontally or if the wave train velocity vector is given but the background is isothermal. For the first case, general analytical solutions are obtained that may be used to model mountain lee waves. For the second case with the additional assumption of horizontal periodicity, upward propagating wave train fronts were found. These wave train fronts modify the mean flow beyond the non-acceleration theorem. Stability analysis reveal that they are intrinsically modulationally unstable. The range of validity for the scale separation parameter was tested with fully nonlinear simulations. Even for large values an excellent agreement with the theory was found.

Limitations of one-dimensional mesoscale PBL parameterizations in reproducing mountain-wave flows

DOE PAGES

Munoz-Esparza, Domingo; Sauer, Jeremy A.; Linn, Rodman R.; ...

2015-12-08

In this study, mesoscale models are considered to be the state of the art in modeling mountain-wave flows. Herein, we investigate the role and accuracy of planetary boundary layer (PBL) parameterizations in handling the interaction between large-scale mountain waves and the atmospheric boundary layer. To that end, we use recent large-eddy simulation (LES) results of mountain waves over a symmetric two-dimensional bell-shaped hill [Sauer et al., J. Atmos. Sci. (2015)], and compare them to four commonly used PBL schemes. We find that one-dimensional PBL parameterizations produce reasonable agreement with the LES results in terms of vertical wavelength, amplitude of velocitymore » and turbulent kinetic energy distribution in the downhill shooting flow region. However, the assumption of horizontal homogeneity in PBL parameterizations does not hold in the context of these complex flow configurations. This inappropriate modeling assumption results in a vertical wavelength shift producing errors of ≈ 10 m s–1 at downstream locations due to the presence of a coherent trapped lee wave that does not mix with the atmospheric boundary layer. In contrast, horizontally-integrated momentum flux derived from these PBL schemes displays a realistic pattern. Therefore results from mesoscale models using ensembles of one-dimensional PBL schemes can still potentially be used to parameterize drag effects in general circulation models. Nonetheless, three-dimensional PBL schemes must be developed in order for mesoscale models to accurately represent complex-terrain and other types of flows where one-dimensional PBL assumptions are violated.« less

Quiescent Prominence Structure and Dynamics: a new View From the Hinode/SOT

NASA Astrophysics Data System (ADS)

Berger, T.; Okamoto, J.; Slater, G.; Magara, T.; Tarbell, T.; Tsuneta, S.; Hurlburt, N.

2008-05-01

To date the Hinode/Solar Optical Telescope (SOT) has produced over a dozen sub-arcsecond, multi-hour movies of quiescent solar prominences in both the Ca II 396.8~nm H-line and the H-alpha 656.3~nm line. These datasets have revealed new details of the structure and dynamics of quiescent prominences including a new form of mass transport in the form of buoyant plume upflows from the chromosphere. We review the SOT prominence datasets to show that quiescent prominences appear in two major morphological categories: "vertically" and "horizontally" structured. The vertically structured prominences all show ubiquitous downflows in 400--700~km wide "streams" with velocities of approximately 10~km~s-1. Most of the vertically structured prominences also show episodic upflows in the form of dark turbulent plumes with typical velocities of 20~km~s-1. Large-scale oscillations are frequently seen in vertical prominences with periods on the order of 10 min and upward propagation speeds of approximately 10~km~s-1. In addition, "bubble" events in which large voids 10--30~Mm across inflate and then burst are seen in some of the vertical prominences. In contrast, the horizontally structured quiescent prominences exhibit only limited flows along the horizontal filaments. We speculate on the origin of the distinction between the vertically and horizontally structured prominences, taking into account viewing angle and the underlying photospheric magnetic flux density. We also discuss the nature of the mysterious dark plumes and bubble expansions and their implications for prominence mass balance in light of recent models of prominence magnetic structure that find vertical flows along some field lines.

Overview of the Frontal Air-Sea Interaction Experiment (FASINEX) - A study of air-sea interaction in a region of strong oceanic gradients

NASA Technical Reports Server (NTRS)

Weller, Robert A.

1991-01-01

From 1984 to 1986 the cooperative Frontal Air-Sea Interaction Experiment (FASINEX) was conducted in the subtropical convergence zone southwest of Bermuda. The overall objective of the experiment was to study air-sea interaction on 1- to 100-km horizontal scales in a region of the open ocean characterized by strong horizontal gradients in upper ocean and sea surface properties. Ocean fronts provided both large spatial gradients in sea surface temperature and strong jetlike flows in the upper ocean. The motivation for and detailed objectives of FASINEX are reviewed. Then the components of the field program are summarized. Finally, selected results are presented in order to provide an overview of the outcome of FASINEX.

Vorticity and divergence in the solar photosphere

NASA Technical Reports Server (NTRS)

Wang, YI; Noyes, Robert W.; Tarbell, Theodore D.; Title, Alan M.

1995-01-01

We have studied an outstanding sequence of continuum images of the solar granulation from Pic du Midi Observatory. We have calculated the horizontal vector flow field using a correlation tracking algorithm, and from this determined three scalar field: the vertical component of the curl; the horizontal divergence; and the horizontal flow speed. The divergence field has substantially longer coherence time and more power than does the curl field. Statistically, curl is better correlated with regions of negative divergence - that is, the vertical vorticity is higher in downflow regions, suggesting excess vorticity in intergranular lanes. The average value of the divergence is largest (i.e., outflow is largest) where the horizontal speed is large; we associate these regions with exploding granules. A numerical simulation of general convection also shows similar statistical differences between curl and divergence. Some individual small bright points in the granulation pattern show large local vorticities.

Trains of large Kelvin-Helmholtz billows observed in the Kuroshio above a seamount

NASA Astrophysics Data System (ADS)

Chang, Ming-Huei; Jheng, Sin-Ya; Lien, Ren-Chieh

2016-08-01

Trains of large Kelvin-Helmholtz (KH) billows within the Kuroshio current at ~230 m depth off southeastern Taiwan and above a seamount were observed by shipboard instruments. The trains of large KH billows were present in a strong shear band along the 0.55 m s-1 isotach within the Kuroshio core; they are presumably produced by flow interactions with the rapidly changing topography. Each individual billow, resembling a cat's eye, had a horizontal length scale of 200 m, a vertical scale of 100 m, and a timescale of 7 min, near the local buoyancy frequency. Overturns were observed frequently in the billow cores and the upper eyelids. The turbulent kinetic energy dissipation rates estimated using the Thorpe scale had an average value of O(10-4) W kg-1 and a maximum value of O(10-3) W kg-1. The turbulence mixing induced by the KH billows may exchange Kuroshio water with the surrounding water masses.

An experimental study of geyser-like flows induced by a pressurized air pocket

NASA Astrophysics Data System (ADS)

Elayeb, I. S.; Leon, A.; Choi, Y.; Alnahit, A. O.

2015-12-01

Previous studies argues that the entrapment of pressurized air pockets within combined sewer systems can produce geyser flows, which is an oscillating jetting of a mixture of gas-liquid flows. To verify that pressurized air pockets can effectively produce geysers, laboratory experiments were conducted. However, past experiments were conducted in relatively small-scale apparatus (i.e. maximum φ2" vertical shaft). This study conducted a set of experiments in a larger apparatus. The experimental setup consists of an upstream head tank, a downstream head tank, a horizontal pipe (46.5ft long, φ6") and a vertical pipe (10ft long, φ6"). The initial condition for the experiments is constant flow discharge through the horizontal pipe. The experiments are initiated by injecting an air pocket with pre-determined volume and pressure at the upstream end of the horizontal pipe. The air pocket propagates through the horizontal pipe until it arrives to the vertical shaft, where it is released producing a geyser-like flow. Three flow rates in the horizontal pipe and three injected air pressures were tested. The variables measured were pressure at two locations in the horizontal pipe and two locations in the vertical pipe. High resolution videos at two regions in the vertical shaft were also recorded. To gain further insights in the physics of air-water interaction, the laboratory experiments were complemented with numerical simulations conducted using a commercial 3D CFD model, previously validated with experiments.

RELATIONSHIPS BETWEEN FLUID VORTICITY, KINETIC HELICITY, AND MAGNETIC FIELD ON SMALL-SCALES (QUIET-NETWORK) ON THE SUN

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sangeetha, C. R.; Rajaguru, S. P., E-mail: crsangeetha@iiap.res.in

We derive horizontal fluid motions on the solar surface over large areas covering the quiet-Sun magnetic network from local correlation tracking of convective granules imaged in continuum intensity and Doppler velocity by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory . From these we calculate the horizontal divergence, the vertical component of vorticity, and the kinetic helicity of fluid motions. We study the correlations between fluid divergence and vorticity, and between vorticity (kinetic helicity) and the magnetic field. We find that the vorticity (kinetic helicity) around small-scale fields exhibits a hemispherical pattern (in sign) similar tomore » that followed by the magnetic helicity of large-scale active regions (containing sunspots). We identify this pattern to be a result of the Coriolis force acting on supergranular-scale flows (both the outflows and inflows), consistent with earlier studies using local helioseismology. Furthermore, we show that the magnetic fields cause transfer of vorticity from supergranular inflow regions to outflow regions, and that they tend to suppress the vortical motions around them when magnetic flux densities exceed about 300 G (from HMI). We also show that such an action of the magnetic fields leads to marked changes in the correlations between fluid divergence and vorticity. These results are speculated to be of importance to local dynamo action (if present) and to the dynamical evolution of magnetic helicity at the small-scale.« less

Subsurface Supergranular Vertical Flows as Measured Using Large Distance Separations in Time-Distance Helioseismology

NASA Technical Reports Server (NTRS)

Duvall, Thomas L., Jr.; Hanasoge, S. M.

2012-01-01

As large-distance rays (say, 10 - 24deg) approach the solar surface approximately vertically, travel times measured from surface pairs for these large separations are mostly sensitive to vertical flows, at least for shallow flows within a few Mm of the solar surface. All previous analyses of supergranulation have used smaller separations and have been hampered by the difficulty of separating the horizontal and vertical flow components. We find that the large-separation travel times associated with supergranulation cannot be studied using the standard phase-speed filters of time-distance helioseismology. These filters, whose use is based upon a refractive model of the perturbations, reduce the resultant travel time signal by at least an order of magnitude at some distances. More effective filters are derived. Modeling suggests that the center-annulus travel-time difference [outward-going time minus inward-going time] in the separation range delta= 10 - 24deg is insensitive to the horizontally diverging flow from the centers of the supergranules and should lead to a constant signal from the vertical flow. Our measurement of this quantity, 5.1+/-0.1 seconds, is constant over the distance range. This magnitude of the signal cannot be caused by the level of upflow at cell centers seen at the photosphere of 10 ms(exp-1) extended in depth. It requires the vertical flow to increase with depth. A simple Gaussian model of the increase with depth implies a peak upward flow of 240 ms(exp-1) at a depth of 2.3 Mm and a peak horizontal flow of 700 ms(exp-1) at a depth of 1.6 Mm.

Comparison of simple, small, full-scale sewage treatment systems in Brazil: UASB-maturation ponds-coarse filter; UASB-horizontal subsurface-flow wetland; vertical-flow wetland (first stage of French system).

PubMed

von Sperling, M

2015-01-01

This paper presents a comparison between three simple sewage treatment lines involving natural processes: (a) upflow anaerobic sludge blanket (UASB) reactor-three maturation ponds in series-coarse rock filter; (b) UASB reactor-horizontal subsurface-flow constructed wetland; and (c) vertical-flow constructed wetlands treating raw sewage (first stage of the French system). The evaluation was based on several years of practical experience with three small full-scale plants receiving the same influent wastewater (population equivalents of 220, 60 and 100 inhabitants) in the city of Belo Horizonte, Brazil. The comparison included interpretation of concentrations and removal efficiencies based on monitoring data (organic matter, solids, nitrogen, phosphorus, coliforms and helminth eggs), together with an evaluation of practical aspects, such as land and volume requirements, sludge production and handling, plant management, clogging and others. Based on an integrated evaluation of all aspects involved, it is worth emphasizing that each system has its own specificities, and no generalization can be made on the best option. The overall conclusion is that the three lines are suitable for sewage treatment in small communities in warm-climate regions.

The impact of resolution on the dynamics of the martian global atmosphere: Varying resolution studies with the MarsWRF GCM

NASA Astrophysics Data System (ADS)

Toigo, Anthony D.; Lee, Christopher; Newman, Claire E.; Richardson, Mark I.

2012-09-01

We investigate the sensitivity of the circulation and thermal structure of the martian atmosphere to numerical model resolution in a general circulation model (GCM) using the martian implementation (MarsWRF) of the planetWRF atmospheric model. We provide a description of the MarsWRF GCM and use it to study the global atmosphere at horizontal resolutions from 7.5° × 9° to 0.5° × 0.5°, encompassing the range from standard Mars GCMs to global mesoscale modeling. We find that while most of the gross-scale features of the circulation (the rough location of jets, the qualitative thermal structure, and the major large-scale features of the surface level winds) are insensitive to horizontal resolution over this range, several major features of the circulation are sensitive in detail. The northern winter polar circulation shows the greatest sensitivity, showing a continuous transition from a smooth polar winter jet at low resolution, to a distinct vertically “split” jet as resolution increases. The separation of the lower and middle atmosphere polar jet occurs at roughly 10 Pa, with the split jet structure developing in concert with the intensification of meridional jets at roughly 10 Pa and above 0.1 Pa. These meridional jets appear to represent the separation of lower and middle atmosphere mean overturning circulations (with the former being consistent with the usual concept of the “Hadley cell”). Further, the transition in polar jet structure is more sensitive to changes in zonal than meridional horizontal resolution, suggesting that representation of small-scale wave-mean flow interactions is more important than fine-scale representation of the meridional thermal gradient across the polar front. Increasing the horizontal resolution improves the match between the modeled thermal structure and the Mars Climate Sounder retrievals for northern winter high latitudes. While increased horizontal resolution also improves the simulation of the northern high latitudes at equinox, even the lowest model resolution considered here appears to do a good job for the southern winter and southern equinoctial pole (although in detail some discrepancies remain). These results suggest that studies of the northern winter jet (e.g., transient waves and cyclogenesis) will be more sensitive to global model resolution that those of the south (e.g., the confining dynamics of the southern polar vortex relevant to studies of argon transport). For surface winds, the major effect of increased horizontal resolution is in the superposition of circulations forced by local-scale topography upon the large-scale surface wind patterns. While passive predictions of dust lifting are generally insensitive to model horizontal resolution when no lifting threshold is considered, increasing the stress threshold produces significantly more lifting in higher resolution simulations with the generation of finer-scale, higher-stress winds due primarily to better-resolved topography. Considering the positive feedbacks expected for radiatively active dust lifting, we expect this bias to increase when such feedbacks are permitted.

Large eddy simulation study of the kinetic energy entrainment by energetic turbulent flow structures in large wind farms

NASA Astrophysics Data System (ADS)

VerHulst, Claire; Meneveau, Charles

2014-02-01

In this study, we address the question of how kinetic energy is entrained into large wind turbine arrays and, in particular, how large-scale flow structures contribute to such entrainment. Previous research has shown this entrainment to be an important limiting factor in the performance of very large arrays where the flow becomes fully developed and there is a balance between the forcing of the atmospheric boundary layer and the resistance of the wind turbines. Given the high Reynolds numbers and domain sizes on the order of kilometers, we rely on wall-modeled large eddy simulation (LES) to simulate turbulent flow within the wind farm. Three-dimensional proper orthogonal decomposition (POD) analysis is then used to identify the most energetic flow structures present in the LES data. We quantify the contribution of each POD mode to the kinetic energy entrainment and its dependence on the layout of the wind turbine array. The primary large-scale structures are found to be streamwise, counter-rotating vortices located above the height of the wind turbines. While the flow is periodic, the geometry is not invariant to all horizontal translations due to the presence of the wind turbines and thus POD modes need not be Fourier modes. Differences of the obtained modes with Fourier modes are documented. Some of the modes are responsible for a large fraction of the kinetic energy flux to the wind turbine region. Surprisingly, more flow structures (POD modes) are needed to capture at least 40% of the turbulent kinetic energy, for which the POD analysis is optimal, than are needed to capture at least 40% of the kinetic energy flux to the turbines. For comparison, we consider the cases of aligned and staggered wind turbine arrays in a neutral atmospheric boundary layer as well as a reference case without wind turbines. While the general characteristics of the flow structures are robust, the net kinetic energy entrainment to the turbines depends on the presence and relative arrangement of the wind turbines in the domain.

TOMOGRAPHY OF PLASMA FLOWS IN THE UPPER SOLAR CONVECTION ZONE USING TIME-DISTANCE INVERSION COMBINING RIDGE AND PHASE-SPEED FILTERING

DOE Office of Scientific and Technical Information (OSTI.GOV)

Svanda, Michal, E-mail: michal@astronomie.cz; Astronomical Institute, Charles University in Prague, Faculty of Mathematics and Physics, V Holesovickach 2, CZ-18000 Prague 8

2013-09-20

The consistency of time-distance inversions for horizontal components of the plasma flow on supergranular scales in the upper solar convection zone is checked by comparing the results derived using two k-{omega} filtering procedures-ridge filtering and phase-speed filtering-commonly used in time-distance helioseismology. I show that both approaches result in similar flow estimates when finite-frequency sensitivity kernels are used. I further demonstrate that the performance of the inversion improves (in terms of a simultaneously better averaging kernel and a lower noise level) when the two approaches are combined together in one inversion. Using the combined inversion, I invert for horizontal flows inmore » the upper 10 Mm of the solar convection zone. The flows connected with supergranulation seem to be coherent only for the top {approx}5 Mm; deeper down there is a hint of change of the convection scales toward structures larger than supergranules.« less

Investigation of the Rocket Induced Flow Field in a Rectangular Duct

NASA Technical Reports Server (NTRS)

Landrum, D. Brian; Thames, Mignon; Parkinson, Doug; Gautney, Serena; Hawk, Clark

1999-01-01

Several tests were performed on a one-sixth scale Rocket Based Combined Cycle (RBCC) engine model at the University of Alabama in Huntsville. The UAH RBCC facility consists of a rectangular duct with a vertical strut mounted in the center. The scaled strut consists of two supersonic rocket nozzles with an embedded vertical turbine between the rocket nozzles. The tests included mass flow, flow visualization and horizontal pressure traverses. The mass flow test indicated a c:hoked condition when the rocket chamber pressure is between 200 psi and 300 psi. The flow visualization tests narrowed the rocket chamber pressure range from, 250 psi to 300 psi. Also, from this t.est, an assumption of a minimum

Convection patterns in a liquid metal under an imposed horizontal magnetic field.

PubMed

Yanagisawa, Takatoshi; Hamano, Yozo; Miyagoshi, Takehiro; Yamagishi, Yasuko; Tasaka, Yuji; Takeda, Yasushi

2013-12-01

We performed laboratory experiments of Rayleigh-Bénard convection with liquid gallium under various intensities of a uniform imposed horizontal magnetic field. An ultrasonic velocity profiling method was used to visualize the spatiotemporal structure of the flows with simultaneous monitoring of the temperature fluctuations in the liquid gallium layer. The explored Rayleigh numbers Ra range from the critical value for onset of convection to 10(5); the Chandrasekhar number Q covers values up to 1100. A regime diagram of the convection patterns was established in relation to the Ra and Q values for a square vessel with aspect ratio 5. We identified five flow regimes: (I) a fluctuating large-scale pattern without rolls, (II) weakly constrained rolls with fluctuations, (III) a continuous oscillation of rolls, (IV) repeated roll number transitions with random reversals of the flow direction, and (V) steady two-dimensional (2D) rolls. These flow regimes are classified by the Ra/Q values, the ratio of the buoyancy to the Lorentz force. Power spectra from the temperature time series indicate that regimes I and II have the features of developed turbulence, while the other regimes do not. The region of steady 2D rolls (Busse balloon) extends to high Ra values in the present setting by a horizontal magnetic field and regime V is located inside the Busse balloon. Concerning the instabilities of the steady 2D rolls, regime III is the traveling wave convection developed from the oscillatory instability. Regime IV can be regarded as a state of phase turbulence, which is induced by intermittent occurrences of the skewed-varicose instability.

Gas-Liquid Flow in Pipelines

DOE Office of Scientific and Technical Information (OSTI.GOV)

Thomas J. Hanratty

A research program was carried out at the University of Illinois in which develops a scientific approach to gas-liquid flows that explains their macroscopic behavior in terms of small scale interactions. For simplicity, fully-developed flows in horizontal and near-horizontal pipes. The difficulty in dealing with these flows is that the phases can assume a variety of configurations. The specific goal was to develop a scientific understanding of transitions from one flow regime to another and a quantitative understanding of how the phases distribute for a give regime. These basic understandings are used to predict macroscopic quantities of interest, such asmore » frictional pressure drop, liquid hold-up, entrainment in annular flow and frequency of slugging in slug flows. A number of scientific issues are addressed. Examples are the rate of atomization of a liquid film, the rate of deposition of drops, the behavior of particles in a turbulent field, the generation and growth of interfacial waves. The use of drag-reducing polymers that change macroscopic behavior by changing small scale interactions was explored.« less

Non-stationary Drainage Flows and Cold Pools in Gentle Terrain

NASA Astrophysics Data System (ADS)

Mahrt, L.

2015-12-01

Previous studies have concentrated on organized topography with well-defined slopes or valleys in an effort to understand the flow dynamics. However, most of the Earth's land surface consists of gentle terrain that is quasi three dimensional. Different scenarios are briefly classified. A network of measurements are analyzed to examine shallow cold pools and drainage flow down the valley which develop for weak ambient wind and relatively clear skies. However, transient modes constantly modulate or intermittently eliminate the cold pool, which makes extraction and analysis of the horizontal structure of the cold pool difficult with traditional analysis methods. Singular value decomposition successfully isolates the effects of large-scale flow from local down-valley cold air drainage within the cold pool in spite of the intermittent nature of this local flow. The traditional concept of a cold pool must be generalized to include cold pool intermittency, complex variation of temperature related to some three-dimensionality and a diffuse cold pool top. Different types of cold pools are classified in terms of the stratification and gradient of potential temperature along the slope. The strength of the cold pool is related to a forcing temperature scale proportional to the net radiative cooling divided by the wind speed above the valley. The scatter is large partly due to nonstationarity of the marginal cold pool in this shallow valley

Black Swans and the Effectiveness of Remediating Groundwater Contamination

NASA Astrophysics Data System (ADS)

Siegel, D. I.; Otz, M. H.; Otz, I.

2013-12-01

Black swans, outliers, dominate science far more than do predictable outcomes. Predictable success constitutes the Black Swan in groundwater remediation. Even the National Research Council concluded that remediating groundwater to drinking water standards has failed in typically complex hydrogeologic settings where heterogeneities and preferential flow paths deflect flow paths obliquely to hydraulic gradients. Natural systems, be they biological or physical, build upon a combination of large-scale regularity coupled to chaos at smaller scales. We show through a review of over 25 case studies that groundwater remediation efforts are best served by coupling parsimonious site characterization to natural and induced geochemical tracer tests to at least know where contamination advects with groundwater in the subsurface. In the majority of our case studies, actual flow paths diverge tens of degrees from anticipated flow paths because of unrecognized heterogeneities in the horizontal direction of transport, let alone the vertical direction. Consequently, regulatory agencies would better serve both the public and the environment by recognizing that long-term groundwater cleanup probably is futile in most hydrogeologic settings except to relaxed standards similar to brownfielding. A Black Swan

Multilevel Cloud Structures above Svalbard

NASA Astrophysics Data System (ADS)

Dörnbrack, Andreas; Pitts, Micheal; Poole, Lamont; Gisinger, Sonja; Maturlli, Marion

2017-04-01

The presentation focusses on the reslts recently published by the authors under the heading "picture of the month" in Monthly Weather Review. The presented picture of the month is a superposition of space-borne lidar observations and high-resolution temperature fields of the ECMWF integrated forecast system (IFS). It displays complex tropospheric and stratospheric clouds in the Arctic winter 2015/16. Near the end of December 2015, the unusual northeastward propagation of warm and humid subtropical air masses as far north as 80°N lifted the tropopause by more than 3 km in 24 h and cooled the stratosphere on a large scale. A widespread formation of thick cirrus clouds near the tropopause and of synoptic-scale polar stratospheric clouds (PSCs) occurred as the temperature dropped below the thresholds for the existence of cloud particles. Additionally, mountain waves were excited by the strong flow at the western edge of the ridge across Svalbard, leading to the formation of mesoscale ice PSCs. The most recent IFS cycle using a horizontal resolution of 8 km globally reproduces the large-scale and mesoscale flow features and leads to a remarkable agreement with the wave structure revealed by the space-borne observations.

Numerical analysis of the wake of a 10kW HAWT

NASA Astrophysics Data System (ADS)

Gong, S. G.; Deng, Y. B.; Xie, G. L.; Zhang, J. P.

2017-01-01

With the rising of wind power industry and the ever-growing scale of wind farm, the research for the wake performance of wind turbine has an important guiding significance for the overall arrangement of wind turbines in the large wind farm. The wake simulation model of 10kW horizontal-axis wind turbine is presented on the basis of Averaged Navier-Stokes (RANS) equations and the RNG k-ε turbulence model for applying to the rotational fluid flow. The sliding mesh technique in ANSYS CFX software is used to solve the coupling equation of velocity and pressure. The characters of the average velocity in the wake zone under rated inlet wind speed and different rotor rotational speeds have been investigated. Based on the analysis results, it is proposed that the horizontal spacing between the wind turbines is less than two times radius of rotor, and its longitudinal spacing is less than five times of radius. And other results have also been obtained, which are of great importance for large wind farms.

Prokaryote genome fluidity: toward a system approach of the mobilome.

PubMed

Toussaint, Ariane; Chandler, Mick

2012-01-01

The importance of horizontal/lateral gene transfer (LGT) in shaping the genomes of prokaryotic organisms has been recognized in recent years as a result of analysis of the increasing number of available genome sequences. LGT is largely due to the transfer and recombination activities of mobile genetic elements (MGEs). Bacterial and archaeal genomes are mosaics of vertically and horizontally transmitted DNA segments. This generates reticulate relationships between members of the prokaryotic world that are better represented by networks than by "classical" phylogenetic trees. In this review we summarize the nature and activities of MGEs, and the problems that presently limit their analysis on a large scale. We propose routes to improve their annotation in the flow of genomic and metagenomic sequences that currently exist and those that become available. We describe network analysis of evolutionary relationships among some MGE categories and sketch out possible developments of this type of approach to get more insight into the role of the mobilome in bacterial adaptation and evolution.

An analytical study on groundwater flow in drainage basins with horizontal wells

NASA Astrophysics Data System (ADS)

Wang, Jun-Zhi; Jiang, Xiao-Wei; Wan, Li; Wang, Xu-Sheng; Li, Hailong

2014-06-01

Analytical studies on release/capture zones are often limited to a uniform background groundwater flow. In fact, for basin-scale problems, the undulating water table would lead to the development of hierarchically nested flow systems, which are more complex than a uniform flow. Under the premise that the water table is a replica of undulating topography and hardly influenced by wells, an analytical solution of hydraulic head is derived for a two-dimensional cross section of a drainage basin with horizontal injection/pumping wells. Based on the analytical solution, distributions of hydraulic head, stagnation points and flow systems (including release/capture zones) are explored. The superposition of injection/pumping wells onto the background flow field leads to the development of new internal stagnation points and new flow systems (including release/capture zones). Generally speaking, the existence of n injection/pumping wells would result in up to n new internal stagnation points and up to 2n new flow systems (including release/capture zones). The analytical study presented, which integrates traditional well hydraulics with the theory of regional groundwater flow, is useful in understanding basin-scale groundwater flow influenced by human activities.

Kelvin-Helmholtz waves in extratropical cyclones passing over mountain ranges: KH Waves in Extratropical Cyclones over Mountain Ranges

DOE Office of Scientific and Technical Information (OSTI.GOV)

Medina, Socorro; Houze, Robert A.

2016-02-19

Kelvin–Helmholtz billows with horizontal scales of 3–4 km have been observed in midlatitude cyclones moving over the Italian Alps and the Oregon Cascades when the atmosphere was mostly statically stable with high amounts of shear and Ri < 0.25. In one case, data from a mobile radar located within a windward facing valley documented a layer in which the shear between down-valley flow below 1.2 km and strong upslope cross-barrier flow above was large. Several episodes of Kelvin–Helmholtz waves were observed within the shear layer. The occurrence of the waves appears to be related to the strength of the shear:more » when the shear attained large values, an episode of billows occurred, followed by a sharp decrease in the shear. The occurrence of large values of shear and Kelvin–Helmholtz billows over two different mountain ranges suggests that they may be important features occurring when extratropical cyclones with statically stable flow pass over mountain ranges.« less

Observations of Seafloor Roughness in a Tidally Modulated Inlet

NASA Astrophysics Data System (ADS)

Lippmann, T. C.; Hunt, J.

2014-12-01

The vertical structure of shallow water flows are influenced by the presence of a bottom boundary layer, which spans the water column for long period waves or mean flows. The nature of the boundary is determined in part by the roughness elements that make up the seafloor, and includes sometimes complex undulations associated with regular and irregular shaped bedforms whose scales range several orders of magnitude from orbital wave ripples (10-1 m) to mega-ripples (100 m) and even larger features (101-103) such as sand waves, bars, and dunes. Modeling efforts often parameterize the effects of roughness elements on flow fields, depending on the complexity of the boundary layer formulations. The problem is exacerbated by the transient nature of bedforms and their large spatial extent and variability. This is particularly important in high flow areas with large sediment transport, such as tidally dominated sandy inlets like New River Inlet, NC. Quantification of small scale seafloor variability over large spatial areas requires the use of mobile platforms that can measure with fine scale (order cm) accuracy in wide swaths. The problem is difficult in shallow water where waves and currents are large, and water clarity is often limited. In this work, we present results from bathymetric surveys obtained with the Coastal Bathymetry Survey System, a personal watercraft equipped with a Imagenex multibeam acoustic echosounder and Applanix POS-MV 320 GPS-aided inertial measurement unit. This system is able to measure shallow water seafloor bathymetry and backscatter intensity with very fine scale (10-1 m) resolution and over relatively large scales (103 m) in the presence of high waves and currents. Wavenumber spectra show that the noise floor of the resolved multibeam bathymetry is on the order of 2.5 - 5 cm in amplitude, depending on water depths ranging 2 - 6 m, and about 30 cm in wavelength. Seafloor roughness elements are estimated from wavenumber spectra across the inlet from bathymetric maps of the seafloor obtained with 10-25 cm horizontal resolution. Implications of the effects of the bottom variability on the vertical structure of the currents will be discussed. This work was supported by ONR and NOAA.

On the scaling features of high-latitude geomagnetic field fluctuations during a large geomagnetic storm

NASA Astrophysics Data System (ADS)

De Michelis, Paola; Federica Marcucci, Maria; Consolini, Giuseppe

2015-04-01

Recently we have investigated the spatial distribution of the scaling features of short-time scale magnetic field fluctuations using measurements from several ground-based geomagnetic observatories distributed in the northern hemisphere. We have found that the scaling features of fluctuations of the horizontal magnetic field component at time scales below 100 minutes are correlated with the geomagnetic activity level and with changes in the currents flowing in the ionosphere. Here, we present a detailed analysis of the dynamical changes of the magnetic field scaling features as a function of the geomagnetic activity level during the well-known large geomagnetic storm occurred on July, 15, 2000 (the Bastille event). The observed dynamical changes are discussed in relationship with the changes of the overall ionospheric polar convection and potential structure as reconstructed using SuperDARN data. This work is supported by the Italian National Program for Antarctic Research (PNRA) - Research Project 2013/AC3.08 and by the European Community's Seventh Framework Programme ([FP7/2007-2013]) under Grant no. 313038/STORM and

Flow and coherent structures around circular cylinders in shallow water

NASA Astrophysics Data System (ADS)

Zeng, Jie; Constantinescu, George

2017-06-01

Eddy-resolving numerical simulations are conducted to investigate the dynamics of the large-scale coherent structures around a circular cylinder in an open channel under very shallow flow conditions where the bed friction significantly affects the wake structure. Results are reported for three test cases, for which the ratio between the cylinder diameter, D, and the channel depth, H, is D/H = 10, 25, and 50, respectively. Simulation results show that a horseshoe vortex system forms in all test cases and the dynamics of the necklace vortices is similar to that during the breakaway sub-regime observed for cases when a laminar horseshoe vortex forms around the base of the cylinder. Given the shallow conditions and turbulence in the incoming channel flow, the necklace vortices occupy a large fraction of the flow depth (they penetrate until the free surface in the shallower cases with D/H = 25 and 50). The oscillations of the necklace vortices become less regular with increasing polar angle magnitude and can induce strong amplification of the bed shear stress beneath their cores. Strong interactions are observed between the legs of the necklace vortices and the eddies shed in the separated shear layers in the cases with D/H = 25 and 50. In these two cases, a vortex-street type wake is formed and strong three-dimensional effects are observed in the near-wake flow. A secondary instability in the form of arrays of co-rotating parallel horizontal vortices develops. Once the roller vortices get away from the cylinder, the horizontal vortices in the array orient themselves along the streamwise direction. This instability is not present for moderately shallow conditions (e.g., D/H ≈ 1) nor for very shallow cases when the wake changes to an unsteady bubble type (e.g., D/H = 50). For cases when this secondary instability is present, the horizontal vortices extend vertically over a large fraction of the flow depth and play an important role in the vertical mixing of fluid situated at the wake edges (e.g., by transporting the near-bed, lower-velocity fluid toward the free surface and vice versa). The largest amplification of the bed shear stress in the near-wake region is observed beneath these horizontal vortices, which means that they would play an important role in promoting bed erosion behind the cylinder in the case of a loose bed. Simulation results suggest that these co-rotating vortices form as a result of the interactions between the legs of the main necklace vortices and the vortical eddies contained into the newly forming roller at the back of the cylinder. The paper also analyzes how D/H affects the separation angle on the cylinder, the size of the recirculation bubble, the bed friction velocity distributions, and turbulence statistics.

From convection rolls to finger convection in double-diffusive turbulence

PubMed Central

Verzicco, Roberto; Lohse, Detlef

2016-01-01

Double-diffusive convection (DDC), which is the buoyancy-driven flow with fluid density depending on two scalar components, is ubiquitous in many natural and engineering environments. Of great interests are scalars' transfer rate and flow structures. Here we systematically investigate DDC flow between two horizontal plates, driven by an unstable salinity gradient and stabilized by a temperature gradient. Counterintuitively, when increasing the stabilizing temperature gradient, the salinity flux first increases, even though the velocity monotonically decreases, before it finally breaks down to the purely diffusive value. The enhanced salinity transport is traced back to a transition in the overall flow pattern, namely from large-scale convection rolls to well-organized vertically oriented salt fingers. We also show and explain that the unifying theory of thermal convection originally developed by Grossmann and Lohse for Rayleigh–Bénard convection can be directly applied to DDC flow for a wide range of control parameters (Lewis number and density ratio), including those which cover the common values relevant for ocean flows. PMID:26699474

Mantle flow tectonics - The influence of a ductile lower crust and implications for the formation of topographic uplands on Venus

NASA Technical Reports Server (NTRS)

Bindschadler, Duane L.; Parmentier, E. Marc

1990-01-01

The crust and mantle of Venus can be represented by a model of a layered structure stratified in both density and viscosity. This structure consists of a brittle-elastic upper crustal layer; a ductile weaker crustal layer; a strong upper mantle layer, about 10 percent denser than the crust; and a weaker substrate, representing the portion of the mantle in which convective flow occurs which is a primary source of large-scale topographic and tectonic features. This paper examines the interactions between these four layers and the mantle flow driven by thermal or compositional variations. Solutions are found for a flow driven by a buoyancy-force distribution within the mantle and by relief at the surface and crust-mantle boundary. It is shown that changes in crustal thickness are driven by vertical normal stresses due to mantle flow and by shear coupling of horizontal mantle flow into the crust.

Modeling the MJO rain rates using parameterized large scale dynamics: vertical structure, radiation, and horizontal advection of dry air

NASA Astrophysics Data System (ADS)

Wang, S.; Sobel, A. H.; Nie, J.

2015-12-01

Two Madden Julian Oscillation (MJO) events were observed during October and November 2011 in the equatorial Indian Ocean during the DYNAMO field campaign. Precipitation rates and large-scale vertical motion profiles derived from the DYNAMO northern sounding array are simulated in a small-domain cloud-resolving model using parameterized large-scale dynamics. Three parameterizations of large-scale dynamics --- the conventional weak temperature gradient (WTG) approximation, vertical mode based spectral WTG (SWTG), and damped gravity wave coupling (DGW) --- are employed. The target temperature profiles and radiative heating rates are taken from a control simulation in which the large-scale vertical motion is imposed (rather than directly from observations), and the model itself is significantly modified from that used in previous work. These methodological changes lead to significant improvement in the results.Simulations using all three methods, with imposed time -dependent radiation and horizontal moisture advection, capture the time variations in precipitation associated with the two MJO events well. The three methods produce significant differences in the large-scale vertical motion profile, however. WTG produces the most top-heavy and noisy profiles, while DGW's is smoother with a peak in midlevels. SWTG produces a smooth profile, somewhere between WTG and DGW, and in better agreement with observations than either of the others. Numerical experiments without horizontal advection of moisture suggest that that process significantly reduces the precipitation and suppresses the top-heaviness of large-scale vertical motion during the MJO active phases, while experiments in which the effect of cloud on radiation are disabled indicate that cloud-radiative interaction significantly amplifies the MJO. Experiments in which interactive radiation is used produce poorer agreement with observation than those with imposed time-varying radiative heating. Our results highlight the importance of both horizontal advection of moisture and cloud-radiative feedback to the dynamics of the MJO, as well as to accurate simulation and prediction of it in models.

Levitation of a drop over a film flow

NASA Astrophysics Data System (ADS)

Sreenivas, K. R.; de, P. K.; Arakeri, Jaywant H.

1999-02-01

A vertical jet of water impinging on a horizontal surface produces a radial film flow followed by a circular hydraulic jump. We report a phenomenon where fairly large (1 ml) drops of liquid levitate just upstream of the jump on a thin air layer between the drop and the film flow. We explain the phenomenon using lubrication theory. Bearing action both in the air film and the water film seems to be necessary to support large drops. Horizontal support is given to the drop by the hydraulic jump. A variety of drop shapes is observed depending on the volume of the drop and liquid properties. We show that interaction of the forces due to gravity, surface tension, viscosity and inertia produces these various shapes.

Lateral fluid flow in a compacting sand-shale sequence: South Caspian basin.

USGS Publications Warehouse

Bredehoeft, J.D.; Djevanshir, R.D.; Belitz, K.R.

1988-01-01

The South Caspian basin contains both sands and shales that have pore-fluid pressures substantially in excess of hydrostatic fluid pressure. Pore-pressure data from the South Caspian basin demonstrate that large differences in excess hydraulic head exist between sand and shale. The data indicate that sands are acting as drains for overlying and underlying compacting shales and that fluid flows laterally through the sand on a regional scale from the basin interior northward to points of discharge. The major driving force for the fluid movement is shale compaction. We present a first- order mathematical analysis in an effort to test if the permeability of the sands required to support a regional flow system is reasonable. The results of the analysis suggest regional sand permeabilities ranging from 1 to 30 md; a range that seems reasonable. This result supports the thesis that lateral fluid flow is occurring on a regional scale within the South Caspian basin. If vertical conduits for flow exist within the basin, they are sufficiently impermeable and do not provide a major outlet for the regional flow system. The lateral fluid flow within the sands implies that the stratigraphic sequence is divided into horizontal units that are hydraulically isolated from one another, a conclusion that has important implications for oil and gas migration.-Authors

Mixing Phenomena in a Bottom Blown Copper Smelter: A Water Model Study

NASA Astrophysics Data System (ADS)

Shui, Lang; Cui, Zhixiang; Ma, Xiaodong; Akbar Rhamdhani, M.; Nguyen, Anh; Zhao, Baojun

2015-03-01

The first commercial bottom blown oxygen copper smelting furnace has been installed and operated at Dongying Fangyuan Nonferrous Metals since 2008. Significant advantages have been demonstrated in this technology mainly due to its bottom blown oxygen-enriched gas. In this study, a scaled-down 1:12 model was set up to simulate the flow behavior for understanding the mixing phenomena in the furnace. A single lance was used in the present study for gas blowing to establish a reliable research technique and quantitative characterisation of the mixing behavior. Operating parameters such as horizontal distance from the blowing lance, detector depth, bath height, and gas flow rate were adjusted to investigate the mixing time under different conditions. It was found that when the horizontal distance between the lance and detector is within an effective stirring range, the mixing time decreases slightly with increasing the horizontal distance. Outside this range, the mixing time was found to increase with increasing the horizontal distance and it is more significant on the surface. The mixing time always decreases with increasing gas flow rate and bath height. An empirical relationship of mixing time as functions of gas flow rate and bath height has been established first time for the horizontal bottom blowing furnace.

Composting plant leachate treatment by a pilot-scale, three-stage, horizontal flow constructed wetland in central Iran.

PubMed

Bakhshoodeh, Reza; Alavi, Nadali; Paydary, Pooya

2017-10-01

Handling and treatment of composting leachate is difficult and poses major burdens on composting facilities. The main goal of this study was to evaluate usage of a three-stage, constructed wetland to treat leachate produced in Isfahan composting facility. A pilot-scale, three-stage, subsurface, horizontal flow constructed wetland, planted with vetiver with a flow rate of 24 L/day and a 15-day hydraulic retention time, was used. Removal of organic matter, ammonia, nitrate, total nitrogen, suspended solids, and several heavy metals from Isfahan composting facility leachate was monitored over a 3-month period. Constructed wetland system was capable of efficiently removing BOD 5 (87.3%), COD (74.5%), ammonia (91.5%), nitrate (87.9%), total nitrogen (87.8%), total suspended solids (85.5%), and heavy metals (ranging from 70 to 90%) from the composting leachate. High contaminant removal efficiencies were achieved, but effluent still failed to meet Iranian standards for treated wastewater. This study shows that although a three-stage horizontal flow constructed wetland planted with vetiver cannot be used alone to treat Isfahan composting facility leachate, but it has the potential to be used as a leachate pre-treatment step, along with another complementary method.

Horizontal heat fluxes over complex terrain computed using a simple mixed-layer model and a numerical model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kimura, Fujio; Kuwagata, Tuneo

1995-02-01

The thermally induced local circulation over a periodic valley is simulated by a two-dimensional numerical model that does-not include condensational processes. During the daytime of a clear, calm day, heat is transported from the mountainous region to the valley area by anabatic wind and its return flow. The specific humidity is, however, transported in an inverse manner. The horizontal exchange rate of sensible heat has a horizontal scale similarity, as long as the horizontal scale is less than a critical width of about 100 km. The sensible heat accumulated in an atmospheric column over an arbitrary point can be estimatedmore » by a simple model termed the uniform mixed-layer model (UML). The model assumes that the potential temperature is both vertically and horizontally uniform in the mixed layer, even over the complex terrain. The UML model is valid only when the horizontal scale of the topography is less than the critical width and the maximum difference in the elevation of the topography is less than about 1500 m. Latent heat is accumulated over the mountainous region while the atmosphere becomes dry over the valley area. When the horizontal scale is close to the critical width, the largest amount of humidity is accumulated during the late afternoon over the mountainous region. 18 refs., 15 figs., 1 tab.« less

Physical and Numerical Model Studies of Cross-flow Turbines Towards Accurate Parameterization in Array Simulations

NASA Astrophysics Data System (ADS)

Wosnik, M.; Bachant, P.

2014-12-01

Cross-flow turbines, often referred to as vertical-axis turbines, show potential for success in marine hydrokinetic (MHK) and wind energy applications, ranging from small- to utility-scale installations in tidal/ocean currents and offshore wind. As turbine designs mature, the research focus is shifting from individual devices to the optimization of turbine arrays. It would be expensive and time-consuming to conduct physical model studies of large arrays at large model scales (to achieve sufficiently high Reynolds numbers), and hence numerical techniques are generally better suited to explore the array design parameter space. However, since the computing power available today is not sufficient to conduct simulations of the flow in and around large arrays of turbines with fully resolved turbine geometries (e.g., grid resolution into the viscous sublayer on turbine blades), the turbines' interaction with the energy resource (water current or wind) needs to be parameterized, or modeled. Models used today--a common model is the actuator disk concept--are not able to predict the unique wake structure generated by cross-flow turbines. This wake structure has been shown to create "constructive" interference in some cases, improving turbine performance in array configurations, in contrast with axial-flow, or horizontal axis devices. Towards a more accurate parameterization of cross-flow turbines, an extensive experimental study was carried out using a high-resolution turbine test bed with wake measurement capability in a large cross-section tow tank. The experimental results were then "interpolated" using high-fidelity Navier--Stokes simulations, to gain insight into the turbine's near-wake. The study was designed to achieve sufficiently high Reynolds numbers for the results to be Reynolds number independent with respect to turbine performance and wake statistics, such that they can be reliably extrapolated to full scale and used for model validation. The end product of this work will be a cross-flow turbine actuator line model to be used as an extension to the OpenFOAM computational fluid dynamics (CFD) software framework, which will likely require modifications to commonly-used dynamic stall models, in consideration of the turbines' high angle of attack excursions during normal operation.

Propagating Waves Transverse to the Magnetic Field in a Solar Prominence

NASA Astrophysics Data System (ADS)

Kucera, Therese A.; Knizhnik, K.; Lopez Ariste, A.; Luna Bennasar, M.; Schmieder, B.; Toot, D.

2013-07-01

We have observed a quiescent prominence with the Hinode Solar Optical Telescope (SOT, in Ca II and H-alpha lines), Sacramento Peak Observatory (in H-alpha, H-beta and Sodium-D lines), and THEMIS/MTR (Télescope Héliographique pour l'Étude du Magnétisme et des Instabilités Solaires/MulTi Raies, providing vector magnetograms), and SDO/AIA (Solar Dynamics Observatory Atmospheric Imaging Assembly, in EUV) over a 4 hour period on 2012 October 10. The small fields of view of SOT, Sac Peak and THEMIS are centered on a large pillar-like prominence footpoint extending towards the surface. This feature appears in the larger field of view of the 304 Å band, as a large, quasi-vertical column with material flowing horizontally on each side. The THEMIS/MTR data indicate that the magnetic field in the pillar is essentially horizontal and the observations in the optical wavelengths show a large number of horizontally aligned features on a much smaller scale than the pillar as a whole. The data are consistent with a model of cool prominence plasma trapped in the dips of horizontal field lines. The SOT and Sac Peak data show what appear to be moving wave pulses. These pulses, which include a Doppler signature, move vertically, perpendicular to the field direction, along quasi-vertical columns. The pulses have a velocity of propagation of about 10 km/s, a period about 260 sec, and a wavelength around 2000 km. We interpret these waves in terms of fast magneto-sonic waves and discuss possible wave drivers.

Nocturnal Reversed Flows Above Parallel Ridges in Perdigão, Portugal

NASA Astrophysics Data System (ADS)

Krishnamurthy, R.; Fernando, H. J.; Leo, L. S.; Vassallo, D.; Hocut, C. M.; Creegan, E.; Rodriguez, C. V.; Palma, J. L.

2017-12-01

Prediction of topographically forced or induced wind events is extremely important for dispersion modeling and wind energy studies in complex terrain. To improve the current understanding of micro-scale processes over complex terrain, a large-scale field experiment was conducted in Perdigão, Portugal from May 1st, 2017 to June 15th, 2017. Measurements over a periodic valley were performed using 52 meteorological met-masts, 30 Doppler Lidars (scanning & vertical profilers), 2 tethered lifting systems and other remote sensing instruments (Sodar-rass, wind profilers & radiometer), and radiosondes were released every 6 hours over the period of study. The observations showed several cases of flow reversals confined to a thin layer of 70 - 100 m above the ridge under stably stratified conditions. These flow reversals were mostly observed during the lee wave formation over the periodic valley. It was observed that the flow reversal occurs predominantly under two atmospheric conditions: a) presence of large recirculation zones on the lee side of the hill causing a pressure gradient between the lee-side floor and the mountain ridge, and b) local change in the horizontal pressure gradient due to differential heating rates of the neighboring valley atmospheres. Microscale flow simulations could capture these observed flow reversals. Based on the network of tower instruments and remote sensing devices, the development, structure and occurrences of the flow reversals are being analyzed and quantified. Since these flow reversals are observed within the rotor swept area of modern wind turbines, they would drastically increase the fatigue loads on wind turbine blades. This presentation will include reversed flow observations from several synchronized scanning Doppler Lidars and meteorological towers and a theoretical framework for reverse flow over parallel valleys.

Observations of Magnetic Evolution and Network Flares Driven by Photospheric Flows in the Quiet Sun

NASA Astrophysics Data System (ADS)

Attie, Raphael; Thompson, Barbara J.

2017-08-01

The quiet Sun may be the biggest laboratory to study physical elementary processes of fundamental importance to space plasma. The advantage is the continuous availability of small-scale events, carrying the hidden microphysics that is responsible for larger-scale phenomena. By small-scale events, we mean spatial dimensions of a few Mm at most, and durations of less than an hour. I present here an attempt to describe and understand the coupling between the photospheric flows, the photospheric magnetic flux, and small-scale energetic transient events. By adapting and improving the highly efficient Balltracking technique for Hinode/SOT data, we relate the fine structures of the supergranular flow fields with the magnetic flux evolution. For studying the dynamics of the latter, and more precisely, the magnetic flux cancellation at sites of energy releases, we applied a new feature tracking algorithm called "Magnetic Balltracking" -- which tracks photospheric magnetic elements -- to high-resolution magnetograms from Hinode/SOT.Using observations of the low corona in soft X-rays with Hinode/XRT, we analyse the triggering mechanism of small-scale network flares. By tracking both the flow fields on the one hand, and the magnetic motions on the other hand, we relate the flows with cancelling magnetic flux. We identify two patterns of horizontal flows that act as catalysts for efficient magnetic reconnection: (i) Funnel-shaped streamlines in which the magnetic flux is carried, and (ii) large-scale vortices (~10 Mm and above) at the network intersections, in which distant magnetic features of opposite polarities seem to be sucked in and ultimately vanish. The excess energy stored in the stressed magnetic field of the vortices is sufficient to power network flares.Prospects for determining the magnetic energy budget in the quiet sun are discussed.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Currie, Laura K., E-mail: lcurrie@astro.ex.ac.uk

Motivated by the significant interaction of convection, rotation, and magnetic field in many astrophysical objects, we investigate the interplay between large-scale flows driven by rotating convection and an imposed magnetic field. We utilize a simple model in two dimensions comprised of a plane layer that is rotating about an axis inclined to gravity. It is known that this setup can result in strong mean flows; we numerically examine the effect of an imposed horizontal magnetic field on such flows. We show that increasing the field strength in general suppresses the time-dependent mean flows, but in some cases it organizes them,more » leading to stronger time-averaged flows. Furthermore, we discuss the effect of the field on the correlations responsible for driving the flows and the competition between Reynolds and Maxwell stresses. A change in behavior is observed when the (fluid and magnetic) Prandtl numbers are decreased. In the smaller Prandtl number regime, it is shown that significant mean flows can persist even when the quenching of the overall flow velocity by the field is relatively strong.« less

Predictive wind turbine simulation with an adaptive lattice Boltzmann method for moving boundaries

NASA Astrophysics Data System (ADS)

Deiterding, Ralf; Wood, Stephen L.

2016-09-01

Operating horizontal axis wind turbines create large-scale turbulent wake structures that affect the power output of downwind turbines considerably. The computational prediction of this phenomenon is challenging as efficient low dissipation schemes are necessary that represent the vorticity production by the moving structures accurately and that are able to transport wakes without significant artificial decay over distances of several rotor diameters. We have developed a parallel adaptive lattice Boltzmann method for large eddy simulation of turbulent weakly compressible flows with embedded moving structures that considers these requirements rather naturally and enables first principle simulations of wake-turbine interaction phenomena at reasonable computational costs. The paper describes the employed computational techniques and presents validation simulations for the Mexnext benchmark experiments as well as simulations of the wake propagation in the Scaled Wind Farm Technology (SWIFT) array consisting of three Vestas V27 turbines in triangular arrangement.

Low-speed wind-tunnel investigation of the longitudinal characteristics of a large-scale variable wing-sweep fighter model in the high-lift configuration

NASA Technical Reports Server (NTRS)

Eckert, W. T.; Maki, R. L.

1973-01-01

The low-speed characteristics of a large-scale model of the U. S. Navy/Grumman F-14A aircraft were studied in tests conducted in the Ames Research Center 40- by 80-Foot Wind Tunnel. The primary purpose of the program was the determination of lift and stability levels and landing approach attitude of the aircraft in its high-lift configuration. Tests were conducted at wing angles of attack between minus 2 deg and 30 deg with zero yaw. Data were taken at Reynolds numbers ranging from 3.48 million to 9.64 million based on a wing mean aerodynamic chord of 7.36 ft. The model configuration was changed as required to show the effects of glove slat, wing slat leading-edge radius, cold flow ducting, flap deflection, direct lift control (spoilers), horizontal tail, speed brake, landing gear and missiles.

The Effect of Large Scale Salinity Gradient on Langmuir Turbulence

NASA Astrophysics Data System (ADS)

Fan, Y.; Jarosz, E.; Yu, Z.; Jensen, T.; Sullivan, P. P.; Liang, J.

2017-12-01

Langmuir circulation (LC) is believed to be one of the leading order causes of turbulent mixing in the upper ocean. It is important for momentum and heat exchange across the mixed layer (ML) and directly impact the dynamics and thermodynamics in the upper ocean and lower atmosphere including the vertical distributions of chemical, biological, optical, and acoustic properties. Based on Craik and Leibovich (1976) theory, large eddy simulation (LES) models have been developed to simulate LC in the upper ocean, yielding new insights that could not be obtained from field observations and turbulent closure models. Due its high computational cost, LES models are usually limited to small domain sizes and cannot resolve large-scale flows. Furthermore, most LES models used in the LC simulations use periodic boundary conditions in the horizontal direction, which assumes the physical properties (i.e. temperature and salinity) and expected flow patterns in the area of interest are of a periodically repeating nature so that the limited small LES domain is representative for the larger area. Using periodic boundary condition can significantly reduce computational effort in problems, and it is a good assumption for isotropic shear turbulence. However, LC is anisotropic (McWilliams et al 1997) and was observed to be modulated by crosswind tidal currents (Kukulka et al 2011). Using symmetrical domains, idealized LES studies also indicate LC could interact with oceanic fronts (Hamlington et al 2014) and standing internal waves (Chini and Leibovich, 2005). The present study expands our previous LES modeling investigations of Langmuir turbulence to the real ocean conditions with large scale environmental motion that features fresh water inflow into the study region. Large scale gradient forcing is introduced to the NCAR LES model through scale separation analysis. The model is applied to a field observation in the Gulf of Mexico in July, 2016 when the measurement site was impacted by large fresh water inflow due to flooding from the Mississippi river. Model results indicate that the strong salinity gradient can reduce the mean flow in the ML and inhibit the turbulence in the planetary boundary layer. The Langmuir cells are also rotated clockwise by the pressure gradient.

Unstructured-grid coastal ocean modelling in Southern Adriatic and Northern Ionian Seas

NASA Astrophysics Data System (ADS)

Federico, Ivan; Pinardi, Nadia; Coppini, Giovanni; Oddo, Paolo

2016-04-01

The Southern Adriatic Northern Ionian coastal Forecasting System (SANIFS) is a short-term forecasting system based on unstructured grid approach. The model component is built on SHYFEM finite element three-dimensional hydrodynamic model. The operational chain exploits a downscaling approach starting from the Mediterranean oceanographic-scale model MFS (Mediterranean Forecasting System, operated by INGV). The implementation set-up has been designed to provide accurate hydrodynamics and active tracer processes in the coastal waters of Southern Eastern Italy (Apulia, Basilicata and Calabria regions), where the model is characterized by a variable resolution in range of 50-500 m. The horizontal resolution is also high in open-sea areas, where the elements size is approximately 3 km. The model is forced: (i) at the lateral open boundaries through a full nesting strategy directly with the MFS (temperature, salinity, non-tidal sea surface height and currents) and OTPS (tidal forcing) fields; (ii) at surface through two alternative atmospheric forcing datasets (ECMWF and COSMOME) via MFS-bulk-formulae. Given that the coastal fields are driven by a combination of both local/coastal and deep ocean forcings propagating along the shelf, the performance of SANIFS was verified first (i) at the large and shelf-coastal scales by comparing with a large scale CTD survey and then (ii) at the coastal-harbour scale by comparison with CTD, ADCP and tide gauge data. Sensitivity tests were performed on initialization conditions (mainly focused on spin-up procedures) and on surface boundary conditions by assessing the reliability of two alternative datasets at different horizontal resolution (12.5 and 7 km). The present work highlights how downscaling could improve the simulation of the flow field going from typical open-ocean scales of the order of several km to the coastal (and harbour) scales of tens to hundreds of meters.

Strongly Stratified Turbulence Wakes and Mixing Produced by Fractal Wakes

NASA Astrophysics Data System (ADS)

Dimitrieva, Natalia; Redondo, Jose Manuel; Chashechkin, Yuli; Fraunie, Philippe; Velascos, David

2017-04-01

This paper describes Shliering and Shadowgraph experiments of the wake induced mixing produced by tranversing a vertical or horizontal fractal grid through the interfase between two miscible fluids at low Atwood and Reynolds numbers. This is a configuration design to models the mixing across isopycnals in stably-stratified flows in many environmental relevant situations (either in the atmosphere or in the ocean. The initial unstable stratification is characterized by a reduced gravity: g' = gΔρ ρ where g is gravity, Δρ being the initial density step and ρ the reference density. Here the Atwood number is A = g' _ 2 g . The topology of the fractal wake within the strong stratification, and the internal wave field produces both a turbulent cascade and a wave cascade, with frecuen parametric resonances, the envelope of the mixing front is found to follow a complex non steady 3rd order polinomial function with a maximum at about 4-5 Brunt-Vaisalla non-dimensional time scales: t/N δ = c1(t/N) + c2g Δρ ρ (t/N)2 -c3(t/N)3. Conductivity probes and Shliering and Shadowgraph visual techniques, including CIV with (Laser induced fluorescence and digitization of the light attenuation across the tank) are used in order to investigate the density gradients and the three-dimensionality of the expanding and contracting wake. Fractal analysis is also used in order to estimate the fastest and slowest growing wavelengths. The large scale structures are observed to increase in wave-length as the mixing progresses, and the processes involved in this increase in scale are also examined.Measurements of the pointwise and horizontally averaged concentrations confirm the picture obtained from past flow visualization studies. They show that the fluid passes through the mixing region with relatively small amounts of molecular mixing,and the molecular effects only dominate on longer time scales when the small scales have penetrated through the large scale structures. The Non-stationary dynamicss and structure of stratified fluid flows around a wedge were also studied based of the fundamental equations set using numerical modeling. Due to breaking of naturally existing background diffusion flux of stratifying agent by an impermeable surface of the wedge a complex multi-level vortex system of compensatory fluid motions is formed around the obstacle. The flow is characterized by a wide range of values of internal scales that are absent in a homogeneous liquid. Numerical solution of the fundamental system with the boundary conditions is constructed using a solver such as stratifiedFoam developed within the frame of the open source computational package OpenFOAM using the finite volume method. The computations were performed in parallel using computing resources of the Scientific Research Supercomputer Complex of MSU (SRCC MSU) and the technological platform UniHUB. The evolution of the flow pattern of the wedge by stratified flow has been demonstrated. The complex structure of the fields of physical quantities and their gradients has been shown. Observed in experiment are multiple flow components, including upstream disturbances, internal waves and the downstream wake with submerged transient vortices well reproduced. Structural elements of flow differ in size and laws of variation in space and time. Rich fine flow structure visualized in vicinity and far from the obstacle. The global efficiency of the mixing process is measured and compared with previous estimates of mixing efficiency.

Wind-tunnel modelling of the tip-speed ratio influence on the wake evolution

NASA Astrophysics Data System (ADS)

Stein, Victor P.; Kaltenbach, Hans-Jakob

2016-09-01

Wind-tunnel measurements on the near-wake evolution of a three bladed horizontal axis wind turbine model (HAWT) in the scale 1:O(350) operating in uniform flow conditions and within a turbulent boundary layer at different tip speed ratios are presented. Operational conditions are chosen to exclude Reynolds number effects regarding the turbulent boundary layer as well as the rotor performance. Triple-wire anemometry is used to measure all three velocity components in the mid-vertical and mid-horizontal plane, covering the range from the near- to the far-wake region. In order to analyse wake properties systematically, power and thrust coefficients of the turbine were measured additionally. It is confirmed that realistic modelling of the wake evolution is not possible in a low-turbulence uniform approach flow. Profiles of mean velocity and turbulence intensity exhibit large deviations between the low-turbulence uniform flow and the turbulent boundary layer, especially in the far-wake region. For nearly constant thrust coefficients differences in the evolution of the near-wake can be identified for tip speed ratios in the range from 6.5 to 10.5. It is shown that with increasing downstream distances mean velocity profiles become indistinguishable whereas for turbulence statistics a subtle dependency on the tip speed ratio is still noticeable in the far-wake region.

Numerical and Experimental Studies of the Natural Convection Flow Within a Horizontal Cylinder Subjected to a Uniformly Cold Wall Boundary Condition. Ph.D. Thesis - Va. Poly. Inst. and State Univ.

NASA Technical Reports Server (NTRS)

Stewart, R. B.

1972-01-01

Numberical solutions are obtained for the quasi-compressible Navier-Stokes equations governing the time dependent natural convection flow within a horizontal cylinder. The early time flow development and wall heat transfer is obtained after imposing a uniformly cold wall boundary condition on the cylinder. Solutions are also obtained for the case of a time varying cold wall boundary condition. Windware explicit differ-encing is used for the numerical solutions. The viscous truncation error associated with this scheme is controlled so that first order accuracy is maintained in time and space. The results encompass a range of Grashof numbers from 8.34 times 10,000 to 7 times 10 to the 7th power which is within the laminar flow regime for gravitationally driven fluid flows. Experiments within a small scale instrumented horizontal cylinder revealed the time development of the temperature distribution across the boundary layer and also the decay of wall heat transfer with time.

Performance comparison and economics analysis of waste stabilization ponds and horizontal subsurface flow constructed wetlands treating domestic wastewater: a case study of the Juja sewage treatment works.

PubMed

Mburu, Njenga; Tebitendwa, Sylvie M; van Bruggen, Johan J A; Rousseau, Diederik P L; Lens, Piet N L

2013-10-15

The performance, effluent quality, land area requirement, investment and operation costs of a full-scale waste stabilization pond (WSP) and a pilot scale horizontal subsurface flow constructed wetland (HSSF-CW) at Jomo Kenyatta University of Agriculture and Technology (JKUAT) were investigated between November 2010 to January 2011. Both systems gave comparable medium to high levels of organic matter and suspended solids removal. However, the WSP showed a better removal for Total Phosphorus (TP) and Ammonium (NH4(+)-N). Based on the population equivalent calculations, the land area requirement per person equivalent of the WSP system was 3 times the area that would be required for the HSSF-CW to treat the same amount of wastewater. The total annual cost estimates consisting of capital, operation and maintenance (O&M) costs were comparable for both systems. However, the evaluation of the capital cost of either system showed that it is largely influenced by the size of the population served, local cost of land and the construction materials involved. Hence, one can select either system in terms of treatment efficiency. When land is available other factor including the volume of wastewater or the investment, and O&M costs determine the technology selection. Copyright © 2013 Elsevier Ltd. All rights reserved.

Three-dimensional evolution of water vapor distributions in the Northern Hemisphere stratosphere as observed by the MLS

NASA Technical Reports Server (NTRS)

Lahoz, W. A.; O'Neill, A.; Carr, E. S.; Harwood, R. S.; Froidevaux, L.; Read, W. G.; Waters, J. W.; Kumer, J. B.; Mergenthaler, J. L.; Roche, A. E.

1994-01-01

The three-dimensional evolution of stratospheric water vapor distributions observed by the Microwave Limb Sounder (MLS) during the period October 1991 - July 1992 is documented. The transport features inferred from the MLS water vapor distributions are corroborated using other dynamical fields, namely, nitrous oxide from the Cryogenic Limb Array Etalon Spectrometer instrument, analyzed winds from the U.K. Meteorological Office (UKMO), UKMO-derived potential vorticity, and the diabatic heating field. By taking a vortex-centered view and an along-track view, the authors observe in great detail the vertical and horizontal structure of the northern winter stratosphere. It is demonstrated that the water vapor distributions show clear signatures of the effects of diabatic descent through isentropic surfaces and quasi-horizontal transport along isentropic surfaces, and that the large-scale winter flow is organized by the interaction between the westerly polar vortex and the Aleutian high.

Wind Tunnel Investigation of the Near-wake Flow Dynamics of a Horizontal Axis Wind Turbine

NASA Astrophysics Data System (ADS)

Hashemi-Tari, P.; Siddiqui, K.; Refan, M.; Hangan, H.

2014-06-01

Experiments conducted in a large wind tunnel set-up investigate the 3D flow dynamics within the near-wake region of a horizontal axis wind turbine. Particle Image Velocimetry (PIV) measurements quantify the mean and turbulent components of the flow field. Measurements are performed in multiple adjacent horizontal planes in order to cover the area behind the rotor in a large radial interval, at several locations downstream of the rotor. The measurements were phase-locked in order to facilitate the re-construction of the threedimensional flow field. The mean velocity and turbulence characteristics clearly correlate with the near-wake vortex dynamics and in particular with the helical structure of the flow, formed immediately behind the turbine rotor. Due to the tip and root vortices, the mean and turbulent characteristics of the flow are highly dependent on the azimuth angle in regions close to the rotor and close to the blade tip and root. Further from the rotor, the characteristics of the flow become phase independent. This can be attributed to the breakdown of the vortical structure of the flow, resulting from the turbulent diffusion. In general, the highest levels of turbulence are observed in shear layer around the tip of the blades, which decrease rapidly downstream. The shear zone grows in the radial direction as the wake moves axially, resulting in velocity recovery toward the centre of the rotor due to momentum transport.

Zonostrophic turbulence

NASA Astrophysics Data System (ADS)

Galperin, Boris; Sukoriansky, Semion; Dikovskaya, Nadejda

2008-12-01

Geostrophic turbulence is a flow regime attained by turbulent, rotating, stably stratified fluids in near-geostrophic balance. When a small-scale forcing is present, flows in this regime may develop an inverse energy cascade. Geostrophic turbulence has been used in geophysical fluid dynamics as a relatively simple model of the large-scale planetary and terrestrial circulations. When the meridional variation of the Coriolis parameter (or a β-effect) is taken into account, the horizontal flow symmetry breaks down giving rise to the emergence of jet flows. In a certain parameter range, a new flow regime comes to life. Its main characteristics include strongly anisotropic kinetic energy spectrum and slowly evolving systems of alternating zonal jets. This regime is a subset of geostrophic turbulence and has been coined zonostrophic turbulence; it can develop both on a β-plane and on the surface of a rotating sphere. This regime was first discovered in computer simulations but later revealed in the laboratory experiments, in the deep terrestrial oceans, and on solar giant planets where it is believed to be the primary physical mechanism responsible for the generation and maintenance of the stable systems of alternating zonal jets. The hallmarks of zonostrophic turbulence are the anisotropic inverse energy cascade and complicated interaction between turbulence and Rossby-Haurwitz waves. Addressing the goals of the conference 'Turbulent Mixing and Beyond' that took place in August 2007 in Trieste, Italy, this paper exposes the regime of zonostrophic turbulence to a wide scientific community, provides a survey of this regime, elaborates its main characteristics, offers novel approaches to describe and understand this phenomenon, and discusses its applicability as a model of the large-scale planetary and terrestrial circulations.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lai, Canhai; Xu, Zhijie; Li, Tingwen

In virtual design and scale up of pilot-scale carbon capture systems, the coupled reactive multiphase flow problem must be solved to predict the adsorber’s performance and capture efficiency under various operation conditions. This paper focuses on the detailed computational fluid dynamics (CFD) modeling of a pilot-scale fluidized bed adsorber equipped with vertical cooling tubes. Multiphase Flow with Interphase eXchanges (MFiX), an open-source multiphase flow CFD solver, is used for the simulations with custom code to simulate the chemical reactions and filtered models to capture the effect of the unresolved details in the coarser mesh for simulations with reasonable simulations andmore » manageable computational effort. Previously developed two filtered models for horizontal cylinder drag, heat transfer, and reaction kinetics have been modified to derive the 2D filtered models representing vertical cylinders in the coarse-grid CFD simulations. The effects of the heat exchanger configurations (i.e., horizontal or vertical) on the adsorber’s hydrodynamics and CO2 capture performance are then examined. The simulation result subsequently is compared and contrasted with another predicted by a one-dimensional three-region process model.« less

Lattice Boltzmann methods applied to large-scale three-dimensional virtual cores constructed from digital optical borehole images of the karst carbonate Biscayne aquifer in southeastern Florida

USGS Publications Warehouse

Michael Sukop,; Cunningham, Kevin J.

2014-01-01

Digital optical borehole images at approximately 2 mm vertical resolution and borehole caliper data were used to create three-dimensional renderings of the distribution of (1) matrix porosity and (2) vuggy megaporosity for the karst carbonate Biscayne aquifer in southeastern Florida. The renderings based on the borehole data were used as input into Lattice Boltzmann methods to obtain intrinsic permeability estimates for this extremely transmissive aquifer, where traditional aquifer test methods may fail due to very small drawdowns and non-Darcian flow that can reduce apparent hydraulic conductivity. Variogram analysis of the borehole data suggests a nearly isotropic rock structure at lag lengths up to the nominal borehole diameter. A strong correlation between the diameter of the borehole and the presence of vuggy megaporosity in the data set led to a bias in the variogram where the computed horizontal spatial autocorrelation is strong at lag distances greater than the nominal borehole size. Lattice Boltzmann simulation of flow across a 0.4 × 0.4 × 17 m (2.72 m3 volume) parallel-walled column of rendered matrix and vuggy megaporosity indicates a high hydraulic conductivity of 53 m s−1. This value is similar to previous Lattice Boltzmann calculations of hydraulic conductivity in smaller limestone samples of the Biscayne aquifer. The development of simulation methods that reproduce dual-porosity systems with higher resolution and fidelity and that consider flow through horizontally longer renderings could provide improved estimates of the hydraulic conductivity and help to address questions about the importance of scale.

Lattice Boltzmann methods applied to large-scale three-dimensional virtual cores constructed from digital optical borehole images of the karst carbonate Biscayne aquifer in southeastern Florida

NASA Astrophysics Data System (ADS)

Sukop, Michael C.; Cunningham, Kevin J.

2014-11-01

Digital optical borehole images at approximately 2 mm vertical resolution and borehole caliper data were used to create three-dimensional renderings of the distribution of (1) matrix porosity and (2) vuggy megaporosity for the karst carbonate Biscayne aquifer in southeastern Florida. The renderings based on the borehole data were used as input into Lattice Boltzmann methods to obtain intrinsic permeability estimates for this extremely transmissive aquifer, where traditional aquifer test methods may fail due to very small drawdowns and non-Darcian flow that can reduce apparent hydraulic conductivity. Variogram analysis of the borehole data suggests a nearly isotropic rock structure at lag lengths up to the nominal borehole diameter. A strong correlation between the diameter of the borehole and the presence of vuggy megaporosity in the data set led to a bias in the variogram where the computed horizontal spatial autocorrelation is strong at lag distances greater than the nominal borehole size. Lattice Boltzmann simulation of flow across a 0.4 × 0.4 × 17 m (2.72 m3 volume) parallel-walled column of rendered matrix and vuggy megaporosity indicates a high hydraulic conductivity of 53 m s-1. This value is similar to previous Lattice Boltzmann calculations of hydraulic conductivity in smaller limestone samples of the Biscayne aquifer. The development of simulation methods that reproduce dual-porosity systems with higher resolution and fidelity and that consider flow through horizontally longer renderings could provide improved estimates of the hydraulic conductivity and help to address questions about the importance of scale.

SUPERGRANULATION AS THE LARGEST BUOYANTLY DRIVEN CONVECTIVE SCALE OF THE SUN

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cossette, Jean-Francois; Rast, Mark P.

The origin of solar supergranulation remains a mystery. Unlike granulation, the size of which is comparable to both the thickness of the radiative boundary layer and local scale-height in the photosphere, supergranulation does not reflect any obvious length scale of the solar convection zone. Moreover, recent observations of flows in the photosphere using Doppler imaging or correlation or feature tracking show a monotonic decrease in horizontal flow power at scales larger than supergranulation. Both local area and global spherical shell simulations of solar convection by contrast show the opposite, an increase in horizontal flow amplitudes to a low wavenumber. Wemore » examine these disparities and investigate how the solar supergranulation may arise as a consequence of nonlocal heat transport by cool diving plumes. Using three-dimensional anelastic simulations with surface driving, we show that the kinetic energy of the largest convective scales in the upper layers of a stratified domain reflects the depth of transition from strong buoyant driving to adiabatic stratification below caused by the dilution of the granular downflows. This depth is quite shallow because of the rapid increase of the mean density below the photosphere. We interpret the observed monotonic decrease in solar convective power at scales larger than supergranulation to be a consequence of this rapid transition, with the supergranular scale the largest buoyantly driven mode of convection in the Sun.« less

Productivity and injectivity of horizontal wells. Quarterly report, October 1--December 31, 1993

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fayers, F.J.; Aziz, K.; Hewett, T.A.

1993-03-10

A number of activities have been carried out in the last three months. A list outlining these efforts is presented below followed by brief description of each activity in the subsequent sections of this report: Progress is being made on the development of a black oil three-phase simulator which will allow the use of a generalized Voronoi grid in the plane perpendicular to a horizontal well. The available analytical solutions in the literature for calculating productivity indices (Inflow Performance) of horizontal wells have been reviewed. The pseudo-steady state analytic model of Goode and Kuchuk has been applied to an examplemore » problem. A general mechanistic two-phase flow model is under development. The model is capable of predicting flow transition boundaries for a horizontal pipe at any inclination angle. It also has the capability of determining pressure drops and holdups for all the flow regimes. A large code incorporating all the features of the model has been programmed and is currently being tested.« less

The nonlinear response in a baroclinic model to stationary forcings

NASA Technical Reports Server (NTRS)

Roads, J. O.

1982-01-01

Wave-zonal-flow interaction studies demonstrated that stratospheric warnings and blocking may be described by the same basic mechanisms although different results occurred for different horizontal scales.

The Role of Rough Topography in Mediating Impacts of Bottom Drag in Eddying Ocean Circulation Models.

PubMed

Trossman, David S; Arbic, Brian K; Straub, David N; Richman, James G; Chassignet, Eric P; Wallcraft, Alan J; Xu, Xiaobiao

2017-08-01

Motivated by the substantial sensitivity of eddies in two-layer quasi-geostrophic (QG) turbulence models to the strength of bottom drag, this study explores the sensitivity of eddies in more realistic ocean general circulation model (OGCM) simulations to bottom drag strength. The OGCM results are interpreted using previous results from horizontally homogeneous, two-layer, flat-bottom, f-plane, doubly periodic QG turbulence simulations and new results from two-layer β -plane QG turbulence simulations run in a basin geometry with both flat and rough bottoms. Baroclinicity in all of the simulations varies greatly with drag strength, with weak drag corresponding to more barotropic flow and strong drag corresponding to more baroclinic flow. The sensitivity of the baroclinicity in the QG basin simulations to bottom drag is considerably reduced, however, when rough topography is used in lieu of a flat bottom. Rough topography reduces the sensitivity of the eddy kinetic energy amplitude and horizontal length scales in the QG basin simulations to bottom drag to an even greater degree. The OGCM simulation behavior is qualitatively similar to that in the QG rough bottom basin simulations in that baroclinicity is more sensitive to bottom drag strength than are eddy amplitudes or horizontal length scales. Rough topography therefore appears to mediate the sensitivity of eddies in models to the strength of bottom drag. The sensitivity of eddies to parameterized topographic internal lee wave drag, which has recently been introduced into some OGCMs, is also briefly discussed. Wave drag acts like a strong bottom drag in that it increases the baroclinicity of the flow, without strongly affecting eddy horizontal length scales.

Use of tandem circulation wells to measure hydraulic conductivity without groundwater extraction

NASA Astrophysics Data System (ADS)

Goltz, Mark N.; Huang, Junqi; Close, Murray E.; Flintoft, Mark J.; Pang, Liping

2008-09-01

Conventional methods to measure the hydraulic conductivity of an aquifer on a relatively large scale (10-100 m) require extraction of significant quantities of groundwater. This can be expensive, and otherwise problematic, when investigating a contaminated aquifer. In this study, innovative approaches that make use of tandem circulation wells to measure hydraulic conductivity are proposed. These approaches measure conductivity on a relatively large scale, but do not require extraction of groundwater. Two basic approaches for using circulation wells to measure hydraulic conductivity are presented; one approach is based upon the dipole-flow test method, while the other approach relies on a tracer test to measure the flow of water between two recirculating wells. The approaches are tested in a relatively homogeneous and isotropic artificial aquifer, where the conductivities measured by both approaches are compared to each other and to the previously measured hydraulic conductivity of the aquifer. It was shown that both approaches have the potential to accurately measure horizontal and vertical hydraulic conductivity for a relatively large subsurface volume without the need to pump groundwater to the surface. Future work is recommended to evaluate the ability of these tandem circulation wells to accurately measure hydraulic conductivity when anisotropy and heterogeneity are greater than in the artificial aquifer used for these studies.

Enhanced Representation of Turbulent Flow Phenomena in Large-Eddy Simulations of the Atmospheric Boundary Layer using Grid Refinement with Pseudo-Spectral Numerics

NASA Astrophysics Data System (ADS)

Torkelson, G. Q.; Stoll, R., II

2017-12-01

Large Eddy Simulation (LES) is a tool commonly used to study the turbulent transport of momentum, heat, and moisture in the Atmospheric Boundary Layer (ABL). For a wide range of ABL LES applications, representing the full range of turbulent length scales in the flow field is a challenge. This is an acute problem in regions of the ABL with strong velocity or scalar gradients, which are typically poorly resolved by standard computational grids (e.g., near the ground surface, in the entrainment zone). Most efforts to address this problem have focused on advanced sub-grid scale (SGS) turbulence model development, or on the use of massive computational resources. While some work exists using embedded meshes, very little has been done on the use of grid refinement. Here, we explore the benefits of grid refinement in a pseudo-spectral LES numerical code. The code utilizes both uniform refinement of the grid in horizontal directions, and stretching of the grid in the vertical direction. Combining the two techniques allows us to refine areas of the flow while maintaining an acceptable grid aspect ratio. In tests that used only refinement of the vertical grid spacing, large grid aspect ratios were found to cause a significant unphysical spike in the stream-wise velocity variance near the ground surface. This was especially problematic in simulations of stably-stratified ABL flows. The use of advanced SGS models was not sufficient to alleviate this issue. The new refinement technique is evaluated using a series of idealized simulation test cases of neutrally and stably stratified ABLs. These test cases illustrate the ability of grid refinement to increase computational efficiency without loss in the representation of statistical features of the flow field.

The influence on response of axial rotation of a six-group local-conductance probe in horizontal oil-water two-phase flow

NASA Astrophysics Data System (ADS)

Weihang, Kong; Lingfu, Kong; Lei, Li; Xingbin, Liu; Tao, Cui

2017-06-01

Water volume fraction is an important parameter of two-phase flow measurement, and it is an urgent task for accurate measurement in horizontal oil field development and optimization of oil production. The previous ring-shaped conductance water-cut meter cannot obtain the response values corresponding to the oil field water conductivity for oil-water two-phase flow in horizontal oil-producing wells characterized by low yield liquid, low velocity and high water cut. Hence, an inserted axisymmetric array structure sensor, i.e. a six-group local-conductance probe (SGLCP), is proposed in this paper. Firstly, the electric field distributions generated by the exciting electrodes of SGLCP are investigated by the finite element method (FEM), and the spatial sensitivity distributions of SGLCP are analyzed from the aspect of different separations between two electrodes and different axial rotation angles respectively. Secondly, the numerical simulation responses of SGLCP in horizontal segregated flow are calculated from the aspect of different water cut and heights of the water layer, respectively. Lastly, an SGLCP-based well logging instrument was developed, and experiments were carried out in a horizontal pipe with an inner diameter of 125 mm on the industrial-scale experimental multiphase flow setup in the Daqing Oilfield, China. In the experiments, the different oil-water two-phase flow, mineralization degree, temperature and pressure were tested. The results obtained from the simulation experiments and simulation well experiments demonstrate that the designed and developed SGLCP-based instrument still has a good response characteristic for measuring water conductivity under the different conditions mentioned above. The validity and reliability of obtaining the response values corresponding to the water conductivity through the designed and developed SGLCP-based instrument are verified by the experimental results. The significance of this work can provide an effective technology for measuring the water volume fraction of oil-water two-phase flow in horizontal oil-producing wells.

Correction of Excessive Precipitation Over Steep and High Mountains in a General Circulation Model

NASA Technical Reports Server (NTRS)

Chao, Winston C.

2012-01-01

Excessive precipitation over steep and high mountains (EPSM) is a well-known problem in GCMs and meso-scale models. This problem impairs simulation and data assimilation products. Among the possible causes investigated in this study, we found that the most important one, by far, is a missing upward transport of heat out of the boundary layer due to the vertical circulations forced by the daytime upslope winds, which are forced by the heated boundary layer on subgrid-scale slopes. These upslope winds are associated with large subgrid-scale topographic variation, which is found over steep and high mountains. Without such subgridscale heat ventilation, the resolvable-scale upslope flow in the boundary layer generated by surface sensible heat flux along the mountain slopes is excessive. Such an excessive resolvablescale upslope flow combined with the high moisture content in the boundary layer results in excessive moisture transport toward mountaintops, which in turn gives rise to EPSM. Other possible causes of EPSM that we have investigated include 1) a poorly-designed horizontal moisture flux in the terrain-following coordinates, 2) the condition for cumulus convection being too easily satisfied at mountaintops, 3) the presence of conditional instability of the computational kind, and 4) the absence of blocked flow drag. These are all minor or inconsequential. We have parameterized the ventilation effects of the subgrid-scale heated-slope-induced vertical circulation (SHVC) by removing heat from the boundary layer and depositing it in layers higher up when the topographic variance exceeds a critical value. Test results using NASA/Goddard's GEOS-S GCM have shown that this largely solved the EPSM problem.

Kinematics and dynamics of Nubia-Somalia divergence along the East African rift

NASA Astrophysics Data System (ADS)

Stamps, Dorothy Sarah

Continental rifting is fundamental to the theory of plate tectonics, yet the force balance driving Earth's largest continental rift system, the East African Rift (EAR), remains debated. The EAR actively diverges the Nubian and Somalian plates spanning ˜5000 km N-S from the Red Sea to the Southwest Indian Ridge and ˜3000 km NW-SE from eastern Congo to eastern Madagascar. Previous studies suggest either lithospheric buoyancy forces or horizontal tractions dominate the force balance acting to rupture East Africa. In this work, we investigate the large-scale dynamics of Nubia-Somalia divergence along the EAR driving present-day kinematics. Because Africa is largely surrounded by spreading ridges, we assume plate-plate interactions are minimal and that the major driving forces are gradients in gravitational potential energy (GPE), which includes the effect of vertical mantle tractions, and horizontal basal tractions arising from viscous coupling to horizontal mantle flow. We quantify a continuous strain rate and velocity field based on kinematic models, an updated GPS velocity solution, and the style of earthquake focal mechanisms, which we use as an observational constraint on surface deformation. We solve the 3D force balance equations and calculate vertically averaged deviatoric stress for a 100 km thick lithosphere constrained by the CRUST2.0 crustal density and thickness model. By comparing vertically integrated deviatoric stress with integrated lithospheric strength we demonstrate forces arising from gradients in gravitational potential energy are insufficient to rupture strong lithosphere, hence weakening mechanisms are required to initiate continental rupture. The next step involves inverting for a stress field boundary condition that is the long-wavelength minimum energy deviatoric stress field required to best-fit the style of our continuous strain rate field in addition to deviatoric stress from gradients in GPE. We infer the stress field boundary condition is an estimate of basal shear stress from viscous coupling to horizontal mantle flow. The stress field boundary condition is small (˜1.6 MPa) compared to deviatoric stress from GPE gradients (8-20 MPa) and does not improve the fit to surface deformation indicators more than 8% when combined with deviatoric stress from GPE gradients. Hence we suggest the style of deformation across the EAR can be explained by forces derived from gradients in GPE. We then calculate dynamic velocities using two types of forward models to solve the instantaneous momentum equations. One method is regional and requires vertically averaged effective viscosity to define lithospheric structure with velocity boundary conditions and a free-slip basal boundary condition. The second is a global model that accounts for a brittle upper crust and viscous mantle lithosphere with velocity boundary conditions imposed at the base of the lithosphere from 5 mantle flow models. With both methods we find deformation driven by internal lithospheric buoyancy forces provides the best-fit to GPS observations of surface velocities on the Somalian plate. We find that any additional contribution from horizontal tractions results in overpredicting surface velocities. This work indicates horizontal mantle flow plays a minimal role in Nubia-Somalia divergence and the EAR is driven largely by gradients in GPE.

Scaling and Numerical Model Evaluation of Snow-Cover Effects on the Generation and Modification of Daytime Mesoscale Circulations.

NASA Astrophysics Data System (ADS)

Segal, M.; Garratt, J. R.; Pielke, R. A.; Ye, Z.

1991-04-01

Consideration of the sensible heat flux characteristics over a snow surface suggests a significant diminution in the magnitude of the flux, compared to that over a snow-free surface under the same environmental conditions. Consequently, the existence of snow-covered mesoscale areas adjacent to snow-free areas produces horizontal thermal gradients in the lower atmosphere during the daytime, possibly resulting in a `snow breeze.' In addition, suppression of the daytime thermally induced upslope flow over snow-covered slopes is likely to occur. The present paper provides scaling and modeling evaluations of these situations, with quantification of the generated and modified circulations. These evaluations suggest that under ideal situations involved with uniform snow cover over large areas, particularly in late winter and early spring, a noticeable `snow breeze' is likely to develop. Additionally: suppression of the daytime thermally induced upslope flow is significant and may even result in a daytime drainage flow. The effects of bare ground patchiness in the snow cover on these circulations are also explored, both for flat terrain and slope-flow situations. A patchiness fraction greater than 0.5 is found to result in a noticeably reduced snow-breeze circulation, while a patchiness fraction of only 0.1 caused the simulated daytime drainage flow over slopes to he reversed.

Guided growth of large-scale, horizontally aligned arrays of single-walled carbon nanotubes and their use in thin-film transistors.

PubMed

Kocabas, Coskun; Hur, Seung-Hyun; Gaur, Anshu; Meitl, Matthew A; Shim, Moonsub; Rogers, John A

2005-11-01

A convenient process for generating large-scale, horizontally aligned arrays of pristine, single-walled carbon nanotubes (SWNTs) is described. The approach uses guided growth, by chemical vapor deposition (CVD), of SWNTs on miscut single-crystal quartz substrates. Studies of the growth reveal important relationships between the density and alignment of the tubes, the CVD conditions, and the morphology of the quartz. Electrodes and dielectrics patterned on top of these arrays yield thin-film transistors that use the SWNTs as effective thin-film semiconductors. The ability to build high-performance devices of this type suggests significant promise for large-scale aligned arrays of SWNTs in electronics, sensors, and other applications.

Fate of estrone in laboratory-scale constructed wetlands

USDA-ARS?s Scientific Manuscript database

A horizontal, subsurface, laboratory-scale constructed wetland (CW) consisting of four cells in series was used to determine the attenuation of the steroid hormone estrone (E1) present in animal wastewater. Liquid swine manure diluted 1:80 with farm pond water and dosed with [14C]E1 flowed through ...

Determination of billows and other turbulent structures, part 4.1A

NASA Technical Reports Server (NTRS)

Rastogi, P. K.

1984-01-01

Billows are regular, wave-like arrays of cross-flow vortices that develop in stratified oceanic or atmospheric flows with large shear. Atmospheric billows can become manifest through condensation. Billows are frequently seen in their characteristic cloud forms in the lower atmosphere. Under suitable viewing conditions, billows can also be seen in noctilucent clouds that form near the polar mesosphere during the summer months. Other turbulent structures -- related to billows -- are the Kelvin-Helmholtz instability (KHI) and cat's eye structures that occur in fully developed turbulent shear flows. Shear flows may contain perturbations at many different horizontal wavelengths and vertical scales. Realistic theoretical models have been constructed to study the stability and growth of these perturbations. The extent to which billows and Kelvin-Helmholtz instability have been observed in the atmosphere with the use of radars is outlined. Most of these observations are confined to the troposphere. Suggestions are made for improved radar experiments that are required to detect these structures at higher altitudes.

Flow of ultra-hot Precambrian orogens and the making of crustal layering in Phanerozoic orogenic plateaux

NASA Astrophysics Data System (ADS)

Chardon, Dominique; Gapais, Denis; Cagnard, Florence; Jayananda, Mudlappa; Peucat, Jean-Jacques

2010-05-01

Reassessment of structural / metamorphic properties of ultra-hot Precambrian orogens and shortening of model weak lithospheres support a syn-convergence flow mode on an orogen scale, with a large component of horizontal finite elongation parallel to the orogen. This orogen-scale flow mode combines distributed shortening, gravity-driven flow, lateral escape, and three-dimensional mass redistribution of buried supracrustal rocks, magmas and migmatites in a thick fluid lower crust. This combination preserves a nearly flat surface and Moho. The upper crust maintains a nearly constant thickness by real-time erosion and near-field clastic sedimentation and by ablation at its base by burial of pop-downs into the lower crust. Steady state regime of these orogens is allowed by activation of an attachment layer that maintains kinematic compatibility between the thin and dominantly plastic upper crust and a thick "water bed" of lower crust. Because very thin lithospheres of orogenic plateaux and Precambrian hot orogens have similar thermomechanical structures, bulk orogenic flow comparable to that governing Precambrian hot orogens should actually operate through today's orogenic plateaux as well. Thus, syn-convergence flow fabrics documented on exposed crustal sections of ancient hot orogens that have not undergone collapse may be used to infer the nature of flow fabrics that are imaged by geophysical techniques beneath orogenic plateaux. We provide a detailed geological perspective on syn-convergence crustal flow in relation to magma emplacement and partial melting on a wide oblique crustal transition of the Neoarchean ultra-hot orogen of Southern India. We document sub-horizontal bulk longitudinal flow of the partially molten lower crust over a protracted period of 60 Ma. Bulk flow results from the interplay of (1) pervasive longitudinal transtensional flow of the partially molten crust, (2) longitudinal coaxial flow on flat fabrics in early plutons, (3) distributed, orogen-normal shortening, (4) emplacement of late prolate shape plutons in the direction of flow, and (5) late, conjugate strike-slip shearing. The macroscopic- to regional scale tectonoplutonic pattern produced by longitudinal flow forms a flat composite anisotropy throughout the lower crust. In the light of GPS data, these results suggest that bulk longitudinal flow accounts for observed deformation of the Tibetan plateau as well as for its seismic structure. This flow mode may be preferred to lateral, east-directed channel flow because it combines both lateral gravity-driven thinning and distributed, orogen-normal shortening of the crust. These results further suggest that lower crustal seismic reflectivity in orogenic belts may not necessarily images fabrics produced by extensional tectonics, as commonly thought, but crustal layering produced by syn-convergence lateral flow.

Sub-grid drag models for horizontal cylinder arrays immersed in gas-particle multiphase flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sarkar, Avik; Sun, Xin; Sundaresan, Sankaran

2013-09-08

Immersed cylindrical tube arrays often are used as heat exchangers in gas-particle fluidized beds. In multiphase computational fluid dynamics (CFD) simulations of large fluidized beds, explicit resolution of small cylinders is computationally infeasible. Instead, the cylinder array may be viewed as an effective porous medium in coarse-grid simulations. The cylinders' influence on the suspension as a whole, manifested as an effective drag force, and on the relative motion between gas and particles, manifested as a correction to the gas-particle drag, must be modeled via suitable sub-grid constitutive relationships. In this work, highly resolved unit-cell simulations of flow around an arraymore » of horizontal cylinders, arranged in a staggered configuration, are filtered to construct sub-grid, or `filtered', drag models, which can be implemented in coarse-grid simulations. The force on the suspension exerted by the cylinders is comprised of, as expected, a buoyancy contribution, and a kinetic component analogous to fluid drag on a single cylinder. Furthermore, the introduction of tubes also is found to enhance segregation at the scale of the cylinder size, which, in turn, leads to a reduction in the filtered gas-particle drag.« less

The stably stratified internal boundary layer for steady and diurnally varying offshore flow

NASA Astrophysics Data System (ADS)

Garratt, J. R.

1987-03-01

A two-dimensional numerical mesoscale model is used to investigate the internal structure and growth of the stably stratified internal boundary layer (IBL) beneath warm, continental air flowing over a cooler sea. Two situations are studied — steady-state and diurnally varying offshore flow. In the steady-state case, vertical profiles of mean quantities and eddy diffusion coefficients ( K) within the IBL show small, but significant, changes with increasing distance from the coast. The top of the IBL is well defined, with large vertical gradients within the layer and a maximum in the coast-normal wind component near the top. Well away from the coast, turbulence, identified by non-zero K, decreases to insignificant levels near the top of the IBL; the IBL itself is characterised by a critical value of the layer-flux Richardson number equal to 0.18. The overall behaviour of the mean profiles is similar to that found in the horizontally homogeneous stable boundary layer over land. A simple physical model is used to relate the depth of the layer h to several relevant physical parameters viz., x, the distance from the coast and U, the large-scale wind (both normal to the coastline) and gδθ/θ, Δθ being the temperature difference between continental mixed-layer air and sea surface, θ is the mean potential temperature and g is the acceleration due to gravity. Excellent agreement with the numerical results is found, with h = 0.014 x 1/2 U ( gδθ/θ)-1/2. In the diurnally varying case, the mean profiles within the IBL show only small differences from the steady-state case, although diurnal variations, particularly in the wind maximum, are evident within a few hundred kilometres of the coast. A mesoscale circulation normal to the coast, and superimposed upon the mean offshore flow, develops seawards of the coastline with maximum vertical velocities about sunset, of depth about 2 km and horizontal scale ≈ 500 km. The circulation is related to the advection, and subsequent decay, of daytime convective turbulence over the sea.

Ocean Research Enabled by Underwater Gliders.

PubMed

Rudnick, Daniel L

2016-01-01

Underwater gliders are autonomous underwater vehicles that profile vertically by changing their buoyancy and use wings to move horizontally. Gliders are useful for sustained observation at relatively fine horizontal scales, especially to connect the coastal and open ocean. In this review, research topics are grouped by time and length scales. Large-scale topics addressed include the eastern and western boundary currents and the regional effects of climate variability. The accessibility of horizontal length scales of order 1 km allows investigation of mesoscale and submesoscale features such as fronts and eddies. Because the submesoscales dominate vertical fluxes in the ocean, gliders have found application in studies of biogeochemical processes. At the finest scales, gliders have been used to measure internal waves and turbulent dissipation. The review summarizes gliders' achievements to date and assesses their future in ocean observation.

Processes governing phytoplankton blooms in estuaries. II: The role of horizontal transport

USGS Publications Warehouse

Lucas, L.V.; Koseff, Jeffrey R.; Monismith, Stephen G.; Cloern, J.E.; Thompson, J.K.

1999-01-01

The development and distribution of phytoplankton blooms in estuaries are functions of both local conditions (i.e. the production-loss balance for a water column at a particular spatial location) and large-scale horizontal transport. In this study, the second of a 2-paper series, we use a depth-averaged hydrodynamic-biological model to identify transport-related mechanisms impacting phytoplankton biomass accumulation and distribution on a system level. We chose South San Francisco Bay as a model domain, since its combination of a deep channel surrounded by broad shoals is typical of drowned-river estuaries. Five general mechanisms involving interaction of horizontal transport with variability in local conditions are discussed. Residual (on the order of days to weeks) transport mechanisms affecting bloom development and location include residence time/export, import, and the role of deep channel regions as conduits for mass transport. Interactions occurring on tidal time scales, i.e. on the order of hours) include the phasing of lateral oscillatory tidal flow relative to temporal changes in local net phytoplankton growth rates, as well as lateral sloshing of shoal-derived biomass into deep channel regions during ebb and back into shallow regions during flood tide. Based on these results, we conclude that: (1) while local conditions control whether a bloom is possible, the combination of transport and spatial-temporal variability in local conditions determines if and where a bloom will actually occur; (2) tidal-time-scale physical-biological interactions provide important mechanisms for bloom development and evolution. As a result of both subtidal and tidal-time-scale transport processes, peak biomass may not be observed where local conditions are most favorable to phytoplankton production, and inherently unproductive areas may be regions of high biomass accumulation.

Vertical length scale selection for pancake vortices in strongly stratified viscous fluids

NASA Astrophysics Data System (ADS)

Godoy-Diana, Ramiro; Chomaz, Jean-Marc; Billant, Paul

2004-04-01

The evolution of pancake dipoles of different aspect ratio is studied in a stratified tank experiment. Two cases are reported here for values of the dipole initial aspect ratio alpha_0 = L_v/L_h (where L_v and L_h are vertical and horizontal length scales, respectively) of alpha_0 = 0.4 (case I) and alpha_0 = 1.2 (case II). In the first case, the usual decay scenario is observed where the dipole diffuses slowly with a growing thickness and a decaying circulation. In case II, we observed a regime where the thickness of the dipole decreases and the circulation in the horizontal mid-plane of the vortices remains constant. We show that this regime where the vertical length scale decreases can be explained by the shedding of two boundary layers at the top and bottom of the dipole that literally peel off vorticity layers. Horizontal advection and vertical diffusion cooperate in this regime and the decrease towards the viscous vertical length scale delta = L_hRe(-1/2) occurs on a time scale alpha_0 Re(1/2) T_A, T_A being the advection time L_h/U. From a scaling analysis of the equations for a stratified viscous fluid in the Boussinesq approximation, two dominant balances depending on the parameter R = ReF_h(2) are discussed, where F_h = U/NL_h is the horizontal Froude number and Re = UL_h/nu is the Reynolds number, U, N and nu being, respectively, the translation speed of the dipole, the Brunt Väisälä frequency and the kinematic viscosity. When R≫ 1 the vertical length scale is determined by buoyancy effects to be of order L_b = U/N. The experiments presented in this paper pertain to the case of small R, where viscous effects govern the selection of the vertical length scale. We show that if initially L_v ≤ delta, the flow diffuses on the vertical (case I), while if L_v ≫ delta (case II), vertically sheared horizontal advection decreases the vertical length scale down to delta. This viscous regime may explain results from experiments and numerical simulations on the late evolution of stratified flows where the decay is observed to be independent of the buoyancy frequency N.

A spatial length scale analysis of turbulent temperature and velocity fluctuations within and above an orchard canopy

USGS Publications Warehouse

Wang, Y.S.; Miller, D.R.; Anderson, D.E.; Cionco, R.M.; Lin, J.D.

1992-01-01

Turbulent flow within and above an almond orchard was measured with three-dimensional wind sensors and fine-wire thermocouple sensors arranged in a horizontal array. The data showed organized turbulent structures as indicated by coherent asymmetric ramp patterns in the time series traces across the sensor array. Space-time correlation analysis indicated that velocity and temperature fluctuations were significantly correlated over a transverse distance more than 4m. Integral length scales of velocity and temperature fluctuations were substantially greater in unstable conditions than those in stable conditions. The coherence spectral analysis indicated that Davenport's geometric similarity hypothesis was satisfied in the lower frequency region. From the geometric similarity hypothesis, the spatial extents of large ramp structures were also estimated with the coherence functions.

Effects of heat exchanger tubes on hydrodynamics and CO 2 capture of a sorbent-based fluidized bed reactor

DOE PAGES

Lai, Canhai; Xu, Zhijie; Li, Tingwen; ...

2017-08-05

In virtual design and scale up of pilot-scale carbon capture systems, the coupled reactive multiphase flow problem must be solved to predict the adsorber's performance and capture efficiency under various operation conditions. This paper focuses on the detailed computational fluid dynamics (CFD) modeling of a pilot-scale fluidized bed adsorber equipped with vertical cooling tubes. Multiphase Flow with Interphase eXchanges (MFiX), an open-source multiphase flow CFD solver, is used for the simulations with custom code to simulate the chemical reactions and filtered sub-grid models to capture the effect of the unresolved details in the coarser mesh for simulations with reasonable accuracymore » and manageable computational effort. Previously developed filtered models for horizontal cylinder drag, heat transfer, and reaction kinetics have been modified to derive the 2D filtered models representing vertical cylinders in the coarse-grid CFD simulations. The effects of the heat exchanger configurations (i.e., horizontal or vertical tubes) on the adsorber's hydrodynamics and CO 2 capture performance are then examined. A one-dimensional three-region process model is briefly introduced for comparison purpose. The CFD model matches reasonably well with the process model while provides additional information about the flow field that is not available with the process model.« less

Noise characteristics of upper surface blown configurations. Experimental program and results

NASA Technical Reports Server (NTRS)

Brown, W. H.; Searle, N.; Blakney, D. F.; Pennock, A. P.; Gibson, J. S.

1977-01-01

An experimental data base was developed from the model upper surface blowing (USB) propulsive lift system hardware. While the emphasis was on far field noise data, a considerable amount of relevant flow field data were also obtained. The data were derived from experiments in four different facilities resulting in: (1) small scale static flow field data; (2) small scale static noise data; (3) small scale simulated forward speed noise and load data; and (4) limited larger-scale static noise flow field and load data. All of the small scale tests used the same USB flap parts. Operational and geometrical variables covered in the test program included jet velocity, nozzle shape, nozzle area, nozzle impingement angle, nozzle vertical and horizontal location, flap length, flap deflection angle, and flap radius of curvature.

Horizontal supergranule-scale motions inferred from TRACE ultraviolet observations of the chromosphere

NASA Astrophysics Data System (ADS)

Tian, H.; Potts, H. E.; Marsch, E.; Attie, R.; He, J.-S.

2010-09-01

Aims: We study horizontal supergranule-scale motions revealed by TRACE observation of the chromospheric emission, and investigate the coupling between the chromosphere and the underlying photosphere. Methods: A highly efficient feature-tracking technique called balltracking has been applied for the first time to the image sequences obtained by TRACE (transition region and coronal explorer) in the passband of white light and the three ultraviolet passbands centered at 1700 Å, 1600 Å, and 1550 Å. The resulting velocity fields have been spatially smoothed and temporally averaged in order to reveal horizontal supergranule-scale motions that may exist at the emission heights of these passbands. Results: We find indeed a high correlation between the horizontal velocities derived in the white-light and ultraviolet passbands. The horizontal velocities derived from the chromospheric and photospheric emission are comparable in magnitude. Conclusions: The horizontal motions derived in the UV passbands might indicate the existence of a supergranule-scale magneto-convection in the chromosphere, which may shed new light on the study of mass and energy supply to the corona and solar wind at the height of the chromosphere. However, it is also possible that the apparent motions reflect the chromospheric brightness evolution as produced by acoustic shocks which might be modulated by the photospheric granular motions in their excitation process, or advected partly by the supergranule-scale flow towards the network while propagating upward from the photosphere. To reach a firm conclusion, it is necessary to investigate the role of granular motions in the excitation of shocks through numerical modeling, and future high-cadence chromospheric magnetograms must be scrutinized.

A dynamic subgrid-scale parameterization of the effective wall stress in atmospheric boundary layer flows over multiscale, fractal-like surfaces

NASA Astrophysics Data System (ADS)

Anderson, William; Meneveau, Charles

2010-05-01

A dynamic subgrid-scale (SGS) parameterization for hydrodynamic surface roughness is developed for large-eddy simulation (LES) of atmospheric boundary layer (ABL) flow over multiscale, fractal-like surfaces. The model consists of two parts. First, a baseline model represents surface roughness at horizontal length-scales that can be resolved in the LES. This model takes the form of a force using a prescribed drag coefficient. This approach is tested in LES of flow over cubes, wavy surfaces, and ellipsoidal roughness elements for which there are detailed experimental data available. Secondly, a dynamic roughness model is built, accounting for SGS surface details of finer resolution than the LES grid width. The SGS boundary condition is based on the logarithmic law of the wall, where the unresolved roughness of the surface is modeled as the product of local root-mean-square (RMS) of the unresolved surface height and an unknown dimensionless model coefficient. This coefficient is evaluated dynamically by comparing the plane-average hydrodynamic drag at two resolutions (grid- and test-filter scale, Germano et al., 1991). The new model is tested on surfaces generated through superposition of random-phase Fourier modes with prescribed, power-law surface-height spectra. The results show that the method yields convergent results and correct trends. Limitations and further challenges are highlighted. Supported by the US National Science Foundation (EAR-0609690).

Study of nitrogen two-phase flow pressure drop in horizontal and vertical orientation

NASA Astrophysics Data System (ADS)

Koettig, T.; Kirsch, H.; Santandrea, D.; Bremer, J.

2017-12-01

The large-scale liquid argon Short Baseline Neutrino Far-detector located at Fermilab is designed to detect neutrinos allowing research in the field of neutrino oscillations. It will be filled with liquid argon and operate at almost ambient pressure. Consequently, its operation temperature is determined at about 87 K. The detector will be surrounded by a thermal shield, which is actively cooled with boiling nitrogen at a pressure of about 2.8 bar absolute, the respective saturation pressure of nitrogen. Due to strict temperature gradient constraints, it is important to study the two-phase flow pressure drop of nitrogen along the cooling circuit of the thermal shield in different orientations of the flow with respect to gravity. An experimental setup has been built in order to determine the two-phase flow pressure drop in nitrogen in horizontal, vertical upward and vertical downward direction. The measurements have been conducted under quasi-adiabatic conditions and at a saturation pressure of 2.8 bar absolute. The mass velocity has been varied in the range of 20 kg·m-2·s-1 to 70 kg·m-2·s-1 and the pressure drop data has been recorded scanning the two-phase region from vapor qualities close to zero up to 0.7. The experimental data will be compared with several established predictions of pressure drop e.g. Mueller-Steinhagen and Heck by using the void fraction correlation of Rouhani.

Shape dependence of slip length on patterned hydrophobic surfaces

NASA Astrophysics Data System (ADS)

Gu, Xiaokun; Chen, Min

2011-08-01

The effects of solid-liquid interfacial shape on the boundary velocity slip of patterned hydrophobic surfaces are investigated. The scaling law in literature is extended to demonstrate the role of such shape, indicating a decrease of the effective slip length with increasing interfacial roughness. A patterned surface with horizontally aligned carbon nanotube arrays reaches an effective slip length of 83 nm, by utilizing large intrinsic slippage of carbon nanotube while keeping away from the negative effects of interfacial curvature through the flow direction. The results emphasize the importance of avoiding the solid-liquid interfacial roughness in low-friction patterned surface design and manufacture.

Fabric Development in a Late-Hercynian Magmatic Strike-Slip Shear Zone in Southern Corsica: Indications of Melt-Supported Large-Scale Deformation Localization

NASA Astrophysics Data System (ADS)

Kruhl, J. H.; Vernon, R. H.

2009-05-01

The calc-alcaline granitoids of the Hercynian Corsica Batholith show a large-scale magmatic flow pattern, outlined by the alignment of large (mm-cm) euhedral feldspar crystals. The trend of the steep magmatic foliation is generally N-S in the northern part of the island, swings to approximately E-W orientation in the central part of the Batholith and back again to approximately N-S orientation in the southern part. This pattern is intensified by large-scale magmatic layering, mainly kilometer long lenses and layers of mafic and intermediate intrusions into the granitoids. On the macro- to micro-scale, magma mingling and mixing are present, reflecting the complex intrusion history and the compositional variability of the Corsica Batholith on different scales. Around the Golf of Valinco, a steep, sinistral magmatic shear zone is represented by E-W trending magmatic layering in mingled dioritic, tonalitic, and granitic magmas - previously misleadingly interpreted as migmatites - and by a magmatic flow foliation formed by the alignment of platy feldspar crystals, as well as amphibole and biotite. Characteristic magmatic structures include multiple thin layering, boudinage, monoclinic folding, melt-injected micro shear zones, and fragmenting and back- veining of dioritic enclaves. The intensity of grain alignment roughly correlates with the thickness of layers. It is low in thick and short boudins and high in cm-thin and cm-m long layers. The monoclinic folds refold the magmatic layering. Flat faces of amphibole and biotite grains are aligned in the axial planes of the folds. The feldspar crystals are locally recrystallized to a few large polygonal grains (up to 1 mm across), and quartz commonly shows chessboard subgrain patterns. No further indications of solid-state deformation are present. Field observations, as well as pattern quantification on variably oriented rock surfaces, indicate variations of crystal alignment and fabric anisotropy in cm- to more than 100m-wide bands parallel to the E-W oriented layering, and various stages of melt-present fragmentation. These variations are interpreted as variations of flow intensity and possibly strain-rate variation. The observations on the macro- as well as the micro-scale point to repeated injection of mafic to felsic magma and crystallization in the presence of a regional stress field. The resulting km-scale sinistral, sub-horizontal synmagmatic shear zone reflects large-scale movements during late-Hercynian crustal reorganization and represents an excellent example of localization of deformation into magma-enriched parts of the continental crust.

Viscoelastic flow in the lower crust after the 1992 landers, california, earthquake

PubMed

Deng; Gurnis; Kanamori; Hauksson

1998-11-27

Space geodesy showed that broad-scale postseismic deformation occurred after the 1992 Landers earthquake. Three-dimensional modeling shows that afterslip can only explain one horizontal component of the postseismic deformation, whereas viscoelastic flow can explain the horizontal and near-vertical displacements. The viscosity of a weak, about 10-km-thick layer, in the lower crust beneath the rupture zone that controls the rebound is about 10(18) pascal seconds. The viscoelastic behavior of the lower crust may help to explain the extensional structures observed in the Basin and Range province and it may be used for the analysis of earthquake hazard.

Mud Flow Characteristics Occurred in Izuoshima Island, 2013

NASA Astrophysics Data System (ADS)

Takebayashi, H.; Egashira, S.; Fujita, M.

2015-12-01

Landslides and mud flows were occurred in the west part of the Izuoshima Island, Japan on 16 October 2013. The Izuoshima Island is a volcanic island and the land surface is covered by the volcanic ash sediment in 1m depth. Hence, the mud flow with high sediment concentration was formed. The laminar layer is formed in the debris flow from the bed to the fluid surface. On the other hand, the laminar flow is restricted near the bed in the mud flow and the turbulence flow is formed on the laminar flow layer. As a result, the equilibrium slope of the mud flow becomes smaller comparing to the debris flow. In this study, the numerical analysis mud flow model considering the effect of turbulence flow on the equilibrium slope of the mud flow is developed. Subsequently, the model is applied to the mud flow occurred in the Izuoshima Island and discussed the applicability of the model and the flow characteristics of the mud flow. The differences of the horizontal flow areas between the simulated results and the field data are compared and it was found that the outline of the horizontal shape of the flow areas is reproduced well. Furthermore, the horizontal distribution of the erosion and deposition area is reproduced by the numerical analysis well except for the residential area (Kandachi area). Kandachi area is judged as the erosion area by the field observation, but the sediment was deposited in the numerical analysis. It is considered that the 1.5hour heavy rain over 100mm/h after the mud flow makes the discrepancy. The difference of the horizontal distribution of the maximum flow surface elevation between the simulated results and the field data are compared and it was found that the simulated flow depth is overestimated slightly, because of the wider erosion area due to the coarse resolution elevation data. The averaged velocity and the depth of the mud flow was enough large to collapse the houses.

Effects of variable thermal diffusivity on the structure of convection

NASA Astrophysics Data System (ADS)

Shcheritsa, O. V.; Getling, A. V.; Mazhorova, O. S.

2018-03-01

The structure of multiscale convection in a thermally stratified plane horizontal fluid layer is investigated by means of numerical simulations. The thermal diffusivity is assumed to produce a thin boundary sublayer convectively much more unstable than the bulk of the layer. The simulated flow is a superposition of cellular structures with three different characteristic scales. In contrast to the largest convection cells, the smaller ones are localised in the upper portion of the layer. The smallest cells are advected by the larger-scale convective flows. The simulated flow pattern qualitatively resembles that observed on the Sun.

Persistent magnetic vortex flow at a supergranular vertex

NASA Astrophysics Data System (ADS)

Requerey, Iker S.; Cobo, Basilio Ruiz; Gošić, Milan; Bellot Rubio, Luis R.

2018-03-01

Context. Photospheric vortex flows are thought to play a key role in the evolution of magnetic fields. Recent studies show that these swirling motions are ubiquitous in the solar surface convection and occur in a wide range of temporal and spatial scales. Their interplay with magnetic fields is poorly characterized, however. Aims: We study the relation between a persistent photospheric vortex flow and the evolution of a network magnetic element at a supergranular vertex. Methods: We used long-duration sequences of continuum intensity images acquired with Hinode and the local correlation-tracking method to derive the horizontal photospheric flows. Supergranular cells are detected as large-scale divergence structures in the flow maps. At their vertices, and cospatial with network magnetic elements, the velocity flows converge on a central point. Results: One of these converging flows is observed as a vortex during the whole 24 h time series. It consists of three consecutive vortices that appear nearly at the same location. At their core, a network magnetic element is also detected. Its evolution is strongly correlated to that of the vortices. The magnetic feature is concentrated and evacuated when it is caught by the vortices and is weakened and fragmented after the whirls disappear. Conclusions: This evolutionary behavior supports the picture presented previously, where a small flux tube becomes stable when it is surrounded by a vortex flow. A movie attached to Fig. 2 is available at http://https://www.aanda.org

Transition from geostrophic turbulence to inertia-gravity waves in the atmospheric energy spectrum.

PubMed

Callies, Jörn; Ferrari, Raffaele; Bühler, Oliver

2014-12-02

Midlatitude fluctuations of the atmospheric winds on scales of thousands of kilometers, the most energetic of such fluctuations, are strongly constrained by the Earth's rotation and the atmosphere's stratification. As a result of these constraints, the flow is quasi-2D and energy is trapped at large scales—nonlinear turbulent interactions transfer energy to larger scales, but not to smaller scales. Aircraft observations of wind and temperature near the tropopause indicate that fluctuations at horizontal scales smaller than about 500 km are more energetic than expected from these quasi-2D dynamics. We present an analysis of the observations that indicates that these smaller-scale motions are due to approximately linear inertia-gravity waves, contrary to recent claims that these scales are strongly turbulent. Specifically, the aircraft velocity and temperature measurements are separated into two components: one due to the quasi-2D dynamics and one due to linear inertia-gravity waves. Quasi-2D dynamics dominate at scales larger than 500 km; inertia-gravity waves dominate at scales smaller than 500 km.

Transition from geostrophic turbulence to inertia–gravity waves in the atmospheric energy spectrum

PubMed Central

Callies, Jörn; Ferrari, Raffaele; Bühler, Oliver

2014-01-01

Midlatitude fluctuations of the atmospheric winds on scales of thousands of kilometers, the most energetic of such fluctuations, are strongly constrained by the Earth’s rotation and the atmosphere’s stratification. As a result of these constraints, the flow is quasi-2D and energy is trapped at large scales—nonlinear turbulent interactions transfer energy to larger scales, but not to smaller scales. Aircraft observations of wind and temperature near the tropopause indicate that fluctuations at horizontal scales smaller than about 500 km are more energetic than expected from these quasi-2D dynamics. We present an analysis of the observations that indicates that these smaller-scale motions are due to approximately linear inertia–gravity waves, contrary to recent claims that these scales are strongly turbulent. Specifically, the aircraft velocity and temperature measurements are separated into two components: one due to the quasi-2D dynamics and one due to linear inertia–gravity waves. Quasi-2D dynamics dominate at scales larger than 500 km; inertia–gravity waves dominate at scales smaller than 500 km. PMID:25404349

Probing Subsurface Flows in NOAA Active Region 12192: Comparison with NOAA 10486

NASA Astrophysics Data System (ADS)

Jain, Kiran; Tripathy, S. C.; Hill, F.

2017-11-01

NOAA Active Region (AR) 12192 is the biggest AR observed in solar cycle 24 so far. This was a long-lived AR that survived for four Carrington rotations (CRs) and exhibited several unusual phenomena. We measure the horizontal subsurface flows in this AR in multiple rotations using the ring-diagram technique of local helioseismology and the Global Oscillation Network Group (GONG+) Dopplergrams, and we investigate how different was the plasma flow in AR 12192 from that in AR 10486. Both regions produced several high M- and X-class flares, but they had different coronal mass ejection (CME) productivity. Our analysis suggests that these ARs had unusually large horizontal flow amplitude with distinctly different directions. While meridional flow in AR 12192 was poleward that supports the flux transport to poles, it was equatorward in AR 10486. Furthermore, there was a sudden increase in the magnitude of estimated zonal flow in shallow layers in AR 12192 during the X3.1 flare; however, it reversed direction in AR 10486 with the X17.2 flare. These flow patterns produced strong twists in horizontal velocity with depth in AR 10486 that persisted throughout the disk passage, as opposed to AR 12192, which produced a twist only after the eruption of the X3.1 flare that disappeared soon after. Our study indicates that the sunspot rotation combined with the reorganization of magnetic field in AR 10486 was not sufficient to decrease the flow energy even after several large flares that might have triggered CMEs. Furthermore, in the absence of sunspot rotation in AR 12192, this reorganization of magnetic field contributed significantly to the substantial release of flow energy after the X3.1 flare.

Sea-town interactions over Marseille: 3D urban boundary layer and thermodynamic fields near the surface

NASA Astrophysics Data System (ADS)

Lemonsu, A.; Pigeon, G.; Masson, V.; Moppert, C.

2006-02-01

3D numerical simulations with the Meso-NH atmospheric model including the Town Energy Balance urban parameterization, are conducted over the south-east of France and the one million inhabitants city of Marseille in the frameworks of the ESCOMPTE-UBL program. The geographic situation of the area is relatively complex, because of the proximity of the Mediterranean Sea and the presence of numerous massifs, inducing complex meteorological flows. The present work is focused on six days of the campaign, characterized by the development of strong summer sea-breeze circulations. A complete evaluation of the model is initially realized at both regional- and city-scales, by using the large available database. The regional evaluation shows a good behavior of the model, during the six days of simulation, either for the parameters near the surface or for the vertical profiles describing the structure of the atmosphere. The urban-scale evaluation indicates that the fine structure of the horizontal fields of air temperature above the city is correctly simulated by the model. A specific attention is then pointed to the 250-m horizontal resolution outputs, focused on the Marseille area, for two days of the campaign. From the study of the vertical structure of the Urban Boundary Layer and the thermodynamic fields near the surface, one underscores the important differences due to the regional and local flows, and the complex interactions that occur between the urban effects and the effects of sea breezes.

Atmospheric energetics in regions of intense convective activity

NASA Technical Reports Server (NTRS)

Fuelberg, H. E.

1977-01-01

Synoptic-scale budgets of kinetic and total potential energy are computed using 3- and 6-h data at nine times from NASA's fourth Atmospheric Variability Experiment (AVE IV). Two intense squall lines occurred during the period. Energy budgets for areas that enclose regions of intense convection are shown to have systematic changes that relate to the life cycles of the convection. Some of the synoptic-scale energy processes associated with the convection are found to be larger than those observed in the vicinity of mature cyclones. Volumes enclosing intense convection are found to have large values of cross-contour conversion of potential to kinetic energy and large horizontal export of kinetic energy. Although small net vertical transport of kinetic energy is observed, values at individual layers indicate large upward transport. Transfer of kinetic energy from grid to subgrid scales of motion occurs in the volumes. Latent heat release is large in the middle and upper troposphere and is thought to be the cause of the observed cyclic changes in the budget terms. Total potential energy is found to be imported horizontally in the lower half of the atmosphere, transported aloft, and then exported horizontally. Although local changes of kinetic energy and total potential energy are small, interaction between volumes enclosing convection with surrounding larger volumes is quite large.

The Østerild Balconies Experiment

NASA Astrophysics Data System (ADS)

Karagali, I.; Mann, J.; Dellwik, E.; Simon, E.; Vasiljevic, N.; Larsen, G. C.

2017-12-01

The New European Wind Atlas (NEWA) is a joint effort of research agencies from eight European countries, co-funded under the ERANET Plus Program. The project is structured around two areas of work: development of dynamical downscaling methodologies and measurement campaigns to validate these methodologies, leading to the creation and publication of a European wind atlas in electronic form. One of the main objectives of the NEWA project is to carry out large-scale field experiments at a high spatial and temporal resolution, and provide a significant upgrade to the experimental databases currently available. The Østerild balconies experiment obtained high resolution measurements over an extended horizontal plane to quantify the effect of unevenly forested terrain on the mean wind field. The experiment was performed, between April 12 and August 12 2016, at the Østerild test station for large wind turbines, where DTU Wind Energy operates two 250 m tall meteorological towers, located in northern Denmark. During the first stage of the measurement campaign, a simple balcony-type structure was installed at a height of 50 m above local ground level (AGL), in each of the towers. A wind scanning lidar was placed on each of the balconies, thus the two instruments were simultaneously scanning an area thereby making it possible to reconstruct the wind field and quantify the wind speed variability in a horizontal plane. During the second phase of the experiment, the balconies were raised at the height of 200 m AGL, in order to verify and study the weaker imprint of surface heterogeneity on mean winds but also to study the large-scale spatial and temporal statistics of winds. The aim of the present study is to present the analyses of the measurements to obtain the horizontal wind field reconstruction and the preliminary results of the average wind flow patterns in relation to the terrain characteristics, during the period of available measurements.

Numerical simulation of turbulent convective flow over wavy terrain

NASA Astrophysics Data System (ADS)

Dörnbrack, A.; Schumann, U.

1993-09-01

By means of a large-eddy simulation, the convective boundary layer is investigated for flows over wavy terrain. The lower surface varies sinusoidally in the downstream direction while remaining constant in the other. Several cases are considered with amplitude δ up to 0.15 H and wavelength λ of H to 8 H, where H is the mean fluid-layer height. At the lower surface, the vertical heat flux is prescribed to be constant and the momentum flux is determined locally from the Monin-Obukhov relationship with a roughness length z o=10-4 H. The mean wind is varied between zero and 5 w *, where w * is the convective velocity scale. After rather long times, the flow structure shows horizontal scales up to 4 H, with a pattern similar to that over flat surfaces at corresponding shear friction. Weak mean wind destroys regular spatial structures induced by the surface undulation at zero mean wind. The surface heating suppresses mean-flow recirculation-regions even for steep surface waves. Short surface waves cause strong drag due to hydrostatic and dynamic pressure forces in addition to frictional drag. The pressure drag increases slowly with the mean velocity, and strongly with δ/ H. The turbulence variances increase mainly in the lower half of the mixed layer for U/w *>2.

Large-scale deformation associated with ridge subduction

USGS Publications Warehouse

Geist, E.L.; Fisher, M.A.; Scholl, D. W.

1993-01-01

Continuum models are used to investigate the large-scale deformation associated with the subduction of aseismic ridges. Formulated in the horizontal plane using thin viscous sheet theory, these models measure the horizontal transmission of stress through the arc lithosphere accompanying ridge subduction. Modelling was used to compare the Tonga arc and Louisville ridge collision with the New Hebrides arc and d'Entrecasteaux ridge collision, which have disparate arc-ridge intersection speeds but otherwise similar characteristics. Models of both systems indicate that diffuse deformation (low values of the effective stress-strain exponent n) are required to explain the observed deformation. -from Authors

Phase-relationships between scales in the perturbed turbulent boundary layer

NASA Astrophysics Data System (ADS)

Jacobi, I.; McKeon, B. J.

2017-12-01

The phase-relationship between large-scale motions and small-scale fluctuations in a non-equilibrium turbulent boundary layer was investigated. A zero-pressure-gradient flat plate turbulent boundary layer was perturbed by a short array of two-dimensional roughness elements, both statically, and under dynamic actuation. Within the compound, dynamic perturbation, the forcing generated a synthetic very-large-scale motion (VLSM) within the flow. The flow was decomposed by phase-locking the flow measurements to the roughness forcing, and the phase-relationship between the synthetic VLSM and remaining fluctuating scales was explored by correlation techniques. The general relationship between large- and small-scale motions in the perturbed flow, without phase-locking, was also examined. The synthetic large scale cohered with smaller scales in the flow via a phase-relationship that is similar to that of natural large scales in an unperturbed flow, but with a much stronger organizing effect. Cospectral techniques were employed to describe the physical implications of the perturbation on the relative orientation of large- and small-scale structures in the flow. The correlation and cospectral techniques provide tools for designing more efficient control strategies that can indirectly control small-scale motions via the large scales.

The effects of divergent and nondivergent winds on the kinetic energy budget of a mid-latitude cyclone - A case study

NASA Technical Reports Server (NTRS)

Chen, T.-C.; Alpert, J. C.; Schlatter, T. W.

1978-01-01

The magnitude of the divergent component of the wind is relatively small compared to that of the nondivergent component in large-scale atmospheric flows; nevertheless, it plays an important role in the case of explosive cyclogenesis examined here. The kinetic energy budget for the life cycle of an intense, developing cyclone over North America is calculated. The principal kinetic energy source is the net horizontal transport across the boundaries of the region enclosing the cyclone. By investigating the relative importance of the divergent and nondivergent wind components in the kinetic energy budget, it was found, as expected, that neglecting the divergent wind component in calculating the magnitude of the kinetic energy is of little consequence, but that the horizontal flux convergence and generation of kinetic energy depend crucially upon the divergent component. Modification of the divergent wind component can result in significant changes in the kinetic energy budget of the synoptic system.

Computational and Experimental Study of the Transient Transport Phenomena in a Full-Scale Twin-Roll Continuous Casting Machine

NASA Astrophysics Data System (ADS)

Xu, Mianguang; Li, Zhongyang; Wang, Zhaohui; Zhu, Miaoyong

2017-02-01

To gain a fundamental understanding of the transient fluid flow in twin-roll continuous casting, the current paper applies both large eddy simulation (LES) and full-scale water modeling experiments to investigate the characteristics of the top free surface, stirring effect of the roll rotation, boundary layer fluctuations, and backflow stability. The results show that, the characteristics of the top free surface and the flow field in the wedge-shaped pool region are quite different with/without the consideration of the roll rotation. The roll rotation decreases the instantaneous fluctuation range of the top free surface, but increases its horizontal velocity. The stirring effect of the roll rotating makes the flow field more homogenous and there exists clear shear flow on the rotating roll surface. The vortex shedding induced by the Kármán Vortex Street from the submerged entry nozzle (SEN) causes the "velocity magnitude wave" and strongly influences the boundary layer stability and the backflow stability. The boundary layer fluctuations or the "velocity magnitude wave" induced by the vortex shedding could give rise to the internal porosity. In strip continuous casting process, the vortex shedding phenomenon indicates that the laminar flow can give rise to instability and that it should be made important in the design of the feeding system and the setting of the operating parameters.

Flow Pattern Identification of Horizontal Two-Phase Refrigerant Flow Using Neural Networks

DTIC Science & Technology

2015-12-31

AFRL-RQ-WP-TP-2016-0079 FLOW PATTERN IDENTIFICATION OF HORIZONTAL TWO-PHASE REFRIGERANT FLOW USING NEURAL NETWORKS (POSTPRINT) Abdeel J...Journal Article Postprint 01 October 2013 – 22 June 2015 4. TITLE AND SUBTITLE FLOW PATTERN IDENTIFICATION OF HORIZONTAL TWO-PHASE REFRIGERANT FLOW USING...networks were used to automatically identify two-phase flow patterns for refrigerant R-134a flowing in a horizontal tube. In laboratory experiments

Impact of freeze-drying, mixing and horizontal transport on water vapor in the upper troposphere and lower stratosphere (UTLS)

NASA Astrophysics Data System (ADS)

Poshyvailo, Liubov; Ploeger, Felix; Müller, Rolf; Tao, Mengchu; Konopka, Paul; Abdoulaye Diallo, Mohamadou; Grooß, Jens-Uwe; Günther, Gebhard; Riese, Martin

2017-04-01

Water vapor in the upper troposphere and lower stratosphere (UTLS) is a key player in the global radiation budget. Therefore, a realistic representation of the water vapor distribution in this region and the involved control processes is critical for climate models, but largely uncertain hitherto. It is known that the extremely low temperatures around the tropical tropopause cause the dominant factor controlling water vapor in the lower stratosphere. Here, we focus on additional processes, such as horizontal transport between tropics and extratropics, small-scale mixing, and freeze-drying. We assess the sensitivities of simulated water vapor in the UTLS from simulations with the Chemical Lagrangian Model of the Stratosphere (CLaMS). CLaMS is a Lagrangian transport model, with a parameterization of small-scale mixing (model diffusion) which is coupled to deformations in the large-scale flow. First, to assess the robustness of water vapor with respect to the meteorological datasets we examine CLaMS driven by ECMWF ERA-Interim and the Japanese 55-year reanalysis. Second, to investigate the effects of small-scale mixing we vary the parameterized mixing strength in the CLaMS model between the reference case with the mixing strength optimized to reproduce atmospheric trace gas observations and a purely advective simulation with parameterized mixing turned off. Also calculation of Lagrangian cold points gives further insight of the processes involved. Third, to assess the effects of horizontal transport between the tropics and extratropics we carry out sensitivity simulations with horizontal transport barriers along latitude circles at the equator, 15°N/S and 35°N/S. Finally, the impact of Antarctic dehydration is estimated from additional sensitivity simulations with switched off freeze-drying in the model at high latitudes of 50°N/S. Our results show that the uncertainty in the tropical tropopause temperatures between current reanalysis datasets causes significant differences in simulated water vapor in the lower stratosphere of about 0.5 ppmv. We further find that small-scale mixing increases troposphere-stratosphere exchange causing moistening of the tropopause region and the tropical stratosphere. Besides, there is an enhancement of water vapor along the subtropical jets, particularly in the Southern hemisphere, and in the Asian monsoon in the UTLS. In the Northern extratropics above about 430K potential temperature, small-scale mixing causes drying by increasing horizontal transport between tropics and extratropics. The negligible effect of a transport barrier along the equator shows that the impact of intrahemispheric exchange on water vapor in the UTLS is very weak. Comparison to simulations with transport barriers in the subtropics, on the other hand, shows the effect of the Asian monsoon in moistening middle and high latitudes and the impact of transported dry air from the tropics towards high latitudes.

An analytical study on nested flow systems in a Tóthian basin with a periodically changing water table

NASA Astrophysics Data System (ADS)

Zhao, Ke-Yu; Jiang, Xiao-Wei; Wang, Xu-Sheng; Wan, Li; Wang, Jun-Zhi; Wang, Heng; Li, Hailong

2018-01-01

Classical understanding on basin-scale groundwater flow patterns is based on Tóth's findings of a single flow system in a unit basin (Tóth, 1962) and nested flow systems in a complex basin (Tóth, 1963), both of which were based on steady state models. Vandenberg (1980) extended Tóth (1962) by deriving a transient solution under a periodically changing water table in a unit basin and examined the flow field distortion under different dimensionless response time, τ∗. Following Vandenberg's (1980) approach, we extended Tóth (1963) by deriving the transient solution under a periodically changing water table in a complex basin and examined the transient behavior of nested flow systems. Due to the effect of specific storage, the flow field is asymmetric with respect to the midline, and the trajectory of internal stagnation points constitutes a non-enclosed loop, whose width decreases when τ∗ decreases. The distribution of the relative magnitude of hydraulic head fluctuation, Δh∗ , is dependent on the horizontal distance away from a divide and the depth below the land surface. In the shallow part, Δh∗ decreases from 1 at the divide to 0 at its neighboring valley under all τ∗, while in the deep part, Δh∗ reaches a threshold, whose value decreases when τ∗ increases. The zones with flowing wells are also found to change periodically. As water table falls, there is a general trend of shrinkage in the area of zones with flowing wells, which has a lag to the declining water table under a large τ∗. Although fluxes have not been assigned in our model, the recharge/discharge flux across the top boundary can be obtained. This study is critical to understand a series of periodically changing hydrogeological phenomena in large-scale basins.

Statistical Analysis of Acoustic Wave Power and Flows around Solar Active Regions

NASA Astrophysics Data System (ADS)

Rabello-Soares, M. Cristina; Bogart, Richard S.; Scherrer, Philip H.

2018-05-01

We analyze the effect of a sunspot in its quiet surroundings applying a helioseismic technique on almost three years of Helioseismic and Magnetic Imager (HMI) observations obtained during solar cycle 24 to further study the sunspot structure below the solar surface. The attenuation of acoustic waves with frequencies lower than 4.2 mHz depends more strongly on the wave direction at a distance of 6°–7° from the sunspot center. The amplification of higher frequency waves is highest 6° away from the active region and is largely independent of the wave’s direction. We observe a mean clockwise flow around active regions, the angular speed of which decreases exponentially with distance and has a coefficient close to ‑0.7 degree‑1. The observed horizontal flow in the direction of the nearby active region agrees with a large-scale circulation around the sunspot in the shape of cylindrical shell. The center of the shell seems to be centered around 7° from the sunspot center, where we observe an inflow close to the surface down to ∼2 Mm, followed by an outflow at deeper layers until at least 7 Mm.

Convective Sedimentation of Colloidal Particles in a Bowl.

PubMed

Stiles; Kagan

1999-08-01

A physical model, which regards a colloidal dispersion as a single fluid continuum, is used to investigate cellular convection accompanying gravitational sedimentation in a hemispherical bowl with a thin cylindrical shaft along its vertical axis of symmetry. We have adapted the stream-function-vorticity form of the Navier-Stokes equations to describe momentum conservation in axially symmetric containers. These hydrodynamic equations have been coupled to the mass balance equation for binary hydrodynamic diffusion in the presence of a vertical gravitational field. Using finite-element software we have solved the equations governing coupled diffusive and hydrodynamic flow. A rapidly intensifying horizontal toroidal vortex develops around the axis of the bowl. This vortex is characterized by downward barycentric flow along the curved surface of the bowl and upward flow in the vicinity of its axis. We find that after a short period of time this large-scale cellular convection associated with the curved boundary of the bowl greatly enhances the rate of sedimentation. Copyright 1999 Academic Press.

Sub-grid drag model for immersed vertical cylinders in fluidized beds

DOE PAGES

Verma, Vikrant; Li, Tingwen; Dietiker, Jean -Francois; ...

2017-01-03

Immersed vertical cylinders are often used as heat exchanger in gas-solid fluidized beds. Computational Fluid Dynamics (CFD) simulations are computationally expensive for large scale systems with bundles of cylinders. Therefore sub-grid models are required to facilitate simulations on a coarse grid, where internal cylinders are treated as a porous medium. The influence of cylinders on the gas-solid flow tends to enhance segregation and affect the gas-solid drag. A correction to gas-solid drag must be modeled using a suitable sub-grid constitutive relationship. In the past, Sarkar et al. have developed a sub-grid drag model for horizontal cylinder arrays based on 2Dmore » simulations. However, the effect of a vertical cylinder arrangement was not considered due to computational complexities. In this study, highly resolved 3D simulations with vertical cylinders were performed in small periodic domains. These simulations were filtered to construct a sub-grid drag model which can then be implemented in coarse-grid simulations. Gas-solid drag was filtered for different solids fractions and a significant reduction in drag was identified when compared with simulation without cylinders and simulation with horizontal cylinders. Slip velocities significantly increase when vertical cylinders are present. Lastly, vertical suspension drag due to vertical cylinders is insignificant however substantial horizontal suspension drag is observed which is consistent to the finding for horizontal cylinders.« less

Initiation of slug flow

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hanratty, T.J.; Woods, B.D.

The initiation of slug flow in a horizontal pipe can be predicted either by considering the stability of a slug or by considering the stability of a stratified flow. Measurements of the shedding rate of slugs are used to define necessary conditions for the existence of a slug. Recent results show that slugs develop from an unstable stratified flow through the evolution of small wavelength waves into large wavelength waves that have the possibility of growing to form a slug. The mechanism appears to be quite different for fluids with viscosities close to water than for fluids with large viscositiesmore » (20 centipoise).« less

Granular flow through an aperture: influence of the packing fraction.

PubMed

Aguirre, M A; De Schant, R; Géminard, J-C

2014-07-01

For the last 50 years, the flow of a granular material through an aperture has been intensely studied in gravity-driven vertical systems (e.g., silos and hoppers). Nevertheless, in many industrial applications, grains are horizontally transported at constant velocity, lying on conveyor belts or floating on the surface of flowing liquids. Unlike fluid flows, that are controlled by the pressure, granular flow is not sensitive to the local pressure but rather to the local velocity of the grains at the outlet. We can also expect the flow rate to depend on the local density of the grains. Indeed, vertical systems are packed in dense configurations by gravity, but, in contrast, in horizontal systems the density can take a large range of values, potentially very small, which may significantly alter the flow rate. In the present article, we study, for different initial packing fractions, the discharge through an orifice of monodisperse grains driven at constant velocity by a horizontal conveyor belt. We report how, during the discharge, the packing fraction is modified by the presence of the outlet, and we analyze how changes in the packing fraction induce variations in the flow rate. We observe that variations of packing fraction do not affect the velocity of the grains at the outlet, and, therefore, we establish that flow-rate variations are directly related to changes in the packing fraction.

Granular flow through an aperture: Influence of the packing fraction

NASA Astrophysics Data System (ADS)

Aguirre, M. A.; De Schant, R.; Géminard, J.-C.

2014-07-01

For the last 50 years, the flow of a granular material through an aperture has been intensely studied in gravity-driven vertical systems (e.g., silos and hoppers). Nevertheless, in many industrial applications, grains are horizontally transported at constant velocity, lying on conveyor belts or floating on the surface of flowing liquids. Unlike fluid flows, that are controlled by the pressure, granular flow is not sensitive to the local pressure but rather to the local velocity of the grains at the outlet. We can also expect the flow rate to depend on the local density of the grains. Indeed, vertical systems are packed in dense configurations by gravity, but, in contrast, in horizontal systems the density can take a large range of values, potentially very small, which may significantly alter the flow rate. In the present article, we study, for different initial packing fractions, the discharge through an orifice of monodisperse grains driven at constant velocity by a horizontal conveyor belt. We report how, during the discharge, the packing fraction is modified by the presence of the outlet, and we analyze how changes in the packing fraction induce variations in the flow rate. We observe that variations of packing fraction do not affect the velocity of the grains at the outlet, and, therefore, we establish that flow-rate variations are directly related to changes in the packing fraction.

Development and performance of an alternative biofilter system.

PubMed

Lee, D H; Lau, A K; Pinder, K L

2001-01-01

Step tracer tests were carried out on lab-scale biofilters to determine the residence time distributions (RTDs) of gases passing through two types of biofilters: a standard biofilter with vertical gas flow and a modified biofilter with horizontal gas flow. Results were used to define the flow patterns in the reactors. "Non-ideal flow" indicates that the flow reactors did not behave like either type of ideal reactor: the perfectly stirred reactor [often called a "continuously stirred tank reactor" (CSTR)] or the plug-flow reactor. The horizontal biofilter with back-mixing was able to accommodate a shorter residence time without the usual requirement of greater biofilter surface area for increased biofiltration efficiency. Experimental results indicated that the first bed of the modified biofilter behaved like two CSTRs in series, while the second bed may be represented by two or three CSTRs in series. Because of the flow baffles used in the horizontal biofilter system, its performance was more similar to completely mixed systems, and hence, it could not be modeled as a plug-flow reactor. For the standard biofilter, the number of CSTRs was found to be between 2 and 9 depending on the airflow rate. In terms of NH3 removal efficiency and elimination capacity, the standard biofilter was not as good as the modified system; moreover, the second bed of the modified biofilter exhibited greater removal efficiency than the first bed. The elimination rate increased as biofilter load increased. An opposite trend was exhibited with respect to removal efficiency.

Bending the law: tidal bending and its effects on ice viscosity and flow

NASA Astrophysics Data System (ADS)

Rosier, S.; Gudmundsson, G. H.

2017-12-01

Many ice shelves are subject to strong ocean tides and, in order to accommodate this vertical motion, the ice must bend within the grounding zone. This tidal bending generates large stresses within the ice, changing its effective viscosity. For a confined ice shelf, this is particularly relevant because the tidal bending stresses occur along the sidewalls, which play an important role in the overall flow regime of the ice shelf. Hence, tidal bending stresses will affect both the mean and time-varying components of ice shelf flow. GPS measurements reveal strong variations in horizontal ice shelf velocities at a variety of tidal frequencies. We show, using full-Stokes viscoelastic modelling, that inclusion of tidal bending within the model accounts for much of the observed tidal modulation of horizontal ice shelf flow. Furthermore, our model shows that in the absence of a vertical tidal forcing, the mean flow of the ice shelf is reduced considerably.

On the role of the Antarctic continent in forcing large-scale circulations in the high southern latitudes

NASA Technical Reports Server (NTRS)

Parish, Thomas R.; Bromwich, David H.; Tzeng, Ren-Yow

1994-01-01

The Antarctic topography and attendant katabatic wind regime appear to play a key role in the climate of the high southern latitudes. During the nonsummer months, persistent and often times intense katabatic winds occur in the lowest few hundred meters of the Antarctic atmosphere. These slope flows transport significant amounts of cold air northward and thereby modify the horizontal pressure field over the high southern latitudes. Three-year seasonal cycle numerical simulations using the NCAR Community Climate Model Version 1 (CCM1) with and without representation of the Antarctic orography were performed to explore the role of the elevated terrain and drainage flows on the distribution and evolution of the horizontal pressure field. The katabatic wind regime is an important part of a clearly defined mean meridional circulation in the high southern latitudes. The position and intensity of the attendant sea level low pressure belt appears to be tied to the Antarctic orography. The seasonal movement of mass in the high southern latitudes is therefore constrained by the presence of the Antarctic ice sheet. The semiannual oscillation of pressure over Antarctica and the high southern latitutdes is well depicted in the CCM1 only when the Antarctic orography is included.

Rheology of the lithosphere and the folding caused by horizontal compression

NASA Astrophysics Data System (ADS)

Birger, B. I.

2015-05-01

The laboratory tests of rock specimens show that transient creep, at which deformations increase with time whereas strain rate decreases occurs when creep strains are sufficiently small. Since plate tectonics only permits small deformations in the lithospheric plates, the creep of the lithosphere is transient (non-steady-state). In this work, we study how the rheology of the lithosphere that possesses elasticity, brittleness (pseudo-plasticity), and creep affects the folding in the Earth's crust. Folding is caused by horizontal compression that results from the collision between the lithospheric plates. The effective viscosity characterizing the transient creep is lower than in the case of a steady-state creep and depends on the characteristic time of the considered process. The allowance for transient creep gives the distribution of the rheological properties of the horizontally compressed lithosphere in which the upper crust is brittle, whereas the lower crust and mantle lithosphere are dominated by transient creep. It is shown that the flows that arise in the lithosphere due to the instability under horizontal compression and cause folding are small-scale. These flows are concentrated in the upper brittle crust, they determine the short-wave Earth's surface topography, penetrate into the lower, creep-dominated crust to a shallow depth, and do not penetrate into the mantle. Therefore, these flows do not deform the Moho.

Rotary balance data for a typical single-engine general aviation design for an angle-of-attack range of 8 deg to 90 deg. 2: Influence of horizontal tail location for Model D

NASA Technical Reports Server (NTRS)

Barnhart, B.

1982-01-01

The influence of horizontal tail location on the rotational flow aerodynamics is discussed for a 1/6-scale general aviation airplane model. The model was tested using various horizontal tail positions, with both a high and a low-wing location and for each of two body lengths. Data were measured, using a rotary balance, over an angle-of-attack range of 8 to 90 deg, and for clockwise and counter-clockwise rotations covering an Omega b/2V range of 0 to 0.9.

Experimental investigation of certain internal condensing and boiling flows: Their sensitivity to pressure fluctuations and heat transfer enhancements

NASA Astrophysics Data System (ADS)

Kivisalu, Michael Toomas

Space-based (satellite, scientific probe, space station, etc.) and millimeter -- to -- micro-scale (such as are used in high power electronics cooling, weapons cooling in aircraft, etc.) condensers and boilers are shear/pressure driven. They are of increasing interest to system engineers for thermal management because flow boilers and flow condensers offer both high fluid flow-rate-specific heat transfer capacity and very low thermal resistance between the fluid and the heat exchange surface, so large amounts of heat may be removed using reasonably-sized devices without the need for excessive temperature differences. However, flow stability issues and degredation of performance of shear/pressure driven condensers and boilers due to non-desireable flow morphology over large portions of their lengths have mostly prevented their use in these applications. This research is part of an ongoing investigation seeking to close the gap between science and engineering by analyzing two key innovations which could help address these problems. First, it is recommended that the condenser and boiler be operated in an innovative flow configuration which provides a non-participating core vapor stream to stabilize the annular flow regime throughout the device length, accomplished in an energy-efficient manner by means of ducted vapor re-circulation. This is demonstrated experimentally.. Second, suitable pulsations applied to the vapor entering the condenser or boiler (from the re-circulating vapor stream) greatly reduce the thermal resistance of the already effective annular flow regime. For experiments reported here, application of pulsations increased time-averaged heat-flux up to 900 % at a location within the flow condenser and up to 200 % at a location within the flow boiler, measured at the heat-exchange surface. Traditional fully condensing flows, reported here for comparison purposes, show similar heat-flux enhancements due to imposed pulsations over a range of frequencies. Shear/pressure driven condensing and boiling flow experiments are carried out in horizontal mm-scale channels with heat exchange through the bottom surface. The sides and top of the flow channel are insulated. The fluid is FC-72 from 3M Corporation.

Multiscale limited penetrable horizontal visibility graph for analyzing nonlinear time series

NASA Astrophysics Data System (ADS)

Gao, Zhong-Ke; Cai, Qing; Yang, Yu-Xuan; Dang, Wei-Dong; Zhang, Shan-Shan

2016-10-01

Visibility graph has established itself as a powerful tool for analyzing time series. We in this paper develop a novel multiscale limited penetrable horizontal visibility graph (MLPHVG). We use nonlinear time series from two typical complex systems, i.e., EEG signals and two-phase flow signals, to demonstrate the effectiveness of our method. Combining MLPHVG and support vector machine, we detect epileptic seizures from the EEG signals recorded from healthy subjects and epilepsy patients and the classification accuracy is 100%. In addition, we derive MLPHVGs from oil-water two-phase flow signals and find that the average clustering coefficient at different scales allows faithfully identifying and characterizing three typical oil-water flow patterns. These findings render our MLPHVG method particularly useful for analyzing nonlinear time series from the perspective of multiscale network analysis.

Statistical properties of gravity-driven granular discharge flow under the influence of an obstacle

NASA Astrophysics Data System (ADS)

Endo, Keita; Katsuragi, Hiroaki

2017-06-01

Two-dimensional granular discharge flow driven by gravity under the influence of an obstacle is experimentally investigated. A horizontal exit of width W is opened at the bottom of vertical Hele-Shaw cell filled with stainless-steel particles to start the discharge flow. In this experiment, a circular obstacle is placed in front of the exit. Thus, the distance between the exit and obstacle L is also an important parameter. During the discharge, granular-flow state is acquired by a high-speed camera. The bulk discharge-flow rate is also measured by load cell sensors. The obtained high-speed-image data are analyzed to clarify the particle-level granular-flow dynamics. Using the measured data, we find that the obstacle above the exit affects the granular- flow field. Specifically, the existence of obstacle results in large horizontal granular temperature and small packing fraction. This tendency becomes significant when L is smaller than approximately 6Dg when W ≃ 4Dg, where Dg is diameter of particles.

Mixed Convection Flow in Horizontal CVD Reactors

NASA Astrophysics Data System (ADS)

Chiu, Wilson K. S.; Richards, Cristy J.; Jaluria, Yogesh

1998-11-01

Increasing demands for high quality films and production rates are challenging current Chemical Vapor Deposition (CVD) technology. Since film quality and deposition rates are strongly dependent on gas flow and heat transfer (W.K.S. Chiu and Y. Jaluria, ASME HTD-Vol. 347, pp. 293-311, 1997.), process improvement is obtained through the study of mixed convection flow and temperature distribution in a CVD reactor. Experimental results are presented for a CVD chamber with a horizontal or inclined resistance heated susceptor. Vaporized glycol solution illuminated by a light sheet is used for flow visualization. Temperature measurements are obtained by inserting thermocouple probes into the gas stream or embedding probes into the reactor walls. Flow visualization and temperature measurements show predominantly two dimensional flow and temperature distributions along the streamwise direction under forced convection conditions. Natural convection dominates under large heating rates and low flow rates. Over the range of parameters studied, several distinct flow regimes, characterized by instability, separation, and turbulence, are evident. Different flow regimes alter the flow pattern and temperature distribution, and in consequence, significantly modify deposition rates and uniformity.

Observations of pockmark flow structure in Belfast Bay, Maine, Part 2: evidence for cavity flow

USGS Publications Warehouse

Fandel, Christina L.; Lippmann, Thomas C.; Foster, Diane L.; Brothers, Laura L.

2017-01-01

Pockmark flow circulation patterns were investigated through current measurements along the rim and center of two pockmarks in Belfast Bay, Maine. Observed time-varying current profiles have a complex vertical and directional structure that rotates significantly with depth and is strongly dependent on the phase of the tide. Observations of the vertical profiles of horizontal velocities in relation to relative geometric parameters of the pockmark are consistent with circulation patterns described qualitatively by cavity flow models (Ashcroft and Zhang 2005). The time-mean behavior of the shear layer is typically used to characterize cavity flow, and was estimated using vorticity thickness to quantify the growth rate of the shear layer horizontally across the pockmark. Estimated positive vorticity thickness spreading rates are consistent with cavity flow predictions, and occur at largely different rates between the two pockmarks. Previously modeled flow (Brothers et al. 2011) and laboratory measurements (Pau et al. 2014) over pockmarks of similar geometry to those examined herein are also qualitatively consistent with cavity flow circulation, suggesting that cavity flow may be a good first-order flow model for pockmarks in general.

Scale dependence of the 200-mb divergence inferred from EOLE data.

NASA Technical Reports Server (NTRS)

Morel, P.; Necco, G.

1973-01-01

The EOLE experiment with 480 constant-volume balloons distributed over the Southern Hemisphere approximately at the 200-mb level, has provided a unique, highly accurate set of tracer trajectories in the general westerly circulation. The trajectories of neighboring balloons are analyzed to estimate the horizontal divergence from the Lagrangian derivative of the area of one cluster. The variance of the divergence estimates results from two almost comparable effects: the true divergence of the horizontal flow and eddy diffusion due to small-scale, two-dimensional turbulence. Taking this into account, the rms divergence is found to be of the order of 0.00001 per sec and decreases logarithmically with cluster size. This scale dependence is shown to be consistent with the quasi-geostrophic turbulence model of the general circulation in midlatitudes.

A study of two subgrid-scale models and their effects on wake breakdown behind a wind turbine in uniform inflow

NASA Astrophysics Data System (ADS)

Martinez, Luis; Meneveau, Charles

2014-11-01

Large Eddy Simulations (LES) of the flow past a single wind turbine with uniform inflow have been performed. A goal of the simulations is to compare two turbulence subgrid-scale models and their effects in predicting the initial breakdown, transition and evolution of the wake behind the turbine. Prior works have often observed negligible sensitivities to subgrid-scale models. The flow is modeled using an in-house LES with pseudo-spectral discretization in horizontal planes and centered finite differencing in the vertical direction. Turbines are represented using the actuator line model. We compare the standard constant-coefficient Smagorinsky subgrid-scale model with the Lagrangian Scale Dependent Dynamic model (LSDM). The LSDM model predicts faster transition to turbulence in the wake, whereas the standard Smagorinsky model predicts significantly delayed transition. The specified Smagorinsky coefficient is larger than the dynamic one on average, increasing diffusion thus delaying transition. A second goal is to compare the resulting near-blade properties such as local aerodynamic forces from the LES with Blade Element Momentum Theory. Results will also be compared with those of the SOWFA package, the wind energy CFD framework from NREL. This work is supported by NSF (IGERT and IIA-1243482) and computations use XSEDE resources, and has benefitted from interactions with Dr. M. Churchfield of NREL.

New Insights into the Morphology of the Galapagos Platform from Lower Crustal Flow Models

NASA Astrophysics Data System (ADS)

Orellana Rovirosa, F.; Richards, M. A.

2014-12-01

The volcanically-active Galapagos Islands are constructed upon a broad platform, with the westernmost islands being the most active, marking the present-day position of the Galapagos hotspot (mantle plume). This volcanic platform overlies relatively young oceanic lithosphere (<15 Myr) and exhibits unique morphologic features along its boundaries. The most spectacular of these features is a system of stepped terraces on the southwestern escarpment, with very large vertical relief (>3 km), and contrasting with relatively gentle slopes off the eastern platform edge toward the Carnegie ridge. Considering the horizontal lithostatic pressure differences associated with this bathymetric relief, along with the high temperatures within this young, hotspot-affected oceanic lithosphere and crust; it is likely that lower crustal flow contributes significantly, perhaps even dominantly, to lithospheric and crustal deformation within the Galapagos Platform. Using 2D numerical models that invoke a thin-sheet approximation for the Stokes' equation for a Newtonian fluid with space- and time-dependent viscosity, and assuming isostatic conditions, we show that: (1) the pronounced bathymetric rim along the Eastern platform region (where gravimetric studies indicate Airy isostasy) near Española Island may be the expression of a mature stage of a lower crustal flow front evolving asymptotically during the last ~3 Myr; (2) the spectacular system of stepped terraces along the southwestern edge of the platform may be explained by lower crustal flow-associated backward tilting of the bathymetric surface that, even with small amounts of angular change (~0.1 deg) and potentially occurring in non-isostatic regimes, effectively hinders the horizontal flow of lava. This process of backward-tilting may have been largely restricted to the last ~1 Myr of platform growth, and hence may be a unique event that may involve horizontal extrusion of large lava sheets from within the southwestward-marching lower-crustal flow front.

Fluid pressure and flow at great depth in the continental crust. A discussion in relation to topography, temperature and salinity distribution using as an example the KTB Fault Zones in connection with the Eger Rift Hot Spot.

NASA Astrophysics Data System (ADS)

Kessels, W.; Kuhlmann, S.; Li, X.

2006-12-01

Hydraulic investigations in and between the two KTB boreholes have shown that groundwater flow is possible at great depth in the crystalline crust. Remarkable permeability was found particularly in the SE1 and SE2 fault zones. The results from a long term pump and injection test, and the related three-dimensional groundwater modelling (Graesle et al., 2006), document the existence of a large-scale (more than 10 km) hydraulic reservoir in the crystalline crust. According to this calculation, an overpressure of 0.4 MPa can be still be expected in KTB-HB in 2009, 4 years after the end of the injection. The good match with the measurement data confirms groundwater pathways at a scale of more than 10 km. The isotopic water composition recovered from the KTB pilot hole indicates a downward water flow along the SE2 fault zone, which is in contact with the Franconian Line. Moreover, there is a deep upward groundwater flow 60 km away in the western Eger Rift Valley as indicated e.g. by the temperature signature and gas flow observations. Therefore, the demand for fluid mass continuity means that water is being supplied by a downstream groundwater flow, probably from the Franconian Line. The question of potential driving processes must be answered to understand and quantify the flow in the deeper crust at a scale of 10 km to 100 km. The processes must result in a sufficient horizontal pressure gradient to allow groundwater flow at great depth. The density variations of groundwater with depth are highly relevant for the calculation of horizontal pressure differences. The two independent potential fields of gravity and pressure have to be considered. Differentiation into 4 relevant driving processes is required: \\bullet The groundwater surface topography related to the groundwater recharge and mean regional distance between neighbouring valleys \\bullet Geothermal gradient and water density depending on temperature and pressure \\bullet Different salt contents in adjacent geological formations \\bullet Gas content in the water and gas dissolution The interpretation of these processes for the Eger Rift Franconian Line area results in horizontal pressure gradients up to 0.5 MPa/km. With these pressure gradients in deep fault zones similar to the KTB fault zones SE1 and SE2, a remarkable groundwater flow is also possible in the deep crystalline crust. For only a 1 MPa pressure difference between the Franconian Line and the Eger Rift Valley, which lie nearly 60 km apart, we get a tracer velocity of 1.0 to 5.0 m/a (using the Darcy relation and porosities for the hydraulic KTB data). The flow system at great depth is determined mainly by the counteractive forces of salinity and temperature with a nonlinear relation to the water density. References GRAESLE, W., KESSELS, W., KUEMPEL, H.-J., LI, XUAN (2006): HYDRAULIC OBSERVATIONS FROM A ONE YEAR FLUID PRODUCTION TEST IN THE 4000 M DEEP KTB PILOT BOREHOLE. GEOFLUIDS, 6, 8 23 KESSELS, W., KUECK, J. (1995): HYDRAULIC COMMUNICATION IN CRYSTALLINE ROCK BETWEEN THE TWO BOREHOLES OF THE CONTINENTAL DEEP DRILLING PROJECT IN GERMANY. INT. J. ROCK MECH. MIN. SCI. &GEOMECH. ABSTR., 32, 37 47

Irreversible transport in the stratosphere by internal waves of short vertical wavelength

NASA Technical Reports Server (NTRS)

Danielsen, Edwin F.; Hipskind, R. S.; Starr, Walter L.; Vedder, James F.; Gaines, Steven E.; Kley, Dieter; Kelley, Ken K.

1991-01-01

Measurements performed during stratospheric flights of the U-2 aircraft confirm that cross-jet transport is dominated by waves, not by large-scale circulations. Monotonic gradients of trace constituents normal to the jet axis, with upper stratospheric tracers increasing poleward and tropospheric tracers increasing equatorward, are augmented by large-scale confluence as the jet intensifies during cyclogenesis. These gradients are rotated, intensified, and significantly increased in areas as their mixing ratio surfaces are folded by the differential transport of a very low frequency transverse wave. The quasi-horizontal transport produces a laminar structure with stable layers rich in upper stratospheric tracers alternating vertically with less stable layers rich in tropospheric tracers. The transport proceeds toward irreversibility at higher frequency, shear-gravity waves extend the folding to smaller horizontal scales.

A Comparative Study of the Electrical Structure of Circum Tibetan Plateau Orogenic Belts and its Tectonic Implications

NASA Astrophysics Data System (ADS)

Jin, Sheng; Zhang, Letian; Wei, Wenbo; Ye, Gaofeng; Jing, Jianen; Dong, Hao; Xie, Chengliang; Yin, Yaotian

2017-04-01

The Tibetan Plateau, as known as "roof of the world", was created through the on-going continent-continent collision between the Indian and Eurasian plates since 55 Ma. As the process continues, the plateau is growing both vertically and horizontally. The horizontal expansion of the plateau is blocked by the Yangtze block in the east, the Tarim block in the north, and the Ordos block in the northeast, and consequently lead to the formation of the circum Tibetan plateau orogenic belts. To better understand the mechanism behind this process, we conducted a comparative study by collecting 7 magnetotelluric (MT) profiles over the margins of the Tibetan plateau, namely, the INDEPTH 100, 700 and 800 lines in the southern Tibet, the INDEPTH 4000 and 5000 lines across the Altyn Tagh fault on the northern margin of the plateau, as well as other two profiles across the Haiyuan fault and the Longmenshan fault on the northeastern and eastern margins of the plateau deployed under the framework of project SinoProbe. The electrical features of the stable blocks surrounding the Tibetan plateau are generally resistive, while crustal conductive layers are found to be wide spread within the plateau. The southern margin of the Tibetan plateau is characterized by large scale underthrust of the Indian lithosphere beneath the plateau. This intense converging process created the thrust fault system distributed along the southern margin of the Tibetan plateau over 1000 km. Crustal conductive layers discovered in southern Tibet are generally associated with the southward crustal flow that originated from the lower crust within the plateau and exhumed along the thrust belts in the Himalayas. On the eastern margin of the Tibetan plateau, the electrical structures suggest that the Yangtze block wedged into the Tibetan lithosphere and caused decoupling between the crust and upper mantel. Large scale conductors discovered beneath the Songpan-Ganze block reflect that the eastward crustal flow was blocked and piled up along the eastern margin of the plateau due to the block of the Sichuan Basin, which further result in the uplift and expansion of the eastern Tibetan plateau. The northeastern and northern margins of the Tibetan plateau is bounded by large scale left-lateral strike-slip Haiyuan and Altyn Tagh faults. In these regions, the plateau interacts with the surrounding stable blocks in a way of oblique strike-slip. The deformation of the northern Tibetan lithosphere is dominated by crustal thickening, where no features of decoupling or large scale underthrusting were seen. Crustal conductors in these regions are generally not very well connected, which suggest the absence of crustal flow. Deep metamorphism fluids could be an alternative interpretation of the crustal conductors in these regions. * This work was jointly supported by the grants from Project SinoProbe-02-04 and National Natural Science Foundation of China (41404060).

Turbulence scalings in pipe flows exhibiting polymer-induced drag reduction

NASA Astrophysics Data System (ADS)

Zadrazil, Ivan; Markides, Christos

2014-11-01

Non-intrusive laser based diagnostics technique, namely Particle Image Velocimetry, was used to in detail characterise polymer induced drag reduction in a turbulent pipe flow. The effect of polymer additives was investigated in a pneumatically-driven flow facility featuring a horizontal pipe test section of inner diameter 25.3 mm and length 8 m. Three high molecular weight polymers (2, 4 and 8 MDa) at concentrations of 5 - 250 wppm were used at Reynolds numbers from 35000 to 210000. The PIV derived results show that the level of drag reduction scales with different normalised turbulence parameters, e.g. streamwise and spanwise velocity fluctuations, vorticity or Reynolds stresses. These scalings are dependent of the distance from the wall, however, are independent of the Reynolds numbers range investigated.

Comparisons between observations and numerical simulations of Japan (East) Sea flow and mass fields in 1999 through 2001

NASA Astrophysics Data System (ADS)

Mooers, Christopher N. K.; Bang, Inkweon; Sandoval, Francisco J.

2005-06-01

The Princeton Ocean Model (POM), as implemented for the Japan (East) Sea (JES) with mesoscale-admitting resolution is driven by seasonal throughflow and synoptic atmospheric forcing for 1999 through 2001. Temperature and salinity profiles from shipborne and PALACE float CTDs, and horizontal velocities at 800 m from PALACE float trajectories, plus horizontal velocities at 15 m from WOCE surface drifters for 1988 through 2001, are used to assess the performance of the numerical simulations for a base case. General agreement exists in the circulation at 15 and 800 m and the horizontal and vertical structure of the upper ocean temperature and salinity fields. The mean observed flow at 15 m defines the two branches of the Tsushima Warm Current and hints at the existence of a large cyclonic gyre over the Japan Basin, which the simulations also produce. The mean observed flow at 800 m defines a large cyclonic recirculation gyre over the Japan Basin that validates the simulated flow pattern. Variances of the observed and simulated flows at 15 and 800 m have similar patterns. The main discrepancies are associated with the strength of the seasonal thermocline and halocline and the location of the Subpolar Front. When smoother topography and smaller lateral friction are used in other cases, the thermocline and halocline strengthen, agreeing better with the observed values, and when 80% of total outflow transport is forced to exit through Soya Strait, the Subpolar Front extends along the coast to the north of Tsugaru Strait, which is an observed feature absent in the base case.

Observations of wave-induced pore pressure gradients and bed level response on a surf zone sandbar

NASA Astrophysics Data System (ADS)

Anderson, Dylan; Cox, Dan; Mieras, Ryan; Puleo, Jack A.; Hsu, Tian-Jian

2017-06-01

Horizontal and vertical pressure gradients may be important physical mechanisms contributing to onshore sediment transport beneath steep, near-breaking waves in the surf zone. A barred beach was constructed in a large-scale laboratory wave flume with a fixed profile containing a mobile sediment layer on the crest of the sandbar. Horizontal and vertical pore pressure gradients were obtained by finite differences of measurements from an array of pressure transducers buried within the upper several centimeters of the bed. Colocated observations of erosion depth were made during asymmetric wave trials with wave heights between 0.10 and 0.98 m, consistently resulting in onshore sheet flow sediment transport. The pore pressure gradient vector within the bed exhibited temporal rotations during each wave cycle, directed predominantly upward under the trough and then rapidly rotating onshore and downward as the wavefront passed. The magnitude of the pore pressure gradient during each phase of rotation was correlated with local wave steepness and relative depth. Momentary bed failures as deep as 20 grain diameters were coincident with sharp increases in the onshore-directed pore pressure gradients, but occurred at horizontal pressure gradients less than theoretical critical values for initiation of the motion for compact beds. An expression combining the effects of both horizontal and vertical pore pressure gradients with bed shear stress and soil stability is used to determine that failure of the bed is initiated at nonnegligible values of both forces.