Domain decomposition algorithms and computation fluid dynamics

In the past several years, domain decomposition was a very popular topic, partly motivated by the potential of parallelization. While a large body of theory and algorithms were developed for model elliptic problems, they are only recently starting to be tested on realistic applications. The application of some of these methods to two model problems in computational fluid dynamics are investigated. Some examples are two dimensional convection-diffusion problems and the incompressible driven cavity flow problem. The construction and analysis of efficient preconditioners for the interface operator to be used in the iterative solution of the interface solution is described. For the convection-diffusion problems, the effect of the convection term and its discretization on the performance of some of the preconditioners is discussed. For the driven cavity problem, the effectiveness of a class of boundary probe preconditioners is discussed

The physics of parallel machines

The idea is considered that architectures for massively parallel computers must be designed to go beyond supporting a particular class of algorithms to supporting the underlying physical processes being modelled. Physical processes modelled by partial differential equations (PDEs) are discussed. Also discussed is the idea that an efficient architecture must go beyond nearest neighbor mesh interconnections and support global and hierarchical communications

ENO-wavelet transforms for piecewise smooth functions

We have designed an adaptive essentially nonoscillatory (ENO)-wavelet transform for approximating discontinuous functions without oscillations near the discontinuities. Our approach is to apply the main idea from ENO schemes for numerical shock capturing to standard wavelet transforms. The crucial point is that the wavelet coefficients are computed without differencing function values across jumps. However, we accomplish this in a different way than in the standard ENO schemes. Whereas in the standard ENO schemes the stencils are adaptively chosen, in the ENO-wavelet transforms we adaptively change the function and use the same uniform stencils. The ENO-wavelet transform retains the essential properties and advantages of standard wavelet transforms such as concentrating the energy to the low frequencies, obtaining maximum accuracy, maintained up to the discontinuities, and having a multiresolution framework and fast algorithms, all without any edge artifacts. We have obtained a rigorous approximation error bound which shows that the error in the ENO-wavelet approximation depends only on the size of the derivative of the function away from the discontinuities. We will show some numerical examples to illustrate this error estimate

Heralded magnetism in non-Hermitian atomic systems

Quantum phase transitions are usually studied in terms of Hermitian Hamiltonians. However, cold-atom experiments are intrinsically non-Hermitian due to spontaneous decay. Here, we show that non-Hermitian systems exhibit quantum phase transitions that are beyond the paradigm of Hermitian physics. We consider the non-Hermitian XY model, which can be implemented using three-level atoms with spontaneous decay. We exactly solve the model in one dimension and show that there is a quantum phase transition from short-range order to quasi-long-range order despite the absence of a continuous symmetry in the Hamiltonian. The ordered phase has a frustrated spin pattern. The critical exponent $\nu$ can be 1 or 1/2. Our results can be seen experimentally with trapped ions, cavity QED, and atoms in optical lattices.Comment: 7 pages + appendi

Arc-Length Continuation and Multigrid Techniques for Nonlinear Elliptic Eigenvalue Problems

We investigate multi-grid methods for solving linear systems arising from arc-length continuation techniques applied to nonlinear elliptic eigenvalue problems. We find that the usual multi-grid methods diverge in the neighborhood of singular points of the solution branches. As a result, the continuation method is unable to continue past a limit point in the Bratu problem. This divergence is analyzed and a modified multi-grid algorithm has been devised based on this analysis. In principle, this new multi-grid algorithm converges for elliptic systems, arbitrarily close to singularity and has been used successfully in conjunction with arc-length continuation procedures on the model problem. In the worst situation, both the storage and the computational work are only about a factor of two more than the unmodified multi-grid methods

Exercise induced collapse: insulin shock

This issue of eMedRef provides information to clinicians on the pathophysiology, diagnosis, and therapeutics of exercise induced collapse caused by insulin shock

Analysis of a parallel multigrid algorithm

The parallel multigrid algorithm of Frederickson and McBryan (1987) is considered. This algorithm uses multiple coarse-grid problems (instead of one problem) in the hope of accelerating convergence and is found to have a close relationship to traditional multigrid methods. Specifically, the parallel coarse-grid correction operator is identical to a traditional multigrid coarse-grid correction operator, except that the mixing of high and low frequencies caused by aliasing error is removed. Appropriate relaxation operators can be chosen to take advantage of this property. Comparisons between the standard multigrid and the new method are made

Entanglement tongue and quantum synchronization of disordered oscillators

We study the synchronization of dissipatively-coupled van der Pol oscillators in the quantum limit, when each oscillator is near its quantum ground state. Two quantum oscillators with different frequencies exhibit an entanglement tongue, which is the quantum analogue of an Arnold tongue. It means that the oscillators are entangled in steady state when the coupling strength is greater than a critical value, and the critical coupling increases with detuning. An ensemble of many oscillators with random frequencies still exhibits a synchronization phase transition in the quantum limit, and we analytically calculate how the critical coupling depends on the frequency disorder. Our results can be experimentally observed with trapped ions or neutral atoms.Comment: 11 pages, 5 figure