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.

Viscous flow computations using a meshless solver, LSFD-U

The Upwind-Least Squares Finite Difference (LSFD-U) scheme has been successfully applied for inviscid flow computations. In the present work, we extend the procedure for computing viscous flows. Different ways of discretizing the viscous fluxes are analysed for the positivity, which determines the robustness of the solution procedure. The scheme which is found to be more positive is employed for viscous flux computation. The numerical results for validating the procedure are presented.

Parallelization of an unstructured data based cell centre finite volume code, HIFUN-3D

This work describes the parallelization of High Resolution flow solver on unstructured meshes, HIFUN-3D, an unstructured data based finite volume solver for 3-D Euler equations. For mesh partitioning, we use METIS, a software based on multilevel graph partitioning. The unstructured graph used for partitioning is associated with weights both on its vertices and edges. The data residing on every processor is split into four layers. Such a novel procedure of handling data helps in maintaining the effectiveness of the serial code. The communication of data across the processors is achieved by explicit message passing using the standard blocking mode feature of Message Passing Interface (MPI). The parallel code is tested on PACE++128 available in CFD Center

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.

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.

Unsteady mixed convection with double diffusion over a horizontal cylinder and a sphere within a porous medium

The combined effects of the permeability of the medium, magnetic field, buoyancy forces and dissipation on the unsteady mixed convection flow over a horizontal cylinder and a sphere embedded in a porous medium have been studied. The nonlinear coupled partial differential equations with three independent variables have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The skin friction, heat transfer and mass transfer increase with the permeability of the medium, magnetic field and buoyancy parameter. The heat and mass transfer continuously decrease with the stream-wise distance, whereas the skin friction increases from zero, attains a maximum and then decreases to zero. The skin friction, heat transfer and mass transfer are significantly affected by the free stream velocity distribution. The effect of dissipation parameter is found to be more pronounced on the heat transfer than on the skin friction and mass transfer

Second-order boundary-layer effects for the unsteady laminar incompressible three-dimensional stagnation-point flow

All the second-order boundary-layer effects on the unsteady laminar incompressible flow at the stagnation-point of a three-dimensional body for both nodal and saddle point regions have been studied. It has been assumed that the free-stream velocity, wall temperature and mass transfer vary arbitrarily with time. The effect of the Prandtl number has been taken into account. The partial differential equations governing the flow have been derived for the first time and then solved numerically unsteady free-stream velocity distributions, the nature of the using an implicit finite-difference scheme. It is found that the stagnation point and the mass transfer strongly affect the skin friction and heat transfer whereas the effects of the Prandtl number and the variation of the wall temperature with time are only on the heat transfer. The skin friction due to the combined effects of first- and second-order boundary layers is less than the skin friction due to, the first-order boundary layers whereas the heat transfer has the opposite behaviour. Suction increases the skin friction and heat transfer but injection does the opposite

Darcian mixed convection plumes along vertical adiabatic surfaces in a saturated porous medium

The mixed convection flow due to a line thermal source embedded at the leading edge of an adiabatic vertical plane surface immersed in a saturated porous medium has been studied. Both weakly and strongly buoyant plume regimes have been considered. The cases of buoyancy assisting and buoyancy opposing flow conditions have been incorporated in the analysis. The results are presented for the entire range of buoyancy parameter from the pure forced convection (xgr=0) to the pure free convection (xgr rarr infin@#@) regimes. For buoyancy-assisting flow, the wall temperature and the velocity at the wall increase as the plume strength increases. However, they all decrease as the free-stream velocity increases. For buoyancyopposing flow, the temperature at the wall increases as the strength of the plume increases but velocity at the wall decreases.

Non-darcy double-diffusive mixed convection from heated vertical and horizontal plates in saturated porous media

The non-darcy mixed convection flows from heated vertical and horizontal plates in saturated porous media have been considered using boundary layer approximations. The flows are considered to be driven by multiple buoyancy forces. The similarity solutions for both vertical and horizontal plates have been obtained. The governing equations have been solved numerically using a shooting method. The heat transfer, mass transfer and skin friction are reduced due to inertial forces. Also, they increase with the buoyancy parameter for aiding flow and decrease for the opposing flow. For aiding flow, the heat and mass transfer coefficients are found to approach asymptotically the forced or free convection values as the buoyancy parameter approaches zero or infinity.

Unsteady mixed convection over two-dimensional and axisymmetric bodies

The unsteady laminar incompressible mixed convection flow over a two-dimensional body (cylinder) and an axisymmetric body (sphere) has been studied when the buboyancy forces arise from both thermal and mass diffusion and the unsteadiness in the flow field is introduced by the time dependent free stream velocity. The nonlinear partial differential equations with three independent variables governing the flow have been solved numerically using an implicit finite-difference scheme in combination with the quasilinearization technique. The results indicate that for the thermally assisting flow the local skin friction, heat transfer and mass diffusion are enhanced when the buoyancy force from mass diffusion assists the thermal buoyancy force. But this trend is opposite for the thermally opposing flow. The point of zero skin friction moves upstream due to unsteadiness. No singularity is observed at the point of zero skin friction for unsteady flow unlike steady flow. The flow reversal is observed after a certain instant of time. The velocity overshoot occurs for assisting flows.

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.

Unsteady mixed convection flow in stagnation region adjacent to a vertical surface

The unsteady two-dimensional laminar mixed convection flow in the stagnation region of a vertical surface has been studied where the buoyancy forces are due to both the temperature and concentration gradients. The unsteadiness in the flow and temperature fields is caused by the time-dependent free stream velocity. Both arbitrary wall temperature and concentration, and arbitrary surface heat and mass flux variations have been considered. The Navier-Stokes equations, the energy equation and the concentration equation, which are coupled nonlinear partial differential equations with three independent variables, have been reduced to a set of nonlinear ordinary differential equations. The analysis has also been done using boundary layer approximations and the difference between the solutions has been discussed. The governing ordinary differential equations for buoyancy assisting and buoyancy opposing regions have been solved numerically using a shooting method. The skin friction, heat transfer and mass transfer coefficients increase with the buoyancy parameter. However, the skin friction coefficient increases with the parameter lambda, which represents the unsteadiness in the free stream velocity, but the heat and mass transfer coefficients decrease. In the case of buoyancy opposed flow, the solution does not exist beyond a certain critical value of the buoyancy parameter. Also, for a certain range of the buoyancy parameter dual solutions exist.

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.