Dans le contexte du stockage des déchets radioactifs en milieu poreux, nous considérons l’équation de convection-diffusion instationnaire et sa discrétisation par des méthodes numériques. La discontinuité des paramètres physiques et la variabilité des échelles d’espace et de temps conduisent à utiliser des discrétisations différentes en temps et en espace dans différentes régions du domaine. Nous choisissons dans cette thèse le schéma volumes finis en dualité discrète (DDFV) et le schéma de Galerkin Discontinu en temps couplés à une méthode de décomposition de domaine de Schwarz de type relaxation d’ondes optimisées (OSWR), ce qui permet de traiter des maillages espace-temps non conformes. La principale difficulté réside dans l’obtention d’une discrétisation amont du flux convectif qui reste locale à un sous-domaine et telle que le schéma monodomaine soit équivalent au schéma multidomaine. Ces difficultés sont appréhendées d’abord en une dimension d’espace où différentes discrétisations sont étudiées. Le schéma retenu introduit une inconnue hybride sur les interfaces entre cellules. L’idée du décentrage amont par rapport à cette inconnue hybride est reprise en dimension deux d’espace, et adaptée au schéma DDFV. Le caractère bien posé de ce schéma et d’un schéma multidomaine équivalent est montré. Ce dernier est résolu par un algorithme OSWR dont la convergence est prouvée. Les paramètres optimisés des conditions de Robin sont obtenus par l'étude du taux de convergence continu ou discret. Différents cas-tests, dont l’un est inspiré du stockage des déchets nucléaires, illustrent ces résultats. / In the context of nuclear waste repositories, we consider the numerical discretization of the non stationary convection diffusion equation. Discontinuous physical parameters and heterogeneous space and time scales lead us to use different space and time discretizations in different parts of the domain. In this work, we choose the discrete duality finite volume (DDFV) scheme and the discontinuous Galerkin scheme in time, coupled by an optimized Scwharz waveform relaxation (OSWR) domain decomposition method, because this allows the use of non-conforming space-time meshes. The main difficulty lies in finding an upwind discretization of the convective flux which remains local to a sub-domain and such that the multidomain scheme is equivalent to the monodomain one. These difficulties are first dealt with in the one-dimensional context, where different discretizations are studied. The chosen scheme introduces a hybrid unknown on the cell interfaces. The idea of upwinding with respect to this hybrid unknown is extended to the DDFV scheme in the two-dimensional setting. The well-posedness of the scheme and of an equivalent multidomain scheme is shown. The latter is solved by an OSWR algorithm, the convergence of which is proved. The optimized parameters in the Robin transmission conditions are obtained by studying the continuous or discrete convergence rates. Several test-cases, one of which inspired by nuclear waste repositories, illustrate these results.

