A new framework for solution of multidimensional population balance equations

A new framework is proposed in this work to solve multidimensional population balance equations (PBEs) using the method of discretization. A continuous PBE is considered as a statement of evolution of one evolving property of particles and conservation of their n internal attributes. Discretization must therefore preserve n + I properties of particles. Continuously distributed population is represented on discrete fixed pivots as in the fixed pivot technique of Kumar and Ramkrishna [1996a. On the solution of population balance equation by discretization-I A fixed pivot technique. Chemical Engineering Science 51(8), 1311-1332] for 1-d PBEs, but instead of the earlier extensions of this technique proposed in the literature which preserve 2(n) properties of non-pivot particles, the new framework requires n + I properties to be preserved. This opens up the use of triangular and tetrahedral elements to solve 2-d and 3-d PBEs, instead of the rectangles and cuboids that are suggested in the literature. Capabilities of computational fluid dynamics and other packages available for generating complex meshes can also be harnessed. The numerical results obtained indeed show the effectiveness of the new framework. It also brings out the hitherto unknown role of directionality of the grid in controlling the accuracy of the numerical solution of multidimensional PBEs. The numerical results obtained show that the quality of the numerical solution can be improved significantly just by altering the directionality of the grid, which does not require any increase in the number of points, or any refinement of the grid, or even redistribution of pivots in space. Directionality of a grid can be altered simply by regrouping of pivots.

Energy identities in water wave theory for free-surface boundary condition with higher-order derivatives

A modified form of Green's integral theorem is employed to derive the energy identity in any water wave diffraction problem in a single-layer fluid for free-surface boundary condition with higher-order derivatives. For a two-layer fluid with free-surface boundary condition involving higher-order derivatives, two forms of energy identities involving transmission and reflection coefficients for any wave diffraction problem are also derived here by the same method. Based on this modified Green's theorem, hydrodynamic relations such as the energy-conservation principle and modified Haskind–Hanaoka relation are derived for radiation and diffraction problems in a single as well as two-layer fluid.

Lightning-Induced Current and Voltage on a Rocket in the Presence of Its Trailing Exhaust Plume

This paper presents time-domain characteristics of induced current and voltage on a rocket in the presence of its exhaust plume when an electromagnetic (EM) wave generated by a nearby lightning discharge is incident on it. For the EM-field interaction with the rocket, the finite-difference time-domain technique has been used. The distributed electrical parameters, such as capacitance and inductance of the rocket and its exhaust plume, are computed using the method of moments technique. For the electrical characterization of the exhaust plume, the computational fluid dynamics technique has been used. The computed peak value of the electrical conductivity of the exhaust plume is 0.12 S/m near the exit plane and it reduces to 0.02 S/m at the downstream end. The relative permittivity varies from 0.91 to 0.99. The exhaust plume behaves as a good conductor for EM fields with frequencies less than 2.285 GHz. It has been observed that the peak value of the induced current on the rocket gets enhanced significantly in the presence of the conducting exhaust plume for the rocket and exhaust plume dimensions and parameters studied. The magnitude of the time-varying induced current at the tail is much more than that of any other section of the rocket.

Lightning-Induced Current and Voltage on a Rocket in the Presence of Its Trailing Exhaust Plume

This paper presents time-domain characteristics of induced current and voltage on a rocket in the presence of its exhaust plume when an electromagnetic (EM) wave generated by a nearby lightning discharge is incident on it. For the EM-field interaction with the rocket, the finite-difference time-domain technique has been used. The distributed electrical parameters, such as capacitance and inductance of the rocket and its exhaust plume, are computed using the method of moments technique. For the electrical characterization of the exhaust plume, the computational fluid dynamics technique has been used. The computed peak value of the electrical conductivity of the exhaust plume is 0.12 S/m near the exit plane and it reduces to 0.02 S/m at the downstream end. The relative permittivity varies from 0.91 to 0.99. The exhaust plume behaves as a good conductor for EM fields with frequencies less than 2.285 GHz. It has been observed that the peak value of the induced current on the rocket gets enhanced significantly in the presence of the conducting exhaust plume for the rocket and exhaust plume dimensions and parameters studied. The magnitude of the time-varying induced current at the tail is much more than that of any other section of the rocket.

