Stochastic perturbations to dynamical systems: a response theory approach

Using the formalism of the Ruelle response theory, we study how the invariant measure of an Axiom A dynamical system changes as a result of adding noise, and describe how the stochastic perturbation can be used to explore the properties of the underlying deterministic dynamics. We first find the expression for the change in the expectation value of a general observable when a white noise forcing is introduced in the system, both in the additive and in the multiplicative case. We also show that the difference between the expectation value of the power spectrum of an observable in the stochastically perturbed case and of the same observable in the unperturbed case is equal to the variance of the noise times the square of the modulus of the linear susceptibility describing the frequency-dependent response of the system to perturbations with the same spatial patterns as the considered stochastic forcing. This provides a conceptual bridge between the change in the fluctuation properties of the system due to the presence of noise and the response of the unperturbed system to deterministic forcings. Using Kramers-Kronig theory, it is then possible to derive the real and imaginary part of the susceptibility and thus deduce the Green function of the system for any desired observable. We then extend our results to rather general patterns of random forcing, from the case of several white noise forcings, to noise terms with memory, up to the case of a space-time random field. Explicit formulas are provided for each relevant case analysed. As a general result, we find, using an argument of positive-definiteness, that the power spectrum of the stochastically perturbed system is larger at all frequencies than the power spectrum of the unperturbed system. We provide an example of application of our results by considering the spatially extended chaotic Lorenz 96 model. These results clarify the property of stochastic stability of SRB measures in Axiom A flows, provide tools for analysing stochastic parameterisations and related closure ansatz to be implemented in modelling studies, and introduce new ways to study the response of a system to external perturbations. Taking into account the chaotic hypothesis, we expect that our results have practical relevance for a more general class of system than those belonging to Axiom A.

A nonparametric ensemble transform method for Bayesian inference

Many applications, such as intermittent data assimilation, lead to a recursive application of Bayesian inference within a Monte Carlo context. Popular data assimilation algorithms include sequential Monte Carlo methods and ensemble Kalman filters (EnKFs). These methods differ in the way Bayesian inference is implemented. Sequential Monte Carlo methods rely on importance sampling combined with a resampling step, while EnKFs utilize a linear transformation of Monte Carlo samples based on the classic Kalman filter. While EnKFs have proven to be quite robust even for small ensemble sizes, they are not consistent since their derivation relies on a linear regression ansatz. In this paper, we propose another transform method, which does not rely on any a priori assumptions on the underlying prior and posterior distributions. The new method is based on solving an optimal transportation problem for discrete random variables. © 2013, Society for Industrial and Applied Mathematics

Basic theory behind parameterizing atmospheric convection

