Analysis of transients in a canal network

A generalised formulation of the mathematical model developed for the analysis of transients in a canal network, under subcritical flow, with any realistic combination of control structures and their multiple operations, has been presented. The model accounts for a large variety of control structures such as weirs, gates, notches etc. discharging under different conditions, namely submerged and unsubmerged. A numerical scheme to compute and approximate steady state flow condition as the initial condition has also been presented. The model can handle complex situations that may arise from multiple gate operations. This has been demonstrated with a problem wherein the boundary conditions change from a gate discharge equation to an energy equation and back to a gate discharge equation. In such a situation the wave strikes a fixed gate and leads to large and rapid fluctuations in both discharge and depth.

The 26 July 2005 heavy rainfall event over Mumbai: numerical modeling aspects

The performance of the Advanced Regional Prediction System (ARPS) in simulating an extreme rainfall event is evaluated, and subsequently the physical mechanisms leading to its initiation and sustenance are explored. As a case study, the heavy precipitation event that led to 65 cm of rainfall accumulation in a span of around 6 h (1430 LT-2030 LT) over Santacruz (Mumbai, India), on 26 July, 2005, is selected. Three sets of numerical experiments have been conducted. The first set of experiments (EXP1) consisted of a four-member ensemble, and was carried out in an idealized mode with a model grid spacing of 1 km. In spite of the idealized framework, signatures of heavy rainfall were seen in two of the ensemble members. The second set (EXP2) consisted of a five-member ensemble, with a four-level one-way nested integration and grid spacing of 54, 18, 6 and 1 km. The model was able to simulate a realistic spatial structure with the 54, 18, and 6 km grids; however, with the 1 km grid, the simulations were dominated by the prescribed boundary conditions. The third and final set of experiments (EXP3) consisted of a five-member ensemble, with a four-level one-way nesting and grid spacing of 54, 18, 6, and 2 km. The Scaled Lagged Average Forecasting (SLAF) methodology was employed to construct the ensemble members. The model simulations in this case were closer to observations, as compared to EXP2. Specifically, among all experiments, the timing of maximum rainfall, the abrupt increase in rainfall intensities, which was a major feature of this event, and the rainfall intensities simulated in EXP3 (at 6 km resolution) were closest to observations. Analysis of the physical mechanisms causing the initiation and sustenance of the event reveals some interesting aspects. Deep convection was found to be initiated by mid-tropospheric convergence that extended to lower levels during the later stage. In addition, there was a high negative vertical gradient of equivalent potential temperature suggesting strong atmospheric instability prior to and during the occurrence of the event. Finally, the presence of a conducive vertical wind shear in the lower and mid-troposphere is thought to be one of the major factors influencing the longevity of the event.

Diffraction of a cylindrical pulse by a metallic half plane

A direct transform technique is applied to the initial and boundary value problem involving diffraction of a cylindrical pulse by a half plane, on which impedance type of boundary conditions must be met by the total field. The solution to the time harmonic incident plane wave is deduced as a particular case of the general time-dependent problem considered here and we avoid the Wiener–Hopf technique which leads to very complicated factorization and which masks the role of the impedance factor Z′ (a small quantity) in the expression for the scattered field.

On the possibility of an alpha-sq omega-type dynamo in a thin layer inside the sun

If the solar dynamo operates in a thin layer of 10,000-km thickness at the interface between the convection zone and the radiative core, using the facts that the dynamo should have a period of 22 years and a half-wavelength of 40 deg in the theta-direction, it is possible to impose restrictions on the values which various dynamo parameters are allowed to have. It is pointed out that the dynamo should be of alpha-sq omega nature, and kinematical calculations are presented for free dynamo waves and for dynamos in thin rectangular slabs with appropriate boundary conditions. An alpha-sq omega dynamo is expected to produce a significant poloidal field which does not leak to the solar surface. It is found that the turbulent diffusity eta and alpha-coefficient are restricted to values within about a factor of 10, the median values being eta of about 10 to the 10th sq cm/sec and alpha of about 10 cm/sec. On the basis of mixing length theory, it is pointed out that such values imply a reasonable turbulent velocity of the order 30 m/s, but rather small turbulent length scales like 300 km.

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.

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.

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.

Risk Minimizing Option Pricing for a Class of Exotic Options in a Markov-Modulated Market

We address risk minimizing option pricing in a regime switching market where the floating interest rate depends on a finite state Markov process. The growth rate and the volatility of the stock also depend on the Markov process. Using the minimal martingale measure, we show that the locally risk minimizing prices for certain exotic options satisfy a system of Black-Scholes partial differential equations with appropriate boundary conditions. We find the corresponding hedging strategies and the residual risk. We develop suitable numerical methods to compute option prices.

Influence of crack tip constraint on void growth in ductile FCC single crystals

Effect of constraint (stress triaxiality) on void growth near a notch tip in a FCC single crystal is investigated. Finite element simulations within the modified boundary layer framework are conducted using crystal plasticity constitutive equations and neglecting elastic anisotropy. Displacement boundary conditions based on model, elastic, two term K-T field are applied on the outer boundary of a large circular domain. A pre-nucleated void is considered ahead of a stationary notch tip. The interaction between the notch tip and the void is studied under different constraints (T-stress levels) and crystal orientations. It is found that negative T-stress retards the mechanisms of ductile fracture. However, the extent of retardation depends on the crystal orientation. Further, it is found that there exists a particular orientation which delays the ductile fracture processes and hence can potentially improve ductility. This optimal orientation depends on the constraint level. (C) 2010 Published by Elsevier B.V.

