Particle-in-cell simulations of EUV-induced hydrogen plasma in the vicinity of a reflective mirror

Particle-In-Cell(PIC)simulations were performed in this work to study the dynamics of the EUVinduced hydrogen plasma.The Monte-Carlo Collision(MCC)model was employed to deal with the collisions between charged particles and background gas molecules.The dynamic evolution of the plasma sheath,as well as the flux and energy distribution of ions impacting on the mirror surface,was discussed.It was found that the emission of secondary electrons under the EUV irradiation on the ruthenium mirror coating creates a positively charged wall and then prevents the ions from impacting on the mirror and therefore changes the flux and energy distribution of ions reaching the mirror.Furthermore,gas pressure has a notable effect on the plasma sheath and the characteristics of the ions impinging on the mirrors.With greater gas pressure,the sheath potential decreases more rapidly.The flux of ions received by the mirror grows approximately linearly and at the same time the energy corresponding to the peak flux decreases slightly.However,the EUV source intensity barely changes the sheath potential,and its influence on the ion impact is mainly limited to the approximate linear increase in ion flux.

Exploration of microscopic physical processes of Z-pinch by a modified electrostatic direct implicit particle-in-cell algorithm

For investigating efficiently the stagnation kinetic-process of Z-pinch,we develop a novel modified electrostatic implicit particle-in-cell algorithm in radial one-dimension for Z-pinch simulation in which a small-angle cumulative binary collision algorithm is used.In our algorithm,the electric field in z-direction is solved by a parallel electrode-plate model,the azimuthal magnetic field is obtained by Ampere’s law,and the term for charged particle gyromotion is approximated by the cross product of the averaged velocity and magnetic field.In simulation results of 2 MA deuterium plasma shell Zpinch,the mass-center implosion trajectory agrees generally with that obtained by one-dimensional MHD simulation,and the plasma current also closely aligns with the external current.The phase space diagrams and radial-velocity probability distributions of ions and electrons are obtained.The main kinetic characteristic of electron motion is thermal equilibrium and oscillation,which should be oscillated around the ions,while that of ion motion is implosion inwards.In the region of stagnation radius,the radial-velocity probability distribution of ions transits from the non-equilibrium to equilibrium state with the current increasing,while of electrons is basically the equilibrium state.When the initial ion density and current peak are not high enough,the ions may not reach their thermal equilibrium state through collisions even in its stagnation phase.

Implementation of a particle-in-cell method for the energy solver in 3D spherical geodynamic modeling

The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.

MULTI-PHASE ACTIVE CONTOUR MODEL FOR IMAGE SEGMENTATION BASED ON LEVEL SETS

A new multi-phase active contour model is proposed for the image segmentation. It is a generalization of the C-V model with the following characteristics： （1） A key technique, called the technique of painting background （TPBG）, is developed to remove the information of the background, which blocks the detection of weak boundaries in the object; （2） The two-phase level set is applied multiple times for getting the multi-phase segmentation model （n-1 times for the n-phase model, n〉1）; （3） A scaling-based method is introduced to improve the basic model. Experimental results show that the proposed model is effective for detecting weak boundaries.

Multi-phase signal setting and capacity of signalized intersection

A capacity model of multi-phase signalized intersections is derived by a stopping-line method. It is simplified with two normal situations： one situation involves one straight lane and one left-turn lane; the other situation involves two straight lanes and one left-turn lane. The results show that the capacity is mainly relative to signal cycle length, phase length, intersection layout and following time. With regard to the vehicles arrival rates, the optimal model is derived based on each phase＇s remaining time balance, and it is solved by Lagrange multipliers. Therefore, the calculation models of the optimal signal cycle length and phase lengths are derived and simplified. Compared to the existing models, the proposed model is more convenient and practical. Finally, a practical intersection is chosen and its signal cycles and phase lengths are calculated by the proposed model.

5千瓦霍尔推力器磁场优化及Particle-in-Cell性能仿真

为了提升兰州空间技术物理研究所研制的一种5 k W霍尔推力器LHT—140的性能,采用ANSOFT软件进行了磁场优化设计,将磁场径向分量的轴向梯度提高了47%,相同励磁激励下放电通道中的磁场强度提高了38%。建立了一个R-Z平面内的二维particle-in-cell(PIC)等离子体模型,对磁场优化后推力器的性能进行了仿真分析,预估其在300~800 V放电电压、10~15 mg/s流率范围内,推力提升了9.6%~22.1%,效率提升了8.7%~19.3%。性能验证试验表明磁场优化后在相同放电电压与流率下,性能提升的测试值高于仿真值。

Numerical simulation of complex multi-phase fluid of casting process and its applications

