Modeling and Simulation of Liquid Pulsed Particulate Fluidized Beds

Pulsed fluidization is of considerable interest in process engineering for improving fluidization quality. Quantitative understanding of the pulsed two-phase flow behaviors is very important for proper design and optimum operation of such contactors. The

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

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.

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by h

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

The Deposition and Burning Characteristics During Slagging Co-Firing Coal and Wood: Modeling and Numerical Simulation

Numerical analysis was used to study the deposition and burning characteristics of combining co-combustion with slagging combustion technologies in this paper. The pyrolysis and burning kinetic models of different fuels were implanted into the WBSF-PCC2 (wall burning and slag flow in pulverized co-combustion) computation code, and then the slagging and co-combustion characteristicsespecially the wall burning mechanism of different solid fuels and their effects on the whole burning behavior in the cylindrical combustor at different mixing ratios under the condition of keeping the heat input samewere simulated numerically. The results showed that adding wood powder at 25% mass fraction can increase the temperature at the initial stage of combustion, which is helpful to utilize the front space of the combustor. Adding wood powder at a 25% mass fraction can increase the reaction rate at the initial combustion stage; also, the coal ignitability is improved, and the burnout efficiency is enhanced by about 5% of suspension and deposition particles, which is helpful for coal particles to burn entirely and for combustion devices to minimize their dimensions or sizes. The results also showed that adding wood powder at a proper ratio is helpful to keep the combustion stability, not only because of the enhancement for the burning characteristics, but also because the running slag layer structure can be changed more continuously, which is very important for avoiding the abnormal slag accumulation in the slagging combustor. The theoretic analysis in this paper proves that unification of co-combustion and slagging combustion technologies is feasible, though more comprehensive and rigorous research is needed.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone). V: Estimation of temperature distribution during cure process

A numerical method to estimate temperature distribution during the cure of epoxy-terminated poly(phenylene ether ketone) (E-PEK)-based composite is suggested. The effect of the temperature distribution on the selection of cure cycle is evaluated using a suggested alternation criterion. The effect of varying heating rate and thickness on the temperature distribution, viscosity distribution and distribution of the extent of cure reaction are discussed based on the combination of the here-established temperature distribution model and the previously established curing kinetics model and chemorheological model. It is found that, for a thin composite (<=10mm) and low heating rate (<=2.5K/min), the effect of temperature distribution on cure cycle and on the processing window for pressure application can be neglected. Low heating rate is of benefit to reduce the temperature gradient. The processing window for pressure application becomes narrower with increasing thicknesses of composite sheets. The validity of the temperature distribution model and the modified processing window is evaluated through the characterization of mechanical and physical properties of E-PEK-based composite fabricated according to different temperature distribution conditions.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone) .3. Determination of the time of pressure application

The curing temperature, pressure, and curing time have significant influence on finished thermosetting composite products. The time of pressure application is one of the most important processing parameters in the manufacture of a thermosetting composite. The determination of the time of pressure application relies on analysis of the viscosity variation of the polymer, associated with curing temperature and curing time. To determine it, the influence of the time of pressure application on the physical properties of epoxy-terminated poly(phenylene ether ketone) (E-PEK)-based continuous carbon fiber composite was studied. It was found that a stepwise temperature cure cycle is more suitable for manufacture of this composite. There are two viscosity valleys, in the case of the E-PEK system, associated with temperature during a stepwise cure cycle. The analysis on the effects of reinforcement fraction and defect content on the composite sheet quality indicates that the width-adjustable second viscosity valley provides a suitable pressing window. The viscosity, ranging from 400 to 1200 Pa . s at the second viscosity valley, is the optimal viscosity range for applying pressure to ensure appropriate resin flow during curing process, which enables one to get a finished composite with optimal fiber volume fraction and low void content. (C) 1997 John Wiley & Sons, Inc.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone) .4. Cure cycle simulation

Epoxy-terminated poly(phenylene ether ketone) (E-PEK) developed in this Institute is a candidate matrix resin for polymer composites as structural materials. Cure cycles for this reaction system were simulated according to the previously established processing model. It is found that for the E-PEK system, the curing process is best completed by a stepwise cure cycle comprising two isothermal processes at different temperatures, T-1 and T-2. The cure cycles over a wide range of processing parameters simulated, based on the established processing model, indicate that the processing window is width-adjustable. Analysis of the mechanical properties of the composite sheets showed that the simulated cure cycles are acceptable and reliable. (C) 1997 John Wiley & Sons, Inc.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone) .1. DSC characterization of curing reaction

