On a class of exact solutions to the Fokker-Planck equations

In this paper we study under which circumstances there exists a general change of gross variables that transforms any FokkerPlanck equation into another of the OrnsteinUhlenbeck class that, therefore, has an exact solution. We find that any FokkerPlanck equation will be exactly solvable by means of a change of gross variables if and only if the curvature tensor and the torsion tensor associated with the diffusion is zero and the transformed drift is linear. We apply our criteria to the Kubo and Gompertz models.

Classical solutions of the leading-logarithm approximation with non-trivial topology

Exact solutions of the classical equations corresponding to the leading-logarithm approximation are obtained. They are classified by an (integer) topological number.

Exact temporal evolution for some non-linear diffusion process

Exact solutions to FokkerPlanck equations with nonlinear drift are considered. Applications of these exact solutions for concrete models are studied. We arrive at the conclusion that for certain drifts we obtain divergent moments (and infinite relaxation time) if the diffusion process can be extended without any obstacle to the whole space. But if we introduce a potential barrier that limits the diffusion process, moments converge with a finite relaxation time.

Application of standard and refined heat balance integral methods to one-dimensional Stefan problems

The work in this paper concerns the study of conventional and refined heat balance integral methods for a number of phase change problems. These include standard test problems, both with one and two phase changes, which have exact solutions to enable us to test the accuracy of the approximate solutions. We also consider situations where no analytical solution is available and compare these to numerical solutions. It is popular to use a quadratic profile as an approximation of the temperature, but we show that a cubic profile, seldom considered in the literature, is far more accurate in most circumstances. In addition, the refined integral method can give greater improvement still and we develop a variation on this method which turns out to be optimal in some cases. We assess which integral method is better for various problems, showing that it is largely dependent on the specified boundary conditions.

Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

An approximate solution method for boundary layer flow of a power law fluid over a flat plate

The work in this paper deals with the development of momentum and thermal boundary layers when a power law fluid flows over a flat plate. At the plate we impose either constant temperature, constant flux or a Newton cooling condition. The problem is analysed using similarity solutions, integral momentum and energy equations and an approximation technique which is a form of the Heat Balance Integral Method. The fluid properties are assumed to be independent of temperature, hence the momentum equation uncouples from the thermal problem. We first derive the similarity equations for the velocity and present exact solutions for the case where the power law index n = 2. The similarity solutions are used to validate the new approximation method. This new technique is then applied to the thermal boundary layer, where a similarity solution can only be obtained for the case n = 1.

Extensional flow of nematic liquid crystal under electric field gradient

Systematic asymptotic methods are used to formulate a model for the extensional flow of a thin sheet of nematic liquid crystal. With no external body forces applied, the model is found to be equivalent to the so-called Trouton model for Newtonian sheets (and fi bers), albeit with a modi fied "Trouton ratio". However, with a symmetry-breaking electric field gradient applied, behavior deviates from the Newtonian case, and the sheet can undergo fi nite-time breakup if a suitable destabilizing field is applied. Some simple exact solutions are presented to illustrate the results in certain idealized limits, as well as sample numerical results to the full model equations.

Systematic weakly nonlinear analysis of radial viscous fingering

We present a weakly nonlinear analysis of the interface dynamics in a radial Hele-Shaw cell driven by both injection and rotation. We extend the systematic expansion introduced in [E. Alvarez-Lacalle et al., Phys. Rev. E 64, 016302 (2001)] to the radial geometry, and compute explicitly the first nonlinear contributions. We also find the necessary and sufficient condition for the uniform convergence of the nonlinear expansion. Within this region of convergence, the analytical predictions at low orders are compared satisfactorily to exact solutions and numerical integration of the problem. This is particularly remarkable in configurations (with no counterpart in the channel geometry) for which the interplay between injection and rotation allows that condition to be satisfied at all times. In the case of the purely centrifugal forcing we demonstrate that nonlinear couplings make the interface more unstable for lower viscosity contrast between the fluids.

