Nonlinear data assimilation in geosciences: an extremely efficient particle filter

Almost all research fields in geosciences use numerical models and observations and combine these using data-assimilation techniques. With ever-increasing resolution and complexity, the numerical models tend to be highly nonlinear and also observations become more complicated and their relation to the models more nonlinear. Standard data-assimilation techniques like (ensemble) Kalman filters and variational methods like 4D-Var rely on linearizations and are likely to fail in one way or another. Nonlinear data-assimilation techniques are available, but are only efficient for small-dimensional problems, hampered by the so-called ‘curse of dimensionality’. Here we present a fully nonlinear particle filter that can be applied to higher dimensional problems by exploiting the freedom of the proposal density inherent in particle filtering. The method is illustrated for the three-dimensional Lorenz model using three particles and the much more complex 40-dimensional Lorenz model using 20 particles. By also applying the method to the 1000-dimensional Lorenz model, again using only 20 particles, we demonstrate the strong scale-invariance of the method, leading to the optimistic conjecture that the method is applicable to realistic geophysical problems. Copyright c 2010 Royal Meteorological Society

The influence of the carbon surface chemical composition on Dubinin-Astakhov equation parameters calculated from SF(6) adsorption data-grand canonical Monte Carlo simulation

Using grand canonical Monte Carlo simulation we show, for the first time, the influence of the carbon porosity and surface oxidation on the parameters of the Dubinin-Astakhov (DA) adsorption isotherm equation. We conclude that upon carbon surface oxidation, the adsorption decreases for all carbons studied. Moreover, the parameters of the DA model depend on the number of surface oxygen groups. That is why in the case of carbons containing surface polar groups, SF(6) adsorption isotherm data cannot be used for characterization of the porosity.

Classifying general nonlinear force laws in cell-based models via the continuum limit

though discrete cell-based frameworks are now commonly used to simulate a whole range of biological phenomena, it is typically not obvious how the numerous different types of model are related to one another, nor which one is most appropriate in a given context. Here we demonstrate how individual cell movement on the discrete scale modeled using nonlinear force laws can be described by nonlinear diffusion coefficients on the continuum scale. A general relationship between nonlinear force laws and their respective diffusion coefficients is derived in one spatial dimension and, subsequently, a range of particular examples is considered. For each case excellent agreement is observed between numerical solutions of the discrete and corresponding continuum models. Three case studies are considered in which we demonstrate how the derived nonlinear diffusion coefficients can be used to (a) relate different discrete models of cell behavior; (b) derive discrete, intercell force laws from previously posed diffusion coefficients, and (c) describe aggregative behavior in discrete simulations.

The Dirichlet-to-Neumann map for the elliptic sine Gordon

We analyse the Dirichlet problem for the elliptic sine Gordon equation in the upper half plane. We express the solution $q(x,y)$ in terms of a Riemann-Hilbert problem whose jump matrix is uniquely defined by a certain function $b(\la)$, $\la\in\R$, explicitly expressed in terms of the given Dirichlet data $g_0(x)=q(x,0)$ and the unknown Neumann boundary value $g_1(x)=q_y(x,0)$, where $g_0(x)$ and $g_1(x)$ are related via the global relation $\{b(\la)=0$, $\la\geq 0\}$. Furthermore, we show that the latter relation can be used to characterise the Dirichlet to Neumann map, i.e. to express $g_1(x)$ in terms of $g_0(x)$. It appears that this provides the first case that such a map is explicitly characterised for a nonlinear integrable {\em elliptic} PDE, as opposed to an {\em evolution} PDE.

Study of the nonlinear optical properties of a crosslinkable chromophore for use in guest host polymer films

A guest/host material system in which the guest molecule is a functionalized, optically nonlinear, chromophore is described. A verification of the crosslinking process, an assessment of the nonlinear properties of the chromophore, using Solvatochromic methods, and an investigation of the electric field induced molecular orientation using second-harmonic generation are included.

A comparison of photoinduced poling and thermal poling of azo‐dye‐doped polymer films for second order nonlinear optical applications

