Shock waves in reactive hydrodynamics

Using the weakly non-linear geometrical acoustics theory, we obtain the small amplitude high frequency asymptotic solution to the basic equations in Eulerian coordinates governing one dimensional unsteady planar, spherically and cylindrically symmetric flow in a reactive hydrodynamic medium. We derive the transport equations for the amplitudes of resonantly interacting waves. The evolutionary behavior of non-resonant wave modes culminating into shock waves is also studied.      

Global Uniqueness of Transonic Shocks in Two-Dimensional Steady Compressible Euler Flows

We prove that for the two-dimensional steady complete compressible Euler system, with given uniform upcoming supersonic flows, the following three fundamental flow patterns (special solutions) in gas dynamics involving transonic shocks are all unique in the class of piecewise C ¹ smooth functions, under appropriate conditions on the downstream subsonic flows: (i) the normal transonic shocks in a straight duct with finite or infinite length, after fixing a point the shock-front passing through; (ii) the oblique transonic shocks attached to an infinite wedge; (iii) a flat Mach configuration containing one supersonic shock, two transonic shocks, and a contact discontinuity, after fixing a point where the four discontinuities intersect. These special solutions are constructed traditionally under the assumption that they are piecewise constant, and they have played important roles in the studies of mathematical gas dynamics. Our results show that the assumption of a piecewise constant can be replaced by some weaker assumptions on the downstream subsonic flows, which are sufficient to uniquely determine these special solutions. Mathematically, these are uniqueness results on solutions of free boundary problems of a quasi-linear system of elliptic-hyperbolic composite-mixed type in bounded or unbounded planar domains, without any assumptions on smallness. The proof relies on an elliptic system of pressure p and the tangent of the flow angle w = v/u obtained by decomposition of the Euler system in Lagrangian coordinates, and a newly developed method for the L ^∞ estimate that is independent of the free boundaries, by combining the maximum principles of elliptic equations, and careful analysis of the shock polar applied on the (maybe curved) shock-fronts.     

Pointwise Estimates and Stability for Dispersive–Diffusive Shock Waves

We study the stability and pointwise behavior of perturbed viscous shock waves for a general scalar conservation law with constant diffusion and dispersion. Along with the usual Lax shocks, such equations are known to admit undercompressive shocks. We unify the treatment of these two cases by introducing a new wave-tracking method based on “instantaneous projection”, giving improved estimates even in the Lax case. Another important feature connected with the introduction of dispersion is the treatment of a non-sectorial operator. An immediate consequence of our pointwise estimates is a simple spectral criterion for stability in all L ^p norms, p≥ 1 for the Lax case and p > 1 for the undercompressive case. Our approach extends immediately to the case of certain scalar equations of higher order, and would also appear suitable for extension to systems. Accepted May 29, 2000?Published online November 16, 2000     

Axisymmetric simulation of the interaction of a rising bubble with a rigid surface in viscous flow

The interaction between a rising deformable gas bubble and a solid wall in viscous liquids is investigated by direct numerical simulation via an arbitrary-Lagrangian–Eulerian (ALE) approach. The flow field is assumed to be axisymmetric. The bubble is driven by gravity only and the motion of the gas inside the bubble is neglected. Deformation of the bubble is tracked by a moving triangular mesh and the liquid motion is obtained by solving the Navier–Stokes equations in a finite element framework. To understand the mechanisms of bubble deformation as it interacts with the wall, the interaction process is studied as a function of two dimensionless parameters, namely, the Morton number (Mo) and Bond number (Bo). We study the range of Bo and Mo from (2, 6.5 × 10⁻⁶) to (16, 0.1). The film drainage process is also considered in this study. It is shown that the deformation of a bubble interacting with a solid wall can be classified into three modes depending on the values of Mo and Bo.     

A mode parabolic equation method in the case of the resonant mode interaction