Pendant ces trois ans, nous nous sommes intéressés à quelques sujets choisis en optique et en électrodynamique quantiques. Le fil rouge de nos interrogations est la fonction d’onde du photon. Les expériences d’optique et d’électrodynamique quantique peuvent-elles être décrites de manière simple, dans l’espace des positions, à l’aide d’une fonction d’onde décrivant le ou les photon(s) impliqués dans l’expérience ? Ce n’est pas entièrement évident :la description usuelle des photons se fait dans l’espace réciproque des vecteurs d’onde. Mais ces expériences gagnent à être décrites par la mécanique ondulatoire en représentation position, comme cela est fait dans les manuels de mécanique quantique pour des situations impliquant des particules massives. De surcroît, une expérience récente[1] a conduit à l’observation de trajectoires de photons uniques à travers un interféromètre à deux fentes d’Young.Pour essayer de décrire formellement ces trajectoires, il est naturel de formuler une mécanique ondulatoire pour les photons. Nous avons donc examiné en détail la construction formelle de la fonction d’onde du photon, un objet qui est resté peu étudié jusqu’aux années 1990. Nous avons également étudié les propriétés de la fonction d’onde du photon en présence de sources, et considéré pour ce faire divers systèmes quantiques ouverts (en interaction). Nous avons vu qu’il existe, en principe, une infinité de possibilités pour le choix de la fonction d’onde du photon.Nous avons mis en évidence un certain nombre de critères sur la base desquels il apparaît que seuls trois choix parmi tous ceux possibles sont intéressants, l’un d’entre eux ramenant à un objet introduit par Glauber [2] pour étudier la détection de la lumière et les corrélations du champ électromagnétique. Nous avons également vu qu’en l’absence de sources l’équation quantique de propagation des photons est formellement identique aux équations de Maxwell.À bas nombre de photons, le formalisme de la fonction d’onde peut se révéler très pratique. Nous avons adapté l’approche aux systèmes en interaction, en nous intéressant dans un premier temps à l’électrodynamique quantique1en cavité [3], en particulier aux expériences réalisées par le groupe de Serge Haroche [4]. Nous avons proposé un modèle simple pour la description des photons dans les cavités d’électrodynamique. À l’aide de ce modèle, et de la fonction d’onde du photon, nous avons étudié la propagation des photons s’échappant de la cavité. Nous avons également construit l’équation maîtresse de Lindblad sans introduire de sauts quantiques non unitaires (voir également [5]). Nous nous sommes enfin intéressés à la question de l’évolution spatiotemporelle d’un photon émis lors d’une désexcitation d’un électron atomique. Après avoir étudié soigneusement la dynamique de la désexcitation de l’électron, notamment aux temps très courts [6, 7], nous nous sommes attachés à décrire, aussi rigoureusement que possible, le champ électromagnétique émis. Celui-ci, de manière surprenante, n’évolue pas causalement. Si cela n’est pas entièrement inattendu au vu du théorème de Hegerfeldt, qui stipule [8] que la causalité est exclue pour les systèmes décrits par un Hamiltonien dont le spectre est borné inférieurement, nous avons identifié [9] deux autres sources de non-causalité, l’une, prédite qualitativement par Shirokov [10], et l’autre, entièrement nouvelle à notre connaissance, et dont la compréhension reste à affiner. / During these three years we focused on several topics in quantum otpics and quantum electrodynamics. A central theme in our investigations is that of the photon wave function. Can quantum optics and quantum electrodynamics experiments be described simply, in position space, with the help of a wave function describing the photon(s) featured in the experiment ? The answer to that question is not quite obvious: the usual description of photons takes place in the reciprocal space of wave vectors. But these experiments call for a wave mechanical description in the position representation, as is done in quantum mechanics textbooks in situations featuring massive particles. Moreover, in a recent experiment [1], single photon trajectories through a Young two-slit setup have been observed. In order to try and describe these trajectories formally, it is natural to build a wave mechanical formalism for photons. We therefore studied in detail the formal construction of the photon wave function, an object which was little studied until the 1990s. We also studied the properties of the photon wave function in the presence of sources.To do that, we considered several open (interacting) quantum systems. We saw that there exists in principle an infinite number of possibilities when defining the photon wave function. We emphasised several criteria on the basis of which it appears that only three choices for the wave function are interesting. One of them coincides with an object introduced and used by Glauber [2] to study light detection andthe correlations of the electromagnetic field in the quantum regime. We also saw that, in the absence of sources, the propagation equation for a single photon is formally equivalent to Maxwell’s equations. At low photon numbers, the wave function formalism can be very useful. We adapted it to interacting systems,first, to cavity quantum electrodynamics (QED) [3], in particular to the experiments carried out by Serge Haroche’s group [4]. We proposed a simple model to describe photons in QED cavities. With this model, and with the helpof the photon wave function, we studied the propagation of photons escaping a cavity. We also constructed the Lindblad master equation without introducing nonunitary quantum jumps (also see [5]). We finally investigated the spacetime evolution of a photon which is emitted during the decay of an atomic electron. After having carefully studied the dynamics of the electronic decay, especially at very short times [6, 7], we set out to describe the emitted electromagnetic field as rigorously as possible. This emitted field, surprisingly, does not evolve causally. Though this is not entirely unexpected in view of Hegerfeldt’s theorem, which states [8] that causality is impossible for quantum systems which are described by a Hamiltonian with a spectrum which is bounded by below, we identified [9] two other sources of non causality. One of them was predicted qualitatively by Shirokov [10], while the other one, which is completely new as far as we can tell, is still to be better understood

