Solution of the Boussinesq equations by means of the finite element method

A finite element method is presented for the computation of flows that are influenced by buoyancy forces. The accuracy of several finite elements is studied by solving the Bénard problem and determining the critical Rayleigh number. It is found that the accuracy is greatly enhanced if the shape functions satisfy a certain requirement that arises from the physical nature of the problem

An asymptotic expansion for product integration applied to Cauchy principal value integrals

Product integration methods for Cauchy principal value integrals based on piecewise Lagrangian interpolation are studied. It is shown that for this class of quadrature methods the truncation error has an asymptotic expansion in integer powers of the step-size, and that a method with an asymptotic expansion in even powers of the step-size does not exist. The relative merits of a quadrature method which employs values of both the integrand and its first derivative and for which the truncation error has an asymptotic expansion in even powers of the step-size are discussed

A Two-Variable Asymptotic Solution for Three - Dimensional, Solar-Powered, Low-Thrust Trajectories in the Vicinity of the Ecliptic Plane

Two-variable asymptotic solution for three-dimensional solar powered low thrust trajectories in vicinity of ecliptic plan

Modeling of droplet breakup in a microfluidic T--shaped junction with a phase--field model

A phase--field method is applied to the modeling of flow and breakup of droplets in a T--shaped junction in the hydrodynamic regime where capillary and viscous stresses dominate over inertial forces, which is characteristic of microfluidic devices. The transport equations are solved numerically in the three--dimensional geometry, and the dependence of the droplet breakup on the flow rates, surface tension and viscosities of the two components is investigated in detail. The model reproduces quite accurately the phase diagram observed in experiments performed with immiscible fluids. The critical capillary number for droplet breakup depends on the viscosity contrast, with a trend which is analogous to that observed for free isolated droplets in hyperbolic flow

Multigrid optimization for space-time discontinuous Galerkin discretizations of advection dominated flows

The goal of this research is to optimize multigrid methods for higher order accurate space-time discontinuous Galerkin discretizations. The main analysis tool is discrete Fourier analysis of two- and three-level multigrid algorithms. This gives the spectral radius of the error transformation operator which predicts the asymptotic rate of convergence of the multigrid algorithm. In the optimization process we therefore choose to minimize the spectral radius of the error transformation operator. We specifically consider optimizing h-multigrid methods with explicit Runge-Kutta type smoothers for second and third order accurate space-time discontinuous Galerkin finite element discretizations of the 2D advection-diffusion equation. The optimized schemes are compared with current h-multigrid techniques employing Runge-Kutta type smoothers. Also, the efficiency of h-, p- and hp-multigrid methods for solving the Euler equations of gas dynamics with a higher order accurate space-time DG method is investigated

A new wireless underground network system for continuous monitoring of soil water contents

A new stand-alone wireless embedded network system has been developed recently for continuous monitoring of soil water contents at multiple depths. This paper presents information on the technical aspects of the system, including the applied sensor technology, the wireless communication protocols, the gateway station for data collection, and data transfer to an end user Web page for disseminating results to targeted audiences. Results from the first test of the network system are presented and discussed, including lessons learned so far and actions to be undertaken in the near future to improve and enhance the operability of this innovative measurement approac

Geodynamo and mantle convection simulations on the Earth Simulator using the Yin-Yang grid

We have developed finite difference codes based on the Yin-Yang grid for the geodynamo simulation and the mantle convection simulation. The Yin-Yang grid is a kind of spherical overset grid that is composed of two identical component grids. The intrinsic simplicity of the mesh configuration of the Yin-Yang grid enables us to develop highly optimized simulation codes on massively parallel supercomputers. The Yin-Yang geodynamo code has achieved 15.2 Tflops with 4096 processors on the Earth Simulator. This represents 46% of the theoretical peak performance. The Yin-Yang mantle code has enabled us to carry out mantle convection simulations in realistic regimes with a Rayleigh number of $10^7$ including strongly temperature-dependent viscosity with spatial contrast up to $10^6$.Comment: Plenary talk at SciDAC 200

Spontaneous Branching of Anode-Directed Streamers between Planar Electrodes

Non-ionized media subject to strong fields can become locally ionized by penetration of finger-shaped streamers. We study negative streamers between planar electrodes in a simple deterministic continuum approximation. We observe that for sufficiently large fields, the streamer tip can split. This happens close to Firsov's limit of `ideal conductivity'. Qualitatively the tip splitting is due to a Laplacian instability quite like in viscous fingering. For future quantitative analytical progress, our stability analysis of planar fronts identifies the screening length as a regularization mechanism.Comment: 4 pages, 6 figures, submitted to PRL on Nov. 16, 2001, revised version of March 10, 200

Nuclear structure studies with the 7Li(e,e'p) reaction

Experimental momentum distributions for the transitions to the ground state and first excited state of 6He have been measured via the reaction 7Li(e,e'p)6He, in the missing momentum range from -70 to 260 MeV/c. They are compared to theoretical distributions calculated with mean-field wave functions and with variational Monte Carlo (VMC) wave functions which include strong state-dependent correlations in both 7Li and 6He. These VMC calculations provide a parameter-free prediction of the momentum distribution that reproduces the measured data, including its normalization. The deduced summed spectroscopic factor for the two transitions is 0.58 +/- 0.05, in perfect agreement with the VMC value of 0.60. This is the first successful comparison of experiment and ab initio theory for spectroscopic factors in p-shell nuclei.Comment: 4 pages, 3 figure