Split operator finite element method for modelling pulmonary gas exchange

A pulmonary gas exchange model which contains six species of interest in three different regions is considered. The model leads to a system of coupled non-linear parabolic partial differential equations and is solved using a transport-reaction split operator scheme. A Galerkin weighted residual finite element method solves the transport algorithm while a simple forward time step is used for the reaction algorithm. Using different time scales for the transport and reaction algorithms, we obtain a reasonable approximation of gas exchange

Using the stochastic Galerkin method as a predictive tool during an epidemic

The ability to accurately predict the course of an epidemic is extremely important. This article looks at an influenza outbreak that spread through a small boarding school. Predictions are made on multiple days throughout the epidemic using the stochastic Galerkin method to consider a range of plausible values for the parameters. These predictions are then compared to known data points. Predictions made before the peak of the epidemic had much larger variances compared to predictions made after the peak of the epidemic. References B. M. Chen-Charpentier, J. C. Cortes, J. V. Romero, and M. D. Rosello. Some recommendations for applying gPC (generalized polynomial chaos) to modeling: An analysis through the Airy random differential equation. Applied Mathematics and Computation, 219(9):4208 – 4218, 2013. doi:10.1016/j.amc.2012.11.007 B. M. Chen-Charpentier and D. Stanescu. Epidemic models with random coefficients. Mathematical and Computer Modelling, 52:1004 – 1010, 2010. doi:10.1016/j.mcm.2010.01.014 D. B. Harman and P. R. Johnston. Applying the stochastic galerkin method to epidemic models with individualised parameter distributions. In Proceedings of the 12th Biennial Engineering Mathematics and Applications Conference, EMAC-2015, volume 57 of ANZIAM J., pages C160–C176, August 2016. doi:10.21914/anziamj.v57i0.10394 D. B. Harman and P. R. Johnston. Applying the stochastic galerkin method to epidemic models with uncertainty in the parameters. Mathematical Biosciences, 277:25 – 37, 2016. doi:10.1016/j.mbs.2016.03.012 D. B. Harman and P. R. Johnston. Boarding house: find border. 2019. doi:10.6084/m9.figshare.7699844.v1 D. B. Harman and P. R. Johnston. SIR uniform equations. 2 2019. doi:10.6084/m9.figshare.7692392.v1 H. W. Hethcote. The mathematics of infectious diseases. SIAM Review, 42(4):599–653, 2000. doi:10.1137/S0036144500371907 R.I. Hickson and M.G. Roberts. How population heterogeneity in susceptibility and infectivity influences epidemic dynamics. Journal of Theoretical Biology, 350(0):70 – 80, 2014. doi:10.1016/j.jtbi.2014.01.014 W. O. Kermack and A. G. McKendrick. A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society of London. Series A, 115(772):700–721, August 1927. doi:10.1098/rspa.1927.0118 M. G. Roberts. A two-strain epidemic model with uncertainty in the interaction. The ANZIAM Journal, 54:108–115, 10 2012. doi:10.1017/S1446181112000326 M. G. Roberts. Epidemic models with uncertainty in the reproduction number. Journal of Mathematical Biology, 66(7):1463–1474, 2013. doi:10.1007/s00285-012-0540-y F. Santonja and B. Chen-Charpentier. Uncertainty quantification in simulations of epidemics using polynomial chaos. Computational and Mathematical Methods in Medicine, 2012:742086, 2012. doi:10.1155/2012/742086 Communicable Disease Surveillance Centre (Public Health Laboratory Service) and Communicable Diseases (Scotland) Unit. Influenza in a boarding school. BMJ, 1(6112):587, 1978. doi:10.1136/bmj.1.6112.586 G. Strang. Linear Algebra and Its Applications. Thomson, Brooks/Cole, 2006. D. Xiu. Numerical Methods for Stochastic Computations: A Spectral Method Approach. Princeton University Press, 2010

Broadband Records of Earthquakes in Deep Gold Mines and a Comparison with Results from SAFOD, California

For one week during September 2007, we deployed a temporary network of field recorders and accelerometers at four sites within two deep, seismically active mines. The ground-motion data, recorded at 200 samples/sec, are well suited to determining source and ground-motion parameters for the mining-induced earthquakes within and adjacent to our network. Four earthquakes with magnitudes close to 2 were recorded with high signal/noise at all four sites. Analysis of seismic moments and peak velocities, in conjunction with the results of laboratory stick-slip friction experiments, were used to estimate source processes that are key to understanding source physics and to assessing underground seismic hazard. The maximum displacements on the rupture surfaces can be estimated from the parameter Rv, where v is the peak ground velocity at a given recording site, and R is the hypocentral distance. For each earthquake, the maximum slip and seismic moment can be combined with results from laboratory friction experiments to estimate the maximum slip rate within the rupture zone. Analysis of the four M 2 earthquakes recorded during our deployment and one of special interest recorded by the in-mine seismic network in 2004 revealed maximum slips ranging from 4 to 27 mm and maximum slip rates from 1.1 to 6:3 m=sec. Applying the same analyses to an M 2.1 earthquake within a cluster of repeating earthquakes near the San Andreas Fault Observatory at Depth site, California, yielded similar results for maximum slip and slip rate, 14 mm and 4:0 m=sec

