Invariant discretizations of partial differential equations

Un algorithme permettant de discrétiser les équations aux dérivées partielles (EDP) tout en préservant leurs symétries de Lie est élaboré. Ceci est rendu possible grâce à l'utilisation de dérivées partielles discrètes se transformant comme les dérivées partielles continues sous l'action de groupes de Lie locaux. Dans les applications, beaucoup d'EDP sont invariantes sous l'action de transformations ponctuelles de Lie de dimension infinie qui font partie de ce que l'on désigne comme des pseudo-groupes de Lie. Afin d'étendre la méthode de discrétisation préservant les symétries à ces équations, une discrétisation des pseudo-groupes est proposée. Cette discrétisation a pour effet de transformer les symétries ponctuelles en symétries généralisées dans l'espace discret. Des schémas invariants sont ensuite créés pour un certain nombre d'EDP. Dans tous les cas, des tests numériques montrent que les schémas invariants approximent mieux leur équivalent continu que les différences finies standard.

A Decomposition and Noise Removal Method Combining Diffusion Equation and Wave Atoms for Textured Images

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Lie group symmetry analysis for granular media stress equations

The Airy stress function, although frequently employed in classical linear elasticity, does not receive similar usage for granular media problems. For plane strain quasi-static deformations of a cohesionless Coulomb–Mohr granular solid, a single nonlinear partial differential equation is formulated for the Airy stress function by combining the equilibrium equations with the yield condition. This has certain advantages from the usual approach, in which two stress invariants and a stress angle are introduced, and a system of two partial differential equations is needed to describe the flow. In the present study, the symmetry analysis of differential equations is utilised for our single partial differential equation, and by computing an optimal system of one-dimensional Lie algebras, a complete set of group-invariant solutions is derived. By this it is meant that any group-invariant solution of the governing partial differential equation (provided it can be derived via the classical symmetries method) may be obtained as a member of this set by a suitable group transformation. For general values of the parameters (angle of internal friction and gravity g) it is found there are three distinct classes of solutions which correspond to granular flows considered previously in the literature. For the two limiting cases of high angle of internal friction and zero gravity, the governing partial differential equation admit larger families of Lie point symmetries, and from these symmetries, further solutions are derived, many of which are new. Furthermore, the majority of these solutions are exact, which is rare for granular flow, especially in the case of gravity driven flows.

Unsteady Incompressible Two-Dimensional and Axisymmetric Turbulent Boundary Layer Flows

The unsteady turbulent incompressible boundary-layer flow over two-dimensional and axisymmetric bodies with pressure gradient has been studied. An eddy-viscosity model has been used to model the Reynolds shear stress. The unsteadiness is due to variations in the free stream velocity with time. The nonlinear partial differential equation with three independent variables governing the flow has been solved using Keller's Box method. The results indicate that the free stram velocity distribution exerts strong influence on the boundary-layer characteristics. The point of zero skin friction is found to move upstream as time increases.

A note on the equivalence of shock manifold equations

The shock manifold equation is a first order nonlinear partial differential equation, which describes the kinematics of a shockfront in an ideal gas with constant specific heats. However, it was found that there was more than one of these shock manifold equations, and the shock surface could be embedded in a one parameter family of surfaces, obtained as a solution of any of these shock manifold equations. Associated with each shock manifold equation is a set of characteristic curves called lsquoshock raysrsquo. This paper investigates the nature of various associated shock ray equations.

基于SVD和能量最小原则的图像自适应降噪算法

基于奇异值分解和能量最小原则,提出了一种自适应图像降噪算法,并给出了基于有界变差的能量降噪模型的代数形式。通过在矩阵范数意义下求能量最小,自适应确定去噪图像重构的奇异值个数。该算法的特点是将能量最小法则和奇异值分解结合起来,在代数空间中建立了一种自适应的图像降噪算法。与基于压缩比和奇异值分解的降噪方法相比,由于该算法避免了图像压缩比函数及其拐点的计算,因此具有快速去噪和简单可行的优点。实验结果证明,该算法是有效的。

Simulations magnétohydrodynamiques en régime idéal

