Numerical modelling of two-dimensional perfect gas flows using RKDG method on unstructured meshes

A two-dimensional, unsteady, transonic, irrotational, inviscid flow of a perfect gas with constant specific heats is considered. The analysis involves perturbations from a uniform sonic isentropic flow. The governing perturbation potential equations are derived for various orders of the ratio of the characteristic time associated with a temporal flow disturbance to the time taken by a sonic disturbance to traverse the transonicregime. The case where this ratio is large compared to one is studied in detail. A similarity solution involving an arbitrary function of time is found and it is shown that this solution corresponds to unsteady chimel flows with either stationary or time-varying wall shapes. Numerical computations are presented showing the temporal changes in flow structure as a disturbance dies out exponentially for the following typical nozzle flows: simple accelerating (Meyer) flow and flow with supersonic pockets (Taylor and limiting Taylor flow).

Download Full-text

Two-dimensional numerical modelling of overtaking shock-wave/moving-body interactions

Shock Waves ◽

10.1007/978-3-540-27009-6_99 ◽

2005 ◽

pp. 659-664

Author(s):

C. Law ◽

B. W. Skews

Keyword(s):

Shock Wave ◽

Numerical Modelling ◽

Two Dimensional ◽

Moving Body

Download Full-text

Derivation of Ice Sliding Properties from The Numerical Modelling of Surging Ice Masses

Journal of Glaciology ◽

10.1017/s0022143000030136 ◽

1979 ◽

Vol 23 (89) ◽

pp. 420-421 ◽

Cited By ~ 2

Author(s):

W. F. Budd ◽

B. J. McInnes ◽

I. Smith

Keyword(s):

Shear Stress ◽

Numerical Modelling ◽

Sliding Friction ◽

Two Dimensional ◽

Dimensional Model ◽

Flow Properties ◽

Two Dimensional Model ◽

Two Parameters ◽

Basal Shear ◽

Best Fit

Abstract It is difficult to deduce sliding properties from the numerical modelling of ordinary glaciers because the flow law of ice is still not known well enough to clearly differentiate sliding from internal deformation of the ice. For glaciers undergoing high-speed surges it appears that the majority of the total speed is due to sliding. Furthermore the average basal shear stress of the ice mass is lowered during the surge. This suggests that surging glaciers can be modelled by incorporating a sliding friction law which has the effective friction coefficient decreasing for high velocities. A relation of this type has been found for ice sliding on granite at −0.5°C by Barnes and others (1971) and has also been obtained for rough slabs with ice at the pressure-melting point by Budd and others (1979). A simple two-dimensional model was developed by Budd and McInnes (1974) and Budd (1975), which was found to exhibit the typical periodic surge-like characteristics of real ice masses. Since the sliding-stress relation for the low velocities and stresses was not known, and was not so important for the surges, it was decided to use the condition of gross equilibrium (i.e. that the ice mass as a whole does not accelerate) together with a single-parameter relation for the way in which the friction decreases with stress and velocity to prescribe the basal shear-stress distribution. The low-stress-velocity relation can thus be obtained as a result. This two-dimensional model has now been parameterized to take account of the three-dimensional aspects of real ice masses. A number of ice masses have since been closely matched by the model including three well-known surging ice masses: Lednik Medvezhiy, Variegated Glacier, and Bruarjökull. Since the flow properties of ice are so poorly known—especially for longitudinal stress and strain-rates—the model has been run with two unknown parameters: one a flow-law parameter (η) and the other a sliding parameter (ø). The model is run over a wide range of these two parameters to see if a good match can be made to the real ice masses and if so what the values of the parameters η and ø are for best fit. The matching of the three above ice masses gave very similar values for each of the two parameters η and ø, the value of η being within the range of values expected for the flow properties of temperate ice as determined by laboratory experiments. Using the same values of η and ø it is found that the ordinary glaciers modelled so far do not develop surging but that they could do if the value of ø were increased or if the mass-balance input were sufficiently increased. For Lednik Medvezhiy a detailed analysis of the friction coefficient with velocity was carried out and it was found that the values required for best fit showed a very close agreement to the sliding friction curve of Barnes and others (1971) at −0.5°C. It is concluded that this type of sliding relation can account for the major features of glacier surge phenomena. Finally it is apparent that the numerical modelling technique can be used very effectively to test any large-scale bulk sliding relation by the analysis of real surges of ice masses and in addition can provide further insight into the sliding relation in association with other stresses in the ice mass.

Download Full-text