A mode parabolic equation method for resonantly interacting modes has been developed. The method ensures acoustic energy flux conservation in accord with multiple scaled analysis accuracy. We carried out testing calculations for the ASA wedge benchmark, that showed excellent agreement with the COUPLE program.     

Mixed convection boundary layer flow over a horizontal circular cylinder with Newtonian heating

The steady mixed convection boundary layer flow over a horizontal circular cylinder, generated by Newtonian heating in which the heat transfer from the surface is proportional to the local surface temperature, is considered in this study. The governing boundary layer equations are first transformed into a system of non-dimensional equations via the non-dimensional variables, and then into non-similar equations before they are solved numerically using a numerical scheme known as the Keller-box method. Numerical solutions are obtained for the skin friction coefficient Re ^1/2 C _f and the local wall temperature θ _w (x) as well as the velocity and temperature profiles with two parameters, namely the mixed convection parameter λ and the Prandtl number Pr.     

Equations of Langmuir turbulence and Zakharov equations: smoothness and approximation

A family of systems parameterized by H > 0,which describes the Langmuir turbulence,is considered.The asymptotic behavior of the solutions(E H,n H) when H goes to zero is studied.The results of convergence of(E H,n H) to the couple(E,n) which is the solution to the Zakharov equations are stated.     

Useful zero friction simulations for assessing MBS codes Pascal's formula giving wheelsets frequency for zero wheel-rail friction

Present study provides a simple analytical formula, the “Klingel-like formula” or “Pascal's Formula” that can be used as a reference to test some results of existing railway codes and specifically those using rigid contact. It develops properly the 3D Newton-Euler equations governing the 6 degrees of freedom (DoF) of unsuspended loaded wheelsets in case of zero wheel-rail friction and constant conicity. Thus, by solving numerically these equations, we got pendulum like harmonic oscillations of which the calculated angular frequency is used for assessing the accuracy of the proposed formula so that it can in turn be used as a fast practical target for testing multi-body system (MBS) railway codes. Due to the harmonic property of these pendulum-like oscillations, the square ω² of their angular frequency can be made in the form of a ratio K/M where K depends on the wheelset geometry and load and M on its inertia. Information on K and M are useful to understand wheelsets behavior. The analytical formula is derived from the first order writing of full trigonometric Newton-Euler equations by setting zero elastic wheel-rail penetration and by assuming small displacements. Full trigonometric equations are numerically solved to assess that the formula provides ω² inside a 1% accuracy for usual wheelsets dimensions. By decreasing the conicity down to 1 × 10^?4 rad, the relative formula accuracy is under 3 × 10^?5. In order to test the formula reliability for rigid contact formulations, the stiffness of elastic contacts can be increased up to practical rigidity (Hertz stiffness × 1000).     

Shocks in local supersonic zones

The results of numerical integration of the Euler equations governing two-dimensional and axisymmetric flows of an ideal (inviscid and non-heat-conducting) gas with local supersonic zones are presented. The subject of the study is the formation of shocks closing local supersonic zones. The flow in the vicinity of the initial point of the closing shock is calculated on embedded, successively refined grids with an accuracy much greater than that previously achieved. The calculations performed, together with the analysis of certain controversial issues, leave no doubt that it is the intersection of C ^?-characteristics proceeding from the sonic line inside the supersonic zone that is responsible for the closing shock formation.     

Nonlinear resonances and wave jumps in media with higher-order dispersion

A new technique for systematically investigating biperiodic (two-wave) steady-state solutions is described with reference to modified Korteweg-de Vries and Schrödinger equations which generalize the conventional model equations for waves on water, in plasmas, and in nonlinear optics [1]. Among these solutions those with ordinary and resonance wave interactions are distinguished. Both singular solutions similar to the solitons of a resonantly interacting wave envelope and solitary waves are found. The soliton-like solutions obtained are used for describing the wave jump structure.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 113–124, July–August, 1996.     

Similarity solution of mixed convection over a horizontal moving plate

