Changes to NDSolve from MMA10.0 to later versions

Questions
1

Hi,

some time ago I posted a question on “Intersection boundary problems when solving PDEs” . It appears that in MMA10.0 I can run the following code without any problems (?):

mr = DiscretizeGraphics[ExampleData[{"Geometry3D", "SpaceShuttle"}]];

uif = NDSolveValue[{Inactive[Laplacian][u[x, y, z], {x, y, z}] == 0, 
    DirichletCondition[u[x, y, z] == 1, z <= -1.3]}, 
   u, {x, y, z} \[Element] mr];

Needs["NDSolve`FEM`"]
ElementMeshSurfacePlot3D[uif, Boxed -> False, ViewPoint -> {0, -4, 2}]

It gives:

ScreenShot2015-11-11at13.46.44.png
ScreenShot2015-11-11at13.46.44.png704×430 43.5 KB

If I run it in 10.1, 10.2 or 10.3 I get:

ScreenShot2015-11-11at13.47.49.png
ScreenShot2015-11-11at13.47.49.png1420×608 68.6 KB

Does anyone know what has happened?

Cheers,

Marco

2

Sorry, a bit unprofessional, but I just noticed that :

mr = BoundaryDiscretizeGraphics[ExampleData[{"Geometry3D", "SpaceShuttle"}]];
uif = NDSolveValue[{Inactive[Laplacian][u[x, y, z], {x, y, z}] == 0, 
    DirichletCondition[u[x, y, z] == 1, z <= -1.3]}, 
   u, {x, y, z} \[Element] mr];

Needs["NDSolve`FEM`"]
ElementMeshSurfacePlot3D[uif, Boxed -> False, ViewPoint -> {0, -4, 2}]

does work in 10.3.

Cheers,

Marco