Consistent thermodynamic derivative estimates for tabular equations of state

Numerical simulations of compressible fluid flows require an equation of state (EOS) to relate the thermodynamic variables of density, internal energy, temperature, and pressure. A valid EOS must satisfy the thermodynamic conditions of consistency (derivation from a free energy) and stability (positive sound speed squared). When phase transitions are significant, the EOS is complicated and can only be specified in a table. For tabular EOS's such as SESAME from Los Alamos National Laboratory, the consistency and stability conditions take the form of a differential equation relating the derivatives of pressure and energy as functions of temperature and density, along with positivity constraints. Typical software interfaces to such tables based on polynomial or rational interpolants compute derivatives of pressure and energy and may enforce the stability conditions, but do not enforce the consistency condition and its derivatives. We describe a new type of table interface based on a constrained local least squares regression technique. It is applied to several SESAME EOS's showing how the consistency condition can be satisfied to round-off while computing first and second derivatives with demonstrated second-order convergence. An improvement of 14 orders of magnitude over conventional derivatives is demonstrated, although the new method is apparently two orders of magnitude slower, due to the fact that every evaluation requires solving an 11-dimensional nonlinear system.Comment: 29 pages, 9 figures, 16 references, submitted to Phys Rev

Solving One Dimensional Scalar Conservation Laws by Particle Management

We present a meshfree numerical solver for scalar conservation laws in one space dimension. Points representing the solution are moved according to their characteristic velocities. Particle interaction is resolved by purely local particle management. Since no global remeshing is required, shocks stay sharp and propagate at the correct speed, while rarefaction waves are created where appropriate. The method is TVD, entropy decreasing, exactly conservative, and has no numerical dissipation. Difficulties involving transonic points do not occur, however inflection points of the flux function pose a slight challenge, which can be overcome by a special treatment. Away from shocks the method is second order accurate, while shocks are resolved with first order accuracy. A postprocessing step can recover the second order accuracy. The method is compared to CLAWPACK in test cases and is found to yield an increase in accuracy for comparable resolutions.Comment: 15 pages, 6 figures. Submitted to proceedings of the Fourth International Workshop Meshfree Methods for Partial Differential Equation

Visuo-spatial cognition in Williams syndrome: Reviewing and accounting for the strengths and weaknesses in performance

Individuals with Williams syndrome typically show relatively poor visuo-spatial abilities in comparison to stronger verbal skills. However, individuals' level of performance is not consistent across all visuo-spatial tasks. The studies assessing visuo-spatial functioning in Williams syndrome are critically reviewed, in order to provide a clear pattern of the relative difficulty of these tasks. This prompts a possible explanation of the variability in performance seen which focuses on the processing demands of some of these tasks. Individuals with Williams syndrome show an atypical processing style on tests of construction, which does not affect tests of perception

Thermochemistry of Alane Complexes for Hydrogen Storage: A Theoretical and Experimental Comparison

Knowledge of the relative stabilities of alane (AlH3) complexes with electron donors is essential for identifying hydrogen storage materials for vehicular applications that can be regenerated by off-board methods; however, almost no thermodynamic data are available to make this assessment. To fill this gap, we employed the G4(MP2) method to determine heats of formation, entropies, and Gibbs free energies of formation for thirty-eight alane complexes with NH3-nRn (R = Me, Et; n = 0-3), pyridine, pyrazine, triethylenediamine (TEDA), quinuclidine, OH2-nRn (R = Me, Et; n = 0-2), dioxane, and tetrahydrofuran (THF). Monomer, bis, and selected dimer complex geometries were considered. Using these data, we computed the thermodynamics of the key formation and dehydrogenation reactions that would occur during hydrogen delivery and alane regeneration, from which trends in complex stability were identified. These predictions were tested by synthesizing six amine-alane complexes involving trimethylamine, triethylamine, dimethylethylamine, TEDA, quinuclidine, and hexamine, and obtaining upper limits of delta G for their formation from metallic aluminum. Combining these computational and experimental results, we establish a criterion for complex stability relevant to hydrogen storage that can be used to assess potential ligands prior to attempting synthesis of the alane complex. Based on this, we conclude that only a subset of the tertiary amine complexes considered and none of the ether complexes can be successfully formed by direct reaction with aluminum and regenerated in an alane-based hydrogen storage system.Comment: Accepted by the Journal of Physical Chemistry

Improving smoothed particle hydrodynamics with an integral approach to calculating gradients

