Direct multiscale coupling of a transport code to gyrokinetic turbulence codes

Direct coupling between a transport solver and local, nonlinear gyrokinetic calculations using the multiscale gyrokinetic code TRINITY [M. Barnes, Ph.D. thesis, arxiv:0901.2868] is described. The coupling of the microscopic and macroscopic physics is done within the framework of multiscale gyrokinetic theory, of which we present the assumptions and key results. An assumption of scale separation in space and time allows for the simulation of turbulence in small regions of the space-time grid, which are embedded in a coarse grid on which the transport equations are implicitly evolved. This leads to a reduction in computational expense of several orders of magnitude, making first-principles simulations of the full fusion device volume over the confinement time feasible on current computing resources. Numerical results from TRINITY simulations are presented and compared with experimental data from JET and ASDEX Upgrade plasmas.Comment: 12 pages, 13 figures, invited paper for 2009 APS-DPP meeting, submitted to Phys. Plasma

Space Shuttle main engine powerhead structural modeling, stress and fatigue life analysis. Volume 1: Gasdynamic environment of the SSME HPFTP and HPOTP turbines

Gasdynamic analysis for the turbine blades and nozzle vanes, HPFTP turbine analysis, and HPOTP turbine analysis are provided

Linearized model Fokker-Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests

A set of key properties for an ideal dissipation scheme in gyrokinetic simulations is proposed, and implementation of a model collision operator satisfying these properties is described. This operator is based on the exact linearized test-particle collision operator, with approximations to the field-particle terms that preserve conservation laws and an H-Theorem. It includes energy diffusion, pitch-angle scattering, and finite Larmor radius effects corresponding to classical (real-space) diffusion. The numerical implementation in the continuum gyrokinetic code GS2 is fully implicit and guarantees exact satisfaction of conservation properties. Numerical results are presented showing that the correct physics is captured over the entire range of collisionalities, from the collisionless to the strongly collisional regimes, without recourse to artificial dissipation.Comment: 13 pages, 8 figures, submitted to Physics of Plasmas; typos fixe

Quantifying alignment effects in 3D coordinate measurement

The use of fixtureless, non-contact coordinate measurement has become increasingly prevalent in manufacturing problem solving. Manufacturers now routinely use measurement systems such as white light area scanners, photogrammetry, laser trackers, and portable laser scanners to conduct studies that require measuring upstream supplier parts, tooling, or in-process subassemblies. For part measurements in these studies, certified fixtures with alignment features such as tooling balls often are not available. Instead, manufacturers rely on ad hoc part-holding fixtures or measure parts without fixtures and perform alignments mathematically. Here, advancements in software are providing operators with numerous alignment options, and users are actively using this functionality. Naturally, these additional capabilities have led to inconsistencies in the alignment method used across measurement studies, often affecting dimensional results. This paper reviews several common alignment or registration methods and provides a metric to assess systematic alignment error. To demonstrate alignment effects, we present a measurement system study of a moderately complex part requiring an over-constrained datum scheme. We first measure the part using a conventional fixture-based method to establish a baseline for static and dynamic repeatability. We then compare these with results from two mathematically-based iterative alignment methods based on fixtureless measurement. Next, we assess the systematic alignment error between the different fixture/alignment alternatives. We show that for the same basic datum scheme provided on engineering drawings, the systematic alignment error is a far more significant issue for problem solving than the repeatability error or equipment accuracy.http://deepblue.lib.umich.edu/bitstream/2027.42/69245/1/99983.pd

Moment Closure - A Brief Review

Moment closure methods appear in myriad scientific disciplines in the modelling of complex systems. The goal is to achieve a closed form of a large, usually even infinite, set of coupled differential (or difference) equations. Each equation describes the evolution of one "moment", a suitable coarse-grained quantity computable from the full state space. If the system is too large for analytical and/or numerical methods, then one aims to reduce it by finding a moment closure relation expressing "higher-order moments" in terms of "lower-order moments". In this brief review, we focus on highlighting how moment closure methods occur in different contexts. We also conjecture via a geometric explanation why it has been difficult to rigorously justify many moment closure approximations although they work very well in practice.Comment: short survey paper (max 20 pages) for a broad audience in mathematics, physics, chemistry and quantitative biolog

A blind detection of a large, complex, Sunyaev--Zel'dovich structure

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100{\mu}Jy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological {\beta}-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 {\mu}K - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.Comment: accepted MNRAS. 12 pages, 9 figure

Simulation of the Magnetothermal Instability

In many magnetized, dilute astrophysical plasmas, thermal conduction occurs almost exclusively parallel to magnetic field lines. In this case, the usual stability criterion for convective stability, the Schwarzschild criterion, which depends on entropy gradients, is modified. In the magnetized long mean free path regime, instability occurs for small wavenumbers when (dP/dz)(dln T/dz) > 0, which we refer to as the Balbus criterion. We refer to the convective-type instability that results as the magnetothermal instability (MTI). We use the equations of MHD with anisotropic electron heat conduction to numerically simulate the linear growth and nonlinear saturation of the MTI in plane-parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured from the simulations are in excellent agreement with the weak field dispersion relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic fields can damp the growth of the MTI and affect the nonlinear regime. The instability saturates when the atmosphere becomes isothermal as the source of free energy is exhausted. By maintaining a fixed temperature difference between the top and bottom boundaries of the simulation domain, sustained convective turbulence can be driven. MTI-stable layers introduced by isotropic conduction are used to prevent the formation of unresolved, thermal boundary layers. We find that the largest component of the time-averaged heat flux is due to advective motions as opposed to the actual thermal conduction itself. Finally, we explore the implications of this instability for a variety of astrophysical systems, such as neutron stars, the hot intracluster medium of galaxy clusters, and the structure of radiatively inefficient accretion flows.Comment: Accepted for publication in Astrophysics and Space Science as proceedings of the 6th High Energy Density Laboratory Astrophysics (HEDLA) Conferenc

On the Origin of Cosmic Magnetic Fields

We review the literature concerning how the cosmic magnetic fields pervading nearly all galaxies actually got started. some observational evidence involves the chemical abundance of the light elements Be and B, while another one is based on strong magnetic fields seen in high red shift galaxies. Seed fields, whose strength is of order 10^{-20} gauss, easily sprung up in the era preceding galaxy formation. Several mechanisms are proposed to amplify these seed fields to microgauss strengths. The standard mechanism is the Alpha-Omega dynamo theory. It has a major difficulty that makes unlikely to provide the sole origin. The difficulty is rooted in the fact that the total flux is constant. This implies that flux must be removed from the galactic discs. This requires that the field and flux be separated, for otherwise interstellar mass must be removed from the deep galactic gravitational and then their strength increased by the alpha omega theory.Comment: 90 pages and 6 figures; accepted for publication in Reports of Progress in Physics as an invited revie