Cette thèse s’intéresse à la modélisation magnétohydrodynamique des écoulements de fluides conducteurs d’électricité multi-échelles en mettant l’emphase sur deux applications particulières de la physique solaire: la modélisation des mécanismes des variations de l’irradiance via la simulation de la dynamo globale et la reconnexion magnétique. Les variations de l’irradiance sur les périodes des jours, des mois et du cycle solaire de 11 ans sont très bien expliquées par le passage des régions actives à la surface du Soleil. Cependant, l’origine ultime des variations se déroulant sur les périodes décadales et multi-décadales demeure un sujet controversé. En particulier, une certaine école de pensée affirme qu’une partie de ces variations à long-terme doit provenir d’une modulation de la structure thermodynamique globale de l’étoile, et que les seuls effets de surface sont incapables d’expliquer la totalité des fluctuations. Nous présentons une simulation globale de la convection solaire produisant un cycle magnétique similaire en plusieurs aspects à celui du Soleil, dans laquelle le flux thermique convectif varie en phase avec l’ ́energie magnétique. La corrélation positive entre le flux convectif et l’énergie magnétique supporte donc l’idée qu’une modulation de la structure thermodynamique puisse contribuer aux variations à long-terme de l’irradiance. Nous analysons cette simulation dans le but d’identifier le mécanisme physique responsable de la corrélation en question et pour prédire de potentiels effets observationnels résultant de la modulation structurelle. La reconnexion magnétique est au coeur du mécanisme de plusieurs phénomènes de la physique solaire dont les éruptions et les éjections de masse, et pourrait expliquer les températures extrêmes caractérisant la couronne. Une correction aux trajectoires du schéma semi-Lagrangien classique est présentée, qui est basée sur la solution à une équation aux dérivées partielles nonlinéaire du second ordre: l’équation de Monge-Ampère. Celle-ci prévient l’intersection des trajectoires et assure la stabilité numérique des simulations de reconnexion magnétique pour un cas de magnéto-fluide relaxant vers un état d’équilibre.

Dissipative three-wave structures in stimulated backscattering. I. A subluminous solitary attractor

We present a solitary solution of the three-wave nonlinear partial differential equation (PDE) model - governing resonant space-time stimulated Brillouin or Raman backscattering - in the presence of a cw pump and dissipative material and Stokes waves. The study is motivated by pulse formation in optical fiber experiments. As a result of the instability any initial bounded Stokes signal is amplified and evolves to a subluminous backscattered Stokes pulse whose shape and velocity are uniquely determined by the damping coefficients and the cw-pump level. This asymptotically stable solitary three-wave structure is an attractor for any initial conditions in a compact support, in contrast to the known superluminous dissipative soliton solution which calls for an unbounded support. The linear asymptotic theory based on the Kolmogorov-Petrovskii-Piskunov assertion allows us to determine analytically the wave-front slope and the subluminous velocity, which are in remarkable agreement with the numerical computation of the nonlinear PDE model when the dynamics attains the asymptotic steady regime. © 1997 The American Physical Society.

Convergence results for stochastic particle systems with social interaction

We consider stochastic individual-based models for social behaviour of groups of animals. In these models the trajectory of each animal is given by a stochastic differential equation with interaction. The social interaction is contained in the drift term of the SDE. We consider a global aggregation force and a short-range repulsion force. The repulsion range and strength gets rescaled with the number of animals N. We show that for N tending to infinity stochastic fluctuations disappear and a smoothed version of the empirical process converges uniformly towards the solution of a nonlinear, nonlocal partial differential equation of advection-reaction-diffusion type. The rescaling of the repulsion in the individual-based model implies that the corresponding term in the limit equation is local while the aggregation term is non-local. Moreover, we discuss the effect of a predator on the system and derive an analogous convergence result. The predator acts as an repulsive force. Different laws of motion for the predator are considered.

Interaction of Kelvin waves and nonlocality of energy transfer in superfluids

We argue that the physics of interacting Kelvin Waves (KWs) is highly nontrivial and cannot be understood on the basis of pure dimensional reasoning. A consistent theory of KW turbulence in superfluids should be based upon explicit knowledge of their interactions. To achieve this, we present a detailed calculation and comprehensive analysis of the interaction coefficients for KW turbuelence, thereby, resolving previous mistakes stemming from unaccounted contributions. As a first application of this analysis, we derive a local nonlinear (partial differential) equation. This equation is much simpler for analysis and numerical simulations of KWs than the Biot-Savart equation, and in contrast to the completely integrable local induction approximation (in which the energy exchange between KWs is absent), describes the nonlinear dynamics of KWs. Second, we show that the previously suggested Kozik-Svistunov energy spectrum for KWs, which has often been used in the analysis of experimental and numerical data in superfluid turbulence, is irrelevant, because it is based upon an erroneous assumption of the locality of the energy transfer through scales. Moreover, we demonstrate the weak nonlocality of the inverse cascade spectrum with a constant particle-number flux and find resulting logarithmic corrections to this spectrum.

Advances in Nonlinear Complexity Analysis for Partial Differential Equations

1 p. -- [Editorial Material]

Initial boundary value problems for the nonlinear Schrödinger equation

A new spectral method for solving initial boundary value problems for linear and integrable nonlinear partial differential equations in two independent variables is applied to the nonlinear Schrödinger equation and to its linearized version in the domain {x≥l(t), t≥0}. We show that there exist two cases: (a) if l″(t)<0, then the solution of the linear or nonlinear equations can be obtained by solving the respective scalar or matrix Riemann-Hilbert problem, which is defined on a time-dependent contour; (b) if l″(t)>0, then the Riemann-Hilbert problem is replaced by a respective scalar or matrix problem on a time-independent domain. In both cases, the solution is expressed in a spectrally decomposed form.

On Modified Method of Simplest Equation for Obtaining Exact Solutions of Nonlinear PDEs: Case of Elliptic Simplest Equation

2010 Mathematics Subject Classification: 74J30, 34L30.