Outward spherical solidification of a superheated melt with time dependent boundary flux

Short-time analytical solutions of solid and liquid temperatures and freezing front have been obtained for the outward radially symmetric spherical solidification of a superheated melt. Although results are presented here only for time dependent boundary flux, the method of solution can be used for other kinds of boundary conditions also. Later, the analytical solution has been compared with the numerical solution obtained with the help of a finite difference numerical scheme in which the grid points change with the freezing front position. An efficient method of execution of the numerical scheme has been discussed in details. Graphs have been drawn for the total solidification times and temperature distributions in the solid.

Flow and heat transfer in steady laminar compressible boundary layer swirling flow in a nozzle

The axisymmetric steady laminar compressible boundary layer swirling flow of a gas with variable properties in a nozzle has been investigated. The partial differential equations governing the non-similar flow have been transformed into new co-ordinates having finite ranges by means of a transformation which maps an infinite range into a finite one. The resulting equations have been solved numerically using an implicit finite-difference scheme. The computations have been carried out for compressible swirling flow through a convergent conical nozzle. The results indicate that the swirl exerts a strong influence on the longitudinal skin friction, but its effect on the tangential skin friction and heat transfer is comparatively small. The effect of the variation of the density-viscosity product across the boundary layer is appreciable only at low-wall temperature. The results are in good agreement with those of the local-similarity method for small values of the longitudinal distance.

Laminar compressible boundary layer flow at a three-dimensional stagnation point with vectored mass transfer

The effect of surface mass transfer velocities having normal, principal and transverse direction components (�vectored� suction and injection) on the steady, laminar, compressible boundary layer at a three-dimensional stagnation point has been investigated both for nodal and saddle points of attachment. The similarity solutions of the boundary layer equations were obtained numerically by the method of parametric differentiation. The principal and transverse direction surface mass transfer velocities significantly affect the skin friction (both in the principal and transverse directions) and the heat transfer. Also the inadequacy of assuming a linear viscosity-temperature relation at low-wall temperatures is shown.

Solution of a Boundary-Value Problem associated with Diffraction of Water Waves by a Nearly Vertical Barrier

An exact solution is derived for a boundary-value problem for Laplace's equation which is a generalization of the one occurring in the course of solution of the problem of diffraction of surface water waves by a nearly vertical submerged barrier. The method of solution involves the use of complex function theory, the Schwarz reflection principle, and reduction to a system of two uncoupled Riemann-Hilbert problems. Known results, representing the reflection and transmission coefficients of the water wave problem involving a nearly vertical barrier, are derived in terms of the shape function.

Radiative interactions in boundary layers

A detailed description of radiative interactions in laminar compressible boundary layers for moderate Mach numbers is presented by way of asymptotic analysis and supporting solutions. The radiation field is described by the differential approximation. While the asymptotic analysis is valid for large N (the ratio of photon mean free path to molecular mean free path) and arbitrary Boltzmann number, Bo (the ratio of convective heat flux to radiation heat flux), the solutions are obtained for Bo [double less-than sign] 1, the case of strong radiative interactions. The asymptotic analysis shows the existence of an optically thin boundary layer for large N and all Bo. For Bo [double less-than sign] 1, two outer regions are observed — one optically thin (at short distances from the leading edge) and the other optically thick (at large distances from the leading edge). An interesting feature not pointed out in the previous literature is the existence of a wall layer at large distances from the leading edge where convective heat flux can be ignored to the leading order of approximation. The radiation field in all cases can be very well approximated by a one-dimensional description. The solutions have been constructed using the ideas of matched asymptotic expansions by approximate analytical procedures and numerical methods. It is shown that, to the leading order of approximation, the radiation slip method yields exactly the same result as the more complicated matching procedure. Both the cases of linear and nonlinear radiation have been considered, the former being of interest in developing approximate methods which are subsequently generalized to handle the nonlinear problem. Detailed results are presented for both cases.

A Direct Approach to the Problem of Diffraction by a Half-Plane under Mixed Boundary Conditions

A general direct technique of solving a mixed boundary value problem in the theory of diffraction by a semi-infinite plane is presented. Taking account of the correct edge-conditions, the unique solution of the problem is derived, by means of Jones' method in the theory of Wiener-Hopf technique, in the case of incident plane wave. The solution of the half-plane problem is found out in exact form. (The far-field is derived by the method of steepest descent.) It is observed that it is not the Wiener-Hopf technique which really needs any modification but a new technique is certainly required to handle the peculiar type of coupled integral equations which the Wiener-Hopf technique leads to. Eine allgemeine direkte Technik zur Lösung eines gemischten Randwertproblems in der Theorie der Beugung an einer halbunendlichen Ebene wird vorgestellt. Unter Berücksichtigung der korrekten Eckbedingungen wird mit der Methode von Jones aus der Theorie der Wiener-Hopf-Technik die eindeutige Lösung für den Fall der einfallenden ebenen Welle hergeleitet. Die Lösung des Halbebenenproblems wird in exakter Form angegeben. (Das Fernfeld wurde mit der Methode des steilsten Abstiegs bestimmt.) Es wurde bemerkt, daß es nicht die Wiener-Hopf-Technik ist, die wirklich irgend welcher Modifikationen bedurfte. Gewiß aber wird eine neue Technik zur Behandlung des besonderen Typs gekoppelter Integralgleichungen benötigt, auf die die Wiener-Hopf-Technik führt.