The need for debugging tools for declarative programming languages has increased due to the rapid development of modeling and simulation tools/programs. Declarative equation-based programming languages have the problem of equation systems being over-, or under-constrained. This means that the system of equations has more equations than variables or more variables than equations respectively, making the system of equations unsolvable. In this study a static debugger is implemented in OpenModelica compiler for the equation-based programming language Modelica to make it easier for the programmer or modeler to locate the equation/s causing the unconstrained system of equations. The debugging techniques used by the debugger are developed by Peter Bunus. Those techniques are able to detect unconstrained systems of equations and give solutions by identifying the minimal set ofequation/s that should be removed or which variable/s should be added to an equation/s to make the system solvable. In this study the debugging techniques for detecting and giving a solution for over-constrained system of equations are shown suitable to be used for the programming language Modelica in the OpenModelica compiler.

L’étude des ondes scélérates est un sujet en plein essor principalement en océanographie mais également dans d’autres domaines. Dans cette thèse, je construis par transformation de Darboux une famille multi-paramétrique de solutions quasi-rationnelles lisses de l’équation de Schödinger non linéaire qui présentent un comportement d’ondes scélérates. Pour un choix générique de paramètres les solutions de deuxième ordre donnent un modèle de "trois sœurs" (une succession de trois vagues plus hautes que prévues) alors que pour un choix particulier de paramètres on obtient les solutions présentées par Akhmediev et al. dans une série d’articles de 2009. Ces solutions me permettent ensuite de construire des solutions rationnelles de l’équation KP-I qui décrit le mouvement des vagues dans une eau peu profonde. / The study of rogue waves is a booming topic mainly in oceanography but also in other fields. In this thesis I construct via Darboux transform a multi-parametric family of smooth quasi-rational solutions of the nonlinear Schödinger equation that present a behavior of rogue waves. For a general choice of parameters the second-order solutions give a model of "three sisters" (three higher than expected waves in a row) while for a particular choice of parameters we obtain the solutions given by Akhmediev et al. in a serie of articles in 2009. Then these solutions allow me to construct rational solutions of the KP-I equation that describe waves in shallow water.

