On weak discontinuities through thermally conducting and dissociating gases

Using singular surface theory, the phenomena associated with the uniform and nonuniform propagation of weak discontinuities through thermally conducting and dissociating gases is studied. The basic differential equations governing the criteria for decay or blow up of these discontinuities is obtained. It turns out that growth and decay of weak discontinuities are derived and solved completely. The the thermal conduction and dissociation allow the existence of a singular surface carrying a weak discontinuity which grows into a shock and the role of dissociation and thermal conduction is to cause rapid damping in the formation of this shock.     

Propagation of weak discontinuities in binary mixtures of ideal gases

The propagation of the weak discontinuities in binary non-reacting mixtures of classical ideal monoatomic gases is analyzed. The normal speeds of propagation are determined and compared with those of a single fluid. The differential equation governing the growth and the decay of the acceleration waves is obtained and the solutions for plane, cylindrical and spherical waves are shown. The influence of the different atomic masses of the constituents is also investigated.     

The propagation of weak discontinuities in quasi-linear symmetric hyperbolic systems

Zusammenfassung Die vorliegende Abhandlung untersucht die Fortpflanzung kleiner Unstetigkeiten in Systemen von nichtlinearen hyperbolischen Differentialgleichungen. Ein Ausdruck wird abgeleitet, der die Änderung in der Intensität der Unstetigkeit angibt, wenn diese sich entlang eines Strahls des hyperbolischen Gleichungssystems fortbewegt. Schliesslich wird als Beispiel mit Hilfe des angegebenen Verfahrens die Fortpflanzung von Schallwellen behandelt.     

On wave propagation in elastic tubes conducting rotational flows

The propagation of perturbations in liquid filled elastic tubes depends on the stream velocity of the basic flow. This phenomenon is currently analyzed with the method of the characteristics which relies upon a basic flow with a rectangular velocity profile. It seems that this one-dimensional flow approximation has not been convincingly validated, which justifies to consider other, more general velocity profiles.In the present analytical study the velocity profile is a quadratic function of the radial coordinate. Small amplitude perturbations are superposed on this inviscid, basic state in which the mean velocity is arbitrarily large. A normal mode analysis shows that the velocity profile and therefore the vorticity of the basic flow influence the more the phenomenon the larger is. For example, a parabolic profile allows countercurrent wave propagation regardless of.This questions the one-dimensional wave propagation theory in compliant tubes and, consequently, the interpretation of several physiological and medical problems mainly in the respiratory and cardio-vascular systems.

---

Resumé La propagation de perturbations dans un tube élastique conduisant un écoulement fluide dépend de la vitesse de l'écoulement de base. Ce phénomène est habituellement étudié avec la méthode des caractéristiques, où l'on suppose que le profil de vitesse de l'écoulement est rectangulaire. Comme cette simplification ne semble pas avoir été bien validée, il paraît indiqué d'étudier l'impact d'autres profils.Dans la présente étude analytique, ce profil de vitesse est une fonction quadratique de la coordonée radiale. A cet écoulement non visqueux, dont la vitesse moyenne est arbitraire, l'on superpose des perturbations de faible amplitude. Une analyse linéarisée montre que le profil de vitesse et donc le rotationel de l'écoulement de base influencent d'autant plus ce phénomène d'ondes que est élévée.Ceci met en question la théorie uni-dimensionelle de la propagation d'ondes dans des tubes compliants et, par là-même, l'interpretation de divers problèmes physiologiques et médicaux, avant tout des systèmes respiratoires et cardio-vasculaires.

     

Growth of weak discontinuities in arbitrarily moving gas

The theory of singular surfaces is combined with the ray theory to the study of anisotropic, non-linear wave-propagation in an arbitrarily moving gas. The governing equation for the strength of the wave along the rays is expressed in an integral form. Use of the analysis is made by working out two examples in detail.     

Evolution of weak discontinuities in shallow water equations

In this paper, we determine the critical time, when a weak discontinuity in the shallow water equations culminates into a bore. Invariance group properties of the governing system of partial differential equations (PDEs), admitting Lie group of point transformations with commuting infinitesimal operators, are presented. Some appropriate canonical variables are characterized that transform equations at hand to an equivalent form, which admits non-constant solutions. The propagation of weak discontinuities is studied in the medium characterized by the particular solution of the governing system.     

On the breakdown of acceleration waves in dissociating gas flows