A one point shock capturing kinetic scheme for hyperbolic conservation laws

We have presented a new low dissipative kinetic scheme based on a modified Courant Splitting of the molecular velocity through a parameter φ. Conditions for the split fluxes derived based on equilibrium determine φ for a one point shock. It turns out that φ is a function of the Left and Right states to the shock and that these states should satisfy the Rankine-Hugoniot Jump condition. Hence φ is utilized in regions where the gradients are sufficiently high, and is switched to unity in smooth regions. Numerical results confirm a discrete shock structure with a single interior point when the shock is aligned with the grid.

Time-dependent analysis of an N2O gasdynamic laser

A simplified two-temperature model is presented for the vibrational energy levels of the N2O and N2 molecules of an N2O-N2-He gasdynamic laser (GDL), and the governing equations for the unsteady flow of the gas mixture in a convergent-divergent contour nozzle are solved using a time-dependent numerical technique. Final steady-state distributions are obtained for vibrational temperatures, population inversion, and the small-signal laser gain along the nozzle. It is demonstrated that, for plenum temperatures lower than 1200 K, an N2O GDL such as the present is more efficient than a CO2 GDL in identical operating conditions

A numerical study of the role of the vertical structure of vorticity during tropical cyclone genesis

An eight-level axisymmetric model with simple parameterizations for clouds and the atmospheric boundary layer was developed to examine the evolution of vortices that are precursors to tropical cyclones. The effect of vertical distributions of vorticity, especially that arising from a merger of mid-level vortices, was studied by us to provide support for a new vortex-merger theory of tropical cyclone genesis. The basic model was validated with the analytical results available for the spin-down of axisymmetric vortices. With the inclusion of the cloud and boundary layer parameterizations, the evolution of deep vortices into hurricanes and the subsequent decay are simulated quite well. The effects of several parameters such as the initial vortex strength, radius of maximum winds, sea-surface temperature and latitude (Coriolis parameter) on the evolution were examined. A new finding is the manner in which mid-level vortices of the same strength decay and how, on simulated merger of these mid-level vortices, the resulting vortex amplifies to hurricane strength in a realistic time frame. The importance of sea-surface temperature on the evolution of full vortices was studied and explained. Also it was found that the strength of the surface vortex determines the time taken by the deep vortex to amplify to hurricane strength.

Induction reheating of A356.2 aluminum alloy and thixocasting as automobile component

The development work for producing an automobile component by thixocasting using A356.2 alloy was introduced. As the first step, the alloy was electromagnetically stirred and solidified to produce a billet with non-dendritic microstructure. The microstructure depended on several process parameters such as stirring intensity, stirring frequency, cooling rate, and melt initial superheat. Through a series of computational studies and controlled experiments, a set of process parameters were identified to produce the best microstructures. Reheating of a billet with non-dendritic microstructure to a semisolid temperature was the next step for thixo-casting of the components. The reheating process was characterized for various reheating cycles using a vertical-type reheating machine. The induction heating cycle was optimized to obtain a near-uniform temperature distribution in radial as well as axial direction of the billet, and the heating was continued until the liquid fraction reached about 50%. These parameters were determined with the help of a computational fluid dynamics (CFD) model of die filling and solidification of the semisolid alloy. The heated billets were subsequently thixo-cast into automobile components using a real-time controlled die casting machine. The results show that the castings are near net shape, free from porosity, good surface finish and have superior mechanical properties compared to those produced by conventional die casting processes using the same alloy.

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.

An efficient marching algorithm for waterhammer analysis by the method of characteristics

A new fast and efficient marching algorithm is introduced to solve the basic quasilinear, hyperbolic partial differential equations describing unsteady, flow in conduits by the method of characteristics. The details of the marching method are presented with an illustration of the waterhammer problem in a simple piping system both for friction and frictionless cases. It is shown that for the same accuracy the new marching method requires fewer computational steps, less computer memory and time.

Exact traveling- wave solution for model geophysical systems