Applicability of Butler's equation in interpreting the thermodynamic behavior of surfaces and adsorption in Fe-S-O melts

Expressions for various second-order derivatives of surface tension with respect to composition at infinite dilution in terms of the interaction parameters of the surface and those of the bulk phases of dilute ternary melts have been presented. A method of deducing the parameters, which consists of repeated differentiation of Butler's equations with subsequent application of the appropriate boundary conditions, has been developed. The present investigation calculates the surface tension and adsorption functions of the Fe-S-O melts at 1873 and 1923 K using the modified form of Butler's equations and the derived values for the surface interaction parameters of the system. The calculated values are found to be in good agreement with those of the experimental data of the system. The present analysis indicates that the energetics of the surface phase are considerably different from those of the bulk phase. The present research investigates a critical compositional range beyond which the surface tension increases with temperature. The observed increase in adsorption of sulfur with consequent desorption of oxygen as a function of temperature above the critical compositional range has been ascribed to the increase of activity ratios of oxygen to sulfur in the surface relative to those in the bulk phase of the system.

Integral treatment for the representation of thermodynamic properties in multicomponent systems using interaction parameters

The partial thermodynamic functions of the solvent component of a ternary system have been deduced in terms of the interaction parameters by integration of several series which emerge from the Maclaurin infinite series based on the integral property of the system and subjected to appropriate boundary conditions. The series integration shows that the resulting partial functions are suitable for interpreting the thermodynamic properties of the system and are independent of compositional paths. In the present analysis, the higher order terms of these series are found to make insignificant contributions.

Breakage of a drop of inviscid fluid due to a pressure fluctuation at its surface

In this work, an attempt is made to gain a better understanding of the breakage of low-viscosity drops in turbulent flows by determining the dynamics of deformation of an inviscid drop in response to a pressure variation acting on the drop surface. Known scaling relationships between wavenumbers and frequencies, and between pressure fluctuations and velocity fluctuations in the inertial subrange are used in characterizing the pressure fluctuation. The existence of a maximum stable drop diameter d(max) follows once scaling laws of turbulent flow are used to correlate the magnitude of the disruptive forces with the duration for which they act. Two undetermined dimensionless quantities, both of order unity, appear in the equations of continuity, motion, and the boundary conditions in terms of pressure fluctuations applied on the surface. One is a constant of proportionality relating root-mean-square values of pressure and velocity differences between two points separated by a distance l. The other is a Weber number based on turbulent stresses acting on the drop and the resisting stresses in the drop due to interfacial tension. The former is set equal to 1, and the latter is determined by studying the interaction of a drop of diameter equal to d(max) with a pressure fluctuation of length scale equal to the drop diameter. The model is then used to study the breakage of drops of diameter greater than d(max) and those with densities different from that of the suspending fluid. It is found that, at least during breakage of a drop of diameter greater than d(max) by interaction with a fluctuation of equal length scale, a satellite drop is always formed between two larger drops. When very large drops are broken by smaller-length-scale fluctuations, highly deformed shapes are produced suggesting the possibility of further fragmentation due to instabilities. The model predicts that as the dispersed-phase density increases, d(max) decreases.

Obliquely incident water waves against a vertical cliff

Utilizing the commutativity property of the Cartesian coordinate differential operators arising in the boundary conditions associated with the propagation of surface water waves against a vertical cliff, under the assumptions of linearized theory, the problem of obliquely incident surface waves is considered for solution. The case of normal incidence, handled by previous workers follow as a particular limiting case of the present problem, which exhibits a source/sink type behavior of the velocity potential at the shore-line. An independent method of attack is also presented to handle the case of normal incidence.

Derivation and consistency of the partial functions of the ternary system involving interaction coefficients

The logarithm of activity coefficients of the components of the ternary system is derived based on the Maclaurin infinite series, which is expressed in terms of the integral property of the system and subjected to appropriate boundary conditions. The derivation of the functions involves extensive summation of various infinite series pertaining to the first-order interaction coefficients that have been shown completely to remove any truncational error. Since the conventional equations involving interaction coefficients are internally inconsistent, a consistent form of the partial functions is developed in the article using the technique just described. The thermodynamic consistency of the functions based on the Maxwell and the Gibbs-Duhem relations has been established. The derived values of the logarithmic activity coefficients of the components have been found to be in agreement with the thermodynamic data of the Fe-Cr-Ni system at 1873 K and have been found to be independent of the compositional paths.

Stability of laminated composite plates with cut-outs

Critical buckling loads of laminated fibre-reinforced plastic square panels have been obtained using the finite element method. Various boundary conditions, lay-up details, fibre orientations, cut-out sizes are considered. A 36 degrees of freedom triangular element, based on the classical lamination theory (CLT) has been used for the analysis. The performance of this element is validated by comparing results with some of those available in literature. New results have been given for several cases of boundary conditions for [0°/ ± 45°/90°]s laminates. The effect of fibre-orientation in the ply on the buckling loads has been investigated by considering [±?]6s laminates.