Adaptive Finite Element Methods with Inexact Solvers for the Nonlinear Poisson-Boltzmann Equation

In this article we study adaptive finite element methods (AFEM) with inexact solvers for a class of semilinear elliptic interface problems. We are particularly interested in nonlinear problems with discontinuous diffusion coefficients, such as the nonlinear Poisson-Boltzmann equation and its regularizations. The algorithm we study consists of the standard SOLVE-ESTIMATE-MARK-REFINE procedure common to many adaptive finite element algorithms, but where the SOLVE step involves only a full solve on the coarsest level, and the remaining levels involve only single Newton updates to the previous approximate solution. We summarize a recently developed AFEM convergence theory for inexact solvers, and present a sequence of numerical experiments that give evidence that the theory does in fact predict the contraction properties of AFEM with inexact solvers. The various routines used are all designed to maintain a linear-time computational complexity.Comment: Submitted to DD20 Proceeding

A tunable radiation source by coupling laser-plasma-generated electrons to a periodic structure

Near-infrared radiation around 1000 nm generated from the interaction of a high-density MeV electron beam, obtained by impinging an intense ultrashort laser pulse on a solid target, with a metal grating is observed experimentally. Theoretical modeling and particle-in-cell simulation suggest that the radiation is caused by the Smith-Purcell mechanism. The results here indicate that tunable terahertz radiation with tens GV=m field strength can be achieved by using appropriate grating parameter

New ventilation cooling for modular PM machines utilizing flux gaps and rotor ducts

This paper investigates a new ventilation cooling technology for modular permanent magnet (PM) machines. Air is pumped into the stator flux gaps and rotor ducts. A significant temperature reduction for the modular machines can be observed by introducing ventilated air. The flux gaps within the modular machines provide extra coolant path, notably modify the flow profile. Moreover, the flux gaps also provide additional contact area between machines and coolant, improving heat transfer rate. Comparative study between different machine configurations, including different inlet and outlet areas, different topologies of shaft, and different stator topologies have been investigated. Their influences on machine cooling have been investigated using computational fluid dynamics (CFD) modelling. Furthermore, since the coolant (air) will pass through the stator-rotor annular gap, i.e., airgap, as well as rotor ducts, the rotor speed will influence the cooling efficiency of the ventilation cooling system. This influence comes from shock loss, friction loss, and combining flow phenomenon. A series of experiments have been conducted to validate the CFD simulations

Immersed boundary-finite element model of fluid-structure interaction in the aortic root

It has long been recognized that aortic root elasticity helps to ensure efficient aortic valve closure, but our understanding of the functional importance of the elasticity and geometry of the aortic root continues to evolve as increasingly detailed in vivo imaging data become available. Herein, we describe fluid-structure interaction models of the aortic root, including the aortic valve leaflets, the sinuses of Valsalva, the aortic annulus, and the sinotubular junction, that employ a version of Peskin's immersed boundary (IB) method with a finite element (FE) description of the structural elasticity. We develop both an idealized model of the root with three-fold symmetry of the aortic sinuses and valve leaflets, and a more realistic model that accounts for the differences in the sizes of the left, right, and noncoronary sinuses and corresponding valve cusps. As in earlier work, we use fiber-based models of the valve leaflets, but this study extends earlier IB models of the aortic root by employing incompressible hyperelastic models of the mechanics of the sinuses and ascending aorta using a constitutive law fit to experimental data from human aortic root tissue. In vivo pressure loading is accounted for by a backwards displacement method that determines the unloaded configurations of the root models. Our models yield realistic cardiac output at physiological pressures, with low transvalvular pressure differences during forward flow, minimal regurgitation during valve closure, and realistic pressure loads when the valve is closed during diastole. Further, results from high-resolution computations demonstrate that IB models of the aortic valve are able to produce essentially grid-converged dynamics at practical grid spacings for the high-Reynolds number flows of the aortic root

Search for the Rare Decays J/Psi --> Ds- e+ nu_e, J/Psi --> D- e+ nu_e, and J/Psi --> D0bar e+ e-

We report on a search for the decays J/Psi --> Ds- e+ nu_e + c.c., J/Psi --> D- e+ nu_e + c.c., and J/Psi --> D0bar e+ e- + c.c. in a sample of 5.8 * 10^7 J/Psi events collected with the BESII detector at the BEPC. No excess of signal above background is observed, and 90% confidence level upper limits on the branching fractions are set: B(J/Psi --> Ds- e+ nu_e + c.c.)<4.8*10^-5, B(J/Psi --> D- e+ nu_e + c.c.) D0bar e+ e- + c.c.)<1.1*10^-5Comment: 10 pages, 4 figure

Measurements of the observed cross sections for exclusive light hadron production in e^+e^- annihilation at \sqrt{s}= 3.773 and 3.650 GeV

By analyzing the data sets of 17.3 pb$^{-1}$ taken at $\sqrt{s}=3.773$ GeV and 6.5 pb$^{-1}$ taken at $\sqrt{s}=3.650$ GeV with the BESII detector at the BEPC collider, we have measured the observed cross sections for 12 exclusive light hadron final states produced in $e^+e^-$ annihilation at the two energy points. We have also set the upper limits on the observed cross sections and the branching fractions for $\psi(3770)$ decay to these final states at 90% C.L.Comment: 8 pages, 5 figur

Measurements of the observed cross sections for e+e−→e^+e^-\to exclusive light hadrons containing π0π0\pi^0\pi^0 at s=3.773\sqrt s= 3.773, 3.650 and 3.6648 GeV

By analyzing the data sets of 17.3, 6.5 and 1.0 pb$^{-1}$ taken, respectively, at $\sqrt s= 3.773$, 3.650 and 3.6648 GeV with the BES-II detector at the BEPC collider, we measure the observed cross sections for $e^+e^-\to \pi^+\pi^-\pi^0\pi^0$, $K^+K^-\pi^0\pi^0$, $2(\pi^+\pi^-\pi^0)$, $K^+K^-\pi^+\pi^-\pi^0\pi^0$ and $3(\pi^+\pi^-)\pi^0\pi^0$ at the three energy points. Based on these cross sections we set the upper limits on the observed cross sections and the branching fractions for $\psi(3770)$ decay into these final states at 90% C.L..Comment: 7 pages, 2 figure

Partial wave analysis of J/\psi \to \gamma \phi \phi

Using $5.8 \times 10^7 J/\psi$ events collected in the BESII detector, the radiative decay $J/\psi \to \gamma \phi \phi \to \gamma K^+ K^- K^0_S K^0_L$ is studied. The $\phi\phi$ invariant mass distribution exhibits a near-threshold enhancement that peaks around 2.24 GeV/$c^{2}$. A partial wave analysis shows that the structure is dominated by a $0^{-+}$ state ($\eta(2225)$) with a mass of $2.24^{+0.03}_{-0.02}{}^{+0.03}_{-0.02}$ GeV/$c^{2}$ and a width of $0.19 \pm 0.03^{+0.06}_{-0.04}$ GeV/$c^{2}$. The product branching fraction is: $Br(J/\psi \to \gamma \eta(2225))\cdot Br(\eta(2225)\to \phi\phi) = (4.4 \pm 0.4 \pm 0.8)\times 10^{-4}$.Comment: 11 pages, 4 figures. corrected proof for journa