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.^

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.

Effect of web holes on web crippling strength of cold-formed steel channel sections under end-one-flange loading condition - Part II: Parametric study and proposed design equations

A parametric study of cold-formed steel sections with web openings subjected to web crippling under end-one-flange (EOF) loading condition is undertaken, using finite element analysis, to investigate the effects of web holes and cross-section sizes. The holes are located either centred above the bearing plates or with a horizontal clear distance to the near edge of the bearing plates. It was demonstrated that the main factors influencing the web crippling strength are the ratio of the hole depth to the depth of the web, the ratio of the length of bearing plates to the flat depth of the web and the location of the holes as defined by the distance of the hole from the edge of the bearing plate divided by the flat depth of web. In this study, design recommendations in the form of web crippling strength reduction factor equations are proposed, which are conservative when compared with the experimental and finite element results.

Projected equations of motion approach to hybrid quantum/classical dynamics in dielectric-metal composites

We introduce a hybrid method for dielectric-metal composites that describes the dynamics of the metallic system classically whilst retaining a quantum description of the dielectric. The time-dependent dipole moment of the classical system is mimicked by the introduction of projected equations of motion (PEOM) and the coupling between the two systems is achieved through an effective dipole-dipole interaction. To benchmark this method, we model a test system (semiconducting quantum dot-metal nanoparticle hybrid). We begin by examining the energy absorption rate, showing agreement between the PEOM method and the analytical rotating wave approximation (RWA) solution. We then investigate population inversion and show that the PEOM method provides an accurate model for the interaction under ultrashort pulse excitation where the traditional RWA breaks down.

Application of Volterra Equations to Solve Unit Commitment Problem of Optimised Energy Storage and Generation

Development of reliable methods for optimised energy storage and generation is one of the most imminent challenges in modern power systems. In this paper an adaptive approach to load leveling problem using novel dynamic models based on the Volterra integral equations of the first kind with piecewise continuous kernels. These integral equations efficiently solve such inverse problem taking into account both the time dependent efficiencies and the availability of generation/storage of each energy storage technology. In this analysis a direct numerical method is employed to find the least-cost dispatch of available storages. The proposed collocation type numerical method has second order accuracy and enjoys self-regularization properties, which is associated with confidence levels of system demand. This adaptive approach is suitable for energy storage optimisation in real time. The efficiency of the proposed methodology is demonstrated on the Single Electricity Market of Republic of Ireland and Northern Ireland.

Population inversion between the sup(3)Hsub(4) and the sup(3)Fsub(4) excited states of Tmsup(3+) investigated by means of numerical solutions of the rate equations system in Tmsup(3+)-doped and Tmsup(3+), Hosup(3+)-codoped fluoride glasses

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Energy transfer rates and population inversion of sup(4)Isub(11/2) excited state of Ersup(3+) investigated by means of numerical solutions of the rate equations system in Er:LiYFsub(4) crystal

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Model order reduction for parameterized nonlinear evolution equations

Otto-von-Guericke-Universität Magdeburg, Fakultät für Mathematik, Dissertation, 2016