The curing reaction process of epoxy-terminated poly(phenylene ether ketone) (E-PEK) with 4,4'-diaminodiphenyl sulfone (DDS) and hexahydrophthalic acid anhydride (Nadic) as curing agents was investigated using isothermal differential scanning calorimetry (IDSC) and nonisothermal differential scanning calorimetry (DDSC) techniques. It was found that the curing reactions of E-PEK/DDS and E-PEK/Nadic are nth-order reactions but not autoaccelerating. The experimental results revealed that the curing reaction kinetics parameters measured from IDSC and DDSC are not equivalent. This means that, in the curing reaction kinetics model for our E-PEK system, both isothermal and nonisothermal reaction kinetics parameters are needed to describe isothermal and nonisothermal curing processes, The isothermal and nonisothermal curing processes were successfully simulated using this model. A new extrapolation method was suggested. On the basis of this method the maximum extent of the curing reaction (A(ult)) that is able to reach a certain temperature can be predicted. The A(ult) for the E-PEK system estimated by the new method agrees well with the results obtained from another procedure reported in the literature. (C) 1997 John Wiley & Sons, Inc.

Cure processing modeling and cure cycle simulation of epoxy-terminated poly(phenylene ether ketone) .2. Chemorheological modeling

Chemorheology and corresponding models for an epoxy-terminated poly(phenylene ether ketone) (E-PEK) and 4,4'-diaminodiphenyl sulfone (DDS) system were investigated using a differential scanning calorimeter (DSC) and a cone-and-plate rheometer. For this system, the reported four-parameter chemorheological model and modified WLF chemorheological model can only be used in an isothermal or nonisothermal process, respectively. In order to predict the resin viscosity variation during a stepwise temperature cure cycle actually used, a new model based on the combination of the four-parameter model and the modified WLF model was developed. The combined model can predict the resin viscosity variation during a stepwise temperature cure cycle more accurately than the above two models. In order to simplify the establishment of this model, a new five-parameter chemorheological model was then developed. The parameters in this five-parameter model can be determined through very few rheology and DSC experiments. This model is practicable to describe the resin viscosity variation for isothermal, nonisothermal, or stepwise temperature cure cycles accurately. The five-parameter chemorheological model has also successfully been used in the E-PEK systems with two other curing agents, i.e., the diamine curing agent with the addition of a boron trifluride monoethylamine (BF3-MEA) accelerator and an anhydride curing agent (hexahydrophthalic acid anhydride). (C) 1997 John Wiley & Sons, Inc.

Modeling and mesoscopic damage constitutive relation of brittle short-fiber-reinforced composites

Aimed at brittle composites reinforced by randomly distributed short-fibers with a relatively large aspect ratio, stiffness modulus and strength, a mesoscopic material model was proposed. Based on the statistical description, damage mechanisms, damage-induced anisotropy, damage rate effect and stress redistribution, the constitutive relation were derived. By taking glass fiber reinforced polypropylene polymers as an example, the effect of initial orientation distribution of fibers, damage-induced anisotropy, and damage-rate effect on macro-behaviors of composites were quantitatively analyzed. The theoretical predictions compared favorably with the experimental results.

Experimental study and simulation principles of an oil-gas multiphase transportation system

Presented is an experimental study on the performance of an oil-gas multiphase transportation system, especially on the multiphase flow patterns, multiphase pumping and multiphase metering of the system. A dynamic simulation analysis is conducted to deduce simulation parameters of the system and similarity criteria under simplified conditions are obtained. The reliability and feasibility of two-phase flow experiment with oil and natural gas simulated by water and air are discussed by using the similarity criteria.

Multiscale Coupling in Complex Mechanical Systems

Multiscale coupling attracts broad interests from mechanics, physics and chemistry to biology. The diversity and coupling of physics at different scales are two essential features of multiscale problems in far-from-equilibrium systems. The two features present fundamental difficulties and are great challenges to multiscale modeling and simulation. The theory of dynamical system and statistical mechanics provide fundamental tools for the multiscale coupling problems. The paper presents some closed multiscale formulations, e.g., the mapping closure approximation, multiscale large-eddy simulation and statistical mesoscopic damage mechanics, for two typical multiscale coupling problems in mechanics, that is, turbulence in fluids and failure in solids. It is pointed that developing a tractable, closed nonequilibrium statistical theory may be an effective approach to deal with the multiscale coupling problems. Some common characteristics of the statistical theory are discussed.

Behaviour analysis of numerical solutions and simulation of coherent structures in compressible mixing layers

For understanding the correctness of simulations the behaviour of numerical solutions is analysed, Tn order to improve the accuracy of solutions three methods are presented. The method with GVC (group velocity control) is used to simulate coherent structures in compressible mixing layers. The effect of initial conditions for the mixing layer with convective Mach number 0.8 on coherent structures is discussed. For the given initial conditions two types of coherent structures in the mixing layer are obtained.