Inverse problem for the wave equation : partial data and novel boundary sources

An inverse problem for the wave equation is a mathematical formulation of the problem to convert measurements of sound waves to information about the wave speed governing the propagation of the waves. This doctoral thesis extends the theory on the inverse problems for the wave equation in cases with partial measurement data and also considers detection of discontinuous interfaces in the wave speed. A possible application of the theory is obstetric sonography in which ultrasound measurements are transformed into an image of the fetus in its mother's uterus. The wave speed inside the body can not be directly observed but sound waves can be produced outside the body and their echoes from the body can be recorded. The present work contains five research articles. In the first and the fifth articles we show that it is possible to determine the wave speed uniquely by using far apart sound sources and receivers. This extends a previously known result which requires the sound waves to be produced and recorded in the same place. Our result is motivated by a possible application to reflection seismology which seeks to create an image of the Earth s crust from recording of echoes stimulated for example by explosions. For this purpose, the receivers can not typically lie near the powerful sound sources. In the second article we present a sound source that allows us to recover many essential features of the wave speed from the echo produced by the source. Moreover, these features are known to determine the wave speed under certain geometric assumptions. Previously known results permitted the same features to be recovered only by sequential measurement of echoes produced by multiple different sources. The reduced number of measurements could increase the number possible applications of acoustic probing. In the third and fourth articles we develop an acoustic probing method to locate discontinuous interfaces in the wave speed. These interfaces typically correspond to interfaces between different materials and their locations are of interest in many applications. There are many previous approaches to this problem but none of them exploits sound sources varying freely in time. Our use of more variable sources could allow more robust implementation of the probing.

Analytical and numerical solutions of inward spherical solidification of a superheated melt with radiative-convective heat transfer and density jump at freezing front

Analytical and numerical solutions of a general problem related to the radially symmetric inward spherical solidification of a superheated melt have been studied in this paper. In the radiation-convection type boundary conditions, the heat transfer coefficient has been taken as time dependent which could be infinite, at time,t=0. This is necessary, for the initiation of instantaneous solidification of superheated melt, over its surface. The analytical solution consists of employing suitable fictitious initial temperatures and fictitious extensions of the original region occupied by the melt. The numerical solution consists of finite difference scheme in which the grid points move with the freezing front. The numerical scheme can handle with ease the density changes in the solid and liquid states and the shrinkage or expansions of volumes due to density changes. In the numerical results, obtained for the moving boundary and temperatures, the effects of several parameters such as latent heat, Boltzmann constant, density ratios, heat transfer coefficients, etc. have been shown. The correctness of numerical results has also been checked by satisfying the integral heat balance at every timestep.

Analytical solution of the laminar mean flow wave equation in a lined or unlined two-dimensional rectangular duct

The prediction of the sound attenuation in lined ducts with sheared mean flow has been a topic of research for many years. This involves solving the sheared mean flow wave equation, satisfying the relevant boundary condition. As far as the authors' knowledge goes, this has always been done using numerical techniques. Here, an analytical solution is presented for the wave propagation in two-dimensional rectangular lined ducts with laminar mean flow. The effect of laminar mean flow is studied for both the downstream and the upstream wave propagation. The attenuation values predicted for the laminar mean flow case are compared with those for the case of uniform mean flow. Analytical expressions are derived for the transfer matrices.

Analysis of heat and mass transfer during bulk hydraircooling of spherical food products

Hydraircooling is a technique used for precooling food products. In this technique chilled water is sprayed over the food products while cold unsaturated air is blown over them. Hydraircooling combines the advantages of both air- and hydrocooling. The present study is concerned with the analysis of bulk hydraircooling as it occurs in a package filled with several layers of spherical food products with chilled water sprayed from the top and cold unsaturated air blown from the bottom. A mathematical model is developed to describe the hydrodynamics and simultaneous heat and mass transfer occurring inside the package. The non-dimensional governing equations are solved using the finite difference numerical methods. The results are presented in the form of time-temperature charts. A correlation is obtained to calculate the process time in terms of the process parameters.

Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

A thermodynamic protocol for explaining subcriticality and other subtle features in self-deflagrating solids

