Optimal velocity functions for car-following models

The integral part of the optimal velocity car-following models is the optimal velocity function (OVF), which can be derived from measured velocity-spacing data. This paper discusses several characteristics of the OVF and presents regression analysis on two classical datasets, the Lincoln and Holland tunnels, with different possible OVFs. The numerical simulation of the formation of traffic congestion is conducted with three different heuristic OVFs, demonstrating that these functions give results similar to those of the famous Bando OVF (Bando et al., 1995). Also an alternative method is present for determining the sensitivity and model parameters based on a single car driving to a fixed barrier.

Left ventricular regional and global diastolic function assessed using Quantitative Tissue velocity Imaging in patients with hypertrophic cardiomyopathy

Objectives The study was performed to assess the left ventricular (LV) regional and global diastolic function、left ventricular wall motion features in patients with Hypertrophic cardiomyopathy by Quantitative Tissue Velocity Imaging (QTVI). Methods 42 patients with hypertrophic cardiomyopathy and 36 age-matched normal subjects underwent QTVI study. Off-line LV regional muscular tissue velocity Imaging along LV apical long-axis view were obtained. Regional diastolic function was assessed in using peak tissue velocities of LV regional muscular tissue during early diastole (Ve)and LA contraction (Va), Ve/Va ratio, derived from Tissue Velocity Imaging. Global diastolic function was reflected by isovolumic relaxation time(IRT) and mitral valve peak flow velocity ( E/A ) calculated with pulsed wave doppler. The end-diastolic interventricular septal thickness (ⅣSt) was measured by conventional 2 - dimension echocardiography. Results ① Ve、 Va、 Ve/Va in the segments of hypertrophic interventricular septum (IVS) reduced wlhile E/A ratio significantly reduced and IRT markedly prolonged in HCM patients than in normal subjects。 ② Ve、 Ve/Va were significant reduced in the segments of hypertrophic interventricular septum compared with other LV segments in HCM patients . ③ There was a correlation between Ve/Va and E/A in HCM patients with abnormal E/A ratio (r = 0. 70). ④ There was a negative correlation between Ve/Va and ⅣSt in non -obstruction HCM patients (B group , r = -0.61 ) Conclusions QTVI offers a newer method in clinical practice which has a higher sensibility and accuracy in evaluating the LV regional and global diastolic function in HCM patients .

MULTI-SCALE COHERENT STRUCTURES IN TURBULENT BOUNDARY LAYER DETECTED BY LOCALLY AVERAGED VELOCITY STRUCTURE FUNCTIONS

The time sequence of longitudinal velocity component at different vertical locations in turbulent boundary layer was finely measured in a wind tunnel. The concept of coarse_grained velocity structure functions, which describes the relative motions of straining and compressing for multi_scale eddy structures in turbulent flows, was put forward based on the theory of locally multi_scale average. Based on the consistency between coarse_grained velocity structure function and Harr wavelet transformation,detecting method was presented, by which the coherent structures and their intermittency was identified by multi_scale flatness factor calculated by locally average structure function. Phase_averaged evolution course for multi_scale coherent eddy structures in wall turbulence were extracted by this conditional sampling to educe scheme. The dynamics course of multi_scale coherent eddy structures and their effects on statistics of turbulent flows were studied.

A multifractal model for linking Lagrangian and Eulerian velocity structure functions

A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian time scale, and it should not be the Eulerian time scale which was widely used in previous studies on Lagrangian statistics. Using the present model, we can obtain the scaling exponents of Lagrangian velocity structure functions from the existing data or models of scaling exponents of Eulerian velocity structure functions. This model is validated by comparing its prediction with the results of experiments, direct numerical simulations, and the previous theoretical models. The comparison shows that the proposed model can better predict the scaling exponents of Lagrangian velocity structure functions, especially for orders larger than 6.

Detailed string stability analysis for bi-directional optimal velocity model