The fluid of casting process is a typical kind of multi-phase flow. Actually, many casting phenomena have close relationship with the multi-phase flow, such as molten metal filling process, air entrapment, slag movement, venting process of die casting, gas escaping of lost foam casting and so on. Obviously, in order to analyze these phenomena accurately, numerical simulation of the multi-phase fluid is necessary. Unfortunately, so far, most of the commercial casting simulation systems do not have the ability of multi-phase flow modeling due to the difficulty in the multi-phase flow calculation. In the paper, Finite Different Method (FDM) technique was adopt to solve the multi-phase fluid model. And a simple object of the muiti-phase fluid was analyzed to obtain the fluid rates of the liquid phase and the entrapped air phase.

A New Method for Hydrocarbon Detection Based on Multi-phase Theory

The hydrocarbon detection techniques used currently are generally based on the theory of single-phase medium, but hydrocarbon reservoir mostly is multi-phase medium, therefore, multisolutions and uncertainties are existed in the result of hydrocarbon detection. This paper presents a fast way to detect hydrocarbon in accordance with BOIT theory and laboratory data. The technique called DHAF technique has been applied to several survey area and obtained good result where the coincidence rate for hydrocarbon detection is higher than other similar techniques. The method shows a good prospect of the application in hydrocarbon detecting at exploration stage and in reservoir monitoring at production stage.

A proposed NMR solution for multi-phase flow fluid detection

In the petroleum industry,detection of multi-phase fluid flow is very important in both surface and down-hole measurements.Accurate measurement of high rate of water or gas multi-phase flow has always been an academic and industrial focus.NMR is an efficient and accurate technique for the detection of fluids;it is widely used in the determination of fluid compositions and properties.This paper is aimed to quantitatively detect multi-phase flow in oil and gas wells and pipelines and to propose an innovative method for online nuclear magnetic resonance(NMR)detection.The online NMR data acquisition,processing and interpretation methods are proposed to fill the blank of traditional methods.A full-bore straight tube design without pressure drop,a Halbach magnet structure design with zero magnetic leakage outside the probe,a separate antenna structure design without flowing effects on NMR measurement and automatic control technology will achieve unattended operation.Through the innovation of this work,the application of NMR for the real-time and quantitative detection of multi-phase flow in oil and gas wells and pipelines can be implemented.

A new analytical model for thermal stresses in multi-phase materials and lifetime prediction methods

Based on the fundamental equations of the mechanics of solid continuum, the paper employs an analytical model for determination of elastic thermal stresses in isotropic continuum represented by periodically distributed spherical particles with different distributions in an infinite matrix, imaginarily divided into identical cells with dimensions equal to inter-particle distances, containing a central spherical particle with or without a spherical envelope on the particle surface. Consequently, the multi-particle-（envelope）- matrix system, as a model system regarding the analytical modelling, is applicable to four types of multi-phase materials. As functions of the particle volume fraction v, the inter-particle distances dl, d2, d3 along three mutually per- pendicular axes, and the particle and envelope radii, R1 and R2, respectively, the thermal stresses within the cell, are originated during a cooling process as a consequence of the difference in thermal expansion coefficients of phases rep- resented by the matrix, envelope and particle. Analytical-（experimental）-computational lifetime prediction methods for multi-phase materials are proposed, which can be used in engineering with appropriate values of parameters of real multi-phase materials.

Phase Transition Analysis Based Quality Prediction for Multi-phase Batch Processes

Batch processes are usually involved with multiple phases in the time domain and many researches on process monitoring as well as quality prediction have been done using phase information. However, few of them consider phase transitions, though they exit widely in batch processes and have non-ignorable impacts on product qualities. In the present work, a phase-based partial least squares (PLS) method utilizing transition information is proposed to give both online and offline quality predictions. First, batch processes are divided into several phases using regression parameters other than prior process knowledge. Then both steady phases and transitions which have great influences on qualities are identified as critical-to-quality phases using statistical methods. Finally, based on the analysis of different characteristics of transitions and steady phases, an integrated algorithm is developed for quality prediction. The application to an injection molding process shows the effectiveness of the proposed algorithm in comparison with the traditional MPLS method and the phase-based PLS method.

Optimization study of a PEM fuel cell performance using 3D multi-phase computational fluid dynamics model

