Dinamikus optimalizálás és a Leontief-modell (Dynamic optimization and the Leontief model)

A tanulmány a variációszámítás gazdasági alkalmazásaiból ismertet hármat. Mindhárom alkalmazás a Leontief-modellen alapszik. Az optimális pályák vizsgálata után arra keressük a választ, hogy az Euler–Lagrange-differenciálegyenlet rendszerrel kapott megoldások valóban optimális megoldásai-e a modelleknek. Arra a következtetésre jut a tanulmány, hogy csak pótlólagos közgazdasági feltételek bevezetésével határozhatók meg az optimális megoldások. Ugyanakkor a megfogalmazott feltételek segítségével az ismertetett modellek egy általánosabb keretbe illeszthetők. A tanulmány végső eredménye az, hogy mind a három modell optimális megoldása a Neumann-sugárnak felel meg. /===/ The study presents three economic applications of variation calculations. All three rely on the Leontief model. After examination of the optimal courses, an answer is sought to whether the solutions to the Euler–Lagrange differential equation system are really opti-mal solutions to the models. The study concludes that the optimal solutions can only be determined by introducing additional economic conditions. At the same time, the models presented can be fitted into a general framework with the help of the conditions outlined. The final conclusion of the study is that the optimal solution of all three models fits into the Neumann band.

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

Theoretical Investigation of Intra- and Inter-cellular Spatiotemporal Calcium Patterns in Microcirculation

Microcirculatory vessels are lined by endothelial cells (ECs) which are surrounded by a single or multiple layer of smooth muscle cells (SMCs). Spontaneous and agonist induced spatiotemporal calcium (Ca2+) events are generated in ECs and SMCs, and regulated by complex bi-directional signaling between the two layers which ultimately determines the vessel tone. The contractile state of microcirculatory vessels is an important factor in the determination of vascular resistance, blood flow and blood pressure. This dissertation presents theoretical insights into some of the important and currently unresolved phenomena in microvascular tone regulation. Compartmental and continuum models of isolated EC and SMC, coupled EC-SMC and a multi-cellular vessel segment with deterministic and stochastic descriptions of the cellular components were developed, and the intra- and inter-cellular spatiotemporal Ca2+ mobilization was examined. Coupled EC-SMC model simulations captured the experimentally observed localized subcellular EC Ca2+ events arising from the opening of EC transient receptor vanilloid 4 (TRPV4) channels and inositol triphosphate receptors (IP3Rs). These localized EC Ca2+ events result in endothelium-derived hyperpolarization (EDH) and Nitric Oxide (NO) production which transmit to the adjacent SMCs to ultimately result in vasodilation. The model examined the effect of heterogeneous distribution of cellular components and channel gating kinetics in determination of the amplitude and spread of the Ca2+ events. The simulations suggested the necessity of co-localization of certain cellular components for modulation of EDH and NO responses. Isolated EC and SMC models captured intracellular Ca2+ wave like activity and predicted the necessity of non-uniform distribution of cellular components for the generation of Ca2+ waves. The simulations also suggested the role of membrane potential dynamics in regulating Ca2+ wave velocity. The multi-cellular vessel segment model examined the underlying mechanisms for the intercellular synchronization of spontaneous oscillatory Ca2+ waves in individual SMC. From local subcellular events to integrated macro-scale behavior at the vessel level, the developed multi-scale models captured basic features of vascular Ca2+ signaling and provide insights for their physiological relevance. The models provide a theoretical framework for assisting investigations on the regulation of vascular tone in health and disease.

Effects of seasonally varying temperature and salinity on the dynamics of sea lice (Lepeophtheirus salmonis)

Sea lice (Lepeophtheirus salmonis) are an economically significant parasite in salmonid aquaculture. They exhibit temperature-dependent development rates and salinity-dependent mortality, which can greatly impact sea lice population dynamics, but no deterministic models have incorporated these seasonal variables. To understand how seasonality affects sea lice population dynamics, I derive a delay differential equation model with temperature and salinity dependence. I find that peak reproductive output in Newfoundland and British Columbia differs by four months. A sensitivity analysis shows sea lice abundance is most sensitive to variation in mean annual water temperature and salinity, whereas it is lease sensitive to infection rate. Additionally, I investigate the effects of production cycle timing on sea lice management and find that optimal production cycle start times are between the 281st and 337th days of the year in Newfoundland. I also demonstrate that adjusting follow-up treatment timing in response to temperature can improve treatment regimes. My results suggest that effective sea lice management requires consideration of local temperature and salinity patterns.

Mathematical Modeling of Perifusion Cell Culture Experiments