Exact traveling-wave solutions of time-dependent nonlinear inhomogeneous PDEs, describing several model systems in geophysical fluid dynamics, are found. The reduced nonlinear ODEs are treated as systems of linear algebraic equations in the derivatives. A variety of solutions are found, depending on the rank of the algebraic systems. The geophysical systems include acoustic gravity waves, inertial waves, and Rossby waves. The solutions describe waves which are, in general, either periodic or monoclinic. The present approach is compared with the earlier one due to Grundland (1974) for finding exact solutions of inhomogeneous systems of nonlinear PDEs.

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.

A bicharacteristic formulation of the ideal MHD equations

On a characteristic surface Omega of a hyperbolic system of first-order equations in multi-dimensions (x, t), there exits a compatibility condition which is in the form of a transport equation along a bicharacteristic on Omega. This result can be interpreted also as a transport equation along rays of the wavefront Omega(t) in x-space associated with Omega. For a system of quasi-linear equations, the ray equations (which has two distinct parts) and the transport equation form a coupled system of underdetermined equations. As an example of this bicharacteristic formulation, we consider two-dimensional unsteady flow of an ideal magnetohydrodynamics gas with a plane aligned magnetic field. For any mode of propagation in this two-dimensional flow, there are three ray equations: two for the spatial coordinates x and y and one for the ray diffraction. In spite of little longer calculations, the final four equations (three ray equations and one transport equation) for the fast magneto-acoustic wave are simple and elegant and cannot be derived in these simple forms by use of a computer program like REDUCE.

Effect of wall emissivities on radiation heat transfer in glass tank forehearths

In an earlier work, we had proposed a two-band, non-grey radiative transfer model for heat transfer in forehearths with simultaneous optically thick and thin approximations for molten glass interiors and at boundaries. Here using the same model, the radiative interaction of the top-crown and bottom-refractory walls with interior layers of shallow molten glass is studied by varying the wall emissivities. The forehearth exit temperature profiles for higher wall emissivities (0.9) show better conditioning of the glass for white flint glasses (optically thin).

Interannual variations of Indian summer monsoon in a GCM: External conditions versus internal feedbacks

The potential predictability of the Indian summer monsoon due to slowly varying sea surface temperature (SST) forcing is examined. Factors responsible for limiting the predictability are also investigated. Three multiyear simulations with the R30 version of the Geophysical Fluid Dynamics Laboratory's climate model are carried out for this purpose, The mean monsoon simulated by this model is realistic including the mean summer precipitation over the Indian continent. The interannual variability of the large-scale component of the monsoon such as the "monsoon shear index" and its teleconnection with Pacific SST is well simulated by the model in a 15-yr integration with observed SST as boundary condition. On regional scales, the skill in simulating the interannual variability of precipitation over the Indian continent by the model is rather modest and its simultaneous correlation with eastern Pacific SST is negative but poor as observed. The poor predictability of precipitation over the Indian region in the model is related to the fact that contribution to the interannual variability over this region due to slow SST variations [El Nino-Southern Oscillation (ENSO) related] is comparable to those due to regional-scale fluctuations unrelated to ENSO SST. The physical mechanism through which ENSO SST tend to produce reduction in precipitation over the Indian continent is also elucidated. A measure of internal variability of the model summer monsoon is obtained from a 20-yr integration of the same model with fixed annual cycle SST as boundary conditions but with predicted soil moisture and snow cover. A comparison of summer monsoon indexes between this run and the observed SST run shows that the internal oscillations can account for a large fraction of the simulated monsoon variability. The regional-scale oscillations in the observed SST run seems to arise from these internal oscillations. It is discovered that most of the interannual internal variability is due to an internal quasi-biennial oscillation (QBO) of the model atmosphere. Such a QBO is also found in the author's third 18-yr simulation in which fixed annual cycle of SST as well as soil moisture and snow cover are prescribed. This shows that the model QBO is not due to land-surface-atmosphere interaction. It is proposed that the model QBO arises due to an interaction between nonlinear intraseasonal oscillations and the annual cycle. Spatial structure of the QBO and its role in limiting the predictability of the Indian summer monsoon is discussed.