A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

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.

Brownian dynamics simulation of dense binary colloidal mixtures. I. Structural evolution and dynamics

We have carried out Brownian dynamics simulations of binary mixtures of charged colloidal suspensions of two different diameter particles with varying volume fractions phi and charged impurity concentrations n(i). For a given phi, the effective temperature is lowered in many steps by reducing n(i) to see how structure and dynamics evolve. The structural quantities studied are the partial and total pair distribution functions g(tau), the static structure factors, the time average g(<(tau)over bar>), and the Wendt-Abraham parameter. The dynamic quantity is the temporal evolution of the total meansquared displacement (MSD). All these parameters show that by lowering the effective temperature at phi = 0.2, liquid freezes into a body-centered-cubic crystal whereas at phi = 0.3, a glassy state is formed. The MSD at intermediate times shows significant subdiffusive behavior whose time span increases with a reduction in the effective temperature. The mean-squared displacements for the supercooled liquid with phi = 0.3 show staircase behavior indicating a strongly cooperative jump motion of the particles.

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.

Geometry of sand-bed channels with seepage

Literature suggests that apart from a significant percentage of water loss from seepage, downward seepage from the channel causes an increase in the mobility of channel bed materials and thus changes the channel stability. Consequently, regime conditions (which provide the relationship among hydraulic parameters) should also be revised by incorporating downward seepage as an additional parameter. In the present work, regime conditions for the prediction of channel geometry in alluvial channels affected by downward seepage have been formulated on the basis of experimental observations. (C) 2011 Elsevier B.V. All rights reserved.

On the dual head-discharge characteristics of a modified chimney weir

This paper is concerned with the dual head-discharge characteristics of a modified chimney weir. It is shown by an optimization procedure that the modified chimney weir having an inward trapezoidal weir over an inverted V-notch can produce discharges proportional to both the logarithm of the head as well as the linear power of the head reckoned independently over two different reference planes, within a prefixed maximum allowable percentage error from the theoretical discharge. A new technique is adopted to optimize the weir parameters, in order to obtain the maximum ranges of measurement under logarithmic, linear as well as combined characteristics. In the case of linear weir it is shown that it is possible to enhance the linearity range of the chimney weir by more than 540% and for a weir with constant indication accuracy by more than 350%. In addition, about 86% of the overall depth of the designed linear weir and over 90% of the logarithmic weir is converted as the corresponding measurable ranges. Experiments with four typical weirs give consistent constant average coefficient of discharge for each weir confirming the theory. The practical application of the weir in minor irrigation, hydraulic and other engineering fields is highlighted.

Performance of a mixed flow pump with varying tip clearance .2.

Measurements in a mixed flow pump of non-dimensional specific speed k = 1.89[N-S = 100 r/min (metric)] are analysed to give loss distribution and local hydraulic efficiencies at different flowrates and values of tip clearance. Fairly close agreement is obtained between the relative flow angles leaving the blading as predicted by simple deviation and slip models and derived from the measurements. The head developed is broken up into two parts: that contributed by Coriolis action and that associated with blade circulation. It is suggested that lift coefficients based on blade circulation are of limited value in selecting blade profiles. The variation of pump efficiency with tip clearance is greater than that reported for centrifugal pumps.

Spin and charge density waves in the extended Hubbard model: A slave-boson approach

We use the extended Hubbard model to investigate the properties of the charge- and spin-density-wave phases in the presence of a nearest-neighbors repulsion term in the framework of the slave-boson technique. We show that, contrary to Hartree-Fock results, an instablity may occur for sufficiently high values of the Hubbard repulsion, both in the spin- and charge-density-wave phase, which makes the system discontinuously jump to a phase with a smaller or zero wave amplitude. The limits of applicability of our approach are discussed and our results are compared with previous numerical analysis. The phase diagram of the model at half-filling is determined.

Nematogens with more flexible chains than aromatic rings in the core

Mesogens containing four rings in the main core can accommodate one terminal and two nearby lateral chains on each outside aromatic ring. These compounds containing six chains present an enantiotropic nematic range which is influenced by the rigidity of the links. The conformational behaviour of the first methyleneoxy group within the chains was investigated by one and two dimensional C-13 NMR. The sign of the jump in chemical shifts when entering the nematic phase indicates the folding of each lateral branch. Dipolar oscillations during cross-polarization contact provide the values of the bond order parameter. The two First lateral fragments do not behave in the same way, demonstrating the influence of the fragment along which the chain is back: folded.

Simplified seepage consolidation test for soft sediments

The consistency of very soft sediments prevents the conventional oedometer test from being applied to study their compressibility and permeability characteristics. The hydraulic consolidation test in existence requires sophisticated instrumentation and testing procedures. The present paper proposes a seepage-force-induced consolidation testing procedure for studying the compressibility and permeability behavior of soft sediments at low effective stress levels. The good agreement that has been observed between the results obtained from the proposed method and the conventional oedometer test at overlapping effective stress levels indicates that the proposed method can be used to study the compressibility and permeability characteristics of soft sediments at low effective stress levels satisfactorily.