In perifusion cell cultures, the culture medium flows continuously through a chamber containing immobilized cells and the effluent is collected at the end. In our main applications, gonadotropin releasing hormone (GnRH) or oxytocin is introduced into the chamber as the input. They stimulate the cells to secrete luteinizing hormone (LH), which is collected in the effluent. To relate the effluent LH concentration to the cellular processes producing it, we develop and analyze a mathematical model consisting of coupled partial differential equations describing the intracellular signaling and the movement of substances in the cell chamber. We analyze three different data sets and give cellular mechanisms that explain the data. Our model indicates that two negative feedback loops, one fast and one slow, are needed to explain the data and we give their biological bases. We demonstrate that different LH outcomes in oxytocin and GnRH stimulations might originate from different receptor dynamics. We analyze the model to understand the influence of parameters, like the rate of the medium flow or the fraction collection time, on the experimental outcomes. We investigate how the rate of binding and dissociation of the input hormone to and from its receptor influence its movement down the chamber. Finally, we formulate and analyze simpler models that allow us to predict the distortion of a square pulse due to hormone-receptor interactions and to estimate parameters using perifusion data. We show that in the limit of high binding and dissociation the square pulse moves as a diffusing Gaussian and in this limit the biological parameters can be estimated.

Essays on Job Search and Labor Market Dynamics

This dissertation consists of three separate essays on job search and labor market dynamics. In the first essay, “The Impact of Labor Market Conditions on Job Creation: Evidence from Firm Level Data”, I study how much changes in labor market conditions reduce employment fluctuations over the business cycle. Changes in labor market conditions make hiring more expensive during expansions and cheaper during recessions, creating counter-cyclical incentives for job creation. I estimate firm level elasticities of labor demand with respect to changes in labor market conditions, considering two margins: changes in labor market tightness and changes in wages. Using employer-employee matched data from Brazil, I find that all firms are more sensitive to changes in wages rather than labor market tightness, and there is substantial heterogeneity in labor demand elasticity across regions. Based on these results, I demonstrate that changes in labor market conditions reduce the variance of employment growth over the business cycle by 20% in a median region, and this effect is equally driven by changes along each margin. Moreover, I show that the magnitude of the effect of labor market conditions on employment growth can be significantly affected by economic policy. In particular, I document that the rapid growth of the national minimum wages in Brazil in 1997-2010 amplified the impact of the change in labor market conditions during local expansions and diminished this impact during local recessions.

In the second essay, “A Framework for Estimating Persistence of Local Labor

Demand Shocks”, I propose a decomposition which allows me to study the persistence of local labor demand shocks. Persistence of labor demand shocks varies across industries, and the incidence of shocks in a region depends on the regional industrial composition. As a result, less diverse regions are more likely to experience deeper shocks, but not necessarily more long lasting shocks. Building on this idea, I propose a decomposition of local labor demand shocks into idiosyncratic location shocks and nationwide industry shocks and estimate the variance and the persistence of these shocks using the Quarterly Census of Employment and Wages (QCEW) in 1990-2013.

In the third essay, “Conditional Choice Probability Estimation of Continuous- Time Job Search Models”, co-authored with Peter Arcidiacono and Arnaud Maurel, we propose a novel, computationally feasible method of estimating non-stationary job search models. Non-stationary job search models arise in many applications, where policy change can be anticipated by the workers. The most prominent example of such policy is the expiration of unemployment benefits. However, estimating these models still poses a considerable computational challenge, because of the need to solve a differential equation numerically at each step of the optimization routine. We overcome this challenge by adopting conditional choice probability methods, widely used in dynamic discrete choice literature, to job search models and show how the hazard rate out of unemployment and the distribution of the accepted wages, which can be estimated in many datasets, can be used to infer the value of unemployment. We demonstrate how to apply our method by analyzing the effect of the unemployment benefit expiration on duration of unemployment using the data from the Survey of Income and Program Participation (SIPP) in 1996-2007.

On thermodynamics in the primary power conversion of oscillating water column wave energy converters

The paper presents an investigation to the thermodynamics of the air flow in the air chamber for the oscillating water column wave energy converters, in which the oscillating water surface in the water column pressurizes or de-pressurises the air in the chamber. To study the thermodynamics and the compressibility of the air in the chamber, a method is developed in this research: the power take-off is replaced with an accepted semi-empirical relationship between the air flow rate and the oscillating water column chamber pressure, and the thermodynamic process is simplified as an isentropic process. This facilitates the use of a direct expression for the work done on the power take-off by the flowing air and the generation of a single differential equation that defines the thermodynamic process occurring inside the air chamber. Solving the differential equation, the chamber pressure can be obtained if the interior water surface motion is known or the chamber volume (thus the interior water surface motion) if the chamber pressure is known. As a result, the effects of the air compressibility can be studied. Examples given in the paper have shown the compressibility, and its effects on the power losses for large oscillating water column devices.

