Overture - Object-oriented tools for overset grid applications

The Overture framework is an object-oriented environment for solving partial differential equations in two and three space dimensions. It is a collection of C++ libraries that enables the use of finite difference and finite volume methods at a level that hides the details of the associated data structures. Overture can be used to solve problems in complicated, moving geometries using the method of overlapping grids. It has support for grid generation, difference operators, boundary conditions, data-base access and graphics. Short sample code segments are presented to show the power of this approach

On the smallest scale for the incompressible Navier-Stokes equations

It is proven that for solutions to the two- and three-dimensional incompressible Navier-Stokes equations the minimum scale is inversely proportional to the square root of the Reynolds number based on the kinematic viscosity and the maximum of the velocity gradients. The bounds on the velocity gradients can be obtained for two-dimensional flows, but have to be assumed to be three-dimensional. Numerical results in two dimensions are given which illustrate and substantiate the features of the proof. Implications of the minimum scale result to the decay rate of the energy spectrum are discussed

Re-introduction of structurally complex wood jams promotes channel and habitat recovery from overwidening: Implications for river conservation

Copyright © 2017 John Wiley & Sons, Ltd. Large wood is a powerful geomorphic agent in rivers, providing important habitat functions for a range of aquatic organisms, but has been subject to a long history of removal. Internationally, approaches to river restoration are increasingly incorporating large wood features, but generally favour simple flow deflectors (e.g. single logs, stripped of branches and anchored in place) over more complex structures that more accurately mimic natural wood jams. This paper explores channel response to wood-based restoration of an overwidened lowland chalk stream that incorporated whole felled trees. Hydraulics, sediment, topography and vegetation data were assessed for a 3year period for two restored reaches: an upstream reach where pre-restoration baseline data were obtained, and a downstream reach restored before data collection. Where pre-restoration data were available, the introduction of wood jams generated sediment deposition within jams leading to the development of vegetated marginal ‘benches’ and bed scour in adjacent areas of flow convergence. Patterns were less clear in the downstream reach, where restoration design was less ambitious and outcomes may have been affected by subsequent restoration work upstream. The results indicate that reintroduction of large wood (whole trees), can promote channel and habitat recovery from overwidening in lowland rivers, creating important ecological benefits through the provision of structurally complex marginal habitat and associated food resources. Longer-term assessments are required to establish whether the trajectories of change are persistent. The work emphasizes the effectiveness of restoration approaches that aim to ‘work with nature’. The ambitious design, incorporating structurally complex wood jams, was also low-cost, using materials available from the river corridor (existing riparian trees). Furthermore, ecosystem engineering effects were amplified by the colonization of wood jams by aquatic vegetation. The approach should, therefore, be transferable to other lowland rivers, subject to wider catchment constraints

Meshfree finite differences for vector Poisson and pressure Poisson equations with electric boundary conditions

We demonstrate how meshfree finite difference methods can be applied to solve vector Poisson problems with electric boundary conditions. In these, the tangential velocity and the incompressibility of the vector field are prescribed at the boundary. Even on irregular domains with only convex corners, canonical nodal-based finite elements may converge to the wrong solution due to a version of the Babuska paradox. In turn, straightforward meshfree finite differences converge to the true solution, and even high-order accuracy can be achieved in a simple fashion. The methodology is then extended to a specific pressure Poisson equation reformulation of the Navier-Stokes equations that possesses the same type of boundary conditions. The resulting numerical approach is second order accurate and allows for a simple switching between an explicit and implicit treatment of the viscosity terms.Comment: 19 pages, 7 figure

Radon and childhood cancer

British Journal of Cancer (2002) 87, 1336–1337. doi:10.1038/sj.bjc.6600671 www.bjcancer.co

New Measurement of Compton Scattering from the Deuteron and an Improved Extraction of the Neutron Electromagnetic Polarizabilities

