Populations of models, Experimental Designs and coverage of parameter space by Latin Hypercube and Orthogonal Sampling

In this paper we have used simulations to make a conjecture about the coverage of a $t$ dimensional subspace of a $d$ dimensional parameter space of size $n$ when performing $k$ trials of Latin Hypercube sampling. This takes the form $P(k,n,d,t)=1-e^{-k/n^{t-1}}$. We suggest that this coverage formula is independent of $d$ and this allows us to make connections between building Populations of Models and Experimental Designs. We also show that Orthogonal sampling is superior to Latin Hypercube sampling in terms of allowing a more uniform coverage of the $t$ dimensional subspace at the sub-block size level.Comment: 9 pages, 5 figure

Convergence of eigenfunction expansions corresponding to nonlinear Sturm-Liouville operators

It is well known that the classical linear Sturm-Liouville eigenvalue problem is self-adjoint and possesses a family of eigenfunctions which form an orthonormal basis for the space L. A natural question is to ask if a similar result holds for nonlinear problems. In the present paper, we examine the basis property for eigenfunctions of nonlinear Sturm-Liouville equations subject to general linear, separated boundary conditions

Uncertainty quantification of coal seam gas production prediction using Polynomial Chaos

A surrogate model approximates a computationally expensive solver. Polynomial Chaos is a method to construct surrogate models by summing combinations of carefully chosen polynomials. The polynomials are chosen to respect the probability distributions of the uncertain input variables (parameters); this allows for both uncertainty quantification and global sensitivity analysis. In this paper we apply these techniques to a commercial solver for the estimation of peak gas rate and cumulative gas extraction from a coal seam gas well. The polynomial expansion is shown to honour the underlying geophysics with low error when compared to a much more complex and computationally slower commercial solver. We make use of advanced numerical integration techniques to achieve this accuracy using relatively small amounts of training data

S13RS SGR No. 23 (TSGDRC)

A RESOLUTION To create the Temporary Student Government Document Revision Committe

Study of the acoustic signature of UHE neutrino interactions in water and ice

The production of acoustic signals from the interactions of ultra-high energy (UHE) cosmic ray neutrinos in water and ice has been studied. A new computationally fast and efficient method of deriving the signal is presented. This method allows the implementation of up to date parameterisations of acoustic attenuation in sea water and ice that now includes the effects of complex attenuation, where appropriate. The methods presented here have been used to compute and study the properties of the acoustic signals which would be expected from such interactions. A matrix method of parameterising the signals, which includes the expected fluctuations, is also presented. These methods are used to generate the expected signals that would be detected in acoustic UHE neutrino telescopes.Comment: 21 pages and 13 figure

Positive solutions for nonlinear m-point Eigenvalue problems

In this paper, we are concerned with determining values of lambda, for which there exist positive solutions of the nonlinear eigenvalue problem [GRAPHICS] where a, b, c, d is an element of [0, infinity), xi(i) is an element of (0, 1), alpha(i), beta(i) is an element of [0 infinity) (for i is an element of {1, ..., m - 2}) are given constants, p, q is an element of C ([0, 1], (0, infinity)), h is an element of C ([0, 1], [0, infinity)), and f is an element of C ([0, infinity), [0, infinity)) satisfying some suitable conditions. Our proofs are based on Guo-Krasnoselskii fixed point theorem. (C) 2004 Elsevier Inc. All rights reserved

Discrete first-order three-point boundary value problem

We study difference equations which arise as discrete approximations to three-point boundary value problems for systems of first-order ordinary differential equations. We obtain new results of the existence of solutions to the discrete problem by employing Euler’s method. The existence of solutions are proven by the contraction mapping theorem and the Brouwer fixed point theorem in Euclidean space. We apply our results to show that solutions to the discrete problem converge to solutions of the continuous problem in an aggregate sense. We also give some examples to illustrate the existence of a unique solution of the contraction mapping theorem

A biomechanical model for fibril recruitment: Evaluation in tendons and arteries

Simulations of soft tissue mechanobiological behaviour are increasingly important for clinical prediction of aneurysm, tendinopathy and other disorders. Mechanical behaviour at low stretches is governed by fibril straightening, transitioning into load-bearing at recruitment stretch, resulting in a tissue stiffening effect. Previous investigations have suggested theoretical relationships between stress-stretch measurements and recruitment probability density function (PDF) but not derived these rigorously nor evaluated these experimentally. Other work has proposed image-based methods for measurement of recruitment but made use of arbitrary fibril critical straightness parameters. The aim of this work was to provide a sound theoretical basis for estimating recruitment PDF from stress-stretch measurements and to evaluate this relationship using image-based methods, clearly motivating the choice of fibril critical straightness parameter in rat tail tendon and porcine artery. Rigorous derivation showed that the recruitment PDF may be estimated from the second stretch derivative of the first Piola-Kirchoff tissue stress. Image-based fibril recruitment identified the fibril straightness parameter that maximised Pearson correlation coefficients (PCC) with estimated PDFs. Using these critical straightness parameters the new method for estimating recruitment PDF showed a PCC with image-based measures of 0.915 and 0.933 for tendons and arteries respectively. This method may be used for accurate estimation of fibril recruitment PDF in mechanobiological simulation where fibril-level mechanical parameters are important for predicting cell behaviour