The class of bi-directional optimal velocity models can describe the bi-directional looking effect that usually exists in the reality and is even enhanced with the development of the connected vehicle technologies. Its combined string stability condition can be obtained through the method of the ring-road based string stability analysis. However, the partial string stability about traffic fluctuation propagated backward or forward was neglected, which will be analyzed in detail in this work by the method of transfer function and its H∞ norm from the viewpoint of control theory. Then, through comparing the conditions of combined and partial string stabilities, their relationships can make traffic flow be divided into three distinguishable regions, displaying various combined and partial string stability performance. Finally, the numerical experiments verify the theoretical results and find that the final displaying string stability or instability performance results from the accumulated and offset effects of traffic fluctuations propagated from different directions.

A control method applied to mixed traffic flow for the coupled-map car-following model

In light of previous work [Phys. Rev. E 60 4000 （1999）], a modified coupled-map car-following model is proposed by considering the headways of two successive vehicles in front of a considered vehicle described by the optimal velocity function. The non-jam conditions are given on the basis of control theory. Through simulation, we find that our model can exhibit a better effect as p = 0.65, which is a parameter in the optimal velocity function. The control scheme, which was proposed by Zhao and Gao, is introduced into the modified model and the feedback gain range is determined. In addition, a modified control method is applied to a mixed traffic system that consists of two types of vehicle. The range of gains is also obtained by theoretical analysis. Comparisons between our method and that of Zhao and Gao are carried out, and the corresponding numerical simulation results demonstrate that the temporal behavior of traffic flow obtained using our method is better than that proposed by Zhao and Gao in mixed traffic systems.

Modified coupled map car-following model and its delayed feedback control scheme

A modified coupled map car-following model is proposed, in which two successive vehicle headways in front of the considering vehicle is incorporated into the optimal velocity function. The steady state under certain conditions is obtained. An error system around the steady state is studied further. Moreover, the condition for the state having no traffic jam is derived. A new control scheme is presented to suppress the traffic jam in the modified coupled map car-following model under the open boundary. A control signal including the velocity differences between the following and the considering vehicles, and between the preceding and the considering vehicles is used. The condition under which the traffic jam can be well suppressed is analysed. The results are compared with that presented by t^onishi et al. （the KKH model）. The simulation results show that the temporal behaviour obtained in our model is better than that in the KKH model. The simulation results are in good agreement with the theoretical analysis.

Evaluation of Phonon-level Density of UO2 by Molecular Dynamics Simulation

Molecular dynamics calculation of UO2 in a wide temperature range are presented and dis- cussed. The calculated lattice parameters, mean square displacements, and dynamic prop- erty of phonon-level density of the velocity auto-correlation functions for UO2 are provided. The Morelon potential and the Basak potential are employed. It confirms that the cal- culated lattice parameters using the Basak potential are in nearly perfect agreement with the reported values. The models successfully predict mean square displacement and Bredig transition. Furthermore, the phonon-level density of uranium dioxide are discussed. The intensity of phonon-level density increases with temperature, and the properties of UO2 are characterized by large thermal vibrations rather than extensive disorder.

The far-field plasma characterization in a 600W Hall thruster plume by laser-induced fluorescence

