Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence

A scale-similarity model for Lagrangian two-point, two-time velocity correlations LVCs in isotropic turbulence is developed from the Kolmogorov similarity hypothesis. It is a second approximation to the isocontours of LVCs, while the Smith-Hay model is only a first approximation. This model expresses the LVC by its space correlation and a dispersion velocity. We derive the analytical expression for the dispersion velocity from the Navier-Stokes equations using the quasinormality assumption. The dispersion velocity is dependent on enstrophy spectra and shown to be smaller than the sweeping velocity for the Eulerian velocity correlation. Therefore, the Lagrangian decorrelation process is slower than the Eulerian decorrelation process. The data from direct numerical simulation of isotropic turbulence support the scale-similarity model: the LVCs for different space separations collapse into a universal form when plotted against the separation axis defined by the model.

Measurement of velocity profiles in a rectangular microchannel with aspect ratio α = 0.35

In this work, we measured 14 horizontal velocity profiles along the vertical direction of a rectangular microchannel with aspect ratio alpha = h/w = 0.35 (h is the height of the channel and w is the width of the channel) using microPIV at Re = 1.8 and 3.6. The experimental velocity profiles are compared with the full 3D theoretical solution, and also with a Poiseuille parabolic profile. It is shown that the experimental velocity profiles in the horizontal and vertical planes are in agreement with the theoretical profiles, except for the planes close to the wall. The discrepancies between the experimental data and 3D theoretical results in the center vertical plane are less than 3.6%. But the deviations between experimental data and Poiseuille's results approaches 5%. It indicates that 2D Poiseuille profile is no longer a perfect theoretical approximation since a = 0.35. The experiments also reveal that, very near the hydrophilic wall (z = 0.5-1 mu m), the measured velocities are significantly larger than the theoretical velocity based on the no-slip assumption. A proper discussion on some physical effects influencing the near wall velocity measurement is given.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

The near wall velocity measurements in microchannels with different diameter particles by microPIV/PTV. 8th International symposium on PIV

