Numerical solution of the Dirichlet initial boundary value problem for the heat equation in exterior 3-dimensional domains using integral equations

A numerical method for the Dirichlet initial boundary value problem for the heat equation in the exterior and unbounded region of a smooth closed simply connected 3-dimensional domain is proposed and investigated. This method is based on a combination of a Laguerre transformation with respect to the time variable and an integral equation approach in the spatial variables. Using the Laguerre transformation in time reduces the parabolic problem to a sequence of stationary elliptic problems which are solved by a boundary layer approach giving a sequence of boundary integral equations of the first kind to solve. Under the assumption that the boundary surface of the solution domain has a one-to-one mapping onto the unit sphere, these integral equations are transformed and rewritten over this sphere. The numerical discretisation and solution are obtained by a discrete projection method involving spherical harmonic functions. Numerical results are included.

Numerical solution of a Cauchy problem for Laplace equation in 3-dimensional domains by integral equations

A numerical method based on integral equations is proposed and investigated for the Cauchy problem for the Laplace equation in 3-dimensional smooth bounded doubly connected domains. To numerically reconstruct a harmonic function from knowledge of the function and its normal derivative on the outer of two closed boundary surfaces, the harmonic function is represented as a single-layer potential. Matching this representation against the given data, a system of boundary integral equations is obtained to be solved for two unknown densities. This system is rewritten over the unit sphere under the assumption that each of the two boundary surfaces can be mapped smoothly and one-to-one to the unit sphere. For the discretization of this system, Weinert’s method (PhD, Göttingen, 1990) is employed, which generates a Galerkin type procedure for the numerical solution, and the densities in the system of integral equations are expressed in terms of spherical harmonics. Tikhonov regularization is incorporated, and numerical results are included showing the efficiency of the proposed procedure.

Computation of rare transitions in the barotropic quasi-geostrophic equations

We investigate the theoretical and numerical computation of rare transitions in simple geophysical turbulent models. We consider the barotropic quasi-geostrophic and two-dimensional Navier–Stokes equations in regimes where bistability between two coexisting large-scale attractors exist. By means of large deviations and instanton theory with the use of an Onsager–Machlup path integral formalism for the transition probability, we show how one can directly compute the most probable transition path between two coexisting attractors analytically in an equilibrium (Langevin) framework and numerically otherWe adapt a class of numerical optimization algorithms known as minimum action methods to simple geophysical turbulent models. We show that by numerically minimizing an appropriate action functional in a large deviation limit, one can predict the most likely transition path for a rare transition between two states. By considering examples where theoretical predictions can be made, we show that the minimum action method successfully predicts the most likely transition path. Finally, we discuss the application and extension of such numerical optimization schemes to the computation of rare transitions observed in direct numerical simulations and experiments and to other, more complex, turbulent systems.

Langevin dynamics, large deviations and instantons for the quasi-geostrophic model and two-dimensional Euler equations

We investigate a class of simple models for Langevin dynamics of turbulent flows, including the one-layer quasi-geostrophic equation and the two-dimensional Euler equations. Starting from a path integral representation of the transition probability, we compute the most probable fluctuation paths from one attractor to any state within its basin of attraction. We prove that such fluctuation paths are the time reversed trajectories of the relaxation paths for a corresponding dual dynamics, which are also within the framework of quasi-geostrophic Langevin dynamics. Cases with or without detailed balance are studied. We discuss a specific example for which the stationary measure displays either a second order (continuous) or a first order (discontinuous) phase transition and a tricritical point. In situations where a first order phase transition is observed, the dynamics are bistable. Then, the transition paths between two coexisting attractors are instantons (fluctuation paths from an attractor to a saddle), which are related to the relaxation paths of the corresponding dual dynamics. For this example, we show how one can analytically determine the instantons and compute the transition probabilities for rare transitions between two attractors.

Lokális energiamódszer kicsi rendben gerjesztett Liénard-egyenletekre (Local energy method for little order forced equations)

Az x''+f(x) x'+g(x) = 0 alakú Liénard-típusú differenciálegyenlet központi szerepet játszik az üzleti ciklusok Káldor-Kalecki-féle [3,4] és Goodwin-féle [2] modelljeiben, sőt egy a munkanélküliség és vállalkozás-ösztönzések ciklikus változásait leíró újabb modellben [1] is. De ugyanez a nemlineáris egyenlettípus a gerjesztett ingák és elektromos rezgőkörök elméletét is felöleli [5]. Az ezzel kapcsolatos irodalom nagyrészt a határciklusok létezését vizsgálja (pl. [5]), pedig az alapvető stabilitási kérdések jóval áttekinthetőbb módon kezelhetők, s a kapott eredmények közvetve a határciklusok létezésének feltételeit is sokkal jobban be tudják határolni. Jelen dolgozatban az egyváltozós analízis hatékony nyelvezetével olyan egyszerűen megfogalmazható eredményekhez jutunk, amelyek képesek kitágítani az üzleti és más közgazdasági ciklusok modelljeinek kereteit, illetve pl. az [1]-beli modellhez újabb szemléltető speciális eseteket is nyerünk. ____ The Liénard type differential equation of the form x00 + f(x) ¢ x0 + g(x) = 0 has a central role in business cycle models by Káldor [3], Kalecki [4] and Goodwin [2], moreover in a new model describing the cyclical behavior of unemployment and entrepreneurship [1]. The same type of nonlinear equation explains the features of forced pendulums and electric circuits [5]. The related literature discusses mainly the existence of limit cycles, although the fundamental stability questions of this topic can be managed much more easily. The achieved results also outline the conditions for the existence of limit cycles. In this work, by the effective language of real valued analysis, we obtain easy-formulated results which may broaden the frames of economic and business cycle models, moreover we may gain new illustrative particular cases for e.g., [1].