Investigation to the mixed convective heat and mass transfer over a horizontal plate has been carried out. By applying transformation group theory to analysis of the governing equations of continuity, momentum, energy and diffusion, we show the existence of similarity solution for the problem provided that the temperature and concentration at the wall are proportional to x ^4/(7-5n) and that the moving speed of the plate is proportional to x ^(3-n)/(7-5n), and further obtain a similarity representation of the problem. The similarity equations have been solved numerically by a fourth-order Runge–Kutta scheme. The numerical results obtained for Pr=0.72 and various values of the parameters Sc, K ₁, K ₂ and K ₃ reveals the influence of the parameters on the flow, heat and mass transfer behavior.     

Mechanical Quadrature Methods and Extrapolation Algorithms for Boundary Integral Equations with Linear Boundary Conditions in Elasticity

By potential theory, elastic problems with linear boundary conditions are converted into boundary integral equations (BIEs) with logarithmic and Cauchy singularity. In this paper, a mechanical quadrature method (MQMs) is presented to deal with the logarithmic and the Cauchy singularity simultaneously for solving the boundary integral equations. The convergence and stability are proved based on Anselone??s collective compact and asymptotical compact theory. Furthermore, an asymptotic expansion with odd powers of errors is presented, which possesses high accuracy order O(h ³). Using h ³?Richardson extrapolation algorithms (EAs), the accuracy order of the approximation can be greatly improved to O(h ⁵), and an a posteriori error estimate can be obtained for constructing a self-adaptive algorithm. The efficiency of the algorithm is illustrated by examples.     

A critical state sand model with elastic–plastic coupling

Soil elastic moduli are highly pressure-dependent. Experimental findings have indicated that the elastic shear modulus of sands depends on p^χ, where p is mean principal effective stress and χ is a non-dimensional parameter. χ practically remains unchanged for shear strains less than 10⁻⁵ where the mechanical behavior is purely elastic. However, experiments have revealed that the emergence of plasticity for shear strains larger than 10⁻⁵ provokes a gradual increase in χ. Technically, this observation is an elastic–plastic coupling effect in which plasticity causes to change the elastic characteristics. Here, this issue is considered in hyper-elasticity framework in conjunction with a critical state compatible bounding surface plasticity platform for granular soils. To this aim, constitutive equations linking χ to a proper kinematic hardening parameter are presented. Then, using the proposed approach, a hyper-elastic theory is modified to consider the mentioned elastic–plastic coupling effect in the whole domain of the elastoplastic behavior. Adopting the improved hyper-elasticity necessitates the modification of a number of basic plasticity platform elements. In this regard, dilatancy and plastic hardening modulus of the bounding surface platform are modified. Successful performance of the modified constitutive model is presented against experimental data of loading/unloading triaxial tests.     

Multiple Timescales Analysis for 1:2 and 1:3 Resonant Hopf Bifurcations

The postcritical behavior of a generaln-dimensional system around a resonant double Hopf bifurcation isanalyzed. Both cases in which the critical eigenvalues are in ratios of1:2 and 1:3 are investigated. The Multiple Scale Method is employedto derive the bifurcation equations systematically in terms of thederivatives of the original vector field evaluated at the criticalstate. Expansions of the n-dimensional vector of state variables andof a three-dimensional vector of control parameters are performed interms of a unique perturbation parameter ε, of the order ofthe amplitude of motion. However, while resonant terms only appear atthe ε³-order in the 1:3 case, they already arise at theε²-order in the 1:2 case. Thus, by truncating theanalysis at the ε³-order in both cases, first orsecond-order bifurcation equations are respectively drawn, the latterrequiring resort to the reconstitution principle. A two-degrees-of-freedom system undergoing resonant double Hopf bifurcations isstudied. The complete postcritical scenario is analyzed in terms of thethree control parameters and the asymptotic results are compared withexact numerical integrations for both resonances. Branches of periodicas well as periodically modulated solutions are found and theirstability analyzed.     