The present paper is devoted to the study of characteristic solution in the neighbourhood of the leading frozen characteristics in dissociating gas flows. It is found that at the cusp of the envelope of intersecting forward characteristics there occurs a breakdown of the wave after a finite critical timet _o. It is observed that there exists a critical value of the initial amplitude of the wave such that all compressive waves with an initial amplitude greater than the critical one will terminate into a shock wave due to non-linear steepening while an initial amplitude less than the critical one will result in a continuous decay. It is also concluded that the breakdown point moves forward along the leading characteristics due to dissociation effects.     

Evolution of weak discontinuities in a non-ideal radiating gas

A method of wavefront analysis is used to study the formation of shock waves in a two-dimensional steady supersonic flow of a non-ideal radiating gas past plane and axisymmetric bodies. The gas is taken to be sufficiently hot for the effect of thermal radiation to be significant, which is, of course, treated by the optically thin approximation to the radiative transfer equation. Transport equations, which lead to the determination of the shock formation distance and also to conditions which insure that no shock will ever evolve on the wavefront, is derived. The influence of the parameter of the non-idealness, upstream flow Mach number in the presence of thermal radiation on the behavior of the wavefront are examined.     

The propagation of weak concentration discontinuities in the isothermal flow of a multicomponent mixture in a porous medium with phase transitions

It is shown that for the general case of a system of non-linear equations, describing multicomponent isothermal flow in a porous medium with phase transitions, as in hyperbolic systems, weak concentration discontinuities propagate with finite velocities, which are determined by solving an eigenvalue problem. If the seeping phases are incompressible and there are no phase transitions, the results obtained for weak discontinuities transfer into the well-known formulae for the Buckley – Leverett model. The results are demonstrated for the case of two-component seepage with phase transitions.     

Anisotropic weak Hardy spaces and interpolation theorems

In this paper, the authors establish the anisotropic weak Hardy spaces associated with very general discrete groups of dilations. Moreover, the atomic decomposition theorem of the anisotropic weak Hardy spaces is also given. As some applications of the above results, the authors prove some interpolation theorems and obtain the boundedness of the singular integral operators on these Hardy spaces.     

Real gas effects in thermally choked nozzle flows

Real gas effects in condensing nozzle flows are discussed by the virial equation of state truncated after the second virial coefficient. The thermal choking conditions in nozzles previously derived for a perfect condensible vapor are generalized to include real gas effects. For these cases it is shown that the critical amount of heat necessary to thermally choke the flow can be defined explicitly only for the expansion of a pure vapor.Alexander von Humboldt Fellow.The flow Mach number is usually taken as the local frozen Mach number.     

Simulation of crack–propagation using embedded discontinuities

The transition from microscale damage phenomena to crack initiation and growth at the macroscale is an important mechanism which constrains the lifetime of concrete structures. Analysing crack growth using the finite element method without enhancement of the shape functions is possible only by continuously updating the corresponding meshes, which constitutes a significant computational effort. But even then the results can be substantially mesh–dependent and hard to interpret. The extended Finite Element Method (XFEM) uses additional discontinuous shape–functions and is one possibility to overcome these problems. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Well-posedness of transonic characteristic discontinuities in two-dimensional steady compressible Euler flows

In our previous work, we have established the existence of transonic characteristic discontinuities separating supersonic flows from a static gas in two-dimensional steady compressible Euler flows under a perturbation with small total variation of the incoming supersonic flow over a solid right wedge. It is a free boundary problem in Eulerian coordinates and, across the free boundary (characteristic discontinuity), the Euler equations are of elliptic–hyperbolic composite-mixed type. In this paper, we further prove that such a transonic characteristic discontinuity solution is unique and L ¹–stable with respect to the small perturbation of the incoming supersonic flow in Lagrangian coordinates.     

Simulation of shallow water flows with shoaling areas and bottom discontinuities

A numerical method based on a second-order accurate Godunov-type scheme is described for solving the shallow water equations on unstructured triangular-quadrilateral meshes. The bottom surface is represented by a piecewise linear approximation with discontinuities, and a new approximate Riemann solver is used to treat the bottom jump. Flows with a dry sloping bottom are computed using a simplified method that admits negative depths and preserves the liquid mass and the equilibrium state. The accuracy and performance of the approach proposed for shallow water flow simulation are illustrated by computing one- and two-dimensional problems.     

Wave propagation and breaking in multidimensional magnetogasdynamic flows

