Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Physical and numerical modeling of horizontal gas-hydrate wells

Methane hydrate, which is usually found under deep seabed or permafrost zones, is a potential energy resource for future years. Depressurization of horizontal wells bored in methane hydrate layer is considered as one possible method for hydrate dissociation and methane extraction from the hosting soil. Since hydrate is likely to behave as a bonding material to sandy soils, supported well construction is necessary to avoid well-collapse due to the loss of the apparent cohesion during depressurization. This paper describes both physical and numerical modeling of such horizontal support wells. The experimental part involves depressurization of small well models in a large pressure cell, while the numerical part simulates the corresponding problem. While the experiment models simulate only gas saturated initial conditions, the numerical analysis simulates both gas-saturated and more realistic water-saturated conditions based on effective stress coupled flow-deformation formulation of these three phases. © 2006 Taylor & Francis Group.

Numerical modeling of 3-D turbulent two-phase flow and coal combustion in a pulverized-coal combustor

In the present paper, a multifluid model of two-phase flows with pulverized-coal combustion, based on a continuum-trajectory model with reacting particle phase, is developed and employed to simulate the 3-D turbulent two-phase hows and combustion in a new type of pulverized-coal combustor with one primary-air jet placed along the wall of the combustor. The results show that: (1) this continuum-trajectory model with reacting particle phase can be used in practical engineering to qualitatively predict the flame stability, concentrations of gas species, possibilities of slag formation and soot deposition, etc.; (2) large recirculation zones can be created in the combustor, which is favorable to the ignition and flame stabilization.

Self-consistent kinetic modeling of electron distribution function in a low-pressure inductively coupled plasma

Based upon the spatially inhomogeneous Boltzmann equation in two-term approximation coupled with electromagnetic and fluid model analysis for the recently developed inductively coupled plasma sources, a self-consistent electron kinetic model is developed. The electron distribution function, spatial distributions of the electron density and ionization rate are calculated and discussed.

Mechanical modeling of grain boundary sliding of polycrystals

Using a variational method, a general three-dimensional solution to the problem of a sliding spherical inclusion embedded in an infinite anisotropic medium is presented in this paper. The inclusion itself is also a general anisotropic elastic medium. The interface is treated as a thin interface layer with interphase anisotropic properties. The displacements in the matrix and the inclusion are expressed as polynomial series of the cartesian coordinate components. Using the virtual work principle, a set of linear algebraic equations about unknown coefficients are obtained. Then the general sliding spherical inclusion problem is accurately solved. Based on this solution, a self-consistent method for sliding polycrystals is proposed. Combining this with a two-dimensional model of an aggregate polycrystal, a systematic analysis of the mechanical behaviour of sliding polycrystals is given in detail. Numerical results are given to show the significant effect of grain boundary sliding on the overall mechanical properties of aggregate polycrystals.

Modeling of a Laminar DC Arc Plasma Torch

In recent years, stable and long laminarplasma jets have been successfully generated, and thus it is possible to achieve low-noise working surroundings, better process repeatability and controllability, and reduced metal-oxidation degree in plasma materials processing. With such a recent development in thermal plasma science and technology as the main research background, modeling studies are performed concerning the DCarcplasmatorch for generating the long laminar argon plasma jet. Two different two-dimensional modeling approaches are employed to deal with the arc-root attachment at the anode surface. The first approach is based on circumferentially uniform arc-root attachment, while the second uses the so-called fictitious anode method. Modeling results show that the highest temperature and maximum axial-velocity at the plasmatorch exit are ～15000 K and ～1100 m/s, respectively, for the case with arc current of 160 A and argon flow rate of 1.95×10{sup}(-4)kg/s.

Physical Modeling of Current-Induced Seabed Scour around a Vibrating Submarine Pipeline

Most of the existing researches either focus on vortex-induced-vibrations (VIV) of a pipeline near a rigid boundary, or on seabed scour around a fixed pipeline. In the fields, pipeline vibration and seabed scour are actually always coupled. Based on the similarity analysis, a series of tests were conducted with a hydro-elastic facility to investigate the influence of pipe vibration on the local scour and the effects of scour process on the pipeline dynamic responses. Experimental results indicate that, there exist two phases in the process of sand scouring around the pipeline with small embedment, i.e. Phase I: scour beneath pipe without VIV, and Phase II: scour with VIV of pipe. It is also found that the gap-to-diameter ratio (e/D) has much effect upon the scour depth for the fixed pipes. For a given value of e/D, the vibrating pipes with close proximity to seabed may induce a deeper scour hole than the fixed ones. Within the examined gap-to-diameter ratio range (425 < e/D < 0.75), the influences of gap-to-diameter ratio on the maximum values of scour-depth for the case of vibrating pipes are not as much as those for the case of fixed pipes.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Establishing National Ocean Service Priorities for Estuarine, Coastal, and Ocean Modeling: Capabilities, Gaps, and Preliminary Prioritization Factors

This report was developed to help establish National Ocean Service priorities and chart new directions for research and development of models for estuarine, coastal and ocean ecosystems based on user-driven requirements and supportive of sound coastal management, stewardship, and an ecosystem approach to management. (PDF contains 63 pages)

A Featureless Approach to 3D Polyhedral Building Modeling from Aerial Images

This paper presents a model-based approach for reconstructing 3D polyhedral building models from aerial images. The proposed approach exploits some geometric and photometric properties resulting from the perspective projection of planar structures. Data are provided by calibrated aerial images. The novelty of the approach lies in its featurelessness and in its use of direct optimization based on image rawbrightness. The proposed framework avoids feature extraction and matching. The 3D polyhedral model is directly estimated by optimizing an objective function that combines an image-based dissimilarity measure and a gradient score over several aerial images. The optimization process is carried out by the Differential Evolution algorithm. The proposed approach is intended to provide more accurate 3D reconstruction than feature-based approaches. Fast 3D model rectification and updating can take advantage of the proposed method. Several results and evaluations of performance from real and synthetic images show the feasibility and robustness of the proposed approach.

Modeling smart grids as complex systems through the implementation of intelligent hubs

ICINCO 2010