Multiple spatial and temporal scales method for numerical simulation of non-classical heat conduction problems: one dimensional case

A multiple spatial and temporal scales method is studied to simulate numerically the phenomenon of non-Fourier heat conduction in periodic heterogeneous materials. The model developed is based on the higher-order homogenization theory with multiple spatial and temporal scales in one dimensional case. The amplified spatial scale and the reduced temporal scale are introduced respectively to account for the fluctuations of non-Fourier heat conduction due to material heterogeneity and non-local effect of the homogenized solution. By separating the governing equations into various scales, the different orders of homogenized non-Fourier heat conduction equations are obtained. The reduced time dependence is thus eliminated and the fourth-order governing differential equations are derived. To avoid the necessity of C¹ continuous finite element implementation, a C⁰ continuous mixed finite element approximation scheme is put forward. Numerical results are shown to demonstrate the efficiency and validity of the proposed method.     

Singular vortices in regular flows

A review of the theory of quasigeostrophic singular vortices embedded in regular flows is presented with emphasis on recent results. The equations governing the joint evolution of singular vortices and regular flow, and the conservation laws (integrals) yielded by these equations are presented. Using these integrals, we prove the nonlinear stability of a vortex pair on the f-plane with respect to any small regular perturbation with finite energy and enstrophy. On the β-plane, a new exact steady-state solution is presented, a hybrid regular-singular modon comprised of a singular vortex and a localized regular component. The unsteady drift of an individual singular β-plane vortex confined to one layer of a two-layer fluid is considered. Analysis of the β-gyres shows that the vortex trajectory is similar to that of a barotropic monopole on the β-plane. Non-stationary behavior of a dipole interacting with a radial flow produced by a point source in a 2D fluid is examined. The dipole always survives after collision with the source and accelerates (decelerates) in a convergent (divergent) radial flow.     

A Geometrical Approach to the Study¶of Unbounded Solutions¶of Quasilinear Parabolic Equations

In this article, we are interested in the existence and uniqueness of solutions for quasilinear parabolic equations set in the whole space ? ^N . We consider, in particular, cases when there is no restriction on the growth or the behavior of these solutions at infinity. Our model equation is the mean-curvature equation for graphs for which Ecker and Huisken have shown the existence of smooth solutions for any locally Lipschitz continuous initial data. We use a geometrical approach which consists in seeing the evolution of the graph of a solution as a geometric motion which is then studied by the so-called “level-set approach”. After determining the right class of quasilinear parabolic PDEs which can be taken into account by this approach, we show how the uniqueness for the original PDE is related to “fattening phenomena” in the level-set approach. Existence of solutions is proved using a local L ^∞ bound obtained by using in an essential way the level-set approach. Finally we apply these results to convex initial data and prove existence and comparison results in full generality, i.e., without restriction on their growth at infinity.     

A method of solving the diffusion problem for a vortex layer in a viscoplastic half-plane

A method is proposed to reduce the classical formulation of the problem to a system of two functional equations whose solution can be found numerically. A number of assertions that characterize the behavior of a rigid zone are proved. In particular, the lower estimate h ⁰(t) = 2b√t for the boundary motion is obtained; an explicit expression for b is given as a boundary stress function.     

An encapsulated dumbbell model for concentrated polymer solutions and melts II. Calculation of material functions and experimental comparisons

The consitutive equation derived in Part I [1] is used to compute steady-state and time-dependent material functions. A procedure is suggested for determining the four independent model parameters (λ_H, b, σ, β) from experimental data.Material functions have been calculated using parameters determined by experimental data for concentrated solutions of nearly monodisperse linear polymers. The theoretical curves for steady shear properties (η and Ψ₁) are in quantitative agreement with the data well into the range of nonlinear behavior. The calculated time-dependent properties (η*, η⁺, Ψ⁺, η^?) are in qualitative agreement with observed behavior. Theoretical curves for elongational viscosity are reasonable, but suitable experimental data are not available for comparison.