Photoinduced poling (PIP) is a new technique which allows the room‐temperature preparation of guest/host polymer films exhibiting significant polar order for nonlinear optical applications. We report a comparison of this novel technique with the conventional electrode poling procedure performed at the glass transition temperature of the polymer using disperse red 1/poly(methylmethacrylate) films. In particular, in situ second harmonic generation measurements show that levels of polar order achieved using these two techniques are similar. In contrast, the stability of the polar order is reduced by up to 20 times in terms of the decay time constant in films prepared using PIP although the stability is very dependent upon the temperature at which the poling was performed.

A stress relaxation model for the viscoelastic solids based on the steady-state creep equation

Creep and stress relaxation are inherent mechanical behaviors of viscoelastic materials. It is considered that both are different performances of one identical physical phenomenon. The relationship between the decay stress and time during stress relaxation has been derived from the power law equation of the steady-state creep. The model was used to analyse the stress relaxation curves of various different viscoelastic materials (such as pure polycrystalline ice, polymers, foods, bones, metal, animal tissues, etc.). The calculated results using the theoretical model agree with the experimental data very well. Here we show that the new mathematical formula is not only simple but its parameters have the clear physical meanings. It is suitable to materials with a very broad scope and has a strong predictive ability.

Stable Nonlinear Receding Horizon Regulator Using RBF Neural Network Models

The general stability theory of nonlinear receding horizon controllers has attracted much attention over the last fifteen years, and many algorithms have been proposed to ensure closed-loop stability. On the other hand many reports exist regarding the use of artificial neural network models in nonlinear receding horizon control. However, little attention has been given to the stability issue of these specific controllers. This paper addresses this problem and proposes to cast the nonlinear receding horizon control based on neural network models within the framework of an existing stabilising algorithm.

Evidence of dispersion relations for the nonlinear response of the Lorenz 63 system

Along the lines of the nonlinear response theory developed by Ruelle, in a previous paper we have proved under rather general conditions that Kramers-Kronig dispersion relations and sum rules apply for a class of susceptibilities describing at any order of perturbation the response of Axiom A non equilibrium steady state systems to weak monochromatic forcings. We present here the first evidence of the validity of these integral relations for the linear and the second harmonic response for the perturbed Lorenz 63 system, by showing that numerical simulations agree up to high degree of accuracy with the theoretical predictions. Some new theoretical results, showing how to derive asymptotic behaviors and how to obtain recursively harmonic generation susceptibilities for general observables, are also presented. Our findings confirm the conceptual validity of the nonlinear response theory, suggest that the theory can be extended for more general non equilibrium steady state systems, and shed new light on the applicability of very general tools, based only upon the principle of causality, for diagnosing the behavior of perturbed chaotic systems and reconstructing their output signals, in situations where the fluctuation-dissipation relation is not of great help.

Exact design manifold control of a class of nonlinear singularly perturbed systems

The integral manifold approach captures from a geometric point of view the intrinsic two-time-scale behavior of singularly perturbed systems. An important class of nonlinear singularly perturbed systems considered in this note are fast actuator-type systems. For a class of fast actuator-type systems, which includes many physical systems, an explicit corrected composite control, the sum of a slow control and a fast control, is derived. This corrected control will steer the system exactly to a required design manifold.

Exact design manifold control of a class of nonlinear singularly perturbed systems

The integral manifold approach captures from a geometric point of view the intrinsic two-time-scale behavior of singularly perturbed systems. An important class of nonlinear singularly perturbed systems considered in this note are fast actuator-type systems. For a class of fast actuator-type systems, which includes many physical systems, an explicit corrected composite control, the sum of a slow control and a fast control, is derived. This corrected control will steer the system exactly to a required design manifold.

Bifurcation Analysis of Eigenstructure Assignment Control in a Simple Nonlinear Aircraft Model

Aircraft systems are highly nonlinear and time varying. High-performance aircraft at high angles of incidence experience undesired coupling of the lateral and longitudinal variables, resulting in departure from normal controlled � ight. The construction of a robust closed-loop control that extends the stable and decoupled � ight envelope as far as possible is pursued. For the study of these systems, nonlinear analysis methods are needed. Previously, bifurcation techniques have been used mainly to analyze open-loop nonlinear aircraft models and to investigate control effects on dynamic behavior. Linear feedback control designs constructed by eigenstructure assignment methods at a � xed � ight condition are investigated for a simple nonlinear aircraft model. Bifurcation analysis, in conjunction with linear control design methods, is shown to aid control law design for the nonlinear system.