PARALLEL SOLUTION METHODS OF PARTIAL DIFFERENTIAL EQUATIONS

Directory of Open Access Journals (Sweden)

Korhan KARABULUT

1998-03-01

Full Text Available Partial differential equations arise in almost all fields of science and engineering. Computer time spent in solving partial differential equations is much more than that of in any other problem class. For this reason, partial differential equations are suitable to be solved on parallel computers that offer great computation power. In this study, parallel solution to partial differential equations with Jacobi, Gauss-Siedel, SOR (Succesive OverRelaxation and SSOR (Symmetric SOR algorithms is studied.

Partial differential equation models in the socio-economic sciences.

Science.gov (United States)

Burger, Martin; Caffarelli, Luis; Markowich, Peter A

2014-11-13

Mathematical models based on partial differential equations (PDEs) have become an integral part of quantitative analysis in most branches of science and engineering, recently expanding also towards biomedicine and socio-economic sciences. The application of PDEs in the latter is a promising field, but widely quite open and leading to a variety of novel mathematical challenges. In this introductory article of the Theme Issue, we will provide an overview of the field and its recent boosting topics. Moreover, we will put the contributions to the Theme Issue in an appropriate perspective. © 2014 The Author(s) Published by the Royal Society. All rights reserved.

Beginning partial differential equations

CERN Document Server

O'Neil, Peter V

2014-01-01

A broad introduction to PDEs with an emphasis on specialized topics and applications occurring in a variety of fields Featuring a thoroughly revised presentation of topics, Beginning Partial Differential Equations, Third Edition provides a challenging, yet accessible,combination of techniques, applications, and introductory theory on the subjectof partial differential equations. The new edition offers nonstandard coverageon material including Burger's equation, the telegraph equation, damped wavemotion, and the use of characteristics to solve nonhomogeneous problems. The Third Edition is or

Hamiltonian partial differential equations and applications

CERN Document Server

Nicholls, David; Sulem, Catherine

2015-01-01

This book is a unique selection of work by world-class experts exploring the latest developments in Hamiltonian partial differential equations and their applications. Topics covered within are representative of the field’s wide scope, including KAM and normal form theories, perturbation and variational methods, integrable systems, stability of nonlinear solutions as well as applications to cosmology, fluid mechanics and water waves. The volume contains both surveys and original research papers and gives a concise overview of the above topics, with results ranging from mathematical modeling to rigorous analysis and numerical simulation. It will be of particular interest to graduate students as well as researchers in mathematics and physics, who wish to learn more about the powerful and elegant analytical techniques for Hamiltonian partial differential equations.

Abstract methods in partial differential equations

CERN Document Server

Carroll, Robert W

2012-01-01

Detailed, self-contained treatment examines modern abstract methods in partial differential equations, especially abstract evolution equations. Suitable for graduate students with some previous exposure to classical partial differential equations. 1969 edition.

Partial differential equations of mathematical physics

CERN Document Server

Sobolev, S L

1964-01-01

Partial Differential Equations of Mathematical Physics emphasizes the study of second-order partial differential equations of mathematical physics, which is deemed as the foundation of investigations into waves, heat conduction, hydrodynamics, and other physical problems. The book discusses in detail a wide spectrum of topics related to partial differential equations, such as the theories of sets and of Lebesgue integration, integral equations, Green's function, and the proof of the Fourier method. Theoretical physicists, experimental physicists, mathematicians engaged in pure and applied math

Numerical Analysis of Partial Differential Equations

CERN Document Server

Lui, S H

2011-01-01

A balanced guide to the essential techniques for solving elliptic partial differential equations Numerical Analysis of Partial Differential Equations provides a comprehensive, self-contained treatment of the quantitative methods used to solve elliptic partial differential equations (PDEs), with a focus on the efficiency as well as the error of the presented methods. The author utilizes coverage of theoretical PDEs, along with the nu merical solution of linear systems and various examples and exercises, to supply readers with an introduction to the essential concepts in the numerical analysis

Parameter Estimation of Partial Differential Equation Models

KAUST Repository

Xun, Xiaolei

2013-09-01

Partial differential equation (PDE) models are commonly used to model complex dynamic systems in applied sciences such as biology and finance. The forms of these PDE models are usually proposed by experts based on their prior knowledge and understanding of the dynamic system. Parameters in PDE models often have interesting scientific interpretations, but their values are often unknown and need to be estimated from the measurements of the dynamic system in the presence of measurement errors. Most PDEs used in practice have no analytic solutions, and can only be solved with numerical methods. Currently, methods for estimating PDE parameters require repeatedly solving PDEs numerically under thousands of candidate parameter values, and thus the computational load is high. In this article, we propose two methods to estimate parameters in PDE models: a parameter cascading method and a Bayesian approach. In both methods, the underlying dynamic process modeled with the PDE model is represented via basis function expansion. For the parameter cascading method, we develop two nested levels of optimization to estimate the PDE parameters. For the Bayesian method, we develop a joint model for data and the PDE and develop a novel hierarchical model allowing us to employ Markov chain Monte Carlo (MCMC) techniques to make posterior inference. Simulation studies show that the Bayesian method and parameter cascading method are comparable, and both outperform other available methods in terms of estimation accuracy. The two methods are demonstrated by estimating parameters in a PDE model from long-range infrared light detection and ranging data. Supplementary materials for this article are available online. © 2013 American Statistical Association.

Elements of partial differential equations

CERN Document Server

Sneddon, Ian Naismith

1957-01-01

Geared toward students of applied rather than pure mathematics, this volume introduces elements of partial differential equations. Its focus is primarily upon finding solutions to particular equations rather than general theory.Topics include ordinary differential equations in more than two variables, partial differential equations of the first and second orders, Laplace's equation, the wave equation, and the diffusion equation. A helpful Appendix offers information on systems of surfaces, and solutions to the odd-numbered problems appear at the end of the book. Readers pursuing independent st

Partial differential equations for scientists and engineers

CERN Document Server

Farlow, Stanley J

1993-01-01

Most physical phenomena, whether in the domain of fluid dynamics, electricity, magnetism, mechanics, optics, or heat flow, can be described in general by partial differential equations. Indeed, such equations are crucial to mathematical physics. Although simplifications can be made that reduce these equations to ordinary differential equations, nevertheless the complete description of physical systems resides in the general area of partial differential equations.This highly useful text shows the reader how to formulate a partial differential equation from the physical problem (constructing th

Modeling Solution of Nonlinear Dispersive Partial Differential Equations using the Marker Method

International Nuclear Information System (INIS)

Lewandowski, Jerome L.V.

2005-01-01

A new method for the solution of nonlinear dispersive partial differential equations is described. The marker method relies on the definition of a convective field associated with the underlying partial differential equation; the information about the approximate solution is associated with the response of an ensemble of markers to this convective field. Some key aspects of the method, such as the selection of the shape function and the initial loading, are discussed in some details

Modeling tree crown dynamics with 3D partial differential equations.

Science.gov (United States)

Beyer, Robert; Letort, Véronique; Cournède, Paul-Henry

2014-01-01

We characterize a tree's spatial foliage distribution by the local leaf area density. Considering this spatially continuous variable allows to describe the spatiotemporal evolution of the tree crown by means of 3D partial differential equations. These offer a framework to rigorously take locally and adaptively acting effects into account, notably the growth toward light. Biomass production through photosynthesis and the allocation to foliage and wood are readily included in this model framework. The system of equations stands out due to its inherent dynamic property of self-organization and spontaneous adaptation, generating complex behavior from even only a few parameters. The density-based approach yields spatially structured tree crowns without relying on detailed geometry. We present the methodological fundamentals of such a modeling approach and discuss further prospects and applications.

Application of partial differential equation modeling of the control/structural dynamics of flexible spacecraft

Science.gov (United States)

Taylor, Lawrence W., Jr.; Rajiyah, H.

1991-01-01

Partial differential equations for modeling the structural dynamics and control systems of flexible spacecraft are applied here in order to facilitate systems analysis and optimization of these spacecraft. Example applications are given, including the structural dynamics of SCOLE, the Solar Array Flight Experiment, the Mini-MAST truss, and the LACE satellite. The development of related software is briefly addressed.

Tracer kinetics: Modelling by partial differential equations of inhomogeneous compartments with age-dependent elimination rates. Pt. 2

International Nuclear Information System (INIS)

Winkler, E.

1991-01-01

The general theory of inhomogeneous compartments with age-dependent elimination rates is illustrated by examples. Mathematically, it turns out that models consisting of partial differential equations include ordinary, delayed and integro-differential equations, a general fact which is treated here in the context of linear tracer kinetics. The examples include standard compartments as a degenerate case, systems of standard compartments (compartment blocks), models resulting in special residence time distributions, models with pipes, and systems with heterogeneous particles. (orig./BBR) [de

Basic linear partial differential equations

CERN Document Server

Treves, Francois

1975-01-01

Focusing on the archetypes of linear partial differential equations, this text for upper-level undergraduates and graduate students features most of the basic classical results. The methods, however, are decidedly nontraditional: in practically every instance, they tend toward a high level of abstraction. This approach recalls classical material to contemporary analysts in a language they can understand, as well as exploiting the field's wealth of examples as an introduction to modern theories.The four-part treatment covers the basic examples of linear partial differential equations and their

Stochastic partial differential equations an introduction

CERN Document Server

Liu, Wei

2015-01-01

This book provides an introduction to the theory of stochastic partial differential equations (SPDEs) of evolutionary type. SPDEs are one of the main research directions in probability theory with several wide ranging applications. Many types of dynamics with stochastic influence in nature or man-made complex systems can be modelled by such equations. The theory of SPDEs is based both on the theory of deterministic partial differential equations, as well as on modern stochastic analysis. Whilst this volume mainly follows the ‘variational approach’, it also contains a short account on the ‘semigroup (or mild solution) approach’. In particular, the volume contains a complete presentation of the main existence and uniqueness results in the case of locally monotone coefficients. Various types of generalized coercivity conditions are shown to guarantee non-explosion, but also a systematic approach to treat SPDEs with explosion in finite time is developed. It is, so far, the only book where the latter and t...

First-order partial differential equations

CERN Document Server

Rhee, Hyun-Ku; Amundson, Neal R

2001-01-01

This first volume of a highly regarded two-volume text is fully usable on its own. After going over some of the preliminaries, the authors discuss mathematical models that yield first-order partial differential equations; motivations, classifications, and some methods of solution; linear and semilinear equations; chromatographic equations with finite rate expressions; homogeneous and nonhomogeneous quasilinear equations; formation and propagation of shocks; conservation equations, weak solutions, and shock layers; nonlinear equations; and variational problems. Exercises appear at the end of mo

Taguchi method for partial differential equations with application in tumor growth.

Science.gov (United States)

Ilea, M; Turnea, M; Rotariu, M; Arotăriţei, D; Popescu, Marilena

2014-01-01

The growth of tumors is a highly complex process. To describe this process, mathematical models are needed. A variety of partial differential mathematical models for tumor growth have been developed and studied. Most of those models are based on the reaction-diffusion equations and mass conservation law. A variety of modeling strategies have been developed, each focusing on tumor growth. Systems of time-dependent partial differential equations occur in many branches of applied mathematics. The vast majority of mathematical models in tumor growth are formulated in terms of partial differential equations. We propose a mathematical model for the interactions between these three cancer cell populations. The Taguchi methods are widely used by quality engineering scientists to compare the effects of multiple variables, together with their interactions, with a simple and manageable experimental design. In Taguchi's design of experiments, variation is more interesting to study than the average. First, Taguchi methods are utilized to search for the significant factors and the optimal level combination of parameters. Except the three parameters levels, other factors levels other factors levels would not be considered. Second, cutting parameters namely, cutting speed, depth of cut, and feed rate are designed using the Taguchi method. Finally, the adequacy of the developed mathematical model is proved by ANOVA. According to the results of ANOVA, since the percentage contribution of the combined error is as small. Many mathematical models can be quantitatively characterized by partial differential equations. The use of MATLAB and Taguchi method in this article illustrates the important role of informatics in research in mathematical modeling. The study of tumor growth cells is an exciting and important topic in cancer research and will profit considerably from theoretical input. Interpret these results to be a permanent collaboration between math's and medical oncologists.

Observability of discretized partial differential equations

Science.gov (United States)

Cohn, Stephen E.; Dee, Dick P.

1988-01-01

It is shown that complete observability of the discrete model used to assimilate data from a linear partial differential equation (PDE) system is necessary and sufficient for asymptotic stability of the data assimilation process. The observability theory for discrete systems is reviewed and applied to obtain simple observability tests for discretized constant-coefficient PDEs. Examples are used to show how numerical dispersion can result in discrete dynamics with multiple eigenvalues, thereby detracting from observability.

Solving Partial Differential Equations Using a New Differential Evolution Algorithm

Directory of Open Access Journals (Sweden)

Natee Panagant

2014-01-01

Full Text Available This paper proposes an alternative meshless approach to solve partial differential equations (PDEs. With a global approximate function being defined, a partial differential equation problem is converted into an optimisation problem with equality constraints from PDE boundary conditions. An evolutionary algorithm (EA is employed to search for the optimum solution. For this approach, the most difficult task is the low convergence rate of EA which consequently results in poor PDE solution approximation. However, its attractiveness remains due to the nature of a soft computing technique in EA. The algorithm can be used to tackle almost any kind of optimisation problem with simple evolutionary operation, which means it is mathematically simpler to use. A new efficient differential evolution (DE is presented and used to solve a number of the partial differential equations. The results obtained are illustrated and compared with exact solutions. It is shown that the proposed method has a potential to be a future meshless tool provided that the search performance of EA is greatly enhanced.

Reduced differential transform method for partial differential equations within local fractional derivative operators

Directory of Open Access Journals (Sweden)

Hossein Jafari

2016-04-01

Full Text Available The non-differentiable solution of the linear and non-linear partial differential equations on Cantor sets is implemented in this article. The reduced differential transform method is considered in the local fractional operator sense. The four illustrative examples are given to show the efficiency and accuracy features of the presented technique to solve local fractional partial differential equations.

Reduced-order modelling of parameter-dependent, linear and nonlinear dynamic partial differential equation models.

Science.gov (United States)

Shah, A A; Xing, W W; Triantafyllidis, V

2017-04-01

In this paper, we develop reduced-order models for dynamic, parameter-dependent, linear and nonlinear partial differential equations using proper orthogonal decomposition (POD). The main challenges are to accurately and efficiently approximate the POD bases for new parameter values and, in the case of nonlinear problems, to efficiently handle the nonlinear terms. We use a Bayesian nonlinear regression approach to learn the snapshots of the solutions and the nonlinearities for new parameter values. Computational efficiency is ensured by using manifold learning to perform the emulation in a low-dimensional space. The accuracy of the method is demonstrated on a linear and a nonlinear example, with comparisons with a global basis approach.

Partial differential equations & boundary value problems with Maple

CERN Document Server

Articolo, George A

2009-01-01

Partial Differential Equations and Boundary Value Problems with Maple presents all of the material normally covered in a standard course on partial differential equations, while focusing on the natural union between this material and the powerful computational software, Maple. The Maple commands are so intuitive and easy to learn, students can learn what they need to know about the software in a matter of hours- an investment that provides substantial returns. Maple''s animation capabilities allow students and practitioners to see real-time displays of the solutions of partial differential equations.  Maple files can be found on the books website. Ancillary list: Maple files- http://www.elsevierdirect.com/companion.jsp?ISBN=9780123747327  Provides a quick overview of the software w/simple commands needed to get startedIncludes review material on linear algebra and Ordinary Differential equations, and their contribution in solving partial differential equationsIncorporates an early introduction to Sturm-L...

A partial differential equation model and its reduction to an ordinary differential equation model for prostate tumor growth under intermittent hormone therapy.

Science.gov (United States)

Tao, Youshan; Guo, Qian; Aihara, Kazuyuki

2014-10-01

Hormonal therapy with androgen suppression is a common treatment for advanced prostate tumors. The emergence of androgen-independent cells, however, leads to a tumor relapse under a condition of long-term androgen deprivation. Clinical trials suggest that intermittent androgen suppression (IAS) with alternating on- and off-treatment periods can delay the relapse when compared with continuous androgen suppression (CAS). In this paper, we propose a mathematical model for prostate tumor growth under IAS therapy. The model elucidates initial hormone sensitivity, an eventual relapse of a tumor under CAS therapy, and a delay of a relapse under IAS therapy, which are due to the coexistence of androgen-dependent cells, androgen-independent cells resulting from reversible changes by adaptation, and androgen-independent cells resulting from irreversible changes by genetic mutations. The model is formulated as a free boundary problem of partial differential equations that describe the evolution of populations of the abovementioned three types of cells during on-treatment periods and off-treatment periods. Moreover, the model can be transformed into a piecewise linear ordinary differential equation model by introducing three new volume variables, and the study of the resulting model may help to devise optimal IAS schedules.

Partial differential equations in action from modelling to theory

CERN Document Server

Salsa, Sandro

2016-01-01

The book is intended as an advanced undergraduate or first-year graduate course for students from various disciplines, including applied mathematics, physics and engineering. It has evolved from courses offered on partial differential equations (PDEs) over the last several years at the Politecnico di Milano. These courses had a twofold purpose: on the one hand, to teach students to appreciate the interplay between theory and modeling in problems arising in the applied sciences, and on the other to provide them with a solid theoretical background in numerical methods, such as finite elements. Accordingly, this textbook is divided into two parts. The first part, chapters 2 to 5, is more elementary in nature and focuses on developing and studying basic problems from the macro-areas of diffusion, propagation and transport, waves and vibrations. In turn the second part, chapters 6 to 11, concentrates on the development of Hilbert spaces methods for the variational formulation and the analysis of (mainly) linear bo...

Partial differential equations in action from modelling to theory

CERN Document Server

Salsa, Sandro

2015-01-01

The book is intended as an advanced undergraduate or first-year graduate course for students from various disciplines, including applied mathematics, physics and engineering. It has evolved from courses offered on partial differential equations (PDEs) over the last several years at the Politecnico di Milano. These courses had a twofold purpose: on the one hand, to teach students to appreciate the interplay between theory and modeling in problems arising in the applied sciences, and on the other to provide them with a solid theoretical background in numerical methods, such as finite elements. Accordingly, this textbook is divided into two parts. The first part, chapters 2 to 5, is more elementary in nature and focuses on developing and studying basic problems from the macro-areas of diffusion, propagation and transport, waves and vibrations. In turn the second part, chapters 6 to 11, concentrates on the development of Hilbert spaces methods for the variational formulation and the analysis of (mainly) linear bo...

A lattice Boltzmann model with an amending function for simulating nonlinear partial differential equations

International Nuclear Information System (INIS)

Lin-Jie, Chen; Chang-Feng, Ma

2010-01-01

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations (NPDEs) in the form u t + αuu x + βu n u x + γu xx + δu xxx + ζu xxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman–Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions. (general)

Parent Ratings of ADHD Symptoms: Generalized Partial Credit Model Analysis of Differential Item Functioning across Gender

Science.gov (United States)

Gomez, Rapson

2012-01-01

Objective: Generalized partial credit model, which is based on item response theory (IRT), was used to test differential item functioning (DIF) for the "Diagnostic and Statistical Manual of Mental Disorders" (4th ed.), inattention (IA), and hyperactivity/impulsivity (HI) symptoms across boys and girls. Method: To accomplish this, parents completed…

An odor interaction model of binary odorant mixtures by a partial differential equation method.

Science.gov (United States)

Yan, Luchun; Liu, Jiemin; Wang, Guihua; Wu, Chuandong

2014-07-09

A novel odor interaction model was proposed for binary mixtures of benzene and substituted benzenes by a partial differential equation (PDE) method. Based on the measurement method (tangent-intercept method) of partial molar volume, original parameters of corresponding formulas were reasonably displaced by perceptual measures. By these substitutions, it was possible to relate a mixture's odor intensity to the individual odorant's relative odor activity value (OAV). Several binary mixtures of benzene and substituted benzenes were respectively tested to establish the PDE models. The obtained results showed that the PDE model provided an easily interpretable method relating individual components to their joint odor intensity. Besides, both predictive performance and feasibility of the PDE model were proved well through a series of odor intensity matching tests. If combining the PDE model with portable gas detectors or on-line monitoring systems, olfactory evaluation of odor intensity will be achieved by instruments instead of odor assessors. Many disadvantages (e.g., expense on a fixed number of odor assessors) also will be successfully avoided. Thus, the PDE model is predicted to be helpful to the monitoring and management of odor pollutions.

An Odor Interaction Model of Binary Odorant Mixtures by a Partial Differential Equation Method

Directory of Open Access Journals (Sweden)

Luchun Yan

2014-07-01

Full Text Available A novel odor interaction model was proposed for binary mixtures of benzene and substituted benzenes by a partial differential equation (PDE method. Based on the measurement method (tangent-intercept method of partial molar volume, original parameters of corresponding formulas were reasonably displaced by perceptual measures. By these substitutions, it was possible to relate a mixture’s odor intensity to the individual odorant’s relative odor activity value (OAV. Several binary mixtures of benzene and substituted benzenes were respectively tested to establish the PDE models. The obtained results showed that the PDE model provided an easily interpretable method relating individual components to their joint odor intensity. Besides, both predictive performance and feasibility of the PDE model were proved well through a series of odor intensity matching tests. If combining the PDE model with portable gas detectors or on-line monitoring systems, olfactory evaluation of odor intensity will be achieved by instruments instead of odor assessors. Many disadvantages (e.g., expense on a fixed number of odor assessors also will be successfully avoided. Thus, the PDE model is predicted to be helpful to the monitoring and management of odor pollutions.

Partial differential operators of elliptic type

CERN Document Server

Shimakura, Norio

1992-01-01

This book, which originally appeared in Japanese, was written for use in an undergraduate course or first year graduate course in partial differential equations and is likely to be of interest to researchers as well. This book presents a comprehensive study of the theory of elliptic partial differential operators. Beginning with the definitions of ellipticity for higher order operators, Shimakura discusses the Laplacian in Euclidean spaces, elementary solutions, smoothness of solutions, Vishik-Sobolev problems, the Schauder theory, and degenerate elliptic operators. The appendix covers such preliminaries as ordinary differential equations, Sobolev spaces, and maximum principles. Because elliptic operators arise in many areas, readers will appreciate this book for the way it brings together a variety of techniques that have arisen in different branches of mathematics.

Application of Stochastic Partial Differential Equations to Reservoir Property Modelling

KAUST Repository

Potsepaev, R.

2010-09-06

Existing algorithms of geostatistics for stochastic modelling of reservoir parameters require a mapping (the \\'uvt-transform\\') into the parametric space and reconstruction of a stratigraphic co-ordinate system. The parametric space can be considered to represent a pre-deformed and pre-faulted depositional environment. Existing approximations of this mapping in many cases cause significant distortions to the correlation distances. In this work we propose a coordinate free approach for modelling stochastic textures through the application of stochastic partial differential equations. By avoiding the construction of a uvt-transform and stratigraphic coordinates, one can generate realizations directly in the physical space in the presence of deformations and faults. In particular the solution of the modified Helmholtz equation driven by Gaussian white noise is a zero mean Gaussian stationary random field with exponential correlation function (in 3-D). This equation can be used to generate realizations in parametric space. In order to sample in physical space we introduce a stochastic elliptic PDE with tensor coefficients, where the tensor is related to correlation anisotropy and its variation is physical space.

Methods for partial differential equations qualitative properties of solutions, phase space analysis, semilinear models

CERN Document Server

Ebert, Marcelo R

2018-01-01

This book provides an overview of different topics related to the theory of partial differential equations. Selected exercises are included at the end of each chapter to prepare readers for the “research project for beginners” proposed at the end of the book. It is a valuable resource for advanced graduates and undergraduate students who are interested in specializing in this area. The book is organized in five parts: In Part 1 the authors review the basics and the mathematical prerequisites, presenting two of the most fundamental results in the theory of partial differential equations: the Cauchy-Kovalevskaja theorem and Holmgren's uniqueness theorem in its classical and abstract form. It also introduces the method of characteristics in detail and applies this method to the study of Burger's equation. Part 2 focuses on qualitative properties of solutions to basic partial differential equations, explaining the usual properties of solutions to elliptic, parabolic and hyperbolic equations for the archetypes...

System Entropy Measurement of Stochastic Partial Differential Systems

Directory of Open Access Journals (Sweden)

Bor-Sen Chen

2016-03-01

Full Text Available System entropy describes the dispersal of a system’s energy and is an indication of the disorder of a physical system. Several system entropy measurement methods have been developed for dynamic systems. However, most real physical systems are always modeled using stochastic partial differential dynamic equations in the spatio-temporal domain. No efficient method currently exists that can calculate the system entropy of stochastic partial differential systems (SPDSs in consideration of the effects of intrinsic random fluctuation and compartment diffusion. In this study, a novel indirect measurement method is proposed for calculating of system entropy of SPDSs using a Hamilton–Jacobi integral inequality (HJII-constrained optimization method. In other words, we solve a nonlinear HJII-constrained optimization problem for measuring the system entropy of nonlinear stochastic partial differential systems (NSPDSs. To simplify the system entropy measurement of NSPDSs, the global linearization technique and finite difference scheme were employed to approximate the nonlinear stochastic spatial state space system. This allows the nonlinear HJII-constrained optimization problem for the system entropy measurement to be transformed to an equivalent linear matrix inequalities (LMIs-constrained optimization problem, which can be easily solved using the MATLAB LMI-toolbox (MATLAB R2014a, version 8.3. Finally, several examples are presented to illustrate the system entropy measurement of SPDSs.

A Line-Tau Collocation Method for Partial Differential Equations ...

African Journals Online (AJOL)

This paper deals with the numerical solution of second order linear partial differential equations with the use of the method of lines coupled with the tau collocation method. The method of lines is used to convert the partial differential equation (PDE) to a sequence of ordinary differential equations (ODEs) which is then ...

Numerical solution of two-dimensional non-linear partial differential ...

African Journals Online (AJOL)

linear partial differential equations using a hybrid method. The solution technique involves discritizing the non-linear system of partial differential equations (PDEs) to obtain a corresponding nonlinear system of algebraic difference equations to be ...

Hyperbolic partial differential equations populations, reactors, tides and waves theory and applications

CERN Document Server

Witten, Matthew

1983-01-01

Hyperbolic Partial Differential Equations, Volume 1: Population, Reactors, Tides and Waves: Theory and Applications covers three general areas of hyperbolic partial differential equation applications. These areas include problems related to the McKendrick/Von Foerster population equations, other hyperbolic form equations, and the numerical solution.This text is composed of 15 chapters and begins with surveys of age specific population interactions, populations models of diffusion, nonlinear age dependent population growth with harvesting, local and global stability for the nonlinear renewal eq

CIME course on Control of Partial Differential Equations

CERN Document Server

Alabau-Boussouira, Fatiha; Glass, Olivier; Le Rousseau, Jérôme; Zuazua, Enrique

2012-01-01

The term “control theory” refers to the body of results - theoretical, numerical and algorithmic - which have been developed to influence the evolution of the state of a given system in order to meet a prescribed performance criterion. Systems of interest to control theory may be of very different natures. This monograph is concerned with models that can be described by partial differential equations of evolution. It contains five major contributions and is connected to the CIME Course on Control of Partial Differential Equations that took place in Cetraro (CS, Italy), July 19 - 23, 2010.  Specifically, it covers the stabilization of evolution equations, control of the Liouville equation, control in fluid mechanics, control and numerics for the wave equation, and Carleman estimates for elliptic and parabolic equations with application to control. We are confident this work will provide an authoritative reference work for all scientists who are interested in this field, representing at the same time a fri...

Ambit processes and stochastic partial differential equations

DEFF Research Database (Denmark)

Barndorff-Nielsen, Ole; Benth, Fred Espen; Veraart, Almut

Ambit processes are general stochastic processes based on stochastic integrals with respect to Lévy bases. Due to their flexible structure, they have great potential for providing realistic models for various applications such as in turbulence and finance. This papers studies the connection betwe...... ambit processes and solutions to stochastic partial differential equations. We investigate this relationship from two angles: from the Walsh theory of martingale measures and from the viewpoint of the Lévy noise analysis....

In silico ordinary differential equation/partial differential equation hemodialysis model estimates methadone removal during dialysis

Science.gov (United States)

Linares, Oscar A; Schiesser, William E; Fudin, Jeffrey; Pham, Thien C; Bettinger, Jeffrey J; Mathew, Roy O; Daly, Annemarie L

2015-01-01

Background There is a need to have a model to study methadone’s losses during hemodialysis to provide informed methadone dose recommendations for the practitioner. Aim To build a one-dimensional (1-D), hollow-fiber geometry, ordinary differential equation (ODE) and partial differential equation (PDE) countercurrent hemodialyzer model (ODE/PDE model). Methodology We conducted a cross-sectional study in silico that evaluated eleven hemodialysis patients. Patients received a ceiling dose of methadone hydrochloride 30 mg/day. Outcome measures included: the total amount of methadone removed during dialysis; methadone’s overall intradialytic mass transfer rate coefficient, km; and, methadone’s removal rate, jME. Each metric was measured at dialysate flow rates of 250 mL/min and 800 mL/min. Results The ODE/PDE model revealed a significant increase in the change of methadone’s mass transfer with increased dialysate flow rate, %Δkm=18.56, P=0.02, N=11. The total amount of methadone mass transferred across the dialyzer membrane with high dialysate flow rate significantly increased (0.042±0.016 versus 0.052±0.019 mg/kg, P=0.02, N=11). This was accompanied by a small significant increase in methadone’s mass transfer rate (0.113±0.002 versus 0.014±0.002 mg/kg/h, P=0.02, N=11). The ODE/PDE model accurately predicted methadone’s removal during dialysis. The absolute value of the prediction errors for methadone’s extraction and throughput were less than 2%. Conclusion ODE/PDE modeling of methadone’s hemodialysis is a new approach to study methadone’s removal, in particular, and opioid removal, in general, in patients with end-stage renal disease on hemodialysis. ODE/PDE modeling accurately quantified the fundamental phenomena of methadone’s mass transfer during hemodialysis. This methodology may lead to development of optimally designed intradialytic opioid treatment protocols, and allow dynamic monitoring of outflow plasma opioid concentrations for model

In silico ordinary differential equation/partial differential equation hemodialysis model estimates methadone removal during dialysis.

Science.gov (United States)

Linares, Oscar A; Schiesser, William E; Fudin, Jeffrey; Pham, Thien C; Bettinger, Jeffrey J; Mathew, Roy O; Daly, Annemarie L

2015-01-01

There is a need to have a model to study methadone's losses during hemodialysis to provide informed methadone dose recommendations for the practitioner. To build a one-dimensional (1-D), hollow-fiber geometry, ordinary differential equation (ODE) and partial differential equation (PDE) countercurrent hemodialyzer model (ODE/PDE model). We conducted a cross-sectional study in silico that evaluated eleven hemodialysis patients. Patients received a ceiling dose of methadone hydrochloride 30 mg/day. Outcome measures included: the total amount of methadone removed during dialysis; methadone's overall intradialytic mass transfer rate coefficient, km ; and, methadone's removal rate, j ME. Each metric was measured at dialysate flow rates of 250 mL/min and 800 mL/min. The ODE/PDE model revealed a significant increase in the change of methadone's mass transfer with increased dialysate flow rate, %Δkm =18.56, P=0.02, N=11. The total amount of methadone mass transferred across the dialyzer membrane with high dialysate flow rate significantly increased (0.042±0.016 versus 0.052±0.019 mg/kg, P=0.02, N=11). This was accompanied by a small significant increase in methadone's mass transfer rate (0.113±0.002 versus 0.014±0.002 mg/kg/h, P=0.02, N=11). The ODE/PDE model accurately predicted methadone's removal during dialysis. The absolute value of the prediction errors for methadone's extraction and throughput were less than 2%. ODE/PDE modeling of methadone's hemodialysis is a new approach to study methadone's removal, in particular, and opioid removal, in general, in patients with end-stage renal disease on hemodialysis. ODE/PDE modeling accurately quantified the fundamental phenomena of methadone's mass transfer during hemodialysis. This methodology may lead to development of optimally designed intradialytic opioid treatment protocols, and allow dynamic monitoring of outflow plasma opioid concentrations for model predictive control during dialysis in humans.

Particle Systems and Partial Differential Equations I

CERN Document Server

Gonçalves, Patricia

2014-01-01

This book presents the proceedings of the international conference Particle Systems and Partial Differential Equations I, which took place at the Centre of Mathematics of the University of Minho, Braga, Portugal, from the 5th to the 7th of December, 2012.  The purpose of the conference was to bring together world leaders to discuss their topics of expertise and to present some of their latest research developments in those fields. Among the participants were researchers in probability, partial differential equations and kinetics theory. The aim of the meeting was to present to a varied public the subject of interacting particle systems, its motivation from the viewpoint of physics and its relation with partial differential equations or kinetics theory, and to stimulate discussions and possibly new collaborations among researchers with different backgrounds.  The book contains lecture notes written by François Golse on the derivation of hydrodynamic equations (compressible and incompressible Euler and Navie...

Introduction to partial differential equations with applications

CERN Document Server

Zachmanoglou, E C

1988-01-01

This text explores the essentials of partial differential equations as applied to engineering and the physical sciences. Discusses ordinary differential equations, integral curves and surfaces of vector fields, the Cauchy-Kovalevsky theory, more. Problems and answers.

On the relation between elementary partial difference equations and partial differential equations

NARCIS (Netherlands)

van den Berg, I.P.

1998-01-01

The nonstandard stroboscopy method links discrete-time ordinary difference equations of first-order and continuous-time, ordinary differential equations of first order. We extend this method to the second order, and also to an elementary, yet general class of partial difference/differential

Partial differential equations

CERN Document Server

Agranovich, M S

2002-01-01

Mark Vishik's Partial Differential Equations seminar held at Moscow State University was one of the world's leading seminars in PDEs for over 40 years. This book celebrates Vishik's eightieth birthday. It comprises new results and survey papers written by many renowned specialists who actively participated over the years in Vishik's seminars. Contributions include original developments and methods in PDEs and related fields, such as mathematical physics, tomography, and symplectic geometry. Papers discuss linear and nonlinear equations, particularly linear elliptic problems in angles and gener

Dynamics of partial differential equations

CERN Document Server

Wayne, C Eugene

2015-01-01

This book contains two review articles on the dynamics of partial differential equations that deal with closely related topics but can be read independently. Wayne reviews recent results on the global dynamics of the two-dimensional Navier-Stokes equations. This system exhibits stable vortex solutions: the topic of Wayne's contribution is how solutions that start from arbitrary initial conditions evolve towards stable vortices. Weinstein considers the dynamics of localized states in nonlinear Schrodinger and Gross-Pitaevskii equations that describe many optical and quantum systems. In this contribution, Weinstein reviews recent bifurcations results of solitary waves, their linear and nonlinear stability properties, and results about radiation damping where waves lose energy through radiation.   The articles, written independently, are combined into one volume to showcase the tools of dynamical systems theory at work in explaining qualitative phenomena associated with two classes of partial differential equ...

Generalized solutions of nonlinear partial differential equations

CERN Document Server

Rosinger, EE

1987-01-01

During the last few years, several fairly systematic nonlinear theories of generalized solutions of rather arbitrary nonlinear partial differential equations have emerged. The aim of this volume is to offer the reader a sufficiently detailed introduction to two of these recent nonlinear theories which have so far contributed most to the study of generalized solutions of nonlinear partial differential equations, bringing the reader to the level of ongoing research.The essence of the two nonlinear theories presented in this volume is the observation that much of the mathematics concernin

Controllability of partial differential equations governed by multiplicative controls

CERN Document Server

Khapalov, Alexander Y

2010-01-01

The goal of this monograph is to address the issue of the global controllability of partial differential equations in the context of multiplicative (or bilinear) controls, which enter the model equations as coefficients. The mathematical models we examine include the linear and nonlinear parabolic and hyperbolic PDE's, the Schrödinger equation, and coupled hybrid nonlinear distributed parameter systems modeling the swimming phenomenon. The book offers a new, high-quality and intrinsically nonlinear methodology to approach the aforementioned highly nonlinear controllability problems.

Applied partial differential equations

CERN Document Server

Logan, J David

2004-01-01

This primer on elementary partial differential equations presents the standard material usually covered in a one-semester, undergraduate course on boundary value problems and PDEs. What makes this book unique is that it is a brief treatment, yet it covers all the major ideas: the wave equation, the diffusion equation, the Laplace equation, and the advection equation on bounded and unbounded domains. Methods include eigenfunction expansions, integral transforms, and characteristics. Mathematical ideas are motivated from physical problems, and the exposition is presented in a concise style accessible to science and engineering students; emphasis is on motivation, concepts, methods, and interpretation, rather than formal theory. This second edition contains new and additional exercises, and it includes a new chapter on the applications of PDEs to biology: age structured models, pattern formation; epidemic wave fronts, and advection-diffusion processes. The student who reads through this book and solves many of t...

Modulating Function-Based Method for Parameter and Source Estimation of Partial Differential Equations

KAUST Repository

Asiri, Sharefa M.

2017-01-01

Partial Differential Equations (PDEs) are commonly used to model complex systems that arise for example in biology, engineering, chemistry, and elsewhere. The parameters (or coefficients) and the source of PDE models are often unknown

Partial differential equations

CERN Document Server

Levine, Harold

1997-01-01

The subject matter, partial differential equations (PDEs), has a long history (dating from the 18th century) and an active contemporary phase. An early phase (with a separate focus on taut string vibrations and heat flow through solid bodies) stimulated developments of great importance for mathematical analysis, such as a wider concept of functions and integration and the existence of trigonometric or Fourier series representations. The direct relevance of PDEs to all manner of mathematical, physical and technical problems continues. This book presents a reasonably broad introductory account of the subject, with due regard for analytical detail, applications and historical matters.

Topics in numerical partial differential equations and scientific computing

CERN Document Server

2016-01-01

Numerical partial differential equations (PDEs) are an important part of numerical simulation, the third component of the modern methodology for science and engineering, besides the traditional theory and experiment. This volume contains papers that originated with the collaborative research of the teams that participated in the IMA Workshop for Women in Applied Mathematics: Numerical Partial Differential Equations and Scientific Computing in August 2014.

Semi-bounded partial differential operators

CERN Document Server

Cialdea, Alberto

2014-01-01

This book examines the conditions for the semi-boundedness of partial differential operators, which are interpreted in different ways. For example, today we know a great deal about L2-semibounded differential and pseudodifferential operators, although their complete characterization in analytic terms still poses difficulties, even for fairly simple operators. In contrast, until recently almost nothing was known about analytic characterizations of semi-boundedness for differential operators in other Hilbert function spaces and in Banach function spaces. This book works to address that gap. As such, various types of semi-boundedness are considered and a number of relevant conditions which are either necessary and sufficient or best possible in a certain sense are presented. The majority of the results reported on are the authors’ own contributions.

Differential equation analysis in biomedical science and engineering partial differential equation applications with R

CERN Document Server

Schiesser, William E

2014-01-01

Features a solid foundation of mathematical and computational tools to formulate and solve real-world PDE problems across various fields With a step-by-step approach to solving partial differential equations (PDEs), Differential Equation Analysis in Biomedical Science and Engineering: Partial Differential Equation Applications with R successfully applies computational techniques for solving real-world PDE problems that are found in a variety of fields, including chemistry, physics, biology, and physiology. The book provides readers with the necessary knowledge to reproduce and extend the com

Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations

International Nuclear Information System (INIS)

Rhebergen, S.; Bokhove, O.; Vegt, J.J.W. van der

2008-01-01

We present space- and space-time discontinuous Galerkin finite element (DGFEM) formulations for systems containing nonconservative products, such as occur in dispersed multiphase flow equations. The main criterium we pose on the weak formulation is that if the system of nonconservative partial differential equations can be transformed into conservative form, then the formulation must reduce to that for conservative systems. Standard DGFEM formulations cannot be applied to nonconservative systems of partial differential equations. We therefore introduce the theory of weak solutions for nonconservative products into the DGFEM formulation leading to the new question how to define the path connecting left and right states across a discontinuity. The effect of different paths on the numerical solution is investigated and found to be small. We also introduce a new numerical flux that is able to deal with nonconservative products. Our scheme is applied to two different systems of partial differential equations. First, we consider the shallow water equations, where topography leads to nonconservative products, in which the known, possibly discontinuous, topography is formally taken as an unknown in the system. Second, we consider a simplification of a depth-averaged two-phase flow model which contains more intrinsic nonconservative products

Entropy and convexity for nonlinear partial differential equations.

Science.gov (United States)

Ball, John M; Chen, Gui-Qiang G

2013-12-28

Partial differential equations are ubiquitous in almost all applications of mathematics, where they provide a natural mathematical description of many phenomena involving change in physical, chemical, biological and social processes. The concept of entropy originated in thermodynamics and statistical physics during the nineteenth century to describe the heat exchanges that occur in the thermal processes in a thermodynamic system, while the original notion of convexity is for sets and functions in mathematics. Since then, entropy and convexity have become two of the most important concepts in mathematics. In particular, nonlinear methods via entropy and convexity have been playing an increasingly important role in the analysis of nonlinear partial differential equations in recent decades. This opening article of the Theme Issue is intended to provide an introduction to entropy, convexity and related nonlinear methods for the analysis of nonlinear partial differential equations. We also provide a brief discussion about the content and contributions of the papers that make up this Theme Issue.

Lagrangian vector field and Lagrangian formulation of partial differential equations

Directory of Open Access Journals (Sweden)

M.Chen

2005-01-01

Full Text Available In this paper we consider the Lagrangian formulation of a system of second order quasilinear partial differential equations. Specifically we construct a Lagrangian vector field such that the flows of the vector field satisfy the original system of partial differential equations.

On the hierarchy of partially invariant submodels of differential equations

Energy Technology Data Exchange (ETDEWEB)

Golovin, Sergey V [Lavrentyev Institute of Hydrodynamics SB RAS, Novosibirsk 630090 (Russian Federation)], E-mail: sergey@hydro.nsc.ru

2008-07-04

It is noted that the partially invariant solution (PIS) of differential equations in many cases can be represented as an invariant reduction of some PISs of the higher rank. This introduces a hierarchic structure in the set of all PISs of a given system of differential equations. An equivalence of the two-step and the direct ways of construction of PISs is proved. The hierarchy simplifies the process of enumeration and analysis of partially invariant submodels to the given system of differential equations. In this framework, the complete classification of regular partially invariant solutions of ideal MHD equations is given.

On the hierarchy of partially invariant submodels of differential equations

Science.gov (United States)

Golovin, Sergey V.

2008-07-01

It is noted that the partially invariant solution (PIS) of differential equations in many cases can be represented as an invariant reduction of some PISs of the higher rank. This introduces a hierarchic structure in the set of all PISs of a given system of differential equations. An equivalence of the two-step and the direct ways of construction of PISs is proved. The hierarchy simplifies the process of enumeration and analysis of partially invariant submodels to the given system of differential equations. In this framework, the complete classification of regular partially invariant solutions of ideal MHD equations is given.

On the hierarchy of partially invariant submodels of differential equations

International Nuclear Information System (INIS)

Golovin, Sergey V

2008-01-01

It is noted that the partially invariant solution (PIS) of differential equations in many cases can be represented as an invariant reduction of some PISs of the higher rank. This introduces a hierarchic structure in the set of all PISs of a given system of differential equations. An equivalence of the two-step and the direct ways of construction of PISs is proved. The hierarchy simplifies the process of enumeration and analysis of partially invariant submodels to the given system of differential equations. In this framework, the complete classification of regular partially invariant solutions of ideal MHD equations is given

Derivation of a macroscale formulation for a class of nonlinear partial differential equations

International Nuclear Information System (INIS)

Pantelis, G.

1995-05-01

A macroscale formulation is constructed from a system of partial differential equations which govern the microscale dependent variables. The construction is based upon the requirement that the solutions of the macroscale partial differential equations satisfy, in some approximate sense, the system of partial differential equations associated with the microscale. These results are restricted to the class of nonlinear partial differential equations which can be expressed as polynomials of the dependent variables and their partial derivatives up to second order. A linear approximation of transformations of second order contact manifolds is employed. 6 refs

From stochastic processes to numerical methods: A new scheme for solving reaction subdiffusion fractional partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Angstmann, C.N.; Donnelly, I.C. [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia); Henry, B.I., E-mail: B.Henry@unsw.edu.au [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia); Jacobs, B.A. [School of Computer Science and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag 3, Wits 2050 (South Africa); DST–NRF Centre of Excellence in Mathematical and Statistical Sciences (CoE-MaSS) (South Africa); Langlands, T.A.M. [Department of Mathematics and Computing, University of Southern Queensland, Toowoomba QLD 4350 (Australia); Nichols, J.A. [School of Mathematics and Statistics, UNSW Australia, Sydney NSW 2052 (Australia)

2016-02-15

We have introduced a new explicit numerical method, based on a discrete stochastic process, for solving a class of fractional partial differential equations that model reaction subdiffusion. The scheme is derived from the master equations for the evolution of the probability density of a sum of discrete time random walks. We show that the diffusion limit of the master equations recovers the fractional partial differential equation of interest. This limiting procedure guarantees the consistency of the numerical scheme. The positivity of the solution and stability results are simply obtained, provided that the underlying process is well posed. We also show that the method can be applied to standard reaction–diffusion equations. This work highlights the broader applicability of using discrete stochastic processes to provide numerical schemes for partial differential equations, including fractional partial differential equations.

Nonlocal symmetry generators and explicit solutions of some partial differential equations

International Nuclear Information System (INIS)

Qin Maochang

2007-01-01

The nonlocal symmetry of a partial differential equation is studied in this paper. The partial differential equation written as a conservation law can be transformed into an equivalent system by introducing a suitable potential. The nonlocal symmetry group generators of original partial differential equations can be obtained through their equivalent system. Further, new explicit solutions can be constructed from the newly obtained symmetry generators. The Burgers equation is chosen as an example; many new valuable explicit solutions and nonlocal symmetry generators are presented

A note on the Lie symmetries of complex partial differential

Indian Academy of Sciences (India)

Folklore suggests that the split Lie-like operators of a complex partial differential equation are symmetries of the split system of real partial differential equations. However, this is not the case generally. We illustrate this by using the complex heat equation, wave equation with dissipation, the nonlinear Burgers equation and ...

ERC Workshop on Geometric Partial Differential Equations

CERN Document Server

Novaga, Matteo; Valdinoci, Enrico

2013-01-01

This book is the outcome of a conference held at the Centro De Giorgi of the Scuola Normale of Pisa in September 2012. The aim of the conference was to discuss recent results on nonlinear partial differential equations, and more specifically geometric evolutions and reaction-diffusion equations. Particular attention was paid to self-similar solutions, such as solitons and travelling waves, asymptotic behaviour, formation of singularities and qualitative properties of solutions. These problems arise in many models from Physics, Biology, Image Processing and Applied Mathematics in general, and have attracted a lot of attention in recent years.

Data-driven discovery of partial differential equations.

Science.gov (United States)

Rudy, Samuel H; Brunton, Steven L; Proctor, Joshua L; Kutz, J Nathan

2017-04-01

We propose a sparse regression method capable of discovering the governing partial differential equation(s) of a given system by time series measurements in the spatial domain. The regression framework relies on sparsity-promoting techniques to select the nonlinear and partial derivative terms of the governing equations that most accurately represent the data, bypassing a combinatorially large search through all possible candidate models. The method balances model complexity and regression accuracy by selecting a parsimonious model via Pareto analysis. Time series measurements can be made in an Eulerian framework, where the sensors are fixed spatially, or in a Lagrangian framework, where the sensors move with the dynamics. The method is computationally efficient, robust, and demonstrated to work on a variety of canonical problems spanning a number of scientific domains including Navier-Stokes, the quantum harmonic oscillator, and the diffusion equation. Moreover, the method is capable of disambiguating between potentially nonunique dynamical terms by using multiple time series taken with different initial data. Thus, for a traveling wave, the method can distinguish between a linear wave equation and the Korteweg-de Vries equation, for instance. The method provides a promising new technique for discovering governing equations and physical laws in parameterized spatiotemporal systems, where first-principles derivations are intractable.

Computational partial differential equations using Matlab

CERN Document Server

Li, Jichun

2008-01-01

Brief Overview of Partial Differential Equations The parabolic equations The wave equations The elliptic equations Differential equations in broader areasA quick review of numerical methods for PDEsFinite Difference Methods for Parabolic Equations Introduction Theoretical issues: stability, consistence, and convergence 1-D parabolic equations2-D and 3-D parabolic equationsNumerical examples with MATLAB codesFinite Difference Methods for Hyperbolic Equations IntroductionSome basic difference schemes Dissipation and dispersion errors Extensions to conservation lawsThe second-order hyperbolic PDE

International Conference on Multiscale Methods and Partial Differential Equations.

Energy Technology Data Exchange (ETDEWEB)

Thomas Hou

2006-12-12

The International Conference on Multiscale Methods and Partial Differential Equations (ICMMPDE for short) was held at IPAM, UCLA on August 26-27, 2005. The conference brought together researchers, students and practitioners with interest in the theoretical, computational and practical aspects of multiscale problems and related partial differential equations. The conference provided a forum to exchange and stimulate new ideas from different disciplines, and to formulate new challenging multiscale problems that will have impact in applications.

Algebraic aspects of evolution partial differential equation arising in the study of constant elasticity of variance model from financial mathematics

Science.gov (United States)

Motsepa, Tanki; Aziz, Taha; Fatima, Aeeman; Khalique, Chaudry Masood

2018-03-01

The optimal investment-consumption problem under the constant elasticity of variance (CEV) model is investigated from the perspective of Lie group analysis. The Lie symmetry group of the evolution partial differential equation describing the CEV model is derived. The Lie point symmetries are then used to obtain an exact solution of the governing model satisfying a standard terminal condition. Finally, we construct conservation laws of the underlying equation using the general theorem on conservation laws.

On the hierarchy of partially invariant submodels of differential equations

OpenAIRE

Golovin, Sergey V.

2007-01-01

It is noticed, that partially invariant solution (PIS) of differential equations in many cases can be represented as an invariant reduction of some PIS of the higher rank. This introduce a hierarchic structure in the set of all PISs of a given system of differential equations. By using this structure one can significantly decrease an amount of calculations required in enumeration of all PISs for a given system of partially differential equations. An equivalence of the two-step and the direct ...

Consistency of direct integral estimator for partially observed systems of ordinary differential equations

NARCIS (Netherlands)

Vujačić, Ivan; Dattner, Itai

In this paper we use the sieve framework to prove consistency of the ‘direct integral estimator’ of parameters for partially observed systems of ordinary differential equations, which are commonly used for modeling dynamic processes.

In silico ordinary differential equation/partial differential equation hemodialysis model estimates methadone removal during dialysis

Directory of Open Access Journals (Sweden)

Linares OA

2015-07-01

Full Text Available Oscar A Linares,1 William E Schiesser,2 Jeffrey Fudin,3–6 Thien C Pham,6 Jeffrey J Bettinger,6 Roy O Mathew,6 Annemarie L Daly7 1Translational Genomic Medicine Lab, Plymouth Pharmacokinetic Modeling Study Group, Plymouth, MI, 2Department of Chemical and Biomolecular Engineering, Lehigh University, Bethlehem, PA, 3University of Connecticut School of Pharmacy, Storrs, CT, 4Western New England College of Pharmacy, Springfield, MA, 5Albany College of Pharmacy and Health Sciences, Albany, NY, 6Stratton VA Medical Center, Albany, NY, 7Grace Hospice of Ann Arbor, Ann Arbor, MI, USA Background: There is a need to have a model to study methadone’s losses during hemodialysis to provide informed methadone dose recommendations for the practitioner. Aim: To build a one-dimensional (1-D, hollow-fiber geometry, ordinary differential equation (ODE and partial differential equation (PDE countercurrent hemodialyzer model (ODE/PDE model. Methodology: We conducted a cross-sectional study in silico that evaluated eleven hemodialysis patients. Patients received a ceiling dose of methadone hydrochloride 30 mg/day. Outcome measures included: the total amount of methadone removed during dialysis; methadone’s overall intradialytic mass transfer rate coefficient, km; and, methadone’s removal rate, jME. Each metric was measured at dialysate flow rates of 250 mL/min and 800 mL/min. Results: The ODE/PDE model revealed a significant increase in the change of methadone’s mass transfer with increased dialysate flow rate, %Δ km=18.56, P=0.02, N=11. The total amount of methadone mass transferred across the dialyzer membrane with high dialysate flow rate significantly increased (0.042±0.016 versus 0.052±0.019 mg/kg, P=0.02, N=11. This was accompanied by a small significant increase in methadone’s mass transfer rate (0.113±0.002 versus 0.014±0.002 mg/kg/h, P=0.02, N=11. The ODE/PDE model accurately predicted methadone’s removal during dialysis. The absolute value

3rd International Conference on Particle Systems and Partial Differential Equations

CERN Document Server

Soares, Ana

2016-01-01

The main focus of this book is on different topics in probability theory, partial differential equations and kinetic theory, presenting some of the latest developments in these fields. It addresses mathematical problems concerning applications in physics, engineering, chemistry and biology that were presented at the Third International Conference on Particle Systems and Partial Differential Equations, held at the University of Minho, Braga, Portugal in December 2014. The purpose of the conference was to bring together prominent researchers working in the fields of particle systems and partial differential equations, providing a venue for them to present their latest findings and discuss their areas of expertise. Further, it was intended to introduce a vast and varied public, including young researchers, to the subject of interacting particle systems, its underlying motivation, and its relation to partial differential equations. This book will appeal to probabilists, analysts and those mathematicians whose wor...

Multigrid for high dimensional elliptic partial differential equations on non-equidistant grids

NARCIS (Netherlands)

bin Zubair, H.; Oosterlee, C.E.; Wienands, R.

2006-01-01

This work presents techniques, theory and numbers for multigrid in a general d-dimensional setting. The main focus is the multigrid convergence for high-dimensional partial differential equations (PDEs). As a model problem we have chosen the anisotropic diffusion equation, on a unit hypercube. We

A model and variance reduction method for computing statistical outputs of stochastic elliptic partial differential equations

International Nuclear Information System (INIS)

Vidal-Codina, F.; Nguyen, N.C.; Giles, M.B.; Peraire, J.

2015-01-01

We present a model and variance reduction method for the fast and reliable computation of statistical outputs of stochastic elliptic partial differential equations. Our method consists of three main ingredients: (1) the hybridizable discontinuous Galerkin (HDG) discretization of elliptic partial differential equations (PDEs), which allows us to obtain high-order accurate solutions of the governing PDE; (2) the reduced basis method for a new HDG discretization of the underlying PDE to enable real-time solution of the parameterized PDE in the presence of stochastic parameters; and (3) a multilevel variance reduction method that exploits the statistical correlation among the different reduced basis approximations and the high-fidelity HDG discretization to accelerate the convergence of the Monte Carlo simulations. The multilevel variance reduction method provides efficient computation of the statistical outputs by shifting most of the computational burden from the high-fidelity HDG approximation to the reduced basis approximations. Furthermore, we develop a posteriori error estimates for our approximations of the statistical outputs. Based on these error estimates, we propose an algorithm for optimally choosing both the dimensions of the reduced basis approximations and the sizes of Monte Carlo samples to achieve a given error tolerance. We provide numerical examples to demonstrate the performance of the proposed method

Canonical coordinates for partial differential equations

Science.gov (United States)

Hunt, L. R.; Villarreal, Ramiro

1988-01-01

Necessary and sufficient conditions are found under which operators of the form Sigma (m, j=1) x (2) sub j + X sub O can be made constant coefficient. In addition, necessary and sufficient conditions are derived which classify those linear partial differential operators that can be moved to the Kolmogorov type.

Canonical coordinates for partial differential equations

Science.gov (United States)

Hunt, L. R.; Villarreal, Ramiro

1987-01-01

Necessary and sufficient conditions are found under which operators of the form Sigma(m, j=1) X(2)sub j + X sub 0 can be made constant coefficient. In addition, necessary and sufficient conditions are derived which classify those linear partial differential operators that can be moved to the Kolmogorov type.

Research on odor interaction between aldehyde compounds via a partial differential equation (PDE) model.

Science.gov (United States)

Yan, Luchun; Liu, Jiemin; Qu, Chen; Gu, Xingye; Zhao, Xia

2015-01-28

In order to explore the odor interaction of binary odor mixtures, a series of odor intensity evaluation tests were performed using both individual components and binary mixtures of aldehydes. Based on the linear relation between the logarithm of odor activity value and odor intensity of individual substances, the relationship between concentrations of individual constituents and their joint odor intensity was investigated by employing a partial differential equation (PDE) model. The obtained results showed that the binary odor interaction was mainly influenced by the mixing ratio of two constituents, but not the concentration level of an odor sample. Besides, an extended PDE model was also proposed on the basis of the above experiments. Through a series of odor intensity matching tests for several different binary odor mixtures, the extended PDE model was proved effective at odor intensity prediction. Furthermore, odorants of the same chemical group and similar odor type exhibited similar characteristics in the binary odor interaction. The overall results suggested that the PDE model is a more interpretable way of demonstrating the odor interactions of binary odor mixtures.

Optimal moving grids for time-dependent partial differential equations

Science.gov (United States)

Wathen, A. J.

1992-01-01

Various adaptive moving grid techniques for the numerical solution of time-dependent partial differential equations were proposed. The precise criterion for grid motion varies, but most techniques will attempt to give grids on which the solution of the partial differential equation can be well represented. Moving grids are investigated on which the solutions of the linear heat conduction and viscous Burgers' equation in one space dimension are optimally approximated. Precisely, the results of numerical calculations of optimal moving grids for piecewise linear finite element approximation of PDE solutions in the least-squares norm are reported.

Constructing general partial differential equations using polynomial and neural networks.

Science.gov (United States)

Zjavka, Ladislav; Pedrycz, Witold

2016-01-01

Sum fraction terms can approximate multi-variable functions on the basis of discrete observations, replacing a partial differential equation definition with polynomial elementary data relation descriptions. Artificial neural networks commonly transform the weighted sum of inputs to describe overall similarity relationships of trained and new testing input patterns. Differential polynomial neural networks form a new class of neural networks, which construct and solve an unknown general partial differential equation of a function of interest with selected substitution relative terms using non-linear multi-variable composite polynomials. The layers of the network generate simple and composite relative substitution terms whose convergent series combinations can describe partial dependent derivative changes of the input variables. This regression is based on trained generalized partial derivative data relations, decomposed into a multi-layer polynomial network structure. The sigmoidal function, commonly used as a nonlinear activation of artificial neurons, may transform some polynomial items together with the parameters with the aim to improve the polynomial derivative term series ability to approximate complicated periodic functions, as simple low order polynomials are not able to fully make up for the complete cycles. The similarity analysis facilitates substitutions for differential equations or can form dimensional units from data samples to describe real-world problems. Copyright © 2015 Elsevier Ltd. All rights reserved.

An Accurate Approximate-Analytical Technique for Solving Time-Fractional Partial Differential Equations

Directory of Open Access Journals (Sweden)

M. Bishehniasar

2017-01-01

Full Text Available The demand of many scientific areas for the usage of fractional partial differential equations (FPDEs to explain their real-world systems has been broadly identified. The solutions may portray dynamical behaviors of various particles such as chemicals and cells. The desire of obtaining approximate solutions to treat these equations aims to overcome the mathematical complexity of modeling the relevant phenomena in nature. This research proposes a promising approximate-analytical scheme that is an accurate technique for solving a variety of noninteger partial differential equations (PDEs. The proposed strategy is based on approximating the derivative of fractional-order and reducing the problem to the corresponding partial differential equation (PDE. Afterwards, the approximating PDE is solved by using a separation-variables technique. The method can be simply applied to nonhomogeneous problems and is proficient to diminish the span of computational cost as well as achieving an approximate-analytical solution that is in excellent concurrence with the exact solution of the original problem. In addition and to demonstrate the efficiency of the method, it compares with two finite difference methods including a nonstandard finite difference (NSFD method and standard finite difference (SFD technique, which are popular in the literature for solving engineering problems.

Numerical Analysis for Stochastic Partial Differential Delay Equations with Jumps

OpenAIRE

Li, Yan; Hu, Junhao

2013-01-01

We investigate the convergence rate of Euler-Maruyama method for a class of stochastic partial differential delay equations driven by both Brownian motion and Poisson point processes. We discretize in space by a Galerkin method and in time by using a stochastic exponential integrator. We generalize some results of Bao et al. (2011) and Jacob et al. (2009) in finite dimensions to a class of stochastic partial differential delay equations with jumps in infinite dimensions.

Constructing and predicting solitary pattern solutions for nonlinear time-fractional dispersive partial differential equations

Science.gov (United States)

Arqub, Omar Abu; El-Ajou, Ahmad; Momani, Shaher

2015-07-01

Building fractional mathematical models for specific phenomena and developing numerical or analytical solutions for these fractional mathematical models are crucial issues in mathematics, physics, and engineering. In this work, a new analytical technique for constructing and predicting solitary pattern solutions of time-fractional dispersive partial differential equations is proposed based on the generalized Taylor series formula and residual error function. The new approach provides solutions in the form of a rapidly convergent series with easily computable components using symbolic computation software. For method evaluation and validation, the proposed technique was applied to three different models and compared with some of the well-known methods. The resultant simulations clearly demonstrate the superiority and potentiality of the proposed technique in terms of the quality performance and accuracy of substructure preservation in the construct, as well as the prediction of solitary pattern solutions for time-fractional dispersive partial differential equations.

Compatible Spatial Discretizations for Partial Differential Equations

Energy Technology Data Exchange (ETDEWEB)

Arnold, Douglas, N, ed.

2004-11-25

From May 11--15, 2004, the Institute for Mathematics and its Applications held a hot topics workshop on Compatible Spatial Discretizations for Partial Differential Equations. The numerical solution of partial differential equations (PDE) is a fundamental task in science and engineering. The goal of the workshop was to bring together a spectrum of scientists at the forefront of the research in the numerical solution of PDEs to discuss compatible spatial discretizations. We define compatible spatial discretizations as those that inherit or mimic fundamental properties of the PDE such as topology, conservation, symmetries, and positivity structures and maximum principles. A wide variety of discretization methods applied across a wide range of scientific and engineering applications have been designed to or found to inherit or mimic intrinsic spatial structure and reproduce fundamental properties of the solution of the continuous PDE model at the finite dimensional level. A profusion of such methods and concepts relevant to understanding them have been developed and explored: mixed finite element methods, mimetic finite differences, support operator methods, control volume methods, discrete differential forms, Whitney forms, conservative differencing, discrete Hodge operators, discrete Helmholtz decomposition, finite integration techniques, staggered grid and dual grid methods, etc. This workshop seeks to foster communication among the diverse groups of researchers designing, applying, and studying such methods as well as researchers involved in practical solution of large scale problems that may benefit from advancements in such discretizations; to help elucidate the relations between the different methods and concepts; and to generally advance our understanding in the area of compatible spatial discretization methods for PDE. Particular points of emphasis included: + Identification of intrinsic properties of PDE models that are critical for the fidelity of numerical

Partial differential equations with variable exponents variational methods and qualitative analysis

CERN Document Server

Radulescu, Vicentiu D

2015-01-01

Partial Differential Equations with Variable Exponents: Variational Methods and Qualitative Analysis provides researchers and graduate students with a thorough introduction to the theory of nonlinear partial differential equations (PDEs) with a variable exponent, particularly those of elliptic type. The book presents the most important variational methods for elliptic PDEs described by nonhomogeneous differential operators and containing one or more power-type nonlinearities with a variable exponent. The authors give a systematic treatment of the basic mathematical theory and constructive meth

Preconditioners based on windowed Fourier frames applied to elliptic partial differential equations

NARCIS (Netherlands)

Bhowmik, S.K.; Stolk, C.C.

2011-01-01

We investigate the application of windowed Fourier frames to the numerical solution of partial differential equations, focussing on elliptic equations. The action of a partial differential operator (PDO) on a windowed plane wave is close to a multiplication, where the multiplication factor is given

Darboux transformations and linear parabolic partial differential equations

International Nuclear Information System (INIS)

Arrigo, Daniel J.; Hickling, Fred

2002-01-01

Solutions for a class of linear parabolic partial differential equation are provided. These solutions are obtained by first solving a system of (n+1) nonlinear partial differential equations. This system arises as the coefficients of a Darboux transformation and is equivalent to a matrix Burgers' equation. This matrix equation is solved using a generalized Hopf-Cole transformation. The solutions for the original equation are given in terms of solutions of the heat equation. These results are applied to the (1+1)-dimensional Schroedinger equation where all bound state solutions are obtained for a 2n-parameter family of potentials. As a special case, the solutions for integral members of the regular and modified Poeschl-Teller potentials are recovered. (author). Letter-to-the-editor

Improved stochastic approximation methods for discretized parabolic partial differential equations

Science.gov (United States)

Guiaş, Flavius

2016-12-01

We present improvements of the stochastic direct simulation method, a known numerical scheme based on Markov jump processes which is used for approximating solutions of ordinary differential equations. This scheme is suited especially for spatial discretizations of evolution partial differential equations (PDEs). By exploiting the full path simulation of the stochastic method, we use this first approximation as a predictor and construct improved approximations by Picard iterations, Runge-Kutta steps, or a combination. This has as consequence an increased order of convergence. We illustrate the features of the improved method at a standard benchmark problem, a reaction-diffusion equation modeling a combustion process in one space dimension (1D) and two space dimensions (2D).

Elliptic partial differential equations

CERN Document Server

Han, Qing

2011-01-01

Elliptic Partial Differential Equations by Qing Han and FangHua Lin is one of the best textbooks I know. It is the perfect introduction to PDE. In 150 pages or so it covers an amazing amount of wonderful and extraordinary useful material. I have used it as a textbook at both graduate and undergraduate levels which is possible since it only requires very little background material yet it covers an enormous amount of material. In my opinion it is a must read for all interested in analysis and geometry, and for all of my own PhD students it is indeed just that. I cannot say enough good things abo

Handbook of differential equations stationary partial differential equations

CERN Document Server

Chipot, Michel

2006-01-01

This handbook is volume III in a series devoted to stationary partial differential quations. Similarly as volumes I and II, it is a collection of self contained state-of-the-art surveys written by well known experts in the field. The topics covered by this handbook include singular and higher order equations, problems near critically, problems with anisotropic nonlinearities, dam problem, T-convergence and Schauder-type estimates. These surveys will be useful for both beginners and experts and speed up the progress of corresponding (rapidly developing and fascinating) areas of mathematics. Ke

FORSIM, Solution of Ordinary or Partial Differential Equation with Initial Conditions

International Nuclear Information System (INIS)

Carver, M.B.

1985-01-01

1 - Description of problem or function: FORSIM is a FORTRAN oriented simulation program which automates the continuous transient solution of systems of ordinary and/or partial differential equations. The user writes his equations in a FORTRAN subroutine, following prescribed rules, and loads this routine along with the executive routines. The executive routines then read in initial data supplied by the user and proceed with the integration. 2 - Method of solution: Partial differential equations are converted to coupled ordinary differential equations by suitable discretization formulae. Integration is done by variable order, variable step-size error controlled algorithms. 3 - Restrictions on the complexity of the problem - Maximum of: 1000 ordinary differential equations

Differential geometry techniques for sets of nonlinear partial differential equations

Science.gov (United States)

Estabrook, Frank B.

1990-01-01

An attempt is made to show that the Cartan theory of partial differential equations can be a useful technique for applied mathematics. Techniques for finding consistent subfamilies of solutions that are generically rich and well-posed and for introducing potentials or other usefully consistent auxiliary fields are introduced. An extended sample calculation involving the Korteweg-de Vries equation is given.

Formulae and Bounds connected to Optimal Design and Homogenization of Partial Differential Operators and Integral Functionals

Energy Technology Data Exchange (ETDEWEB)

Lukkassen, D.

1996-12-31

When partial differential equations are set up to model physical processes in strongly heterogeneous materials, effective parameters for heat transfer, electric conductivity etc. are usually required. Averaging methods often lead to convergence problems and in homogenization theory one is therefore led to study how certain integral functionals behave asymptotically. This mathematical doctoral thesis discusses (1) means and bounds connected to homogenization of integral functionals, (2) reiterated homogenization of integral functionals, (3) bounds and homogenization of some particular partial differential operators, (4) applications and further results. 154 refs., 11 figs., 8 tabs.

A Numerical Method for Partial Differential Algebraic Equations Based on Differential Transform Method

Directory of Open Access Journals (Sweden)

Murat Osmanoglu

2013-01-01

Full Text Available We have considered linear partial differential algebraic equations (LPDAEs of the form , which has at least one singular matrix of . We have first introduced a uniform differential time index and a differential space index. The initial conditions and boundary conditions of the given system cannot be prescribed for all components of the solution vector here. To overcome this, we introduced these indexes. Furthermore, differential transform method has been given to solve LPDAEs. We have applied this method to a test problem, and numerical solution of the problem has been compared with analytical solution.

A model reduction approach to numerical inversion for a parabolic partial differential equation

International Nuclear Information System (INIS)

Borcea, Liliana; Druskin, Vladimir; Zaslavsky, Mikhail; Mamonov, Alexander V

2014-01-01

We propose a novel numerical inversion algorithm for the coefficients of parabolic partial differential equations, based on model reduction. The study is motivated by the application of controlled source electromagnetic exploration, where the unknown is the subsurface electrical resistivity and the data are time resolved surface measurements of the magnetic field. The algorithm presented in this paper considers inversion in one and two dimensions. The reduced model is obtained with rational interpolation in the frequency (Laplace) domain and a rational Krylov subspace projection method. It amounts to a nonlinear mapping from the function space of the unknown resistivity to the small dimensional space of the parameters of the reduced model. We use this mapping as a nonlinear preconditioner for the Gauss–Newton iterative solution of the inverse problem. The advantage of the inversion algorithm is twofold. First, the nonlinear preconditioner resolves most of the nonlinearity of the problem. Thus the iterations are less likely to get stuck in local minima and the convergence is fast. Second, the inversion is computationally efficient because it avoids repeated accurate simulations of the time-domain response. We study the stability of the inversion algorithm for various rational Krylov subspaces, and assess its performance with numerical experiments. (paper)

A model reduction approach to numerical inversion for a parabolic partial differential equation

Science.gov (United States)

Borcea, Liliana; Druskin, Vladimir; Mamonov, Alexander V.; Zaslavsky, Mikhail

2014-12-01

We propose a novel numerical inversion algorithm for the coefficients of parabolic partial differential equations, based on model reduction. The study is motivated by the application of controlled source electromagnetic exploration, where the unknown is the subsurface electrical resistivity and the data are time resolved surface measurements of the magnetic field. The algorithm presented in this paper considers inversion in one and two dimensions. The reduced model is obtained with rational interpolation in the frequency (Laplace) domain and a rational Krylov subspace projection method. It amounts to a nonlinear mapping from the function space of the unknown resistivity to the small dimensional space of the parameters of the reduced model. We use this mapping as a nonlinear preconditioner for the Gauss-Newton iterative solution of the inverse problem. The advantage of the inversion algorithm is twofold. First, the nonlinear preconditioner resolves most of the nonlinearity of the problem. Thus the iterations are less likely to get stuck in local minima and the convergence is fast. Second, the inversion is computationally efficient because it avoids repeated accurate simulations of the time-domain response. We study the stability of the inversion algorithm for various rational Krylov subspaces, and assess its performance with numerical experiments.

Partial differential equations an introduction

CERN Document Server

Colton, David

2004-01-01

Intended for a college senior or first-year graduate-level course in partial differential equations, this text offers students in mathematics, engineering, and the applied sciences a solid foundation for advanced studies in mathematics. Classical topics presented in a modern context include coverage of integral equations and basic scattering theory. This complete and accessible treatment includes a variety of examples of inverse problems arising from improperly posed applications. Exercises at the ends of chapters, many with answers, offer a clear progression in developing an understanding of

Partial differential equations in several complex variables

CERN Document Server

Chen, So-Chin

2001-01-01

This book is intended both as an introductory text and as a reference book for those interested in studying several complex variables in the context of partial differential equations. In the last few decades, significant progress has been made in the fields of Cauchy-Riemann and tangential Cauchy-Riemann operators. This book gives an up-to-date account of the theories for these equations and their applications. The background material in several complex variables is developed in the first three chapters, leading to the Levi problem. The next three chapters are devoted to the solvability and regularity of the Cauchy-Riemann equations using Hilbert space techniques. The authors provide a systematic study of the Cauchy-Riemann equations and the \\bar\\partial-Neumann problem, including L^2 existence theorems on pseudoconvex domains, \\frac 12-subelliptic estimates for the \\bar\\partial-Neumann problems on strongly pseudoconvex domains, global regularity of \\bar\\partial on more general pseudoconvex domains, boundary ...

Analytic continuation of solutions of some nonlinear convolution partial differential equations

Directory of Open Access Journals (Sweden)

Hidetoshi Tahara

2015-01-01

Full Text Available The paper considers a problem of analytic continuation of solutions of some nonlinear convolution partial differential equations which naturally appear in the summability theory of formal solutions of nonlinear partial differential equations. Under a suitable assumption it is proved that any local holomorphic solution has an analytic extension to a certain sector and its extension has exponential growth when the variable goes to infinity in the sector.

Convergence criteria for systems of nonlinear elliptic partial differential equations

International Nuclear Information System (INIS)

Sharma, R.K.

1986-01-01

This thesis deals with convergence criteria for a special system of nonlinear elliptic partial differential equations. A fixed-point algorithm is used, which iteratively solves one linearized elliptic partial differential equation at a time. Conditions are established that help foresee the convergence of the algorithm. Under reasonable hypotheses it is proved that the algorithm converges for such nonlinear elliptic systems. Extensive experimental results are reported and they show the algorithm converges in a wide variety of cases and the convergence is well correlated with the theoretical conditions introduced in this thesis

Parameter Estimation for Partial Differential Equations by Collage-Based Numerical Approximation

Directory of Open Access Journals (Sweden)

Xiaoyan Deng

2009-01-01

into a minimization problem of a function of several variables after the partial differential equation is approximated by a differential dynamical system. Then numerical schemes for solving this minimization problem are proposed, including grid approximation and ant colony optimization. The proposed schemes are applied to a parameter estimation problem for the Belousov-Zhabotinskii equation, and the results show that the proposed approximation method is efficient for both linear and nonlinear partial differential equations with respect to unknown parameters. At worst, the presented method provides an excellent starting point for traditional inversion methods that must first select a good starting point.

Algebraic dynamics solutions and algebraic dynamics algorithm for nonlinear partial differential evolution equations of dynamical systems

Institute of Scientific and Technical Information of China (English)

2008-01-01

Using functional derivative technique in quantum field theory, the algebraic dy-namics approach for solution of ordinary differential evolution equations was gen-eralized to treat partial differential evolution equations. The partial differential evo-lution equations were lifted to the corresponding functional partial differential equations in functional space by introducing the time translation operator. The functional partial differential evolution equations were solved by algebraic dynam-ics. The algebraic dynamics solutions are analytical in Taylor series in terms of both initial functions and time. Based on the exact analytical solutions, a new nu-merical algorithm—algebraic dynamics algorithm was proposed for partial differ-ential evolution equations. The difficulty of and the way out for the algorithm were discussed. The application of the approach to and computer numerical experi-ments on the nonlinear Burgers equation and meteorological advection equation indicate that the algebraic dynamics approach and algebraic dynamics algorithm are effective to the solution of nonlinear partial differential evolution equations both analytically and numerically.

A Priori Regularity of Parabolic Partial Differential Equations

KAUST Repository

Berkemeier, Francisco

2018-01-01

In this thesis, we consider parabolic partial differential equations such as the heat equation, the Fokker-Planck equation, and the porous media equation. Our aim is to develop methods that provide a priori estimates for solutions with singular

Auxiliary equation method for solving nonlinear partial differential equations

International Nuclear Information System (INIS)

Sirendaoreji,; Jiong, Sun

2003-01-01

By using the solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct several kinds of exact travelling wave solutions for some nonlinear partial differential equations. By this method some physically important nonlinear equations are investigated and new exact travelling wave solutions are explicitly obtained with the aid of symbolic computation

Image denoising using new pixon representation based on fuzzy filtering and partial differential equations

DEFF Research Database (Denmark)

Nadernejad, Ehsan; Nikpour, Mohsen

2012-01-01

In this paper, we have proposed two extensions to pixon-based image modeling. The first one is using bicubic interpolation instead of bilinear interpolation and the second one is using fuzzy filtering method, aiming to improve the quality of the pixonal image. Finally, partial differential...

Functional Determinants for Radially Separable Partial Differential Operators

Directory of Open Access Journals (Sweden)

G. V. Dunne

2007-01-01

Full Text Available Functional determinants of differential operators play a prominent role in many fields of theoretical and mathematical physics, ranging from condensed matter physics, to atomic, molecular and particle physics. They are, however, difficult to compute reliably in non-trivial cases. In one dimensional problems (i.e. functional determinants of ordinary differential operators, a classic result of Gel’fand and Yaglom greatly simplifies the computation of functional determinants. Here I report some recent progress in extending this approach to higher dimensions (i.e., functional determinants of partial differential operators, with applications in quantum field theory. 

Discrete variational derivative method a structure-preserving numerical method for partial differential equations

CERN Document Server

Furihata, Daisuke

2010-01-01

Nonlinear Partial Differential Equations (PDEs) have become increasingly important in the description of physical phenomena. Unlike Ordinary Differential Equations, PDEs can be used to effectively model multidimensional systems. The methods put forward in Discrete Variational Derivative Method concentrate on a new class of ""structure-preserving numerical equations"" which improves the qualitative behaviour of the PDE solutions and allows for stable computing. The authors have also taken care to present their methods in an accessible manner, which means that the book will be useful to engineer

Stability test for a parabolic partial differential equation

NARCIS (Netherlands)

Vajta, Miklos

2001-01-01

The paper describes a stability test applied to coupled parabolic partial differential equations. The PDE's describe the temperature distribution of composite structures with linear inner heat sources. The distributed transfer functions are developed based on the transmission matrix of each layer.

Soliton solution for nonlinear partial differential equations by cosine-function method

International Nuclear Information System (INIS)

Ali, A.H.A.; Soliman, A.A.; Raslan, K.R.

2007-01-01

In this Letter, we established a traveling wave solution by using Cosine-function algorithm for nonlinear partial differential equations. The method is used to obtain the exact solutions for five different types of nonlinear partial differential equations such as, general equal width wave equation (GEWE), general regularized long wave equation (GRLW), general Korteweg-de Vries equation (GKdV), general improved Korteweg-de Vries equation (GIKdV), and Coupled equal width wave equations (CEWE), which are the important soliton equations

Representations of Lie algebras and partial differential equations

CERN Document Server

Xu, Xiaoping

2017-01-01

This book provides explicit representations of finite-dimensional simple Lie algebras, related partial differential equations, linear orthogonal algebraic codes, combinatorics and algebraic varieties, summarizing the author’s works and his joint works with his former students.  Further, it presents various oscillator generalizations of the classical representation theorem on harmonic polynomials, and highlights new functors from the representation category of a simple Lie algebra to that of another simple Lie algebra. Partial differential equations play a key role in solving certain representation problems. The weight matrices of the minimal and adjoint representations over the simple Lie algebras of types E and F are proved to generate ternary orthogonal linear codes with large minimal distances. New multi-variable hypergeometric functions related to the root systems of simple Lie algebras are introduced in connection with quantum many-body systems in one dimension. In addition, the book identifies certai...

Statistical mechanics of normal grain growth in one dimension: A partial integro-differential equation model

International Nuclear Information System (INIS)

Ng, Felix S.L.

2016-01-01

We develop a statistical-mechanical model of one-dimensional normal grain growth that does not require any drift-velocity parameterization for grain size, such as used in the continuity equation of traditional mean-field theories. The model tracks the population by considering grain sizes in neighbour pairs; the probability of a pair having neighbours of certain sizes is determined by the size-frequency distribution of all pairs. Accordingly, the evolution obeys a partial integro-differential equation (PIDE) over ‘grain size versus neighbour grain size’ space, so that the grain-size distribution is a projection of the PIDE's solution. This model, which is applicable before as well as after statistically self-similar grain growth has been reached, shows that the traditional continuity equation is invalid outside this state. During statistically self-similar growth, the PIDE correctly predicts the coarsening rate, invariant grain-size distribution and spatial grain size correlations observed in direct simulations. The PIDE is then reducible to the standard continuity equation, and we derive an explicit expression for the drift velocity. It should be possible to formulate similar parameterization-free models of normal grain growth in two and three dimensions.

Non-linear partial differential equations an algebraic view of generalized solutions

CERN Document Server

Rosinger, Elemer E

1990-01-01

A massive transition of interest from solving linear partial differential equations to solving nonlinear ones has taken place during the last two or three decades. The availability of better computers has often made numerical experimentations progress faster than the theoretical understanding of nonlinear partial differential equations. The three most important nonlinear phenomena observed so far both experimentally and numerically, and studied theoretically in connection with such equations have been the solitons, shock waves and turbulence or chaotical processes. In many ways, these phenomen

Optimal Control Problems for Partial Differential Equations on Reticulated Domains

CERN Document Server

Kogut, Peter I

2011-01-01

In the development of optimal control, the complexity of the systems to which it is applied has increased significantly, becoming an issue in scientific computing. In order to carry out model-reduction on these systems, the authors of this work have developed a method based on asymptotic analysis. Moving from abstract explanations to examples and applications with a focus on structural network problems, they aim at combining techniques of homogenization and approximation. Optimal Control Problems for Partial Differential Equations on Reticulated Domains is an excellent reference tool for gradu

Using Partial Differential Equations for Pricing of Goods and Services

Directory of Open Access Journals (Sweden)

Traykov Metodi

2016-06-01

Full Text Available This article is based on the methodology of comparative analysis, using an innovative approach for pricing of various goods and services. Benchmarking is the continuous search to find and adapt better pricing methods that leading to increased profits. We will consider the numerical solution of partial differential equations, based on Black-Scholes model for pricing of goods and services within European option. Also, we will present formulation and numerical behavior of explicit and implicit methods that can be use in pricing for company assets within European option.

Higher order multi-term time-fractional partial differential equations involving Caputo-Fabrizio derivative

OpenAIRE

Erkinjon Karimov; Sardor Pirnafasov

2017-01-01

In this work we discuss higher order multi-term partial differential equation (PDE) with the Caputo-Fabrizio fractional derivative in time. Using method of separation of variables, we reduce fractional order partial differential equation to the integer order. We represent explicit solution of formulated problem in particular case by Fourier series.

Exact solutions of some nonlinear partial differential equations using ...

Indian Academy of Sciences (India)

Nonlinear partial differential equations (NPDEs) are encountered in various ... such as physics, mechanics, chemistry, biology, mathematics and engineering. ... In §3, this method is applied to the generalized forms of Klein–Gordon equation,.

Advanced Topics in Computational Partial Differential Equations: Numerical Methods and Diffpack Programming

International Nuclear Information System (INIS)

Katsaounis, T D

2005-01-01

The scope of this book is to present well known simple and advanced numerical methods for solving partial differential equations (PDEs) and how to implement these methods using the programming environment of the software package Diffpack. A basic background in PDEs and numerical methods is required by the potential reader. Further, a basic knowledge of the finite element method and its implementation in one and two space dimensions is required. The authors claim that no prior knowledge of the package Diffpack is required, which is true, but the reader should be at least familiar with an object oriented programming language like C++ in order to better comprehend the programming environment of Diffpack. Certainly, a prior knowledge or usage of Diffpack would be a great advantage to the reader. The book consists of 15 chapters, each one written by one or more authors. Each chapter is basically divided into two parts: the first part is about mathematical models described by PDEs and numerical methods to solve these models and the second part describes how to implement the numerical methods using the programming environment of Diffpack. Each chapter closes with a list of references on its subject. The first nine chapters cover well known numerical methods for solving the basic types of PDEs. Further, programming techniques on the serial as well as on the parallel implementation of numerical methods are also included in these chapters. The last five chapters are dedicated to applications, modelled by PDEs, in a variety of fields. In summary, the book focuses on the computational and implementational issues involved in solving partial differential equations. The potential reader should have a basic knowledge of PDEs and the finite difference and finite element methods. The examples presented are solved within the programming framework of Diffpack and the reader should have prior experience with the particular software in order to take full advantage of the book. Overall

Higher order multi-term time-fractional partial differential equations involving Caputo-Fabrizio derivative

Directory of Open Access Journals (Sweden)

Erkinjon Karimov

2017-10-01

Full Text Available In this work we discuss higher order multi-term partial differential equation (PDE with the Caputo-Fabrizio fractional derivative in time. Using method of separation of variables, we reduce fractional order partial differential equation to the integer order. We represent explicit solution of formulated problem in particular case by Fourier series.

Unified algorithm for partial differential equations and examples of numerical computation

International Nuclear Information System (INIS)

Watanabe, Tsuguhiro

1999-01-01

A new unified algorithm is proposed to solve partial differential equations which describe nonlinear boundary value problems, eigenvalue problems and time developing boundary value problems. The algorithm is composed of implicit difference scheme and multiple shooting scheme and is named as HIDM (Higher order Implicit Difference Method). A new prototype computer programs for 2-dimensional partial differential equations is constructed and tested successfully to several problems. Extension of the computer programs to 3 or more higher order dimension problems will be easy due to the direct product type difference scheme. (author)

Construction and accuracy of partial differential equation approximations to the chemical master equation.

Science.gov (United States)

Grima, Ramon

2011-11-01

The mesoscopic description of chemical kinetics, the chemical master equation, can be exactly solved in only a few simple cases. The analytical intractability stems from the discrete character of the equation, and hence considerable effort has been invested in the development of Fokker-Planck equations, second-order partial differential equation approximations to the master equation. We here consider two different types of higher-order partial differential approximations, one derived from the system-size expansion and the other from the Kramers-Moyal expansion, and derive the accuracy of their predictions for chemical reactive networks composed of arbitrary numbers of unimolecular and bimolecular reactions. In particular, we show that the partial differential equation approximation of order Q from the Kramers-Moyal expansion leads to estimates of the mean number of molecules accurate to order Ω(-(2Q-3)/2), of the variance of the fluctuations in the number of molecules accurate to order Ω(-(2Q-5)/2), and of skewness accurate to order Ω(-(Q-2)). We also show that for large Q, the accuracy in the estimates can be matched only by a partial differential equation approximation from the system-size expansion of approximate order 2Q. Hence, we conclude that partial differential approximations based on the Kramers-Moyal expansion generally lead to considerably more accurate estimates in the mean, variance, and skewness than approximations of the same order derived from the system-size expansion.

Tracer kinetics: Modelling by partial differential equations of inhomogeneous compartments with age-dependent elimination rates. Pt. 1

International Nuclear Information System (INIS)

Winkler, E.

1991-01-01

Mathematical models in tracer kinetics are usually based on ordinary differential equations which correspond to a system of kinetically homogeneous compartments (standard compartments). A generalization is possible by the admission of inhomogeneities in the behaviour of the elements belonging to a compartment. The important special case of the age-dependence of elimination rates is treated in its deterministic version. It leads to partial different equations (i.e., systems with distributed coefficients) with the 'age' or the 'residence time' of an element of the compartment as a variable additional to 'time'. The basic equations for one generalized compartment and for systems of such compartments are given together with their general solutions. (orig.) [de

Partial differential equations mathematical techniques for engineers

CERN Document Server

Epstein, Marcelo

2017-01-01

This monograph presents a graduate-level treatment of partial differential equations (PDEs) for engineers. The book begins with a review of the geometrical interpretation of systems of ODEs, the appearance of PDEs in engineering is motivated by the general form of balance laws in continuum physics. Four chapters are devoted to a detailed treatment of the single first-order PDE, including shock waves and genuinely non-linear models, with applications to traffic design and gas dynamics. The rest of the book deals with second-order equations. In the treatment of hyperbolic equations, geometric arguments are used whenever possible and the analogy with discrete vibrating systems is emphasized. The diffusion and potential equations afford the opportunity of dealing with questions of uniqueness and continuous dependence on the data, the Fourier integral, generalized functions (distributions), Duhamel's principle, Green's functions and Dirichlet and Neumann problems. The target audience primarily comprises graduate s...

Book review: Partial Differential Equations and Fluid Mechanics

NARCIS (Netherlands)

Muntean, A.

2011-01-01

The baak is the result of the workshop Partial Differential Equations and Fluid Dynamics that look place at the Mathematics Institute of the University of Warwick. May 21st - 23rd, 2007. It contains ten review and research papers which provide an accessible summary of a wide range of active research

Introduction to partial differential equations and Hilbert space methods

CERN Document Server

Gustafson, Karl E

1997-01-01

Easy-to-use text examines principal method of solving partial differential equations, 1st-order systems, computation methods, and much more. Over 600 exercises, with answers for many. Ideal for a 1-semester or full-year course.

Calculation of similarity solutions of partial differential equations

International Nuclear Information System (INIS)

Dresner, L.

1980-08-01

When a partial differential equation in two independent variables is invariant to a group G of stretching transformations, it has similarity solutions that can be found by solving an ordinary differential equation. Under broad conditions, this ordinary differential equation is also invariant to another stretching group G', related to G. The invariance of the ordinary differential equation to G' can be used to simplify its solution, particularly if it is of second order. Then a method of Lie's can be used to reduce it to a first-order equation, the study of which is greatly facilitated by analysis of its direction field. The method developed here is applied to three examples: Blasius's equation for boundary layer flow over a flat plate and two nonlinear diffusion equations, cc/sub t/ = c/sub zz/ and c/sub t/ = (cc/sub z/)/sub z/

Optimal Control Method of Parabolic Partial Differential Equations and Its Application to Heat Transfer Model in Continuous Cast Secondary Cooling Zone

Directory of Open Access Journals (Sweden)

Yuan Wang

2015-01-01

Full Text Available Our work is devoted to a class of optimal control problems of parabolic partial differential equations. Because of the partial differential equations constraints, it is rather difficult to solve the optimization problem. The gradient of the cost function can be found by the adjoint problem approach. Based on the adjoint problem approach, the gradient of cost function is proved to be Lipschitz continuous. An improved conjugate method is applied to solve this optimization problem and this algorithm is proved to be convergent. This method is applied to set-point values in continuous cast secondary cooling zone. Based on the real data in a plant, the simulation experiments show that the method can ensure the steel billet quality. From these experiment results, it is concluded that the improved conjugate gradient algorithm is convergent and the method is effective in optimal control problem of partial differential equations.

Numerical models for differential problems

CERN Document Server

Quarteroni, Alfio

2017-01-01

In this text, we introduce the basic concepts for the numerical modelling of partial differential equations. We consider the classical elliptic, parabolic and hyperbolic linear equations, but also the diffusion, transport, and Navier-Stokes equations, as well as equations representing conservation laws, saddle-point problems and optimal control problems. Furthermore, we provide numerous physical examples which underline such equations. We then analyze numerical solution methods based on finite elements, finite differences, finite volumes, spectral methods and domain decomposition methods, and reduced basis methods. In particular, we discuss the algorithmic and computer implementation aspects and provide a number of easy-to-use programs. The text does not require any previous advanced mathematical knowledge of partial differential equations: the absolutely essential concepts are reported in a preliminary chapter. It is therefore suitable for students of bachelor and master courses in scientific disciplines, an...

Mobile point sensors and actuators in the controllability theory of partial differential equations

CERN Document Server

Khapalov, Alexander Y

2017-01-01

This book presents a concise study of controllability theory of partial differential equations when they are equipped with actuators and/or sensors that are finite dimensional at every moment of time. Based on the author’s extensive research in the area of controllability theory, this monograph specifically focuses on the issues of controllability, observability, and stabilizability for parabolic and hyperbolic partial differential equations. The topics in this book also cover related applied questions such as the problem of localization of unknown pollution sources based on information obtained from point sensors that arise in environmental monitoring. Researchers and graduate students interested in controllability theory of partial differential equations and its applications will find this book to be an invaluable resource to their studies.

Numerical Methods for Partial Differential Equations

CERN Document Server

Guo, Ben-yu

1987-01-01

These Proceedings of the first Chinese Conference on Numerical Methods for Partial Differential Equations covers topics such as difference methods, finite element methods, spectral methods, splitting methods, parallel algorithm etc., their theoretical foundation and applications to engineering. Numerical methods both for boundary value problems of elliptic equations and for initial-boundary value problems of evolution equations, such as hyperbolic systems and parabolic equations, are involved. The 16 papers of this volume present recent or new unpublished results and provide a good overview of current research being done in this field in China.

Function spaces and partial differential equations 2 volume set

CERN Document Server

Taheri, Ali

2015-01-01

This is a book written primarily for graduate students and early researchers in the fields of Analysis and Partial Differential Equations (PDEs). Coverage of the material is essentially self-contained, extensive and novel with great attention to details and rigour.

The convergence of the order sequence and the solution function sequence on fractional partial differential equation

Science.gov (United States)

Rusyaman, E.; Parmikanti, K.; Chaerani, D.; Asefan; Irianingsih, I.

2018-03-01

One of the application of fractional ordinary differential equation is related to the viscoelasticity, i.e., a correlation between the viscosity of fluids and the elasticity of solids. If the solution function develops into function with two or more variables, then its differential equation must be changed into fractional partial differential equation. As the preliminary study for two variables viscoelasticity problem, this paper discusses about convergence analysis of function sequence which is the solution of the homogenous fractional partial differential equation. The method used to solve the problem is Homotopy Analysis Method. The results show that if given two real number sequences (αn) and (βn) which converge to α and β respectively, then the solution function sequences of fractional partial differential equation with order (αn, βn) will also converge to the solution function of fractional partial differential equation with order (α, β).

Stochastic partial differential equations a modeling, white noise functional approach

CERN Document Server

Holden, Helge; Ubøe, Jan; Zhang, Tusheng

1996-01-01

This book is based on research that, to a large extent, started around 1990, when a research project on fluid flow in stochastic reservoirs was initiated by a group including some of us with the support of VISTA, a research coopera­ tion between the Norwegian Academy of Science and Letters and Den norske stats oljeselskap A.S. (Statoil). The purpose of the project was to use stochastic partial differential equations (SPDEs) to describe the flow of fluid in a medium where some of the parameters, e.g., the permeability, were stochastic or "noisy". We soon realized that the theory of SPDEs at the time was insufficient to handle such equations. Therefore it became our aim to develop a new mathematically rigorous theory that satisfied the following conditions. 1) The theory should be physically meaningful and realistic, and the corre­ sponding solutions should make sense physically and should be useful in applications. 2) The theory should be general enough to handle many of the interesting SPDEs that occur in r...

Partial differential equations methods, applications and theories

CERN Document Server

Hattori, Harumi

2013-01-01

This volume is an introductory level textbook for partial differential equations (PDE's) and suitable for a one-semester undergraduate level or two-semester graduate level course in PDE's or applied mathematics. Chapters One to Five are organized according to the equations and the basic PDE's are introduced in an easy to understand manner. They include the first-order equations and the three fundamental second-order equations, i.e. the heat, wave and Laplace equations. Through these equations we learn the types of problems, how we pose the problems, and the methods of solutions such as the separation of variables and the method of characteristics. The modeling aspects are explained as well. The methods introduced in earlier chapters are developed further in Chapters Six to Twelve. They include the Fourier series, the Fourier and the Laplace transforms, and the Green's functions. The equations in higher dimensions are also discussed in detail. This volume is application-oriented and rich in examples. Going thr...

Advances in nonlinear partial differential equations and stochastics

CERN Document Server

Kawashima, S

1998-01-01

In the past two decades, there has been great progress in the theory of nonlinear partial differential equations. This book describes the progress, focusing on interesting topics in gas dynamics, fluid dynamics, elastodynamics etc. It contains ten articles, each of which discusses a very recent result obtained by the author. Some of these articles review related results.

Reciprocal links among differential parenting, perceived partiality, and self-worth: a three-wave longitudinal study.

Science.gov (United States)

Shebloski, Barbara; Conger, Katherine J; Widaman, Keith F

2005-12-01

This study examined reciprocal links between parental differential treatment, siblings' perception of partiality, and self-worth with 3 waves of data from 384 adolescent sibling dyads. Results suggest that birth-order status was significantly associated with self-worth and perception of maternal and paternal differential treatment. There was a consistent across-time effect of self-worth on perception of parental partiality for later born siblings, but not earlier born siblings, and a consistent effect of differential treatment on perception of partiality for earlier born but not later born siblings. The results contribute new insight into the associations between perception of differential parenting and adolescents' adjustment and the role of birth order. Copyright 2006 APA, all rights reserved).

Energy preserving integration of bi-Hamiltonian partial differential equations

NARCIS (Netherlands)

Karasozen, B.; Simsek, G.

2013-01-01

The energy preserving average vector field (AVF) integrator is applied to evolutionary partial differential equations (PDEs) in bi-Hamiltonian form with nonconstant Poisson structures. Numerical results for the Korteweg de Vries (KdV) equation and for the Ito type coupled KdV equation confirm the

Reconsidering harmonic and anharmonic coherent states: Partial differential equations approach

Energy Technology Data Exchange (ETDEWEB)

Toutounji, Mohamad, E-mail: Mtoutounji@uaeu.ac.ae

2015-02-15

This article presents a new approach to dealing with time dependent quantities such as autocorrelation function of harmonic and anharmonic systems using coherent states and partial differential equations. The approach that is normally used to evaluate dynamical quantities involves formidable operator algebra. That operator algebra becomes insurmountable when employing Morse oscillator coherent states. This problem becomes even more complicated in case of Morse oscillator as it tends to exhibit divergent dynamics. This approach employs linear partial differential equations, some of which may be solved exactly and analytically, thereby avoiding the cumbersome noncommutative algebra required to manipulate coherent states of Morse oscillator. Additionally, the arising integrals while using the herein presented method feature stability and high numerical efficiency. The correctness, applicability, and utility of the above approach are tested by reproducing the partition and optical autocorrelation function of the harmonic oscillator. A closed-form expression for the equilibrium canonical partition function of the Morse oscillator is derived using its coherent states and partial differential equations. Also, a nonequilibrium autocorrelation function expression for weak electron–phonon coupling in condensed systems is derived for displaced Morse oscillator in electronic state. Finally, the utility of the method is demonstrated through further simplifying the Morse oscillator partition function or autocorrelation function expressions reported by other researchers in unevaluated form of second-order derivative exponential. Comparison with exact dynamics shows identical results.

Partially nested designs in psychotherapy trials: A review of modeling developments.

Science.gov (United States)

Sterba, Sonya K

2017-07-01

Individually-randomized psychotherapy trials are often partially nested. For instance, individuals assigned to a treatment arm may be clustered into therapy groups for purposes of treatment administration, whereas individuals assigned to a wait-list control are unclustered. The past several years have seen rapid expansion and investigation of methods for analyzing partially nested data. Yet partial nesting often remains ignored in psychotherapy trials. This review integrates and disseminates developments in the analysis of partially nested data that are particularly relevant for psychotherapy researchers. First, we differentiate among alternative partially nested designs. Then, we present adaptations of multilevel model specifications that accommodate each design. Next, we address how moderation by treatment as well as mediation of the treatment effect can be investigated in partially nested designs. Model fitting results, annotated software syntax, and illustrative data sets are provided and key methodological issues are discussed. We emphasize that cluster-level variability in the treatment arm need not be considered a nuisance; it can be modeled to yield insights about the treatment process.

A deterministic partial differential equation model for dose calculation in electron radiotherapy.

Science.gov (United States)

Duclous, R; Dubroca, B; Frank, M

2010-07-07

High-energy ionizing radiation is a prominent modality for the treatment of many cancers. The approaches to electron dose calculation can be categorized into semi-empirical models (e.g. Fermi-Eyges, convolution-superposition) and probabilistic methods (e.g.Monte Carlo). A third approach to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. We derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on Berthon et al (2007 J. Sci. Comput. 31 347-89) that exactly preserves the key properties of the analytical solution on the discrete level. We discuss several test cases taken from the medical physics literature. A test case with an academic Henyey-Greenstein scattering kernel is considered. We compare our model to a benchmark discrete ordinate solution. A simplified model of electron interactions with tissue is employed to compute the dose of an electron beam in a water phantom, and a case of irradiation of the vertebral column. Here our model is compared to the PENELOPE Monte Carlo code. In the academic example, the fluences computed with the new model and a benchmark result differ by less than 1%. The depths at half maximum differ by less than 0.6%. In the two comparisons with Monte Carlo, our model gives qualitatively reasonable dose distributions. Due to the crude interaction model, these so far do not have the accuracy needed in clinical practice. However, the new model has a computational cost that is less than one-tenth of the cost of a Monte Carlo simulation. In addition, simulations can be set up in a similar way as a Monte Carlo simulation. If more detailed effects such as coupled electron-photon transport, bremsstrahlung

A deterministic partial differential equation model for dose calculation in electron radiotherapy

Science.gov (United States)

Duclous, R.; Dubroca, B.; Frank, M.

2010-07-01

High-energy ionizing radiation is a prominent modality for the treatment of many cancers. The approaches to electron dose calculation can be categorized into semi-empirical models (e.g. Fermi-Eyges, convolution-superposition) and probabilistic methods (e.g. Monte Carlo). A third approach to dose calculation has only recently attracted attention in the medical physics community. This approach is based on the deterministic kinetic equations of radiative transfer. We derive a macroscopic partial differential equation model for electron transport in tissue. This model involves an angular closure in the phase space. It is exact for the free streaming and the isotropic regime. We solve it numerically by a newly developed HLLC scheme based on Berthon et al (2007 J. Sci. Comput. 31 347-89) that exactly preserves the key properties of the analytical solution on the discrete level. We discuss several test cases taken from the medical physics literature. A test case with an academic Henyey-Greenstein scattering kernel is considered. We compare our model to a benchmark discrete ordinate solution. A simplified model of electron interactions with tissue is employed to compute the dose of an electron beam in a water phantom, and a case of irradiation of the vertebral column. Here our model is compared to the PENELOPE Monte Carlo code. In the academic example, the fluences computed with the new model and a benchmark result differ by less than 1%. The depths at half maximum differ by less than 0.6%. In the two comparisons with Monte Carlo, our model gives qualitatively reasonable dose distributions. Due to the crude interaction model, these so far do not have the accuracy needed in clinical practice. However, the new model has a computational cost that is less than one-tenth of the cost of a Monte Carlo simulation. In addition, simulations can be set up in a similar way as a Monte Carlo simulation. If more detailed effects such as coupled electron-photon transport, bremsstrahlung

Partial differential equations and their applications

International Nuclear Information System (INIS)

Gauthier-Villars

1998-01-01

This book is dedicated to the French mathematician J.L.Lions. It represents a compilation of articles from about 80 authors. The topics treated are diverse but the more or less commune matter is the study of the characteristics of some partial differential equations. Stability, optimal approximation, numerical resolution, particular applications are among the subjects reviewed. An article deals with the MHD stability of fusion plasmas in tokamaks, another presents the scientific and technical challenges of nuclear energy in France. The latter that contains no equations can be considered as an enjoyable break in a sea of about 40 mathematical articles. (A.C.)

Exact solutions of some nonlinear partial differential equations using ...

Indian Academy of Sciences (India)

The functional variable method is a powerful solution method for obtaining exact solutions of some nonlinear partial differential equations. In this paper, the functional variable method is used to establish exact solutions of the generalized forms of Klein–Gordon equation, the (2 + 1)-dimensional Camassa–Holm ...

Workload Characterization of CFD Applications Using Partial Differential Equation Solvers

Science.gov (United States)

Waheed, Abdul; Yan, Jerry; Saini, Subhash (Technical Monitor)

1998-01-01

Workload characterization is used for modeling and evaluating of computing systems at different levels of detail. We present workload characterization for a class of Computational Fluid Dynamics (CFD) applications that solve Partial Differential Equations (PDEs). This workload characterization focuses on three high performance computing platforms: SGI Origin2000, EBM SP-2, a cluster of Intel Pentium Pro bases PCs. We execute extensive measurement-based experiments on these platforms to gather statistics of system resource usage, which results in workload characterization. Our workload characterization approach yields a coarse-grain resource utilization behavior that is being applied for performance modeling and evaluation of distributed high performance metacomputing systems. In addition, this study enhances our understanding of interactions between PDE solver workloads and high performance computing platforms and is useful for tuning these applications.

Implementation of geomechanical models for engineered clay barriers in multi-physic partial differential equation solvers

International Nuclear Information System (INIS)

Navarro, V.; Alonso, J.; Asensio, L.; Yustres, A.; Pintado, X.

2012-01-01

Document available in extended abstract form only. The use of numerical methods, especially the Finite Element Method (FEM), for solving boundary problems in Unsaturated Soil Mechanics has experienced significant progress. Several codes, both built mainly for research purposes and commercial software, are now available. In the last years, Multi-physic Partial Differentiation Equation Solvers (MPDES) have turned out to be an interesting proposal. In this family of solvers, the user defines the governing equations and the behaviour models, generally using a computer algebra environment. The code automatically assembles and solves the equation systems, saving the user having to redefine the structures of memory storage or to implement solver algorithms. The user can focus on the definition of the physics of the problem, while it is possible to couple virtually any physical or chemical process that can be described by a PDE. This can be done, for instance, in COMSOL Multiphysics (CM). Nonetheless, the versatility of CM is compromised by the impossibility to implement models with variables defined by implicit functions. Elasto-plastic models involve an implicit coupling among stress increments, plastic strains and plastic variables increments. For this reason, they cannot be implemented in CM in a straightforward way. This means a very relevant limitation for the use of this tool in the analysis of geomechanical boundary value problems. In this work, a strategy to overcome this problem using the multi-physics concept is presented. A mixed method is proposed, considering the constitutive stresses, the pre-consolidation pressure and the plastic variables as main unknowns of the model. Mixed methods usually present stability problems. However, the algorithmics present in CM include several numerical strategies to minimise this kind of problems. Besides, CM is based on the application of the FEM with Lagrange multipliers, an approach that significantly contributes stability

A practical course in differential equations and mathematical modeling

CERN Document Server

Ibragimov , Nail H

2009-01-01

A Practical Course in Differential Equations and Mathematical Modelling is a unique blend of the traditional methods of ordinary and partial differential equations with Lie group analysis enriched by the author's own theoretical developments. The book which aims to present new mathematical curricula based on symmetry and invariance principles is tailored to develop analytic skills and working knowledge in both classical and Lie's methods for solving linear and nonlinear equations. This approach helps to make courses in differential equations, mathematical modelling, distributions and fundame

Partial differential equations for self-organization in cellular and developmental biology

International Nuclear Information System (INIS)

Baker, R E; Gaffney, E A; Maini, P K

2008-01-01

Understanding the mechanisms governing and regulating the emergence of structure and heterogeneity within cellular systems, such as the developing embryo, represents a multiscale challenge typifying current integrative biology research, namely, explaining the macroscale behaviour of a system from microscale dynamics. This review will focus upon modelling how cell-based dynamics orchestrate the emergence of higher level structure. After surveying representative biological examples and the models used to describe them, we will assess how developments at the scale of molecular biology have impacted on current theoretical frameworks, and the new modelling opportunities that are emerging as a result. We shall restrict our survey of mathematical approaches to partial differential equations and the tools required for their analysis. We will discuss the gap between the modelling abstraction and biological reality, the challenges this presents and highlight some open problems in the field. (invited article)

Some overdetermined systems of complex partial differential equations

International Nuclear Information System (INIS)

Le Hung Son.

1990-01-01

In this paper we extend some properties of analytic functions on several complex variables to solutions of overdetermined systems of complex partial differential equations. It is proved that many global properties of analytic functions are true for solutions of the Vekua system in special cases. The relation between analytic functions and solutions of quasi-linear systems is discussed in the paper. (author). 8 refs

A concise course on stochastic partial differential equations

CERN Document Server

Prévôt, Claudia

2007-01-01

These lectures concentrate on (nonlinear) stochastic partial differential equations (SPDE) of evolutionary type. All kinds of dynamics with stochastic influence in nature or man-made complex systems can be modelled by such equations. To keep the technicalities minimal we confine ourselves to the case where the noise term is given by a stochastic integral w.r.t. a cylindrical Wiener process.But all results can be easily generalized to SPDE with more general noises such as, for instance, stochastic integral w.r.t. a continuous local martingale. There are basically three approaches to analyze SPDE: the "martingale measure approach", the "mild solution approach" and the "variational approach". The purpose of these notes is to give a concise and as self-contained as possible an introduction to the "variational approach". A large part of necessary background material, such as definitions and results from the theory of Hilbert spaces, are included in appendices.

Partially Observed Mixtures of IRT Models: An Extension of the Generalized Partial-Credit Model

Science.gov (United States)

Von Davier, Matthias; Yamamoto, Kentaro

2004-01-01

The generalized partial-credit model (GPCM) is used frequently in educational testing and in large-scale assessments for analyzing polytomous data. Special cases of the generalized partial-credit model are the partial-credit model--or Rasch model for ordinal data--and the two parameter logistic (2PL) model. This article extends the GPCM to the…

Multigrid methods for partial differential equations - a short introduction

International Nuclear Information System (INIS)

Linden, J.; Stueben, K.

1993-01-01

These notes summarize the multigrid methods and emphasis is laid on the algorithmic concepts of multigrid for solving linear and non-linear partial differential equations. In this paper there is brief description of the basic structure of multigrid methods. Detailed introduction is also contained with applications to VLSI process simulation. (A.B.)

Workshop on Recent Trends in Complex Methods for Partial Differential Equations

CERN Document Server

Celebi, A; Tutschke, Wolfgang

1999-01-01

This volume is a collection of manscripts mainly originating from talks and lectures given at the Workshop on Recent Trends in Complex Methods for Par­ tial Differential Equations held from July 6 to 10, 1998 at the Middle East Technical University in Ankara, Turkey, sponsored by The Scientific and Tech­ nical Research Council of Turkey and the Middle East Technical University. This workshop is a continuation oftwo workshops from 1988 and 1993 at the In­ ternational Centre for Theoretical Physics in Trieste, Italy entitled Functional analytic Methods in Complex Analysis and Applications to Partial Differential Equations. Since classical complex analysis of one and several variables has a long tra­ dition it is of high level. But most of its basic problems are solved nowadays so that within the last few decades it has lost more and more attention. The area of complex and functional analytic methods in partial differential equations, however, is still a growing and flourishing field, in particular as these ...

A partial differential equation for pseudocontact shift.

Science.gov (United States)

Charnock, G T P; Kuprov, Ilya

2014-10-07

It is demonstrated that pseudocontact shift (PCS), viewed as a scalar or a tensor field in three dimensions, obeys an elliptic partial differential equation with a source term that depends on the Hessian of the unpaired electron probability density. The equation enables straightforward PCS prediction and analysis in systems with delocalized unpaired electrons, particularly for the nuclei located in their immediate vicinity. It is also shown that the probability density of the unpaired electron may be extracted, using a regularization procedure, from PCS data.

New Solutions of Three Nonlinear Space- and Time-Fractional Partial Differential Equations in Mathematical Physics

International Nuclear Information System (INIS)

Yao Ruo-Xia; Wang Wei; Chen Ting-Hua

2014-01-01

Motivated by the widely used ansätz method and starting from the modified Riemann—Liouville derivative together with a fractional complex transformation that can be utilized to transform nonlinear fractional partial differential equations to nonlinear ordinary differential equations, new types of exact traveling wave solutions to three important nonlinear space- and time-fractional partial differential equations are obtained simultaneously in terms of solutions of a Riccati equation. The results are new and first reported in this paper. (general)

Oscillation of certain higher-order neutral partial functional differential equations.

Science.gov (United States)

Li, Wei Nian; Sheng, Weihong

2016-01-01

In this paper, we study the oscillation of certain higher-order neutral partial functional differential equations with the Robin boundary conditions. Some oscillation criteria are established. Two examples are given to illustrate the main results in the end of this paper.

Function spaces and partial differential equations volume 2 : contemporary analysis

CERN Document Server

Taheri, Ali

2015-01-01

This is a book written primarily for graduate students and early researchers in the fields of Analysis and Partial Differential Equations (PDEs). Coverage of the material is essentially self-contained, extensive and novel with great attention to details and rigour.

Partial Differential Equations in General Relativity

International Nuclear Information System (INIS)

Choquet-Bruhat, Yvonne

2008-01-01

General relativity is a physical theory basic in the modeling of the universe at the large and small scales. Its mathematical formulation, the Einstein partial differential equations, are geometrically simple, but intricate for the analyst, involving both hyperbolic and elliptic PDE, with local and global problems. Many problems remain open though remarkable progress has been made recently towards their solutions. Alan Rendall's book states, in a down-to-earth form, fundamental results used to solve different types of equations. In each case he gives applications to special models as well as to general properties of Einsteinian spacetimes. A chapter on ODE contains, in particular, a detailed discussion of Bianchi spacetimes. A chapter entitled 'Elliptic systems' treats the Einstein constraints. A chapter entitled 'Hyperbolic systems' is followed by a chapter on the Cauchy problem and a chapter 'Global results' which contains recently proved theorems. A chapter is dedicated to the Einstein-Vlasov system, of which the author is a specialist. On the whole, the book surveys, in a concise though precise way, many essential results of recent interest in mathematical general relativity, and it is very clearly written. Each chapter is followed by an up to date bibliography. In conclusion, this book will be a valuable asset to relativists who wish to learn clearly-stated mathematical results and to mathematicians who want to penetrate into the subtleties of general relativity, as a mathematical and physical theory. (book review)

Partial Differential Equations in General Relativity

Energy Technology Data Exchange (ETDEWEB)

Choquet-Bruhat, Yvonne

2008-09-07

General relativity is a physical theory basic in the modeling of the universe at the large and small scales. Its mathematical formulation, the Einstein partial differential equations, are geometrically simple, but intricate for the analyst, involving both hyperbolic and elliptic PDE, with local and global problems. Many problems remain open though remarkable progress has been made recently towards their solutions. Alan Rendall's book states, in a down-to-earth form, fundamental results used to solve different types of equations. In each case he gives applications to special models as well as to general properties of Einsteinian spacetimes. A chapter on ODE contains, in particular, a detailed discussion of Bianchi spacetimes. A chapter entitled 'Elliptic systems' treats the Einstein constraints. A chapter entitled 'Hyperbolic systems' is followed by a chapter on the Cauchy problem and a chapter 'Global results' which contains recently proved theorems. A chapter is dedicated to the Einstein-Vlasov system, of which the author is a specialist. On the whole, the book surveys, in a concise though precise way, many essential results of recent interest in mathematical general relativity, and it is very clearly written. Each chapter is followed by an up to date bibliography. In conclusion, this book will be a valuable asset to relativists who wish to learn clearly-stated mathematical results and to mathematicians who want to penetrate into the subtleties of general relativity, as a mathematical and physical theory. (book review)

Introductory Applications of Partial Differential Equations With Emphasis on Wave Propagation and Diffusion

CERN Document Server

Lamb, George L

1995-01-01

INTRODUCTORY APPLICATIONS OF PARTIAL DIFFERENTIAL EQUATIONS. With Emphasis on Wave Propagation and Diffusion. This is the ideal text for students and professionals who have some familiarity with partial differential equations, and who now wish to consolidate and expand their knowledge. Unlike most other texts on this topic, it interweaves prior knowledge of mathematics and physics, especially heat conduction and wave motion, into a presentation that demonstrates their interdependence. The result is a superb teaching text that reinforces the reader's understanding of both mathematics and physic

RECTC/RECTCF, 2. Order Elliptical Partial Differential Equation, Arbitrary Boundary Conditions

International Nuclear Information System (INIS)

Hackbusch, W.

1983-01-01

1 - Description of problem or function: A general linear elliptical second order partial differential equation on a rectangle with arbitrary boundary conditions is solved. 2 - Method of solution: Multi-grid iteration

Boundary value problems and partial differential equations

CERN Document Server

Powers, David L

2005-01-01

Boundary Value Problems is the leading text on boundary value problems and Fourier series. The author, David Powers, (Clarkson) has written a thorough, theoretical overview of solving boundary value problems involving partial differential equations by the methods of separation of variables. Professors and students agree that the author is a master at creating linear problems that adroitly illustrate the techniques of separation of variables used to solve science and engineering.* CD with animations and graphics of solutions, additional exercises and chapter review questions* Nearly 900 exercises ranging in difficulty* Many fully worked examples

Nonlinear elliptic partial differential equations an introduction

CERN Document Server

Le Dret, Hervé

2018-01-01

This textbook presents the essential parts of the modern theory of nonlinear partial differential equations, including the calculus of variations. After a short review of results in real and functional analysis, the author introduces the main mathematical techniques for solving both semilinear and quasilinear elliptic PDEs, and the associated boundary value problems. Key topics include infinite dimensional fixed point methods, the Galerkin method, the maximum principle, elliptic regularity, and the calculus of variations. Aimed at graduate students and researchers, this textbook contains numerous examples and exercises and provides several comments and suggestions for further study.

Mild Solutions of Neutral Stochastic Partial Functional Differential Equations

Directory of Open Access Journals (Sweden)

T. E. Govindan

2011-01-01

Full Text Available This paper studies the existence and uniqueness of a mild solution for a neutral stochastic partial functional differential equation using a local Lipschitz condition. When the neutral term is zero and even in the deterministic special case, the result obtained here appears to be new. An example is included to illustrate the theory.

Solution of partial differential equations by agent-based simulation

International Nuclear Information System (INIS)

Szilagyi, Miklos N

2014-01-01

The purpose of this short note is to demonstrate that partial differential equations can be quickly solved by agent-based simulation with high accuracy. There is no need for the solution of large systems of algebraic equations. This method is especially useful for quick determination of potential distributions and demonstration purposes in teaching electromagnetism. (letters and comments)

Second-order oriented partial-differential equations for denoising in electronic-speckle-pattern interferometry fringes.

Science.gov (United States)

Tang, Chen; Han, Lin; Ren, Hongwei; Zhou, Dongjian; Chang, Yiming; Wang, Xiaohang; Cui, Xiaolong

2008-10-01

We derive the second-order oriented partial-differential equations (PDEs) for denoising in electronic-speckle-pattern interferometry fringe patterns from two points of view. The first is based on variational methods, and the second is based on controlling diffusion direction. Our oriented PDE models make the diffusion along only the fringe orientation. The main advantage of our filtering method, based on oriented PDE models, is that it is very easy to implement compared with the published filtering methods along the fringe orientation. We demonstrate the performance of our oriented PDE models via application to two computer-simulated and experimentally obtained speckle fringes and compare with related PDE models.

A New Numerical Technique for Solving Systems Of Nonlinear Fractional Partial Differential Equations

Directory of Open Access Journals (Sweden)

Mountassir Hamdi Cherif

2017-11-01

Full Text Available In this paper, we apply an efficient method called the Aboodh decomposition method to solve systems of nonlinear fractional partial differential equations. This method is a combined form of Aboodh transform with Adomian decomposition method. The theoretical analysis of this investigated for systems of nonlinear fractional partial differential equations is calculated in the explicit form of a power series with easily computable terms. Some examples are given to shows that this method is very efficient and accurate. This method can be applied to solve others nonlinear systems problems.

Spectral methods for time dependent partial differential equations

Science.gov (United States)

Gottlieb, D.; Turkel, E.

1983-01-01

The theory of spectral methods for time dependent partial differential equations is reviewed. When the domain is periodic Fourier methods are presented while for nonperiodic problems both Chebyshev and Legendre methods are discussed. The theory is presented for both hyperbolic and parabolic systems using both Galerkin and collocation procedures. While most of the review considers problems with constant coefficients the extension to nonlinear problems is also discussed. Some results for problems with shocks are presented.

A Novel Method for Analytical Solutions of Fractional Partial Differential Equations

Directory of Open Access Journals (Sweden)

Mehmet Ali Akinlar

2013-01-01

Full Text Available A new solution technique for analytical solutions of fractional partial differential equations (FPDEs is presented. The solutions are expressed as a finite sum of a vector type functional. By employing MAPLE software, it is shown that the solutions might be extended to an arbitrary degree which makes the present method not only different from the others in the literature but also quite efficient. The method is applied to special Bagley-Torvik and Diethelm fractional differential equations as well as a more general fractional differential equation.

Plane waves and spherical means applied to partial differential equations

CERN Document Server

John, Fritz

2004-01-01

Elementary and self-contained, this heterogeneous collection of results on partial differential equations employs certain elementary identities for plane and spherical integrals of an arbitrary function, showing how a variety of results on fairly general differential equations follow from those identities. The first chapter deals with the decomposition of arbitrary functions into functions of the type of plane waves. Succeeding chapters introduce the first application of the Radon transformation and examine the solution of the initial value problem for homogeneous hyperbolic equations with con

Exp-function method for solving fractional partial differential equations.

Science.gov (United States)

Zheng, Bin

2013-01-01

We extend the Exp-function method to fractional partial differential equations in the sense of modified Riemann-Liouville derivative based on nonlinear fractional complex transformation. For illustrating the validity of this method, we apply it to the space-time fractional Fokas equation and the nonlinear fractional Sharma-Tasso-Olver (STO) equation. As a result, some new exact solutions for them are successfully established.

Differential Equations Models to Study Quorum Sensing.

Science.gov (United States)

Pérez-Velázquez, Judith; Hense, Burkhard A

2018-01-01

Mathematical models to study quorum sensing (QS) have become an important tool to explore all aspects of this type of bacterial communication. A wide spectrum of mathematical tools and methods such as dynamical systems, stochastics, and spatial models can be employed. In this chapter, we focus on giving an overview of models consisting of differential equations (DE), which can be used to describe changing quantities, for example, the dynamics of one or more signaling molecule in time and space, often in conjunction with bacterial growth dynamics. The chapter is divided into two sections: ordinary differential equations (ODE) and partial differential equations (PDE) models of QS. Rates of change are represented mathematically by derivatives, i.e., in terms of DE. ODE models allow describing changes in one independent variable, for example, time. PDE models can be used to follow changes in more than one independent variable, for example, time and space. Both types of models often consist of systems (i.e., more than one equation) of equations, such as equations for bacterial growth and autoinducer concentration dynamics. Almost from the onset, mathematical modeling of QS using differential equations has been an interdisciplinary endeavor and many of the works we revised here will be placed into their biological context.

The fifth-order partial differential equation for the description of the α + β Fermi-Pasta-Ulam model

Science.gov (United States)

Kudryashov, Nikolay A.; Volkov, Alexandr K.

2017-01-01

We study a new nonlinear partial differential equation of the fifth order for the description of perturbations in the Fermi-Pasta-Ulam mass chain. This fifth-order equation is an expansion of the Gardner equation for the description of the Fermi-Pasta-Ulam model. We use the potential of interaction between neighbouring masses with both quadratic and cubic terms. The equation is derived using the continuous limit. Unlike the previous works, we take into account higher order terms in the Taylor series expansions. We investigate the equation using the Painlevé approach. We show that the equation does not pass the Painlevé test and can not be integrated by the inverse scattering transform. We use the logistic function method and the Laurent expansion method to find travelling wave solutions of the fifth-order equation. We use the pseudospectral method for the numerical simulation of wave processes, described by the equation.

Path integral solution of linear second order partial differential equations I: the general construction

International Nuclear Information System (INIS)

LaChapelle, J.

2004-01-01

A path integral is presented that solves a general class of linear second order partial differential equations with Dirichlet/Neumann boundary conditions. Elementary kernels are constructed for both Dirichlet and Neumann boundary conditions. The general solution can be specialized to solve elliptic, parabolic, and hyperbolic partial differential equations with boundary conditions. This extends the well-known path integral solution of the Schroedinger/diffusion equation in unbounded space. The construction is based on a framework for functional integration introduced by Cartier/DeWitt-Morette

Nonlinear partial differential equations and their applications

CERN Document Server

Lions, Jacques Louis

2002-01-01

This book contains the written versions of lectures delivered since 1997 in the well-known weekly seminar on Applied Mathematics at the Collège de France in Paris, directed by Jacques-Louis Lions. It is the 14th and last of the series, due to the recent and untimely death of Professor Lions. The texts in this volume deal mostly with various aspects of the theory of nonlinear partial differential equations. They present both theoretical and applied results in many fields of growing importance such as Calculus of variations and optimal control, optimization, system theory and control, op

ON PARTIAL DIFFERENTIAL AND DIFFERENCE EQUATIONS WITH SYMMETRIES DEPENDING ON ARBITRARY FUNCTIONS

Directory of Open Access Journals (Sweden)

Giorgio Gubbiotti

2016-06-01

Full Text Available In this note we present some ideas on when Lie symmetries, both point and generalized, can depend on arbitrary functions. We show a few examples, both in partial differential and partial difference equations where this happens. Moreover we show that the infinitesimal generators of generalized symmetries depending on arbitrary functions, both for continuous and discrete equations, effectively play the role of master symmetries.

Maillet type theorem for singular first order nonlinear partial differential equations of totally characteristic type. Part II

Directory of Open Access Journals (Sweden)

Akira Shirai

2015-01-01

Full Text Available In this paper, we study the following nonlinear first order partial differential equation: \\[f(t,x,u,\\partial_t u,\\partial_x u=0\\quad\\text{with}\\quad u(0,x\\equiv 0.\\] The purpose of this paper is to determine the estimate of Gevrey order under the condition that the equation is singular of a totally characteristic type. The Gevrey order is indicated by the rate of divergence of a formal power series. This paper is a continuation of the previous papers [Convergence of formal solutions of singular first order nonlinear partial differential equations of totally characteristic type, Funkcial. Ekvac. 45 (2002, 187-208] and [Maillet type theorem for singular first order nonlinear partial differential equations of totally characteristic type, Surikaiseki Kenkyujo Kokyuroku, Kyoto University 1431 (2005, 94-106]. Especially the last-mentioned paper is regarded as part I of this paper.

Computing rates of Markov models of voltage-gated ion channels by inverting partial differential equations governing the probability density functions of the conducting and non-conducting states.

Science.gov (United States)

Tveito, Aslak; Lines, Glenn T; Edwards, Andrew G; McCulloch, Andrew

2016-07-01

Markov models are ubiquitously used to represent the function of single ion channels. However, solving the inverse problem to construct a Markov model of single channel dynamics from bilayer or patch-clamp recordings remains challenging, particularly for channels involving complex gating processes. Methods for solving the inverse problem are generally based on data from voltage clamp measurements. Here, we describe an alternative approach to this problem based on measurements of voltage traces. The voltage traces define probability density functions of the functional states of an ion channel. These probability density functions can also be computed by solving a deterministic system of partial differential equations. The inversion is based on tuning the rates of the Markov models used in the deterministic system of partial differential equations such that the solution mimics the properties of the probability density function gathered from (pseudo) experimental data as well as possible. The optimization is done by defining a cost function to measure the difference between the deterministic solution and the solution based on experimental data. By evoking the properties of this function, it is possible to infer whether the rates of the Markov model are identifiable by our method. We present applications to Markov model well-known from the literature. Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.

The Cousin problems in the viewpoint of partial differential equations

International Nuclear Information System (INIS)

Le Hung Son.

1990-01-01

In this paper we consider the Cousin problems for overdetermined systems of partial differential equations, which are generalizations of the Cauchy-Riemann system. The general methods for solving these problems are given. Applying the given methods we can solve the Cousin problems for many important systems in theoretical physics. (author). 19 refs

A Novel Method for Analytical Solutions of Fractional Partial Differential Equations

OpenAIRE

Mehmet Ali Akinlar; Muhammet Kurulay

2013-01-01

A new solution technique for analytical solutions of fractional partial differential equations (FPDEs) is presented. The solutions are expressed as a finite sum of a vector type functional. By employing MAPLE software, it is shown that the solutions might be extended to an arbitrary degree which makes the present method not only different from the others in the literature but also quite efficient. The method is applied to special Bagley-Torvik and Diethelm fractional differential equations as...

Isostable reduction with applications to time-dependent partial differential equations.

Science.gov (United States)

Wilson, Dan; Moehlis, Jeff

2016-07-01

Isostables and isostable reduction, analogous to isochrons and phase reduction for oscillatory systems, are useful in the study of nonlinear equations which asymptotically approach a stationary solution. In this work, we present a general method for isostable reduction of partial differential equations, with the potential power to reduce the dimensionality of a nonlinear system from infinity to 1. We illustrate the utility of this reduction by applying it to two different models with biological relevance. In the first example, isostable reduction of the Fokker-Planck equation provides the necessary framework to design a simple control strategy to desynchronize a population of pathologically synchronized oscillatory neurons, as might be relevant to Parkinson's disease. Another example analyzes a nonlinear reaction-diffusion equation with relevance to action potential propagation in a cardiac system.

Differential invariants in nonclassical models of hydrodynamics

Science.gov (United States)

Bublik, Vasily V.

2017-10-01

In this paper, differential invariants are used to construct solutions for equations of the dynamics of a viscous heat-conducting gas and the dynamics of a viscous incompressible fluid modified by nanopowder inoculators. To describe the dynamics of a viscous heat-conducting gas, we use the complete system of Navier—Stokes equations with allowance for heat fluxes. Mathematical description of the dynamics of liquid metals under high-energy external influences (laser radiation or plasma flow) includes, in addition to the Navier—Stokes system of an incompressible viscous fluid, also heat fluxes and processes of nonequilibrium crystallization of a deformable fluid. Differentially invariant solutions are a generalization of partially invariant solutions, and their active study for various models of continuous medium mechanics is just beginning. Differentially invariant solutions can also be considered as solutions with differential constraints; therefore, when developing them, the approaches and methods developed by the science schools of academicians N. N. Yanenko and A. F. Sidorov will be actively used. In the construction of partially invariant and differentially invariant solutions, there are overdetermined systems of differential equations that require a compatibility analysis. The algorithms for reducing such systems to involution in a finite number of steps are described by Cartan, Finikov, Kuranishi, and other authors. However, the difficultly foreseeable volume of intermediate calculations complicates their practical application. Therefore, the methods of computer algebra are actively used here, which largely helps in solving this difficult problem. It is proposed to use the constructed exact solutions as tests for formulas, algorithms and their software implementations when developing and creating numerical methods and computational program complexes. This combination of effective numerical methods, capable of solving a wide class of problems, with

Numerical solutions of ordinary and partial differential equations in the frequency domain

International Nuclear Information System (INIS)

Hazi, G.; Por, G.

1997-01-01

Numerical problems during the noise simulation in a nuclear power plant are discussed. The solutions of ordinary and partial differential equations are studied in the frequency domain. Numerical methods by the transfer function method are applied. It is shown that the correctness of the numerical methods is limited for ordinary differential equations in the frequency domain. To overcome the difficulties, step-size selection is suggested. (author)

Nonclassical Symmetries for Nonlinear Partial Differential Equations via Compatibility

International Nuclear Information System (INIS)

El-Sabbagh, Mostafa F.; Ahmad, Ali T.

2011-01-01

The determining equations for the nonclassical symmetry reductions of nonlinear partial differential equations with arbitrary order can be obtained by requiring the compatibility between the original equations and the invariant surface conditions. The (2+1)-dimensional shallow water wave equation, Boussinesq equation, and the dispersive wave equations in shallow water serve as examples illustrating how compatibility leads quickly and easily to the determining equations for their nonclassical symmetries. (general)

Using the phase-space imager to analyze partially coherent imaging systems: bright-field, phase contrast, differential interference contrast, differential phase contrast, and spiral phase contrast

Science.gov (United States)

Mehta, Shalin B.; Sheppard, Colin J. R.

2010-05-01

Various methods that use large illumination aperture (i.e. partially coherent illumination) have been developed for making transparent (i.e. phase) specimens visible. These methods were developed to provide qualitative contrast rather than quantitative measurement-coherent illumination has been relied upon for quantitative phase analysis. Partially coherent illumination has some important advantages over coherent illumination and can be used for measurement of the specimen's phase distribution. However, quantitative analysis and image computation in partially coherent systems have not been explored fully due to the lack of a general, physically insightful and computationally efficient model of image formation. We have developed a phase-space model that satisfies these requirements. In this paper, we employ this model (called the phase-space imager) to elucidate five different partially coherent systems mentioned in the title. We compute images of an optical fiber under these systems and verify some of them with experimental images. These results and simulated images of a general phase profile are used to compare the contrast and the resolution of the imaging systems. We show that, for quantitative phase imaging of a thin specimen with matched illumination, differential phase contrast offers linear transfer of specimen information to the image. We also show that the edge enhancement properties of spiral phase contrast are compromised significantly as the coherence of illumination is reduced. The results demonstrate that the phase-space imager model provides a useful framework for analysis, calibration, and design of partially coherent imaging methods.

A three operator split-step method covering a larger set of non-linear partial differential equations

Science.gov (United States)

Zia, Haider

2017-06-01

This paper describes an updated exponential Fourier based split-step method that can be applied to a greater class of partial differential equations than previous methods would allow. These equations arise in physics and engineering, a notable example being the generalized derivative non-linear Schrödinger equation that arises in non-linear optics with self-steepening terms. These differential equations feature terms that were previously inaccessible to model accurately with low computational resources. The new method maintains a 3rd order error even with these additional terms and models the equation in all three spatial dimensions and time. The class of non-linear differential equations that this method applies to is shown. The method is fully derived and implementation of the method in the split-step architecture is shown. This paper lays the mathematical ground work for an upcoming paper employing this method in white-light generation simulations in bulk material.

A Novel Partial Differential Algebraic Equation (PDAE) Solver

DEFF Research Database (Denmark)

Lim, Young-il; Chang, Sin-Chung; Jørgensen, Sten Bay

2004-01-01

For solving partial differential algebraic equations (PDAEs), the space-time conservation element/solution element (CE/SE) method is addressed in this study. The method of lines (MOL) using an implicit time integrator is compared with the CE/SE method in terms of computational efficiency, solution...... or nonlinear adsorption isotherm are solved by the two methods. The CE/SE method enforces both local and global flux conservation in space and time, and uses a simple stencil structure (two points at the previous time level and one point at the present time level). Thus, accurate and computationally...

Partial Differential Equations

CERN Document Server

1988-01-01

The volume contains a selection of papers presented at the 7th Symposium on differential geometry and differential equations (DD7) held at the Nankai Institute of Mathematics, Tianjin, China, in 1986. Most of the contributions are original research papers on topics including elliptic equations, hyperbolic equations, evolution equations, non-linear equations from differential geometry and mechanics, micro-local analysis.

RBSDE's with jumps and the related obstacle problems for integral-partial differential equations

Institute of Scientific and Technical Information of China (English)

FAN; Yulian

2006-01-01

The author proves, when the noise is driven by a Brownian motion and an independent Poisson random measure, the one-dimensional reflected backward stochastic differential equation with a stopping time terminal has a unique solution. And in a Markovian framework, the solution can provide a probabilistic interpretation for the obstacle problem for the integral-partial differential equation.

A Priori Regularity of Parabolic Partial Differential Equations

KAUST Repository

Berkemeier, Francisco

2018-05-13

In this thesis, we consider parabolic partial differential equations such as the heat equation, the Fokker-Planck equation, and the porous media equation. Our aim is to develop methods that provide a priori estimates for solutions with singular initial data. These estimates are obtained by understanding the time decay of norms of solutions. First, we derive regularity results for the heat equation by estimating the decay of Lebesgue norms. Then, we apply similar methods to the Fokker-Planck equation with suitable assumptions on the advection and diffusion. Finally, we conclude by extending our techniques to the porous media equation. The sharpness of our results is confirmed by examining known solutions of these equations. The main contribution of this thesis is the use of functional inequalities to express decay of norms as differential inequalities. These are then combined with ODE methods to deduce estimates for the norms of solutions and their derivatives.

Partial wave analysis for folded differential cross sections

Science.gov (United States)

Machacek, J. R.; McEachran, R. P.

2018-03-01

The value of modified effective range theory (MERT) and the connection between differential cross sections and phase shifts in low-energy electron scattering has long been recognized. Recent experimental techniques involving magnetically confined beams have introduced the concept of folded differential cross sections (FDCS) where the forward (θ ≤ π/2) and backward scattered (θ ≥ π/2) projectiles are unresolved, that is the value measured at the angle θ is the sum of the signal for particles scattered into the angles θ and π - θ. We have developed an alternative approach to MERT in order to analyse low-energy folded differential cross sections for positrons and electrons. This results in a simplified expression for the FDCS when it is expressed in terms of partial waves and thereby enables one to extract the first few phase shifts from a fit to an experimental FDCS at low energies. Thus, this method predicts forward and backward angle scattering (0 to π) using only experimental FDCS data and can be used to determine the total elastic cross section solely from experimental results at low-energy, which are limited in angular range.

Partial differential equations II elements of the modern theory equations with constant coefficients

CERN Document Server

Shubin, M

1994-01-01

This book, the first printing of which was published as Volume 31 of the Encyclopaedia of Mathematical Sciences, contains a survey of the modern theory of general linear partial differential equations and a detailed review of equations with constant coefficients. Readers will be interested in an introduction to microlocal analysis and its applications including singular integral operators, pseudodifferential operators, Fourier integral operators and wavefronts, a survey of the most important results about the mixed problem for hyperbolic equations, a review of asymptotic methods including short wave asymptotics, the Maslov canonical operator and spectral asymptotics, a detailed description of the applications of distribution theory to partial differential equations with constant coefficients including numerous interesting special topics.

4th International Conference on Particle Systems and Partial Differential Equations

CERN Document Server

Soares, Ana

2017-01-01

'This book addresses mathematical problems motivated by various applications in physics, engineering, chemistry and biology. It gathers the lecture notes from the mini-course presented by Jean-Christophe Mourrat on the construction of the various stochastic “basic” terms involved in the formulation of the dynamic Ö4  theory in three space dimensions, as well as selected contributions presented at the fourth meeting on Particle Systems and PDEs, which was held at the University of Minho’s Centre of Mathematics in December 2015. The purpose of the conference was to bring together prominent researchers working in the fields of particle systems and partial differential equations, offering them a forum to present their recent results and discuss their topics of expertise. The meeting was also intended to present to a vast and varied public, including young researchers, the area of interacting particle systems, its underlying motivation, and its relation to partial differential equations.  The book w...

A boundary value approach for solving three-dimensional elliptic and hyperbolic partial differential equations.

Science.gov (United States)

Biala, T A; Jator, S N

2015-01-01

In this article, the boundary value method is applied to solve three dimensional elliptic and hyperbolic partial differential equations. The partial derivatives with respect to two of the spatial variables (y, z) are discretized using finite difference approximations to obtain a large system of ordinary differential equations (ODEs) in the third spatial variable (x). Using interpolation and collocation techniques, a continuous scheme is developed and used to obtain discrete methods which are applied via the Block unification approach to obtain approximations to the resulting large system of ODEs. Several test problems are investigated to elucidate the solution process.

Collage-based approaches for elliptic partial differential equations inverse problems

Science.gov (United States)

Yodzis, Michael; Kunze, Herb

2017-01-01

The collage method for inverse problems has become well-established in the literature in recent years. Initial work developed a collage theorem, based upon Banach's fixed point theorem, for treating inverse problems for ordinary differential equations (ODEs). Amongst the subsequent work was a generalized collage theorem, based upon the Lax-Milgram representation theorem, useful for treating inverse problems for elliptic partial differential equations (PDEs). Each of these two different approaches can be applied to elliptic PDEs in one space dimension. In this paper, we explore and compare how the two different approaches perform for the estimation of the diffusivity for a steady-state heat equation.

Nonlinear partial differential equations of second order

CERN Document Server

Dong, Guangchang

1991-01-01

This book addresses a class of equations central to many areas of mathematics and its applications. Although there is no routine way of solving nonlinear partial differential equations, effective approaches that apply to a wide variety of problems are available. This book addresses a general approach that consists of the following: Choose an appropriate function space, define a family of mappings, prove this family has a fixed point, and study various properties of the solution. The author emphasizes the derivation of various estimates, including a priori estimates. By focusing on a particular approach that has proven useful in solving a broad range of equations, this book makes a useful contribution to the literature.

Analytical solutions for coupling fractional partial differential equations with Dirichlet boundary conditions

Science.gov (United States)

Ding, Xiao-Li; Nieto, Juan J.

2017-11-01

In this paper, we consider the analytical solutions of coupling fractional partial differential equations (FPDEs) with Dirichlet boundary conditions on a finite domain. Firstly, the method of successive approximations is used to obtain the analytical solutions of coupling multi-term time fractional ordinary differential equations. Then, the technique of spectral representation of the fractional Laplacian operator is used to convert the coupling FPDEs to the coupling multi-term time fractional ordinary differential equations. By applying the obtained analytical solutions to the resulting multi-term time fractional ordinary differential equations, the desired analytical solutions of the coupling FPDEs are given. Our results are applied to derive the analytical solutions of some special cases to demonstrate their applicability.

Unconditionally stable difference methods for delay partial differential equations

OpenAIRE

Huang, Chengming; Vandewalle, Stefan

2012-01-01

This paper is concerned with the numerical solution of parabolic partial differential equations with time-delay. We focus in particular on the delay dependent stability analysis of difference methods that use a non-constrained mesh, i.e., the time step-size is not required to be a submultiple of the delay. We prove that the fully discrete system unconditionally preserves the delay dependent asymptotic stability of the linear test problem under consideration, when the following discretizati...

Multivariate Padé Approximation for Solving Nonlinear Partial Differential Equations of Fractional Order

Directory of Open Access Journals (Sweden)

Veyis Turut

2013-01-01

Full Text Available Two tecHniques were implemented, the Adomian decomposition method (ADM and multivariate Padé approximation (MPA, for solving nonlinear partial differential equations of fractional order. The fractional derivatives are described in Caputo sense. First, the fractional differential equation has been solved and converted to power series by Adomian decomposition method (ADM, then power series solution of fractional differential equation was put into multivariate Padé series. Finally, numerical results were compared and presented in tables and figures.

Variational and potential formulation for stochastic partial differential equations

International Nuclear Information System (INIS)

Munoz S, A G; Ojeda, J; Sierra D, P; Soldovieri, T

2006-01-01

Recently there has been interest in finding a potential formulation for stochastic partial differential equations (SPDEs). The rationale behind this idea lies in obtaining all the dynamical information of the system under study from one single expression. In this letter we formally provide a general Lagrangian formalism for SPDEs using the Hojman et al method. We show that it is possible to write the corresponding effective potential starting from an s-equivalent Lagrangian, and that this potential is able to reproduce all the dynamics of the system once a special differential operator has been applied. This procedure can be used to study the complete time evolution and spatial inhomogeneities of the system under consideration, and is also suitable for the statistical mechanics description of the problem. (letter to the editor)

Inverse problems for partial differential equations

CERN Document Server

Isakov, Victor

2017-01-01

This third edition expands upon the earlier edition by adding nearly 40 pages of new material reflecting the analytical and numerical progress in inverse problems in last 10 years. As in the second edition, the emphasis is on new ideas and methods rather than technical improvements. These new ideas include use of the stationary phase method in the two-dimensional elliptic problems and of multi frequencies\\temporal data to improve stability and numerical resolution. There are also numerous corrections and improvements of the exposition throughout. This book is intended for mathematicians working with partial differential equations and their applications, physicists, geophysicists, and financial, electrical, and mechanical engineers involved with nondestructive evaluation, seismic exploration, remote sensing, and various kinds of tomography. Review of the second edition: "The first edition of this excellent book appeared in 1998 and became a standard reference for everyone interested in analysis and numerics of...

The pseudo-compartment method for coupling partial differential equation and compartment-based models of diffusion.

Science.gov (United States)

Yates, Christian A; Flegg, Mark B

2015-05-06

Spatial reaction-diffusion models have been employed to describe many emergent phenomena in biological systems. The modelling technique most commonly adopted in the literature implements systems of partial differential equations (PDEs), which assumes there are sufficient densities of particles that a continuum approximation is valid. However, owing to recent advances in computational power, the simulation and therefore postulation, of computationally intensive individual-based models has become a popular way to investigate the effects of noise in reaction-diffusion systems in which regions of low copy numbers exist. The specific stochastic models with which we shall be concerned in this manuscript are referred to as 'compartment-based' or 'on-lattice'. These models are characterized by a discretization of the computational domain into a grid/lattice of 'compartments'. Within each compartment, particles are assumed to be well mixed and are permitted to react with other particles within their compartment or to transfer between neighbouring compartments. Stochastic models provide accuracy, but at the cost of significant computational resources. For models that have regions of both low and high concentrations, it is often desirable, for reasons of efficiency, to employ coupled multi-scale modelling paradigms. In this work, we develop two hybrid algorithms in which a PDE in one region of the domain is coupled to a compartment-based model in the other. Rather than attempting to balance average fluxes, our algorithms answer a more fundamental question: 'how are individual particles transported between the vastly different model descriptions?' First, we present an algorithm derived by carefully redefining the continuous PDE concentration as a probability distribution. While this first algorithm shows very strong convergence to analytical solutions of test problems, it can be cumbersome to simulate. Our second algorithm is a simplified and more efficient implementation of

A new RBF-Trefftz meshless method for partial differential equations

International Nuclear Information System (INIS)

Cao Leilei; Zhao Ning; Qin Qinghua

2010-01-01

Based on the radial basis functions (RBF) and T-Trefftz solution, this paper presents a new meshless method for numerically solving various partial differential equation systems. First, the analog equation method (AEM) is used to convert the original patial differential equation to an equivalent Poisson's equation. Then, the radial basis functions (RBF) are employed to approxiamate the inhomogeneous term, while the homogeneous solution is obtained by linear combination of a set of T-Trefftz solutions. The present scheme, named RBF-Trefftz has the advantage over the fundamental solution (MFS) method due to the use of nonsingular T-Trefftz solution rather than singular fundamental solutions, so it does not require the artificial boundary. The application and efficiency of the proposed method are validated through several examples which include different type of differential equations, such as Laplace equation, Hellmholtz equation, convectin-diffusion equation and time-dependent equation.

The generalized tanh method to obtain exact solutions of nonlinear partial differential equation

OpenAIRE

Gómez, César

2007-01-01

In this paper, we present the generalized tanh method to obtain exact solutions of nonlinear partial differential equations, and we obtain solitons and exact solutions of some important equations of the mathematical physics.

Rasch-Master's Partial Credit Model in the assessment of children's creativity in drawings.

Science.gov (United States)

Nakano, Tatiana de Cássia; Primi, Ricardo

2014-01-01

The purpose of the present study was to use the Partial Credit Model to study the factors of the Test of Creativity in Children and identify which characteristics of the creative person would be more effective to differentiate subjects according to their ability level. A sample of 1426 students from first to eighth grades answered the instrument. The Partial Credits model was used to estimate the ability of the subjects and item difficulties on a common scale for each of the four factors, indicating which items required a higher level of creativity to be scored and will differentiate the more creative individuals. The results demonstrated that the greater part of the characteristics showed good fit indices, with values between 0.80 and 1.30 both infit and outfit, indicating a response pattern consistent with the model. The characteristics of Unusual Perspective, Expression of Emotion and Originality have been identified as better predictors of creative performance because requires greater ability level (usually above two standard deviation). These results may be used in the future development of an instrument's reduced form or simplification of the current correction model.

Partial differential equations and boundary-value problems with applications

CERN Document Server

Pinsky, Mark A

2011-01-01

Building on the basic techniques of separation of variables and Fourier series, the book presents the solution of boundary-value problems for basic partial differential equations: the heat equation, wave equation, and Laplace equation, considered in various standard coordinate systems-rectangular, cylindrical, and spherical. Each of the equations is derived in the three-dimensional context; the solutions are organized according to the geometry of the coordinate system, which makes the mathematics especially transparent. Bessel and Legendre functions are studied and used whenever appropriate th

Asymptotic behavior of solutions of diffusion-like partial differential equations invariant to a family of affine groups

International Nuclear Information System (INIS)

Dresner, L.

1990-07-01

This report deals with the asymptotic behavior of certain solutions of partial differential equations in one dependent and two independent variables (call them c, z, and t, respectively). The partial differential equations are invariant to one-parameter families of one-parameter affine groups of the form: c' = λ α c, t' = λ β t, z' = λz, where λ is the group parameter that labels the individual transformations and α and β are parameters that label groups of the family. The parameters α and β are connected by a linear relation, Mα + Nβ = L, where M, N, and L are numbers determined by the structure of the partial differential equation. It is shown that when L/M and N/M are L/M t -N/M for large z or small t. Some practical applications of this result are discussed. 8 refs

Functional analytic methods in complex analysis and applications to partial differential equations

International Nuclear Information System (INIS)

Mshimba, A.S.A.; Tutschke, W.

1990-01-01

The volume contains 24 lectures given at the Workshop on Functional Analytic Methods in Complex Analysis and Applications to Partial Differential Equations held in Trieste, Italy, between 8-19 February 1988, at the ICTP. A separate abstract was prepared for each of these lectures. Refs and figs

Constrained Optimization and Optimal Control for Partial Differential Equations

CERN Document Server

Leugering, Günter; Griewank, Andreas

2012-01-01

This special volume focuses on optimization and control of processes governed by partial differential equations. The contributors are mostly participants of the DFG-priority program 1253: Optimization with PDE-constraints which is active since 2006. The book is organized in sections which cover almost the entire spectrum of modern research in this emerging field. Indeed, even though the field of optimal control and optimization for PDE-constrained problems has undergone a dramatic increase of interest during the last four decades, a full theory for nonlinear problems is still lacking. The cont

Superdiffusions and positive solutions of nonlinear partial differential equations

CERN Document Server

Dynkin, E B

2004-01-01

This book is devoted to the applications of probability theory to the theory of nonlinear partial differential equations. More precisely, it is shown that all positive solutions for a class of nonlinear elliptic equations in a domain are described in terms of their traces on the boundary of the domain. The main probabilistic tool is the theory of superdiffusions, which describes a random evolution of a cloud of particles. A substantial enhancement of this theory is presented that can be of interest for everybody who works on applications of probabilistic methods to mathematical analysis.

Malliavin Calculus With Applications to Stochastic Partial Differential Equations

CERN Document Server

Sanz-Solé, Marta

2005-01-01

Developed in the 1970s to study the existence and smoothness of density for the probability laws of random vectors, Malliavin calculus--a stochastic calculus of variation on the Wiener space--has proven fruitful in many problems in probability theory, particularly in probabilistic numerical methods in financial mathematics.This book presents applications of Malliavin calculus to the analysis of probability laws of solutions to stochastic partial differential equations driven by Gaussian noises that are white in time and coloured in space. The first five chapters introduce the calculus itself

Nonlinear perturbations of systems of partial differential equations with constant coefficients

Directory of Open Access Journals (Sweden)

Carmen J. Vanegas

2000-01-01

Full Text Available In this article, we show the existence of solutions to boundary-value problems, consisting of nonlinear systems of partial differential equations with constant coefficients. For this purpose, we use the right inverse of an associated operator and a fix point argument. As illustrations, we apply this method to Helmholtz equations and to second order systems of elliptic equations.

Magnetic Evidence for a Partially Differentiated Carbonaceous Chondrite Parent Body and Possible Implications for Asteroid 21 Lutetia

Science.gov (United States)

Weiss, Benjamin; Carporzen, L.; Elkins-Tanton, L.; Shuster, D. L.; Ebel, D. S.; Gattacceca, J.; Binzel, R. P.

2010-10-01

The origin of remanent magnetization in the CV carbonaceous chondrite Allende has been a longstanding mystery. The possibility of a core dynamo like that known for achondrite parent bodies has been discounted because chondrite parent bodies are assumed to be undifferentiated. Here we report that Allende's magnetization was acquired over several million years (Ma) during metasomatism on the parent planetesimal in a > 20 microtesla field 8-9 Ma after solar system formation. This field was present too recently and directionally stable for too long to have been the generated by the protoplanetary disk or young Sun. The field intensity is in the range expected for planetesimal core dynamos (Weiss et al. 2010), suggesting that CV chondrites are derived from the outer, unmelted layer of a partially differentiated body with a convecting metallic core (Elkins-Tanton et al. 2010). This suggests that asteroids with differentiated interiors could be present today but masked under chondritic surfaces. In fact, CV chondrites are spectrally similar to many members of the Eos asteroid family whose spectral diversity has been interpreted as evidence for a partially differentiated parent asteroid (Mothe-Diniz et al. 2008). CV chondrite spectral and polarimetric data also resemble those of asteroid 21 Lutetia (e.g., Belskaya et al. 2010), recently encountered by the Rosetta spacecraft. Ground-based measurements of Lutetia indicate a high density of 2.4-5.1 g cm-3 (Drummond et al. 2010), while radar data seem to rule out a metallic surface composition (Shepard et al. 2008). If Rosetta spacecraft measurements confirm a high density and a CV-like surface composition for Lutetia, then we propose Lutetia may be an example of a partially differentiated carbonaceous chondrite parent body. Regardless, the very existence of primitive achondrites, which contain evidence of both relict chondrules and partial melting, are prima facie evidence for the formation of partially differentiated bodies.

A Simple Differential Modulation Scheme for Quasi-Orthogonal Space-Time Block Codes with Partial Transmit Diversity

Directory of Open Access Journals (Sweden)

Lingyang Song

2007-04-01

Full Text Available We report a simple differential modulation scheme for quasi-orthogonal space-time block codes. A new class of quasi-orthogonal coding structures that can provide partial transmit diversity is presented for various numbers of transmit antennas. Differential encoding and decoding can be simplified for differential Alamouti-like codes by grouping the signals in the transmitted matrix and decoupling the detection of data symbols, respectively. The new scheme can achieve constant amplitude of transmitted signals, and avoid signal constellation expansion; in addition it has a linear signal detector with very low complexity. Simulation results show that these partial-diversity codes can provide very useful results at low SNR for current communication systems. Extension to more than four transmit antennas is also considered.

ICM: an Integrated Compartment Method for numerically solving partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Yeh, G.T.

1981-05-01

An integrated compartment method (ICM) is proposed to construct a set of algebraic equations from a system of partial differential equations. The ICM combines the utility of integral formulation of finite element approach, the simplicity of interpolation of finite difference approximation, and the flexibility of compartment analyses. The integral formulation eases the treatment of boundary conditions, in particular, the Neumann-type boundary conditions. The simplicity of interpolation provides great economy in computation. The flexibility of discretization with irregular compartments of various shapes and sizes offers advantages in resolving complex boundaries enclosing compound regions of interest. The basic procedures of ICM are first to discretize the region of interest into compartments, then to apply three integral theorems of vectors to transform the volume integral to the surface integral, and finally to use interpolation to relate the interfacial values in terms of compartment values to close the system. The Navier-Stokes equations are used as an example of how to derive the corresponding ICM alogrithm for a given set of partial differential equations. Because of the structure of the algorithm, the basic computer program remains the same for cases in one-, two-, or three-dimensional problems.

Analytical Solutions for Systems of Singular Partial Differential-Algebraic Equations

Directory of Open Access Journals (Sweden)

U. Filobello-Nino

2015-01-01

Full Text Available This paper proposes power series method (PSM in order to find solutions for singular partial differential-algebraic equations (SPDAEs. We will solve three examples to show that PSM method can be used to search for analytical solutions of SPDAEs. What is more, we will see that, in some cases, Padé posttreatment, besides enlarging the domain of convergence, may be employed in order to get the exact solution from the truncated series solutions of PSM.

Solutions to an advanced functional partial differential equation of the pantograph type.

Science.gov (United States)

Zaidi, Ali A; Van Brunt, B; Wake, G C

2015-07-08

A model for cells structured by size undergoing growth and division leads to an initial boundary value problem that involves a first-order linear partial differential equation with a functional term. Here, size can be interpreted as DNA content or mass. It has been observed experimentally and shown analytically that solutions for arbitrary initial cell distributions are asymptotic as time goes to infinity to a certain solution called the steady size distribution. The full solution to the problem for arbitrary initial distributions, however, is elusive owing to the presence of the functional term and the paucity of solution techniques for such problems. In this paper, we derive a solution to the problem for arbitrary initial cell distributions. The method employed exploits the hyperbolic character of the underlying differential operator, and the advanced nature of the functional argument to reduce the problem to a sequence of simple Cauchy problems. The existence of solutions for arbitrary initial distributions is established along with uniqueness. The asymptotic relationship with the steady size distribution is established, and because the solution is known explicitly, higher-order terms in the asymptotics can be readily obtained.

Introduction to partial differential equations from Fourier series to boundary-value problems

CERN Document Server

Broman, Arne

2010-01-01

This well-written, advanced-level text introduces students to Fourier analysis and some of its applications. The self-contained treatment covers Fourier series, orthogonal systems, Fourier and Laplace transforms, Bessel functions, and partial differential equations of the first and second orders. Over 260 exercises with solutions reinforce students' grasp of the material. 1970 edition.

A direct algebraic method applied to obtain complex solutions of some nonlinear partial differential equations

International Nuclear Information System (INIS)

Zhang Huiqun

2009-01-01

By using some exact solutions of an auxiliary ordinary differential equation, a direct algebraic method is described to construct the exact complex solutions for nonlinear partial differential equations. The method is implemented for the NLS equation, a new Hamiltonian amplitude equation, the coupled Schrodinger-KdV equations and the Hirota-Maccari equations. New exact complex solutions are obtained.

On new classes of solutions of nonlinear partial differential equations in the form of convergent special series

Science.gov (United States)

Filimonov, M. Yu.

2017-12-01

The method of special series with recursively calculated coefficients is used to solve nonlinear partial differential equations. The recurrence of finding the coefficients of the series is achieved due to a special choice of functions, in powers of which the solution is expanded in a series. We obtain a sequence of linear partial differential equations to find the coefficients of the series constructed. In many cases, one can deal with a sequence of linear ordinary differential equations. We construct classes of solutions in the form of convergent series for a certain class of nonlinear evolution equations. A new class of solutions of generalized Boussinesque equation with an arbitrary function in the form of a convergent series is constructed.

Asymptotic problems for stochastic partial differential equations

Science.gov (United States)

Salins, Michael

Stochastic partial differential equations (SPDEs) can be used to model systems in a wide variety of fields including physics, chemistry, and engineering. The main SPDEs of interest in this dissertation are the semilinear stochastic wave equations which model the movement of a material with constant mass density that is exposed to both determinstic and random forcing. Cerrai and Freidlin have shown that on fixed time intervals, as the mass density of the material approaches zero, the solutions of the stochastic wave equation converge uniformly to the solutions of a stochastic heat equation, in probability. This is called the Smoluchowski-Kramers approximation. In Chapter 2, we investigate some of the multi-scale behaviors that these wave equations exhibit. In particular, we show that the Freidlin-Wentzell exit place and exit time asymptotics for the stochastic wave equation in the small noise regime can be approximated by the exit place and exit time asymptotics for the stochastic heat equation. We prove that the exit time and exit place asymptotics are characterized by quantities called quasipotentials and we prove that the quasipotentials converge. We then investigate the special case where the equation has a gradient structure and show that we can explicitly solve for the quasipotentials, and that the quasipotentials for the heat equation and wave equation are equal. In Chapter 3, we study the Smoluchowski-Kramers approximation in the case where the material is electrically charged and exposed to a magnetic field. Interestingly, if the system is frictionless, then the Smoluchowski-Kramers approximation does not hold. We prove that the Smoluchowski-Kramers approximation is valid for systems exposed to both a magnetic field and friction. Notably, we prove that the solutions to the second-order equations converge to the solutions of the first-order equation in an Lp sense. This strengthens previous results where convergence was proved in probability.

Use of fast Fourier transforms for solving partial differential equations in physics

CERN Document Server

Le Bail, R C

1972-01-01

The use of fast Fourier techniques for the direct solution of an important class of elliptic, parabolic, and hyperbolic partial differential equations in two dimensions is described. Extensions to higher-order and higher-dimension equations as well as to integrodifferential equations are presented, and several numerical examples with their resulting precision and timing are reported. (12 refs).

BOOK REVIEW: Partial Differential Equations in General Relativity

Science.gov (United States)

Halburd, Rodney G.

2008-11-01

Although many books on general relativity contain an overview of the relevant background material from differential geometry, very little attention is usually paid to background material from the theory of differential equations. This is understandable in a first course on relativity but it often limits the kinds of problems that can be studied rigorously. Einstein's field equations lie at the heart of general relativity. They are a system of partial differential equations (PDEs) relating the curvature of spacetime to properties of matter. A central part of most problems in general relativity is to extract information about solutions of these equations. Most standard texts achieve this by studying exact solutions or numerical and analytical approximations. In the book under review, Alan Rendall emphasises the role of rigorous qualitative methods in general relativity. There has long been a need for such a book, giving a broad overview of the relevant background from the theory of partial differential equations, and not just from differential geometry. It should be noted that the book also covers the basic theory of ordinary differential equations. Although there are many good books on the rigorous theory of PDEs, methods related to the Einstein equations deserve special attention, not only because of the complexity and importance of these equations, but because these equations do not fit into any of the standard classes of equations (elliptic, parabolic, hyperbolic) that one typically encounters in a course on PDEs. Even specifying exactly what ones means by a Cauchy problem in general relativity requires considerable care. The main problem here is that the manifold on which the solution is defined is determined by the solution itself. This means that one does not simply define data on a submanifold. Rendall's book gives a good overview of applications and results from the qualitative theory of PDEs to general relativity. It would be impossible to give detailed

The Spectral/hp-Finite Element Method for Partial Differential Equations

DEFF Research Database (Denmark)

Engsig-Karup, Allan Peter

2009-01-01

dimensions. In the course the chosen programming environment is Matlab, however, this is by no means a necessary requirement. The mathematical level needed to grasp the details of this set of notes requires an elementary background in mathematical analysis and linear algebra. Each chapter is supplemented......This set of lecture notes provides an elementary introduction to both the classical Finite Element Method (FEM) and the extended Spectral/$hp$-Finite Element Method for solving Partial Differential Equations (PDEs). Many problems in science and engineering can be formulated mathematically...

Stochastic Partial Differential Equation Solver for Hydroacoustic Modeling: Improvements to Paracousti Sound Propagation Solver

Science.gov (United States)

Preston, L. A.

2017-12-01

Marine hydrokinetic (MHK) devices offer a clean, renewable alternative energy source for the future. Responsible utilization of MHK devices, however, requires that the effects of acoustic noise produced by these devices on marine life and marine-related human activities be well understood. Paracousti is a 3-D full waveform acoustic modeling suite that can accurately propagate MHK noise signals in the complex bathymetry found in the near-shore to open ocean environment and considers real properties of the seabed, water column, and air-surface interface. However, this is a deterministic simulation that assumes the environment and source are exactly known. In reality, environmental and source characteristics are often only known in a statistical sense. Thus, to fully characterize the expected noise levels within the marine environment, this uncertainty in environmental and source factors should be incorporated into the acoustic simulations. One method is to use Monte Carlo (MC) techniques where simulation results from a large number of deterministic solutions are aggregated to provide statistical properties of the output signal. However, MC methods can be computationally prohibitive since they can require tens of thousands or more simulations to build up an accurate representation of those statistical properties. An alternative method, using the technique of stochastic partial differential equations (SPDE), allows computation of the statistical properties of output signals at a small fraction of the computational cost of MC. We are developing a SPDE solver for the 3-D acoustic wave propagation problem called Paracousti-UQ to help regulators and operators assess the statistical properties of environmental noise produced by MHK devices. In this presentation, we present the SPDE method and compare statistical distributions of simulated acoustic signals in simple models to MC simulations to show the accuracy and efficiency of the SPDE method. Sandia National Laboratories

Rail-to-rail low-power fully differential OTA utilizing adaptive biasing and partial feedback

DEFF Research Database (Denmark)

Tuan Vu, Cao; Wisland, Dag T.; Lande, Tor Sverre

A fully differential rail-to-rail Operational Transconductance Amplifier (OTA) with improved DC-gain and reduced power consumption is proposed in this paper. By using the adaptive biasing circuit and two differential inputs, a low stand-by current can be obtained together with reduced power...... consumption. The DC-gain of the proposed OTA is improved by adding a partial feedback loop. A Common-Mode Feedback (CMFB) circuit is required for fully differential rail-to-rail operation. Simulations show that the OTA topology has a low stand-by power consumption of 96μW and a high FoM of 3.84 [(V...

Analytical solutions for systems of partial differential-algebraic equations.

Science.gov (United States)

Benhammouda, Brahim; Vazquez-Leal, Hector

2014-01-01

This work presents the application of the power series method (PSM) to find solutions of partial differential-algebraic equations (PDAEs). Two systems of index-one and index-three are solved to show that PSM can provide analytical solutions of PDAEs in convergent series form. What is more, we present the post-treatment of the power series solutions with the Laplace-Padé (LP) resummation method as a useful strategy to find exact solutions. The main advantage of the proposed methodology is that the procedure is based on a few straightforward steps and it does not generate secular terms or depends of a perturbation parameter.

Partially composite Higgs models

DEFF Research Database (Denmark)

Alanne, Tommi; Buarque Franzosi, Diogo; Frandsen, Mads T.

2018-01-01

We study the phenomenology of partially composite-Higgs models where electroweak symmetry breaking is dynamically induced, and the Higgs is a mixture of a composite and an elementary state. The models considered have explicit realizations in terms of gauge-Yukawa theories with new strongly...... interacting fermions coupled to elementary scalars and allow for a very SM-like Higgs state. We study constraints on their parameter spaces from vacuum stability and perturbativity as well as from LHC results and find that requiring vacuum stability up to the compositeness scale already imposes relevant...... constraints. A small part of parameter space around the classically conformal limit is stable up to the Planck scale. This is however already strongly disfavored by LHC results. in different limits, the models realize both (partially) composite-Higgs and (bosonic) technicolor models and a dynamical extension...

A comparison of numerical solutions of partial differential equations with probabilistic and possibilistic parameters for the quantification of uncertainty in subsurface solute transport.

Science.gov (United States)

Zhang, Kejiang; Achari, Gopal; Li, Hua

2009-11-03

Traditionally, uncertainty in parameters are represented as probabilistic distributions and incorporated into groundwater flow and contaminant transport models. With the advent of newer uncertainty theories, it is now understood that stochastic methods cannot properly represent non random uncertainties. In the groundwater flow and contaminant transport equations, uncertainty in some parameters may be random, whereas those of others may be non random. The objective of this paper is to develop a fuzzy-stochastic partial differential equation (FSPDE) model to simulate conditions where both random and non random uncertainties are involved in groundwater flow and solute transport. Three potential solution techniques namely, (a) transforming a probability distribution to a possibility distribution (Method I) then a FSPDE becomes a fuzzy partial differential equation (FPDE), (b) transforming a possibility distribution to a probability distribution (Method II) and then a FSPDE becomes a stochastic partial differential equation (SPDE), and (c) the combination of Monte Carlo methods and FPDE solution techniques (Method III) are proposed and compared. The effects of these three methods on the predictive results are investigated by using two case studies. The results show that the predictions obtained from Method II is a specific case of that got from Method I. When an exact probabilistic result is needed, Method II is suggested. As the loss or gain of information during a probability-possibility (or vice versa) transformation cannot be quantified, their influences on the predictive results is not known. Thus, Method III should probably be preferred for risk assessments.

Paleomagnetic evidence for a partially differentiated H chondrite parent planetesimal

Science.gov (United States)

Bryson, J. F. J.; Weiss, B. P.; Scholl, A.; Getzin, B. L.; Abrahams, J. N. H.; Nimmo, F.

2016-12-01

The texture, composition and ages of chondrites have all been used to argue that the parent bodies of these meteorites did not undergo planetary differentiation. Without a core, these planetesimals could not have generated planetary magnetic fields, hence chondrites are predicted to be unmagnetized. Here, we test this hypothesis by applying synchrotron x-ray microscopy to the metallic melt veins in the metamorphosed H chondrite breccia Portales Valley. We find that tetrataenite nanostructures in these veins are uniformly magnetized, suggesting that the H chondrite parent body generated a stable, 10 µT ancient field. We also performed alternating field (AF) demagnetization on bulk silicate-rich portions of Portales Valley, finding that both the large grain size of the metal in these subsamples and the presence of tetrataenite hinder the reliable interpretation of these measurements. Based on 40Ar/39Ar dating and the metallographic cooling rate, we propose that this field inferred from x-ray microscopy was generated 100 Myr after solar system formation and lasted >5 Myr. These properties are consistent with a dynamo field generated by core solidification, implying that the H chondrite parent body was partially differentiated. This conclusion is supported by our analyses of the H4 chondrite Forest Vale, which show that H chondrite magnetization is unlikely to be a relic signature of early nebular or solar wind fields (Getzin et al., this meeting; Oran et al., this meeting). We propose that partial differentiation could result form prolonged accretion over millions of years, possibly in two stages. In this scenario, the earliest accreted material melted from the radioactive decay of abundant 26Al, forming a core and rocky achondritic mantle, while the later accreted material was less metamorphosed, forming an undifferentiated crust. We demonstrate that, with the inclusion of an insulating regolith, the thermal evolution of such a body is consistent with the measured

Convergence of method of lines approximations to partial differential equations

International Nuclear Information System (INIS)

Verwer, J.G.; Sanz-Serna, J.M.

1984-01-01

Many existing numerical schemes for evolutionary problems in partial differential equations (PDEs) can be viewed as method of lines (MOL) schemes. This paper treats the convergence of one-step MOL schemes. The main purpose is to set up a general framework for a convergence analysis applicable to nonlinear problems. The stability materials for this framework are taken from the field of nonlinear stiff ODEs. In this connection, important concepts are the logarithmic matrix norm and C-stability. A nonlinear parabolic equation and the cubic Schroedinger equation are used for illustrating the ideas. (Auth.)

Numerical methods for stochastic partial differential equations with white noise

CERN Document Server

Zhang, Zhongqiang

2017-01-01

This book covers numerical methods for stochastic partial differential equations with white noise using the framework of Wong-Zakai approximation. The book begins with some motivational and background material in the introductory chapters and is divided into three parts. Part I covers numerical stochastic ordinary differential equations. Here the authors start with numerical methods for SDEs with delay using the Wong-Zakai approximation and finite difference in time. Part II covers temporal white noise. Here the authors consider SPDEs as PDEs driven by white noise, where discretization of white noise (Brownian motion) leads to PDEs with smooth noise, which can then be treated by numerical methods for PDEs. In this part, recursive algorithms based on Wiener chaos expansion and stochastic collocation methods are presented for linear stochastic advection-diffusion-reaction equations. In addition, stochastic Euler equations are exploited as an application of stochastic collocation methods, where a numerical compa...

The Dif Identification in Constructed Response Items Using Partial Credit Model

OpenAIRE

Heri Retnawati

2017-01-01

The study was to identify the load, the type and the significance of differential item functioning (DIF) in constructed response item using the partial credit model (PCM). The data in the study were the students’ instruments and the students’ responses toward the PISA-like test items that had been completed by 386 ninth grade students and 460 tenth grade students who had been about 15 years old in the Province of Yogyakarta Special Region in Indonesia. The analysis toward the item characteris...

Elliptic–hyperbolic partial differential equations a mini-course in geometric and quasilinear methods

CERN Document Server

Otway, Thomas H

2015-01-01

This text is a concise introduction to the partial differential equations which change from elliptic to hyperbolic type across a smooth hypersurface of their domain. These are becoming increasingly important in diverse sub-fields of both applied mathematics and engineering, for example:   • The heating of fusion plasmas by electromagnetic waves • The behaviour of light near a caustic • Extremal surfaces in the space of special relativity • The formation of rapids; transonic and multiphase fluid flow • The dynamics of certain models for elastic structures • The shape of industrial surfaces such as windshields and airfoils • Pathologies of traffic flow • Harmonic fields in extended projective space   They also arise in models for the early universe, for cosmic acceleration, and for possible violation of causality in the interiors of certain compact stars. Within the past 25 years, they have become central to the isometric embedding of Riemannian manifolds and the prescription of Gauss curvatur...

Partial Differential Equations and Solitary Waves Theory

CERN Document Server

Wazwaz, Abdul-Majid

2009-01-01

"Partial Differential Equations and Solitary Waves Theory" is a self-contained book divided into two parts: Part I is a coherent survey bringing together newly developed methods for solving PDEs. While some traditional techniques are presented, this part does not require thorough understanding of abstract theories or compact concepts. Well-selected worked examples and exercises shall guide the reader through the text. Part II provides an extensive exposition of the solitary waves theory. This part handles nonlinear evolution equations by methods such as Hirota’s bilinear method or the tanh-coth method. A self-contained treatment is presented to discuss complete integrability of a wide class of nonlinear equations. This part presents in an accessible manner a systematic presentation of solitons, multi-soliton solutions, kinks, peakons, cuspons, and compactons. While the whole book can be used as a text for advanced undergraduate and graduate students in applied mathematics, physics and engineering, Part II w...

A higher order numerical method for time fractional partial differential equations with nonsmooth data

Science.gov (United States)

Xing, Yanyuan; Yan, Yubin

2018-03-01

Gao et al. [11] (2014) introduced a numerical scheme to approximate the Caputo fractional derivative with the convergence rate O (k 3 - α), 0 equation is sufficiently smooth, Lv and Xu [20] (2016) proved by using energy method that the corresponding numerical method for solving time fractional partial differential equation has the convergence rate O (k 3 - α), 0 equation has low regularity and in this case the numerical method fails to have the convergence rate O (k 3 - α), 0 quadratic interpolation polynomials. Based on this scheme, we introduce a time discretization scheme to approximate the time fractional partial differential equation and show by using Laplace transform methods that the time discretization scheme has the convergence rate O (k 3 - α), 0 0 for smooth and nonsmooth data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the theoretical results are consistent with the numerical results.

Modulating functions-based method for parameters and source estimation in one-dimensional partial differential equations

KAUST Repository

Asiri, Sharefa M.; Laleg-Kirati, Taous-Meriem

2016-01-01

In this paper, modulating functions-based method is proposed for estimating space–time-dependent unknowns in one-dimensional partial differential equations. The proposed method simplifies the problem into a system of algebraic equations linear

Partial differential equation techniques for analysing animal movement: A comparison of different methods.

Science.gov (United States)

Wang, Yi-Shan; Potts, Jonathan R

2017-03-07

Recent advances in animal tracking have allowed us to uncover the drivers of movement in unprecedented detail. This has enabled modellers to construct ever more realistic models of animal movement, which aid in uncovering detailed patterns of space use in animal populations. Partial differential equations (PDEs) provide a popular tool for mathematically analysing such models. However, their construction often relies on simplifying assumptions which may greatly affect the model outcomes. Here, we analyse the effect of various PDE approximations on the analysis of some simple movement models, including a biased random walk, central-place foraging processes and movement in heterogeneous landscapes. Perhaps the most commonly-used PDE method dates back to a seminal paper of Patlak from 1953. However, our results show that this can be a very poor approximation in even quite simple models. On the other hand, more recent methods, based on transport equation formalisms, can provide more accurate results, as long as the kernel describing the animal's movement is sufficiently smooth. When the movement kernel is not smooth, we show that both the older and newer methods can lead to quantitatively misleading results. Our detailed analysis will aid future researchers in the appropriate choice of PDE approximation for analysing models of animal movement. Copyright © 2017 Elsevier Ltd. All rights reserved.

Partially ordered models

NARCIS (Netherlands)

Fernandez, R.; Deveaux, V.

2010-01-01

We provide a formal definition and study the basic properties of partially ordered chains (POC). These systems were proposed to model textures in image processing and to represent independence relations between random variables in statistics (in the later case they are known as Bayesian networks).

A linearizing transformation for the Korteweg-de Vries equation; generalizations to higher-dimensional nonlinear partial differential equations

NARCIS (Netherlands)

Dorren, H.J.S.

1998-01-01

It is shown that the Korteweg–de Vries (KdV) equation can be transformed into an ordinary linear partial differential equation in the wave number domain. Explicit solutions of the KdV equation can be obtained by subsequently solving this linear differential equation and by applying a cascade of

Atmospheres of partially differentiated super-Earth exoplanets

Science.gov (United States)

Schaefer, Laura; Sasselov, Dimitar

2015-11-01

Terrestrial exoplanets have been discovered in a range of sizes, densities and orbital locations that defy our expectations based upon the Solar System. Planets discovered to date with radii less than ~1.5-1.6 Earth radii all seem to fall on an iso-density curve with the Earth [1]. However, mass and radius determinations, which depend on the known properties of the host star, are not accurate enough to distinguish between a fully differentiated three-layer planet (core, mantle, ocean/atmosphere) and an incompletely differentiated planet [2]. Full differentiation of a planet will depend upon the conditions at the time of accretion, including the abundance of short-lived radioisotopes, which will vary from system to system, as well as the number of giant impacts the planet experiences. Furthermore, separation of metal and silicates at the much larger pressures found inside super-Earths will depend on how the chemistry of these materials change at high pressures. There are therefore hints emerging that not all super-Earths will be fully differentiated. Incomplete differentiation will result in a more reduced mantle oxidation state and may have implications for the composition of an outgassed atmosphere. Here we will present the first results from a chemical equilibrium model of the composition of such an outgassed atmosphere and discuss the possibility of distinguishing between fully and incompletely differentiated planets through atmospheric observations.[1] Rogers, L. 2015. ApJ, 801, 41. [2] Zeng, L. & Sasselov, D. 2013. PASP, 125, 227.

Bäcklund transformation of fractional Riccati equation and its applications to nonlinear fractional partial differential equations

International Nuclear Information System (INIS)

Lu, Bin

2012-01-01

In this Letter, the fractional derivatives in the sense of modified Riemann–Liouville derivative and the Bäcklund transformation of fractional Riccati equation are employed for constructing the exact solutions of nonlinear fractional partial differential equations. The power of this manageable method is presented by applying it to several examples. This approach can also be applied to other nonlinear fractional differential equations. -- Highlights: ► Backlund transformation of fractional Riccati equation is presented. ► A new method for solving nonlinear fractional differential equations is proposed. ► Three important fractional differential equations are solved successfully. ► Some new exact solutions of the fractional differential equations are obtained.

Analytical approach to linear fractional partial differential equations arising in fluid mechanics

International Nuclear Information System (INIS)

Momani, Shaher; Odibat, Zaid

2006-01-01

In this Letter, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving linear fractional partial differential equations arising in fluid mechanics. The fractional derivatives are described in the Caputo sense. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these methods, the solution takes the form of a convergent series with easily computable components. The corresponding solutions of the integer order equations are found to follow as special cases of those of fractional order equations. Some numerical examples are presented to illustrate the efficiency and reliability of the two methods

Parallels between control PDE's (Partial Differential Equations) and systems of ODE's (Ordinary Differential Equations)

Science.gov (United States)

Hunt, L. R.; Villarreal, Ramiro

1987-01-01

System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differentail equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralleled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.

Optimized difference schemes for multidimensional hyperbolic partial differential equations

Directory of Open Access Journals (Sweden)

Adrian Sescu

2009-04-01

Full Text Available In numerical solutions to hyperbolic partial differential equations in multidimensions, in addition to dispersion and dissipation errors, there is a grid-related error (referred to as isotropy error or numerical anisotropy that affects the directional dependence of the wave propagation. Difference schemes are mostly analyzed and optimized in one dimension, wherein the anisotropy correction may not be effective enough. In this work, optimized multidimensional difference schemes with arbitrary order of accuracy are designed to have improved isotropy compared to conventional schemes. The derivation is performed based on Taylor series expansion and Fourier analysis. The schemes are restricted to equally-spaced Cartesian grids, so the generalized curvilinear transformation method and Cartesian grid methods are good candidates.

On the strong solution of a class of partial differential equations that arise in the pricing of mortgage backed securities

KAUST Repository

Parshad, Rana; Bayazit, Derviş; Barlow, Nathaniel S.; Prasad, V. Ramchandra

2011-01-01

We consider a reduced form pricing model for mortgage backed securities, formulated as a non-linear partial differential equation. We prove that the model possesses a weak solution. We then show that under additional regularity assumptions on the initial data, we also have a mild solution. This mild solution is shown to be a strong solution via further regularity arguments. We also numerically solve the reduced model via a Fourier spectral method. Lastly, we compare our numerical solution to real market data. We observe interestingly that the reduced model captures a number of recent market trends in this data, that have escaped previous models.

A dynamical regularization algorithm for solving inverse source problems of elliptic partial differential equations

Science.gov (United States)

Zhang, Ye; Gong, Rongfang; Cheng, Xiaoliang; Gulliksson, Mårten

2018-06-01

This study considers the inverse source problem for elliptic partial differential equations with both Dirichlet and Neumann boundary data. The unknown source term is to be determined by additional boundary conditions. Unlike the existing methods found in the literature, which usually employ the first-order in time gradient-like system (such as the steepest descent methods) for numerically solving the regularized optimization problem with a fixed regularization parameter, we propose a novel method with a second-order in time dissipative gradient-like system and a dynamical selected regularization parameter. A damped symplectic scheme is proposed for the numerical solution. Theoretical analysis is given for both the continuous model and the numerical algorithm. Several numerical examples are provided to show the robustness of the proposed algorithm.

2013 CIME Course Vector-valued Partial Differential Equations and Applications

CERN Document Server

Marcellini, Paolo

2017-01-01

Collating different aspects of Vector-valued Partial Differential Equations and Applications, this volume is based on the 2013 CIME Course with the same name which took place at Cetraro, Italy, under the scientific direction of John Ball and Paolo Marcellini. It contains the following contributions: The pullback equation (Bernard Dacorogna), The stability of the isoperimetric inequality (Nicola Fusco), Mathematical problems in thin elastic sheets: scaling limits, packing, crumpling and singularities (Stefan Müller), and Aspects of PDEs related to fluid flows (Vladimir Sverák). These lectures are addressed to graduate students and researchers in the field.

Study of coupled nonlinear partial differential equations for finding exact analytical solutions.

Science.gov (United States)

Khan, Kamruzzaman; Akbar, M Ali; Koppelaar, H

2015-07-01

Exact solutions of nonlinear partial differential equations (NPDEs) are obtained via the enhanced (G'/G)-expansion method. The method is subsequently applied to find exact solutions of the Drinfel'd-Sokolov-Wilson (DSW) equation and the (2+1)-dimensional Painlevé integrable Burgers (PIB) equation. The efficiency of this method for finding these exact solutions is demonstrated. The method is effective and applicable for many other NPDEs in mathematical physics.

Cellular automata for spatiotemporal pattern formation from reaction–diffusion partial differential equations

International Nuclear Information System (INIS)

Ohmori, Shousuke; Yamazaki, Yoshihiro

2016-01-01

Ultradiscrete equations are derived from a set of reaction–diffusion partial differential equations, and cellular automaton rules are obtained on the basis of the ultradiscrete equations. Some rules reproduce the dynamical properties of the original reaction–diffusion equations, namely, bistability and pulse annihilation. Furthermore, other rules bring about soliton-like preservation and periodic pulse generation with a pacemaker, which are not obtained from the original reaction–diffusion equations. (author)

Essential partial differential equations analytical and computational aspects

CERN Document Server

Griffiths, David F; Silvester, David J

2015-01-01

This volume provides an introduction to the analytical and numerical aspects of partial differential equations (PDEs). It unifies an analytical and computational approach for these; the qualitative behaviour of solutions being established using classical concepts: maximum principles and energy methods.   Notable inclusions are the treatment of irregularly shaped boundaries, polar coordinates and the use of flux-limiters when approximating hyperbolic conservation laws. The numerical analysis of difference schemes is rigorously developed using discrete maximum principles and discrete Fourier analysis. A novel feature is the inclusion of a chapter containing projects, intended for either individual or group study, that cover a range of topics such as parabolic smoothing, travelling waves, isospectral matrices, and the approximation of multidimensional advection–diffusion problems.   The underlying theory is illustrated by numerous examples and there are around 300 exercises, designed to promote and test unde...

Partially molten magma ocean model

International Nuclear Information System (INIS)

Shirley, D.N.

1983-01-01

The properties of the lunar crust and upper mantle can be explained if the outer 300-400 km of the moon was initially only partially molten rather than fully molten. The top of the partially molten region contained about 20% melt and decreased to 0% at 300-400 km depth. Nuclei of anorthositic crust formed over localized bodies of magma segregated from the partial melt, then grew peripherally until they coverd the moon. Throughout most of its growth period the anorthosite crust floated on a layer of magma a few km thick. The thickness of this layer is regulated by the opposing forces of loss of material by fractional crystallization and addition of magma from the partial melt below. Concentrations of Sr, Eu, and Sm in pristine ferroan anorthosites are found to be consistent with this model, as are trends for the ferroan anorthosites and Mg-rich suites on a diagram of An in plagioclase vs. mg in mafics. Clustering of Eu, Sr, and mg values found among pristine ferroan anorthosites are predicted by this model

Mathematical Modelling of Intraretinal Oxygen Partial Pressure ...

African Journals Online (AJOL)

Purpose: The aim of our present work is to develop a simple steady state model for intraretinal oxygen partial pressure distribution and to investigate the effect of various model parameters on the partial pressure distribution under adapted conditions of light and darkness.. Method: A simple eight-layered mathematical model ...

Simple equation method for nonlinear partial differential equations and its applications

Directory of Open Access Journals (Sweden)

Taher A. Nofal

2016-04-01

Full Text Available In this article, we focus on the exact solution of the some nonlinear partial differential equations (NLPDEs such as, Kodomtsev–Petviashvili (KP equation, the (2 + 1-dimensional breaking soliton equation and the modified generalized Vakhnenko equation by using the simple equation method. In the simple equation method the trial condition is the Bernoulli equation or the Riccati equation. It has been shown that the method provides a powerful mathematical tool for solving nonlinear wave equations in mathematical physics and engineering problems.

Student Solutions Manual to Boundary Value Problems and Partial Differential Equations

CERN Document Server

Powers, David L

2005-01-01

This student solutions manual accompanies the text, Boundary Value Problems and Partial Differential Equations, 5e. The SSM is available in print via PDF or electronically, and provides the student with the detailed solutions of the odd-numbered problems contained throughout the book.Provides students with exercises that skillfully illustrate the techniques used in the text to solve science and engineering problemsNearly 900 exercises ranging in difficulty from basic drills to advanced problem-solving exercisesMany exercises based on current engineering applications

Differential-discrete mathematical model of two phase flow heat exchanger

International Nuclear Information System (INIS)

Debeljkovic, D.Lj.; Zitek, Pavel; Simeunovic, G.; Inard, Christian

2007-01-01

A dynamic thermal-hydraulic mathematical model of evaporator dynamics of a once - through sub critical steam generator is derived and presented. This model allows the investigation of evaporator dynamics including its transients responses. The evaporator was considered as a part of three-section (economizer, evaporator and super-heater) model with time varying phase boundaries and is described by a set of linearized discrete - difference equations which, with some other algebraic equations, constitutes a closed system of equations possible for exact computer solution. This model has been derived upon the fundamental equations of mass, energy and momentum balance. For the first time, a discrete differential approach has been applied in order to investigate such complex, two phase processes. Namely, this approach allows one to escape from the model of this process usually described by a set of partial differential equations and enables one, using this method, to simulate evaporators dynamics in an extraordinarily simple way. In current literature this approach is sometimes called physical discretization. (author)

Numerical and computational analysis of the partial differential equations in hydrocodes and wavecodes

International Nuclear Information System (INIS)

Hicks, D.L.; Walsh, R.T.

1976-06-01

Discrete methods for the solution of the partial differential equations arising in hydrocodes and wavecodes are presented in a tutorial fashion. By discrete methods is meant, for example, the methods of finite differences, finite elements, discretized characteristics, etc. The concepts of stability, consistency, convergence, order of accuracy, true accuracy, etc., and their relevance to the hydrocodes and wavecodes are discussed

On the solution of elliptic partial differential equations on regions with corners

International Nuclear Information System (INIS)

Serkh, Kirill; Rokhlin, Vladimir

2016-01-01

In this paper we investigate the solution of boundary value problems on polygonal domains for elliptic partial differential equations. We observe that when the problems are formulated as the boundary integral equations of classical potential theory, the solutions are representable by series of elementary functions. In addition to being analytically perspicuous, the resulting expressions lend themselves to the construction of accurate and efficient numerical algorithms. The results are illustrated by a number of numerical examples.

Baecklund transformations and zero-curvature representations of systems of partial differential equations

International Nuclear Information System (INIS)

Brandt, F.

1993-01-01

It is shown that Baecklund transformations (BTs) and zero-curvature representations (ZCRs) of systems of partial differential equations (PDEs) are closely related. The connection is established by nonlinear representations of the symmetry group underlying the ZCR which induce gauge transformations relating different BTs. This connection is used to construct BTs from ZCRs (and vice versa). Furthermore a procedure is outlined which allows a systematic search for ZCRs of a given system of PDEs. (orig.)

Study of coupled nonlinear partial differential equations for finding exact analytical solutions

Science.gov (United States)

Khan, Kamruzzaman; Akbar, M. Ali; Koppelaar, H.

2015-01-01

Exact solutions of nonlinear partial differential equations (NPDEs) are obtained via the enhanced (G′/G)-expansion method. The method is subsequently applied to find exact solutions of the Drinfel'd–Sokolov–Wilson (DSW) equation and the (2+1)-dimensional Painlevé integrable Burgers (PIB) equation. The efficiency of this method for finding these exact solutions is demonstrated. The method is effective and applicable for many other NPDEs in mathematical physics. PMID:26587256

Uniqueness of global quasi-classical solutions of the Cauchy problems for first-order nonlinear partial differential equations

International Nuclear Information System (INIS)

Tran Duc Van

1994-01-01

The notion of global quasi-classical solutions of the Cauchy problems for first-order nonlinear partial differential equations is presented, some uniqueness theorems and a stability result are established by the method based on the theory of differential inclusions. In particular, the answer to an open problem of S.N. Kruzhkov is given. (author). 10 refs, 1 fig

Nonlinear partial differential equations for scientists and engineers

CERN Document Server

Debnath, Lokenath

1997-01-01

"An exceptionally complete overview. There are numerous examples and the emphasis is on applications to almost all areas of science and engineering. There is truly something for everyone here. This reviewer feels that it is a very hard act to follow, and recommends it strongly. [This book] is a jewel." ---Applied Mechanics Review (Review of First Edition) This expanded and revised second edition is a comprehensive and systematic treatment of linear and nonlinear partial differential equations and their varied applications. Building upon the successful material of the first book, this edition contains updated modern examples and applications from areas of fluid dynamics, gas dynamics, plasma physics, nonlinear dynamics, quantum mechanics, nonlinear optics, acoustics, and wave propagation. Methods and properties of solutions are presented, along with their physical significance, making the book more useful for a diverse readership. Topics and key features: * Thorough coverage of derivation and methods of soluti...

Reversible effects of oxygen partial pressure on genes associated with placental angiogenesis and differentiation in primary-term cytotrophoblast cell culture.

Science.gov (United States)

Debiève, F; Depoix, C; Gruson, D; Hubinont, C

2013-09-01

Timely regulated changes in oxygen partial pressure are important for placental formation. Disturbances could be responsible for pregnancy-related diseases like preeclampsia and intrauterine growth restriction. We aimed to (i) determine the effect of oxygen partial pressure on cytotrophoblast differentiation; (ii) measure mRNA expression and protein secretion from genes associated with placental angiogenesis; and (iii) determine the reversibility of these effects at different oxygen partial pressures. Term cytotrophoblasts were incubated at 21% and 2.5% O2 for 96 hr, or were switched between the two oxygen concentrations after 48 hr. Real-time PCR and enzyme-linked immunosorbent assays (ELISAs) were used to evaluate cell fusion and differentiation, measuring transcript levels for those genes involved in cell fusion and placental angiogenesis, including VEGF, PlGF, VEGFR1, sVEGFR1, sENG, INHA, and GCM1. Cytotrophoblasts underwent fusion and differentiation in 2.5% O2 . PlGF expression was inhibited while sVEGFR1 expression increased. VEGF and sENG mRNA expressions increased in 2.5% compared to 21% O2 , but no protein was detected in the cell supernatants. Finally, GCM1 mRNA expression increased during trophoblast differentiation at 21% O2 , but was inhibited at 2.5% O2 . These mRNA expression effects were reversed by returning the cells to 21% O2 . Thus, low-oxygen partial pressure does not inhibit term-cytotrophoblast cell fusion and differentiation in vitro. Lowering the oxygen partial pressure from 21% to 2.5% caused normal-term trophoblasts to reversibly modify their expression of genes associated with placental angiogenesis. This suggests that modifications observed in pregnancy diseases such as preeclampsia or growth retardation are probably due to an extrinsic effect on trophoblasts. Copyright © 2013 Wiley Periodicals, Inc.

Parametric Borel summability for some semilinear system of partial differential equations

Directory of Open Access Journals (Sweden)

Hiroshi Yamazawa

2015-01-01

Full Text Available In this paper we study the Borel summability of formal solutions with a parameter of first order semilinear system of partial differential equations with \\(n\\ independent variables. In [Singular perturbation of linear systems with a regular singularity, J. Dynam. Control. Syst. 8 (2002, 313-322], Balser and Kostov proved the Borel summability of formal solutions with respect to a singular perturbation parameter for a linear equation with one independent variable. We shall extend their results to a semilinear system of equations with general independent variables.

"Real-Time Optical Laboratory Linear Algebra Solution Of Partial Differential Equations"

Science.gov (United States)

Casasent, David; Jackson, James

1986-03-01

A Space Integrating (SI) Optical Linear Algebra Processor (OLAP) employing space and frequency-multiplexing, new partitioning and data flow, and achieving high accuracy performance with a non base-2 number system is described. Laboratory data on the performance of this system and the solution of parabolic Partial Differential Equations (PDEs) is provided. A multi-processor OLAP system is also described for the first time. It use in the solution of multiple banded matrices that frequently arise is then discussed. The utility and flexibility of this processor compared to digital systolic architectures should be apparent.

Lp Theory for Super-Parabolic Backward Stochastic Partial Differential Equations in the Whole Space

International Nuclear Information System (INIS)

Du Kai; Qiu, Jinniao; Tang Shanjian

2012-01-01

This paper is concerned with semi-linear backward stochastic partial differential equations (BSPDEs for short) of super-parabolic type. An L p -theory is given for the Cauchy problem of BSPDEs, separately for the case of p∈(1,2] and for the case of p∈(2,∞). A comparison theorem is also addressed.

The application of Legendre-tau approximation to parameter identification for delay and partial differential equations

Science.gov (United States)

Ito, K.

1983-01-01

Approximation schemes based on Legendre-tau approximation are developed for application to parameter identification problem for delay and partial differential equations. The tau method is based on representing the approximate solution as a truncated series of orthonormal functions. The characteristic feature of the Legendre-tau approach is that when the solution to a problem is infinitely differentiable, the rate of convergence is faster than any finite power of 1/N; higher accuracy is thus achieved, making the approach suitable for small N.

Conservation laws for certain time fractional nonlinear systems of partial differential equations

Science.gov (United States)

Singla, Komal; Gupta, R. K.

2017-12-01

In this study, an extension of the concept of nonlinear self-adjointness and Noether operators is proposed for calculating conserved vectors of the time fractional nonlinear systems of partial differential equations. In our recent work (J Math Phys 2016; 57: 101504), by proposing the symmetry approach for time fractional systems, the Lie symmetries for some fractional nonlinear systems have been derived. In this paper, the obtained infinitesimal generators are used to find conservation laws for the corresponding fractional systems.

Solution of Fractional Partial Differential Equations in Fluid Mechanics by Extension of Some Iterative Method

Directory of Open Access Journals (Sweden)

A. A. Hemeda

2013-01-01

Full Text Available An extension of the so-called new iterative method (NIM has been used to handle linear and nonlinear fractional partial differential equations. The main property of the method lies in its flexibility and ability to solve nonlinear equations accurately and conveniently. Therefore, a general framework of the NIM is presented for analytical treatment of fractional partial differential equations in fluid mechanics. The fractional derivatives are described in the Caputo sense. Numerical illustrations that include the fractional wave equation, fractional Burgers equation, fractional KdV equation, fractional Klein-Gordon equation, and fractional Boussinesq-like equation are investigated to show the pertinent features of the technique. Comparison of the results obtained by the NIM with those obtained by both Adomian decomposition method (ADM and the variational iteration method (VIM reveals that the NIM is very effective and convenient. The basic idea described in this paper is expected to be further employed to solve other similar linear and nonlinear problems in fractional calculus.

Parabolic partial differential equations with discrete state-dependent delay: Classical solutions and solution manifold

Czech Academy of Sciences Publication Activity Database

Krisztin, T.; Rezunenko, Oleksandr

2016-01-01

Roč. 260, č. 5 (2016), s. 4454-4472 ISSN 0022-0396 R&D Projects: GA ČR GAP103/12/2431 Institutional support: RVO:67985556 Keywords : Parabolic partial differential equations * State dependent delay * Solution manifold Subject RIV: BC - Control Systems Theory Impact factor: 1.988, year: 2016 http://library.utia.cas.cz/separaty/2016/AS/rezunenko-0457879.pdf

Rail-to-rail low-power fully differential OTA utilizing adaptive biasing and partial feedback

DEFF Research Database (Denmark)

Tuan Vu, Cao; Wisland, Dag T.; Lande, Tor Sverre

consumption. The DC-gain of the proposed OTA is improved by adding a partial feedback loop. A Common-Mode Feedback (CMFB) circuit is required for fully differential rail-to-rail operation. Simulations show that the OTA topology has a low stand-by power consumption of 96μW and a high FoM of 3.84 [(V...

New finite volume methods for approximating partial differential equations on arbitrary meshes

International Nuclear Information System (INIS)

Hermeline, F.

2008-12-01

This dissertation presents some new methods of finite volume type for approximating partial differential equations on arbitrary meshes. The main idea lies in solving twice the problem to be dealt with. One addresses the elliptic equations with variable (anisotropic, antisymmetric, discontinuous) coefficients, the parabolic linear or non linear equations (heat equation, radiative diffusion, magnetic diffusion with Hall effect), the wave type equations (Maxwell, acoustics), the elasticity and Stokes'equations. Numerous numerical experiments show the good behaviour of this type of method. (author)

XMDS2: Fast, scalable simulation of coupled stochastic partial differential equations

Science.gov (United States)

Dennis, Graham R.; Hope, Joseph J.; Johnsson, Mattias T.

2013-01-01

XMDS2 is a cross-platform, GPL-licensed, open source package for numerically integrating initial value problems that range from a single ordinary differential equation up to systems of coupled stochastic partial differential equations. The equations are described in a high-level XML-based script, and the package generates low-level optionally parallelised C++ code for the efficient solution of those equations. It combines the advantages of high-level simulations, namely fast and low-error development, with the speed, portability and scalability of hand-written code. XMDS2 is a complete redesign of the XMDS package, and features support for a much wider problem space while also producing faster code. Program summaryProgram title: XMDS2 Catalogue identifier: AENK_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AENK_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 2 No. of lines in distributed program, including test data, etc.: 872490 No. of bytes in distributed program, including test data, etc.: 45522370 Distribution format: tar.gz Programming language: Python and C++. Computer: Any computer with a Unix-like system, a C++ compiler and Python. Operating system: Any Unix-like system; developed under Mac OS X and GNU/Linux. RAM: Problem dependent (roughly 50 bytes per grid point) Classification: 4.3, 6.5. External routines: The external libraries required are problem-dependent. Uses FFTW3 Fourier transforms (used only for FFT-based spectral methods), dSFMT random number generation (used only for stochastic problems), MPI message-passing interface (used only for distributed problems), HDF5, GNU Scientific Library (used only for Bessel-based spectral methods) and a BLAS implementation (used only for non-FFT-based spectral methods). Nature of problem: General coupled initial-value stochastic partial differential equations. Solution method: Spectral method

Graphical Rasch models

DEFF Research Database (Denmark)

Kreiner, Svend; Christensen, Karl Bang

Rasch models; Partial Credit models; Rating Scale models; Item bias; Differential item functioning; Local independence; Graphical models......Rasch models; Partial Credit models; Rating Scale models; Item bias; Differential item functioning; Local independence; Graphical models...

Differential models in ecology

International Nuclear Information System (INIS)

Barco Gomez, Carlos; Barco Gomez, German

2002-01-01

The models mathematical writings with differential equations are used to describe the populational behavior through the time of the animal species. These models can be lineal or no lineal. The differential models for unique specie include the exponential pattern of Malthus and the logistical pattern of Verlhust. The lineal differential models to describe the interaction between two species include the competition relationships, predation and symbiosis

A Model for Partial Kantian Cooperation

OpenAIRE

Kordonis, Ioannis

2016-01-01

This work presents a game theoretic model to describe game situations in which there is a partial cooperation among the players. Specifically, we assume that the players partially follow Kant's "Categorical Imperative". The model is stated for games with a continuum of players and the basic assumption made is that the participants consider that they belong to virtual groups in which they optimize their actions as if they were bound to follow the same strategy. The relation with the Nash, (Ben...

The modified simplest equation method to look for exact solutions of nonlinear partial differential equations

OpenAIRE

Efimova, Olga Yu.

2010-01-01

The modification of simplest equation method to look for exact solutions of nonlinear partial differential equations is presented. Using this method we obtain exact solutions of generalized Korteweg-de Vries equation with cubic source and exact solutions of third-order Kudryashov-Sinelshchikov equation describing nonlinear waves in liquids with gas bubbles.

Partially Hidden Markov Models

DEFF Research Database (Denmark)

Forchhammer, Søren Otto; Rissanen, Jorma

1996-01-01

Partially Hidden Markov Models (PHMM) are introduced. They differ from the ordinary HMM's in that both the transition probabilities of the hidden states and the output probabilities are conditioned on past observations. As an illustration they are applied to black and white image compression where...

On k-summability of formal solutions for certain partial differential operators with polynomial coefficients

Directory of Open Access Journals (Sweden)

Kunio Ichinobe

2015-01-01

Full Text Available We study the \\(k\\-summability of divergent formal solutions for the Cauchy problem of certain linear partial differential operators with coefficients which are polynomial in \\(t\\. We employ the method of successive approximation in order to construct the formal solutions and to obtain the properties of analytic continuation of the solutions of convolution equations and their exponential growth estimates.

An approximation theory for nonlinear partial differential equations with applications to identification and control

Science.gov (United States)

Banks, H. T.; Kunisch, K.

1982-01-01

Approximation results from linear semigroup theory are used to develop a general framework for convergence of approximation schemes in parameter estimation and optimal control problems for nonlinear partial differential equations. These ideas are used to establish theoretical convergence results for parameter identification using modal (eigenfunction) approximation techniques. Results from numerical investigations of these schemes for both hyperbolic and parabolic systems are given.

Equivalent construction of the infinitesimal time translation operator in algebraic dynamics algorithm for partial differential evolution equation

Institute of Scientific and Technical Information of China (English)

无

2010-01-01

We give an equivalent construction of the infinitesimal time translation operator for partial differential evolution equation in the algebraic dynamics algorithm proposed by Shun-Jin Wang and his students. Our construction involves only simple partial differentials and avoids the derivative terms of δ function which appear in the course of computation by means of Wang-Zhang operator. We prove Wang’s equivalent theorem which says that our construction and Wang-Zhang’s are equivalent. We use our construction to deal with several typical equations such as nonlinear advection equation, Burgers equation, nonlinear Schrodinger equation, KdV equation and sine-Gordon equation, and obtain at least second order approximate solutions to them. These equations include the cases of real and complex field variables and the cases of the first and the second order time derivatives.

Reduced basis methods for partial differential equations an introduction

CERN Document Server

Quarteroni, Alfio; Negri, Federico

2016-01-01

This book provides a basic introduction to reduced basis (RB) methods for problems involving the repeated solution of partial differential equations (PDEs) arising from engineering and applied sciences, such as PDEs depending on several parameters and PDE-constrained optimization.  The book presents a general mathematical formulation of RB methods, analyzes their fundamental theoretical properties, discusses the related algorithmic and implementation aspects, and highlights their built-in algebraic and geometric structures.  More specifically, the authors discuss alternative strategies for constructing accurate RB spaces using greedy algorithms and proper orthogonal decomposition techniques, investigate their approximation properties and analyze offline-online decomposition strategies aimed at the reduction of computational complexity. Furthermore, they carry out both a priori and a posteriori error analysis.  The whole mathematical presentation is made more stimulating by the use of representative examp...

The solids-flux theory--confirmation and extension by using partial differential equations.

Science.gov (United States)

Diehl, Stefan

2008-12-01

The solids-flux theory has been used for half a century as a tool for estimating concentration and fluxes in the design and operation of secondary settling tanks during stationary conditions. The flux theory means that the conservation of mass is used in one dimension together with the batch-settling flux function according to the Kynch assumption. The flux theory results correspond to stationary solutions of a partial differential equation, a conservation law, with discontinuous coefficients modelling the continuous-sedimentation process in one dimension. The mathematical analysis of such an equation is intricate, partly since it cannot be interpreted in the classical sense. Recent results, however, make it possible to partly confirm and extend the previous flux theory statements, partly draw new conclusions also on the dynamic behaviour and the possibilities and limitations for control. We use here a single example of an ideal settling tank and a given batch-settling flux in a whole series of calculations. The mathematical results are adapted towards the application and many of them are conveniently presented in terms of operating charts.

Semi-groups of operators and some of their applications to partial differential equations

International Nuclear Information System (INIS)

Kisynski, J.

1976-01-01

Basic notions and theorems of the theory of one-parameter semi-groups of linear operators are given, illustrated by some examples concerned with linear partial differential operators. For brevity, some important and widely developed parts of the semi-group theory such as the general theory of holomorphic semi-groups or the theory of temporally inhomogeneous evolution equations are omitted. This omission includes also the very important application of semi-groups to investigating stochastic processes. (author)

Total, partial and differential ionization cross sections in proton-hydrogen collisions at low energy

Energy Technology Data Exchange (ETDEWEB)

Zou, Shiyang [Graduate University for Advanced Studies, School of Mathematical and Physical Science, Toki, Gifu (Japan); Pichl, Lukas [University of Aizu, Foundation of Computer Science Laboratory, Aizuwakamatsu, Fukushima (Japan); Kimura, Mineo [Yamaguchi Univ., Graduate School of Science and Engineering, Ube, Yamaguchi (Japan); Kato, Takako [National Inst. for Fusion Science, Toki, Gifu (Japan)

2003-01-01

Single-differential, partial and total ionization cross sections for the proton-hydrogen collision system at low energy range (0.1-10 keV/amu) are determined by using the electron translation factor corrected molecular-orbital close-coupling method. Full convergence of ionization cross sections as a function of H{sub 2}{sup +} molecular basis size is achieved by including up to 10 bound states, and 11 continuum partial waves. The present cross sections are in an excellent agreement with the recent experiments of Shah et al., but decrease more rapidly than the cross sections measured by Pieksma et al. with decreasing energy. The calculated cross section data are included in this report. (author)

PDASAC, Partial Differential Sensitivity Analysis of Stiff System

International Nuclear Information System (INIS)

Caracotsios, M.; Stewart, W.E.

2001-01-01

1 - Description of program or function: PDASAC solves stiff, nonlinear initial-boundary-value problems in a timelike dimension t and a space dimension x. Plane, circular cylindrical or spherical boundaries can be handled. Mixed-order systems of partial differential and algebraic equations can be analyzed with members of order or 0 or 1 in t, 0, 1 or 2 in x. Parametric sensitivities of the calculated states are computed simultaneously on request, via the Jacobian of the state equations. Initial and boundary conditions are efficiently reconciled. Local error control (in the max-norm or the 2-norm) is provided for the state vector and can include the parametric sensitivities if desired. 2 - Method of solution: The method of lines is used, with a user- selected x-grid and a minimum-bandwidth finite-difference approximations of the x-derivatives. Starting conditions are reconciled with a damped Newton algorithm adapted from Bain and Stewart (1991). Initial step selection is done by the first-order algorithms of Shampine (1987), extended here to differential- algebraic equation systems. The solution is continued with the DASSL predictor-corrector algorithm (Petzold 1983, Brenan et al. 1989) with the initial acceleration phase deleted and with row scaling of the Jacobian added. The predictor and corrector are expressed in divided-difference form, with the fixed-leading-coefficient form of corrector (Jackson and Sacks-Davis 1989; Brenan et al. 1989). Weights for the error tests are updated in each step with the user's tolerances at the predicted state. Sensitivity analysis is performed directly on the corrector equations of Caracotsios and Stewart (1985) and is extended here to the initialization when needed. 3 - Restrictions on the complexity of the problem: This algorithm, like DASSL, performs well on differential-algebraic equation systems of index 0 and 1 but not on higher-index systems; see Brenan et al. (1989). The user assigned the work array lengths and the output

A hybrid algorithm for coupling partial differential equation and compartment-based dynamics.

Science.gov (United States)

Harrison, Jonathan U; Yates, Christian A

2016-09-01

Stochastic simulation methods can be applied successfully to model exact spatio-temporally resolved reaction-diffusion systems. However, in many cases, these methods can quickly become extremely computationally intensive with increasing particle numbers. An alternative description of many of these systems can be derived in the diffusive limit as a deterministic, continuum system of partial differential equations (PDEs). Although the numerical solution of such PDEs is, in general, much more efficient than the full stochastic simulation, the deterministic continuum description is generally not valid when copy numbers are low and stochastic effects dominate. Therefore, to take advantage of the benefits of both of these types of models, each of which may be appropriate in different parts of a spatial domain, we have developed an algorithm that can be used to couple these two types of model together. This hybrid coupling algorithm uses an overlap region between the two modelling regimes. By coupling fluxes at one end of the interface and using a concentration-matching condition at the other end, we ensure that mass is appropriately transferred between PDE- and compartment-based regimes. Our methodology gives notable reductions in simulation time in comparison with using a fully stochastic model, while maintaining the important stochastic features of the system and providing detail in appropriate areas of the domain. We test our hybrid methodology robustly by applying it to several biologically motivated problems including diffusion and morphogen gradient formation. Our analysis shows that the resulting error is small, unbiased and does not grow over time. © 2016 The Authors.

Variance Function Partially Linear Single-Index Models1.

Science.gov (United States)

Lian, Heng; Liang, Hua; Carroll, Raymond J

2015-01-01

We consider heteroscedastic regression models where the mean function is a partially linear single index model and the variance function depends upon a generalized partially linear single index model. We do not insist that the variance function depend only upon the mean function, as happens in the classical generalized partially linear single index model. We develop efficient and practical estimation methods for the variance function and for the mean function. Asymptotic theory for the parametric and nonparametric parts of the model is developed. Simulations illustrate the results. An empirical example involving ozone levels is used to further illustrate the results, and is shown to be a case where the variance function does not depend upon the mean function.

Mixed problem with integral boundary condition for a high order mixed type partial differential equation

OpenAIRE

M. Denche; A. L. Marhoune

2003-01-01

In this paper, we study a mixed problem with integral boundary conditions for a high order partial differential equation of mixed type. We prove the existence and uniqueness of the solution. The proof is based on energy inequality, and on the density of the range of the operator generated by the considered problem.

Computation of Green function of the Schroedinger-like partial differential equations by the numerical functional integration

International Nuclear Information System (INIS)

Lobanov, Yu.Yu.; Shahbagian, R.R.; Zhidkov, E.P.

1991-01-01

A new method for numerical solution of the boundary problem for Schroedinger-like partial differential equations in R n is elaborated. The method is based on representation of multidimensional Green function in the form of multiple functional integral and on the use of approximation formulas which are constructed for such integrals. The convergence of approximations to the exact value is proved, the remainder of the formulas is estimated. Method reduces the initial differential problem to quadratures. 16 refs.; 7 tabs

Learning partial differential equations via data discovery and sparse optimization.

Science.gov (United States)

Schaeffer, Hayden

2017-01-01

We investigate the problem of learning an evolution equation directly from some given data. This work develops a learning algorithm to identify the terms in the underlying partial differential equations and to approximate the coefficients of the terms only using data. The algorithm uses sparse optimization in order to perform feature selection and parameter estimation. The features are data driven in the sense that they are constructed using nonlinear algebraic equations on the spatial derivatives of the data. Several numerical experiments show the proposed method's robustness to data noise and size, its ability to capture the true features of the data, and its capability of performing additional analytics. Examples include shock equations, pattern formation, fluid flow and turbulence, and oscillatory convection.

Solution of Nonlinear Partial Differential Equations by New Laplace Variational Iteration Method

Directory of Open Access Journals (Sweden)

Eman M. A. Hilal

2014-01-01

Full Text Available The aim of this study is to give a good strategy for solving some linear and nonlinear partial differential equations in engineering and physics fields, by combining Laplace transform and the modified variational iteration method. This method is based on the variational iteration method, Laplace transforms, and convolution integral, introducing an alternative Laplace correction functional and expressing the integral as a convolution. Some examples in physical engineering are provided to illustrate the simplicity and reliability of this method. The solutions of these examples are contingent only on the initial conditions.

Existence of pseudo almost periodic solutions for a class of partial functional differential equations

Directory of Open Access Journals (Sweden)

Hui-Sheng Ding

2013-04-01

Full Text Available In this paper, we first introduce a new class of pseudo almost periodic type functions and investigate some properties of pseudo almost periodic type functions; and then we discuss the existence of pseudo almost periodic solutions to the class of abstract partial functional differential equations $x'(t=Ax(t+f(t,x_t$ with finite delay in a Banach space X.

Simulation, optimal control and parametric sensitivity analysis of a molten carbonate fuel cell using a partial differential algebraic dynamic equation system; Simulation, Optimale Steuerung und Sensitivitaetsanalyse einer Schmelzkarbonat-Brennstoffzelle mithilfe eines partiellen differential-algebraischen dynamischen Gleichungssystems

Energy Technology Data Exchange (ETDEWEB)

Sternberg, K

2007-02-08

Molten carbonate fuel cells (MCFCs) allow an efficient and environmentally friendly energy production by converting the chemical energy contained in the fuel gas in virtue of electro-chemical reactions. In order to predict the effect of the electro-chemical reactions and to control the dynamical behavior of the fuel cell a mathematical model has to be found. The molten carbonate fuel cell (MCFC) can indeed be described by a highly complex,large scale, semi-linear system of partial differential algebraic equations. This system includes a reaction-diffusion-equation of parabolic type, several reaction-transport-equations of hyperbolic type, several ordinary differential equations and finally a system of integro-differential algebraic equations which describes the nonlinear non-standard boundary conditions for the entire partial differential algebraic equation system (PDAE-system). The existence of an analytical or the computability of a numerical solution for this high-dimensional PDAE-system depends on the kind of the differential equations and their special characteristics. Apart from theoretical investigations, the real process has to be controlled, more precisely optimally controlled. Hence, on the basis of the PDAE-system an optimal control problem is set up, whose analytical and numerical solvability is closely linked to the solvability of the PDAE-system. Moreover the solution of that optimal control problem is made more difficult by inaccuracies in the underlying database, which does not supply sufficiently accurate values for the model parameters. Therefore the optimal control problem must also be investigated with respect to small disturbances of model parameters. The aim of this work is to analyze the relevant dynamic behavior of MCFCs and to develop concepts for their optimal process control. Therefore this work is concerned with the simulation, the optimal control and the sensitivity analysis of a mathematical model for MCDCs, which can be characterized

Computationally efficient statistical differential equation modeling using homogenization

Science.gov (United States)

Hooten, Mevin B.; Garlick, Martha J.; Powell, James A.

2013-01-01

Statistical models using partial differential equations (PDEs) to describe dynamically evolving natural systems are appearing in the scientific literature with some regularity in recent years. Often such studies seek to characterize the dynamics of temporal or spatio-temporal phenomena such as invasive species, consumer-resource interactions, community evolution, and resource selection. Specifically, in the spatial setting, data are often available at varying spatial and temporal scales. Additionally, the necessary numerical integration of a PDE may be computationally infeasible over the spatial support of interest. We present an approach to impose computationally advantageous changes of support in statistical implementations of PDE models and demonstrate its utility through simulation using a form of PDE known as “ecological diffusion.” We also apply a statistical ecological diffusion model to a data set involving the spread of mountain pine beetle (Dendroctonus ponderosae) in Idaho, USA.

The orthogonal gradients method: A radial basis functions method for solving partial differential equations on arbitrary surfaces

KAUST Repository

Piret, Cé cile

2012-01-01

Much work has been done on reconstructing arbitrary surfaces using the radial basis function (RBF) method, but one can hardly find any work done on the use of RBFs to solve partial differential equations (PDEs) on arbitrary surfaces. In this paper

Constraints on the rheology of the partially molten mantle from numerical models of laboratory experiments

Science.gov (United States)

Rudge, J. F.; Alisic Jewell, L.; Rhebergen, S.; Katz, R. F.; Wells, G. N.

2015-12-01

One of the fundamental components in any dynamical model of melt transport is the rheology of partially molten rock. This rheology is poorly understood, and one way in which a better understanding can be obtained is by comparing the results of laboratory deformation experiments to numerical models. Here we present a comparison between numerical models and the laboratory setup of Qi et al. 2013 (EPSL), where a cylinder of partially molten rock containing rigid spherical inclusions was placed under torsion. We have replicated this setup in a finite element model which solves the partial differential equations describing the mechanical process of compaction. These computationally-demanding 3D simulations are only possible due to the recent development of a new preconditioning method for the equations of magma dynamics. The experiments show a distinct pattern of melt-rich and melt-depleted regions around the inclusions. In our numerical models, the pattern of melt varies with key rheological parameters, such as the ratio of bulk to shear viscosity, and the porosity- and strain-rate-dependence of the shear viscosity. These observed melt patterns therefore have the potential to constrain rheological properties. While there are many similarities between the experiments and the numerical models, there are also important differences, which highlight the need for better models of the physics of two-phase mantle/magma dynamics. In particular, the laboratory experiments display more pervasive melt-rich bands than is seen in our numerics.

The Reduction of Chazy Classes and Other Third-Order Differential Equations Related to Boundary Layer Flow Models

International Nuclear Information System (INIS)

Fakhar, K.; Kara, A. H.

2012-01-01

We study the symmetries, conservation laws and reduction of third-order equations that evolve from a prior reduction of models that arise in fluid phenomena. These could be the ordinary differential equations (ODEs) that are reductions of partial differential equations (PDEs) or, alternatively, PDEs related to given ODEs. In this class, the analysis includes the well-known Blasius, Chazy, and other associated third-order ODEs. (general)

Gröbner Bases and Generation of Difference Schemes for Partial Differential Equations

Directory of Open Access Journals (Sweden)

Vladimir P. Gerdt

2006-05-01

Full Text Available In this paper we present an algorithmic approach to the generation of fully conservative difference schemes for linear partial differential equations. The approach is based on enlargement of the equations in their integral conservation law form by extra integral relations between unknown functions and their derivatives, and on discretization of the obtained system. The structure of the discrete system depends on numerical approximation methods for the integrals occurring in the enlarged system. As a result of the discretization, a system of linear polynomial difference equations is derived for the unknown functions and their partial derivatives. A difference scheme is constructed by elimination of all the partial derivatives. The elimination can be achieved by selecting a proper elimination ranking and by computing a Gröbner basis of the linear difference ideal generated by the polynomials in the discrete system. For these purposes we use the difference form of Janet-like Gröbner bases and their implementation in Maple. As illustration of the described methods and algorithms, we construct a number of difference schemes for Burgers and Falkowich-Karman equations and discuss their numerical properties.

Effective action for stochastic partial differential equations

International Nuclear Information System (INIS)

Hochberg, David; Molina-Paris, Carmen; Perez-Mercader, Juan; Visser, Matt

1999-01-01

Stochastic partial differential equations (SPDEs) are the basic tool for modeling systems where noise is important. SPDEs are used for models of turbulence, pattern formation, and the structural development of the universe itself. It is reasonably well known that certain SPDEs can be manipulated to be equivalent to (nonquantum) field theories that nevertheless exhibit deep and important relationships with quantum field theory. In this paper we systematically extend these ideas: We set up a functional integral formalism and demonstrate how to extract all the one-loop physics for an arbitrary SPDE subject to arbitrary Gaussian noise. It is extremely important to realize that Gaussian noise does not imply that the field variables undergo Gaussian fluctuations, and that these nonquantum field theories are fully interacting. The limitation to one loop is not as serious as might be supposed: Experience with quantum field theories (QFTs) has taught us that one-loop physics is often quite adequate to give a good description of the salient issues. The limitation to one loop does, however, offer marked technical advantages: Because at one loop almost any field theory can be rendered finite using zeta function technology, we can sidestep the complications inherent in the Martin-Siggia-Rose formalism (the SPDE analog of the Becchi-Rouet-Stora-Tyutin formalism used in QFT) and instead focus attention on a minimalist approach that uses only the physical fields (this ''direct approach'' is the SPDE analog of canonical quantization using physical fields). After setting up the general formalism for the characteristic functional (partition function), we show how to define the effective action to all loops, and then focus on the one-loop effective action and its specialization to constant fields: the effective potential. The physical interpretation of the effective action and effective potential for SPDEs is addressed and we show that key features carry over from QFT to the case of

Effective action for stochastic partial differential equations

Energy Technology Data Exchange (ETDEWEB)

Hochberg, David [Laboratorio de Astrofisica Espacial y Fisica Fundamental, Apartado 50727, 28080 Madrid, (Spain); Centro de Astrobiologia, INTA, Carratera Ajalvir, Km. 4, 28850 Torrejon, Madrid, (Spain); Molina-Paris, Carmen [Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 (United States); Perez-Mercader, Juan [Laboratorio de Astrofisica Espacial y Fisica Fundamental, Apartado 50727, 28080 Madrid, (Spain); Visser, Matt [Physics Department, Washington University, Saint Louis, Missouri 63130-4899 (United States)

1999-12-01

Stochastic partial differential equations (SPDEs) are the basic tool for modeling systems where noise is important. SPDEs are used for models of turbulence, pattern formation, and the structural development of the universe itself. It is reasonably well known that certain SPDEs can be manipulated to be equivalent to (nonquantum) field theories that nevertheless exhibit deep and important relationships with quantum field theory. In this paper we systematically extend these ideas: We set up a functional integral formalism and demonstrate how to extract all the one-loop physics for an arbitrary SPDE subject to arbitrary Gaussian noise. It is extremely important to realize that Gaussian noise does not imply that the field variables undergo Gaussian fluctuations, and that these nonquantum field theories are fully interacting. The limitation to one loop is not as serious as might be supposed: Experience with quantum field theories (QFTs) has taught us that one-loop physics is often quite adequate to give a good description of the salient issues. The limitation to one loop does, however, offer marked technical advantages: Because at one loop almost any field theory can be rendered finite using zeta function technology, we can sidestep the complications inherent in the Martin-Siggia-Rose formalism (the SPDE analog of the Becchi-Rouet-Stora-Tyutin formalism used in QFT) and instead focus attention on a minimalist approach that uses only the physical fields (this ''direct approach'' is the SPDE analog of canonical quantization using physical fields). After setting up the general formalism for the characteristic functional (partition function), we show how to define the effective action to all loops, and then focus on the one-loop effective action and its specialization to constant fields: the effective potential. The physical interpretation of the effective action and effective potential for SPDEs is addressed and we show that key features carry over from

NATO Advanced Research Workshop on Approximation by Solutions of Partial Differential Equations, Quadrature Formulae, and Related Topics

CERN Document Server

Goldstein, M; Haussmann, W; Hayman, W; Rogge, L

1992-01-01

This volume consists of the proceedings of the NATO Advanced Research Workshop on Approximation by Solutions of Partial Differential Equations, Quadrature Formulae, and Related Topics, which was held at Hanstholm, Denmark. These proceedings include the main invited talks and contributed papers given during the workshop. The aim of these lectures was to present a selection of results of the latest research in the field. In addition to covering topics in approximation by solutions of partial differential equations and quadrature formulae, this volume is also concerned with related areas, such as Gaussian quadratures, the Pompelu problem, rational approximation to the Fresnel integral, boundary correspondence of univalent harmonic mappings, the application of the Hilbert transform in two dimensional aerodynamics, finely open sets in the limit set of a finitely generated Kleinian group, scattering theory, harmonic and maximal measures for rational functions and the solution of the classical Dirichlet problem. In ...

Implicit particle filtering for models with partial noise, and an application to geomagnetic data assimilation

Directory of Open Access Journals (Sweden)

M. Morzfeld

2012-06-01

Full Text Available Implicit particle filtering is a sequential Monte Carlo method for data assimilation, designed to keep the number of particles manageable by focussing attention on regions of large probability. These regions are found by minimizing, for each particle, a scalar function F of the state variables. Some previous implementations of the implicit filter rely on finding the Hessians of these functions. The calculation of the Hessians can be cumbersome if the state dimension is large or if the underlying physics are such that derivatives of F are difficult to calculate, as happens in many geophysical applications, in particular in models with partial noise, i.e. with a singular state covariance matrix. Examples of models with partial noise include models where uncertain dynamic equations are supplemented by conservation laws with zero uncertainty, or with higher order (in time stochastic partial differential equations (PDE or with PDEs driven by spatially smooth noise processes. We make the implicit particle filter applicable to such situations by combining gradient descent minimization with random maps and show that the filter is efficient, accurate and reliable because it operates in a subspace of the state space. As an example, we consider a system of nonlinear stochastic PDEs that is of importance in geomagnetic data assimilation.

An ansatz for solving nonlinear partial differential equations in mathematical physics.

Science.gov (United States)

Akbar, M Ali; Ali, Norhashidah Hj Mohd

2016-01-01

In this article, we introduce an ansatz involving exact traveling wave solutions to nonlinear partial differential equations. To obtain wave solutions using direct method, the choice of an appropriate ansatz is of great importance. We apply this ansatz to examine new and further general traveling wave solutions to the (1+1)-dimensional modified Benjamin-Bona-Mahony equation. Abundant traveling wave solutions are derived including solitons, singular solitons, periodic solutions and general solitary wave solutions. The solutions emphasize the nobility of this ansatz in providing distinct solutions to various tangible phenomena in nonlinear science and engineering. The ansatz could be more efficient tool to deal with higher dimensional nonlinear evolution equations which frequently arise in many real world physical problems.

Local Fractional Laplace Variational Iteration Method for Solving Linear Partial Differential Equations with Local Fractional Derivative

Directory of Open Access Journals (Sweden)

Ai-Min Yang

2014-01-01

Full Text Available The local fractional Laplace variational iteration method was applied to solve the linear local fractional partial differential equations. The local fractional Laplace variational iteration method is coupled by the local fractional variational iteration method and Laplace transform. The nondifferentiable approximate solutions are obtained and their graphs are also shown.

Mixed problem with nonlocal boundary conditions for a third-order partial differential equation of mixed type

OpenAIRE

Denche, M.; Marhoune, A. L.

2001-01-01

We study a mixed problem with integral boundary conditions for a third-order partial differential equation of mixed type. We prove the existence and uniqueness of the solution. The proof is based on two-sided a priori estimates and on the density of the range of the operator generated by the considered problem.

A Posteriori Finite Element Bounds for Sensitivity Derivatives of Partial-Differential-Equation Outputs. Revised

Science.gov (United States)

Lewis, Robert Michael; Patera, Anthony T.; Peraire, Jaume

1998-01-01

We present a Neumann-subproblem a posteriori finite element procedure for the efficient and accurate calculation of rigorous, 'constant-free' upper and lower bounds for sensitivity derivatives of functionals of the solutions of partial differential equations. The design motivation for sensitivity derivative error control is discussed; the a posteriori finite element procedure is described; the asymptotic bounding properties and computational complexity of the method are summarized; and illustrative numerical results are presented.

On mixed derivatives type high dimensional multi-term fractional partial differential equations approximate solutions

Science.gov (United States)

Talib, Imran; Belgacem, Fethi Bin Muhammad; Asif, Naseer Ahmad; Khalil, Hammad

2017-01-01

In this research article, we derive and analyze an efficient spectral method based on the operational matrices of three dimensional orthogonal Jacobi polynomials to solve numerically the mixed partial derivatives type multi-terms high dimensions generalized class of fractional order partial differential equations. We transform the considered fractional order problem to an easily solvable algebraic equations with the aid of the operational matrices. Being easily solvable, the associated algebraic system leads to finding the solution of the problem. Some test problems are considered to confirm the accuracy and validity of the proposed numerical method. The convergence of the method is ensured by comparing our Matlab software simulations based obtained results with the exact solutions in the literature, yielding negligible errors. Moreover, comparative results discussed in the literature are extended and improved in this study.

Advances in phase space analysis of partial differential equations in honor of Ferruccio Colombini's 60th birthday

CERN Document Server

Bove, Antonio; Murthy, MK Venkatesha

2009-01-01

This collection of original articles and surveys addresses the recent advances in linear and nonlinear aspects of the theory of partial differential equations. The key topics include operators as "sums of squares" of real and complex vector fields, nonlinear evolution equations, local solvability, and hyperbolic questions.

Towards information-optimal simulation of partial differential equations.

Science.gov (United States)

Leike, Reimar H; Enßlin, Torsten A

2018-03-01

Most simulation schemes for partial differential equations (PDEs) focus on minimizing a simple error norm of a discretized version of a field. This paper takes a fundamentally different approach; the discretized field is interpreted as data providing information about a real physical field that is unknown. This information is sought to be conserved by the scheme as the field evolves in time. Such an information theoretic approach to simulation was pursued before by information field dynamics (IFD). In this paper we work out the theory of IFD for nonlinear PDEs in a noiseless Gaussian approximation. The result is an action that can be minimized to obtain an information-optimal simulation scheme. It can be brought into a closed form using field operators to calculate the appearing Gaussian integrals. The resulting simulation schemes are tested numerically in two instances for the Burgers equation. Their accuracy surpasses finite-difference schemes on the same resolution. The IFD scheme, however, has to be correctly informed on the subgrid correlation structure. In certain limiting cases we recover well-known simulation schemes like spectral Fourier-Galerkin methods. We discuss implications of the approximations made.

A New Fractional Projective Riccati Equation Method for Solving Fractional Partial Differential Equations

International Nuclear Information System (INIS)

Feng Qing-Hua

2014-01-01

In this paper, a new fractional projective Riccati equation method is proposed to establish exact solutions for fractional partial differential equations in the sense of modified Riemann—Liouville derivative. This method can be seen as the fractional version of the known projective Riccati equation method. For illustrating the validity of this method, we apply this method to solve the space-time fractional Whitham—Broer—Kaup (WBK) equations and the nonlinear fractional Sharma—Tasso—Olever (STO) equation, and as a result, some new exact solutions for them are obtained. (general)

Automating Embedded Analysis Capabilities and Managing Software Complexity in Multiphysics Simulation, Part II: Application to Partial Differential Equations

Directory of Open Access Journals (Sweden)

Roger P. Pawlowski

2012-01-01

Full Text Available A template-based generic programming approach was presented in Part I of this series of papers [Sci. Program. 20 (2012, 197–219] that separates the development effort of programming a physical model from that of computing additional quantities, such as derivatives, needed for embedded analysis algorithms. In this paper, we describe the implementation details for using the template-based generic programming approach for simulation and analysis of partial differential equations (PDEs. We detail several of the hurdles that we have encountered, and some of the software infrastructure developed to overcome them. We end with a demonstration where we present shape optimization and uncertainty quantification results for a 3D PDE application.

Adaptive Partially Hidden Markov Models

DEFF Research Database (Denmark)

Forchhammer, Søren Otto; Rasmussen, Tage

1996-01-01

Partially Hidden Markov Models (PHMM) have recently been introduced. The transition and emission probabilities are conditioned on the past. In this report, the PHMM is extended with a multiple token version. The different versions of the PHMM are applied to bi-level image coding....

MRI differential diagnosis of complete and partial tears of the anterior cruciate ligament of the knee: the usefulness of oblique coronal T2-weighted image

International Nuclear Information System (INIS)

Lee, Seo Young; Shim, Jae Chan; Lee, Ghi Jai; Bang, Sun Woo; Ryu, Seok Jong; Kim, Ho Kyun; Kim, Jeong Seok

2002-01-01

To assess the usefulness of T2-weighted oblique coronal MR imaging (T2OCI) in the differential diagnosis of complete and partial tears of the anterior cruciate ligament (ACL) of the knee. Thirty-three patients with ACL tear (16 complete and 17 partial tears), comfirmed by arthroscopy, were included in this study. Conventional MR imaging and T2OCI were performed, and the findings were retrospectively reviewed by two radiologists in terms of continuity, shape, axis and internal signal intensity of the ligament. Each finding was tested if there were stastistically significant differences in its prevalence between partial and complete tears. The diagnostic accuracy of T2OCI and conventional MR imaging in the detection of partial and complete tears of the ACL were compared. Conventional MR imaging revealed no statistically significant finding for differential diagnosis of complete and partial ACL tears. The reliable and statistically significant (p<0.001) findings of T2OCI were complete discontinuity of the ligament in cases involving complete ACL tears (14 of 16 complete tears and 2 of 17 partial tears) and the preservation of the band form for partial ACL tears (2 of 16 complete tears and 15 of 17 partial tears). The accuracy of T2OCI and conventional MR imaging was 88% and 70%, respectively. When ACL injury is vague on conventional MR images, a modality which is more useful in the differential diagnosis of partial and complete tears of the ACL, and in predicting the site of a tear, is T2-weighted oblique coronal imaging

Survey of the status of finite element methods for partial differential equations

Science.gov (United States)

Temam, Roger

1986-01-01

The finite element methods (FEM) have proved to be a powerful technique for the solution of boundary value problems associated with partial differential equations of either elliptic, parabolic, or hyperbolic type. They also have a good potential for utilization on parallel computers particularly in relation to the concept of domain decomposition. This report is intended as an introduction to the FEM for the nonspecialist. It contains a survey which is totally nonexhaustive, and it also contains as an illustration, a report on some new results concerning two specific applications, namely a free boundary fluid-structure interaction problem and the Euler equations for inviscid flows.

Series expansion solutions for the multi-term time and space fractional partial differential equations in two- and three-dimensions

Science.gov (United States)

Ye, H.; Liu, F.; Turner, I.; Anh, V.; Burrage, K.

2013-09-01

Fractional partial differential equations with more than one fractional derivative in time describe some important physical phenomena, such as the telegraph equation, the power law wave equation, or the Szabo wave equation. In this paper, we consider two- and three-dimensional multi-term time and space fractional partial differential equations. The multi-term time-fractional derivative is defined in the Caputo sense, whose order belongs to the interval (1,2],(2,3],(3,4] or (0, m], and the space-fractional derivative is referred to as the fractional Laplacian form. We derive series expansion solutions based on a spectral representation of the Laplacian operator on a bounded region. Some applications are given for the two- and three-dimensional telegraph equation, power law wave equation and Szabo wave equation.

Mesh refinement and numerical sensitivity analysis for parameter calibration of partial differential equations

Science.gov (United States)

Becker, Roland; Vexler, Boris

2005-06-01

We consider the calibration of parameters in physical models described by partial differential equations. This task is formulated as a constrained optimization problem with a cost functional of least squares type using information obtained from measurements. An important issue in the numerical solution of this type of problem is the control of the errors introduced, first, by discretization of the equations describing the physical model, and second, by measurement errors or other perturbations. Our strategy is as follows: we suppose that the user defines an interest functional I, which might depend on both the state variable and the parameters and which represents the goal of the computation. First, we propose an a posteriori error estimator which measures the error with respect to this functional. This error estimator is used in an adaptive algorithm to construct economic meshes by local mesh refinement. The proposed estimator requires the solution of an auxiliary linear equation. Second, we address the question of sensitivity. Applying similar techniques as before, we derive quantities which describe the influence of small changes in the measurements on the value of the interest functional. These numbers, which we call relative condition numbers, give additional information on the problem under consideration. They can be computed by means of the solution of the auxiliary problem determined before. Finally, we demonstrate our approach at hand of a parameter calibration problem for a model flow problem.

The Adomian decomposition method for solving partial differential equations of fractal order in finite domains

Energy Technology Data Exchange (ETDEWEB)

El-Sayed, A.M.A. [Faculty of Science University of Alexandria (Egypt)]. E-mail: amasyed@hotmail.com; Gaber, M. [Faculty of Education Al-Arish, Suez Canal University (Egypt)]. E-mail: mghf408@hotmail.com

2006-11-20

The Adomian decomposition method has been successively used to find the explicit and numerical solutions of the time fractional partial differential equations. A different examples of special interest with fractional time and space derivatives of order {alpha}, 0<{alpha}=<1 are considered and solved by means of Adomian decomposition method. The behaviour of Adomian solutions and the effects of different values of {alpha} are shown graphically for some examples.

Partial Orders and Fully Abstract Models for Concurrency

DEFF Research Database (Denmark)

Engberg, Uffe Henrik

1990-01-01

In this thesis sets of labelled partial orders are employed as fundamental mathematical entities for modelling nondeterministic and concurrent processes thereby obtaining so-called noninterleaving semantics. Based on different closures of sets of labelled partial orders, simple algebraic language...

Explicit finite difference predictor and convex corrector with applications to hyperbolic partial differential equations

Science.gov (United States)

Dey, C.; Dey, S. K.

1983-01-01

An explicit finite difference scheme consisting of a predictor and a corrector has been developed and applied to solve some hyperbolic partial differential equations (PDEs). The corrector is a convex-type function which is applied at each time level and at each mesh point. It consists of a parameter which may be estimated such that for larger time steps the algorithm should remain stable and generate a fast speed of convergence to the steady-state solution. Some examples have been given.

Distribution of the Discretization and Algebraic Error in Numerical Solution of Partial Differential Equations

Czech Academy of Sciences Publication Activity Database

Papež, Jan; Liesen, J.; Strakoš, Z.

2014-01-01

Roč. 449, 15 May (2014), s. 89-114 ISSN 0024-3795 R&D Projects: GA AV ČR IAA100300802; GA ČR GA201/09/0917 Grant - others:GA MŠk(CZ) LL1202; GA UK(CZ) 695612 Institutional support: RVO:67985807 Keywords : numerical solution of partial differential equations * finite element method * adaptivity * a posteriori error analysis * discretization error * algebra ic error * spatial distribution of the error Subject RIV: BA - General Mathematics Impact factor: 0.939, year: 2014

A partially ionized plasma modeling; Un modele de plasma partiellement ionise

Energy Technology Data Exchange (ETDEWEB)

Le Thanh, K.C.; Raviart, P.A

2003-07-01

We propose a model for the partially ionized plasma sheaths near the anode of an anodic spot electric arc where the cathode is considered as an electron emitter. A fluid description takes into account the heating and the ionization of the plasma induced by the electron beam. As physical hypothesis we assume that the condition of charge neutrality is valid. According that the electron mass can be neglected compared to the ion mass, we can assume that ions and atoms have the same velocity and the same temperature. Electrons and heavy particles are then regarded as two separate fluids coexisting in the plasma. Governing equations are then multi-fluid equations with relaxation correction to the local thermodynamic equilibrium (LTE) and heating by Joule effect. Equations are solved by an operator splitting procedure. That is we first discretize the homogeneous conservation laws (i.e. without source terms) by a finite volume method. The second step is to solve the ordinary differential system (i.e, governing equation without transport terms) with an implicit scheme. (authors)

Grid generation for the solution of partial differential equations

Science.gov (United States)

Eiseman, Peter R.; Erlebacher, Gordon

1989-01-01

A general survey of grid generators is presented with a concern for understanding why grids are necessary, how they are applied, and how they are generated. After an examination of the need for meshes, the overall applications setting is established with a categorization of the various connectivity patterns. This is split between structured grids and unstructured meshes. Altogether, the categorization establishes the foundation upon which grid generation techniques are developed. The two primary categories are algebraic techniques and partial differential equation techniques. These are each split into basic parts, and accordingly are individually examined in some detail. In the process, the interrelations between the various parts are accented. From the established background in the primary techniques, consideration is shifted to the topic of interactive grid generation and then to adaptive meshes. The setting for adaptivity is established with a suitable means to monitor severe solution behavior. Adaptive grids are considered first and are followed by adaptive triangular meshes. Then the consideration shifts to the temporal coupling between grid generators and PDE-solvers. To conclude, a reflection upon the discussion, herein, is given.

Testing the generalized partial credit model

NARCIS (Netherlands)

Glas, Cornelis A.W.

1996-01-01

The partial credit model (PCM) (G.N. Masters, 1982) can be viewed as a generalization of the Rasch model for dichotomous items to the case of polytomous items. In many cases, the PCM is too restrictive to fit the data. Several generalizations of the PCM have been proposed. In this paper, a

Partial differential equation-based localization of a monopole source from a circular array.

Science.gov (United States)

Ando, Shigeru; Nara, Takaaki; Levy, Tsukassa

2013-10-01

Wave source localization from a sensor array has long been the most active research topics in both theory and application. In this paper, an explicit and time-domain inversion method for the direction and distance of a monopole source from a circular array is proposed. The approach is based on a mathematical technique, the weighted integral method, for signal/source parameter estimation. It begins with an exact form of the source-constraint partial differential equation that describes the unilateral propagation of wide-band waves from a single source, and leads to exact algebraic equations that include circular Fourier coefficients (phase mode measurements) as their coefficients. From them, nearly closed-form, single-shot and multishot algorithms are obtained that is suitable for use with band-pass/differential filter banks. Numerical evaluation and several experimental results obtained using a 16-element circular microphone array are presented to verify the validity of the proposed method.

Differentiation of mucosal disease from partial development of the paranasal sinuses in pediatric patients

International Nuclear Information System (INIS)

Duerinckx, A.J.; Whyte, A.M.; Lufkin, R.B.; Hall, T.R.; Kangarloo, H.

1988-01-01

On magnetic resonance (MR) images of pediatric patients, sinus mucosal disease may have an appearance similar to that of the normal partially developed sinus, leading to an increase in the number of patients labeled as having incidental sinusitis. The paranasal sinuses were retrospectively evaluated in 27 infants and children aged 0-11 years undergoing brain MR imaging for indications both unrelated and related to sinus disease. The authors developed criteria for grading paranasal sinus development and mucosal disease. Incidental mucosal disease is not uncommon, occurring in 28% of patients aged 0-7 years. In children under 3 years of age, inflammatory mucosal thickening and marrow surrounding the partially developed sinus have a high signal on many MR sequences and may be confused. Recognition of the low-intensity peripheral cortical margin of the sinus and awareness of the stages of normal sinus development allow differentiation

On a higher order multi-term time-fractional partial differential equation involving Caputo-Fabrizio derivative

OpenAIRE

Pirnapasov, Sardor; Karimov, Erkinjon

2017-01-01

In the present work we discuss higher order multi-term partial differential equation (PDE) with the Caputo-Fabrizio fractional derivative in time. We investigate a boundary value problem for fractional heat equation involving higher order Caputo-Fabrizio derivatives in time-variable. Using method of separation of variables and integration by parts, we reduce fractional order PDE to the integer order. We represent explicit solution of formulated problem in particular case by Fourier series.

Building bridges connections and challenges in modern approaches to numerical partial differential equations

CERN Document Server

Brezzi, Franco; Cangiani, Andrea; Georgoulis, Emmanuil

2016-01-01

This volume contains contributed survey papers from the main speakers at the LMS/EPSRC Symposium “Building bridges: connections and challenges in modern approaches to numerical partial differential equations”. This meeting took place in July 8-16, 2014, and its main purpose was to gather specialists in emerging areas of numerical PDEs, and explore the connections between the different approaches. The type of contributions ranges from the theoretical foundations of these new techniques, to the applications of them, to new general frameworks and unified approaches that can cover one, or more than one, of these emerging techniques.

Stable multi-domain spectral penalty methods for fractional partial differential equations

Science.gov (United States)

Xu, Qinwu; Hesthaven, Jan S.

2014-01-01

We propose stable multi-domain spectral penalty methods suitable for solving fractional partial differential equations with fractional derivatives of any order. First, a high order discretization is proposed to approximate fractional derivatives of any order on any given grids based on orthogonal polynomials. The approximation order is analyzed and verified through numerical examples. Based on the discrete fractional derivative, we introduce stable multi-domain spectral penalty methods for solving fractional advection and diffusion equations. The equations are discretized in each sub-domain separately and the global schemes are obtained by weakly imposed boundary and interface conditions through a penalty term. Stability of the schemes are analyzed and numerical examples based on both uniform and nonuniform grids are considered to highlight the flexibility and high accuracy of the proposed schemes.

Application of Generalized Hukuhara derivative approach in an economic production quantity model with partial trade credit policy under fuzzy environment

Directory of Open Access Journals (Sweden)

Pinki Majumder

2016-01-01

Full Text Available In this present study, a production inventory model with partial trade credit is formulated and solved in fuzzy environment via Generalized Hukuhara derivative approach. To capture the market, a supplier offers a trade credit period to its retailers. Due to this facility, retailer also offers a partial trade credit period to his/her customer to boost the demand of the item. In practical life situation, demands are generally dependent upon time. Constant demand of an item varies time to time. In this vague situation, demands are taken as time dependent, where its constant part is taken as Left Right - type fuzzy number. In this paper, Generalized Hukuhara derivative approach is used to solve the fuzzy inventory model. Four different cases are considered by using Generalized Hukuhara-(i differentiability and Generalized Hukuhara-(ii differentiability. The objective of this paper is to find out the optimal time so as the total inventory cost is minimum. Finally the model is solved by generalized reduced gradient method. The proposed model and technique are illustrated by numerical examples. Some sensitivity analyses both in tabular and graphical forms are presented and the effects of minimum cost with respect to various inventory parameters are discussed.

Maximum profile likelihood estimation of differential equation parameters through model based smoothing state estimates.

Science.gov (United States)

Campbell, D A; Chkrebtii, O

2013-12-01

Statistical inference for biochemical models often faces a variety of characteristic challenges. In this paper we examine state and parameter estimation for the JAK-STAT intracellular signalling mechanism, which exemplifies the implementation intricacies common in many biochemical inference problems. We introduce an extension to the Generalized Smoothing approach for estimating delay differential equation models, addressing selection of complexity parameters, choice of the basis system, and appropriate optimization strategies. Motivated by the JAK-STAT system, we further extend the generalized smoothing approach to consider a nonlinear observation process with additional unknown parameters, and highlight how the approach handles unobserved states and unevenly spaced observations. The methodology developed is generally applicable to problems of estimation for differential equation models with delays, unobserved states, nonlinear observation processes, and partially observed histories. Crown Copyright © 2013. Published by Elsevier Inc. All rights reserved.

Matrix-oriented implementation for the numerical solution of the partial differential equations governing flows and transport in porous media

KAUST Repository

Sun, Shuyu; Salama, Amgad; El-Amin, Mohamed

2012-01-01

In this paper we introduce a new technique for the numerical solution of the various partial differential equations governing flow and transport phenomena in porous media. This method is proposed to be used in high level programming languages like

Modeling Ability Differentiation in the Second-Order Factor Model

Science.gov (United States)

Molenaar, Dylan; Dolan, Conor V.; van der Maas, Han L. J.

2011-01-01

In this article we present factor models to test for ability differentiation. Ability differentiation predicts that the size of IQ subtest correlations decreases as a function of the general intelligence factor. In the Schmid-Leiman decomposition of the second-order factor model, we model differentiation by introducing heteroscedastic residuals,…

A Generalized Partial Credit Model: Application of an EM Algorithm.

Science.gov (United States)

Muraki, Eiji

1992-01-01

The partial credit model with a varying slope parameter is developed and called the generalized partial credit model (GPCM). Analysis results for simulated data by this and other polytomous item-response models demonstrate that the rating formulation of the GPCM is adaptable to the analysis of polytomous item responses. (SLD)

Testing a model of codependency for college students in Taiwan based on Bowen's concept of differentiation.

Science.gov (United States)

Chang, Shih-Hua

2018-04-01

The purpose of this study was to test a model of codependency based on Bowen's concept of differentiation for college students in Taiwan. The relations between family-of-origin dysfunction, differentiation of self, codependency traits and related symptoms including low self-esteem, relationship distress and psychological adjustment problems were examined. Data were collected from 567 college students from 2 large, urban universities in northern Taiwan. Results indicated a significantly negative relationship between levels of codependency and self-differentiation and that self-differentiation partially mediated the relationship between family-of-origin dysfunction and codependency. The implications of these findings for counselling Taiwanese college students who experience codependency traits and related symptoms as well as suggestions for future research are discussed. © 2016 International Union of Psychological Science.

Preconditioning for partial differential equation constrained optimization with control constraints

KAUST Repository

Stoll, Martin

2011-10-18

Optimal control problems with partial differential equations play an important role in many applications. The inclusion of bound constraints for the control poses a significant additional challenge for optimization methods. In this paper, we propose preconditioners for the saddle point problems that arise when a primal-dual active set method is used. We also show for this method that the same saddle point system can be derived when the method is considered as a semismooth Newton method. In addition, the projected gradient method can be employed to solve optimization problems with simple bounds, and we discuss the efficient solution of the linear systems in question. In the case when an acceleration technique is employed for the projected gradient method, this again yields a semismooth Newton method that is equivalent to the primal-dual active set method. We also consider the Moreau-Yosida regularization method for control constraints and efficient preconditioners for this technique. Numerical results illustrate the competitiveness of these approaches. © 2011 John Wiley & Sons, Ltd.

Partial differential equations an accessible route through theory and applications

CERN Document Server

Vasy, András

2015-01-01

This text on partial differential equations is intended for readers who want to understand the theoretical underpinnings of modern PDEs in settings that are important for the applications without using extensive analytic tools required by most advanced texts. The assumed mathematical background is at the level of multivariable calculus and basic metric space material, but the latter is recalled as relevant as the text progresses. The key goal of this book is to be mathematically complete without overwhelming the reader, and to develop PDE theory in a manner that reflects how researchers would think about the material. A concrete example is that distribution theory and the concept of weak solutions are introduced early because while these ideas take some time for the students to get used to, they are fundamentally easy and, on the other hand, play a central role in the field. Then, Hilbert spaces that are quite important in the later development are introduced via completions which give essentially all the fea...

Preconditioning for partial differential equation constrained optimization with control constraints

KAUST Repository

Stoll, Martin; Wathen, Andy

2011-01-01

Optimal control problems with partial differential equations play an important role in many applications. The inclusion of bound constraints for the control poses a significant additional challenge for optimization methods. In this paper, we propose preconditioners for the saddle point problems that arise when a primal-dual active set method is used. We also show for this method that the same saddle point system can be derived when the method is considered as a semismooth Newton method. In addition, the projected gradient method can be employed to solve optimization problems with simple bounds, and we discuss the efficient solution of the linear systems in question. In the case when an acceleration technique is employed for the projected gradient method, this again yields a semismooth Newton method that is equivalent to the primal-dual active set method. We also consider the Moreau-Yosida regularization method for control constraints and efficient preconditioners for this technique. Numerical results illustrate the competitiveness of these approaches. © 2011 John Wiley & Sons, Ltd.

Some Considerations on the Partial Credit Model

Science.gov (United States)

Verhelst, N. D.; Verstralen, H. H. F. M.

2008-01-01

The Partial Credit Model (PCM) is sometimes interpreted as a model for stepwise solution of polytomously scored items, where the item parameters are interpreted as difficulties of the steps. It is argued that this interpretation is not justified. A model for stepwise solution is discussed. It is shown that the PCM is suited to model sums of binary…

Modulating functions-based method for parameters and source estimation in one-dimensional partial differential equations

KAUST Repository

Asiri, Sharefa M.

2016-10-20

In this paper, modulating functions-based method is proposed for estimating space–time-dependent unknowns in one-dimensional partial differential equations. The proposed method simplifies the problem into a system of algebraic equations linear in unknown parameters. The well-posedness of the modulating functions-based solution is proved. The wave and the fifth-order KdV equations are used as examples to show the effectiveness of the proposed method in both noise-free and noisy cases.

From ordinary to partial differential equations

CERN Document Server

Esposito, Giampiero

2017-01-01

This book is addressed to mathematics and physics students who want to develop an interdisciplinary view of mathematics, from the age of Riemann, Poincaré and Darboux to basic tools of modern mathematics. It enables them to acquire the sensibility necessary for the formulation and solution of difficult problems, with an emphasis on concepts, rigour and creativity. It consists of eight self-contained parts: ordinary differential equations; linear elliptic equations; calculus of variations; linear and non-linear hyperbolic equations; parabolic equations; Fuchsian functions and non-linear equations; the functional equations of number theory; pseudo-differential operators and pseudo-differential equations. The author leads readers through the original papers and introduces new concepts, with a selection of topics and examples that are of high pedagogical value.

Analytical solutions to time-fractional partial differential equations in a two-dimensional multilayer annulus

Science.gov (United States)

Chen, Shanzhen; Jiang, Xiaoyun

2012-08-01

In this paper, analytical solutions to time-fractional partial differential equations in a multi-layer annulus are presented. The final solutions are obtained in terms of Mittag-Leffler function by using the finite integral transform technique and Laplace transform technique. In addition, the classical diffusion equation (α=1), the Helmholtz equation (α→0) and the wave equation (α=2) are discussed as special cases. Finally, an illustrative example problem for the three-layer semi-circular annular region is solved and numerical results are presented graphically for various kind of order of fractional derivative.

Testing the generalized partial credit model

OpenAIRE

Glas, Cornelis A.W.

1996-01-01

The partial credit model (PCM) (G.N. Masters, 1982) can be viewed as a generalization of the Rasch model for dichotomous items to the case of polytomous items. In many cases, the PCM is too restrictive to fit the data. Several generalizations of the PCM have been proposed. In this paper, a generalization of the PCM (GPCM), a further generalization of the one-parameter logistic model, is discussed. The model is defined and the conditional maximum likelihood procedure for the method is describe...

Polynomial chaos methods for hyperbolic partial differential equations numerical techniques for fluid dynamics problems in the presence of uncertainties

CERN Document Server

Pettersson, Mass Per; Nordström, Jan

2015-01-01

This monograph presents computational techniques and numerical analysis to study conservation laws under uncertainty using the stochastic Galerkin formulation. With the continual growth of computer power, these methods are becoming increasingly popular as an alternative to more classical sampling-based techniques. The approach described in the text takes advantage of stochastic Galerkin projections applied to the original conservation laws to produce a large system of modified partial differential equations, the solutions to which directly provide a full statistical characterization of the effect of uncertainties. Polynomial Chaos Methods of Hyperbolic Partial Differential Equations focuses on the analysis of stochastic Galerkin systems obtained for linear and non-linear convection-diffusion equations and for a systems of conservation laws; a detailed well-posedness and accuracy analysis is presented to enable the design of robust and stable numerical methods. The exposition is restricted to one spatial dime...

An Efficient Numerical Approach for Solving Nonlinear Coupled Hyperbolic Partial Differential Equations with Nonlocal Conditions

Directory of Open Access Journals (Sweden)

A. H. Bhrawy

2014-01-01

Full Text Available One of the most important advantages of collocation method is the possibility of dealing with nonlinear partial differential equations (PDEs as well as PDEs with variable coefficients. A numerical solution based on a Jacobi collocation method is extended to solve nonlinear coupled hyperbolic PDEs with variable coefficients subject to initial-boundary nonlocal conservation conditions. This approach, based on Jacobi polynomials and Gauss-Lobatto quadrature integration, reduces solving the nonlinear coupled hyperbolic PDEs with variable coefficients to a system of nonlinear ordinary differential equation which is far easier to solve. In fact, we deal with initial-boundary coupled hyperbolic PDEs with variable coefficients as well as initial-nonlocal conditions. Using triangular, soliton, and exponential-triangular solutions as exact solutions, the obtained results show that the proposed numerical algorithm is efficient and very accurate.

Mathematical and numerical methods for partial differential equations applications for engineering sciences

CERN Document Server

Chaskalovic, Joël

2014-01-01

This self-tutorial offers a concise yet thorough introduction into the mathematical analysis of approximation methods for partial differential equation. A particular emphasis is put on finite element methods. The unique approach first summarizes and outlines the finite-element mathematics in general and then, in the second and major part, formulates problem examples that clearly demonstrate the techniques of functional analysis via numerous and diverse exercises. The solutions of the problems are given directly afterwards. Using this approach, the author motivates and encourages the reader to actively acquire the knowledge of finite- element methods instead of passively absorbing the material, as in most standard textbooks. This English edition is based on the Finite Element Methods for Engineering Sciences by Joel Chaskalovic

CPDS3, Coupled 3-D Partial Differential Equation Solution

International Nuclear Information System (INIS)

Anderson, D.V.; Koniges, A.E.; Shumaker, D.E.

1992-01-01

1 - Description of program or function: CPDES3 solves the linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils, permits general couplings between all of the component PDE's, and automatically generates the matrix structures needed to perform the algorithm. 2 - Method of solution: The resulting sparse matrix equation with a complicated sub-band structure and generally asymmetric is solved by either the preconditioned conjugate gradient (CG) method or the preconditioned bi-conjugate gradient (BCG) algorithm. BCG enjoys faster convergence in most cases but in rare instances diverges. Then, CG iterations must be used. 3 - Restrictions on the complexity of the problem: The discretization of the coupled three-dimensional PDE's and their boundary conditions must result in an operator stencil which fits in the Cray2 memory. In addition, the matrix must possess a reasonable amount of diagonal dominance for the preconditioning technique to be effective

CPDES2, Coupled 2-D Partial Differential Equation Solution

International Nuclear Information System (INIS)

Anderson, D.V.; Koniges, A.E.; Shumaker, D.E.

1992-01-01

1 - Description of program or function: CPDES2 solves the linear asymmetric equations arising from coupled partial differential equations in two dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES2 allows each spatial operator to have 5 or 9 point stencils, permits general coupling between all of the component PDE's, and automatically generates the matrix structures needed to perform the algorithm. 2 - Method of solution: The resulting sparse matrix equation with a complicated sub-band structure and generally asymmetric is solved by either the preconditioned conjugate gradient (CG) method or the preconditioned bi-conjugate gradient (BCG) algorithm. BCG enjoys faster convergence in most cases but in rare instances diverges. Then, CG iterations must be used. 3 - Restrictions on the complexity of the problem: The discretization of the coupled two-dimensional PDE's and their boundary conditions must result in an operator stencil which fits in the Cray2 memory. In addition, the matrix must possess a reasonable amount of diagonal dominance for the preconditioning technique to be effective

A hybrid perturbation-Galerkin technique for partial differential equations

Science.gov (United States)

Geer, James F.; Anderson, Carl M.

1990-01-01

A two-step hybrid perturbation-Galerkin technique for improving the usefulness of perturbation solutions to partial differential equations which contain a parameter is presented and discussed. In the first step of the method, the leading terms in the asymptotic expansion(s) of the solution about one or more values of the perturbation parameter are obtained using standard perturbation methods. In the second step, the perturbation functions obtained in the first step are used as trial functions in a Bubnov-Galerkin approximation. This semi-analytical, semi-numerical hybrid technique appears to overcome some of the drawbacks of the perturbation and Galerkin methods when they are applied by themselves, while combining some of the good features of each. The technique is illustrated first by a simple example. It is then applied to the problem of determining the flow of a slightly compressible fluid past a circular cylinder and to the problem of determining the shape of a free surface due to a sink above the surface. Solutions obtained by the hybrid method are compared with other approximate solutions, and its possible application to certain problems associated with domain decomposition is discussed.

Scale-invariant solutions to partial differential equations of fractional order with a moving boundary condition

International Nuclear Information System (INIS)

Li Xicheng; Xu Mingyu; Wang Shaowei

2008-01-01

In this paper, we give similarity solutions of partial differential equations of fractional order with a moving boundary condition. The solutions are given in terms of a generalized Wright function. The time-fractional Caputo derivative and two types of space-fractional derivatives are considered. The scale-invariant variable and the form of the solution of the moving boundary are obtained by the Lie group analysis. A comparison between the solutions corresponding to two types of fractional derivative is also given

High-order fractional partial differential equation transform for molecular surface construction.

Science.gov (United States)

Hu, Langhua; Chen, Duan; Wei, Guo-Wei

2013-01-01

Fractional derivative or fractional calculus plays a significant role in theoretical modeling of scientific and engineering problems. However, only relatively low order fractional derivatives are used at present. In general, it is not obvious what role a high fractional derivative can play and how to make use of arbitrarily high-order fractional derivatives. This work introduces arbitrarily high-order fractional partial differential equations (PDEs) to describe fractional hyperdiffusions. The fractional PDEs are constructed via fractional variational principle. A fast fractional Fourier transform (FFFT) is proposed to numerically integrate the high-order fractional PDEs so as to avoid stringent stability constraints in solving high-order evolution PDEs. The proposed high-order fractional PDEs are applied to the surface generation of proteins. We first validate the proposed method with a variety of test examples in two and three-dimensional settings. The impact of high-order fractional derivatives to surface analysis is examined. We also construct fractional PDE transform based on arbitrarily high-order fractional PDEs. We demonstrate that the use of arbitrarily high-order derivatives gives rise to time-frequency localization, the control of the spectral distribution, and the regulation of the spatial resolution in the fractional PDE transform. Consequently, the fractional PDE transform enables the mode decomposition of images, signals, and surfaces. The effect of the propagation time on the quality of resulting molecular surfaces is also studied. Computational efficiency of the present surface generation method is compared with the MSMS approach in Cartesian representation. We further validate the present method by examining some benchmark indicators of macromolecular surfaces, i.e., surface area, surface enclosed volume, surface electrostatic potential and solvation free energy. Extensive numerical experiments and comparison with an established surface model

ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY

OpenAIRE

Enrique Gonzalo Reyes Garcia

2004-01-01

ON DIFFERENTIAL EQUATIONS, INTEGRABLE SYSTEMS, AND GEOMETRY Equations in partial derivatives appeared in the 18th century as essential tools for the analytic study of physical models and, later, they proved to be fundamental for the progress of mathematics. For example, fundamental results of modern differential geometry are based on deep theorems on differential equations. Reciprocally, it is possible to study differential equations through geometrical means just like it was done by o...

Survival of partially differentiated mouse embryonic stem cells in the scala media of the guinea pig cochlea.

Science.gov (United States)

Hildebrand, Michael S; Dahl, Hans-Henrik M; Hardman, Jennifer; Coleman, Bryony; Shepherd, Robert K; de Silva, Michelle G

2005-12-01

The low regenerative capacity of the hair cells of the mammalian inner ear is a major obstacle for functional recovery following sensorineural hearing loss. A potential treatment is to replace damaged tissue by transplantation of stem cells. To test this approach, undifferentiated and partially differentiated mouse embryonic stem (ES) cells were delivered into the scala media of the deafened guinea pig cochlea. Transplanted cells survived in the scala media for a postoperative period of at least nine weeks, evidenced by histochemical and direct fluorescent detection of enhanced green fluorescent protein (EGFP). Transplanted cells were discovered near the spiral ligament and stria vascularis in the endolymph fluid of the scala media. In some cases, cells were observed close to the damaged organ of Corti structure. There was no evidence of significant immunological rejection of the implanted ES cells despite the absence of immunosuppression. Our surgical approach allowed efficient delivery of ES cells to the scala media while preserving the delicate structures of the cochlea. This is the first report of the survival of partially differentiated ES cells in the scala media of the mammalian cochlea, and it provides support for the potential of cell-based therapies for sensorineural hearing impairment.

Partial-Order Reduction for GPU Model Checking

NARCIS (Netherlands)

Neele, T.; Wijs, A.; Bosnacki, D.; van de Pol, Jan Cornelis; Artho, C; Legay, A.; Peled, D.

2016-01-01

Model checking using GPUs has seen increased popularity over the last years. Because GPUs have a limited amount of memory, only small to medium-sized systems can be verified. For on-the-fly explicit-state model checking, we improve memory efficiency by applying partial-order reduction. We propose

Three tesla magnetic resonance imaging of the anterior cruciate ligament of the knee: can we differentiate complete from partial tears?

Energy Technology Data Exchange (ETDEWEB)

Dyck, Pieter van; Gielen, Jan L.; Parizel, Paul M. [University Hospital Antwerp and University of Antwerp, Department of Radiology, Antwerp (Edegem) (Belgium); Vanhoenacker, Filip M. [University Hospital Antwerp and University of Antwerp, Department of Radiology, Antwerp (Edegem) (Belgium); AZ St-Maarten Duffel/Mechelen, Department of Radiology, Duffel (Belgium); Dossche, Lieven; Gestel, Jozef van [University Hospital Antwerp and University of Antwerp, Department of Orthopedics, Antwerp (Edegem) (Belgium); Wouters, Kristien [University Hospital Antwerp and University of Antwerp, Department of Scientific Coordination and Biostatistics, Antwerp (Edegem) (Belgium)

2011-06-15

To determine the ability of 3.0T magnetic resonance (MR) imaging to identify partial tears of the anterior cruciate ligament (ACL) and to allow distinction of complete from partial ACL tears. One hundred seventy-two patients were prospectively studied by 3.0T MR imaging and arthroscopy in our institution. MR images were interpreted in consensus by two experienced reviewers, and the ACL was diagnosed as being normal, partially torn, or completely torn. Diagnostic accuracy of 3.0T MR for the detection of both complete and partial tears of the ACL was calculated using arthroscopy as the standard of reference. There were 132 patients with an intact ACL, 17 had a partial, and 23 had a complete tear of the ACL seen at arthroscopy. Sensitivity, specificity, and accuracy of 3.0T MR for complete ACL tears were 83, 99, and 97%, respectively, and, for partial ACL tears, 77, 97, and 95%, respectively. Five of 40 ACL lesions (13%) could not correctly be identified as complete or partial ACL tears. MR imaging at 3.0T represents a highly accurate method for identifying tears of the ACL. However, differentiation between complete and partial ACL tears and identification of partial tears of this ligament remains difficult, even at 3.0T. (orig.)

Three tesla magnetic resonance imaging of the anterior cruciate ligament of the knee: can we differentiate complete from partial tears?

International Nuclear Information System (INIS)

Dyck, Pieter van; Gielen, Jan L.; Parizel, Paul M.; Vanhoenacker, Filip M.; Dossche, Lieven; Gestel, Jozef van; Wouters, Kristien

2011-01-01

To determine the ability of 3.0T magnetic resonance (MR) imaging to identify partial tears of the anterior cruciate ligament (ACL) and to allow distinction of complete from partial ACL tears. One hundred seventy-two patients were prospectively studied by 3.0T MR imaging and arthroscopy in our institution. MR images were interpreted in consensus by two experienced reviewers, and the ACL was diagnosed as being normal, partially torn, or completely torn. Diagnostic accuracy of 3.0T MR for the detection of both complete and partial tears of the ACL was calculated using arthroscopy as the standard of reference. There were 132 patients with an intact ACL, 17 had a partial, and 23 had a complete tear of the ACL seen at arthroscopy. Sensitivity, specificity, and accuracy of 3.0T MR for complete ACL tears were 83, 99, and 97%, respectively, and, for partial ACL tears, 77, 97, and 95%, respectively. Five of 40 ACL lesions (13%) could not correctly be identified as complete or partial ACL tears. MR imaging at 3.0T represents a highly accurate method for identifying tears of the ACL. However, differentiation between complete and partial ACL tears and identification of partial tears of this ligament remains difficult, even at 3.0T. (orig.)

Flow visualization via partial differential equations

NARCIS (Netherlands)

Preusser, T.; Rumpf, M.; Telea, A.C.; Möller, T.; Hamann, B.; Russell, R.D.

2009-01-01

The visualization of stationary and time-dependent flow is an important and chaltenging topic in scientific visualization. lts aim is 10 represent transport phenomena govemed by vector fjelds in an intuitively understandable way. In this paper. we review the use of methods based on partial

ADI splitting schemes for a fourth-order nonlinear partial differential equation from image processing

KAUST Repository

Calatroni, Luca

2013-08-01

We present directional operator splitting schemes for the numerical solution of a fourth-order, nonlinear partial differential evolution equation which arises in image processing. This equation constitutes the H -1-gradient flow of the total variation and represents a prototype of higher-order equations of similar type which are popular in imaging for denoising, deblurring and inpainting problems. The efficient numerical solution of this equation is very challenging due to the stiffness of most numerical schemes. We show that the combination of directional splitting schemes with implicit time-stepping provides a stable and computationally cheap numerical realisation of the equation.

ADI splitting schemes for a fourth-order nonlinear partial differential equation from image processing

KAUST Repository

Calatroni, Luca; Dü ring, Bertram; Schö nlieb, Carola-Bibiane

2013-01-01

We present directional operator splitting schemes for the numerical solution of a fourth-order, nonlinear partial differential evolution equation which arises in image processing. This equation constitutes the H -1-gradient flow of the total variation and represents a prototype of higher-order equations of similar type which are popular in imaging for denoising, deblurring and inpainting problems. The efficient numerical solution of this equation is very challenging due to the stiffness of most numerical schemes. We show that the combination of directional splitting schemes with implicit time-stepping provides a stable and computationally cheap numerical realisation of the equation.

First-order systems of linear partial differential equations: normal forms, canonical systems, transform methods

Directory of Open Access Journals (Sweden)

Heinz Toparkus

2014-04-01

Full Text Available In this paper we consider first-order systems with constant coefficients for two real-valued functions of two real variables. This is both a problem in itself, as well as an alternative view of the classical linear partial differential equations of second order with constant coefficients. The classification of the systems is done using elementary methods of linear algebra. Each type presents its special canonical form in the associated characteristic coordinate system. Then you can formulate initial value problems in appropriate basic areas, and you can try to achieve a solution of these problems by means of transform methods.

New Traveling Wave Solutions of the Higher Dimensional Nonlinear Partial Differential Equation by the Exp-Function Method

Directory of Open Access Journals (Sweden)

Hasibun Naher

2012-01-01

Full Text Available We construct new analytical solutions of the (3+1-dimensional modified KdV-Zakharov-Kuznetsev equation by the Exp-function method. Plentiful exact traveling wave solutions with arbitrary parameters are effectively obtained by the method. The obtained results show that the Exp-function method is effective and straightforward mathematical tool for searching analytical solutions with arbitrary parameters of higher-dimensional nonlinear partial differential equation.

A Differential Evolution Based MPPT Method for Photovoltaic Modules under Partial Shading Conditions

Directory of Open Access Journals (Sweden)

Kok Soon Tey

2014-01-01

Full Text Available Partially shaded photovoltaic (PV modules have multiple peaks in the power-voltage (P-V characteristic curve and conventional maximum power point tracking (MPPT algorithm, such as perturbation and observation (P&O, which is unable to track the global maximum power point (GMPP accurately due to its localized search space. Therefore, this paper proposes a differential evolution (DE based optimization algorithm to provide the globalized search space to track the GMPP. The direction of mutation in the DE algorithm is modified to ensure that the mutation always converges to the best solution among all the particles in the generation. This helps to provide the rapid convergence of the algorithm. Simulation of the proposed PV system is carried out in PSIM and the results are compared to P&O algorithm. In the hardware implementation, a high step-up DC-DC converter is employed to verify the proposed algorithm experimentally on partial shading conditions, load variation, and solar intensity variation. The experimental results show that the proposed algorithm is able to converge to the GMPP within 1.2 seconds with higher efficiency than P&O.

Issues in developing parallel iterative algorithms for solving partial differential equations on a (transputer-based) distributed parallel computing system

International Nuclear Information System (INIS)

Rajagopalan, S.; Jethra, A.; Khare, A.N.; Ghodgaonkar, M.D.; Srivenkateshan, R.; Menon, S.V.G.

1990-01-01

Issues relating to implementing iterative procedures, for numerical solution of elliptic partial differential equations, on a distributed parallel computing system are discussed. Preliminary investigations show that a speed-up of about 3.85 is achievable on a four transputer pipeline network. (author). 2 figs., 3 a ppendixes., 7 refs

Partial chord diagrams and matrix models

DEFF Research Database (Denmark)

Andersen, Jørgen Ellegaard; Fuji, Hiroyuki; Manabe, Masahide

In this article, the enumeration of partial chord diagrams is discussed via matrix model techniques. In addition to the basic data such as the number of backbones and chords, we also consider the Euler characteristic, the backbone spectrum, the boundary point spectrum, and the boundary length spe...

Patch Similarity Modulus and Difference Curvature Based Fourth-Order Partial Differential Equation for Image Denoising

Directory of Open Access Journals (Sweden)

Yunjiao Bai

2015-01-01

Full Text Available The traditional fourth-order nonlinear diffusion denoising model suffers the isolated speckles and the loss of fine details in the processed image. For this reason, a new fourth-order partial differential equation based on the patch similarity modulus and the difference curvature is proposed for image denoising. First, based on the intensity similarity of neighbor pixels, this paper presents a new edge indicator called patch similarity modulus, which is strongly robust to noise. Furthermore, the difference curvature which can effectively distinguish between edges and noise is incorporated into the denoising algorithm to determine the diffusion process by adaptively adjusting the size of the diffusion coefficient. The experimental results show that the proposed algorithm can not only preserve edges and texture details, but also avoid isolated speckles and staircase effect while filtering out noise. And the proposed algorithm has a better performance for the images with abundant details. Additionally, the subjective visual quality and objective evaluation index of the denoised image obtained by the proposed algorithm are higher than the ones from the related methods.

Ultrasound speckle reduction based on fractional order differentiation.

Science.gov (United States)

Shao, Dangguo; Zhou, Ting; Liu, Fan; Yi, Sanli; Xiang, Yan; Ma, Lei; Xiong, Xin; He, Jianfeng

2017-07-01

Ultrasound images show a granular pattern of noise known as speckle that diminishes their quality and results in difficulties in diagnosis. To preserve edges and features, this paper proposes a fractional differentiation-based image operator to reduce speckle in ultrasound. An image de-noising model based on fractional partial differential equations with balance relation between k (gradient modulus threshold that controls the conduction) and v (the order of fractional differentiation) was constructed by the effective combination of fractional calculus theory and a partial differential equation, and the numerical algorithm of it was achieved using a fractional differential mask operator. The proposed algorithm has better speckle reduction and structure preservation than the three existing methods [P-M model, the speckle reducing anisotropic diffusion (SRAD) technique, and the detail preserving anisotropic diffusion (DPAD) technique]. And it is significantly faster than bilateral filtering (BF) in producing virtually the same experimental results. Ultrasound phantom testing and in vivo imaging show that the proposed method can improve the quality of an ultrasound image in terms of tissue SNR, CNR, and FOM values.

Harmonic analysis, partial differential equations and applications in honor of Richard L. Wheeden

CERN Document Server

Franchi, Bruno; Lu, Guozhen; Perez, Carlos; Sawyer, Eric

2017-01-01

This is a collection of contributed papers by many eminent Harmonic Analysts and specialists of Partial Differential equations. The papers focus on weighted norm equalities for singular integrals, focusing wave equations, degenerate elliptic equations, Navier-Stokes flow in two dimensions and Poincare-Sobolev inequalities in the setting of metric spaces equipped with measures among others. Many topics considered in this volume stem from the interests of Richard L. Wheeden whose contributions to Potential Theory, singular integral theory and degenerate elliptic PDE theory this volume honors. Luis Caffarelli, Sagun Chanillo, Bruno Franchi, Cristian Guttierez, Xiaojun Huang, Carlos Kenig, Ermanno Lanconelli, Eric Sawyer and Alexander Volberg, are some of the many contributors to this volume. .

On Direct Transformation Approach to Asymptotical Analytical Solutions of Perturbed Partial Differential Equation

International Nuclear Information System (INIS)

Liu Hongzhun; Pan Zuliang; Li Peng

2006-01-01

In this article, we will derive an equality, where the Taylor series expansion around ε = 0 for any asymptotical analytical solution of the perturbed partial differential equation (PDE) with perturbing parameter ε must be admitted. By making use of the equality, we may obtain a transformation, which directly map the analytical solutions of a given unperturbed PDE to the asymptotical analytical solutions of the corresponding perturbed one. The notion of Lie-Baecklund symmetries is introduced in order to obtain more transformations. Hence, we can directly create more transformations in virtue of known Lie-Baecklund symmetries and recursion operators of corresponding unperturbed equation. The perturbed Burgers equation and the perturbed Korteweg-de Vries (KdV) equation are used as examples.

Bayesian inference for partially identified models exploring the limits of limited data

CERN Document Server

Gustafson, Paul

2015-01-01

Introduction Identification What Is against Us? What Is for Us? Some Simple Examples of Partially Identified ModelsThe Road Ahead The Structure of Inference in Partially Identified Models Bayesian Inference The Structure of Posterior Distributions in PIMs Computational Strategies Strength of Bayesian Updating, Revisited Posterior MomentsCredible Intervals Evaluating the Worth of Inference Partial Identification versus Model Misspecification The Siren Call of Identification Comp

Coenzyme Q10 partially restores pathological alterations in a macrophage model of Gaucher disease.

Science.gov (United States)

de la Mata, Mario; Cotán, David; Oropesa-Ávila, Manuel; Villanueva-Paz, Marina; de Lavera, Isabel; Álvarez-Córdoba, Mónica; Luzón-Hidalgo, Raquel; Suárez-Rivero, Juan M; Tiscornia, Gustavo; Sánchez-Alcázar, José A

2017-02-06

Gaucher disease (GD) is caused by mutations in the GBA1 gene which encodes lysosomal β-glucocerebrosidase (GCase). In GD, partial or complete loss of GCase activity causes the accumulation of the glycolipids glucosylceramide (GlcCer) and glucosylsphingosine in the lysosomes of macrophages. In this manuscript, we investigated the effects of glycolipids accumulation on lysosomal and mitochondrial function, inflammasome activation and efferocytosis capacity in a THP-1 macrophage model of Gaucher disease. In addition, the beneficial effects of coenzyme Q 10 (CoQ) supplementation on cellular alterations were evaluated. Chemically-induced Gaucher macrophages were developed by differentiateing THP-1 monocytes to macrophages by treatment with phorbol 12-myristate 13-acetate (PMA) and then inhibiting intracellular GCase with conduritol B-epoxide (CBE), a specific irreversible inhibitor of GCase activity, and supplementing the medium with exogenous GlcCer. This cell model accumulated up to 16-fold more GlcCer compared with control THP-1 cells. Chemically-induced Gaucher macrophages showed impaired autophagy flux associated with mitochondrial dysfunction and increased oxidative stress, inflammasome activation and impaired efferocytosis. All abnormalities were partially restored by supplementation with CoQ. These data suggest that targeting mitochondria function and oxidative stress by CoQ can ameliorate the pathological phenotype of Gaucher cells. Chemically-induced Gaucher macrophages provide cellular models that can be used to investigate disease pathogenesis and explore new therapeutics for GD.

THREE-POINT BACKWARD FINITE DIFFERENCE METHOD FOR SOLVING A SYSTEM OF MIXED HYPERBOLIC-PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS. (R825549C019)

Science.gov (United States)

A three-point backward finite-difference method has been derived for a system of mixed hyperbolic¯¯parabolic (convection¯¯diffusion) partial differential equations (mixed PDEs). The method resorts to the three-point backward differenci...

Efficient solution of ordinary differential equations modeling electrical activity in cardiac cells.

Science.gov (United States)

Sundnes, J; Lines, G T; Tveito, A

2001-08-01

The contraction of the heart is preceded and caused by a cellular electro-chemical reaction, causing an electrical field to be generated. Performing realistic computer simulations of this process involves solving a set of partial differential equations, as well as a large number of ordinary differential equations (ODEs) characterizing the reactive behavior of the cardiac tissue. Experiments have shown that the solution of the ODEs contribute significantly to the total work of a simulation, and there is thus a strong need to utilize efficient solution methods for this part of the problem. This paper presents how an efficient implicit Runge-Kutta method may be adapted to solve a complicated cardiac cell model consisting of 31 ODEs, and how this solver may be coupled to a set of PDE solvers to provide complete simulations of the electrical activity.

Partial differential equations of first order and their applications to physics

CERN Document Server

López, Gustavo

2012-01-01

This book tries to point out the mathematical importance of the Partial Differential Equations of First Order (PDEFO) in Physics and Applied Sciences. The intention is to provide mathematicians with a wide view of the applications of this branch in physics, and to give physicists and applied scientists a powerful tool for solving some problems appearing in Classical Mechanics, Quantum Mechanics, Optics, and General Relativity. This book is intended for senior or first year graduate students in mathematics, physics, or engineering curricula. This book is unique in the sense that it covers the applications of PDEFO in several branches of applied mathematics, and fills the theoretical gap between the formal mathematical presentation of the theory and the pure applied tool to physical problems that are contained in other books. Improvements made in this second edition include corrected typographical errors; rewritten text to improve the flow and enrich the material; added exercises in all chapters; new applicati...

Sparse dynamics for partial differential equations.

Science.gov (United States)

Schaeffer, Hayden; Caflisch, Russel; Hauck, Cory D; Osher, Stanley

2013-04-23

We investigate the approximate dynamics of several differential equations when the solutions are restricted to a sparse subset of a given basis. The restriction is enforced at every time step by simply applying soft thresholding to the coefficients of the basis approximation. By reducing or compressing the information needed to represent the solution at every step, only the essential dynamics are represented. In many cases, there are natural bases derived from the differential equations, which promote sparsity. We find that our method successfully reduces the dynamics of convection equations, diffusion equations, weak shocks, and vorticity equations with high-frequency source terms.

Some Considerations on the Partial Credit Model

OpenAIRE

H.H.F.M. Verstralen; N.D. Verhelst

2008-01-01

The Partial Credit Model (PCM) is sometimes interpreted as a model for stepwise solution of polytomously scored items, where the item parameters are interpreted as di culties of the steps. It is argued that this interpretation is not justi ed. A model for stepwise solution is discussed. It is shown that the PCM is suited to model sums of binary responses which are not supposed to be stochastically independent. As a practical result, a statistical test of sto...

Mathematical analysis of partial differential equations modeling electrostatic MEMS

CERN Document Server

Esposito, Pierpaolo; Guo, Yujin

2010-01-01

Micro- and nanoelectromechanical systems (MEMS and NEMS), which combine electronics with miniature-size mechanical devices, are essential components of modern technology. It is the mathematical model describing "electrostatically actuated" MEMS that is addressed in this monograph. Even the simplified models that the authors deal with still lead to very interesting second- and fourth-order nonlinear elliptic equations (in the stationary case) and to nonlinear parabolic equations (in the dynamic case). While nonlinear eigenvalue problems-where the stationary MEMS models fit-are a well-developed

Differential geometry based solvation model II: Lagrangian formulation.

Science.gov (United States)

Chen, Zhan; Baker, Nathan A; Wei, G W

2011-12-01

computation, thanks to the equivalence of the Laplace-Beltrami operator in the two representations. The coupled partial differential equations (PDEs) are solved with an iterative procedure to reach a steady state, which delivers desired solvent-solute interface and electrostatic potential for problems of interest. These quantities are utilized to evaluate the solvation free energies and protein-protein binding affinities. A number of computational methods and algorithms are described for the interconversion of Lagrangian and Eulerian representations, and for the solution of the coupled PDE system. The proposed approaches have been extensively validated. We also verify that the mean curvature flow indeed gives rise to the minimal molecular surface and the proposed variational procedure indeed offers minimal total free energy. Solvation analysis and applications are considered for a set of 17 small compounds and a set of 23 proteins. The salt effect on protein-protein binding affinity is investigated with two protein complexes by using the present model. Numerical results are compared to the experimental measurements and to those obtained by using other theoretical methods in the literature. © Springer-Verlag 2011

Calculus for cognitive scientists partial differential equation models

CERN Document Server

Peterson, James K

2016-01-01

This book shows cognitive scientists in training how mathematics, computer science and science can be usefully and seamlessly intertwined. It is a follow-up to the first two volumes on mathematics for cognitive scientists, and includes the mathematics and computational tools needed to understand how to compute the terms in the Fourier series expansions that solve the cable equation. The latter is derived from first principles by going back to cellular biology and the relevant biophysics.  A detailed discussion of ion movement through cellular membranes, and an explanation of how the equations that govern such ion movement leading to the standard transient cable equation are included. There are also solutions for the cable model using separation of variables, as well an explanation of why Fourier series converge and a description of the implementation of MatLab tools to compute the solutions. Finally, the standard Hodgkin - Huxley model is developed for an excitable neuron and is solved using MatLab.

Dimensional analysis to transform the differential equations in partial derivates in the theory of heat transmission into ordinary ones

International Nuclear Information System (INIS)

Diaz Sanchidrian, C.

1989-01-01

The present paper applies dimensional analysis with spatial discrimination to transform the differential equations in partial derivatives developed in the theory of heat transmission into ordinary ones. The effectivity of the method is comparable to that methods based in transformations of uni or multiparametric groups, with the advantage of being more direct and simple. (Author)

Combination of oriented partial differential equation and shearlet transform for denoising in electronic speckle pattern interferometry fringe patterns.

Science.gov (United States)

Xu, Wenjun; Tang, Chen; Gu, Fan; Cheng, Jiajia

2017-04-01

It is a key step to remove the massive speckle noise in electronic speckle pattern interferometry (ESPI) fringe patterns. In the spatial-domain filtering methods, oriented partial differential equations have been demonstrated to be a powerful tool. In the transform-domain filtering methods, the shearlet transform is a state-of-the-art method. In this paper, we propose a filtering method for ESPI fringe patterns denoising, which is a combination of second-order oriented partial differential equation (SOOPDE) and the shearlet transform, named SOOPDE-Shearlet. Here, the shearlet transform is introduced into the ESPI fringe patterns denoising for the first time. This combination takes advantage of the fact that the spatial-domain filtering method SOOPDE and the transform-domain filtering method shearlet transform benefit from each other. We test the proposed SOOPDE-Shearlet on five experimentally obtained ESPI fringe patterns with poor quality and compare our method with SOOPDE, shearlet transform, windowed Fourier filtering (WFF), and coherence-enhancing diffusion (CEDPDE). Among them, WFF and CEDPDE are the state-of-the-art methods for ESPI fringe patterns denoising in transform domain and spatial domain, respectively. The experimental results have demonstrated the good performance of the proposed SOOPDE-Shearlet.

Response Styles in the Partial Credit Model

OpenAIRE

Tutz, Gerhard; Schauberger, Gunther; Berger, Moritz

2016-01-01

In the modelling of ordinal responses in psychological measurement and survey- based research, response styles that represent specific answering patterns of respondents are typically ignored. One consequence is that estimates of item parameters can be poor and considerably biased. The focus here is on the modelling of a tendency to extreme or middle categories. An extension of the Partial Credit Model is proposed that explicitly accounts for this specific response style. In contrast to exi...

MODEL PENSKORAN PARTIAL CREDIT PADA BUTIR MULTIPLE TRUE-FALSE BIDANG FISIKA

Directory of Open Access Journals (Sweden)

Wasis Wasis

2013-01-01

Full Text Available Tujuan penelitian ini menghasilkan model penskoran politomus untuk respons butir multiple true-false, sehingga dapat mengestimasi secara lebih akurat kemampuan di bidang fisika. Pengembangan penskoran menggunakan Four-D model dan diuji akurasinya melalui penelitian empiris dan simulasi. Penelitian empiris menggunakan 15 butir multiple true-false yang diambil dari soal UMPTN tahun 1996-2006 dan dikenakan pada 410 mahasiswa baru FMIPA Universitas Negeri Surabaya angkatan tahun 2007. Respons peserta tes diskor dengan tiga model partial credit (PCM I; II; dan III dan secara dikotomus. Hasil penskoran dianalisis dengan program Quest untuk mendapat-kan estimasi tingkat kesukaran butir (δ dan estimasi ke-mampuan peserta (θ untuk menentukan nilai fungsi informasi tes dan kesalahan baku estimasi. Penelitian simulasi mengguna-kan data bangkitan berdasarkan parameter empiris (δ dan θ memakai program statistik SAS dan akurasi estimasinya di-analisis dengan metode root mean squared error (RMSE. Hasil penelitian ini menunjukkan: (i Penskoran PCM dengan pem-bobotan mampu mengestimasi kemampuan lebih akurat di-bandingkan tanpa pembobotan maupun secara dikotomus; (ii Semakin banyak jumlah kategori dalam penskoran partial credit, semakin akurat. Kata kunci: model penskoran partial credit, butir multiple true-false ____________________________________________________________ THE PARTIAL CREDIT SCORING MODEL FOR THE MULTIPLE TRUE-FALSE BUTIRS IN PHYSICS Abstract This study is an attempt to overcome the weaknesses. This study aims to produce a polytomous scoring model for responses to multiple true-false butirs in order to get a more accurate estimation of abilities in physics. It adopts the Four-D model and its accuracy is assessed through empirical and simulation studies. The empirical study employed 15 multiple true-false butirs taken from the New Students Entrance Test of State University the year of 1996–2006. It administered to 410 new students enrolled

Some Considerations on the Partial Credit Model

Directory of Open Access Journals (Sweden)

H.H.F.M. Verstralen

2008-01-01

Full Text Available The Partial Credit Model (PCM is sometimes interpreted as a model for stepwise solution of polytomously scored items, where the item parameters are interpreted as di culties of the steps. It is argued that this interpretation is not justi ed. A model for stepwise solution is discussed. It is shown that the PCM is suited to model sums of binary responses which are not supposed to be stochastically independent. As a practical result, a statistical test of stochastic independence in the Rasch model is derived

Exploring the interplay of resilience and energy consumption for a task-based partial differential equations preconditioner

KAUST Repository

Rizzi, F.

2017-05-25

We discuss algorithm-based resilience to silent data corruptions (SDCs) in a task-based domain-decomposition preconditioner for partial differential equations (PDEs). The algorithm exploits a reformulation of the PDE as a sampling problem, followed by a solution update through data manipulation that is resilient to SDCs. The implementation is based on a server-client model where all state information is held by the servers, while clients are designed solely as computational units. Scalability tests run up to ∼ 51K cores show a parallel efficiency greater than 90%. We use a 2D elliptic PDE and a fault model based on random single and double bit-flip to demonstrate the resilience of the application to synthetically injected SDC. We discuss two fault scenarios: one based on the corruption of all data of a target task, and the other involving the corruption of a single data point. We show that for our application, given the test problem considered, a four-fold increase in the number of faults only yields a 2% change in the overhead to overcome their presence, from 7% to 9%. We then discuss potential savings in energy consumption via dynamic voltage/frequency scaling, and its interplay with fault-rates, and application overhead.

Exploring the interplay of resilience and energy consumption for a task-based partial differential equations preconditioner

KAUST Repository

Rizzi, F.; Morris, K.; Sargsyan, K.; Mycek, P.; Safta, C.; Le Maî tre, O.; Knio, Omar; Debusschere, B.J.

2017-01-01

We discuss algorithm-based resilience to silent data corruptions (SDCs) in a task-based domain-decomposition preconditioner for partial differential equations (PDEs). The algorithm exploits a reformulation of the PDE as a sampling problem, followed by a solution update through data manipulation that is resilient to SDCs. The implementation is based on a server-client model where all state information is held by the servers, while clients are designed solely as computational units. Scalability tests run up to ∼ 51K cores show a parallel efficiency greater than 90%. We use a 2D elliptic PDE and a fault model based on random single and double bit-flip to demonstrate the resilience of the application to synthetically injected SDC. We discuss two fault scenarios: one based on the corruption of all data of a target task, and the other involving the corruption of a single data point. We show that for our application, given the test problem considered, a four-fold increase in the number of faults only yields a 2% change in the overhead to overcome their presence, from 7% to 9%. We then discuss potential savings in energy consumption via dynamic voltage/frequency scaling, and its interplay with fault-rates, and application overhead.

Quantifying uncertainty, variability and likelihood for ordinary differential equation models

LENUS (Irish Health Repository)

Weisse, Andrea Y

2010-10-28

Abstract Background In many applications, ordinary differential equation (ODE) models are subject to uncertainty or variability in initial conditions and parameters. Both, uncertainty and variability can be quantified in terms of a probability density function on the state and parameter space. Results The partial differential equation that describes the evolution of this probability density function has a form that is particularly amenable to application of the well-known method of characteristics. The value of the density at some point in time is directly accessible by the solution of the original ODE extended by a single extra dimension (for the value of the density). This leads to simple methods for studying uncertainty, variability and likelihood, with significant advantages over more traditional Monte Carlo and related approaches especially when studying regions with low probability. Conclusions While such approaches based on the method of characteristics are common practice in other disciplines, their advantages for the study of biological systems have so far remained unrecognized. Several examples illustrate performance and accuracy of the approach and its limitations.

Reactor modeling and process analysis for partial oxidation of natural gas

NARCIS (Netherlands)

Albrecht, B.A.

2004-01-01

This thesis analyses a novel process of partial oxidation of natural gas and develops a numerical tool for the partial oxidation reactor modeling. The proposed process generates syngas in an integrated plant of a partial oxidation reactor, a syngas turbine and an air separation unit. This is called

Phase behavior and reactive transport of partial melt in heterogeneous mantle model

Science.gov (United States)

Jordan, J.; Hesse, M. A.

2013-12-01

The reactive transport of partial melt is the key process that leads to the chemical and physical differentiation of terrestrial planets and smaller celestial bodies. The essential role of the lithological heterogeneities during partial melting of the mantle is increasingly recognized. How far can enriched melts propagate while interacting with the ambient mantle? Can the melt flow emanating from a fertile heterogeneity be localized through a reactive infiltration feedback in a model without exogenous factors or contrived initial conditions? A full understanding of the role of heterogeneities requires reactive melt transport models that account for the phase behavior of major elements. Previous work on reactive transport in the mantle focuses on trace element partitioning; we present the first nonlinear chromatographic analysis of reactive melt transport in systems with binary solid solution. Our analysis shows that reactive melt transport in systems with binary solid solution leads to the formation of two separate reaction fronts: a slow melting/freezing front along which enthalpy change is dominant and a fast dissolution/precipitation front along which compositional changes are dominated by an ion-exchange process over enthalpy change. An intermediate state forms between these two fronts with a bulk-rock composition and enthalpy that are not necessarily bounded by the bulk-rock composition and enthalpy of either the enriched heterogeneity or the depleted ambient mantle. The formation of this intermediate state makes it difficult to anticipate the porosity changes and hence the stability of reaction fronts. Therefore, we develop a graphical representation for the solution that allows identification of the intermediate state by inspection, for all possible bulk-rock compositions and enthalpies of the heterogeneity and the ambient mantle. We apply the analysis to the partial melting of an enriched heterogeneity. This leads to the formation of moving precipitation

New Generalized Hyperbolic Functions to Find New Exact Solutions of the Nonlinear Partial Differential Equations

Directory of Open Access Journals (Sweden)

Yusuf Pandir

2013-01-01

Full Text Available We firstly give some new functions called generalized hyperbolic functions. By the using of the generalized hyperbolic functions, new kinds of transformations are defined to discover the exact approximate solutions of nonlinear partial differential equations. Based on the generalized hyperbolic function transformation of the generalized KdV equation and the coupled equal width wave equations (CEWE, we find new exact solutions of two equations and analyze the properties of them by taking different parameter values of the generalized hyperbolic functions. We think that these solutions are very important to explain some physical phenomena.

Mathematical Methods for Engineers and Scientists 3 Fourier Analysis, Partial Differential Equations and Variational Methods

CERN Document Server

Tang, Kwong-Tin

2007-01-01

Pedagogical insights gained through 30 years of teaching applied mathematics led the author to write this set of student oriented books. Topics such as complex analysis, matrix theory, vector and tensor analysis, Fourier analysis, integral transforms, ordinary and partial differential equations are presented in a discursive style that is readable and easy to follow. Numerous clearly stated, completely worked out examples together with carefully selected problem sets with answers are used to enhance students' understanding and manipulative skill. The goal is to make students comfortable and confident in using advanced mathematical tools in junior, senior, and beginning graduate courses.

Latent Partially Ordered Classification Models and Normal Mixtures

Science.gov (United States)

Tatsuoka, Curtis; Varadi, Ferenc; Jaeger, Judith

2013-01-01

Latent partially ordered sets (posets) can be employed in modeling cognitive functioning, such as in the analysis of neuropsychological (NP) and educational test data. Posets are cognitively diagnostic in the sense that classification states in these models are associated with detailed profiles of cognitive functioning. These profiles allow for…

The orthogonal gradients method: A radial basis functions method for solving partial differential equations on arbitrary surfaces

KAUST Repository

Piret, Cécile

2012-05-01

Much work has been done on reconstructing arbitrary surfaces using the radial basis function (RBF) method, but one can hardly find any work done on the use of RBFs to solve partial differential equations (PDEs) on arbitrary surfaces. In this paper, we investigate methods to solve PDEs on arbitrary stationary surfaces embedded in . R3 using the RBF method. We present three RBF-based methods that easily discretize surface differential operators. We take advantage of the meshfree character of RBFs, which give us a high accuracy and the flexibility to represent the most complex geometries in any dimension. Two out of the three methods, which we call the orthogonal gradients (OGr) methods are the result of our work and are hereby presented for the first time. © 2012 Elsevier Inc.

Nonlinear self-adjointness, conservation laws, and the construction of solutions of partial differential equations using conservation laws

International Nuclear Information System (INIS)

Ibragimov, N Kh; Avdonina, E D

2013-01-01

The method of nonlinear self-adjointness, which was recently developed by the first author, gives a generalization of Noether's theorem. This new method significantly extends approaches to constructing conservation laws associated with symmetries, since it does not require the existence of a Lagrangian. In particular, it can be applied to any linear equations and any nonlinear equations that possess at least one local conservation law. The present paper provides a brief survey of results on conservation laws which have been obtained by this method and published mostly in recent preprints of the authors, along with a method for constructing exact solutions of systems of partial differential equations with the use of conservation laws. In most cases the solutions obtained by the method of conservation laws cannot be found as invariant or partially invariant solutions. Bibliography: 23 titles

Melt migration modeling in partially molten upper mantle

Science.gov (United States)

Ghods, Abdolreza

The objective of this thesis is to investigate the importance of melt migration in shaping major characteristics of geological features associated with the partial melting of the upper mantle, such as sea-floor spreading, continental flood basalts and rifting. The partial melting produces permeable partially molten rocks and a buoyant low viscosity melt. Melt migrates through the partially molten rocks, and transfers mass and heat. Due to its much faster velocity and appreciable buoyancy, melt migration has the potential to modify dynamics of the upwelling partially molten plumes. I develop a 2-D, two-phase flow model and apply it to investigate effects of melt migration on the dynamics and melt generation of upwelling mantle plumes and focusing of melt migration beneath mid-ocean ridges. Melt migration changes distribution of the melt-retention buoyancy force and therefore affects the dynamics of the upwelling plume. This is investigated by modeling a plume with a constant initial melt of 10% where no further melting is considered. Melt migration polarizes melt-retention buoyancy force into high and low melt fraction regions at the top and bottom portions of the plume and therefore results in formation of a more slender and faster upwelling plume. Allowing the plume to melt as it ascends through the upper mantle also produces a slender and faster plume. It is shown that melt produced by decompressional melting of the plume migrates to the upper horizons of the plume, increases the upwelling velocity and thus, the volume of melt generated by the plume. Melt migration produces a plume which lacks the mushroom shape observed for the plume models without melt migration. Melt migration forms a high melt fraction layer beneath the sloping base of the impermeable oceanic lithosphere. Using realistic conditions of melting, freezing and melt extraction, I examine whether the high melt fraction layer is able to focus melt from a wide partial melting zone to a narrow region

Fitting and Testing Conditional Multinormal Partial Credit Models

Science.gov (United States)

Hessen, David J.

2012-01-01

A multinormal partial credit model for factor analysis of polytomously scored items with ordered response categories is derived using an extension of the Dutch Identity (Holland in "Psychometrika" 55:5-18, 1990). In the model, latent variables are assumed to have a multivariate normal distribution conditional on unweighted sums of item…

High-order asynchrony-tolerant finite difference schemes for partial differential equations

Science.gov (United States)

Aditya, Konduri; Donzis, Diego A.

2017-12-01

Synchronizations of processing elements (PEs) in massively parallel simulations, which arise due to communication or load imbalances between PEs, significantly affect the scalability of scientific applications. We have recently proposed a method based on finite-difference schemes to solve partial differential equations in an asynchronous fashion - synchronization between PEs is relaxed at a mathematical level. While standard schemes can maintain their stability in the presence of asynchrony, their accuracy is drastically affected. In this work, we present a general methodology to derive asynchrony-tolerant (AT) finite difference schemes of arbitrary order of accuracy, which can maintain their accuracy when synchronizations are relaxed. We show that there are several choices available in selecting a stencil to derive these schemes and discuss their effect on numerical and computational performance. We provide a simple classification of schemes based on the stencil and derive schemes that are representative of different classes. Their numerical error is rigorously analyzed within a statistical framework to obtain the overall accuracy of the solution. Results from numerical experiments are used to validate the performance of the schemes.

Error propagation of partial least squares for parameters optimization in NIR modeling

Science.gov (United States)

Du, Chenzhao; Dai, Shengyun; Qiao, Yanjiang; Wu, Zhisheng

2018-03-01

A novel methodology is proposed to determine the error propagation of partial least-square (PLS) for parameters optimization in near-infrared (NIR) modeling. The parameters include spectral pretreatment, latent variables and variable selection. In this paper, an open source dataset (corn) and a complicated dataset (Gardenia) were used to establish PLS models under different modeling parameters. And error propagation of modeling parameters for water quantity in corn and geniposide quantity in Gardenia were presented by both type І and type II error. For example, when variable importance in the projection (VIP), interval partial least square (iPLS) and backward interval partial least square (BiPLS) variable selection algorithms were used for geniposide in Gardenia, compared with synergy interval partial least squares (SiPLS), the error weight varied from 5% to 65%, 55% and 15%. The results demonstrated how and what extent the different modeling parameters affect error propagation of PLS for parameters optimization in NIR modeling. The larger the error weight, the worse the model. Finally, our trials finished a powerful process in developing robust PLS models for corn and Gardenia under the optimal modeling parameters. Furthermore, it could provide a significant guidance for the selection of modeling parameters of other multivariate calibration models.

Error propagation of partial least squares for parameters optimization in NIR modeling.

Science.gov (United States)

Du, Chenzhao; Dai, Shengyun; Qiao, Yanjiang; Wu, Zhisheng

2018-03-05

A novel methodology is proposed to determine the error propagation of partial least-square (PLS) for parameters optimization in near-infrared (NIR) modeling. The parameters include spectral pretreatment, latent variables and variable selection. In this paper, an open source dataset (corn) and a complicated dataset (Gardenia) were used to establish PLS models under different modeling parameters. And error propagation of modeling parameters for water quantity in corn and geniposide quantity in Gardenia were presented by both type І and type II error. For example, when variable importance in the projection (VIP), interval partial least square (iPLS) and backward interval partial least square (BiPLS) variable selection algorithms were used for geniposide in Gardenia, compared with synergy interval partial least squares (SiPLS), the error weight varied from 5% to 65%, 55% and 15%. The results demonstrated how and what extent the different modeling parameters affect error propagation of PLS for parameters optimization in NIR modeling. The larger the error weight, the worse the model. Finally, our trials finished a powerful process in developing robust PLS models for corn and Gardenia under the optimal modeling parameters. Furthermore, it could provide a significant guidance for the selection of modeling parameters of other multivariate calibration models. Copyright © 2017. Published by Elsevier B.V.

An analytic, approximate method for modeling steady, three-dimensional flow to partially penetrating wells

Science.gov (United States)

Bakker, Mark

2001-05-01

An analytic, approximate solution is derived for the modeling of three-dimensional flow to partially penetrating wells. The solution is written in terms of a correction on the solution for a fully penetrating well and is obtained by dividing the aquifer up, locally, in a number of aquifer layers. The resulting system of differential equations is solved by application of the theory for multiaquifer flow. The presented approach has three major benefits. First, the solution may be applied to any groundwater model that can simulate flow to a fully penetrating well; the solution may be superimposed onto the solution for the fully penetrating well to simulate the local three-dimensional drawdown and flow field. Second, the approach is applicable to isotropic, anisotropic, and stratified aquifers and to both confined and unconfined flow. Third, the solution extends over a small area around the well only; outside this area the three-dimensional effect of the partially penetrating well is negligible, and no correction to the fully penetrating well is needed. A number of comparisons are made to existing three-dimensional, analytic solutions, including radial confined and unconfined flow and a well in a uniform flow field. It is shown that a subdivision in three layers is accurate for many practical cases; very accurate solutions are obtained with more layers.

The strategic value of partial vertical integration

OpenAIRE

Fiocco, Raffaele

2014-01-01

We investigate the strategic incentives for partial vertical integration, namely, partial ownership agreements between manufacturers and retailers, when retailers privately know their costs and engage in differentiated good price competition. The partial misalignment between the profit objectives within a partially integrated manufacturer-retailer hierarchy entails a higher retail price than under full integration. This `information vertical effect' translates into an opposite ...

A fast direct method for block triangular Toeplitz-like with tri-diagonal block systems from time-fractional partial differential equations

Science.gov (United States)

Ke, Rihuan; Ng, Michael K.; Sun, Hai-Wei

2015-12-01

In this paper, we study the block lower triangular Toeplitz-like with tri-diagonal blocks system which arises from the time-fractional partial differential equation. Existing fast numerical solver (e.g., fast approximate inversion method) cannot handle such linear system as the main diagonal blocks are different. The main contribution of this paper is to propose a fast direct method for solving this linear system, and to illustrate that the proposed method is much faster than the classical block forward substitution method for solving this linear system. Our idea is based on the divide-and-conquer strategy and together with the fast Fourier transforms for calculating Toeplitz matrix-vector multiplication. The complexity needs O (MNlog2 ⁡ M) arithmetic operations, where M is the number of blocks (the number of time steps) in the system and N is the size (number of spatial grid points) of each block. Numerical examples from the finite difference discretization of time-fractional partial differential equations are also given to demonstrate the efficiency of the proposed method.

Substrate stiffness and oxygen as regulators of stem cell differentiation during skeletal tissue regeneration: a mechanobiological model.

Directory of Open Access Journals (Sweden)

Darren Paul Burke

Full Text Available Extrinsic mechanical signals have been implicated as key regulators of mesenchymal stem cell (MSC differentiation. It has been possible to test different hypotheses for mechano-regulated MSC differentiation by attempting to simulate regenerative events such as bone fracture repair, where repeatable spatial and temporal patterns of tissue differentiation occur. More recently, in vitro studies have identified other environmental cues such as substrate stiffness and oxygen tension as key regulators of MSC differentiation; however it remains unclear if and how such cues determine stem cell fate in vivo. As part of this study, a computational model was developed to test the hypothesis that substrate stiffness and oxygen tension regulate stem cell differentiation during fracture healing. Rather than assuming mechanical signals act directly on stem cells to determine their differentiation pathway, it is postulated that they act indirectly to regulate angiogenesis and hence partially determine the local oxygen environment within a regenerating tissue. Chondrogenesis of MSCs was hypothesized to occur in low oxygen regions, while in well vascularised regions of the regenerating tissue a soft local substrate was hypothesised to facilitate adipogenesis while a stiff substrate facilitated osteogenesis. Predictions from the model were compared to both experimental data and to predictions of a well established computational mechanobiological model where tissue differentiation is assumed to be regulated directly by the local mechanical environment. The model predicted all the major events of fracture repair, including cartilaginous bridging, endosteal and periosteal bony bridging and bone remodelling. It therefore provides support for the hypothesis that substrate stiffness and oxygen play a key role in regulating MSC fate during regenerative events such as fracture healing.

Sparse grid spectral methods for the numerical solution of partial differential equations with periodic boundary conditions

International Nuclear Information System (INIS)

Kupka, F.

1997-11-01

This thesis deals with the extension of sparse grid techniques to spectral methods for the solution of partial differential equations with periodic boundary conditions. A review on boundary and initial-boundary value problems and a discussion on numerical resolution is used to motivate this research. Spectral methods are introduced by projection techniques, and by three model problems: the stationary and the transient Helmholtz equations, and the linear advection equation. The approximation theory on the hyperbolic cross is reviewed and its close relation to sparse grids is demonstrated. This approach extends to non-periodic problems. Various Sobolev spaces with dominant mixed derivative are introduced to provide error estimates for Fourier approximation and interpolation on the hyperbolic cross and on sparse grids by means of Sobolev norms. The theorems are immediately applicable to the stability and convergence analysis of sparse grid spectral methods. This is explicitly demonstrated for the three model problems. A variant of the von Neumann condition is introduced to simplify the stability analysis of the time-dependent model problems. The discrete Fourier transformation on sparse grids is discussed together with its software implementation. Results on numerical experiments are used to illustrate the performance of the new method with respect to the smoothness properties of each example. The potential of the method in mathematical modelling is estimated and generalizations to other sparse grid methods are suggested. The appendix includes a complete Fortran90 program to solve the linear advection equation by the sparse grid Fourier collocation method and a third-order Runge-Kutta routine for integration in time. (author)

Exponentially Convergent Algorithms for Abstract Differential Equations

CERN Document Server

Gavrilyuk, Ivan; Vasylyk, Vitalii

2011-01-01

This book presents new accurate and efficient exponentially convergent methods for abstract differential equations with unbounded operator coefficients in Banach space. These methods are highly relevant for the practical scientific computing since the equations under consideration can be seen as the meta-models of systems of ordinary differential equations (ODE) as well as the partial differential equations (PDEs) describing various applied problems. The framework of functional analysis allows one to obtain very general but at the same time transparent algorithms and mathematical results which

Differential-difference model for textile engineering

International Nuclear Information System (INIS)

Wu Guocheng; Zhao Ling; He Jihuan

2009-01-01

Woven fabric is manifestly not a continuum and therefore Darcy's law or its modifications, or any other differential models are invalid theoretically. A differential-difference model for air transport in discontinuous media is introduced using conservation of mass, conservation of energy, and the equation of state in discrete space and continuous time, capillary pressure is obtained by dimensional analysis.

A Table Lookup Method for Exact Analytical Solutions of Nonlinear Fractional Partial Differential Equations

Directory of Open Access Journals (Sweden)

Ji Juan-Juan

2017-01-01

Full Text Available A table lookup method for solving nonlinear fractional partial differential equations (fPDEs is proposed in this paper. Looking up the corresponding tables, we can quickly obtain the exact analytical solutions of fPDEs by using this method. To illustrate the validity of the method, we apply it to construct the exact analytical solutions of four nonlinear fPDEs, namely, the time fractional simplified MCH equation, the space-time fractional combined KdV-mKdV equation, the (2+1-dimensional time fractional Zoomeron equation, and the space-time fractional ZKBBM equation. As a result, many new types of exact analytical solutions are obtained including triangular periodic solution, hyperbolic function solution, singular solution, multiple solitary wave solution, and Jacobi elliptic function solution.

Mathematical Modelling of Intraretinal Oxygen Partial Pressure

African Journals Online (AJOL)

Erah

The system of non-linear differential equations was solved numerically using Runge-kutta. Nystroms method. ... artery occlusion. Keywords: Mathematical modeling, Intraretinal oxygen pressure, Retinal capillaries, Oxygen ..... Mass transfer,.

Numerical Solution of the Fractional Partial Differential Equations by the Two-Dimensional Fractional-Order Legendre Functions

Directory of Open Access Journals (Sweden)

Fukang Yin

2013-01-01

Full Text Available A numerical method is presented to obtain the approximate solutions of the fractional partial differential equations (FPDEs. The basic idea of this method is to achieve the approximate solutions in a generalized expansion form of two-dimensional fractional-order Legendre functions (2D-FLFs. The operational matrices of integration and derivative for 2D-FLFs are first derived. Then, by these matrices, a system of algebraic equations is obtained from FPDEs. Hence, by solving this system, the unknown 2D-FLFs coefficients can be computed. Three examples are discussed to demonstrate the validity and applicability of the proposed method.

Automatic differentiation algorithms in model analysis

NARCIS (Netherlands)

Huiskes, M.J.

2002-01-01

Title: Automatic differentiation algorithms in model analysis
Author: M.J. Huiskes
Date: 19 March, 2002

In this thesis automatic differentiation algorithms and derivative-based methods

Coding with partially hidden Markov models

DEFF Research Database (Denmark)

Forchhammer, Søren; Rissanen, J.

1995-01-01

Partially hidden Markov models (PHMM) are introduced. They are a variation of the hidden Markov models (HMM) combining the power of explicit conditioning on past observations and the power of using hidden states. (P)HMM may be combined with arithmetic coding for lossless data compression. A general...... 2-part coding scheme for given model order but unknown parameters based on PHMM is presented. A forward-backward reestimation of parameters with a redefined backward variable is given for these models and used for estimating the unknown parameters. Proof of convergence of this reestimation is given....... The PHMM structure and the conditions of the convergence proof allows for application of the PHMM to image coding. Relations between the PHMM and hidden Markov models (HMM) are treated. Results of coding bi-level images with the PHMM coding scheme is given. The results indicate that the PHMM can adapt...

On the removal of boundary errors caused by Runge-Kutta integration of non-linear partial differential equations

Science.gov (United States)

Abarbanel, Saul; Gottlieb, David; Carpenter, Mark H.

1994-01-01

It has been previously shown that the temporal integration of hyperbolic partial differential equations (PDE's) may, because of boundary conditions, lead to deterioration of accuracy of the solution. A procedure for removal of this error in the linear case has been established previously. In the present paper we consider hyperbolic (PDE's) (linear and non-linear) whose boundary treatment is done via the SAT-procedure. A methodology is present for recovery of the full order of accuracy, and has been applied to the case of a 4th order explicit finite difference scheme.

Mathematical Modelling of Continuous Biotechnological Processes

Science.gov (United States)

Pencheva, T.; Hristozov, I.; Shannon, A. G.

2003-01-01

Biotechnological processes (BTP) are characterized by a complicated structure of organization and interdependent characteristics. Partial differential equations or systems of partial differential equations are used for their behavioural description as objects with distributed parameters. Modelling of substrate without regard to dispersion…

3D early embryogenesis image filtering by nonlinear partial differential equations.

Science.gov (United States)

Krivá, Z; Mikula, K; Peyriéras, N; Rizzi, B; Sarti, A; Stasová, O

2010-08-01

We present nonlinear diffusion equations, numerical schemes to solve them and their application for filtering 3D images obtained from laser scanning microscopy (LSM) of living zebrafish embryos, with a goal to identify the optimal filtering method and its parameters. In the large scale applications dealing with analysis of 3D+time embryogenesis images, an important objective is a correct detection of the number and position of cell nuclei yielding the spatio-temporal cell lineage tree of embryogenesis. The filtering is the first and necessary step of the image analysis chain and must lead to correct results, removing the noise, sharpening the nuclei edges and correcting the acquisition errors related to spuriously connected subregions. In this paper we study such properties for the regularized Perona-Malik model and for the generalized mean curvature flow equations in the level-set formulation. A comparison with other nonlinear diffusion filters, like tensor anisotropic diffusion and Beltrami flow, is also included. All numerical schemes are based on the same discretization principles, i.e. finite volume method in space and semi-implicit scheme in time, for solving nonlinear partial differential equations. These numerical schemes are unconditionally stable, fast and naturally parallelizable. The filtering results are evaluated and compared first using the Mean Hausdorff distance between a gold standard and different isosurfaces of original and filtered data. Then, the number of isosurface connected components in a region of interest (ROI) detected in original and after the filtering is compared with the corresponding correct number of nuclei in the gold standard. Such analysis proves the robustness and reliability of the edge preserving nonlinear diffusion filtering for this type of data and lead to finding the optimal filtering parameters for the studied models and numerical schemes. Further comparisons consist in ability of splitting the very close objects which

Basis adaptation and domain decomposition for steady-state partial differential equations with random coefficients

Energy Technology Data Exchange (ETDEWEB)

Tipireddy, R.; Stinis, P.; Tartakovsky, A. M.

2017-12-01

We present a novel approach for solving steady-state stochastic partial differential equations (PDEs) with high-dimensional random parameter space. The proposed approach combines spatial domain decomposition with basis adaptation for each subdomain. The basis adaptation is used to address the curse of dimensionality by constructing an accurate low-dimensional representation of the stochastic PDE solution (probability density function and/or its leading statistical moments) in each subdomain. Restricting the basis adaptation to a specific subdomain affords finding a locally accurate solution. Then, the solutions from all of the subdomains are stitched together to provide a global solution. We support our construction with numerical experiments for a steady-state diffusion equation with a random spatially dependent coefficient. Our results show that highly accurate global solutions can be obtained with significantly reduced computational costs.

Final Report: Symposium on Adaptive Methods for Partial Differential Equations

Energy Technology Data Exchange (ETDEWEB)

Pernice, M.; Johnson, C.R.; Smith, P.J.; Fogelson, A.

1998-12-10

OAK-B135 Final Report: Symposium on Adaptive Methods for Partial Differential Equations. Complex physical phenomena often include features that span a wide range of spatial and temporal scales. Accurate simulation of such phenomena can be difficult to obtain, and computations that are under-resolved can even exhibit spurious features. While it is possible to resolve small scale features by increasing the number of grid points, global grid refinement can quickly lead to problems that are intractable, even on the largest available computing facilities. These constraints are particularly severe for three dimensional problems that involve complex physics. One way to achieve the needed resolution is to refine the computational mesh locally, in only those regions where enhanced resolution is required. Adaptive solution methods concentrate computational effort in regions where it is most needed. These methods have been successfully applied to a wide variety of problems in computational science and engineering. Adaptive methods can be difficult to implement, prompting the development of tools and environments to facilitate their use. To ensure that the results of their efforts are useful, algorithm and tool developers must maintain close communication with application specialists. Conversely it remains difficult for application specialists who are unfamiliar with the methods to evaluate the trade-offs between the benefits of enhanced local resolution and the effort needed to implement an adaptive solution method.

Effect of partial heating at mid of vertical plate adjacent to porous medium

Science.gov (United States)

Mulla, Mohammed Fahimuddin; Pallan, Khalid. M.; Al-Rashed, A. A. A. A.

2018-05-01

Heat and mass transfer in porous medium due to heating of vertical plate at mid-section is analyzed for various physical parameters. The heat and mass transfer in porous medium is modeled with the help of momentum, energy and concentration equations in terms of non-dimensional partial differential equations. The partial differential equations are converted into simpler form of algebraic equations with the help of finite element method. A computer code is developed to assemble the matrix form of algebraic equations into global matrices and then to solve them in an iterative manner to obtain the temperature, concentration and streamline distribution inside the porous medium. It is found that the heat transfer behavior of porous medium heated at middle section is considerably different from other cases.

Hopf bifurcation in a partial dependent predator-prey system with delay

International Nuclear Information System (INIS)

Zhao Huitao; Lin Yiping

2009-01-01

In this paper, a partial dependent predator-prey model with time delay is studied by using the theory of functional differential equation and Hassard's method, the condition on which positive equilibrium exists and Hopf bifurcation occurs are given. Finally, numerical simulations are performed to support the analytical results, and the chaotic behaviors are observed.

Introduction to computation and modeling for differential equations

CERN Document Server

Edsberg, Lennart

2008-01-01

An introduction to scientific computing for differential equationsIntroduction to Computation and Modeling for Differential Equations provides a unified and integrated view of numerical analysis, mathematical modeling in applications, and programming to solve differential equations, which is essential in problem-solving across many disciplines, such as engineering, physics, and economics. This book successfully introduces readers to the subject through a unique ""Five-M"" approach: Modeling, Mathematics, Methods, MATLAB, and Multiphysics. This approach facilitates a thorough understanding of h

OUTLIER DETECTION IN PARTIAL ERRORS-IN-VARIABLES MODEL

Directory of Open Access Journals (Sweden)

JUN ZHAO

Full Text Available The weighed total least square (WTLS estimate is very sensitive to the outliers in the partial EIV model. A new procedure for detecting outliers based on the data-snooping is presented in this paper. Firstly, a two-step iterated method of computing the WTLS estimates for the partial EIV model based on the standard LS theory is proposed. Secondly, the corresponding w-test statistics are constructed to detect outliers while the observations and coefficient matrix are contaminated with outliers, and a specific algorithm for detecting outliers is suggested. When the variance factor is unknown, it may be estimated by the least median squares (LMS method. At last, the simulated data and real data about two-dimensional affine transformation are analyzed. The numerical results show that the new test procedure is able to judge that the outliers locate in x component, y component or both components in coordinates while the observations and coefficient matrix are contaminated with outliers

Differential equation models for sharp threshold dynamics.

Science.gov (United States)

Schramm, Harrison C; Dimitrov, Nedialko B

2014-01-01

We develop an extension to differential equation models of dynamical systems to allow us to analyze probabilistic threshold dynamics that fundamentally and globally change system behavior. We apply our novel modeling approach to two cases of interest: a model of infectious disease modified for malware where a detection event drastically changes dynamics by introducing a new class in competition with the original infection; and the Lanchester model of armed conflict, where the loss of a key capability drastically changes the effectiveness of one of the sides. We derive and demonstrate a step-by-step, repeatable method for applying our novel modeling approach to an arbitrary system, and we compare the resulting differential equations to simulations of the system's random progression. Our work leads to a simple and easily implemented method for analyzing probabilistic threshold dynamics using differential equations. Published by Elsevier Inc.

Understanding Rasch Measurement: Partial Credit Model and Pivot Anchoring.

Science.gov (United States)

Bode, Rita K.

2001-01-01

Describes the Rasch measurement partial credit model, what it is, how it differs from other Rasch models, and when and how to use it. Also describes the calibration of instruments with increasingly complex items. Explains pivot anchoring and illustrates its use and describes the effect of pivot anchoring on step calibrations, item hierarchy, and…

An Introduction to the Partial Credit Model for Developing Nursing Assessments.

Science.gov (United States)

Fox, Christine

1999-01-01

Demonstrates how the partial credit model, a variation of the Rasch Measurement Model, can be used to develop performance-based assessments for nursing education. Applies the model using the Practical Knowledge Inventory for Nurses. (SK)

A Sequential, Implicit, Wavelet-Based Solver for Multi-Scale Time-Dependent Partial Differential Equations

Directory of Open Access Journals (Sweden)

Donald A. McLaren

2013-04-01

Full Text Available This paper describes and tests a wavelet-based implicit numerical method for solving partial differential equations. Intended for problems with localized small-scale interactions, the method exploits the form of the wavelet decomposition to divide the implicit system created by the time-discretization into multiple smaller systems that can be solved sequentially. Included is a test on a basic non-linear problem, with both the results of the test, and the time required to calculate them, compared with control results based on a single system with fine resolution. The method is then tested on a non-trivial problem, its computational time and accuracy checked against control results. In both tests, it was found that the method requires less computational expense than the control. Furthermore, the method showed convergence towards the fine resolution control results.

Optimal Designs for the Generalized Partial Credit Model

OpenAIRE

Bürkner, Paul-Christian; Schwabe, Rainer; Holling, Heinz

2018-01-01

Analyzing ordinal data becomes increasingly important in psychology, especially in the context of item response theory. The generalized partial credit model (GPCM) is probably the most widely used ordinal model and finds application in many large scale educational assessment studies such as PISA. In the present paper, optimal test designs are investigated for estimating persons' abilities with the GPCM for calibrated tests when item parameters are known from previous studies. We will derive t...

Controllability and Stabilization of Bilinear and Semilinear Partial Differential Equations

DEFF Research Database (Denmark)

Krishnaswamy, Vijayaraghavan

The topic of the thesis is the investigation of the question of controllability of weakly nonlinear partial differntial equations. The method is based on the Hilbert Uniqueness Method.......The topic of the thesis is the investigation of the question of controllability of weakly nonlinear partial differntial equations. The method is based on the Hilbert Uniqueness Method....

A Solution Space for a System of Null-State Partial Differential Equations: Part 1

Science.gov (United States)

Flores, Steven M.; Kleban, Peter

2015-01-01

This article is the first of four that completely and rigorously characterize a solution space for a homogeneous system of 2 N + 3 linear partial differential equations (PDEs) in 2 N variables that arises in conformal field theory (CFT) and multiple Schramm-Löwner evolution (SLE). In CFT, these are null-state equations and conformal Ward identities. They govern partition functions for the continuum limit of a statistical cluster or loop-gas model, such as percolation, or more generally the Potts models and O( n) models, at the statistical mechanical critical point. (SLE partition functions also satisfy these equations.) For such a lattice model in a polygon with its 2 N sides exhibiting a free/fixed side-alternating boundary condition , this partition function is proportional to the CFT correlation function where the w i are the vertices of and where is a one-leg corner operator. (Partition functions for "crossing events" in which clusters join the fixed sides of in some specified connectivity are linear combinations of such correlation functions.) When conformally mapped onto the upper half-plane, methods of CFT show that this correlation function satisfies the system of PDEs that we consider. In this first article, we use methods of analysis to prove that the dimension of this solution space is no more than C N , the Nth Catalan number. While our motivations are based in CFT, our proofs are completely rigorous. This proof is contained entirely within this article, except for the proof of Lemma 14, which constitutes the second article (Flores and Kleban, in Commun Math Phys, arXiv:1404.0035, 2014). In the third article (Flores and Kleban, in Commun Math Phys, arXiv:1303.7182, 2013), we use the results of this article to prove that the solution space of this system of PDEs has dimension C N and is spanned by solutions constructed with the CFT Coulomb gas (contour integral) formalism. In the fourth article (Flores and Kleban, in Commun Math Phys, arXiv:1405

Modeling and Prediction Using Stochastic Differential Equations

DEFF Research Database (Denmark)

Juhl, Rune; Møller, Jan Kloppenborg; Jørgensen, John Bagterp

2016-01-01

Pharmacokinetic/pharmakodynamic (PK/PD) modeling for a single subject is most often performed using nonlinear models based on deterministic ordinary differential equations (ODEs), and the variation between subjects in a population of subjects is described using a population (mixed effects) setup...... deterministic and can predict the future perfectly. A more realistic approach would be to allow for randomness in the model due to e.g., the model be too simple or errors in input. We describe a modeling and prediction setup which better reflects reality and suggests stochastic differential equations (SDEs...

Optimal Retail Price Model for Partial Consignment to Multiple Retailers

Directory of Open Access Journals (Sweden)

Po-Yu Chen

2017-01-01

Full Text Available This paper investigates the product pricing decision-making problem under a consignment stock policy in a two-level supply chain composed of one supplier and multiple retailers. The effects of the supplier’s wholesale prices and its partial inventory cost absorption of the retail prices of retailers with different market shares are investigated. In the partial product consignment model this paper proposes, the seller and the retailers each absorb part of the inventory costs. This model also provides general solutions for the complete product consignment and the traditional policy that adopts no product consignment. In other words, both the complete consignment and nonconsignment models are extensions of the proposed model (i.e., special cases. Research results indicated that the optimal retail price must be between 1/2 (50% and 2/3 (66.67% times the upper limit of the gross profit. This study also explored the results and influence of parameter variations on optimal retail price in the model.

Ionic diffusion through confined geometries: from Langevin equations to partial differential equations

International Nuclear Information System (INIS)

Nadler, Boaz; Schuss, Zeev; Singer, Amit; Eisenberg, R S

2004-01-01

Ionic diffusion through and near small domains is of considerable importance in molecular biophysics in applications such as permeation through protein channels and diffusion near the charged active sites of macromolecules. The motion of the ions in these settings depends on the specific nanoscale geometry and charge distribution in and near the domain, so standard continuum type approaches have obvious limitations. The standard machinery of equilibrium statistical mechanics includes microscopic details, but is also not applicable, because these systems are usually not in equilibrium due to concentration gradients and to the presence of an external applied potential, which drive a non-vanishing stationary current through the system. We present a stochastic molecular model for the diffusive motion of interacting particles in an external field of force and a derivation of effective partial differential equations and their boundary conditions that describe the stationary non-equilibrium system. The interactions can include electrostatic, Lennard-Jones and other pairwise forces. The analysis yields a new type of Poisson-Nernst-Planck equations, that involves conditional and unconditional charge densities and potentials. The conditional charge densities are the non-equilibrium analogues of the well studied pair correlation functions of equilibrium statistical physics. Our proposed theory is an extension of equilibrium statistical mechanics of simple fluids to stationary non-equilibrium problems. The proposed system of equations differs from the standard Poisson-Nernst-Planck system in two important aspects. First, the force term depends on conditional densities and thus on the finite size of ions, and second, it contains the dielectric boundary force on a discrete ion near dielectric interfaces. Recently, various authors have shown that both of these terms are important for diffusion through confined geometries in the context of ion channels

Fractal-Based Methods and Inverse Problems for Differential Equations: Current State of the Art

Directory of Open Access Journals (Sweden)

Herb E. Kunze

2014-01-01

Full Text Available We illustrate, in this short survey, the current state of the art of fractal-based techniques and their application to the solution of inverse problems for ordinary and partial differential equations. We review several methods based on the Collage Theorem and its extensions. We also discuss two innovative applications: the first one is related to a vibrating string model while the second one considers a collage-based approach for solving inverse problems for partial differential equations on a perforated domain.

DISPL-1, 2. Order Nonlinear Partial Differential Equation System Solution for Kinetics Diffusion Problems

International Nuclear Information System (INIS)

Leaf, G.K.; Minkoff, M.

1982-01-01

1 - Description of problem or function: DISPL1 is a software package for solving second-order nonlinear systems of partial differential equations including parabolic, elliptic, hyperbolic, and some mixed types. The package is designed primarily for chemical kinetics- diffusion problems, although not limited to these problems. Fairly general nonlinear boundary conditions are allowed as well as inter- face conditions for problems in an inhomogeneous medium. The spatial domain is one- or two-dimensional with rectangular Cartesian, cylindrical, or spherical (in one dimension only) geometry. 2 - Method of solution: The numerical method is based on the use of Galerkin's procedure combined with the use of B-Splines (C.W.R. de-Boor's B-spline package) to generate a system of ordinary differential equations. These equations are solved by a sophisticated ODE software package which is a modified version of Hindmarsh's GEAR package, NESC Abstract 592. 3 - Restrictions on the complexity of the problem: The spatial domain must be rectangular with sides parallel to the coordinate geometry. Cross derivative terms are not permitted in the PDE. The order of the B-Splines is at most 12. Other parameters such as the number of mesh points in each coordinate direction, the number of PDE's etc. are set in a macro table used by the MORTRAn2 preprocessor in generating the object code

Parameter estimation in stochastic differential equations

CERN Document Server

Bishwal, Jaya P N

2008-01-01

Parameter estimation in stochastic differential equations and stochastic partial differential equations is the science, art and technology of modelling complex phenomena and making beautiful decisions. The subject has attracted researchers from several areas of mathematics and other related fields like economics and finance. This volume presents the estimation of the unknown parameters in the corresponding continuous models based on continuous and discrete observations and examines extensively maximum likelihood, minimum contrast and Bayesian methods. Useful because of the current availability of high frequency data is the study of refined asymptotic properties of several estimators when the observation time length is large and the observation time interval is small. Also space time white noise driven models, useful for spatial data, and more sophisticated non-Markovian and non-semimartingale models like fractional diffusions that model the long memory phenomena are examined in this volume.

Computable Error Estimates for Finite Element Approximations of Elliptic Partial Differential Equations with Rough Stochastic Data

KAUST Repository

Hall, Eric Joseph

2016-12-08

We derive computable error estimates for finite element approximations of linear elliptic partial differential equations with rough stochastic coefficients. In this setting, the exact solutions contain high frequency content that standard a posteriori error estimates fail to capture. We propose goal-oriented estimates, based on local error indicators, for the pathwise Galerkin and expected quadrature errors committed in standard, continuous, piecewise linear finite element approximations. Derived using easily validated assumptions, these novel estimates can be computed at a relatively low cost and have applications to subsurface flow problems in geophysics where the conductivities are assumed to have lognormal distributions with low regularity. Our theory is supported by numerical experiments on test problems in one and two dimensions.

Solution of differential equations by application of transformation groups

Science.gov (United States)

Driskell, C. N., Jr.; Gallaher, L. J.; Martin, R. H., Jr.

1968-01-01

Report applies transformation groups to the solution of systems of ordinary differential equations and partial differential equations. Lies theorem finds an integrating factor for appropriate invariance group or groups can be found and can be extended to partial differential equations.

A partial hearing animal model for chronic electro-acoustic stimulation

Science.gov (United States)

Irving, S.; Wise, A. K.; Millard, R. E.; Shepherd, R. K.; Fallon, J. B.

2014-08-01

Objective. Cochlear implants (CIs) have provided some auditory function to hundreds of thousands of people around the world. Although traditionally carried out only in profoundly deaf patients, the eligibility criteria for implantation have recently been relaxed to include many partially-deaf patients with useful levels of hearing. These patients receive both electrical stimulation from their implant and acoustic stimulation via their residual hearing (electro-acoustic stimulation; EAS) and perform very well. It is unclear how EAS improves speech perception over electrical stimulation alone, and little evidence exists about the nature of the interactions between electric and acoustic stimuli. Furthermore, clinical results suggest that some patients that undergo cochlear implantation lose some, if not all, of their residual hearing, reducing the advantages of EAS over electrical stimulation alone. A reliable animal model with clinically-relevant partial deafness combined with clinical CIs is important to enable these issues to be studied. This paper outlines such a model that has been successfully used in our laboratory. Approach. This paper outlines a battery of techniques used in our laboratory to generate, validate and examine an animal model of partial deafness and chronic CI use. Main results. Ototoxic deafening produced bilaterally symmetrical hearing thresholds in neonatal and adult animals. Electrical activation of the auditory system was confirmed, and all animals were chronically stimulated via adapted clinical CIs. Acoustic compound action potentials (CAPs) were obtained from partially-hearing cochleae, using the CI amplifier. Immunohistochemical analysis allows the effects of deafness and electrical stimulation on cell survival to be studied. Significance. This animal model has applications in EAS research, including investigating the functional interactions between electric and acoustic stimulation, and the development of techniques to maintain residual

A discrete model of a modified Burgers' partial differential equation

Science.gov (United States)

Mickens, R. E.; Shoosmith, J. N.

1990-01-01

A new finite-difference scheme is constructed for a modified Burger's equation. Three special cases of the equation are considered, and the 'exact' difference schemes for the space- and time-independent forms of the equation are presented, along with the diffusion-free case of Burger's equation modeled by a difference equation. The desired difference scheme is then obtained by imposing on any difference model of the initial equation the requirement that, in the appropriate limits, its difference scheme must reduce the results of the obtained equations.

Pricing derivatives under Lévy models modern finite-difference and pseudo-differential operators approach

CERN Document Server

Itkin, Andrey

2017-01-01

This monograph presents a novel numerical approach to solving partial integro-differential equations arising in asset pricing models with jumps, which greatly exceeds the efficiency of existing approaches. The method, based on pseudo-differential operators and several original contributions to the theory of finite-difference schemes, is new as applied to the Lévy processes in finance, and is herein presented for the first time in a single volume. The results within, developed in a series of research papers, are collected and arranged together with the necessary background material from Lévy processes, the modern theory of finite-difference schemes, the theory of M-matrices and EM-matrices, etc., thus forming a self-contained work that gives the reader a smooth introduction to the subject. For readers with no knowledge of finance, a short explanation of the main financial terms and notions used in the book is given in the glossary. The latter part of the book demonstrates the efficacy of the method by solvin...

Normal differential renal function does not indicate a normal kidney after partial ureteropelvic obstruction and subsequent relief in 2-week-old piglets

Energy Technology Data Exchange (ETDEWEB)

Dissing, Thomas H.; Mikkelsen, Mette Marie; Pedersen, Michael; Froekiaer, Joergen; Djurhuus, Jens Christian [University of Aarhus, Institute of Clinical Medicine, Aarhus (Denmark); Eskild-Jensen, Anni [Aarhus University Hospital, Department of Nuclear Medicine, Aarhus Sygehus, Aarhus (Denmark); Gordon, Isky [University College London, Institute of Child Health, London (United Kingdom); University College London, Radiology and Physics Unit, Institute of Child Health, London (United Kingdom)

2008-09-15

We investigated the functional consequences of relieving ureteric obstruction in young pigs with experimental hydronephrosis (HN) induced by partial unilateral ureteropelvic obstruction. Three groups of animals were followed from the age of 2 weeks to the age of 14 weeks: Eight animals had severe or grades 3-4 HN throughout the study. Six animals had relief of the obstruction after 4 weeks. Six animals received sham operations at both ages. Morphological and functional examinations were performed at age 6 weeks and again at age 14 weeks and consisted of magnetic resonance imaging (MRI), technetium-diethylenetriaminepentaaceticacid ({sup 99m}Tc-DTPA) renography, renal technetium-dimercaptosuccinicacid ({sup 99m}Tc-DMSA) scintigraphy, and glomerular filtration rate (GFR) measurement. After relief of the partial obstruction, there was reduction of the pelvic diameter and improvement of urinary drainage. Global and relative kidney function was not significantly affected by either obstruction or its relief. Renal {sup 99m}Tc-DMSA scintigraphy showed a change in both the appearance of the kidney and a change in the distribution within kidneys even after relief of obstruction. This study shows that partial ureteric obstruction in young pigs may be associated with little effect on global and differential kidney function. However, even after relief of HN, the distribution of {sup 99m}Tc-DMSA in the kidney remains abnormal suggesting that a normal differential renal function may not represent a normal kidney. (orig.)

Reduced basis ANOVA methods for partial differential equations with high-dimensional random inputs

Energy Technology Data Exchange (ETDEWEB)

Liao, Qifeng, E-mail: liaoqf@shanghaitech.edu.cn [School of Information Science and Technology, ShanghaiTech University, Shanghai 200031 (China); Lin, Guang, E-mail: guanglin@purdue.edu [Department of Mathematics & School of Mechanical Engineering, Purdue University, West Lafayette, IN 47907 (United States)

2016-07-15

In this paper we present a reduced basis ANOVA approach for partial deferential equations (PDEs) with random inputs. The ANOVA method combined with stochastic collocation methods provides model reduction in high-dimensional parameter space through decomposing high-dimensional inputs into unions of low-dimensional inputs. In this work, to further reduce the computational cost, we investigate spatial low-rank structures in the ANOVA-collocation method, and develop efficient spatial model reduction techniques using hierarchically generated reduced bases. We present a general mathematical framework of the methodology, validate its accuracy and demonstrate its efficiency with numerical experiments.

Iterative Observer-based Estimation Algorithms for Steady-State Elliptic Partial Differential Equation Systems

KAUST Repository

Majeed, Muhammad Usman

2017-07-19

Steady-state elliptic partial differential equations (PDEs) are frequently used to model a diverse range of physical phenomena. The source and boundary data estimation problems for such PDE systems are of prime interest in various engineering disciplines including biomedical engineering, mechanics of materials and earth sciences. Almost all existing solution strategies for such problems can be broadly classified as optimization-based techniques, which are computationally heavy especially when the problems are formulated on higher dimensional space domains. However, in this dissertation, feedback based state estimation algorithms, known as state observers, are developed to solve such steady-state problems using one of the space variables as time-like. In this regard, first, an iterative observer algorithm is developed that sweeps over regular-shaped domains and solves boundary estimation problems for steady-state Laplace equation. It is well-known that source and boundary estimation problems for the elliptic PDEs are highly sensitive to noise in the data. For this, an optimal iterative observer algorithm, which is a robust counterpart of the iterative observer, is presented to tackle the ill-posedness due to noise. The iterative observer algorithm and the optimal iterative algorithm are then used to solve source localization and estimation problems for Poisson equation for noise-free and noisy data cases respectively. Next, a divide and conquer approach is developed for three-dimensional domains with two congruent parallel surfaces to solve the boundary and the source data estimation problems for the steady-state Laplace and Poisson kind of systems respectively. Theoretical results are shown using a functional analysis framework, and consistent numerical simulation results are presented for several test cases using finite difference discretization schemes.

Automatic differentiation bibliography

Energy Technology Data Exchange (ETDEWEB)

Corliss, G.F. [comp.

1992-07-01

This is a bibliography of work related to automatic differentiation. Automatic differentiation is a technique for the fast, accurate propagation of derivative values using the chain rule. It is neither symbolic nor numeric. Automatic differentiation is a fundamental tool for scientific computation, with applications in optimization, nonlinear equations, nonlinear least squares approximation, stiff ordinary differential equation, partial differential equations, continuation methods, and sensitivity analysis. This report is an updated version of the bibliography which originally appeared in Automatic Differentiation of Algorithms: Theory, Implementation, and Application.

Exploiting partial knowledge for efficient model analysis