An optimization study using a comprehensive 3D, multi-phase, non-isothermal model of a PEM （proton exchange membrane） fuel cell that incorporates significant physical processes and key parameters affecting fuel cell performance is presented and discussed in detail. The model accounts for both gas and liquid phase in the same computational domain, and thus allows for the implementation of phase change inside the gas diffusion layers. The model includes the transport of gaseous species, liquid water, protons, energy, and water dissolved in the ion-conducting polymer. Water is assumed to be exchanged among three phases： liquid, vapottr, and dissolved, with equilibrium among these phases being assumed. This model also takes into account convection and diffusion of different species in the channels as well as in the porous gas diffusion layer, heat transfer in the solids as well as in the gases, and electrochemical reactions. The results showed that the present multi-phase model is capable of identifying important parameters for the wetting behaviour of the gas diffusion layers and can be used to identify conditions that might lead to the onset of pore plugging, which has a detrimental effect on the fuel cell performance. This model is used to study the effects of several operating, design, and material parameters on fuel cell performance. Detailed analyses of the fuel cell performance under various operating conditions have been conducted and examined.

Multi-phase field simulation of grain growth in multiple phase transformations of a binary alloy

This work establishes a temperature-controlled sequence function, and a new multi-phase-field model, for liquid- solid-solid multi-phase transformation by coupling the liquid-solid phase transformation model with the solid-solid phase transformation model. Taking an Fe-C alloy as an example, the continuous evolution of a multi-phase transformation is simulated by using this new model. In addition, the growth of grains affected by the grain orientation of the parent phase （generated in liquid-solid phase transformation） in the solid-solid phase transformation is studied. The results show that the morphology of ferrite grains which nucleate at the boundaries of the austenite grains is influenced by the orientation of the parent austenite grains. The growth rate of ferrite grains which nucleate at small-angle austenite grain boundaries is faster than those that nucleate at large-angle austenite grain boundaries. The difference of the growth rate of ferrites grains in different parent phase that nucleate at large-angle austenite grain boundaries, on both sides of the boundaries, is greater than that of ferrites nucleating at small-angle austenite grain boundaries.

On the energy conservation electrostatic particle-in-cell/Monte Carlo simulation: Benchmark and application to the radio frequency discharges

We benchmark and analyze the error of energy conservation （EC） scheme in particle-in-cell/Monte Carlo （PIC/MC） algorithms by simulating the radio frequency discharge. The plasma heating behaviors and electron distributing functions obtained by one-dimensional （1D） simulation are analyzed. Both explicit and implicit algorithms are checked. The results showed that the EC scheme can eliminated the self-heating with wide grid spacing in both cases with a small reduction of the accuracies. In typical parameters, the EC implicit scheme has higher precision than EC explicit scheme. Some ＂numerical cooling＂ behaviors are observed and analyzed. Some other errors are also analyzed. The analysis showed that the EC implicit scheme can be used to qualitative estimation of some discharge problems with much less computational resource cost without much loss of accuracies.

Mechanical properties of irradiated multi-phase polycrystalline BCC materials

Structure materials under severe irradiations in nuclear environments are known to degrade because of irradiation hardening and loss of ductility,resulting from irradiation-induced defects such as vacancies,interstitials and dislocation loops,etc.In this paper,we develop an elastic-viscoplastic model for irradiated multi-phase polycrystalline BCC materials in which the mechanical behaviors of individual grains and polycrystalline aggregates are both explored.At the microscopic grain scale,we use the internal variable model and propose a new tensorial damage descriptor to represent the geometry character of the defect loop,which facilitates the analysis of the defect loop evolutions and dislocation-defect interactions.At the macroscopic polycrystal scale,the self-consistent scheme is extended to consider the multiphase problem and used to bridge the individual grain behavior to polycrystal properties.Based on the proposed model,we found that the work-hardening coefficient decreases with the increase of irradiation-induced defect loops,and the orientation/loading dependence of mechanical properties is mainly attributed to the different Schmid factors.At the polycrystalline scale,numerical results for pure Fe match well with the irradiation experiment data.The model is further extended to predict the hardening effect of dispersoids in oxide-dispersed strengthened steels by the considering the Orowan bowing.The influences of grain size and irradiation are found to compete to dominate the strengthening behaviors of materials.

THE RATE-INDEPENDENT CONSTITUTIVE MODELING FOR POROUS AND MULTI-PHASE NANOCRYSTALLINE MATERIALS

To determine the time-independent constitutive modeling for porous and multi- phase nanocrystalline materials and understand the effects of grain size and porosity on their mechanical behavior, each phase was treated as a mixture of grain interior and grain bound- ary, and pores were taken as a single phase, then Budiansky＇s self-consistent method was used to calculate the Young＇s modulus of porous, possible multi-phase, nanocrystalline materials, the prediction being in good agreement with the results in the literature. Further, the established method is extended to simulate the constitutive relations of porous and possible multi-phase nanocrystalline materials with small plastic deformation in conjunction with the secant-moduli approach and iso-strain assumption. Comparisons between the experimental grain size and porosity dependent mechanical data and the corresponding predictions using the established model show that it appears to be capable of describing the time-independent mechanical behaviors for porous and multi-phase nanocrystalline materials in a small plastic strain range. Further discussion on the modification factor, the advantages and limitations of the model developed were present.