Diffraction of a cylindrical pulse by a half plane under mixed boundary conditions

The general time dependent source problem has been solved by the method of transforms (Laplace, Lebedev–Kontorovich in succession) and the solution is obtained in the form of an infinite series involving Legendre functions. The solutions in the case of harmonic time dependence and the incident plane wave have been derived from the above solution and are presented in the form of an infinite series. In the case of an incident plane wave, the series has been summed and the final solution involves an improper integral which behaves like a complementary error function for large values of the argument. Finally, the far field evaluation has been shown. The results are compared with those of Sommerfeld's half-plane diffraction problem with unmixed boundary conditions.

A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

Boundary layer swirling flow through a conical hydrocyclone

The paper deals with the flow and heat-transfer problem of a steady axisymmetric laminar incompressible boundary layer swirling flow of a fluid through a conical hydrocyclone. The implicit finitedifference scheme is used to solve the partial differential equations governing the flow. The effect of swirl is found to be more pronounced on the longitudinal skin friction than on the tangential skin friction and heat transfer. The skin friction and heat transfer increase with swirl or with longitudinal distance. Swirl also gives rise to velocity overshoot in the longitudinal velocity profiles and the magnitude of the velocity overshoot increases as the swirl parameter increases. The results are found to be in good agreement with those of the local nonsimilarity and momentum integral methods but they differ appreciably from those of the local similarity method except for the longitudinal skin friction which is fairly in good agreement with that of the local similarity method.Die Arbeit beschäftigt sich mit der Strömung und dem Wärmeübergang in einem konischen Zyklon unter der Voraussetzung stationärer, achsensymmetrischer, laminarer, inkompressibler Grenzschichtströmung. Ein implizites Differenzenverfahren wird benutzt, um die partiellen Differentialgleichungen zu lösen. Der Einfluß des Dralls ist besonders ausgeprägt auf die longitudinale Komponente der Oberflächenreibung, weniger dagegen bei der tangentialen Komponente und beim Wärmeübergang. Die Oberflächenreibung und der Wärmeübergang nehmen zu mit dem Drall, sowie mit dem longitudinalen Abstand. Der Drall erzeugt ein Überschießen der Geschwindigkeit in der longitudinalen Abstand. Der Drall erzeugt ein Überschießen der Geschwindigkeit in der Längsrichtung. Die Größe des Überschusses nimmt mit wachsendem Drallparameter zu. Die Resultate stimmen gut mit den Ergebnissen der Theorie der lokalen Nichtähnlichkeit und der Impulsintegralmethode überein. Dagegen weichen sie mit Ausnahme der longitudinalen Komponente der Oberflächenreibung beträchtlich von den Resultaten der Theorie der lokalen Ähnlichkeit ab.

Nonsimilar boundary layers for non-Darcy mixed convection flow about a horizontal surface in a saturated porous medium

The nonsimilar non-Darcy mixed convection flow about a heated horizontal surface in a saturated porous medium has been studied when the surface temperature is a power function of distance (Tw = T∞ ± Axλ). The analysis is performed for the cases of parallel and stagnation flows with favourable induced pressure gradient. The partial differential equations governing the flow have been solved numerically using the Keller box method. The heat transfer is enhanced due to the buoyancy parameter and wall temperature, but the non-Darcy parameter reduces it. For non-Darcy flow, the similarity solution exists only for the case of parallel flow.

A sufficient criterion for Kelvin–Helmholtz instability in the magnetopause boundary‐layer region

The Kelvin–Helmholtz instability has been investigated for the magnetopause boundary‐layer region by the linearized method. The plasma in magnetosheath and magnetopause is assumed to be semi‐infinitely extended homogeneous, nondissipative, and incompressible. It is observed that, if one relation of two plasma speeds on the two sides of the magnetopause, wave number, and boundary‐layer thickness exceeds a certain threshold, the instability sets in. This new analytically sufficient criterion for excitation of instability in the three‐layer plasma flow generalizes the corresponding Chandrasekhar’s instability criterion for two‐layer plasma flow. The known results have been recovered and modified, the new results have been discovered. It is proved that the velocity threshold for the onset of instability is low when the magnitude of the magnetosheath and boundary‐layer region magnetic field and the angle between them are small. Also the threshold depends on the direction of plasma flow. The following results are observed numerically. The growth of the instability is sensitive to the magnetic field direction in the magnetosheath. A slight variation in the magnetic field direction in the second region can substantially change the relative velocity threshold for instability. When the ratio of the density of the second and third layer (magnetosphere) increases or that of the first and third layer decreases, the threshold decreases. Apart from this a necessary criterion for instability is obtained for a particular case.