Energy transfer efficiency distributions in polymers in solution during folding and unfolding

Distribution of fluorescence resonance energy transfer (FRET) efficiency between the two ends of a Lennard-Jones polymer chain both at equilibrium and during folding and unfolding has been calculated, for the first time, by Brownian dynamics simulations. The distribution of FRET efficiency becomes bimodal during folding of the extended state subsequent to a temperature quench, with the width of the distribution for the extended state broader than that for the folded state. The reverse process of unfolding subsequent to a upward temperature jump shows different characteristics. The distributions show significant viscosity dependence which can be tested against experiments.

Phase transformation and semiconductor-metal transition in thin films of VO2 deposited by low-pressure metalorganic chemical vapor deposition

Thin films of the semiconducting, monoclinic vanadium dioxide, VO2(M) have been prepared on ordinary glass by two methods: directly by low-pressure metalorganic chemical vapor deposition (MOCVD), and by argon-annealing films of the VO2(B) phase deposited by MOCVD. The composition and microstructure of the films have been examined by x-ray diffraction, scanning electron microscopy, and transmission electron microscopy. Films made predominantly of either the B or the M phase, as deposited, can only be obtained over a narrow range of deposition temperatures. At the lower end of this temperature range, the as-deposited films are strongly oriented, although the substrate is glass. This can be understood from the drive to minimize surface energy. Films of the B phase have a platelet morphology, which leads to an unusual microstructure at the lower-deposition temperatures. Those grown at similar to370 degreesC convert to the metallic, rutile (R) phase when annealed at 550 degreesC, whereas those deposited at 420 degreesC transform to the R phase only at 580 degreesC. (When cooled to room temperature, the annealed films convert reversibly from the R phase to the M phase.) Electron microscopy shows that annealing leads to disintegration of the single crystalline VO2(B) platelets into small crystallites of VO2(R), although the platelet morphology is retained. When the annealing temperature is relatively low, these crystallites are nanometer sized. At a higher-annealing temperature, the transformation leads to well-connected and similarly oriented large grains of VO2(R), enveloped in the original platelet. The semiconductor-metal transition near 68 degreesC leads to a large jump in resistivity in all the VO2(M) films, nearly as large as in epitaxial films on single-crystal substrates. When the annealed films contain well-connected large grains, the transition is very sharp. Even when preferred orientation is present, the transition is not as sharp in as-deposited VO2(M), because the crystallites are not densely packed as in annealed VO2(B). However, the high degree of orientation in these films leads to a narrow temperature hysteresis. (C) 2002 American Institute of Physics.

Time Evolution of the Microstructure of VO2(B) Films Deposited on Glass by MOCVD

Thin films of VO2(B), a metastable polymorph of vanadium dioxide, have been grown on glass by low-pressure metalorganic chemical vapor deposition (MOCVD). The films grown for 90 minutes have atypical microstructure, comprising micrometer-sized, island-like entities made up of numerous small, single-crystalline platelets (≅1 μm) emerging orthogonally from larger ones at the center. Microstructure evolution as a function of deposition time has been examined by X-ray diffraction, scanning electron microscopy, and transmission electron microscopy. The metastable VO2(B) transforms to the stable rutile (R) phase at 550°C in inert ambient, which on cooling convert reversibly to M phase. Electron microscopy shows that annealing leads to the disintegration of the VO2(B) platelets into small crystallites of the rutile phase VO2(R), although the platelet morphology is retained. The magnitude of the jump in resistance at the semiconductor-to-metal, VO2(M)→VO2(R) phase transition depends on the arrangement of polycrystalline platelets in the films.

Analytical solutions for protective filters based on soil-retention and permeability criteria with respect to the phenomenon of soil boiling

For the successful performance of a granular filter medium, existing design guidelines, which are based on the particle size distribution (PSD) characteristics of the base soil and filter medium, require two contradictory conditions to be satisfied, viz., soil retention and permeability. In spite of the wider applicability of these guidelines, it is well recognized that (i) they are applicable to a particular range of soils tested in the laboratory, (ii) the design procedures do not include performance-based selection criteria, and (iii) there are no means to establish the sensitivity of the important variables influencing performance. In the present work, analytical solutions are developed to obtain a factor of safety with respect to soil-retention and permeability criteria for a base soil - filter medium system subjected to a soil boiling condition. The proposed analytical solutions take into consideration relevant geotechnical properties such as void ratio, permeability, dry unit weight, effective friction angle, shape and size of soil particles, seepage discharge, and existing hydraulic gradient. The solution is validated through example applications and experimental results, and it is established that it can be used successfully in the selection as well as design of granular filters and can be applied to all types of base soils.