A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

A fuzzy dynamic flood routing model for natural channels

A fuzzy dynamic flood routing model (FDFRM) for natural channels is presented, wherein the flood wave can be approximated to a monoclinal wave. This study is based on modification of an earlier published work by the same authors, where the nature of the wave was of gravity type. Momentum equation of the dynamic wave model is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. Hence, the FDFRM gets rid of the assumptions associated with the momentum equation. Also, it overcomes the necessity of calculating friction slope (S-f) in flood routing and hence the associated uncertainties are eliminated. The fuzzy rule based model is developed on an equation for wave velocity, which is obtained in terms of discontinuities in the gradient of flow parameters. The channel reach is divided into a number of approximately uniform sub-reaches. Training set required for development of the fuzzy rule based model for each sub-reach is obtained from discharge-area relationship at its mean section. For highly heterogeneous sub-reaches, optimized fuzzy rule based models are obtained by means of a neuro-fuzzy algorithm. For demonstration, the FDFRM is applied to flood routing problems in a fictitious channel with single uniform reach, in a fictitious channel with two uniform sub-reaches and also in a natural channel with a number of approximately uniform sub-reaches. It is observed that in cases of the fictitious channels, the FDFRM outputs match well with those of an implicit numerical model (INM), which solves the dynamic wave equations using an implicit numerical scheme. For the natural channel, the FDFRM Outputs are comparable to those of the HEC-RAS model.

A numerical study of division of flow in open channels

A two-dimensional numerical model which employs the depth-averaged forms of continuity and momentum equations along with k-e turbulence closure scheme is used to simulate the flow at the open channel divisions. The model is generalised to flows of arbitrary geometries and MacCormack finite volume method is used for solving governing equations. Application of cartesian version of the model to analyse the flow at right-angled junction is presented. The numerical predictions are compared with experimental data of earlier investigators and measurements made as part of the present study. Performance of the model in predicting discharge distribution, surface profiles, separation zone parameters and energy losses is evaluated and discussed in detail. To illustrate the application of the numerical model to analyse the flow in acute angled offtakes and streamlined branch entries, a few computational results are presented.

Numerical modeling of heat and mass transfer in laser surface alloying: Non-equilibrium solidification effects

A transient macroscopic model is developed for studying heat and mass transfer in a single-pass laser surface alloying process, with particular emphasis on non-equilibrium solidification considerations. The solution for species concentration distribution requires suitable treatment of non-equilibrium mass transfer conditions. In this context, microscopic features pertaining to non-equilibrium effects on account of solutal undercooling are incorporated through the formulation of a modified partition-coefficient. The effective partition-coefficient is numerically modeled by Means of a number of macroscopically observable parameters related to the solidifying domain. The numerical model is so developed that the modifications on account of non-equilibrium solidification considerations can be conveniently implemented in existing numerical codes based on equilibrium solidification considerations.

Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Nonlinear cyclic load analysis for lateral response of batter piles in soft clay with a rigorous degradation model

Nonlinear analysis of batter piles in soft clay is performed using the finite element technique. As the batter piles are not only governed by lateral load but also axial load, the effect of P- Delta moment and geometric stiffness matrix is included in the analysis. For implementing the nonlinear soil behavior, reduction in soil strength (degradation), and formation of gap with number of load cycles, a numerical model is developed where a hyperbolic relation is adopted for the soil in static condition and hyperbolic relation considering degradation and gap for cyclic load condition. The numerical model is validated with published experimental results for cyclic lateral loading and the hysteresis loops are developed to predict the load-deflection behavior and soil resistance behavior during consecutive cycles of loading. This paper highlights the importance of a rigorous degradation model for subsequent cycles of loading on the pile-soil system by a hysteretic representation.

The upper Bay of Bengal salinity structure in a high-resolution model