A novel universal approach to understand the self-deflagration in solids has been attempted by using basic thermodynamic equation of partial differentiation, where burning mte depends on the initial temperature and pressure of the system. Self-deflagrating solids are rare and are reported only in few compounds like ammonium perchlorate (AP), polystyrene peroxide and tetrazole. This approach has led us to understand the unique characteristics of AP, viz. the existence of low pressure deflagration limit (LPL 20 atm), hitherto not understood sufficiently. This analysis infers that the overall surface activation energy comprises of two components governed by the condensed phase and gas phase processes. The most attractive feature of the model is the identification of a new subcritical regime I' below LPL where AP does not burn. The model is aptly supported by the thermochemical computations and temperature-profile analyses of the combustion train. The thermodynamic model is further corroborated from the kinetic analysis of the high pressure (1-30 atm) DTA thermograms which affords distinct empirical decomposition rate laws in regimes I' and 1 (20-60 atm). Using Fourier-Kirchoff one dimensional heat transfer differential equation, the phase transition thickness and the melt-layer thickness have been computed which conform to the experimental data.

Casting and heat treatment variables of Al–7Si–Mg alloy

A review of the research work that has been carried out thus far relating the casting and heat treatment variables to the structure and mechanical properties of Al–7Si–Mg (wt-%) is presented here. Although specifications recommend a wide range of magnesium contents and a fairly high content of iron, a narrow range of magnesium contents, closer to either the upper or lower specified limits depending on the properties desired, and a low iron content will have to be maintained to obtain optimum and consistent mechanical properties. A few studies have revealed that the modification of eutectic silicon slightly increases ductility and fracture toughness and also that the effect of modification is predominant at low iron content. Generally, higher solidification rates give superior mechanical properties. Delayed aging (the time elapsed between quenching and artificial aging during precipitation hardening) severely affects the strength of the alloy. The mechanism of delayed aging can be explained on the basis of Pashley's kinetic model. It has been reported that certain trace additions (cadmium, indium, tin, etc.) neutralise the detrimental effect of delayed aging. In particular, it should be noted that delayed aging is not mentioned in any of the specifications. With reference to the mechanism by which trace additions neutralise the detrimental effect of delayed aging, various hypotheses have been postulated, of which impurity–vacancy interaction appears to be the most widely accepted.

Reynolds stress and heat flux in spherical shell convection

Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Aims. We compute turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime in order to study whether restricted geometry introduces artefacts into the results. Methods. We employ direct numerical simulations of turbulent convection in spherical and Cartesian geometries. In order to alleviate the computational cost in the spherical runs and to reach as high spatial resolution as possible, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number or inverse Rossby number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in accordance with earlier numerical results, but in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong 'banana cells'. Their effect in the spherical case does not appear to be as large. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. Longitudinal heat flux is always in the retrograde direction. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor-Proudman balance.

Exact analysis of Burgers equation on semiline with flux condition at the origin

The initial boundary value problem for the Burgers equation in the domain x greater-or-equal, slanted 0, t > 0 with flux boundary condition at x = 0 has been solved exactly. The behaviour of the solution as t tends to infinity is studied and the “asymptotic profile at infinity” is obtained. In addition, the uniqueness of the solution of the initial boundary value problem is proved and its inviscid limit as var epsilon → 0 is obtained.

Heat capacity measurements on Ge20Se80−xBix glasses

The specific heat Cp of glassy Ge20Se80−xBix (0 ≤ × ≤ 12) samples is investigated. The Cp at 323K and the ΔCp at glass transition temperature Tg1 show anomalous features around x = 8 at.%, where p−n conduction type inversion also take place. These features are discussed in the light of Phillips model of phase separation in these glasses at the microscopic level.

Nonaxisymmetric unsteady motion over a rotating disk in the presence of free convection and magnetic field