Summary The theory of shocks evolving from weak discontinuities is applied to magnetogasdynamics and worked out in detail for two situations: a spherical disturbance propagating into a uniform magnetic field, and a planar disturbance propagating in the field of a current-carrying wire. For both cases, shock formation is delayed as the magnetoacoustic speed of propagation approaches the Alfvén speed. This agrees with the fact that Alfvén waves do not steepen to form shocks.

---

Résumé On applique à la magnétogasdynamique la théorie des ondes de choc produites par l'évolution des discontinuités faibles. On présente des résultats sur la formation de chocs par des ondes magnétoacoustiques. Lorsque la vitesse de propagation de ces ondes se rapproche de la vitesse d'Alfvén, l'apparition des chocs est retardée, en accord avec le résultat bien connu que les ondes d'Alfvén ne se transforment pas pour former des chocs.

---

---

Research supported in part by Grant No. AFOSR-73-2561.     

A difference scheme with anti-diffusion terms to improve the computation of contact discontinuities and the study of weak discontinuities in compressible flow

The investigation of Mach reflection formed after the impingement of a weak plane shock wave on a wedge with shock Mach number M_s near 1, is still an open problem[12]. It's difficult for shock tube experiments with interferometer to detect contact discontinuities if it is too weak; also difficult to catch with due accuracy the transition condition between Mach reflection and regular reflection. The interest to this phenomenon is continuing, especially for weak shocks, because there was systematic discrepancy between simplified three shock theory of von Neumann [8] and shock tube results [15] which was named by G. Birkhoff as “von Neumann Paradox on three shock theory” [18].In 1972, K.O.Friedrichs called for more computational efforts on this problem. Recently it is known that for weak impinging shocks it's still difficult to get contact discontinuities and curved Mach stem with satisfactory accuracy. Recent numerical computation sometimes even fails to show reflected shock wave[6]. These explain why von Neumann paradox of the three shock theory in case of weak discontinuities is still a problem of interesting [9,12,14]. In this paper, on one hand, we investigate the numerical methods for Euler's equation for compressible inviscid flow, aiming at improving the computation of contact discontinuities, on the other hand, a methodology is suggested to correctly plot flow data from the massive information in storage. On this basis, all the reflected shock wave , contact discontinuities and the curved Mach stem are determined. We get Mach reflection under the condition when over-simplified shock theory predicts no such configuration[5].     

Anisotropic weak Hardy spaces of Musielak-Orlicz type and their applications

Anisotropy is a common attribute of the nature, which shows different characterizations in different directions of all or part of the physical or chemical properties of an object. The anisotropic property, in mathematics, can be expressed by a fairly general discrete group of dilations {A ^k : k ∈ ?}, where A is a real n × n matrix with all its eigenvalues λ satisfy |λ| > 1. The aim of this article is to study a general class of anisotropic function spaces, some properties and applications of these spaces. Let φ: ? ⁿ ×[0,∞) → [0,∞) be an anisotropic p-growth function with p ∈ (0, 1]. The purpose of this article is to find an appropriate general space which includes weak Hardy space of Fefferman and Soria, weighted weak Hardy space of Quek and Yang, and anisotropic weak Hardy space of Ding and Lan. For this reason, we introduce the anisotropic weak Hardy space of Musielak-Orlicz type H _A ^φ,∞ (? ⁿ ) and obtain its atomic characterization. As applications, we further obtain an interpolation theorem adapted to H _A ^φ,∞ (? ⁿ ) and the boundedness of the anisotropic Calderón-Zygmund operator from H _A ^φ,∞ (? ⁿ ) to L _A ^φ,∞ (? ⁿ ). It is worth mentioning that the superposition principle adapted to the weak Musielak-Orlicz function space, which is an extension of a result of E. M. Stein, M. Taibleson and G. Weiss, plays an important role in the proofs of the atomic decomposition of H _A ^φ,∞ (? ⁿ ) and the interpolation theorem.     

Finite element analysis of thermally coupled nonlinear Darcy flows

We consider a coupled system describing nonlinear Darcy flows with temperature dependent viscosity and with viscous heating. We first establish existence, uniqueness, and regularity of the weak solution of the system of equations. Next, we decouple the coupled system by a fixed point algorithm and propose its finite element approximation. Finally, we present convergence analysis with an error estimate between continuous solution and its iterative finite element approximation.© 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010