Les problèmes inverses liés à des processus physiques sont d'une grande importance dans la plupart des domaines liés au traitement du signal, tels que la tomographie, l'acoustique, les communications sans fil, le radar, l'imagerie médicale, pour n'en nommer que quelques uns. Dans le même temps, beaucoup de ces problèmes soulèvent des défis en raison de leur nature mal posée. Par ailleurs, les signaux émanant de phénomènes physiques sont souvent gouvernées par des lois s'exprimant sous la forme d'équations aux dérivées partielles (EDP) linéaires, ou, de manière équivalente, par des équations intégrales et leurs fonctions de Green associées. De plus, ces phénomènes sont habituellement induits par des singularités, apparaissant comme des sources ou des puits d'un champ vectoriel. Dans cette thèse, nous étudions en premier lieu le couplage entre de telles lois physiques et une hypothèse initiale de parcimonie des origines du phénomène physique. Ceci donne naissance à un concept de dualité des régularisations, formulées soit comme un problème d'analyse coparcimonieuse (menant à la représentation en EDP), soit comme une parcimonie à la synthèse équivalente à la précédente (lorsqu'on fait plutôt usage des fonctions de Green). Nous dédions une part significative de notre travail à la comparaison entre les approches de synthèse et d'analyse. Nous défendons l'idée qu'en dépit de leur équivalence formelle, leurs propriétés computationnelles sont très différentes. En effet, en raison de la parcimonie héritée par la version discrétisée de l'EDP (incarnée par l'opérateur d'analyse), l'approche coparcimonieuse passe bien plus favorablement à l'échelle que le problème équivalent régularisé par parcimonie à la synthèse. Nos constatations sont illustrées dans le cadre de deux applications : la localisation de sources acoustiques, et la localisation de sources de crises épileptiques à partir de signaux électro-encéphalographiques. Dans les deux cas, nous vérifions que l'approche coparcimonieuse démontre de meilleures capacités de passage à l'échelle, au point qu'elle permet même une interpolation complète du champ de pression dans le temps et en trois dimensions. De plus, dans le cas des sources acoustiques, l'optimisation fondée sur le modèle d'analyse \emph{bénéficie} d'une augmentation du nombre de données observées, ce qui débouche sur une accélération du temps de traitement, plus rapide que l'approche de synthèse dans des proportions de plusieurs ordres de grandeur. Nos simulations numériques montrent que les méthodes développées pour les deux applications sont compétitives face à des algorithmes de localisation constituant l'état de l'art. Pour finir, nous présentons deux méthodes fondées sur la parcimonie à l'analyse pour l'estimation aveugle de la célérité du son et de l'impédance acoustique, simultanément à l'interpolation du champ sonore. Ceci constitue une étape importante en direction de la mise en œuvre de nos méthodes en en situation réelle. / Inverse problems related to physical processes are of great importance in practically every field related to signal processing, such as tomography, acoustics, wireless communications, medical and radar imaging, to name only a few. At the same time, many of these problems are quite challenging due to their ill-posed nature. On the other hand, signals originating from physical phenomena are often governed by laws expressible through linear Partial Differential Equations (PDE), or equivalently, integral equations and the associated Green’s functions. In addition, these phenomena are usually induced by sparse singularities, appearing as sources or sinks of a vector field. In this thesis we primarily investigate the coupling of such physical laws with a prior assumption on the sparse origin of a physical process. This gives rise to a “dual” regularization concept, formulated either as sparse analysis (cosparse), yielded by a PDE representation, or equivalent sparse synthesis regularization, if the Green’s functions are used instead. We devote a significant part of the thesis to the comparison of these two approaches. We argue that, despite nominal equivalence, their computational properties are very different. Indeed, due to the inherited sparsity of the discretized PDE (embodied in the analysis operator), the analysis approach scales much more favorably than the equivalent problem regularized by the synthesis approach. Our findings are demonstrated on two applications: acoustic source localization and epileptic source localization in electroencephalography. In both cases, we verify that cosparse approach exhibits superior scalability, even allowing for full (time domain) wavefield interpolation in three spatial dimensions. Moreover, in the acoustic setting, the analysis-based optimization benefits from the increased amount of observation data, resulting in a speedup in processing time that is orders of magnitude faster than the synthesis approach. Numerical simulations show that the developed methods in both applications are competitive to state-of-the-art localization algorithms in their corresponding areas. Finally, we present two sparse analysis methods for blind estimation of the speed of sound and acoustic impedance, simultaneously with wavefield interpolation. This is an important step toward practical implementation, where most physical parameters are unknown beforehand. The versatility of the approach is demonstrated on the “hearing behind walls” scenario, in which the traditional localization methods necessarily fail. Additionally, by means of a novel algorithmic framework, we challenge the audio declipping problemregularized by sparsity or cosparsity. Our method is highly competitive against stateof-the-art, and, in the cosparse setting, allows for an efficient (even real-time) implementation.

Cette thèse porte sur la diffraction d’une onde plane électromagnétique par une surface lisse parfaitement conductrice (PEC). Elle présente l’analyse des propriétés d’une nouvelle formulation des trois principales équations intégrales de frontières de la théorie de la diffraction électromagnétique (EFIE, MFIE et CFIE). L’idée est d’adapter les équations intégrales conventionnelles à la diffraction d’une onde plane en supposant que la fonction de phase de l’onde plane incidente détermine la fonction de phase de la distribution de courant induit sur la surface.L’idée d’utiliser la phase dans la diffraction d’ondes planes a déjà été étudiée pour les hautes fréquences, notamment dans les thèses de Zhou (1995) et Darrigrand (2002) qui adaptèrent les espaces d’approximation des éléments finis. Dans cette thèse, cependant, nous suivons une formulation plus récente, donnée par Herberthson (2008), où la fonction de phase est incorporée dans la distribution du noyau des opérateurs intégraux.En présentant les versions modifiées de l’EFIE et de la MFIE (dénommées HEFIE et HMFIE)dans des espaces fonctionnels appropriés, nous prouvons ici l’existence d’une solution unique à cette formulation spécifique et présentons une mise en oeuvre pratique originale qui tire parti de l’expérience acquise sur l’EFIE/MFIE. Par la suite, nous explorons une propriété importante offerte par ces nouvelles formulations: la possibilité de réduire le nombre de degrés de liberté requis pour obtenir une solution précise du problème. / This thesis is about the scattering of an electromagnetic plane wave incidenton a perfectly conducting smooth surface. It presents the analysis of the properties of a newformulation of the three principal boundary integral equations of electromagnetic scattering theory(EFIE, MFIE and CFIE). The basic idea is to adapt the conventional integral equations toplane-wave scattering by supposing that the phase function of an incident plane wave determinesthe phase function of the induced boundary current distribution.This idea of using the phase in plane wave scattering has previously been studied in highfrequencyscattering, in particular in the theses by Zhou (1995) and Darrigrand (2002) whoadapt the finite element approximation spaces. In this thesis, though, we follow a more recentformulation, given by Herberthson (2008), where the phase function is incorporated in the kerneldistribution of the integral operators.Presenting the modified version of the EFIE and the MFIE (denoted HEFIE and HMFIE) inappropriate function spaces, we prove the existence of a unique solution to this specific formulationand developp an original practical implementation which takes advantage of the gainedexperience on the EFIE/MFIE. Then, we explore another important property provided by thenew formulations: the possibility to reduce the number of degrees of freedom required to get anaccurate solution of the problem.

We consider mathematical modeling of thin film flow between two rough surfaces which are in relative motion. For example such flows take place in different kinds of bearings and gears when a lubricant is used to reduce friction and wear between the surfaces. The mathematical foundations of lubrication theory is given by the Navier--Stokes equation, which describes the motion of viscous fluids. In thin domains several approximations are possible which lead to the so called Reynolds equation. This equation is crucial to describe the pressure in the lubricant film. When the pressure is found it is possible to predict vorous important physical quantities such as friction (stresses on the bounding surfaces), load carrying capacity and velocity field. In hydrodynamic lubrication the effect of surface roughness is not negligible, because in practical situations the amplitude of the surface roughness are of the same order as the film thickness. Moreover, a perfectly smooth surface does not exist in reality due to imperfections in the manufacturing process. Therefore, any realistic lubrication model should account for the effects of surface roughness. This implies that the mathematical modeling leads to partial differential equations with coefficients that will oscillate rapidly in space and time. A direct numerical computation is therefore very difficult, since an extremely dense mesh is needed to resolve the oscillations due to the surface roughness. A natural approach is to do some type of averaging. In this PhD thesis we use and develop modern homogenization theory to be able to handle the questions above. Especially, we use, develop and apply the method based on the multiple scale expansions and two-scale convergence. The thesis is based on five papers (A-E), with an appendix to paper A, and an extensive introduction, which puts these publications in a larger context. In Paper A the connection between the Stokes equation and the Reynolds equation is investigated. More precisely, the asymptotic behavior as both the film thickness <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cepsilon" /> and wavelength <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Cmu" /> of the roughness tend to zero is analyzed and described. Three different limit equations are derived. Time-dependent equations of Reynolds type are obtained in all three cases (Stokes roughness, Reynolds roughness and high frequency roughness regime). In paper C we extend the work done in Paper A where we compare the roughness regimes by numeric computations for the stationary case. In paper B we present a mathematical model that takes into account cavitation, surfaces roughness and compressibility of the fluid. We compute the homogenized coefficients in the case of unidirectional roughness.In the paper D we derive a mathematical model of thin film flow between two close rough surfaces, which takes into account cavitation, surface roughness and pressure dependent density. Moreover, we use two-scale convergence to homogenize the model. Finally, in paper E we prove the existence of solutions to a frequently used mathematical model of thin film flow, which takes cavitation into account.

Ce manuscrit présente des résultats d’analyse mathématique autour de deux exemples de problèmes singuliers d’équations aux dérivées partielles issus de la modélisation. I. Diffusion croisée en dynamique des populations. En dynamique des populations, les systèmes de réaction –diffusion croisée modélisent l’évolution de populations d’espèces en compétition avec un effet répulsif entre les individus. Pour ces systèmes fortement couplés, souvent non linéaires, une question aussi fondamentale que l’existence de solutions se révèle extrêmement complexe. Dans ce manuscrit, on introduit une approche basée sur des extensions récentes de lemmes de dualité et sur des méthodes d’entropie. On démontre l’existence de solutions faibles dans un cadre général de systèmes de réaction-diffusion croisée, ainsi que certaines propriétés qualitatives des solutions. II. Équation de Boltzmann en domaine borné. L’équation de Boltzmann, introduite en 1872, modélise la dynamique des gaz raréfiés hors équilibre. Malgré les nombreux résultats autour de la question de l’existence de solutions fortes proches de l’équilibre, très peu concernent le cas d’un domaine borné, situation pourtant fréquente dans les applications. Une raison de la difficulté du problème est l’irruption des singularités le long des trajectoires rasant le bord du domaine. Dans ce manuscrit, on présente une théorie de la régulation de l’équation de Boltzmann en domaine borné. Grâce à l’introduction d’une distance cinétique qui compense les singularités au bord, on montre des résultats de propagation de normes de Sobolev et de propagation C^1 en domaine convexe. En domaine non convexe, on montre un résultat de propagation de régularité BV. / This manuscript presents results of mathematical analysis concerning two singular problems of partial differential equations coming from the modeling. I. Cross-diffusion in Population dynamics. In Population dynamics, reaction-cross diffusion systems model the evolution of the populations of competing species with a repulsive effect between individuals. For these strongly coupled, often non linear systems, a question as basic as the existence of solutions appears to be extremely complex. In this manuscript, we introduce an approach based on the most recent extensions of duality lemmas and on entropy methods. We prove the existence of weak solutions in a general setting of reaction-cross diffusion systems, as well as some qualitative properties of the solutions. II. Boltzmann equation in bounded domains The Boltzmann equation, introduced in 1872, model the evolution of a rarefied gas out of equilibrium. Despite the numerous results concerning the existence of strong solutions close to equilibrium, very few concern the case of bounded domain, though this situation is very useful in applications. A crucial reason of the difficulty of this problem is the formation of a singularity on the trajectories grazing the boundary. In this manuscript, we present a theory of the regularity of the Boltzmann equation in bounded domains. Thanks to the introduction of a kinetic distance which balances the singularity, we obtain results of propagation of Sobolev norms and C^1 propagation in convex domains. In non convex domains, we obtain the propagation of BV regularity.

This thesis deals with problems arising in the study of nonlinear partial differential equations arising from many-body problems. It is divided into two parts: The first part concerns the derivation of a nonlinear diffusion equation from a microscopic stochastic process. We give a new method to show that in the hydrodynamic limit, the particle densities of a one-dimensional zero range process on a periodic lattice converge to the solution of a nonlinear diffusion equation. This method allows for the first time an explicit uniform-in-time bound on the rate of convergence in the hydrodynamic limit. We also discuss how to extend this method to the multi-dimensional case. Furthermore we present an argument, which seems to be new in the context of hydrodynamic limits, how to deduce the convergence of the microscopic entropy and Fisher information towards the corresponding macroscopic quantities from the validity of the hydrodynamic limit and the initial convergence of the entropy. The second part deals with problems arising in the analysis of nonlinear Schrödinger equations of Gross-Pitaevskii type. First, we consider the Cauchy problem for (energy-subcritical) nonlinear Schrödinger equations with sub-quadratic external potentials and an additional angular momentum rotation term. This equation is a well-known model for superfluid quantum gases in rotating traps. We prove global existence (in the energy space) for defocusing nonlinearities without any restriction on the rotation frequency, generalizing earlier results given in the literature. Moreover, we find that the rotation term has a considerable influence in proving finite time blow-up in the focusing case. Finally, a mathematical framework for optimal bilinear control of nonlinear Schrödinger equations arising in the description of Bose-Einstein condensates is presented. The obtained results generalize earlier efforts found in the literature in several aspects. In particular, the cost induced by the physical work load over the control process is taken into account rather then often used L^2- or H^1-norms for the cost of the control action. We prove well-posedness of the problem and existence of an optimal control. In addition, the first order optimality system is rigorously derived. Also a numerical solution method is proposed, which is based on a Newton type iteration, and used to solve several coherent quantum control problems.

Mathematical models displaying travelling waves appear in a variety of domains. These waves are often faced with different kinds of perturbations. In some cases, these perturbations result in propagation failure, also known as wave-blocking. Wave-blocking has been studied in the case of several specific models, often with the help of numerical tools. In this thesis, we will display a technique that uses symmetry and a center manifold reduction to find a criterion which defines regions in parameter space where a wave will be blocked. We focus on waves with low velocity and small symmetry-breaking perturbations, which is where the blocking initiates; the organising center. The range of the tools used makes the technique easily generalizable to higher dimensions. In order to demonstrate this technique, we apply it to the bistable equation. This allows us to do calculations explicitly. As a result, we show that wave-blocking occurs inside a wedge originating from the organising center and derive an expression for this wedge to leading order. We verify our results with some numerical simulations.