Last fall, a network of the European Cooperation in Science and Technology (COST), called “Basic Concepts for Convection Parameterization in Weather Forecast and Climate Models” (COST Action ES0905; see http://w3.cost.esf.org/index.php?id=205&action_number=ES0905), organized a 10-day training course on atmospheric convection and its parameterization. The aim of the workshop, held on the island of Brac, Croatia, was to help young scientists develop an in-depth understanding of the core theory underpinning convection parameterizations. The speakers also sought to impart an appreciation of the various approximations, compromises, and ansatz necessary to translate theory into operational practice for numerical models.

EXACT AND QUASI-EXACT MODELS OF STRANGE STARS

We construct and compare in this work a variety of simple models for strange stars, namely, hypothetical self-bound objects made of a cold stable version of the quark-gluon plasma. Exact, quasi-exact and numerical models are examined to find the most economical description for these objects. A simple and successful parametrization of them is given in terms of the central density, and the differences among the models are explicitly shown and discussed. In particular, we present a model starting with a Gaussian ansatz for the density profile that provides a very accurate and almost complete analytical integration of the problem, modulo a small difference for one of the metric potentials.

A Hopping Mechanism for Cargo Transport by Molecular Motors on Crowded Microtubules

Most models designed to study the bidirectional movement of cargos as they are driven by molecular motors rely on the idea that motors of different polarities can be coordinated by external agents if arranged into a motor-cargo complex to perform the necessary work Gross, Hither and yon: a review of bidirectional microtubule-based transport (Gross in Phys. Biol. 1:R1-R11, 2004). Although these models have provided us with important insights into these phenomena, there are still many unanswered questions regarding the mechanisms through which the movement of the complex takes place on crowded microtubules. For example (i) how does cargo-binding affect motor motility? and in connection with that-(ii) how does the presence of other motors (and also other cargos) on the microtubule affect the motility of the motor-cargo complex? We discuss these questions from a different perspective. The movement of a cargo is conceived here as a hopping process resulting from the transference of cargo between neighboring motors. In the light of this, we examine the conditions under which cargo might display bidirectional movement even if directed by motors of a single polarity. The global properties of the model in the long-time regime are obtained by mapping the dynamics of the collection of interacting motors and cargos into an asymmetric simple exclusion process (ASEP) which can be resolved using the matrix ansatz introduced by Derrida (Derrida and Evans in Nonequilibrium Statistical Mechanics in One Dimension, pp. 277-304, 1997; Derrida et al. in J. Phys. A 26: 1493-1517, 1993).

Self-dual hopfions

We construct static and time-dependent exact soliton solutions with nontrivial Hopf topological charge for a field theory in 3 + 1 dimensions with the target space being the two dimensional sphere S(2). The model considered is a reduction of the so-called extended Skyrme-Faddeev theory by the removal of the quadratic term in derivatives of the fields. The solutions are constructed using an ansatz based on the conformal and target space symmetries. The solutions are said self-dual because they solve first order differential equations which together with some conditions on the coupling constants, imply the second order equations of motion. The solutions belong to a sub-sector of the theory with an infinite number of local conserved currents. The equation for the profile function of the ansatz corresponds to the Bogomolny equation for the sine-Gordon model.

Static Hopfions in the extended Skyrme-Faddeev model

We construct static soliton solutions with non-zero Hopf topological charges to a theory which is an extension of the Skyrme-Faddeev model by the addition of a further quartic term in derivatives. We use an axially symmetric ansatz based on toroidal coordinates, and solve the resulting two coupled non-linear partial differential equations in two variables by a successive over-relaxation (SOR) method. We construct numerical solutions with Hopf charge up to four, and calculate their analytical behavior in some limiting cases. The solutions present an interesting behavior under the changes of a special combination of the coupling constants of the quartic terms. Their energies and sizes tend to zero as that combination approaches a particular special value. We calculate the equivalent of the Vakulenko and Kapitanskii energy bound for the theory and find that it vanishes at that same special value of the coupling constants. In addition, the model presents an integrable sector with an in finite number of local conserved currents which apparently are not related to symmetries of the action. In the intersection of those two special sectors the theory possesses exact vortex solutions (static and time dependent) which were constructed in a previous paper by one of the authors. It is believed that such model describes some aspects of the low energy limit of the pure SU(2) Yang-Mills theory, and our results may be important in identifying important structures in that strong coupling regime.

Axially symmetric soliton solutions in a Skyrme-Faddeev-type model with Gies`s extension

We construct static soliton solutions with non-zero Hopf topological charges to a theory which is the extended Skyrme-Faddeev model with a further quartic term in derivatives. We use an axially symmetric ansatz based on toroidal coordinates, and solve the resulting two coupled nonlinear partial differential equations in two variables by a successive over-relaxation method. We construct numerical solutions with the Hopf charge up to 4. The solutions present an interesting behavior under the changes of a special combination of the coupling constants of the quartic terms.

Simulação numérica da combustão turbulenta de gás natural em câmara cilíndrica

O presente trabalho apresenta uma modelagem detalhada de processos de combustao turbulentos para um jato concentrico de combustıvel e ar. A modelagem é fundamentada nas equacões de conservacão de massa, de quantidade de movimento, de energia e de espécies quımicas. A turbulencia é resolvida pela utilizacão do modelo k- padrão. Dois modelos de reacões quımicas são apresentados. O modelo SCRS – Simple Chemically-Reacting Systems, que assume taxas instantâneas de reacões quımicas. Também é abordado o modelo E-A – Eddy Breakup - Arrhenius, que assume taxas finitas de reacões quımicas. A radiacão térmica, fenômeno de grande importância devido as altas temperaturas alcancadas em processos de combustão,é modelada através do Método das Zonas. O modelo da soma ponderada de gases cinzas – WSGGM,é usado para determinar o espectro de emissão e absorcão dos gases no processo. Para a solucão destas equacões diferenciais, juntamente com os modelos de turbulência, de reaçõoes químicas e radiação térmica, faz-se o uso do Método dos Volumes Finitos. Para validar a modelagem apresentada resolve-se o processo de combustão em uma câmara cilíndrica. A câmara de combustão usada áa mesma abordada no First Workshop on Aerodynamics of Steady State Combustion Chambers and Furnaces, organizado pela ERCORTAC - European Research Community On Flow Turbulence And Combustion, em outubro de 1994, que apresenta dados experimentais de temperatura e concentração das espécies químicas para várias posições de interesse no interior da câmara. Utiliza-se o gás natural como combustível e o ar atmosférico como oxidante. O processo de combustão sem pré-mistura é resolvido para a condição de excesso de combustível de 5 % para ambos os modelos, onde o gás natural é injetado por um duto circular central, e o ar atmosférico por um orifício anular externo a esse duto, no mesmo plano Uma reação química não estagiada é assumida para o modelo SCRS. Para o modelo E-A duas situações são resolvidas: combustão não estagiada, com uma etapa global de reação química; e reação quımica estagiada, com duas etapas globais. Os resultados obtidos com o modelo SCRS para a distribuição de temperaturas, em termos de tendências gerais, são razoáveis. Já as concentrações de espécies químicas não apresentam dados satisfatórios para este modelo. Para o modelo E-A os resultados apresentam boa concordância com os dados experimentais, principalmente para a situação em que o processo de combustão é assumido em duas etapas globais. ´E analisado em detalhe o papel desempenhado pela transferencia de calor por radiacao, com meio participante. Para melhor verificar as trocas de calor, assume-se uma camara de combustao cilındrica com paredes d’agua. A injecao do combustıvel e do oxidante e feita atraves de um queimador central, semelhante ao usado para validar a modelagem, porem com dois orifıcios concentricos para injecao de combustıvel. Nesta situação o efeito do turbilhonamento (swril), assumido como 20 % da velocidade axial de entrada, sobre a injecao de ar e computado atraves da condicao contorno da equacao de conservacao da quantidade de movimento angular. Nesta fase apenas o modelo E-A, com duas etapas globais de reacoes quımicas, e considerado, ja que o mesmo apresenta os melhores resultados. O processo de combustao e simulado com e sem a presenca da radiacao termica. Verifica-se que a presenca da radiacao termica homogeneiza a temperatura dos gases no interior da camara. Com isso verifica-se tambem alterações nas taxas de reacoes quımicas, modificando a magnitude das fracoes das especies quımicas Quando a radiacao termica e considerada efeitos de extinção local da chama sao verificados nas regioes de temperaturas mais altas, diminuindo o consumo de oxigenio e aumentando a producao de monoxido de carbono, caracterizando assim uma combustao incompleta. Em algumas situacoes tem-se uma variacao de temperatura de ate 500 K, a montante da chama. A radiacao termica tambem aumenta a taxa de transferencia de calor dos gases quentes para as paredes da camara, e desta para o seu exterior. Com os resultados obtidos a partir desta modelagem e possıvel determinar o perfil da zona de combustao, a distribuicao de concentracoes de especies quımicas, o campo de velocidades e as taxas de transferencia de calor para as paredes da camara de combustao, total, por conveccao superficial e por radiacao. Estes resultados sao de extrema importancia para prever a performance de camaras de combustao, assim como auxiliar na sua otimizacao.

Uma solução híbrida para a equação de Fokker-Planck dependente da energia usando a técnica da transformada de Laplace e diferenças finitas

Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Estabilidade não linear de um feixe de partículas carregadas sujeito a um campo magnético focalizador

Uma análise da estabilidade não-linear de feixes com seção transversal circular considerando perturbações sem simetria axial é executada. Mostra-se que as oscilações simétricas oficialmente circulares de um feixe podem induzir oscilações não-lineares do tipo anti-simétrica(elípticas), com um conseqüente aumento do tamanho do feixe ao logo de uma direção preferencial. O mecanismo da instabilidade e sua relevância às perdas de partícula no feixe são discutidos.

Correlação linear entre componentes da produção e produtividade do arroz de terras altas em sistema plantio direto

The rice is one of the main sources of the humanity's feeding. During the agricultural year 2009/2010, in Selviria County, Mato Grosso do Sul State, in the Brazilian Savannah, an experiment was installed with rice upland in a Dystropherric Red Latosol (Typic Acrustox) under no-tillage, irrigated by central pivot, with the purpose of selecting the best components production to explain the variability the irrigated rice yield upland. The geostatistical grid was installed, to collect the data, with 120 sampling points, in an area of 3.0 ha and and homogeneous slope of 0.055 m m(-1). The medium rice yield was of the 5980 kg ha(-1). For the simple lineal regressions, the number of spikelets grenades for panicle presented the best direct potential correlation with the yield rice, given for: PGO = 115,5.NEG(0,770). However, for the multiple lineal regressions, the equation equacao PGO = 2754,30-411,55.NEG-461,07. NEC+436,59. NET it was the one that better she came to esteem the yield rice. However, spatial, it was not possible to establish correlation between the yield rice and the components production, once none of those it presented spatial dependence in their data.

Variational supersymmetric approach to evaluate Fokker-Planck probability

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Strength of the trilinear Higgs boson coupling in technicolor models

We discuss the strength of the trilinear Higgs boson coupling in technicolor (or composite) models in a model independent way. The coupling is determined as a function of a very general ansatz for the technicolor self-energy, and turns out to be equal or smaller than the one of the Standard Model Higgs boson depending on the dynamics of the theory. (c) 2006 Published by Elsevier B.V.