基于LabVIEW的纳米颗粒动态散射光模拟

采用实验室虚拟仪器工程平台LabVIEW实现了纳米颗粒动态光散射信号的计算机模拟。用G语言(图形语言)设计了信号模拟的框图程序,给出了5nm, 15 nm两种粒径颗粒的模拟动态散射光信号及自相关函数。对模拟信号的粒度分析表明,这两种模拟信号产生的测量偏差分别为0.4nm和-0.6nm。     

Large eddy simulation of particle deposition and resuspension in turbulent duct flows

Particle deposition and resuspension in a horizontal, fully developed turbulent square duct flow at four flow bulk Reynolds numbers (10,320, 83k, 215k and 250k) is simulated by applying large eddy simulation coupled with Lagrangian particle tracking technique. Forces acting on particles includes drag, lift, buoyancy and gravity. Four particle sizes are considered with the diameters of 5?μm, 50?μm, 100?μm and 500?μm. Results obtained for the fluid phase are in good agreement with the available experimental and numerical data. Predictions for particles show that particle size, flow Reynolds number and the duct (celling, floor and vertical) walls play important roles in near-wall particle deposition and resuspension. For the smallest particle (5?μm), the particle deposition rates in duct ceiling, floor and vertical walls are found to be similar with each other and all increase with the flow Reynolds number, while the particle resuspension tends to occur in the middle wall regions and corners of the duct with less influenced by the flow Reynolds number. The ceiling deposition rate gradually decreases with particle size while the floor and vertical wall deposition rates both increase with particle size. The ceiling particle deposition rate increases with Reynolds number while the floor deposition rate decreases with it. The vertical deposition rate for the small particles (5–50?μm) increases with the flow Reynolds number obviously, while for the large particles (100–500?μm) that becomes insensitive. In addition, the flow Reynolds number is found to have an obvious effect on particle resuspension while the effect of particle size on particle resuspension decreases with Reynolds number. Eventually, a dynamic analysis was conducted for particles deposition and resuspension in turbulent duct flows.     

空分板翅式换热器热力性能动态仿真研究

本文主要通过建立空分板翅式换热器的物理模型和数学模型 ,运用数值计算手段对其热力性能控制方程进行了仿真 ,并对一具体的板翅式换热器给出了一些典型的仿真结果。在空分装置设计和生产实验过程中 ,板翅式换热器热力性能仿真是提高设计效率、缩短设计和试验时间、降低成本的有效手段     

Three-dimensional simulation of a plasma jet with transverse particle and carrier gas injection

In the present paper, three-dimensional modeling results are presented concerning the turbulent thermal plasma jet with transversely injected carrier gas and metal particles at atmospheric pressure. The standard K− model is employed for the numerical simulation of the turbulent plasma flow in coupling with the variable-property continuity, momentum and energy equations. For predicting the motion of the injected particles in the turbulent flow field, an improved particle stochastic-trajectory model is adopted in the calculation. The heating histories of the injected particles are also calculated in their moving processes. The modeling results show that including the effect of carrier gas on jet and particle behavior is very important. The plasma jet is deflected from its geometrical axis due to the transverse injection of carrier gas, and the particle trajectories are also appreciably changed by the carrier gas injection. The particles disperse around their average trajectories in the turbulent flow field.     

Inline method of droplet and particle size distribution analysis in dilute disperse systems

In the present article a measurement method of particle size distributions (PSD) in industrial installations which use a dispersed phase of low concentration (like spray dryers or spray scrubbers) is introduced. A new type of inline-measurement system has been developed and designed to work in spray drying conditions. A standard digital camera is used to record shadows of flowing particles inside the spray drying chamber. Collected images were analyzed by a newly developed software which recognizes particles only in the focus area and eliminates several types of artifacts. The constructed prototype of the PSD inline-analyzer was installed and used to monitor large laboratory scale spray dryer. All data collected by the designed system during the spray drying experiments were compared with data measured with an offline reference system to show accuracy of the new measurement technique.     

Numerical simulation of dispersed gas-liquid flows

