Oscillons are extremely long lived, oscillatory, spatially localized field configurations that arise from generic initial conditions in a large number of non-linear field theories. With an eye towards their cosmological implications, we investigate their properties in an expanding universe. We (1) provide an analytic solution for one dimensional oscillons (for the models under consideration) and discuss their generalization to 3 dimensions, (2) discuss their stability against long wavelength perturbations and (3) estimate the effects of expansion on their shapes and life-times. In particular, we discuss a new, extended class of oscillons with surprisingly flat tops. We show that these flat topped oscillons are more robust against collapse instabilities in (3+1) dimensions than their usual counterparts. Unlike the solutions found in the small amplitude analysis, the width of these configurations is a non-monotonic function of their amplitudes.Comment: v2-matches version published in Phys. Rev D. Updated references and minor modification to section 4.

We study the thermalization of the classical Klein-Gordon equation under a u 4 interaction. We numerically show that even in the presence of strong nonlinearities, the local thermodynamic equilibrium state exhibits a weakly nonlinear behavior in a renormalized wave basis. The renormalized basis is defined locally in time by a linear transformation and the requirement of vanishing wavewave correlations. We show that the renormalized waves oscillate around one frequency, and that the frequency dispersion relation undergoes a nonlinear shift proportional to the mean square field. In addition, the renormalized waves exhibit a Planck like spectrum. Namely, there is equipartition of energy in the low frequency modes described by a Boltzmann distribution, followed by a linear exponential decay in the high frequency modes.

Common efficient schemes for the incompressible Navier-Stokes equations, such as projection or fractional step methods, have limited temporal accuracy as a result of matrix splitting errors, or introduce errors near the domain boundaries (which destroy uniform convergence to the solution). In this paper we recast the incompressible (constant density) Navier-Stokes equations (with the velocity prescribed at the boundary) as an equivalent system, for the primary variables velocity and pressure. We do this in the usual way away from the boundaries, by replacing the incompressibility condition on the velocity by a Poisson equation for the pressure. The key difference from the usual approaches occurs at the boundaries, where we use boundary conditions that unequivocally allow the pressure to be recovered from knowledge of the velocity at any fixed time. This avoids the common difficulty of an, apparently, over-determined Poisson problem. Since in this alternative formulation the pressure can be accurately and efficiently recovered from the velocity, the recast equations are ideal for numerical marching methods. The new system can be discretized using a variety of methods, including semiimplicit treatments of viscosity, and in principle to any desired order of accuracy. In this work we illustrate the approach with a 2-D second order finite difference scheme on a Cartesian grid, and devise an algorithm to solve the equations on domains with curved (non-conforming) boundaries, including a case with a non-trivial topology (a circular obstruction inside the domain).This algorithm achieves second order accuracy in the L ∞ norm, for both the velocity and the pressure. The scheme has a natural extension to 3-D.

The volume penalty method provides a simple, efficient approach for solving the incompressible Navier-Stokes equations in domains with boundaries or in the presence of moving objects. Despite the simplicity, the method is typically limited to first order spatial accuracy. We demonstrate that one may achieve high order accuracy by introducing an active penalty term. One key difference from other works is that we use a sharp, unregularized mask function. We discuss how to construct the active penalty term, and provide numerical examples, in dimensions one and two. We demonstrate second and third order convergence for the heat equation, and second order convergence for the Navier-Stokes equations. In addition, we show that modifying the penalty term does not significantly alter the time step restriction from that of the conventional penalty method.

This paper presents a new class of high order linear ImEx multistep schemes with large regions of unconditional stability. Unconditional stability is a desirable property of a time stepping scheme, as it allows the choice of time step solely based on accuracy considerations. Of particular interest are problems for which both the implicit and explicit parts of the ImEx splitting are stiff. Such splittings can arise, for example, in variable-coefficient problems, or the incompressible Navier-Stokes equations. To characterize the new ImEx schemes, an unconditional stability region is introduced, which plays a role analogous to that of the stability region in conventional multistep methods. Moreover, computable quantities (such as a numerical range) are provided that guarantee an unconditionally stable scheme for a proposed implicit-explicit matrix splitting. The new approach is illustrated with several examples. Coefficients of the new schemes up to fifth order are provided.