Discrete multicomponent model for biodiesel spray combustion simulation

An important first step in spray combustion simulation is an accurate determination of the fuel properties which affects the modelling of spray formation and reaction. In a practical combustion simulation, the implementation of a multicomponent model is important in capturing the relative volatility of different fuel components. A Discrete Multicomponent (DM) model is deemed to be an appropriate candidate to model a composite fuel like biodiesel which consists of four components of fatty acid methyl esters (FAME). In this paper, the DM model is compared with the traditional Continuous Thermodynamics (CTM) model for both diesel and biodiesel. The CTM model is formulated based on mixing rules that incorporate the physical and thermophysical properties of pure components into a single continuous surrogate for the composite fuel. The models are implemented within the open-source CFD code OpenFOAM, and a semi-quantitative comparison is made between the predicted spray-combustion characteristics and optical measurements of a swirl-stabilised flame of diesel and biodiesel. The DM model performs better than the CTM model in predicting a higher magnitude of heat release rate in the top flame brush region of the biodiesel flame compared to that of the diesel flame. Using both the DM and CTM models, the simulation successfully reproduces the droplet size, volume flux, and droplet density profiles of diesel and biodiesel. The DM model predicts a longer spray penetration length for biodiesel compared to that of diesel, as seen in the experimental data. Also, the DM model reproduces a segregated biodiesel fuel vapour field and spray in which the most abundant FAME component has the longest vapour penetration. In the biodiesel flame, the relative abundance of each fuel component is found to dominate over the relative volatility in terms of the vapour species distribution and vice versa in the liquid species distribution. © 2014 Elsevier Ltd. All rights reserved.

Minimalistic models of an energy-efficient vertical-hopping robot

The use of free vibration in elastic structure can lead to energy-efficient robot locomotion, since it significantly reduces the energy expenditure if properly designed and controlled. However, it is not well understood how to harness the dynamics of free vibration for the robot locomotion, because of the complex dynamics originated in discrete events and energy dissipation during locomotion. From this perspective, the goals of this paper are to propose a design strategy of hopping robot based on elastic curved beams and actuated rotating masses and to identify the minimalistic model that can characterize the basic principle of robot locomotion. Since the robot mainly exhibits vertical hopping, three 1-D models are examined that contain different configurations of simple spring-damper-mass components. The real-world and simulation experiments show that one of the models best characterizes the robot hopping, through analyzing the basic kinematics and negative works in actuation. Based on this model, the self-stability of hopping motion under disturbances is investigated, and design and control parameters are analyzed for the energy-efficient hopping. In addition, further analyses show that this robot can achieve the energy-efficient hopping with the variation in payload, and the source of energy dissipation of the robot hopping is investigated. © 1982-2012 IEEE.

Minimalistic models of an energy efficient vertical hopping robot

The use of free vibration in elastic structure can lead to energy efficient robot locomotion, since it significantly reduces the energy expenditure if properly designed and controlled. However, it is not well understood how to harness the dynamics of free vibration for the robot locomotion, because of the complex dynamics originated in discrete events and energy dissipation during locomotion. From this perspective, this paper explores three minimalistic models of free vibration that can characterize the basic principle of robot locomotion. Since the robot mainly exhibits vertical hopping, three one-dimensional models are examined that contain different configurations of simple spring-damper-mass components. The self-stability of these models are also investigated in simulation. The real-world and simulation experiments show that one of the models best characterizes the robot hopping, through analyzing the basic kinematics and negative works in actuation. Based on this model, the control parameters are analyzed for the energy efficient hopping. © 2013 IEEE.

Simulation on Gain Recovery of Quantum Dot Semiconductor Optical Amplifiers by Rate Equation

The gain recoveries in quantum dot semiconductor optical amplifiers are numerically studied by rate equation models. Similar to the optical pump-probe experiment, the injection of double optical pulses is used to simulate the gain recovery of a weak continuous signal for the QD SOAs. The gain recoveries are fitted by a response function with multiple exponential terms. For the pulses duration of 10 ps, the gain recovery can be described by three exponential terms with the time constants, and for the pulse with the width of 150 fs, the gain recovery can be described by two exponential terms, the reason is that the short pulse does not consume lot of carriers.

simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.

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.

SIMULATION OF LATERAL CONFINEMENT IN VERY NARROW CHANNELS

A theoretical study is presented of the lateral confinement potential (CP) in the very narrow mesa channels fabricated in the conventional two-dimensional (2D) electron gas in GaAs-AlxGa1-xAs heterostructures. The ID electronic structures are calculated in the framework of the confinement potential: V(x) = m* omega0(2)x2/2 for Absolute value of x models consisting of flat bottoms and soft walls are more realistic than this complex potential. The experimental result for a wire of W(str) = 250 nm cannot be explained within the composite-well model including these two cases. Thus how to explain the experiment for this wire remains an open question.

THE STUDY OF MICROSTRUCTURE IN IRON-DOPED AND SULFUR-DOPED INP CRYSTALS BY MEANS OF HREM AND COMPUTER-SIMULATION

The microstructures in iron- and sulphur-doped InP crystals were studied using both electron microscopy and electron diffraction. A modulated structure has been found in S-doped InP crystal, where the commensurate modulations corresponded to periodicities of 0.68 nm and 0.7 nm in real space and were related to the reflections of the cubic lattice in [111] and [113BAR] directions; they were indexed as q111* = 1/2(a* + b* + c*) and q113BAR* = 1/4(-a* - b* + 3c*), respectively. Single atomic layers of iron precipitate were observed, with preferred orientations along which precipitates are formed. Simulated calculations by means of the dynamical theory of electron diffraction using models for the precipitate structure were in good agreement with our experimental results. The relation between the modulated structure and the precipitates is also discussed.