This paper reviews the available information on numerical simulation of dispersed gas-liquid flows. Emphasis is on informing the reader about various aspects of constructing simulation models rather than giving an exhaustive literature review. The information is organised in a way so as to provide answers to the following questions: how to formulate model equations? how to select suitable algorithms and numerical techniques to solve these model equations? how to translate these into workable computer codes? how to use such codes for simulating flows in industrial equipment? Though greater emphasis is given to dispersed gas-liquid flows, the methodology can be applied to any multi-phase problem. Special features of multi-phase flow simulation over single-phase simulation are highlighted. A case of gas-liquid flow in a bubble column is presented as an illustration for the general methodology. The simulated mean flow field agrees reasonably with the experimental data. Properly validatedcfd codes thus can serve as a useful tool for design engineers of multi-phase systems. Some of the common pitfalls in using simulation codes for design are also discussed. This review is expected to give an overall idea about the present state-of-art of two-phase simulation in industrial equipment.     

CFD modelling of moisture evaporation in an industrial dispersed system

A methodology is presented here for the mathematical modelling of moisture evaporation in a dispersed system in an industrial tower. An empirical model using characteristic drying curves was applied to calculate moisture evaporation from a droplet and particle. A simple method was developed to calculate the agglomeration process of the dispersed phase in the drying towers, using transient functions between the initial and final particle size distributions, at the nozzle and in the final product. The developed model and simulation results were validated on the basis of industrial spray tower experiments. High instability of the airflow due to the geometry of the dryer and the construction of the air inlets was observed. The general methodology applied within this CFD model is universal, and can be applied to the scaling-up of installations for dewatering in dispersed systems in order to determine configurations of feeding systems and control the product quality and safety of the process.     

Effect of particle size distribution on hydrodynamics of pneumatic conveying system based on CPFD simulation

The effect of particle size distribution on the hydrodynamics of dilute-phase pneumatic conveying system was analyzed using computational particle fluid dynamics (CPFD) simulation. The influence of a simulation parameter, i.e., correction factor of drag coefficient (k), on the hydrodynamics of pneumatic conveying system was determined via CPFD simulation. When results of simulation were compared with experimental data of previous studies, the average error of pressure drop per length predicted by the CPFD approach with the correction factor was below 4.4%. Saltation velocity and the pressure drop per unit length declined as the drag force coefficient increased. Simulation results also revealed that the pressure drop per length and the saltation velocity were decreased when the fine powder fraction in the particle size distribution was increased, although the width of particle size distribution was widened, and the standard deviation was increased. Finally, the Relative Standard Deviation (RSD) of pressure drop per length was measured and compared with median diameter (d₅₀), Sauter mean diameter, geometric mean diameter, and arithmetic mean diameter. The RSD of the Sauter mean diameter was 5.8%, approximately twice less than the RSD value of d₅₀ commonly used in pneumatic conveying.     

Monte Carlo simulation for model-based fault diagnosis in dynamic systems

Fault diagnosis requires the accurate estimation of the dynamic state of the system in real time. This can be pursued starting from a model of the system dynamics and on measurements related to the state of the system. In real applications, the nonlinearity of the model and non-Gaussianity of the noise typically affecting the measurement challenge the classical approximate approaches, e.g. the extended-Kalman, Gaussian-sum and grid-based filters, which often turn out to be inaccurate and/or too computationally expensive for real-time applications. On the contrary, Monte Carlo estimation methods, also called particle filters, can be very effective. Based on sequential importance sampling and on a Bayesian formulation of the estimation problem, these methods recursively approximate the relevant probability distributions of the system state by random measures composed of particles (sampled values of the unknown state variables) and associated weights.The present paper aims at demonstrating the power of particle filtering for fault diagnosis. This is done by applying an estimation procedure called sampling importance resampling (SIR) to a case study of literature.     

电力推进船舶中央冷却系统动态特性建模及控制仿真

通过对电力推进船舶中央冷却系统的综合分析，建立电力推进船舶中冷系统的稳态换热模型和动态换热模型。对中冷系统的冷却水出口温度的控制方案，进行了传统的PID控制和基于功率的前馈、反馈控制研究，并在此基础上对控制方案计算机仿真，比较两种控制方案的区别。     

Self-consistent particle simulation of radio frequency SF6 discharges: Role of ion recombination