An investigation into equations for estimating water requirements and the development of new equations for predicting total water intake

The primary purpose of this study was to investigate agreement among five equations by which clinicians estimate water requirements (EWR) and to determine how well these equations predict total water intake (TWI). The Institute of Medicine has used TWI as a measure of water requirements. A secondary goal of this study was to develop practical equations to predict TWI. These equations could then be considered accurate predictors of an individual’s water requirement. ^ Regressions were performed to determine agreement between the five equations and between the five equations and TWI using NHANES 1999–2004. The criteria for agreement was (1) strong correlation coefficients between all comparisons and (2) regression line that was not significantly different when compared to the line of equality (x=y) i.e., the 95% CI of the slope and intercept must include one and zero, respectively. Correlations were performed to determine association between fat-free mass (FFM) and TWI. Clinically significant variables were selected to build equations for predicting TWI. All analyses were performed with SAS software and were weighted to account for the complex survey design and for oversampling. ^ Results showed that the five EWR equations were strongly correlated but did not agree with each other. Further, the EWR equations were all weakly associated to TWI and lacked agreement with TWI. The strongest agreement between the NRC equation and TWI explained only 8.1% of the variability of TWI. Fat-free mass was positively correlated to TWI. Two models were created to predict TWI. Both models included the variables, race/ethnicity, kcals, age, and height, but one model also included FFM and gender. The other model included BMI and osmolality. Neither model accounted for more than 28% of the variability of TWI. These results provide evidence that estimates of water requirements would vary depending upon which EWR equation was selected by the clinician. None of the existing EWR equations predicted TWI, nor could a prediction equation be created which explained a satisfactory amount of variance in TWI. A good estimate of water requirements may not be predicted by TWI. Future research should focus on using more valid measures to predict water requirements.^

Simultaneous high-resolution pH and spectrophotometric recordings of oxygen binding in native blood microvolumes, link to supplementary material

Oxygen equilibrium curves have been widely used to understand oxygen transport in numerous organisms. A major challenge has been to monitor oxygen binding characteristics and concomitant pH changes as they occur in vivo, in limited sample volumes. Here we report a technique allowing highly resolved and simultaneous monitoring of pH and blood pigment saturation in minute blood volumes. We equipped a gas diffusion chamber with a broad range fibre optic spectrophotometer and a micro-pH optode and recorded changes of pigment oxygenation along PO2 and pH gradients to test the setup. Oxygen binding parameters derived from measurements in only 15 µl of haemolymph from the cephalopod Octopus vulgaris showed low instrumental error (0.93%) and good agreement with published data. Broad range spectra, each resolving 2048 data points, provided detailed insight into the complex absorbance characteristics of diverse blood types. After consideration of photobleaching and intrinsic fluorescence, pH optodes yielded accurate recordings and resolved a sigmoidal shift of 0.03 pH units in response to changing PO2 from 0-21 kPa. Highly resolved continuous recordings along pH gradients conformed to stepwise measurements at low rates of pH changes. In this study we showed that a diffusion chamber upgraded with a broad range spectrophotometer and an optical pH sensor accurately characterizes oxygen binding with minimal sample consumption and manipulation. We conclude that the modified diffusion chamber is highly suitable for experimental biologists who demand high flexibility, detailed insight into oxygen binding as well as experimental and biological accuracy combined in a single set up.

Spectrophotometric study in the near-ir of a sample of MIPS selected galaxies at Z~2

Our main objective is to determine what kind of galaxies dominate the cosmic SFR density at z~2. Our sample consists of 24 galaxies in Chandra Deep Field South, a unique field for the study of galaxy evolution (12 observed with GNIRS/GEMINI and 12 with ISAAC/VLT). We use H alpha together with the already merged X-ray, ultraviolet, optical, near and mid-infrared imaging data to obtain estimations of SFRs, metallicities, stellar and dynamical masses, AGN activity, and extinction properties. We have obtained 15 Hα detections, 4 rotation curves, and SFR relationship for 7 galaxies. The metallicities obtained for 8 galaxies of the sample are compatible with the metallicities of local galaxies.

Exterior Differential Systems and Euler-Lagrange Partial Differential Equations

The book also covers the Second Variation, Euler-Lagrange PDE systems, and higher-order conservation laws.