In this paper we develop and test a fully conservative SPH scheme based on a tensor formulation that can be applied to simulate astrophysical systems. In the proposed scheme, derivatives are calculated from an integral expression that leads to a tensor (instead of a vectorial) estimation of gradients and reduces to the standard formulation in the continuum limit. The new formulation improves the interpolation of physical magnitudes, leading to a set of conservative equations that resembles those of standard SPH. The resulting scheme is verified using a variety of well-known tests, all of them simulated in two dimensions. We also discuss an application of the proposed tensor method to astrophysics by simulating the stability of a Sun-like polytrope calculated in three dimensions.Comment: 15 pages, 12 figures, accepted for publication in Astronomy and Astrophysics; minor changes due to language editin

A trial on unruptured intracranial aneurysms (the TEAM trial): results, lessons from a failure and the necessity for clinical care trials

The trial on endovascular management of unruptured intracranial aneurysms (TEAM), a prospective randomized trial comparing coiling and conservative management, initiated in September 2006, was stopped in June 2009 because of poor recruitment (80 patients). Aspects of the trial design that may have contributed to this failure are reviewed in the hope of identifying better ways to successfully complete this special type of pragmatic trial which seeks to test two strategies that are in routine clinical use. Cultural, conceptual and bureaucratic hurdles and difficulties obstruct all trials. These obstacles are however particularly misplaced when the trial aims to identify what a good medical practice should be. A clean separation between research and practice, with diverging ethical and scientific requirements, has been enforced for decades, but it cannot work when care needs to be provided in the presence of pervasive uncertainty. Hence valid and robust scientific methods need to be legitimately re-integrated into clinical practice when reliable knowledge is in want

Integrating cancer survivors' experiences into UK cancer registries: design and development of the ePOCS system (electronic Patient-reported Outcomes from Cancer Survivors)

BACKGROUND: Understanding the psychosocial challenges of cancer survivorship, and identifying which patients experience ongoing difficulties, is a key priority. The ePOCS (electronic patient-reported outcomes from cancer survivors) project aims to develop and evaluate a cost-efficient, UK-scalable electronic system for collecting patient-reported outcome measures (PROMs), at regular post-diagnostic timepoints, and linking these with clinical data in cancer registries. METHODS: A multidisciplinary team developed the system using agile methods. Design entailed process mapping the system's constituent parts, data flows and involved human activities, and undertaking usability testing. Informatics specialists built new technical components, including a web-based questionnaire tool and tracking database, and established component-connecting data flows. Development challenges were overcome, including patient usability and data linkage and security. RESULTS: We have developed a system in which PROMs are completed online, using a secure questionnaire administration tool, accessed via a public-facing website, and the responses are linked and stored with clinical registry data. Patient monitoring and communications are semiautomated via a tracker database, and patient correspondence is primarily Email-based. The system is currently honed for clinician-led hospital-based patient recruitment. CONCLUSIONS: A feasibility test study is underway. Although there are possible challenges to sustaining and scaling up ePOCS, the system has potential to support UK epidemiological PROMs collection and clinical data linkage

Episodic Occurrence of Favourable Weather Constrains Recovery of a Cold Desert Shrubland After Fire

Key to the long-term resilience of dryland ecosystems is the recovery of foundation plant species following disturbance. In ecosystems with high interannual weather variability, understanding the influence of short-term environmental conditions on establishment of foundation species is essential for identifying vulnerable landscapes and developing restoration strategies. We asked how annual environmental conditions affect post-fire establishment of Artemisia tridentata, a shrub species that dominates landscapes across much of the western United States, and evaluated the influence of episodic establishment on population recovery. We collected A. tridentata stem samples from 33 plots in 12 prescribed fire sites that burned 8–11 years before sampling. We determined individual establishment years using annual growth rings. We measured seasonal soil environmental conditions at the study sites and asked if these conditions predicted annual establishment density. We then evaluated whether establishment patterns could be predicted by site-level climate or dominant subspecies. Finally, we tested the effect of the magnitude and frequency of post-fire establishment episodes on long-term population recovery. Annual post-fire recruitment of A. tridentata was driven by the episodic availability of spring soil moisture. Annual establishment was highest with wetter spring soils (relative influence [RI] = 19.4%) and later seasonal dry-down (RI = 11.8%) in the year of establishment. Establishment density declined greatly 4 to 5 years after fire (RI = 17.1%). Post-fire establishment patterns were poorly predicted by site-level mean climate (marginal R2 ≤ 0.18) and dominant subspecies (marginal R2 ≤ 0.43). Population recovery reflected the magnitude, but not the frequency, of early post-fire establishment pulses. Post-fire A. tridentata density and cover (measured 8–11 years after fire) were more strongly related to the magnitude of the largest establishment pulse than to establishment frequency, suggesting that population recovery may occur with a single favourable establishment year. Synthesis and applications. This study demonstrates the importance of episodic periods of favourable weather for long-term plant population recovery following disturbance. Management strategies that increase opportunities for seed availability to coincide with favourable weather conditions, such as retaining unburned patches or repeated seeding treatments, can improve restoration outcomes in high-priority areas