Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

On some results of A. E. Livingston and coefficient problems for concave univalent functions

We consider functions that map the open unit disc conformally onto the complement of a bounded convex set. We call these functions concave univalent functions. In 1994, Livingston presented a characterization for these functions. In this paper, we observe that there is a minor flaw with this characterization. We obtain certain sharp estimates and the exact set of variability involving Laurent and Taylor coefficients for concave functions. We also present the exact set of variability of the linear combination of certain successive Taylor coefficients of concave functions.

Nature of Cl center dot center dot center dot Cl Intermolecular Interactions via Experimental and Theoretical Charge Density Analysis:Correlation of Polar Flattening Effects with Geometry

The experimental charge density distribution in three compounds, 2-chloro-3-quinolinyl methanol, 2-chloro-3-hydroxypyridine, and 2-chloro-3-chloromethyl-8-methylquinoline, has been obtained using high-resolution X-ray diffraction data collected at 100 K based on the aspherical multipole modeling of electron density. These compounds represent type I (cis), type I (trans), and type II geometries, respectively, as defined for short Cl center dot center dot center dot Cl interactions. The experimental results are compared with the theoretical charge densities using theoretical structure factors obtained from a periodic quantum calculation at the B3LYP/6-31G** level. The topological features derived from the Bader's theory of atoms in molecules (AIM) approach unequivocally suggest that both cis and trans type I geometries show decreased repulsion, whereas type II geometry is attractive based on the nature of polar flattening of the electron density around the Cl atom.

Numerical and experimental study of the thermal behaviour of coining and upsetting processes

A numerical simulation technique has been employed to study the thermal behavior of hot-forging type forming processes. Experiments on the coining and upsetting of an aluminum billet were conducted to validate the numerical predictions. Typical forming conditions for both the coining and upsetting processes were then studied in detail. an electrical analogy scheme was used to determine the thermal contact resistance. This scheme can conviniently provide the interface characteristics for typical processing conditions, which normally involve high pressures and temperatures. A single forging cycle was first considered, and then a batch of twenty-five forgings was studied. Each forging cycle includes the billet mounting, ascent, loading, dwelling, unloading, descent, and billet removal stages. The temperature distribution in the first forging to be formed is found to be significantly different from that at the end of the batch. In industry, forging is essentially a batch operation. The influence of forming speed and reduction on thermal characteristics was investigated also. The variations that can occur in the process design by considering differences in temperature characteristics are discussed also.

On extremal solutions to stochastic control problems, II

The set of attainable laws of the joint state-control process of a controlled diffusion is analyzed from a convex analytic viewpoint. Various equivalence relations depending on one-dimensional marginals thereof are defined on this set and the corresponding equivalence classes are studied.

A monolithic strategy for fluid-structure interaction problems

In this work, we present a new monolithic strategy for solving fluid-structure interaction problems involving incompressible fluids, within the context of the finite element method. This strategy, similar to the continuum dynamics, conserves certain properties, and thus provides a rational basis for the design of the time-stepping strategy; detailed proofs of the conservation of these properties are provided. The proposed algorithm works with displacement and velocity variables for the structure and fluid, respectively, and introduces no new variables to enforce velocity or traction continuity. Any existing structural dynamics algorithm can be used without change in the proposed method. Use of the exact tangent stiffness matrix ensures that the algorithm converges quadratically within each time step. An analytical solution is presented for one of the benchmark problems used in the literature, namely, the piston problem. A number of benchmark problems including problems involving free surfaces such as sloshing and the breaking dam problem are used to demonstrate the good performance of the proposed method. Copyright (C) 2010 John Wiley & Sons, Ltd.

Studies on helium flow in experimental loops with different impedances operated with fountain effect pumps

An experimental flow loop with He II flow driven by fountain effect pumps (FEPs) is studied with respect to operation at different flow impedances and with thermal loads applied at different positions. The measured values of temperature, flow rate and pressure drop are compared with calculations resulting from a simplified model which assumes ideal performance of the porous plug and of the heat exchangers and which does not take into account Gorter-Mellink (GM) conduction. The main features of the loop are shown to be well described by this model. Refined calculations with a more complex model, including GM conduction of the He II, are only required for predicting the temperature distribution in some discrete regions of the loop.

Time-domain-finite-wave analysis of the engine exhaust system by means of the stationary-frame method of characteristics. Part II. Computed results and experimental corroboration thereof

Time-domain-finite-wave analysis of engine exhaust systems is usually carried out by means of the method of characteristics. The theory and the computational details of the stationary-frame method have been worked out in the accompanying paper (part I). In this paper (part II), typical computed results are given and discussed. A setup designed for experimental corroboration is described. The results obtained from the simulation are found to be in good agreement with experimental observations.

Experimental Test of the Quantum No-Hiding Theorem

The no-hiding theorem says that if any physical process leads to bleaching of quantum information from the original system, then it must reside in the rest of the Universe with no information being hidden in the correlation between these two subsystems. Here, we report an experimental test of the no-hiding theorem with the technique of nuclear magnetic resonance. We use the quantum state randomization of a qubit as one example of the bleaching process and show that the missing information can be fully recovered up to local unitary transformations in the ancilla qubits.

Selectivity in Garratt-Braverman Cyclization: An Experimental and Computational Study

Bispropargyl sulfones equipped with aromatic rings of dissimilar nature were synthesized. Under basic conditions, these sulfones isomerized to the bisallenic sulfones, creating a competitive scenario between two alternate Garratt-Braverman (GB) cyclization pathways. The observed product distribution ruled out the involvement of any ionic intermediate and supported the diradical mechanism with greater involvement of the electron-rich aromatic ring via the more nucleophilic radical. DFT-based calculations supported the diradical mechanism along with the observed selectivity.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Diastereospecific Coupling of Imines by Low-Valent Titanium: An Experimental and Computational Study

The reaction of phenylsilane (PhSiH3) and titanium(IV) isopropoxide Ti(OiPr)(4)] generates low-valent titanium alkoxides that reduce and reductively couple imines. The C-C coupling reaction is diastereospecific, with exclusive formation of the (+/-)-isomer. The yield is dependent on the concentration of titanium(IV) isopropoxide used. Using imines with varying electronic demand revealed that the coupling is most efficient for unsubstituted imines. The involvement of a trimeric titanium biradical species is invoked to explain the observed concentration dependence and diastereospecificity in the reaction. Multilayer ONIOM (''our own n-layered integrated molecular orbital and molecular mechanics'') calculations were carried out on the plausible intermediates involved by using the Guassian 09 suite of programs. Formation of the trimeric biradical species leading to the formation of the (+/-)-isomer is more favored than formation of the trimeric biradical species, which leads to the meso-isomer. Similar calculations with dimeric intermediates do not predict the (+/-)-selectivity. Mass spectral analysis of the reaction mixture indicates the presence of a trimeric titanium alkoxide species.

Particle filters for structural system identification using multiple test and sensor data: A combined computational and experimental study

The problem of structural system identification when measurements originate from multiple tests and multiple sensors is considered. An offline solution to this problem using bootstrap particle filtering is proposed. The central idea of the proposed method is the introduction of a dummy independent variable that allows for simultaneous assimilation of multiple measurements in a sequential manner. The method can treat linear/nonlinear structural models and allows for measurements on strains and displacements under static/dynamic loads. Illustrative examples consider measurement data from numerical models and also from laboratory experiments. The results from the proposed method are compared with those from a Kalman filter-based approach and the superior performance of the proposed method is demonstrated. Copyright (C) 2009 John Wiley & Sons, Ltd.