Interface dynamics in Hele-Shaw flows with centrifugal forces. Preventing cusp singularities with rotation

A class of exact solutions of Hele-Shaw flows without surface tension in a rotating cell is reported. We show that the interplay between injection and rotation modifies the scenario of formation of finite-time cusp singularities. For a subclass of solutions, we show that, for any given initial condition, there exists a critical rotation rate above which cusp formation is suppressed. We also find an exact sufficient condition to avoid cusps simultaneously for all initial conditions within the above subclass.

Dynamical Systems approach to Saffman-Taylor fingering. A Dynamical Solvability Scenario

A dynamical systems approach to competition of Saffman-Taylor fingers in a Hele-Shaw channel is developed. This is based on global analysis of the phase space flow of the low-dimensional ordinary-differential-equation sets associated with the classes of exact solutions of the problem without surface tension. Some simple examples are studied in detail. A general proof of the existence of finite-time singularities for broad classes of solutions is given. Solutions leading to finite-time interface pinchoff are also identified. The existence of a continuum of multifinger fixed points and its dynamical implications are discussed. We conclude that exact zero-surface tension solutions taken in a global sense as families of trajectories in phase space are unphysical because the multifinger fixed points are nonhyperbolic, and an unfolding does not exist within the same class of solutions. Hyperbolicity (saddle-point structure) of the multifinger fixed points is argued to be essential to the physically correct qualitative description of finger competition. The restoring of hyperbolicity by surface tension is proposed as the key point to formulate a generic dynamical solvability scenario for interfacial pattern selection.

Surface tension and dynamics of fingering patterns

We study the minimal class of exact solutions of the Saffman-Taylor problem with zero surface tension, which contains the physical fixed points of the regularized (nonzero surface tension) problem. New fixed points are found and the basin of attraction of the Saffman-Taylor finger is determined within that class. Specific features of the physics of finger competition are identified and quantitatively defined, which are absent in the zero surface tension case. This has dramatic consequences for the long-time asymptotics, revealing a fundamental role of surface tension in the dynamics of the problem. A multifinger extension of microscopic solvability theory is proposed to elucidate the interplay between finger widths, screening and surface tension.

Black diholes

We present and analyze exact solutions of the Einstein-Maxwell and Einstein-Maxwell-dilaton equations that describe static pairs of oppositely charged extremal black holes, i.e., black diholes. The holes are suspended in equilibrium in an external magnetic field, or held apart by cosmic strings. We comment as well on the relation of these solutions to brane-antibrane configurations in string and M theory.

Solution to telegrapher's equation in the presence of reflecting and partly reflecting broundaries

We show that the reflecting boundary condition for a one-dimensional telegraphers equation is the same as that for the diffusion equation, in contrast to what is found for the absorbing boundary condition. The radiation boundary condition is found to have a quite complicated form. We also obtain exact solutions of the telegraphers equation in the presence of these boundaries.

Higher-dimensional rotating charged black holes

Using the blackfold approach, we study new classes of higher-dimensional rotating black holes with electric charges and string dipoles, in theories of gravity coupled to a 2-form or 3-form field strength and to a dilaton with arbitrary coupling. The method allows to describe not only black holes with large angular momenta, but also other regimes that include charged black holes near extremality with slow rotation. We construct explicit examples of electric rotating black holes of dilatonic and non-dilatonic Einstein-Maxwell theory, with horizons of spherical and non-spherical topology. We also find new families of solutions with string dipoles, including a new class of prolate black rings. Whenever there are exact solutions that we can compare to, their properties in the appropriate regime are reproduced precisely by our solutions. The analysis of blackfolds with string charges requires the formulation of the dynamics of anisotropic fluids with conserved string-number currents, which is new, and is carried out in detail for perfect fluids. Finally, our results indicate new instabilities of near-extremal, slowly rotating charged black holes, and motivate conjectures about topological constraints on dipole hair.