Numerical simulation of thermal loading produced by shaped high power laser onto engine parts

Recently a new method for simulating the thermal loading on pistons of diesel engines was reported. The spatially shaped high power laser is employed as the heat source, and some preliminary experimental and numerical work was carried out. In this paper, a further effort was made to extend this simulation method to some other important engine parts such as cylinder heads. The incident Gaussian beam was transformed into concentric multi-circular patterns of specific intensity distributions, with the aid of diffractive optical elements (DOEs). By incorporating the appropriate repetitive laser pulses, the designed transient temperature fields and thermal loadings in the engine parts could be simulated. Thermal-structural numerical models for pistons and cylinder heads were built to predict the transient temperature and thermal stress. The models were also employed to find the optimal intensity distributions of the transformed laser beam that could produce the target transient temperature fields. Comparison of experimental and numerical results demonstrated that this systematic approach is effective in simulating the thermal loading on the engine parts. (C) 2009 Elsevier Ltd. All rights reserved.

CFD study of a sudden-expanding coal combustor using Euler-Euler and Euler-Lagrange models

In this study, the Euler-Euler (E-E) and Euler-Lagrange (E-L) models designed for the same chemical mechanism of heterogeneous reactions were used to predict the performance of a typical sudden-expanding coal combustor. The results showed that the current E-E model underestimated the coal burnout rate because the particle temperature fluctuation on char combustion is not adequately considered. A comparison of the E-E and E-L simulations showed the underestimation of heterogeneous chemical reaction rates by the E-E model. (C) 2010 Elsevier Ltd. All rights reserved.

CFD-DEM simulation of three-dimensional aeolian sand movement

A three-dimensional CFD-DEM model is proposed to investigate the aeolian sand movement. The results show that the mean particle horizontal velocity can be expressed by a power function of heights. The probability distribution of the impact and lift-off velocities of particles can be described by a log-normal function, and that of the impact and lift-off angles can be expressed by an exponential function. The probability distribution of particle horizontal velocity at different heights can be described as a lognormal function, while the probability distribution of longitudinal and vertical velocity can be described as a normal function. The comparison with previous two-dimensional calculations shows that the variations of mean particle horizontal velocity along the heights in two-dimensional and three-dimensional models are similar. However, the mean particle density of the two-dimensional model is larger than that in reality, which will result in the overestimation of sand transportation rate in the two-dimensional calculation. The study also shows that the predicted probability distributions of particle velocities are in good agreement with the experimental results.

Simulation study of China's net primary production

Spatial and temporal distribution of vegetation net primary production (NPP) in China was studied using three light-use efficiency models (CASA, GLOPEM and GEOLUE) and two mechanistic ecological process models (CEVSA, GEOPRO). Based on spatial and temporal analysis (e.g. monthly, seasonally and annually) of simulated results from ecological process mechanism models of CASA, GLOPEM and CEVSA, the following conclusions could be made: (1) during the last 20 years, NPP change in China followed closely the seasonal change of climate affected by monsoon with an overall trend of increasing; (2) simulated average seasonal NPP was: 0.571 +/- 0.2 GtC in spring, 1.573 +/- 0.4 GtC in summer, 0.6 +/- 0.2 GtC in autumn, and 0.12 +/- 0.1 GtC in winter. Average annual NPP in China was 2.864 +/- 1 GtC. All the five models were able to simulate seasonal and spatial features of biomass for different ecological types in China. This paper provides a baseline for China's total biomass production. It also offers a means of estimating the NPP change due to afforestation, reforestation, conservation and other human activities and could aid people in using for-mentioned carbon sinks to fulfill China's commitment of reducing greenhouse gases.

Computer simulation of the damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects

The damage evolution of fiber-reinforced polypropylene-matrix composites with matrix defects was studied via a Monte Carlo technique combined with a finite element method. A finite element model was constructed to predict the effects of various matrix defect shapes on the stress distributions. The results indicated that a small matrix defect had almost no effect on fiber stress distributions other than interfacial shear stress distributions. Then, a finite element model with a statistical distribution of the fiber strength was constructed to investigate the influences of the spatial distribution and the volume fraction of matrix defects on composite failure. The results showed that it was accurate to use the shear-lag models and Green's function methods to predict the tensile strength of composites even though the axial stresses in the matrix were neglected.

Effects of chain flexibility on polymer conformation in dilute solution studied by lattice Monte Carlo simulation

Effects of chain flexibility on the conformation of homopolymers in good solvents have been investigated by Monte Carlo simulation. Bond angle constraint coupled with persistence length of polymer chains has been introduced in the modified eight-site bond fluctuation simulation model. The study about the effects of chain flexibility on polymer sizes reveals that the orientation of polymer chains under confinement is driven by the loss of conformation entropy. The conformation of polymer chains undergoing a gradual change from spherical iso-diametric ellipsoid to rodlike iso-diametric ellipsoid with the decrease of polymer chain flexibility in a wide region has been clearly illustrated from several aspects. Furthermore, a comparison of the freely jointed chain (FJC) model and the wormlike chain (WLC) model has also been made to describe the polymer sizes in terms of chain flexibility and quasi-quantitative boundary toward the suitability of the models.