The nonaxisymmetric unsteady motion produced by a buoyancy-induced cross-flow of an electrically conducting fluid over an infinite rotating disk in a vertical plane and in the presence of an applied magnetic field normal to the disk has been studied. Both constant wall and constant heat flux conditions have been considered. It has been found that if the angular velocity of the disk and the applied magnetic field squared vary inversely as a linear function of time (i.e. as (1??t*)?1, the governing Navier-Stokes equation and the energy equation admit a locally self-similar solution. The resulting set of ordinary differential equations has been solved using a shooting method with a generalized Newton's correction procedure for guessed boundary conditions. It is observed that in a certain region near the disk the buoyancy induced cross-flow dominates the primary von Karman flow. The shear stresses induced by the cross-flow are found to be more than these of the primary flow and they increase with magnetic parameter or the parameter ? characterizing the unsteadiness. The velocity profiles in the x- and y-directions for the primary flow at any two values of the unsteady parameter ? cross each other towards the edge of the boundary layer. The heat transfer increases with the Prandtl number but reduces with the magnetic parameter.

Simultaneous heat and mass transfer under unsteady mixed convection along a vertical slender cylinder embedded in a porous medium

Unsteady laminar mixed convection flow (combined free and forced convection flow) along a vertical slender cylinder embedded in a porous medium under the combined buoyancy effect of thermal and species diffusion has been studied. The effect of the permeability of the medium as well as the magnetic field has been included in the analysis. The partial differential equations with three independent variables governing the flow have been solved numerically using a implicit finite difference scheme in combination with the quasilinearization technique. Computations have been carried out for accelerating, decelerating and oscillatory free stream velocity distributions. The effects of the permeability of the medium, buoyancy forces, transverse curvature and magnetic field on skin friction, heat transfer and mass transfer have been studied. It is found that the effect of free stream velocity distribution is more pronounced on the skin friction than on the heat and mass transfer. The permeability and magnetic parameters increase the skin friction, but reduce the heat and mass transfer. The skin friction, heat transfer and mass transfer are enhanced due to the buoyancy forces and curvature parameter. The heat transfer is strongly dependent on the viscous dissipation parameter and the Prandtl number, and the mass transfer on the Schmidt number. Untersucht wurde die instationäre laminare Mischkonvektion längs eines vertikalen, in einem porösen Medium eingebetteten Zylinders unter kombinierten Auftriebseffekten von thermischer und spezieller Diffusion. Der Einfluß der Permeabilität des Mediums sowie des magnetischen Feldes wurden in die Betrachtung einbezogen. Die partiellen Differentialgleichungen mit drei unabhängigen Variablen, welche die Strömung beschreiben, wurde numerisch anhand des Schemas der endlichen Differenzen in Verbindung mit der Technik der Quasilinearisation gelöst. Berechnungen für die beschleunigte, verzögerte und oszillierende Geschwindigkeitsverteilung der freien Strömung sind durchgeführt worden. Untersucht wurden ebenfalls die Effekte der Permeabilität des Mediums, der Auftriebskräfte, der transversalen Krümmung, des magnetischen Feldes auf die Oberflächenreibung sowie die Wärmeund Stoffübertragung. Es wurde festgestellt, daß die Geschwindigkeit mehr Einfluß auf die Oberflächenreibung als auf die Wärmeund Stoffübertragung hat. Die Oberflächenreibung sowie die Wärme- und Stoffübertragung werden durch die Auftriebskräfte und die Krümmungsparameter verbessert. Die Wärmeübertragung ist stark abhängig von den Parametern der viskosen Dissipation und der Prandtl-Zahl; die Stoffübertragung von der Schmidt-Zahl.

Solubility of sulfur dioxide at low partial pressures in dilute sulfuric acid solutions

Equilibrium of dissolution of sulfur dioxide at ppm levels in aqueous solutions of dilute sulfuric acid is analyzed, and a general expression is derived relating the total concentration of sulfur dioxide in the liquid phase to the partial pressure of SO2 in the gas and to the concentration of sulfuric acid in the solution. The equation is simplified for zero and high concentrations of the acid. Experiments at high concentrations of sulfuric acid have enabled the direct determination of Henry’s constant and its dependency on temperature. Heat of dissolution is -31.47 kJ/mol. Experiments in the absence of sulfuric acid and the related simplified expression have led to the determination of the equilibrium constant of the hydrolysis of aqueous sulfur dioxide and its temperature dependency.The heat of hydrolysis is 15.69 kJ/mol. The model equation with these parameters predicts the experimental data of the present work as well as the reported data very well.

Specific-heat measurements during cooling through the glass-transition region

In this paper we report the measurements of specific heats of five glass formers as they are cooled through the glass-transition region. The measurements are compared with other specific-heat measurements such as adiabatic-calorimetry and ac-calorimetry measurements. The data are then analyzed using a model of enthalpy relaxation and nonequilibrium cooling, which can track the nonequilibrium relaxation time tau(S). The relevant parameters that describe tau(S) are obtained, allowing us to compare the enthalpy-relaxation times obtained from this method with other methods. We display the clear connection of the unrelaxed enthalpy with the nonequilibrium relaxation time and also show the role played by the delayed heat release from the unrelaxed enthalpy in the glass-transition region. We have also made certain definite observations regarding the equilibrium configurational specific heat and the Vogel-Fulcher law, which describes tau(S).