Multi-phase computer modeling and laboratory study of dust capture by an inertial Vortecone scrubber

Dust generated in mining and tunneling activities is hazardous to health of persons and safety of operations. These projects employ pick-milling machines to extract minerals and rock by mechanical breakage.The machines are equipped with flooded-bed scrubbers that encase dust particles within fine water films as particles encounter a flooded wire-mesh screen. A major disadvantage is that the screen gets clogged when particles become trapped within the wire mesh, reducing airflow through the scrubber and increasing ambient dust concentrations. Thus, the system requires frequent maintenance or replacement. The application of a Vortecone scrubber as an improved alternative to conventional fibrous type scrubbers is investigated. A Vortecone forces dust-laden air and water to follow a complex, rapidly swirling motion.The momentum drives dust particles towards the periphery where they are captured by the water film.The operating characteristics of a reduced-scale physical model of a Vortecone, with its primary axis mounted in the horizontal orientation, was analyzed numerically and experimentally. Computational fluid dynamics(CFD) models depicting the spraying action and multi-phase air/water flows using the volume of fraction(VOF) approach, are presented. Experimental results, utilizing an optical particle counting technique to establish the dust-cleaning capabilities of the model, are also described.

Electromagnetic Particle-in-Cell Simulations of Electron Holes Formed During the Electron Two-Stream Instability

Previous electrostatic particle-in-cell （PIC） simulations have pointed out that elec- tron phase-space holes （electron holes） can be formed during the nonlinear evolution of the electron two-stream instability. The parallel cuts of the parallel and perpendicular electric field have bipolar and unipolar structures in these electron holes, respectively. In this study, two-dimensional （2D） electromagnetic PIC simulations are performed in the x - y plane to investigate the evolution of the electron two-stream instability, with the emphasis on the magnetic structures associated with these electron holes in different plasma conditions. In the simulations, the background magnetic field （Bo = Boer） is along the x direction. In weakly magnetized plasma （Ωe 〈ωpe, where Ωe and ωpe are the electron gyrofrequency and electron plasma frequency, respectively）, several 2D electron holes are formed. In these 2D electron holes, the parallel cut of the fluctuating magnetic field δBx and δBz has unipolar structures, while the fluctuating magnetic field δBy has bipolar structures. In strongly magnetized plasma （Ωe 〉 ωpe）, several quasi-lD electron holes are formed. The electrostatic whistler waves with streaked structures of Ey are excited. The fluctuating mag- netic field δBx and δBz also have streaked structures. The fluctuating magnetic field δBx and δBy are produced by the current in the z direction due to the electric field drift of the trapped elec- trons, while the fluctuating magnetic field δBz can be explained by the Lorentz transformation of a moving quasielectrostatic structure. The influences of the initial temperature anisotropy on the magnetic structures of the electron holes are also analyzed. The electromagnetic whistler waves are found to be excited in weakly magnetized plasma. However, they do not have any significant effects on the electrostatic structures of the electron holes.

Implicit electrostatic particle-in-cell/Monte Carlo simulation for the magnetized plasma:Algorithms and application in gas-inductive breakdown

An implicit electrostatic particle-in-cell/Monte Carlo （PIC/MC） algorithm is developed for the magnetized discharging device simulation. The inductive driving force can be considered. The direct implicit PIC algorithm （DIPIC） and energy conservation scheme are applied together and the grid heating can be eliminated in most cases. A tensor-susceptibility Poisson equation is constructed. Its discrete form is made up by a hybrid scheme in one-dimensional （1D） and two- dimensional （2D） cylindrical systems. A semi-coarsening multigrid method is used to solve the discrete system. The algorithm is applied to simulate the cylindrical magnetized target fusion （MTF） pre-ionization process and get qualitatively correct results. The potential application of the algorithm is discussed briefly.

Numerical Simulation of Shaped Charge Jet Using Multi-Phase SPH Method

Since the jets and detonation gaseous products are separated by sharp interfaces, the traditional smoothed particle hydrodynamics (SPH) method is difficult to avoid the computational instability at interfaces. The multi-phase SPH (MSPH) method was applied to improving the stabil-ity, which smoothes the particle density and makes pressure continuous at interfaces. Numericalexamples of jet forming process were used to test capability of the MSPH method. The results show that the method remains algorithm stability for large density gradient between the jets and gaseous products and has potential application to both the explosion and the jet problems. The effect of initiation ways of the shaped charge was discussed as well.