A Lagrangian stochastic model for the trajectories of particle pairs and its application to the prediction of concentration variance within plant canopies

A simple Lagrangian stochastic model for the trajectories of particle pairs in high Reynolds-number turbulent flows is presented. In this model, the velocities of particle pairs are initially correlated but subsequently each particle moves independently. The independent single-particle trajectories are simulated using Thomson's model [J. Fluid Mech. 180, 529–556, 1987]. This two-particle model exactly satisfies the well-mixed condition for Gaussian turbulence when length scales, characterizing the two-point Eulerian velocity correlation function, vanish. Temperature variances, due to heat released as a passive scalar from an elevated plane source, within a model plant canopy (Coppin et al. Boundary Layer Meteorol. 35, 167–191, 1986) are shown to be well predicted by the model. It is suggested that for strongly inhomogeneous flows, the two-point Eulerian velocity function is of secondary importance in determining the simulated trajectories of particle pairs compared to the importance of ensuring satisfaction of the two-to-one constraint (Borgas and Sawford. J. Fluid Mech. 279, 69–99, 1994); i.e ensuring that one-particle statistics obtained from the two-particle model are the same as those obtained from the corresponding one-particle model. Limitations of this modelling approach are discussed.     

Mode identification of terrestrial ULF waves observed by Cluster: A case study

Using multipoint measurements from the Cluster mission wave identification techniques are applied to observations of ULF waves made in the terrestrial foreshock with the aim of identifying the modes and properties of the waves taking into account the effects of a high beta plasma. The wave properties in the spacecraft and plasma rest frames are experimentally derived using minimum variance analysis. Two waves with periods of 30 and 3 s dominate the dynamic frequency spectrum. The results indicate that these waves propagate in the fast magnetosonic and Alfvén/Ion Cyclotron modes, respectively. Both waves propagate in the upstream direction in the plasma rest frame but are convected downstream in the spacecraft frame. The measured wave properties in the plasma rest frame are in good agreement with those obtained from the theoretical kinetic dispersion relation taking into account the effects of different plasma beta. The dispersion results show a rather significant deviation from fluid model, especially when high beta plasma conditions occur. These experimentally derived foreshock ULF wave properties are in good agreement with previous results but when the effects of a high beta plasma are considered it is not as straight forward to choose the correct wave mode branch.     

近50a东北冷涡异常特征及其与淮河流域降水的关系

利用美国国家环境预报/大气研究中心(NCEP/NCAR)再分析月平均数据分析了1960—2009 年夏季东北冷涡的异常特征,研究了夏季东北冷涡与同期淮河流域降水的关系,发现:东北冷涡偏强时,淮河流域的降水很可能偏多,东北冷涡偏弱时,淮河流域的降水很可能偏少。东北冷涡异常强年,淮河流域高低层环流具有斜压性,且低层有显著的正涡度发展,促进了上升对流运动活跃发展。而东北冷涡活动异常频繁,有利于引导潮湿阴冷的北方气流南下,与东亚夏季盛行的西南暖湿气流在淮河流域上空交汇,在上升运动的触发下,导致淮河流域降水明显增多;东北冷涡弱年的情况正好相反。     

Tropical Cyclone Energy Dispersion and Possible Self Maintenance Mechanism of the Synoptic Scale Wave Trains

The interaction between tropical cyclone (TC) and the large-scale mean flows such as the inter-tropical convergence zone (ITCZ) is investigated using a three-dimensional primitive equation model. Once a TC develops in the vicinity of the ITCZ region where satisfies both barotropic and baroclinic instabilities, the southeastward energy dispersion from the TC may disturb the ITCZ and thus help its breakdown. Cumulus convection can be organized in the region of cyclonic circulation, and the interaction between convective heating and the perturbation circulation may enhance the development of the waves, leading to the generation of a new tropical cyclone to the east. While the TC moves to the high latitude, the ITCZ will reform. Though repeating of this process, a synoptic-scale wave train oriented in the northwest-southeast direction can be generated and self-maintained. The results suggest that the mutual interaction among the low-frequency background flow, wave train pattern and TCs provides a possible mechanism for the origin of the summer synoptic scale wave train pattern over the western North Pacific.     

Observed winds,turbulence, and dispersion in built-up downtown areas of Oklahoma City and Manhattan

