An integral equation method for a boundary value problem arising in unsteady water wave problems

In this paper we consider the 2D Dirichlet boundary value problem for Laplace’s equation in a non-locally perturbed half-plane, with data in the space of bounded and continuous functions. We show uniqueness of solution, using standard Phragmen-Lindelof arguments. The main result is to propose a boundary integral equation formulation, to prove equivalence with the boundary value problem, and to show that the integral equation is well posed by applying a recent partial generalisation of the Fredholm alternative in Arens et al [J. Int. Equ. Appl. 15 (2003) pp. 1-35]. This then leads to an existence proof for the boundary value problem. Keywords. Boundary integral equation method, Water waves, Laplace’s

Hydraulic modeling of runoff over a rough surface under partial inundation

This study presents a numerical method to derive the Darcy- Weisbach friction coefficient for overland flow under partial inundation of surface roughness. To better account for the variable influence of roughness with varying levels of emergence, we model the flow over a network which evolves as the free surface rises. This network is constructed using a height numerical map, based on surface roughness data, and a discrete geometry skeletonization algorithm. By applying a hydraulic model to the flows through this network, local heads, velocities, and Froude and Reynolds numbers over the surface can be estimated. These quantities enable us to analyze the flow and ultimately to derive a bulk friction factor for flow over the entire surface which takes into account local variations in flow quantities. Results demonstrate that although the flow is laminar, head losses are chiefly inertial because of local flow disturbances. The results also emphasize that for conditions of partial inundation, flow resistance varies nonmonotonically but does generally increase with progressive roughness inundation.

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold. (c) 2005 IMACS. Published by Elsevier B.V. All rights reserved.

Inclusion of nonlinear demand-supply relationships within large-scale partial equilibrium linear programming models

This paper presents a new method for the inclusion of nonlinear demand and supply relationships within a linear programming model. An existing method for this purpose is described first and its shortcomings are pointed out before showing how the new approach overcomes those difficulties and how it provides a more accurate and 'smooth' (rather than a kinked) approximation of the nonlinear functions as well as dealing with equilibrium under perfect competition instead of handling just the monopolistic situation. The workings of the proposed method are illustrated by extending a previously available sectoral model for the UK agriculture.

Assays for enhanced activity of low efficacy partial agonists at the D-2 dopamine receptor

Background and purpose: Low efficacy partial agonists at the D-2 dopamine receptor may be useful for treating schizophrenia. In this report we describe a method for assessing the efficacy of these compounds based on stimulation of [S-35]GTP gamma S binding. Experimental approach: Agonist efficacy was assessed from [S-35]GTP gamma S binding to membranes of CHO cells expressing D2 dopamine receptors in buffers with and without Na+. Effects of Na+ on receptor/G protein coupling were assessed using agonist/[H-3] spiperone competition binding assays. Key results: When [S-35]GTP gamma S binding assays were performed in buffers containing Na+, some agonists (aripiprazole, AJ-76, UH-232) exhibited very low efficacy whereas other agonists exhibited measurable efficacy. When Na+ was substituted by N-methyl D-glucamine, the efficacy of all agonists increased (relative to that of dopamine) but particularly for aripiprazole, aplindore, AJ-76, (-)-3-PPP and UH-232. In ligand binding assays, substitution of Na+ by N-methyl D-glucamine increased receptor/G protein coupling for some agonists -. aplindore, dopamine and (-)-3-PPP-but for aripiprazole, AJ-76 and UH-232 there was little effect on receptor/G protein coupling. Conclusions and implications: Substitution of Na+ by NMDG increases sensitivity in [S-35] GTPgS binding assays so that very low efficacy agonists were detected clearly. For some agonists the effect seems to be mediated via enhanced receptor/G protein coupling whereas for others the effect is mediated at another point in the G protein activation cycle. AJ-76, aripiprazole and UH-232 seem particularly sensitive to this change in assay conditions. This work provides a new method to discover these very low efficacy agonists.

Velocity-based moving mesh methods for nonlinear partial differential equations

This article describes a number of velocity-based moving mesh numerical methods formultidimensional nonlinear time-dependent partial differential equations (PDEs). It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation. Finite element algorithms are derived for both mass-conserving and non mass-conserving problems, and results shown for a number of multidimensional nonlinear test problems, including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem. Further applications and extensions are referenced.

GSHMC: an efficient method for molecular simulation

The hybrid Monte Carlo (HMC) method is a popular and rigorous method for sampling from a canonical ensemble. The HMC method is based on classical molecular dynamics simulations combined with a Metropolis acceptance criterion and a momentum resampling step. While the HMC method completely resamples the momentum after each Monte Carlo step, the generalized hybrid Monte Carlo (GHMC) method can be implemented with a partial momentum refreshment step. This property seems desirable for keeping some of the dynamic information throughout the sampling process similar to stochastic Langevin and Brownian dynamics simulations. It is, however, ultimate to the success of the GHMC method that the rejection rate in the molecular dynamics part is kept at a minimum. Otherwise an undesirable Zitterbewegung in the Monte Carlo samples is observed. In this paper, we describe a method to achieve very low rejection rates by using a modified energy, which is preserved to high-order along molecular dynamics trajectories. The modified energy is based on backward error results for symplectic time-stepping methods. The proposed generalized shadow hybrid Monte Carlo (GSHMC) method is applicable to NVT as well as NPT ensemble simulations.

A modified PNN algorithm with optimal PD modeling using the orthogonal least squares method

In this paper a modified algorithm is suggested for developing polynomial neural network (PNN) models. Optimal partial description (PD) modeling is introduced at each layer of the PNN expansion, a task accomplished using the orthogonal least squares (OLS) method. Based on the initial PD models determined by the polynomial order and the number of PD inputs, OLS selects the most significant regressor terms reducing the output error variance. The method produces PNN models exhibiting a high level of accuracy and superior generalization capabilities. Additionally, parsimonious models are obtained comprising a considerably smaller number of parameters compared to the ones generated by means of the conventional PNN algorithm. Three benchmark examples are elaborated, including modeling of the gas furnace process as well as the iris and wine classification problems. Extensive simulation results and comparison with other methods in the literature, demonstrate the effectiveness of the suggested modeling approach.

Identification of partial resetting using De as a function of illumination time

Modern age samples from various depositional environments were examined for signal resetting. For 19 modern aeolian/beach samples all De values obtained were View the MathML source, with ∼70% having View the MathML source. For 21 fluvial/colluvial samples, all De values were View the MathML source with ∼80% being View the MathML source. De as a function of illumination (OSL measurement) time (De(t)) plots were examined for all samples. Based on previous laboratory experiments, increases in De(t) were expected for partially reset samples, and constant De(t) for fully reset samples. All aeolian samples, both modern age and additional ‘young’ samples (<1000 years), showed constant (flat) De(t) while all modern, non-zero De, fluvial/colluvial samples showed increasing De(t). ‘Replacement plots’, where a regenerated signal is substituted for the natural, yielded constant (flat) De(t). These findings support strongly the use of De(t) as a method of identifying incomplete resetting in fluvial samples. Potential complicating factors, such as illumination (bleaching) spectrum, thermal instability and component composition are discussed and a series of internal checks on the applicability of the De(t) for each individual aliquot/grain level are outlined.