The role of positive−negative ion recombination in simulating radio frequency (rf) SF₆ discharges is investigated based on Nanbu and Denpoh's theory [J. Phys. Soc. Jpn. 1998; 67: 1288]. The model of ion recombination is combined with the Particle-in-Cell/Monte Carlo (PIC/MC) simulation. Results show that the ion recombination is a dominant factor in simulating a steady rf discharge of SF₆. A criterion determining whether a periodic steady state is reached is proposed based on the ion recombination rate. It is found that the ion recombination rapidly increases as the gas pressure increases. From 25 to 50 m Torr, the loss of positive ions caused by ion recombination becomes comparable with that lost from collisions with electrodes.     

Numerical simulation of the guided Lamb wave propagation in particle reinforced composites

This paper deals with the investigation of the Lamb wave propagation in particle reinforced composites excited by piezoelectric patch actuators. A three-dimensional finite element method (FEM) modeling approach is set up to perform parameter studies in order to better understand how the Lamb wave propagation in particle reinforced composite plates is affected by change of central frequency of excitation signal, volume fraction of particles, size of particles and stiffness to density ratio of particles. Furthermore, the influence of different arrangements is investigated. Finally, the results of simplified models using material data obtained from numerical homogenization are compared to the results of models with heterogeneous build-up. The results show that the Lamb wave propagation properties are mainly affected by the volume fraction and ratio of stiffness to density of particles, whereas the particle size does not affect the Lamb wave propagation in the considered range. As the contribution of the stiffer material increases, the group velocity and the wave length also increase while the energy transmission reduces. Simplified models based on homogenization technique enabled a tremendous drop in computational costs and show reasonable agreement in terms of group velocity and wave length.     

Dynamic behaviour of dry and water-saturated sand under planar shock conditions

Plane shock-wave experiments were performed on dry and partially water-saturated sand, using three water contents, in order to validate predictive models of material behaviour at stress levels between 1 and 10 GPa. Gas and powder guns were used to load the sample under uni-axial strain conditions at low and high stress levels, respectively. Wave motions were detected by piezoelectric pins in the samples and a VISAR (Velocity Interferometer System for Any Reflector) recorded the free-surface velocity on the back target. This study presents both experimental and simulated results. Experimental data are used to determine shock Hugoniot states. Significant differences are observed in the dynamic response of the materials under various water-saturated conditions, and are reproduced with good agreement by numerical simulations using the ARMORS (A Rheological MOdel of Rocks Saturated) model.     

Preparation and characterization of highly dispersed nanocrystalline rutile powders

A highly dispersed nanocrystalline rutile powder (NRP) having ultrafine and narrow-distributed diameters was successfully prepared via an improved conventional liquid one-step method by introducing dioctyl sulfosuccinate sodium salt (AOT) surfactant in the peptization reaction anaphase. The average grain size of the product is 10-20 nm. The effects of the titanium tetrabutoxide (Ti(OBu)₄) to ethanol volumetric ratio, the precipitation temperature, the peptization temperature and the presence of dioctyl sulfosuccinate sodium salt (AOT) on the preparation of NRP were studied. A new procedure for synthesizing NRP is proposed, and the optimum conditions for preparing NRP were obtained. The TEM, XRD, SEM and FT-IR results were also investigated in this work.     

Simulation of microstructural evolution in materials with misfitting precipitates

To predict the behavior of directional coarsening and the temporal evolution of the shape of coherent precipitates in two-phase materials, a dislocation-free model is proposed, based on a combination of statistical mechanics and linear elasticity. This model takes elastic anisotropy and isotropic interfacial energy into account. Based on an example of isolated precipitates under plane strain condition, the influence of particle size, inhomogeneity, direction and sign of external loads on the equilibrium shape will be discussed in terms of a generalized thermodynamic force acting on the interface. To simulate the morphological diffusion process of typical microstructures with several random distributed misfitting inclusions, a computational technique in form of a finite element Monte Carlo simulation is presented. Within this numerical technique, no restrictions on the particle shape or the elastic anisotropy of both phases are made.     

Influence of particle size of the dispersoid on compressibility and sinterability of TiO2 dispersed Al 7075 alloy composites prepared by mechanical milling