Evaporation maps for ternary non-ideal liquid mixtures

This thesis deals with the evaporation of non-ideal liquid mixtures using a multicomponent mass transfer approach. It develops the concept of evaporation maps as a convenient way of representing the dynamic composition changes of ternary mixtures during an evaporation process. Evaporation maps represent the residual composition of evaporating ternary non-ideal mixtures over the full range of composition, and are analogous to the commonly-used residue curve maps of simple distillation processes. The evaporation process initially considered in this work involves gas-phase limited evaporation from a liquid or wetted-solid surface, over which a gas flows at known conditions. Evaporation may occur into a pure inert gas, or into one pre-loaded with a known fraction of one of the ternary components. To explore multicomponent masstransfer effects, a model is developed that uses an exact solution to the Maxwell-Stefan equations for mass transfer in the gas film, with a lumped approach applied to the liquid phase. Solutions to the evaporation model take the form of trajectories in temperaturecomposition space, which are then projected onto a ternary diagram to form the map. Novel algorithms are developed for computation of pseudo-azeotropes in the evaporating mixture, and for calculation of the multicomponent wet-bulb temperature at a given liquid composition. A numerical continuation method is used to track the bifurcations which occur in the evaporation maps, where the composition of one component of the pre-loaded gas is the bifurcation parameter. The bifurcation diagrams can in principle be used to determine the required gas composition to produce a specific terminal composition in the liquid. A simple homotopy method is developed to track the locations of the various possible pseudo-azeotropes in the mixture. The stability of pseudo-azeotropes in the gas-phase limited case is examined using a linearized analysis of the governing equations. Algorithms for the calculation of separation boundaries in the evaporation maps are developed using an optimization-based method, as well as a method employing eigenvectors derived from the linearized analysis. The flexure of the wet-bulb temperature surface is explored, and it is shown how evaporation trajectories cross ridges and valleys, so that ridges and valleys of the surface do not coincide with separation boundaries. Finally, the assumption of gas-phase limited mass transfer is relaxed, by employing a model that includes diffusion in the liquid phase. A finite-volume method is used to solve the system of partial differential equations that results. The evaporation trajectories for the distributed model reduce to those of the lumped (gas-phase limited) model as the diffusivity in the liquid increases; under the same gas-phase conditions the permissible terminal compositions of the distributed and lumped models are the same.

C-MEMS Based Micro Enzymatic Biofuel Cells

Miniaturized, self-sufficient bioelectronics powered by unconventional micropower may lead to a new generation of implantable, wireless, minimally invasive medical devices, such as pacemakers, defibrillators, drug-delivering pumps, sensor transmitters, and neurostimulators. Studies have shown that micro-enzymatic biofuel cells (EBFCs) are among the most intuitive candidates for in vivo micropower. In the fisrt part of this thesis, the prototype design of an EBFC chip, having 3D intedigitated microelectrode arrays was proposed to obtain an optimum design of 3D microelectrode arrays for carbon microelectromechanical systems (C-MEMS) based EBFCs. A detailed modeling solving partial differential equations (PDEs) by finite element techniques has been developed on the effect of 1) dimensions of microelectrodes, 2) spatial arrangement of 3D microelectrode arrays, 3) geometry of microelectrode on the EBFC performance based on COMSOL Multiphysics. In the second part of this thesis, in order to investigate the performance of an EBFC, behavior of an EBFC chip performance inside an artery has been studied. COMSOL Multiphysics software has also been applied to analyze mass transport for different orientations of an EBFC chip inside a blood artery. Two orientations: horizontal position (HP) and vertical position (VP) have been analyzed. The third part of this thesis has been focused on experimental work towards high performance EBFC. This work has integrated graphene/enzyme onto three-dimensional (3D) micropillar arrays in order to obtain efficient enzyme immobilization, enhanced enzyme loading and facilitate direct electron transfer. The developed 3D graphene/enzyme network based EBFC generated a maximum power density of 136.3 μWcm-2 at 0.59 V, which is almost 7 times of the maximum power density of the bare 3D carbon micropillar arrays based EBFC. To further improve the EBFC performance, reduced graphene oxide (rGO)/carbon nanotubes (CNTs) has been integrated onto 3D mciropillar arrays to further increase EBFC performance in the fourth part of this thesisThe developed rGO/CNTs based EBFC generated twice the maximum power density of rGO based EBFC. Through a comparison of experimental and theoretical results, the cell performance efficiency is noted to be 67%.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Funzioni di Sobolev

