Reduction of nonholonomic variation problems to isoperimetric ones and connections in principal bundles

Translated from Matematicheskie Zametki, Vol. 49, No. 5, pp. 37–44, May, 1991.     

Calculation of the Radiating Characteristics of a Thermodynamically Nonuniform Gas Medium

Passive methods of remote analysis of waste gases of industrial enterprises and aircrafts engines and gases ejected in eruptions of volcanos, in fires, and in other processes disturbing the atmosphere are now of limited utility because of the lack of well-developed mathematical techniques for solving the inverse problems of gas analysis. However, the advantage of these methods is obvious, since the equipment is limited to a receiving device and the gas medium serves as a radiator. At the same time, in modern techniques, the interpretation of measurements of emissivities of gases is performed using a phenomenological approach such that experimental data on radiation characteristics of various objects are approximated by simple polynomial relations whose coefficients are determined by minimization methods. This approach does not allow one to investigate the mechanism of the processes occurring in a radiating medium and gives no way of solving the inverse problem of monitoring the composition and thermodynamic parameters of a medium which disturbs the atmosphere by thermal effects. In elaborating theoretical methods of gas analysis, the elaboration of the direct problem of calculation of radiating characteristics of gases is of importance, since it is not clear how the emissivity of a gas medium depends on its thermodynamic parameters. This paper is devoted to the above problem. Here, the error is estimated that arises in modeling a nonuniform gas medium with an actual temperature distribution, which linearly depends on coordinate, by a uniform gas layer. It is shown that the difference between the radiation intensities calculated for uniform and nonuniform media can be significant in the practically important case where the recording device is in a zone of normal temperatures. This implies that the error of reconstruction of the spectroscopic and molecular parameters of gaseous compounds from high-temperature measurements of the radiating and absorbing characteristics should take into account the thermodynamic nonuniformity of the medium.     

Research problems

This section appears from time to time. Contributions are invited, and should be submitted to R. C. Thompson, Mathematics Department, University of California, Santa Barbara, CA 93106.     

Numerical study of simultaneous natural convection heat transfer from both surfaces of a uniformly heated thin plate with arbitrary inclination

 Numerical studies were conducted to investigate the natural convection heat transfer around a uniformly heated thin plate with arbitrary inclination in an infinite space. The numerical approach was based on the finite volume technique with a nonstaggered grid arrangement. For handling the pressure–velocity coupling the SIMPLE-algorithm was used. QUICK scheme and first order upwind scheme were employed for discretization of the momentum and energy convective terms respectively. Plate width and heating rate were used to vary the modified Rayleigh number over the range of 4.8×10⁶ to 1.87×10⁸. Local and average heat transfer characteristics were compared with regarding to the inclination angle. The empirical expressions for local and average Nusselt number were correlated. It has been found that for inclination angle less than 10^∘, the flow and heat transfer characteristics are complicated and the average Nusselt number can not be correlated by one equation while for inclination angle larger than 10^∘, the average Nusselt number can be correlated into an elegant correlation. Received on 18 April 2001 / Published online: 29 November 2001     

Finite element method in the solution of the Euler and Navier-Stokes equations for internal flow

Finite element solutions of the Euler and Navier-Stokes equations are presented, using a simple dissipation model. The discretization is based on the weak-Galerkin weighted residual method and equal interpolation functions for all the unknowns are permitted. The nonlinearity is iterated upon using a Newton method and at each iteration the linear algebraic system is solved by a direct solver with all unknowns fully coupled. Results are presented for two-dimensional transonic inviscid flows and two- and three-dimensional incompressible viscous flows. Convergence of the algorithm is shown to be quadratic, reaching machine accuracy in very few iterations. The inviscid results demonstrate the existence of nonunique numerical solutions to the steady Euler equations.     

A note on dephasing and pair-breaking

It is shown that the dephasing which suppresses the weak localization correction to the conductivity has the same physical origin than pair-breaking in superconductors and thus-following de Gennes-may be expressed in terms of the correlation function of the operator for time reversal.     

Theoretical and experimental investigation of the linear optical response of YBa2Cu4O8 The relevance of LDA calculations

We report results of ellipsometric measurements of the dielectric tensor of YBa₂Cu₄O₈ as well as calculations employing density functional theory in the local density approximation (LDA). For the computation of the one-electron eigenvalues and states the linear-muffin-tin-orbital method (LMTO) was used. Measured and calculated results show good agreement. In particular, the calculation predicts different band structures for YBa₂Cu₃O₇ and YBa₂Cu₄O₈ in the vicinity of the Fermi energy which lead to significant differences in the optical spectra. One such difference is a region of small ε₂ in the near infrared which is a fingerprint of the high (relative to the top of the highest fully occupied band) Fermi energy as compared to the situation in YBa₂Cu₃O₇. Our experimental results confirm the theoretical findings.     

A globally convergent multi-grid algorithm for moving boundary problems of two-phase Stefan type

In this paper we deal with the numerical solution of movingboundary problems of two-phase Stefan type. Based on an implicitdiscretization in time and the use of continuous, piecewiselinear finite elements in the space variables with respect tothe weak formulation of the problem, a globally convergent multi-gridalgorithm is developed. That algorithm strongly relies on thevariational characterization of the fully discretized problemas the unconstrained minimization of a subdifferentiable convexobjective functional. Numerical results indicate a significantimprovement in efficiency compared with previous multi-gridapproaches.