Particulate TiO₂ (with varying particle size produced by mechanical milling) dispersed AA7075 composites are synthesised by short duration milling (10 min) followed by room temperature unidirectional compaction (with varying pressure) and sintering. Apparent and relative density of the alloy powder and composites are measured. The effect of reinforcement particle size on the compressibility behaviour of the composites is demonstrated. Mechanically milled (for 25 h) alloy powder shows lower relative density than coarse alloy powder. In addition, compressibility of the alloy composites decreases with decreasing particle size of the reinforcement. In contrast, the sinterability of the composites increases with decreasing dispersoid’s size due to easy filling up of finer pores and particle induced precipitation.     

A Monte Carlo simulation approach to the availability assessment of multi-state systems with operational dependencies

This paper deals with multi-state systems (MSS), whose performance can settle on different levels, e.g. 100%, 80%, 50% of the nominal capacity, depending on the operative conditions of the constitutive multi-state elements. Examples are manufacturing, production, power generation and gas and oil transportation systems. Often in practice, MSS are such that operational dependencies exist between the system state and the state of its components. For example, in a production line of nodal series structure, with no buffers between the nodes, if one of the nodes throughput changes (e.g. switches from 100% to 50% due to a deterministic or stochastic transition of one of its components), the other nodes must be reconfigured (i.e. their components must deterministically change their states) so as to provide the same throughput.In this paper, we present a Monte Carlo simulation technique which allows modelling the complex dynamics of multi-state components subject to operational dependencies with the system overall state. A correlation method is tailored to model the automatic change of state of the relevant components following a change in one of the system nodes. The proposed technique is verified on a simple case study of literature.     

A rapid route for the synthesis of submicron Se and Te rod-like crystals

Submicron Se and Te rod-like crystals were successfully synthesized via a rapid polyol process by refluxing Na₂EO₃ (ESe, Te) and (NH₃)₂S₂O₃ in ethylene glycol at 180 °C for 0.5 h under vigorous stirring. The obtained samples were characterized by X-ray diffraction (XRD), scanning electron microscope (SEM), transmission electron microscopy (TEM), and XPS procedures. Studies found that reaction time and temperature have great influences on the morphology of the final products.     

Dynamic simulation of particle packing influenced by size, aspect ratio and surface energy

The multi-sphere method and JKR model are used in the discrete element method simulation to investigate the effect of the particle size, aspect ratios, and cohesiveness on the packing structure, characterized by porosity, radial distribution function (RDF), coordination number and contact geometry. In the absence of cohesive force, the porosity is nearly invariable with fixed aspect ratio, regardless of the size of the particles. In contrast, as surface energy increases, the porosity increases with decreasing particle size and increasing aspect ratio. The RDF results show that the number of peaks for different aspect ratios changes and show trends similar to the relaxation algorithm, expected for the finer particles. In the case of finer, cohesive particles, the most novel outcome of contact analysis is the existence of single contact, attributed to the formation of a cage structure, which has not been previously reported. The peak position and the width of the contact distributions are affected by higher surface energy because fewer contacts are required to achieve the mechanical equilibrium. Another interesting observation is that higher porosity does not always imply fewer contacts for particles with non-zero aspect ratios and high surface energies. The analysis of the distribution of the contact vector angles is found to better explain increased porosity in spite of higher coordination numbers. The results presented shed light on the packing density and structure, revealing features not easily discerned via experiments, and confirming the important role of the cohesion and aspect ratio in packing.     

Dynamic simulation of one-stage Oxford split-Stirling cryocooler and comparison with experiment

In this paper, a mathematical model of a one-stage Oxford split-Stirling cryocooler has been made to simulate its dynamic performance. A set of governing equations has been solved by a numerical method. The present paper emphasizes the validation of this model against experimental results for the Oxford split-Stirling cryocooler. The comparison shows acceptable agreement between the simulation and experimental results. The simulated results on compression work have relative discrepancies of 6.8% (cold-tip temperature 80 K) and 6.5% (cold-tip temperature 54 K) with their experimental results respectively. The greatest deviations for the minimum pressure and net cooling power are about 6.0% and 8.0% respectively.