Non-intrusive characterization of the singly ionized xenon velocity in Hall thruster plume using laser induced fluorescence(LIF)is critical for constructing a complete picture of plume plasma,deeply understanding the ion dynamics in the plume,and providing validation data for numerical simulation.This work presents LIF measurements of singly ionized xenon axial velocity on a grid ranging from 100 to 300 mm in axial direction and from 0 to 50 mm in radial direction for a600 W Hall thruster operating at the nominal condition of discharge voltage 300 V and discharge current 2 A,the influence of discharge voltage is investigated as well.The ion velocity distribution function(IVDF)results in the far-field plume demonstrate a profile of bimodal IVDFs,especially prominent at radial distances greater than channel inner radius of 22 mm at axial position of 100 mm,which is quite different from that of the near-field plume where bimodal IVDFs occur in the central core region for the same power Hall thruster when compared to previous LIF measurements of BHT-600 by Hargus(2010 J.Propulsion Power 26135).Beyond 100 mm,only single-peak IVDFs are measured.The two-dimensional ion velocity vector field indicates the bimodal axial IVDF is merely a geometry effect for the annular discharge channel in the far-field plume.Results about the IVDF,the most probable velocity and the accelerating potential profile along the centerline all indicate that ions are still accelerating at axial distances greater than 100 mm,and the maximum most probable velocity measured at300 mm downstream of the exit plane is about 19 km s-1.In addition,the most probable velocity of ions along radial direction changes a little except the lower velocity ion populations in the bimodal IVDF cases.The ion temperature at axial distances of 10 and 300 mm oscillates along the radial direction,while the ion temperature first increases,and then decreases for the 200 mm case.Finally,the axial position for the ion peak axial velocity on the thruster centerline is shifted upstream for higher discharge voltages,and the velocity curve is becoming steeper with the discharge voltage before reaching the maximum.This observation can be used as a criterion to optimize the thruster performance.

Influence of Cardiac Structure, Blood Blow, Velocity and Heart Function on Circulating Atrial Natriuretic Peptide and Reninangiotension System

Heart function and plasma atrial natriuretic peptide (ANP).plasma renin actiity(PRA) andangiotensionⅡ(Ang Ⅱ) were examined with echocardiography and radioimmunoassay in patients with dilated cardiomyopathy(DCM),

Quantum twisted double-slits experiments:confirming wavefunctions'physical reality

Are quantum states real? This most fundamental question in quantum mechanics has not yet been satisfactorily resolved, although its realistic interpretation seems to have been rejected by various delayedchoice experiments. Here, to address this long-standing issue, we present a quantum twisted double-slit experiment. By exploiting the subluminal feature of twisted photons, the real nature of a photon during its time in flight is revealed for the first time. We found that photons' arrival times were inconsistent with the states obtained in measurements but agreed with the states during propagation. Our results demonstrate that wavefunctions describe the realistic existence and evolution of quantum entities rather than a pure mathematical abstraction providing a probability list of measurement outcomes. This finding clarifies the long-held misunderstanding of the role of wavefunctions and their collapse in the evolution of quantum entities.

Evolution of vortex structure around a wall-mounted rough hemisphere

In this paper,the large eddy simulation(LES)method is used to calculate Reynolds numbers the range in the range of 3000–30000(based on the diameter of the hemisphere D),and then the statistical analysis method is used to study the influence of Re on the characteristics of the flow field around the rough hemisphere and the evolution law of the tail vortex.The effect of the Reynolds number on the structure of the rough hemisphere flow field and the evolution of vortex structure is studied by vorticity distribution,POD method,and correlation analysis.It can be observed that with the increase of Reynolds number,the shear in the recirculating region of the rough hemisphere is enhanced,the large-scale structure is reduced,the small-scale structure is increased,and the angle between the hairpin vortex and the wall is reduced.λ_(ci)vortex identification and the Q-criterion are used to capture the evolution of a periodic vortex structure.Q-criterion gives the three-dimensional vortex structure behind the hemispherical flow.It is found that the horseshoe vortex and arch vortex generated by the incoming flow after passing through the hemisphere will interact with the wall surface to form a hairpin vortex,and the hairpin vortex structure will continue to evolve downstream to form a new vortex structure.

Convergence proof of the DSMC method and the Gas-Kinetic Unified Algorithm for the Boltzmann equation

This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.

Numerical study on the gas-kinetic high-order schemes for solving Boltzmann model equation

