A Parallel SSOR Preconditioner for Lattice QCD

A parallelizable SSOR preconditioning scheme for Krylov subspace iterative solvers in lattice QCD applications involving Wilson fermions is presented. In actual Hybrid Monte Carlo and quark propagator calculations it helps to reduce the number of iterations by a factor of 2 compared to conventional odd-even preconditioning. This corresponds to a gain in cpu-time of 30\% - 70\% over odd-even preconditioning.Comment: Talk presented at LATTICE96(algorithms), 3 pages, LaTeX file, 3 epsf-files include

Hitting Time of Quantum Walks with Perturbation

The hitting time is the required minimum time for a Markov chain-based walk (classical or quantum) to reach a target state in the state space. We investigate the effect of the perturbation on the hitting time of a quantum walk. We obtain an upper bound for the perturbed quantum walk hitting time by applying Szegedy's work and the perturbation bounds with Weyl's perturbation theorem on classical matrix. Based on the definition of quantum hitting time given in MNRS algorithm, we further compute the delayed perturbed hitting time (DPHT) and delayed perturbed quantum hitting time (DPQHT). We show that the upper bound for DPQHT is actually greater than the difference between the square root of the upper bound for a perturbed random walk and the square root of the lower bound for a random walk.Comment: 9 page

CFD investigation of a complete floating offshore wind turbine

This chapter presents numerical computations for floating offshore wind turbines for a machine of 10-MW rated power. The rotors were computed using the Helicopter Multi-Block flow solver of the University of Glasgow that solves the Navier-Stokes equations in integral form using the arbitrary Lagrangian-Eulerian formulation for time-dependent domains with moving boundaries. Hydrodynamic loads on the support platform were computed using the Smoothed Particle Hydrodynamics method. This method is mesh-free, and represents the fluid by a set of discrete particles. The motion of the floating offshore wind turbine is computed using a Multi-Body Dynamic Model of rigid bodies and frictionless joints. Mooring cables are modelled as a set of springs and dampers. All solvers were validated separately before coupling, and the loosely coupled algorithm used is described in detail alongside the obtained results

A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection-diffusion PDE's coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton-Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that on be able to solve chemical equilibrium problems (and compute derivatives), without having to know the solution method. An additional advantage of the Newton-Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.Comment: Computational Geosciences (2009) http://www.springerlink.com/content/933p55085742m203/?p=db14bb8c399b49979ba8389a3cae1b0f&pi=1

SESAM and TXL Results for Wilson Action--A Status Report

Results from two studies of full QCD with two flavours of dynamical Wilson fermions are presented. At beta=5.6, the region 0.83 > m_pi/m_rho > 0.56 at m_pia > 0.23 L^{-1} is explored. The SESAM collaboration has generated ensembles of about 200 statistically independent configurations on a 16^3 x 32-lattice at three different kappa-values and is entering the final phase of data analysis. The TXL simulation on a 24^3 x 40-lattice at two kappa-values has reached half statistics and data analysis has started recently, hence most results presented here are preliminary. The focus of this report is fourfold: we demonstrate that algorithmic improvements like fast Krylov solvers and parallel preconditioning recently introduced can be put into practise in full QCD simulations, we present encouraging observations as to the critical dynamics of the Hybrid Monte Carlo algorithm in the approach to the chiral limit, we mention signal improvements of noisy estimator techniques for disconnected diagrams to the pi-N sigma term, and we report on SESAM's results for light hadron spectrum, light quark masses, and heavy quarkonia.Comment: 24 pages, tex + postscript figures, to appear in Proceedings of Int. Workshop "Lattice QCD on Parallel Computers", University of Tsukuba, Japa

Grid-Adapted FUN3D Computations for the Second High Lift Prediction Workshop

Contributions of the unstructured Reynolds-averaged Navier-Stokes code FUN3D to the 2nd AIAA CFD High Lift Prediction Workshop are described, and detailed comparisons are made with experimental data. Using workshop-supplied grids, results for the clean wing configuration are compared with results from the structured code CFL3D Using the same turbulence model, both codes compare reasonably well in terms of total forces and moments, and the maximum lift is similarly over-predicted for both codes compared to experiment. By including more representative geometry features such as slat and flap brackets and slat pressure tube bundles, FUN3D captures the general effects of the Reynolds number variation, but under-predicts maximum lift on workshop-supplied grids in comparison with the experimental data, due to excessive separation. However, when output-based, off-body grid adaptation in FUN3D is employed, results improve considerably. In particular, when the geometry includes both brackets and the pressure tube bundles, grid adaptation results in a more accurate prediction of lift near stall in comparison with the wind-tunnel data. Furthermore, a rotation-corrected turbulence model shows improved pressure predictions on the outboard span when using adapted grids

Therapeutic impact of cytoreductive surgery and irradiation of posterior fossa ependymoma in the molecular era: a retrospective multicohort analysis

PURPOSE: Posterior fossa ependymoma comprises two distinct molecular variants termed EPN_PFA and EPN_PFB that have a distinct biology and natural history. The therapeutic value of cytoreductive surgery and radiation therapy for posterior fossa ependymoma after accounting for molecular subgroup is not known. METHODS: Four independent nonoverlapping retrospective cohorts of posterior fossa ependymomas (n = 820) were profiled using genome-wide methylation arrays. Risk stratification models were designed based on known clinical and newly described molecular biomarkers identified by multivariable Cox proportional hazards analyses. RESULTS: Molecular subgroup is a powerful independent predictor of outcome even when accounting for age or treatment regimen. Incompletely resected EPN_PFA ependymomas have a dismal prognosis, with a 5-year progression-free survival ranging from 26.1% to 56.8% across all four cohorts. Although first-line (adjuvant) radiation is clearly beneficial for completely resected EPN_PFA, a substantial proportion of patients with EPN_PFB can be cured with surgery alone, and patients with relapsed EPN_PFB can often be treated successfully with delayed external-beam irradiation. CONCLUSION: The most impactful biomarker for posterior fossa ependymoma is molecular subgroup affiliation, independent of other demographic or treatment variables. However, both EPN_PFA and EPN_PFB still benefit from increased extent of resection, with the survival rates being particularly poor for subtotally resected EPN_PFA, even with adjuvant radiation therapy. Patients with EPN_PFB who undergo gross total resection are at lower risk for relapse and should be considered for inclusion in a randomized clinical trial of observation alone with radiation reserved for those who experience recurrence