Wind and tracer data from the Oklahoma City Joint Urban 2003 (JU2003) and the Manhattan Madison Square Garden 2005 (MSG05) urban field experiments are being analyzed to aid in understanding air flow and dispersion near street-level in built-up downtown areas. The mean winds are separately calculated for groups of anemometers having similar exposures such as “near street level” and “on building top”. Several general results are found, such as the scalar wind speed at street level is about 1/3 of that at building top. Turbulent standard deviations of wind speed components and temperature, and vertical fluxes of momentum and sensible heat, are calculated from sonic anemometers near street level at 20 locations in JU2003 and five locations in MSG05, and from two rooftop locations in MSG05. The turbulence observations are consistent with observations in the literature at other cities, although the JU2003 and MSG05 data are unique in that many data are available near street level. For example, it is found that the local (i.e., at the measuring height) averages about 1.5 and the local averages about 0.25 in the two cities, where is the standard deviation of vertical velocity fluctuations, is the friction velocity, and u is the wind speed. The ratio of temperature fluctuations to temperature scale, , averages about −3 in both cities, consistent with similarity theory for slightly unstable conditions, where is the standard deviation of temperature fluctuations, and is the temperature scale. The calculated Obukhov length, L, is also consistent with slightly unstable conditions near street level, even at night during JU2003. The SF₆ tracer concentration observations from JU2003 are analyzed. Values of for the continuous releases are calculated for each release and arc distance, where is the 30-min average arc maximum concentration, Q is the continuous source emission rate, and u is the spatial-averaged wind speed in the downtown area. The basic characteristics of the JU2003 plot of averaged agree reasonably well with similar plots for other urban experiments in Salt Lake City and London (i.e., at . A is found to be about 3 during the day and about 10 during the night.     

Multi-resolution adaptive modeling of groundwater flow and transport problems

Many groundwater flow and transport problems, especially those with sharp fronts, narrow transition zones, layers and fingers, require extensive computational resources. In this paper, we present a novel multi-resolution adaptive Fup approach to solve the above mentioned problems. Our numerical procedure is the Adaptive Fup Collocation Method (AFCM), based on Fup basis functions and designed through a method of lines (MOL). Fup basis functions are localized and infinitely differentiable functions with compact support and are related to more standard choices such as splines or wavelets. This method enables the adaptive multi-resolution approach to solve problems with different spatial and temporal scales with a desired level of accuracy using the entire family of Fup basis functions. In addition, the utilized collocation algorithm enables the mesh free approach with consistent velocity approximation and flux continuity due to properties of the Fup basis functions. The introduced numerical procedure was tested and verified by a few characteristic groundwater flow and transport problems, the Buckley–Leverett multiphase flow problem, the 1-D vertical density driven problem and the standard 2-D seawater intrusion benchmark–Henry problem. The results demonstrate that the method is robust and efficient particularly when describing sharp fronts and narrow transition zones changing in space and time.     

http://www.sciencedirect.com/science/article/pii/S1674987110000058

<正>The thermo-electric coefficients of twenty-six magnetite samples,formed either by magmatism or metamorphism,were tested by the thermo-electric instrument BHET—06.Results showed that the coefficient is of a constant value of about -0.05 mV/℃.It is emphasized that because every magnetite grain was tested randomly,the coefficient is independent of the crystallographic direction.This fact means the thermal voltage generated from a single magnetite crystal can be accumulated,and as a result a new thermo-electric field can arise when a gradient thermal field exists and is active within the earth's crust.Because magnetite is widespread in the earth's crust(generally appearing more in the middle-lower crust),there is more-than-random probability that the additional thermo-electric field can be generated when certain thermal conditions are fulfilled.We,therefore,used the thermo-electric effect of magnetite to study the mechanism responsible for the presence of abnormal geo-electric fields during earthquake formation and occurrence, because gradient thermal fields always exist before earthquakes.The possible presence of additional thermo-electric fields was calculated under theoretical seismological conditions,using the following calcu-lation formula:E= - 0.159(σ×△T×φ×ρ_2×[(h~2-2x~2)cosα+ 3hxsinα]/ρ_1(h~2 +x~2)~(5/2)).In the above formula,σis thermo-electric coefficient of magnetite,△T is the temperature difference acting on it,φis a sectional area on a block of magnetite vertically perpendicular to the direction of the thermal current.ρ_1 andρ_2 are the respective resistivities of magnetite and the crust,and h,α,and x,respectively,h is the depth of embedded magnetite block,αmeans the angle created by the horizontal line and ligature of the two poles of magnetite block,and x is the distance from observation point to projective center point of the magnetite block on earth surface.According to simulations calculated with this formula,additional thermo-electric field intensity may reach as high as n to n×10~2 mV/km.This field is strong enough to cause obvious anomalies in the background geo-electric field,and can be easy probed by earthquake monitoring equipment. Therefore,we hypothesize that geo-electric abnormalities which occur during earthquakes may be caused by the thermo-electric effect of magnetite.     

A fast inversion method for ocean parameters based on dispersion curves with a single hydrophone

The dispersion characteristics of shallow water can be described by the dispersion curves, which contain substantial ocean parameter information. A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper. The method is achieved through Bayesian theory. Several sets of dispersion curves extracted from measured data are used as the input function. The inversion is performed by matching a replica calculated with a dispersion formula. The bottom characteristics can be described by the bottom reflection phase shift parameter P. The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P. The inversion results improve the inversion efficiency of the seabed parameters. Consequently, the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased. The inversion results have lower error than the reference values, and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data; thus, the effectiveness of the inversion method is demonstrated.     

Distribution of stars perpendicular to the plane of the galaxy

We present here rigorous analytical solutions for the Boltzmann-Poisson equation concerning the distribution of stars above the galactic plane. The number density of stars is considered to follow a behaviour n(m,0) ∼H(m - m₀)m^−x, wherem is the mass of a star andx an arbitrary exponent greater than 2 and also the velocity dispersion of the stars is assumed to behave as < v²(m)> ∼ m^−θ the exponent θ being arbitrary and positive. It is shown that an analytic expression can be found for the gravitational field K_z, in terms of confluent hypergeometric functions, the limiting trends being K_z∼z for z →0, while K_z → constant for z → infinity. We also study the behaviour of < |z(m)|²>,i.e. the dispersion of the distance from the galactic disc for the stars of massm. It is seen that the quantity < |z(m)|²>∼ m^t-θ, for m→ t, while it departs significantly from this harmonic oscillator behaviour for stars of lighter masses. It is suggested that observation of < |z(m)|²> can be used as a probe to findx and hence obtain information about the mass spectrum.