A numerical simulation model for predicting residual stresses which arise during the solidification process of pressed glass bulb panel was developed. The solidification of a molten layer of glass between cooled parallel plates was used to model the mechanics of the buildup of residual stresses in the forming process. A thermorheologically simple thermoviscoelastic model was assumed for the material. The finite element method employed was based on the theory of shells as an assembly of flat elements. This approach calculates residual stresses layer by layer like a truly three-dimensional calculation, which is well suited for thin pressed products of complex shape. An experimental comparison was employed to verify the proposed models and methods.

In this paper, high-order Discontinuous Galerkin (DG) method is used to solve the two-dimensional Euler equations. A shock-capturing method based on the artificial viscosity technique is employed to handle physical discontinuities. Numerical tests show that the shocks can be captured within one element even on very coarse grids. The thickness of the shocks is dominated by the local mesh size and the local order of the basis functions. In order to obtain better shock resolution, a straightforwardhp-adaptivity strategy is introduced, which is based on the high-order contribution calculated using hierarchical basis. Numerical results indicate that thehp-adaptivity method is easy to implement and better shock resolution can be obtained with smaller local mesh size and higher local order.