What is the Nariai lambdavacuum?
Well, I have no idea what New Scientist is on about this time

but I can tell you what the Nariai lambdavacuum is in gtr and be pretty confident that your UCL mathematician must be talking about the same thing. (Do you recall a name? And is that University College London?)
But first: you should know that most real math/physics types have canceled their subscriptions to NS, which a dozen years ago was a fairly respectable publication, due to an unending spate of sensationalized and highly misleading articles in recent years, such as the now notorious article on an alleged device which would violate conservation of momentum if it really operated as its inventor claims; see
this discussion.
OK, enough of that!
Many Lorentzian manifolds are Cartesian products of a signature (-1,1) manifold (some would say signature zero) with a Riemanian two-manifold. Very very few of these are also solutions of the Einstein field equation. The Nariai lambdavacuum is one of the exceptions. Physically speaking it is a model spacetime in which the only source of the gravitational field is a Lambda term, i.e. that mysterious thing which goes by the name of "dark energy". As you probably know, Lorentzian manifolds can be described using many distinct coordinate charts, and when written in terms of different charts, the metric tensor may appear very different. Nonetheless such different ways of writing a solution really are describing the same geometry, and thus physically speaking (in terms of gtr), the same physical scenario or model universe.
So here is one way of writing down the Nariai solution:
<br />
ds^2 = -\sinh(kx)^2 \, dt^2 + dx^2 + \sinh(kz)^2 \, dy^2 + dz^2,<br />
-\infty < t, \, x, \, y, \, z < \infty<br />
On the domain
<br />
-\infty < t,y < \infty, <br />
\; 0 < x < 1/2 \, \operatorname{arccosh}(\sqrt{2}), <br />
\; 0 < z < \infty<br />
we can adopt the frame field
<br />
\vec{e}_1 = <br />
\frac{1}{\sinh(kx)^2} \, \partial_t <br />
+ \frac{\sqrt{2-\cosh(kx)^2}}{\sinh(kx)} \, \partial_x<br />
<br />
\vec{e}_2 = <br />
\frac{\sqrt{2-\cosh(kx)^2}} {\sinh(kx)^2} \, \partial_t<br />
\frac{1}{\sinh(kx)} \, \partial_x<br />
<br />
\vec{e}_3 =<br />
\frac{1}{\sinh(kz)} \, \partial_y<br />
<br />
\vec{e}_4 = \partial_z<br />
Here, \vec{e}_1 is a timelike unit vector field, whose integral curves represent the world lines of a family of observers, while the other three are spacelike unit vector fields and represent the "spatial directions" defining a "local Lorentz frame" for each of these observers, at each event on his world line (in the domain where the frame field is valid). These observers happen to be inertial observers (their world lines have vanishing path curvature, i.e. vanishing acceleration vector, so they are freely falling in the gravitational field modeled by our spacetime geometry). Furthermore, their world lines form a
vorticity-free timelike geodesic congruence, so they define a unique family of spatial hyperslices orthogonal to the world lines, and our frame is nonspinning, meaning that the spatial vectors are not rotating with respect to the spin axes of gyros carried by our observers. Such nonspinning inertial frames are the closest we can come in general relativity to the global Lorentz frames of special relativity (which are also nonspinning inertial).
There is a way of writing down a dual coframe
\sigma^1, \; \sigma^2, \; \sigma^3, \; \sigma^4
which consists of mutually orthogonal unit covector fields and gives an equivalent description of the same family of observers, but I won't bother to write it out here. As an exercise in computing curvature following the method of Cartan (connection one-forms, curvature two-forms), the Nariai lambda vacuum is good practice. The answer is remarkably simple:
<br />
{\Omega^{\hat{1}}}_{\hat{2}} = -k^2 \, \sigma^1 \wedge \sigma^2, \;<br />
{\Omega^{\hat{3}}}_{\hat{4}} = -k^2 \, \sigma^3 \wedge \sigma^4<br />
which gives the Riemann tensor (dropping the hats sometimes used to stress that we are computing "physical components", i.e. components wrt a frame field, not a coordinate basis--- in this post, all components are taken wrt our frame field)
<br />
R_{1212} = k^2, \; R_{3434} = -k^2<br />
In other words, we have here the constant curvature Lorentzian manifold H^{1,1} \times H^2 of "radius" 1/k. (This is analogous to the Riemannian manifold H^2 \times H^2, the Cartesian product of two hyperbolic planes.) From the Riemann tensor we can compute the Ricci tensor and then the Einstein tensor, which turns out to be
<br />
G^{ab} = k^2 \; \operatorname{diag} \left(-1,1,1,1 \right)<br />
which is diagnostic of a "Lambda" term giving constant negative energy and equal but positive isotropic pressure (in any frame!). The Riemann tensor (considered as a linear operator on bivectors) has one double eigenvalue, -a^2, and one quadruple eigenvalue, zero, which is again a reflection of the highly unusual nature of the Nariai lambdavacuum (a product of two two-dimensional manifolds).
Whenever one is confronted with a timelike congruence in gtr (such as the world lines of our nonspinning inertial observers), one should compute its
kinematic decomposition (acceleration vector, expansion tensor, and vorticity vector). In this case we already noted that the acceleration and vorticity vanish and the expansion tensor is very simple:
<br />
{\theta \left[ \vec{e}_1 \right]}_{ab} = \frac{-k \, \cosh(kx)}{\sqrt{2 - \cosh(kx)^2}} \;<br />
\operatorname{diag} \left( 1,0,0 \right)<br />
The tidal tensor is
<br />
{E \left[ \vec{e}_1 \right]}_{ab} = k^2 \, \operatorname{diag} \left( 1,0,0 \right)<br />
(The kinematic decomposition and the Bel decomposition produce
three-dimensional tensors; the brackets are meant to stress that these decompositions are taken with respect to a particular timelike congruence. Such a procedure is analogous to the way that the EM field tensor is decomposed into electric and magnetic vectors, which are
three-dimensional vectors.) This means that our observers are drawing apart along the \vec{e}_2 direction, due to a constant tidal tension along that direction which each of our ideal observers can in principle measure with sensitive strain gauges attached to the skin of his spaceship. The magnetogravitic tensor vanishes (no spin-spin forces on gyroscopes carried by our observers!), and the Riemannian curvature tensor of our spatial hyperslices shows that these manifolds are all isometric to the Cartesian product R^1 \times H^2 where the spherical factor has "radius" 1/k.
I seem to be spending more effort than usual warning various PF posters to be cautious about phrases like "space is expanding". As I keep repeating, almost always what this really means is that some family of observers is expanding in the technical sense of expansion tensor. Indeed, we see very clearly from the above analysis that our inertial observers are drawing apart along a preferred direction. (I left as an exercise the task of determing the Lie algebra of Killing vectors, but anyone who knows what that means can probably write down generators from the information that as a Lorentzian manifold we are dealing here with H^{1,1} \times H^2.)
Without more information (such as a citation of the paper by your UCL mathematician) I cannot answer your question about how pieces of such a Nariai lambdavacuum are supposed to be lurking inside black holes, but I can say this: suggestions that the interior of black holes might contain regions which can considered "phenomonologically" (there's that word again!) as something like pieces of the Nariai lambdavacuum (or pieces of the de Sitter lambdavacuum) are nothing new: this rather vague idea has been kicked around for at least a decade. As I said, New Scientist is known for not being very scrupulous in its reporting