In near wall measurements with microPIV/PTV, whether seeding particles can be effectively used to detect local fluid velocity is a crucial problem. This talk presents our recent measurements in microchannels [1][2]. Based on measured velocity profiles with 200nm and 50nm in pure water, we found that the measured velocity profiles are agreed with the theoretical values in the middle of channel, but large deviations between measured data and theoretical prediction appear close to wall (0.25mm velocity deviation appeared in shear flows.

Group velocity of probe light pulse in an open lambda system

The group velocity of the probe light pulse (GVPLP) propagating through an open Lambda-type atomic system with a spontaneously generated coherence is investigated when the weak probe and strong driving light fields have different frequencies. It is found that adjusting the detuning or Rabi frequency of the probe light field can realize switching of the GVPLP from subluminal to superluminal. Changing the relative phase between the probe and driving light. elds or atomic exit and injection rates can lead to GVPLP varying in a wider range, but cannot induce transformation of the property of the GVPLP. The absolute value of the GVPLP always increases with Rabi frequency of the driving light field increasing. For subluminal and superluminal propagation, the system always exhibits the probe absorption, and GVPLP is mainly determined by the slope of the steep dispersion.

Controllable group velocity of the probe light in a Lambda-type system with two fold levels

The group velocities of the probe laser field are studied in a A-type system where one lower state has two fold levels coupled by a control field. It is found that the interaction of double dark states leads to controllable group velocity of the probe field in this system. It can be easily realized, due to the interacting double dark resonances, that one of the group velocities at transparency positions is much slower than the other by tuning the control field to be off resonance. In particular, when the control field is on resonance. we can obtain two equal slow group velocities with a broader EIT width, which provides potential applications in quantum storage and retrieval of light. (c) 2005 Elsevier B.V. All rights reserved.

Influence of Coulomb potential for photoionization of H atoms in an elliptically polarized laser field: Velocity gauge versus length gauge

We theoretically study the influence of Coulomb potential for photoionization of hydrogen atoms in an intense laser field with elliptical polarization. The total ionization rates, photoelectron energy spectra, and photoelectron angular distributions are calculated with the Coulomb-Volkov wave functions in the velocity gauge and compared with those calculated in the length gauge as well as those calculated with the Volkov wave functions. By comparing the results obtained by the Coulomb-Volkov and Volkov wave functions, we find that for linear polarization the influence of Coulomb potential is obvious for low-energy photoelectrons, and as the photoelectron energy and/or the laser intensity increase, its influence becomes smaller. This trend, however, is not so clear for the case of elliptical polarization. We also find that the twofold symmetry in the photoelectron angular distributions for elliptical polarization is caused by the cooperation of Coulomb potential and interference of multiple transition channels. About the gauge issue, we show that the difference in the photoelectron angular distributions obtained by the velocity and length gauges becomes rather obvious for elliptical polarization, while the difference is generally smaller for linear polarization.

Teleseismic array analysis of upper mantle compressional velocity structure

Large quantities of teleseismic short-period seismograms recorded at SCARLET provide travel time, apparent velocity and waveform data for study of upper mantle compressional velocity structure. Relative array analysis of arrival times from distant (30° < Δ < 95°) earthquakes at all azimuths constrains lateral velocity variations beneath southern California. We compare dT/dΔ back azimuth and averaged arrival time estimates from the entire network for 154 events to the same parameters derived from small subsets of SCARLET. Patterns of mislocation vectors for over 100 overlapping subarrays delimit the spatial extent of an east-west striking, high-velocity anomaly beneath the Transverse Ranges. Thin lens analysis of the averaged arrival time differences, called 'net delay' data, requires the mean depth of the corresponding lens to be more than 100 km. Our results are consistent with the PKP-delay times of Hadley and Kanamori (1977), who first proposed the high-velocity feature, but we place the anomalous material at substantially greater depths than their 40-100 km estimate.

Detailed analysis of travel time, ray parameter and waveform data from 29 events occurring in the distance range 9° to 40° reveals the upper mantle structure beneath an oceanic ridge to depths of over 900 km. More than 1400 digital seismograms from earthquakes in Mexico and Central America yield 1753 travel times and 58 dT/dΔ measurements as well as high-quality, stable waveforms for investigation of the deep structure of the Gulf of California. The result of a travel time inversion with the tau method (Bessonova et al., 1976) is adjusted to fit the p(Δ) data, then further refined by incorporation of relative amplitude information through synthetic seismogram modeling. The application of a modified wave field continuation method (Clayton and McMechan, 1981) to the data with the final model confirms that GCA is consistent with the entire data set and also provides an estimate of the data resolution in velocity-depth space. We discover that the upper mantle under this spreading center has anomalously slow velocities to depths of 350 km, and place new constraints on the shape of the 660 km discontinuity.

Seismograms from 22 earthquakes along the northeast Pacific rim recorded in southern California form the data set for a comparative investigation of the upper mantle beneath the Cascade Ranges-Juan de Fuca region, an ocean-continent transit ion. These data consist of 853 seismograms (6° < Δ < 42°) which produce 1068 travel times and 40 ray parameter estimates. We use the spreading center model initially in synthetic seismogram modeling, and perturb GCA until the Cascade Ranges data are matched. Wave field continuation of both data sets with a common reference model confirms that real differences exist between the two suites of seismograms, implying lateral variation in the upper mantle. The ocean-continent transition model, CJF, features velocities from 200 and 350 km that are intermediate between GCA and T7 (Burdick and Helmberger, 1978), a model for the inland western United States. Models of continental shield regions (e.g., King and Calcagnile, 1976) have higher velocities in this depth range, but all four model types are similar below 400 km. This variation in rate of velocity increase with tectonic regime suggests an inverse relationship between velocity gradient and lithospheric age above 400 km depth.

Part I. A study of the velocity structure of the earth by the use of core phases. Part II. The 1971 San Fernando earthquake series focal mechanisms and tectonics

The initial objective of Part I was to determine the nature of upper mantle discontinuities, the average velocities through the mantle, and differences between mantle structure under continents and oceans by the use of P'dP', the seismic core phase P'P' (PKPPKP) that reflects at depth d in the mantle. In order to accomplish this, it was found necessary to also investigate core phases themselves and their inferences on core structure. P'dP' at both single stations and at the LASA array in Montana indicates that the following zones are candidates for discontinuities with varying degrees of confidence: 800-950 km, weak; 630-670 km, strongest; 500-600 km, strong but interpretation in doubt; 350-415 km, fair; 280-300 km, strong, varying in depth; 100-200 km, strong, varying in depth, may be the bottom of the low-velocity zone. It is estimated that a single station cannot easily discriminate between asymmetric P'P' and P'dP' for lead times of about 30 sec from the main P'P' phase, but the LASA array reduces this uncertainty range to less than 10 sec. The problems of scatter of P'P' main-phase times, mainly due to asymmetric P'P', incorrect identification of the branch, and lack of the proper velocity structure at the velocity point, are avoided and the analysis shows that one-way travel of P waves through oceanic mantle is delayed by 0.65 to 0.95 sec relative to United States mid-continental mantle.

A new P-wave velocity core model is constructed from observed times, dt/dΔ's, and relative amplitudes of P'; the observed times of SKS, SKKS, and PKiKP; and a new mantle-velocity determination by Jordan and Anderson. The new core model is smooth except for a discontinuity at the inner-core boundary determined to be at a radius of 1215 km. Short-period amplitude data do not require the inner core Q to be significantly lower than that of the outer core. Several lines of evidence show that most, if not all, of the arrivals preceding the DF branch of P' at distances shorter than 143° are due to scattering as proposed by Haddon and not due to spherically symmetric discontinuities just above the inner core as previously believed. Calculation of the travel-time distribution of scattered phases and comparison with published data show that the strongest scattering takes place at or near the core-mantle boundary close to the seismic station.

In Part II, the largest events in the San Fernando earthquake series, initiated by the main shock at 14 00 41.8 GMT on February 9, 1971, were chosen for analysis from the first three months of activity, 87 events in all. The initial rupture location coincides with the lower, northernmost edge of the main north-dipping thrust fault and the aftershock distribution. The best focal mechanism fit to the main shock P-wave first motions constrains the fault plane parameters to: strike, N 67° (± 6°) W; dip, 52° (± 3°) NE; rake, 72° (67°-95°) left lateral. Focal mechanisms of the aftershocks clearly outline a downstep of the western edge of the main thrust fault surface along a northeast-trending flexure. Faulting on this downstep is left-lateral strike-slip and dominates the strain release of the aftershock series, which indicates that the downstep limited the main event rupture on the west. The main thrust fault surface dips at about 35° to the northeast at shallow depths and probably steepens to 50° below a depth of 8 km. This steep dip at depth is a characteristic of other thrust faults in the Transverse Ranges and indicates the presence at depth of laterally-varying vertical forces that are probably due to buckling or overriding that causes some upward redirection of a dominant north-south horizontal compression. Two sets of events exhibit normal dip-slip motion with shallow hypocenters and correlate with areas of ground subsidence deduced from gravity data. Several lines of evidence indicate that a horizontal compressional stress in a north or north-northwest direction was added to the stresses in the aftershock area 12 days after the main shock. After this change, events were contained in bursts along the downstep and sequencing within the bursts provides evidence for an earthquake-triggering phenomenon that propagates with speeds of 5 to 15 km/day. Seismicity before the San Fernando series and the mapped structure of the area suggest that the downstep of the main fault surface is not a localized discontinuity but is part of a zone of weakness extending from Point Dume, near Malibu, to Palmdale on the San Andreas fault. This zone is interpreted as a decoupling boundary between crustal blocks that permits them to deform separately in the prevalent crustal-shortening mode of the Transverse Ranges region.

Study of noncollinear quasi-phase matching optical parametric amplification with group-velocity matching based on periodically poled crystals

The properties of noncollinear optical parametric amplification (NOPA) based on quasi-phase matching of periodically poled crystals are investigated, under the condition that the group velocity matching (GVM) of the signal and idler pulses is satisfied. Our study focuses on the dependence of the gain spectrum upon the noncollinear angle, crystal temperature, and crystal angle with periodically poled KTiOPO4 (PPKTP), periodically poled LiNbO3 (PPLN), and periodically poled LiTaO3 (PPLT), and the NOPA gain properties of the three crystals are compared. Broad gain bandwidth exists above 85 nm at a signal wavelength of 800 nm with a 532 nm pump pulse, with proper noncollinear angle and grating period at a fixed temperature for GVM. Deviation from the group-velocity-matched noncollinear angle can be compensated by accurately tuning the crystal angle or temperature with a fixed grating period for phase matching. Moreover, there is a large capability of crystal angle tuning.

Group velocity of the light pulse in an open V-type system

We investigate the group velocity of the probe light pulse in an open V-type system with spontaneously generated coherence. We find that, not only varying the relative phase between the probe and driving pulses can but varying the atomic exit rate or incoherent pumping rate also can manipulate dramatically the group velocity, even make the pulse propagation switching from subluminal to superluminal; the subliminal propagation can be companied with gain or absorption, but the superluminal propagation is always companied with absorption. (c) 2006 Elsevier GmbH. All rights reserved.

Velocity resolved - scalar modeled simulations of high Schmidt number turbulent transport

The objective of this thesis is to develop a framework to conduct velocity resolved - scalar modeled (VR-SM) simulations, which will enable accurate simulations at higher Reynolds and Schmidt (Sc) numbers than are currently feasible. The framework established will serve as a first step to enable future simulation studies for practical applications. To achieve this goal, in-depth analyses of the physical, numerical, and modeling aspects related to Sc>>1 are presented, specifically when modeling in the viscous-convective subrange. Transport characteristics are scrutinized by examining scalar-velocity Fourier mode interactions in Direct Numerical Simulation (DNS) datasets and suggest that scalar modes in the viscous-convective subrange do not directly affect large-scale transport for high Sc. Further observations confirm that discretization errors inherent in numerical schemes can be sufficiently large to wipe out any meaningful contribution from subfilter models. This provides strong incentive to develop more effective numerical schemes to support high Sc simulations. To lower numerical dissipation while maintaining physically and mathematically appropriate scalar bounds during the convection step, a novel method of enforcing bounds is formulated, specifically for use with cubic Hermite polynomials. Boundedness of the scalar being transported is effected by applying derivative limiting techniques, and physically plausible single sub-cell extrema are allowed to exist to help minimize numerical dissipation. The proposed bounding algorithm results in significant performance gain in DNS of turbulent mixing layers and of homogeneous isotropic turbulence. Next, the combined physical/mathematical behavior of the subfilter scalar-flux vector is analyzed in homogeneous isotropic turbulence, by examining vector orientation in the strain-rate eigenframe. The results indicate no discernible dependence on the modeled scalar field, and lead to the identification of the tensor-diffusivity model as a good representation of the subfilter flux. Velocity resolved - scalar modeled simulations of homogeneous isotropic turbulence are conducted to confirm the behavior theorized in these a priori analyses, and suggest that the tensor-diffusivity model is ideal for use in the viscous-convective subrange. Simulations of a turbulent mixing layer are also discussed, with the partial objective of analyzing Schmidt number dependence of a variety of scalar statistics. Large-scale statistics are confirmed to be relatively independent of the Schmidt number for Sc>>1, which is explained by the dominance of subfilter dissipation over resolved molecular dissipation in the simulations. Overall, the VR-SM framework presented is quite effective in predicting large-scale transport characteristics of high Schmidt number scalars, however, it is determined that prediction of subfilter quantities would entail additional modeling intended specifically for this purpose. The VR-SM simulations presented in this thesis provide us with the opportunity to overlap with experimental studies, while at the same time creating an assortment of baseline datasets for future validation of LES models, thereby satisfying the objectives outlined for this work.