Numerical simulation of laminar separated flows on adaptive tri-tree grids with the finite volume method

In this work, a code has been developed that solves the Navier-Stokes equations using the finite volume method with unstructured triangular grids. A cell-centred, finite volume method is used and the pressure-velocity coupling is treated using both the SMTLE and the MAC algorithms. The major advantage of using triangular grids is their applicability to complex geometry. A special treatment is developed to ensure good quality triangular elements around the boundaries. The numerical simulation of incompressible flow at low Reynolds number is studied in this thesis. A code for generating triangular grids using the tri-tree algorithm has been written and an adaptive finite volume method developed for calculating laminar fluid flow. The grid is locally adapted at each time step, with grid refinement and derefinement dependent on the vorticity magnitude. The resulting grids have fine local resolution and are economical in reducing the numerical simulation time. The discretised equations are solved by using an iterative point by point Gauss-Seidel solver. For calculating the values of velocity and pressure at vertices of triangular grids, special interpolation schemes (averaged linear-interpolation and scattered interpolation) are used to increase the accuracy. To avoid the well known checkerboard error problems, i. e., the oscillations occurring in the pressure field, third derivative terms in pressure, first introduced by Rhie-chow (1983), are added to the mass flux velocity. Convective terms are approximated using a QUICK (Quadratic Upstream Interpolation for Convective Kinematics) differencing scheme which has been developed here in for unstructured grids. Three cases of two-dimensional viscous incompressible fluid flow have been investigated: the first is channel flow, in which the numerical results are compared with the analytical solution; the second case is the backward-facing step flow; and the third case is flow past circular cylinders at low Reynolds number (Re). The numerical results obtained for the last two cases are compared with published data. The evolution of vortex shedding is presented for the case of unidirectional flow past a circular cylinder at Re=200. In addition, drag and lift force coefficients are calculated and compared for single and multiple cylinders in unidirectional flow

Observation of Two New N* Peaks in J/psi -> ppi−nˉp pi^- \bar n and pˉπ+n\bar p\pi^+n Decays

The $\pi N$ system in decays of $J/\psi\to\bar NN\pi$ is limited to be isospin 1/2 by isospin conservation. This provides a big advantage in studying $N^*\to \pi N$ compared with $\pi N$ and $\gamma N$ experiments which mix isospin 1/2 and 3/2 for the $\pi N$ system. Using 58 million $J/\psi$ decays collected with the Beijing Electron Positron Collider, more than 100 thousand $J/\psi \to p \pi^- \bar n + c.c.$ events are obtained. Besides two well known $N^*$ peaks at 1500 MeV and 1670 MeV, there are two new, clear $N^*$ peaks in the $p\pi$ invariant mass spectrum around 1360 MeV and 2030 MeV. They are the first direct observation of the $N^*(1440)$ peak and a long-sought "missing" $N^*$ peak above 2 GeV in the $\pi N$ invariant mass spectrum. A simple Breit-Wigner fit gives the mass and width for the $N^*(1440)$ peak as $1358\pm 6 \pm 16$ MeV and $179\pm 26\pm 50$ MeV, and for the new $N^*$ peak above 2 GeV as $2068\pm 3^{+15}_{-40}$ MeV and $165\pm 14\pm 40$ MeV, respectively

Partial Wave Analysis of J/ψ→γ(K+K−π+π−)J/\psi \to \gamma (K^+K^-\pi^+\pi^-)

BES data on $J/\psi \to \gamma (K^+K^-\pi^+\pi^-)$ are presented. The $K^*\bar K^*$ contribution peaks strongly near threshold. It is fitted with a broad $0^{-+}$ resonance with mass $M = 1800 \pm 100$ MeV, width $\Gamma = 500 \pm 200$ MeV. A broad $2^{++}$ resonance peaking at 2020 MeV is also required with width $\sim 500$ MeV. There is further evidence for a $2^{-+}$ component peaking at 2.55 GeV. The non-$K^*\bar K^*$ contribution is close to phase space; it peaks at 2.6 GeV and is very different from $K^{*}\bar{K^{*}}$.Comment: 15 pages, 6 figures, 1 table, Submitted to PL

Direct Measurements of the Branching Fractions for D0→K−e+νeD^0 \to K^-e^+\nu_e and D0→π−e+νeD^0 \to \pi^-e^+\nu_e and Determinations of the Form Factors f+K(0)f_{+}^{K}(0) and f+π(0)f^{\pi}_{+}(0)

