Holder Continuous Solutions of Active Scalar Equations

We consider active scalar equations $\partial_t \theta + \nabla \cdot (u \, \theta) = 0$, where $u = T[\theta]$ is a divergence-free velocity field, and $T$ is a Fourier multiplier operator with symbol $m$. We prove that when $m$ is not an odd function of frequency, there are nontrivial, compactly supported solutions weak solutions, with H\"older regularity $C^{1/9-}_{t,x}$. In fact, every integral conserving scalar field can be approximated in ${\cal D}'$ by such solutions, and these weak solutions may be obtained from arbitrary initial data. We also show that when the multiplier $m$ is odd, weak limits of solutions are solutions, so that the $h$-principle for odd active scalars may not be expected.Comment: 61 page

A methodology for airplane parameter estimation and confidence interval determination in nonlinear estimation problems

An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. With the fitted surface, sensitivity information can be updated at each iteration with less computational effort than that required by either a finite-difference method or integration of the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, and thus provides flexibility to use model equations in any convenient format. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. The degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels and to predict the degree of agreement between CR bounds and search estimates

A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection-diffusion PDE's coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton-Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that on be able to solve chemical equilibrium problems (and compute derivatives), without having to know the solution method. An additional advantage of the Newton-Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.Comment: Computational Geosciences (2009) http://www.springerlink.com/content/933p55085742m203/?p=db14bb8c399b49979ba8389a3cae1b0f&pi=1

Steady and Stable: Numerical Investigations of Nonlinear Partial Differential Equations

Excerpt: Mathematics is a language which can describe patterns in everyday life as well as abstract concepts existing only in our minds. Patterns exist in data, functions, and sets constructed around a common theme, but the most tangible patterns are visual. Visual demonstrations can help undergraduate students connect to abstract concepts in advanced mathematical courses. The study of partial differential equations, in particular, benefits from numerical analysis and simulation

An optimal three-point eighth-order iterative method without memory for solving nonlinear equations with its dynamics

We present a three-point iterative method without memory for solving nonlinear equations in one variable. The proposed method provides convergence order eight with four function evaluations per iteration. Hence, it possesses a very high computational efficiency and supports Kung and Traub's conjecture. The construction, the convergence analysis, and the numerical implementation of the method will be presented. Using several test problems, the proposed method will be compared with existing methods of convergence order eight concerning accuracy and basin of attraction. Furthermore, some measures are used to judge methods with respect to their performance in finding the basin of attraction.Comment: arXiv admin note: substantial text overlap with arXiv:1508.0174