Salinity in the Bay of Bengal is highly heterogeneous, with extremely fresh waters found at the surface in the Northern part of the basin, and saltier waters at subsurface as well as to the south. This paper investigates the seasonal structure of sea surface salinity of the Bay in a regional high-resolution model forced by ERA-Interim reanalysis and various precipitation products. Surface circulation is believed to drive the spreading of northern Bay of Bengal fresh waters to the rest of the Indian Ocean. We first present a series of experiments to infer the sensitivity of modeled circulation to various numerical choices. Surface circulation is found to be sensitive to the horizontal resolution of the model, with the 1/12 degrees version appearing much more realistic than the 1/4 degrees version. The sidewall boundary condition is also drastically influencing the characteristics of the western boundary current simulated. We then investigate the sensitivity of the salinity response to the various precipitation products. We observe that ERA-Interim excess precipitation induces a fresh bias in the surface salinity response. Spaceborne precipitation products are more satisfactory. We then identify the pathways of the northern Bay freshwater mass, based on passive tracers experiments. Our model suggests that over timescales of a few months, vertical exchanges between the upper fresh layer and the underlying saltier layer appear to be the main export pathway for the freshwater. The horizontal circulation within the mixed layer also acts to convey fresh waters out of the Bay at these timescales, but in a lesser quantity compared to the vertical export. Beyond its intrinsic interest for the understanding of Bay of Bengal physics, this study highlights the need for a careful design of any realistic numerical model, in three key aspects: the choice of the resolution of the model, the choice of the sub-grid scale parameterizations, and the choice of the forcing fluxes. (C) 2013 Elsevier Ltd. All rights reserved.

Numerical study of solidification of A356 aluminum alloy flowing on an oblique plate with experimental validation

Molten A356 aluminum alloy flowing on an oblique plate is water cooled from underneath. The melt partially solidifies on plate wall with continuous formation of columnar dendrites. These dendrites are continuously sheared off into equiaxed/fragmented grains and carried away with the melt by producing semisolid slurry collected at plate exit. Melt pouring temperature provides required solidification whereas plate inclination enables necessary shear for producing slurry of desired solid fraction. A numerical model concerning transport equations of mass, momentum, energy and species is developed for predicting velocity, temperature, macrosegregation and solid fraction. The model uses FVM with phase change algorithm, VOF and variable viscosity. The model introduces solid phase movement with gravity effect as well. Effects of melt pouring temperature and plate inclination on hydrodynamic and thermo-solutal behaviors are studied subsequently. Slurry solid fractions at plate exit are 27%, 22%, 16%, and 10% for pouring temperatures of 620 degrees C, 625 degrees C, 630 degrees C, and 635 degrees C, respectively. And, are 27%, 25%, 22%, and 18% for plate inclinations of 30, 45, 60, and 75, respectively. Melt pouring temperature of 625 degrees C with plate inclination of 60 generates appropriate quality of slurry and is the optimum. Both numerical and experimental results are in good agreement with each other. (C) 2015 Taiwan Institute of Chemical Engineers. Published by Elsevier B.V. All rights reserved.

Three-dimensional numerical analysis of geocell-reinforced soft clay beds by considering the actual geometry of geocell pockets

Due to its complex honeycomb structure, the numerical modeling of the geocell has always been a big challenge. Generally, the equivalent composite approach is used to model the geocells. In the equivalent composite approach, the geocellsoil composite is treated as the soil layer with improved strength and stiffness values. Though this approach is very simple, it is unrealistic to model the geocells as the soil layer. This paper presents a more realistic approach of modeling the geocells in three-dimensional (3D) framework by considering the actual curvature of the geocell pocket. A square footing resting on geocell reinforced soft clay bed was modeled using the ``fast Lagrangian analysis of continua in 3D'' (FLAC(3D)) finite difference package. Three different material models, namely modified Cam-clay, Mohr-Coulomb, and linear elastic were used to simulate the behaviour of foundation soil, infill soil and the geocell, respectively. It was found that the geocells distribute the load laterally to the wider area below the footing as compared to the unreinforced case. More than 50% reduction in the stress was observed in the clay bed in the presence of geocells. In addition to geocells, two other cases, namely, only geogrid and geocell with additional basal geogrid cases were also simulated. The numerical model was systematically validated with the results of the physical model tests. Using the validated numerical model, parametric studies were conducted to evaluate the influence of various geocell properties on the performance of reinforced clay beds.

Effects of the filling process on the evolution of the mushy zone and macrosegregation in alloy casting

A numerical model of the entire casting process starting from the mould filling stage to complete solidification is presented. The model takes into consideration any phase change taking place during the filling process. A volume of fluid method is used for tracking the metal–air interface during filling and an enthalpy based macro-scale solidification model is used for the phase change process. The model is demonstrated for the case of filling and solidification of Pb–15 wt%Sn alloy in a side-cooled two-dimensional rectangular cavity, and the resulting evolution of a mushy region and macrosegregation are studied. The effects of process parameters related to filling, namely degree of melt superheat and filling velocity on macrosegregation in the cavity, are also investigated. Results show significant differences in the progress of the mushy zone and macrosegregation pattern between this analysis and conventional analysis without the filling effect.

Studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring

Numerical and experimental studies on transport phenomena during solidification of an aluminum alloy in the presence of linear electromagnetic stirring are performed. The alloy is electromagnetically stirred to produce semisolid slurry in a cylindrical graphite mould placed in the annulus of a linear electromagnetic stirrer. The mould is cooled at the bottom, such that solidification progresses from the bottom to the top of the cylindrical mould. A numerical model is developed for simulating the transport phenomena associated with the solidification process using a set of single-phase governing equations of mass. momentum, energy. and species conservation. The viscosity variation of the slurry, used in the model, is determined experimentally using a rotary viscometer. The set of governing equations is solved using a pressure-based finite volume technique, along with an enthalpy based phase change algorithm. The numerical study involves prediction of temperature, velocity, species and solid fraction distribution in the mould. Corresponding solidification experiments are performed, with time-temperature history recorded at key locations. The microstructures at various temperature measurement locations in the solidified billet are analyzed. The numerical predictions of temperature variations are in good agreement with experiments, and the predicted flow field evolution correlates well with the microstructures observed at various locations.

Effect of reinforcement form on the bearing capacity of square footings on sand

This paper presents the results of laboratory model loading tests and numerical studies carried out on square footings supported on geosynthetic reinforced sand beds. The relative performance of different forms of geosynthetic reinforcement (i.e. geocell, planar layers and randomly distributed mesh elements) in foundation beds is compared; using same quantity of reinforcement in each test. A biaxial geogrid and a geonet are used for reinforcing the sand beds. Geonet is used in two forms of reinforcement, viz. Planar layers and geocell, while the biaxial geogrid was used in three forms of reinforcement, viz. planar layers, geocell and randomly distributed mesh elements. Laboratory load tests on unreinforced and reinforced footings are simulated in a numerical model and the results are analyzed to understand the distribution of displacements and stresses below the footing better. Both the experimental and numerical studies demonstrated that the geocell is the most advantageous form of soil reinforcement technique of those investigated, provided there is no rupture of the material during loading. Geogrid used in the form of randomly distributed mesh elements is found to be inferior to the other two forms. Some significant observations on the difference in reinforcement mechanism for different forms of reinforcement are presented in this paper.

Centrifugal pump transient characteristics and analysis using the method of characteristics

A theoretical and experimental study has been carried out on the transient characteristics of a centrifugal pump during starting and stopping periods. Experiments have been conducted on a volute pump with different valve openings to study the dynamic behaviour of the pump during normal start up and stopping, when a small length of discharge pipe line is connected to discharge flange of the pump. Similar experiments have also been conducted when the test pump was part of a hydraulic system to study the system effect on the transient characteristics. Instantaneous rotational speed, flowrate, and delivery and suction pressures of the pump are recorded and it is observed in ail the tested cases that the change of pump behaviour during the transient period is quasi-steady. The dynamic characteristics of the pump have been analysed by a numerical model using the method of characteristics. The model is presented and the results are compared with the experimental data. As the model contains speed acceleration and unsteady discharge terms, the model can be applied for analyses of purely unsteady cases where the pump dynamic characteristics show considerable departure from their steady-state characteristics.

The energy balance in the Ramdas layer

On calm clear nights, air at a height of a few decimetres above bare soil can be cooler than the surface by several degrees in what we shall call the Ramdas layer (Ramdas and Atmanathan, 1932). The authors have recently offered a logical explanation for such a lifted temperature minimum, together with a detailed numerical model. In this paper, we provide physical insight into the phenomenon by a detailed discussion of the energy budget in four typical cases, including one with a lifted minimum. It is shown that the net cooling rate near ground is the small difference between two dominant terms, representing respectively radiative upflux from the ground and from the air layers just above ground. The delicate energy balance that leads to the lifted minimum is upset by turbulent transport, by surface emissivity approaching unity, or by high ground cooling rates. The rapid variation of the flux emissivity of humid air is shown to dominate radiative transport near the ground.

Three-dimensional computational modeling of momentum, heat, and mass transfer in a laser surface alloying process

A three- dimensional, transient model is developed for studying heat transfer, fluid flow, and mass transfer for the case of a single- pass laser surface alloying process. The coupled momentum, energy, and species conservation equations are solved using a finite volume procedure. Phase change processes are modeled using a fixed-grid enthalpy-porosity technique, which is capable of predicting the continuously evolving solid- liquid interface. The three- dimensional model is able to predict the species concentration distribution inside the molten pool during alloying, as well as in the entire cross section of the solidified alloy. The model is simulated for different values of various significant processing parameters such as laser power, scanning speed, and powder feedrate in order to assess their influences on geometry and dynamics of the pool, cooling rates, as well as species concentration distribution inside the substrate. Effects of incorporating property variations in the numerical model are also discussed.