The absolute branching fractions for the decays $D^0 \to K^-e ^+\nu_e$ and $D^0 \to \pi^-e^+\nu_e$ are determined using $7584\pm 198 \pm 341$ singly tagged $\bar D^0$ sample from the data collected around 3.773 GeV with the BES-II detector at the BEPC. In the system recoiling against the singly tagged $\bar D^0$ meson, $104.0\pm 10.9$ events for $D^0 \to K^-e ^+\nu_e$ and $9.0 \pm 3.6$ events for $D^0 \to \pi^-e^+\nu_e$ decays are observed. Those yield the absolute branching fractions to be $BF(D^0 \to K^-e^+\nu_e)=(3.82 \pm 0.40\pm 0.27)$ and $BF(D^0 \to \pi^-e^+\nu_e)=(0.33 \pm 0.13\pm 0.03)$. The vector form factors are determined to be $|f^K_+(0)| = 0.78 \pm 0.04 \pm 0.03$ and $|f^{\pi}_+(0)| = 0.73 \pm 0.14 \pm 0.06$. The ratio of the two form factors is measured to be $|f^{\pi}_+(0)/f^K_+(0)|= 0.93 \pm 0.19 \pm 0.07$.Comment: 6 pages, 5 figure

Study of J/psi decays to Lambda Lambdabar and Sigma0 Sigma0bar

The branching ratios and Angular distributions for J/psi decays to Lambda Lambdabar and Sigma0 Sigma0bar are measured using BESII 58 million J/psi.Comment: 11 pages, 5 figure

Measurements of J/psi Decays into 2(pi+pi-)eta and 3(pi+pi-)eta

Based on a sample of 5.8X 10^7 J/psi events taken with the BESII detector, the branching fractions of J/psi--> 2(pi+pi-)eta and J/psi-->3(pi+pi-)eta are measured for the first time to be (2.26+-0.08+-0.27)X10^{-3} and (7.24+-0.96+-1.11)X10^{-4}, respectively.Comment: 11 pages, 6 figure

BESII Detector Simulation

A Monte Carlo program based on Geant3 has been developed for BESII detector simulation. The organization of the program is outlined, and the digitization procedure for simulating the response of various sub-detectors is described. Comparisons with data show that the performance of the program is generally satisfactory.Comment: 17 pages, 14 figures, uses elsart.cls, to be submitted to NIM

Measurement of branching fractions for the inclusive Cabibbo-favored ~K*0(892) and Cabibbo-suppressed K*0(892) decays of neutral and charged D mesons

The branching fractions for the inclusive Cabibbo-favored ~K*0 and Cabibbo-suppressed K*0 decays of D mesons are measured based on a data sample of 33 pb-1 collected at and around the center-of-mass energy of 3.773 GeV with the BES-II detector at the BEPC collider. The branching fractions for the decays D+(0) -> ~K*0(892)X and D0 -> K*0(892)X are determined to be BF(D0 -> \~K*0X) = (8.7 +/- 4.0 +/- 1.2)%, BF(D+ -> ~K*0X) = (23.2 +/- 4.5 +/- 3.0)% and BF(D0 -> K*0X) = (2.8 +/- 1.2 +/- 0.4)%. An upper limit on the branching fraction at 90% C.L. for the decay D+ -> K*0(892)X is set to be BF(D+ -> K*0X) < 6.6%

Measurements of the Mass and Full-Width of the ηc\eta_c Meson

In a sample of 58 million $J/\psi$ events collected with the BES II detector, the process J/$\psi\to\gamma\eta_c$ is observed in five different decay channels: $\gamma K^+K^-\pi^+\pi^-$, $\gamma\pi^+\pi^-\pi^+\pi^-$, $\gamma K^\pm K^0_S \pi^\mp$ (with $K^0_S\to\pi^+\pi^-$), $\gamma \phi\phi$ (with $\phi\to K^+K^-$) and $\gamma p\bar{p}$. From a combined fit of all five channels, we determine the mass and full-width of $\eta_c$ to be $m_{\eta_c}=2977.5\pm1.0 ({stat.})\pm1.2 ({syst.})$ MeV/$c^2$ and $\Gamma_{\eta_c} = 17.0\pm3.7 ({stat.})\pm7.4 ({syst.})$ MeV/$c^2$.Comment: 9 pages, 2 figures and 4 table. Submitted to Phys. Lett.