The high-order compact finite difference technique is introduced to solve the Boltzmann model equation, and the gas-kinetic high-order schemes are developed to simulate the different kinetic model equations such as the BGK model, the Shakhov model and the Ellipsoidal Statistical (ES) model in this paper. The methods are tested for the one-dimensional unsteady shock-tube problems with various Knudsen numbers, the inner flows of normal shock wave for different Mach numbers, and the two-dimensional flows past a circular cylinder and a NACA 002 airfoil to verify the reliability of the present high-order algorithm and simulate gas transport phenomena covering various flow regimes. The computed results are found in good agreement both with the theoretical prediction from continuum to rarefied gas dynamics, the related DSMC solutions, and with the experimental results. The numerical effect of the schemes with the different precision and the different types of Boltzmann collision models on the computational efficiency and computed results is investigated and analyzed. The numerical experience indicates that an approach developing and applying the gas-kinetic high-order algorithm is feasible for directly solving the Boltzmann model equation.

A High-Order Accurate Gas-Kinetic Scheme for One-and Two-Dimensional Flow Simulation

This paper develops a high-order accurate gas-kinetic scheme in the framework of the finite volume method for the one-and two-dimensional flow simulations,which is an extension of the third-order accurate gas-kinetic scheme[Q.B.Li,K.Xu,and S.Fu,J.Comput.Phys.,229(2010),6715-6731]and the second-order accurate gas-kinetic scheme[K.Xu,J.Comput.Phys.,171(2001),289-335].It is formed by two parts:quartic polynomial reconstruction of the macroscopic variables and fourth-order accurate flux evolution.The first part reconstructs a piecewise cell-center based quartic polynomial and a cell-vertex based quartic polynomial according to the“initial”cell average approximation of macroscopic variables to recover locally the non-equilibrium and equilibrium single particle velocity distribution functions around the cell interface.It is in view of the fact that all macroscopic variables become moments of a single particle velocity distribution function in the gas-kinetic theory.The generalized moment limiter is employed there to suppress the possible numerical oscillation.In the second part,the macroscopic flux at the cell interface is evolved in fourth-order accuracy by means of the simple particle transport mechanism in the microscopic level,i.e.free transport and the Bhatnagar-Gross-Krook(BGK)collisions.In other words,the fourth-order flux evolution is based on the solution(i.e.the particle velocity distribution function)of the BGK model for the Boltzmann equation.Several 1D and 2D test problems are numerically solved by using the proposed high-order accurate gas-kinetic scheme.By comparing with the exact solutions or the numerical solutions obtained the secondorder or third-order accurate gas-kinetic scheme,the computations demonstrate that our scheme is effective and accurate for simulating invisid and viscous fluid flows,and the accuracy of the high-order GKS depends on the choice of the(numerical)collision time.

Analysis of the Microphysical Properties of a Stratiform Rain Event Using an L-Band Profiler Radar

This paper investigates spatial and temporal distributions of the microphysical properties of precipitating stratiform clouds based on Doppler spectra of rain particles observed by an L-band profiler radar.The retrieval of raindrop size distributions（RSDs） is accomplished through eliminating vertical air motion and isolating the terminal fall velocity of raindrops in the observed Doppler velocity spectrum.The microphysical properties of raindrops in a broad stratiform region with weak convective cells are studied using data collected from a 1320-MHz wind profiler radar in Huayin,Shaanxi Province on 14 May 2009.RSDs and gamma function parameters are retrieved at altitudes between 700 and 3000 m above the surface,below a melting layer.It is found that the altitude of the maximum number of raindrops was closely related to the surface rain rate.The maximum number of large drops was observed at lower altitudes earlier in the precipitation event but at higher altitudes in later periods,suggesting decreases in the numbers of large and medium size raindrops.These decreases may have been caused by the breakup of larger drops and evaporation of smaller drops as they fell.The number of medium size drops decreased with increasing altitude.The relationship between reflectivity and liquid water content during this precipitation event was Z = 1.69×10~4M^（1.5）,and the relationship between reflectivity and rain intensity was Z = 256I^（1.4）.