Molt Migration of Scoters at Cape Peirce, Alaska

There is presently little specific information on the molt migrations of scoters in the Nearctic. We conducted migration watches from 21 June to 31 July 1984 (total of 96 h) and from 5 to 15 July 1985 (total of 36 h) during daylight hours to estimate abundance and species composition of scoters engaged in a molt migration at Cape Peirce, southeast Bering Sea, Alaska. We counted 22,897 scoters moving west past the observation site in 1984, the year in which we probably made observations over most of the migration period. Allowing for daylight hours without watches, an estimated 66,500 scoters passed in 1984. Peak passage occurred on 11 July in both years. Species composition of migrants in 1984 was approximately 77% white-winged scoters (Melanitta fusca), 12% surf scoters (M. perspicillata), and 11% black scoters (M. nigra). Most of the migrants were adult males, probably migrating from breeding grounds in Interior Alaska. Our findings suggest that a large and presently undescribed molting area of white-winged scoters exists somewhere in the waters of western Alaska or eastern Siberia.Key words: scoters, seaducks, molt migration, Cape Peirce, Togiak National Wildlife Refuge, Bering Sea, AlaskaMots clés: macreuses, canards de mer, migration d’oiseaux ayant effectué la mue, Cape Peirce, Togiak National Wildlife Refuge, mer de Béring, Alask

Randomised controlled feasibility trial of an evidence-informed behavioural intervention for obese adults with additional risk factors

Peer reviewedPublisher PD

Phase equilibria constraints on Archean crustal genesis from crystallization experiments on trondhjemite with water at 10-17 kbar

The formation of continental crust during the Archean and early Proterozoic occurred through a different mechanisms than the currently active processes of calc-alkaline volcanism in orogenic regions. In view that most crustal growth models imply that by the end of the Archean a continental mass equivalent to 75% or more of the current crust had evolved, it seems highly relevant to study early crustal genesis

HST-NICMOS Observations of Terzan 5: Stellar Content and Structure of the Core

We report results from HST-NICMOS imaging of the extremely dense core of the globular cluster Terzan 5. This highly obscured bulge cluster has been estimated to have one of the highest collision rates of any galactic globular cluster, making its core a particularly conducive environment for the production of interacting binary systems. We have reconstructed high-resolution images of the central 19"x19" region of Terzan 5 by application of the drizzle algorithm to dithered NIC2 images in the F110W, F187W, and F187N near-infrared filters. We have used a DAOPHOT/ALLSTAR analysis of these images to produce the deepest color-magnitude diagram (CMD) yet obtained for the core of Terzan 5. We have also analyzed the parallel 11"X11" NIC1 field, centered 30" from the cluster center and imaged in F110W and F160W, and an additional NIC2 field that is immediately adjacent to the central field. This imaging results in a clean detection of the red-giant branch and horizontal branch in the central NIC2 field, and the detection of these plus the main-sequence turnoff and the upper main sequence in the NIC1 field. We have constructed an H versus J-H CMD for the NIC1 field. We obtain a new distance estimate of 8.7 kpc, which places Terzan 5 within less than 1 kpc of the galactic center. We have also determined a central surface-density profile which results in a maximum likelihood estimate of 7.9" +/- 0.6" for the cluster core radius. We discuss the implications of these results for the dynamical state of Terzan 5.Comment: 17 pages, 9 figures, accepted for publication in ApJ, for May 20, 200

A Two Micron All-Sky Survey View of the Sagittarius Dwarf Galaxy: II. Swope Telescope Spectroscopy of M Giant Stars in the Dynamically Cold Sagittarius Tidal Stream

We present moderate resolution (~6 km/s) spectroscopy of 284 M giant candidates selected from the Two Micron All Sky Survey photometry. Radial velocities (RVs) are presented for stars mainly in the south, with a number having positions consistent with association to the trailing tidal tail of the Sagittarius (Sgr) dwarf galaxy. The latter show a clear RV trend with orbital longitude, as expected from models of the orbit and destruction of Sgr. A minimum 8 kpc width of the trailing stream about the Sgr orbital midplane is implied by verified RV members. The coldness of this stream (dispersion ~10 km/s) provides upper limits on the combined contributions of stream heating by a lumpy Galactic halo and the intrinsic dispersion of released stars, which is a function of the Sgr core mass. The Sgr trailing arm is consistent with a Galactic halo containing one dominant, LMC-like lump, however some lumpier halos are not ruled out. An upper limit to the total M/L of the Sgr core is 21 in solar units. A second structure that roughly mimics expectations for wrapped, leading Sgr arm debris crosses the trailing arm in the Southern Hemisphere; however, this may also be an unrelated tidal feature. Among the <13 kpc M giants toward the South Galactic Pole are some with large RVs that identify them as halo stars, perhaps part of the Sgr leading arm near the Sun. The positions and RVs of Southern Hemisphere M giants are compared with those of southern globular clusters potentially stripped from the Sgr system and support for association of Pal 2 and Pal 12 with Sgr debris is found. Our discussion includes description of a masked-filtered cross-correlation methodology that achieves better than 1/20 of a resolution element RVs in moderate resolution spectra.Comment: 41 pages, 6 figures, Astronomical Journal, in press (submitted Nov. 24, 2003; tentatively scheduled for July 2004 issue