Compressible MHD boundary layer in the stagnation region of a sphere

The steady laminar compressible boundary layer flow of an electrically conducting fluid in the stagnation region of a sphere with an applied magnetic field has been studied. The effects of the induced magnetic field, mass transfer, and viscous dissipation have been taken into account. Both isothermal and adiabatic wall conditions have been considered. The governing equations have been solved numerically using a shooting method. The skin friction and heat transfer are found to be strongly affected by the magnetic field, mass transfer, wall temperature and Mach number. It is found that the magnetic field reduces the heat transfer. This is a significant result which can be used in controlling the heat transfer rate. The boundary layer solutions break down as the magnetic parameter tends to a certain critical value.

Biorthogonal series method for Oseen type flows

A biorthogonal series method is developed to solve Oseen type flow problems. The theory leads to a new set of eigenfunctions for a specific class of linear non-selfadjoint operators containing the biharmonic one. These eigenfunctions differ from those given earlier in the literature for the biharmonic operator. The method is applied to the problem of thermocapillary flow in a cylindrical liquid bridge of finite length with axial through flow. Flow and temperature distributions are obtained at leading order of an expansion for small surface tension Reynolds number and Prandtl number. Another related problem considered is that of cylindrical cavity flow. Solutions for both cases are presented in terms of biorthogonal series. The effect of axial through flow on velocity and temperature fields is discussed by numerical evaluation of the truncated analytical series. The presence of axial through flow not only convectively shifts the vortices induced by surface forces in the direction of the through flow, but also moves their centers toward the outer cylindrical boundary. This process can lead to significantly asymmetric flow structures.

Time-domain-finite-wave analysis of the engine exhaust system by means of the stationary-frame method of characteristics. Part I. Theory

Time-domain-finite-wave analysis of the engine exhaust system is usually done using the method of characteristics. This makes use of either the moving frame method, or the stationary frame method. The stationary frame method is more convenient than its counterpart inasmuch as it avoids the tedium of graphical computations. In this paper (part I), the stationary-frame computational scheme along with the boundary conditions has been implemented. The analysis of a uniform tube, cavity-pipe junction including the engine and the radiation ends, and also the simple area discontinuities has been presented. The analysis has been done accounting for wall friction and heat-transfer for a one-dimensional unsteady flow. In the process, a few inconsistencies in the formulations reported in the literature have been pointed out and corrected. In the accompanying paper (part II) results obtained from the simulation are shown to be in good agreement with the experimental observations.

Accelerated gradient based diffuse optical tomographic image reconstruction

Purpose: Fast reconstruction of interior optical parameter distribution using a new approach called Broyden-based model iterative image reconstruction (BMOBIIR) and adjoint Broyden-based MOBIIR (ABMOBIIR) of a tissue and a tissue mimicking phantom from boundary measurement data in diffuse optical tomography (DOT). Methods: DOT is a nonlinear and ill-posed inverse problem. Newton-based MOBIIR algorithm, which is generally used, requires repeated evaluation of the Jacobian which consumes bulk of the computation time for reconstruction. In this study, we propose a Broyden approach-based accelerated scheme for Jacobian computation and it is combined with conjugate gradient scheme (CGS) for fast reconstruction. The method makes explicit use of secant and adjoint information that can be obtained from forward solution of the diffusion equation. This approach reduces the computational time many fold by approximating the system Jacobian successively through low-rank updates. Results: Simulation studies have been carried out with single as well as multiple inhomogeneities. Algorithms are validated using an experimental study carried out on a pork tissue with fat acting as an inhomogeneity. The results obtained through the proposed BMOBIIR and ABMOBIIR approaches are compared with those of Newton-based MOBIIR algorithm. The mean squared error and execution time are used as metrics for comparing the results of reconstruction. Conclusions: We have shown through experimental and simulation studies that Broyden-based MOBIIR and adjoint Broyden-based methods are capable of reconstructing single as well as multiple inhomogeneities in tissue and a tissue-mimicking phantom. Broyden MOBIIR and adjoint Broyden MOBIIR methods are computationally simple and they result in much faster implementations because they avoid direct evaluation of Jacobian. The image reconstructions have been carried out with different initial values using Newton, Broyden, and adjoint Broyden approaches. These algorithms work well when the initial guess is close to the true solution. However, when initial guess is far away from true solution, Newton-based MOBIIR gives better reconstructed images. The proposed methods are found to be stable with noisy measurement data. (C) 2011 American Association of Physicists in Medicine. DOI: 10.1118/1.3531572]