In questa tesi cercherò di analizzare le funzioni di Sobolev su R}^{n}, seguendo le trattazioni Measure Theory and Fine Properties of Functions di L.C. Evans e R.F.Gariepy e l'elaborato Functional Analysis, Sobolev Spaces and Partial Differential Equations di H. Brezis. Le funzioni di Sobolev si caratterizzano per essere funzioni con le derivate prime deboli appartenenti a qualche spazio L^{p}. I vari spazi di Sobolev hanno buone proprietà di completezza e compattezza e conseguentemente sono spesso i giusti spazi per le applicazioni di analisi funzionale. Ora, come vedremo, per definizione, l'integrazione per parti è valida per le funzioni di Sobolev. È, invece, meno ovvio che altre regole di calcolo siano allo stesso modo valide. Così, ho inteso chiarire questa questione di carattere generale, con particolare attenzione alle proprietà puntuali delle funzioni di Sobolev. Abbiamo suddiviso il lavoro svolto in cinque capitoli. Il capitolo 1 contiene le definizioni di base necessarie per la trattazione svolta; nel secondo capitolo sono stati derivati vari modi di approssimazione delle funzioni di Sobolev con funzioni lisce e sono state fornite alcune regole di calcolo per tali funzioni. Il capitolo 3 darà un' interpretazione dei valori al bordo delle funzioni di Sobolev utilizzando l'operatore Traccia, mentre il capitolo 4 discute l' estensione su tutto R^{n} di tali funzioni. Proveremo infine le principali disuguaglianze di Sobolev nel Capitolo 5.

Nonlinear Dynamics Near the Cell Membrane of Chara Corallina

The phenomenon of patterned distribution of pH near the cell membrane of the algae Chara corallina upon illumination is well-known. In this paper, we develop a mathematical model, based on the detailed kinetic analysis of proton fluxes across the cell membrane, to explain this phenomenon. The model yields two coupled nonlinear partial differential equations which describe the spatial dynamics of proton concentration changes and transmembrane potential generation. The experimental observation of pH pattern formation, its period and amplitude of oscillation, and also its hysteresis in response to changing illumination, are all reproduced by our model. A comparison of experimental results and predictions of our theory is made. Finally, a mechanism for pattern formation in Chara corallina is proposed.

Applications of a Collaborative Open Object Oriented Learning Environment to Electronics, Computing and Mathematics Education

A lightweight Java application suite has been developed and deployed allowing collaborative learning between students and tutors at remote locations. Students can engage in group activities online and also collaborate with tutors. A generic Java framework has been developed and applied to electronics, computing and mathematics education. The applications are respectively: (a) a digital circuit simulator, which allows students to collaborate in building simple or complex electronic circuits; (b) a Java programming environment where the paradigm is behavioural-based robotics, and (c) a differential equation solver useful in modelling of any complex and nonlinear dynamic system. Each student sees a common shared window on which may be added text or graphical objects and which can then be shared online. A built-in chat room supports collaborative dialogue. Students can work either in collaborative groups or else in teams as directed by the tutor. This paper summarises the technical architecture of the system as well as the pedagogical implications of the suite. A report of student evaluation is also presented distilled from use over a period of twelve months. We intend this suite to facilitate learning between groups at one or many institutions and to facilitate international collaboration. We also intend to use the suite as a tool to research the establishment and behaviour of collaborative learning groups. We shall make our software freely available to interested researchers.

Eigenfunctions and fundamental solutions of the fractional Laplace and Dirac operators: the Riemann-Liouville case

In this paper we study eigenfunctions and fundamental solutions for the three parameter fractional Laplace operator $\Delta_+^{(\alpha,\beta,\gamma)}:= D_{x_0^+}^{1+\alpha} +D_{y_0^+}^{1+\beta} +D_{z_0^+}^{1+\gamma},$ where $(\alpha, \beta, \gamma) \in \,]0,1]^3$, and the fractional derivatives $D_{x_0^+}^{1+\alpha}$, $D_{y_0^+}^{1+\beta}$, $D_{z_0^+}^{1+\gamma}$ are in the Riemann-Liouville sense. Applying operational techniques via two-dimensional Laplace transform we describe a complete family of eigenfunctions and fundamental solutions of the operator $\Delta_+^{(\alpha,\beta,\gamma)}$ in classes of functions admitting a summable fractional derivative. Making use of the Mittag-Leffler function, a symbolic operational form of the solutions is presented. From the obtained family of fundamental solutions we deduce a family of fundamental solutions of the fractional Dirac operator, which factorizes the fractional Laplace operator. We apply also the method of separation of variables to obtain eigenfunctions and fundamental solutions.

A Discretization of the Hadamard fractional derivative

We present a new discretization for the Hadamard fractional derivative, that simplifies the computations. We then apply the method to solve a fractional differential equation and a fractional variational problem with dependence on the Hadamard fractional derivative.