The electromagnetic polarizabilities of the nucleon are fundamental properties that describe its response to external electric and magnetic fields. They can be extracted from Compton-scattering data --- and have been, with good accuracy, in the case of the proton. In contradistinction, information for the neutron requires the use of Compton scattering from nuclear targets. Here we report a new measurement of elastic photon scattering from deuterium using quasimonoenergetic tagged photons at the MAX IV Laboratory in Lund, Sweden. These first new data in more than a decade effectively double the world dataset. Their energy range overlaps with previous experiments and extends it by 20 MeV to higher energies. An analysis using Chiral Effective Field Theory with dynamical \Delta(1232) degrees of freedom shows the data are consistent with and within the world dataset. After demonstrating that the fit is consistent with the Baldin sum rule, extracting values for the isoscalar nucleon polarizabilities and combining them with a recent result for the proton, we obtain the neutron polarizabilities as \alpha_n = [11.55 +/- 1.25(stat) +/- 0.2(BSR) +/- 0.8(th)] X 10^{-4} fm^3 and \beta_n = [3.65 -/+ 1.25(stat) +/- 0.2(BSR) -/+ 0.8(th)] X 10^{-4} fm3, with \chi^2 = 45.2 for 44 degrees of freedom.Comment: 6 pages, 3 figures, comments from Physical Review Letters Referees addresse

The complex multiscale structure in simulated and observed emission maps of the proto-cluster cloud G0.253+0.016 (\u27the Brick\u27)

The Central Molecular Zone (the central ∼500 pc of the Milky Way) hosts molecular clouds in an extreme environment of strong shear, high gas pressure and density, and complex chemistry. G0.253+0.016, also known as \u27the Brick\u27, is the densest, most compact, and quiescent of these clouds. High-resolution observations with the Atacama Large Millimetre/submillimetre Array (ALMA) have revealed its complex, hierarchical structure. In this paper we compare the properties of recent hydrodynamical simulations of the Brick to those of the ALMA observations. To facilitate the comparison, we post-process the simulations and create synthetic ALMA maps of molecular line emission from eight molecules. We correlate the line emission maps to each other and to the mass column density and find that HNCO is the best mass tracer of the eight emission lines within the simulations. Additionally, we characterize the spatial structure of the observed and simulated cloud using the density probability distribution function (PDF), spatial power spectrum, fractal dimension, and moments of inertia. While we find good agreement between the observed and simulated data in terms of power spectra and fractal dimensions, there are key differences in the density PDFs and moments of inertia, which we attribute to the omission of magnetic fields in the simulations. This demonstrates that the presence of the Galactic potential can reproduce many cloud properties, but additional physical processes are needed to fully explain the gas structure

The three-dimensional structure of Galactic molecular cloud complexes out to 2.5 kpc

Knowledge of the three-dimensional structure of Galactic molecular clouds is important for understanding how clouds are affected by processes such as turbulence and magnetic fields and how this structure effects star formation within them. Great progress has been made in this field with the arrival of the Gaia mission, which provides accurate distances to $\sim10^{9}$ stars. Combining these distances with extinctions inferred from optical-IR, we recover the three-dimensional structure of 16 Galactic molecular cloud complexes at $\sim1$pc resolution using our novel three-dimensional dust mapping algorithm \texttt{Dustribution}. Using \texttt{astrodendro} we derive a catalogue of physical parameters for each complex. We recover structures with aspect ratios between 1 and 11, i.e.\ everything from near-spherical to very elongated shapes. We find a large variation in cloud environments that is not apparent when studying them in two-dimensions. For example, the nearby California and Orion A clouds look similar on-sky, but we find California to be more sheet-like, and massive, which could explain their different star-formation rates. In Carina, our most distant complex, we observe evidence for dust sputtering, which explains its measured low dust mass. By calculating the total mass of these individual clouds, we demonstrate that it is necessary to define cloud boundaries in three-dimensions in order to obtain an accurate mass; simply integrating the extinction overestimates masses. We find that Larson's relationship on mass vs radius holds true whether you assume a spherical shape for the cloud or take their true extents.Comment: accepted for publication by MNRAS, 23 pages, 9 figures, 3 table

Experience with custom processors in space flight applications

The Applied Physics Laboratory (APL) has developed a magnetometer instrument for a swedish satellite named Freja with launch scheduled for August 1992 on a Chinese Long March rocket. The magnetometer controller utilized a custom microprocessor designed at APL with the Genesil silicon compiler. The processor evolved from our experience with an older bit-slice design and two prior single chip efforts. The architecture of our microprocessor greatly lowered software development costs because it was optimized to provide an interactive and extensible programming environment hosted by the target hardware. Radiation tolerance of the microprocessor was also tested and was adequate for Freja's mission -- 20 kRad(Si) total dose and very infrequent latch-up and single event upset events