Monte Carlo simulation of phase separation of A/B/A-B ternary mixtures

Monte Carlo simulations were used to model A/B/A-B ternary mixtures with different AB diblock copolymer volume fractions for which both the dispersed and continuous phase volume fractions were kept constant. For concentrations of the diblock copolymer below a critical value, the domain size increment of the dispersed phase decreases linearly with the copolymer concentration. This is in agreement with the predictions of Noolandi and Hong. The dependence of the domain size as a function of the copolymer volume fraction can also be fitted by the equation of Tang and Huang. Our simulations indicate, for the first time, that the micelles form before saturation of the interface occurs. This means that the formation of the micelles is not a result of the saturation of the interface.

Size distributions of quantum islands on stepped substrates

The size distributions of self-assembled quantum islands on stepped substrates are studied using kinetic Monte Carlo simulations. It is found that the energy barrier E-SW between the step and the terrace region is the key factor in affecting the size distribution of islands. With small E-SW (<= 0.1 eV), lines of uniform islands can be obtained at relative low surface coverage. As the surface coverage is increased, wirelike islands can be obtained. Scaling behavior is obeyed for the size distributions of the wirelike islands. When the size distributions are separated into their width and length components, however, scaling is only observed in the length distribution of the wirelike islands. With larger E-SW, the size distribution of islands shows a clear bimodal size distribution and anomalous growth temperature dependent island size evolutions are observed. The simulation results reproduce qualitatively the phenomena observed in the cases of InAs islands grown on stepped GaAs substrates. (c) 2009 American Institute of Physics. [doi:10.1063/1.3248367]

Quantum-dot growth simulation on periodic stress of substrate

InAs quantum dots (QDs) are grown on the cleaved edge of an InxGa1-xAs/GaAs supperlattice experimentally and a good linear alignment of these QDs on the surface of an InxGa1-xAs layer has been realized. The modulation effects of periodic strain on the substrate are investigated theoretically using a kinetic Monte Carlo method. Our results show that a good alignment of QDs can be achieved when the strain energy reaches 2% of the atomic binding energy. The simulation results are in excellent qualitative agreement with our experiments. (C) 2005 American Institute of Physics.

Formation trends of ordered self-assembled nanoislands on stepped substrates

The growth of ordered self-assembled nanoislands on stepped substrates is studied systematically by kinetic Monte Carlo simulations. As the terrace width is small, the formation of nanoislands is confined in the steps and nanoislands ordered in lines or nanowires can be obtained. The Schwoebel barrier at the step edges has a great influence on the evolution of both the size and space distributions of the islands. When the terrace width is relatively large, self-ordering of nanoislands in the center regions of the terraces happens. An unexpected trend of the nanoisland self-ordering is found as the deposition thickness is larger than 0.2 ML, which can be related to the attractive migrations between nearby islands.

Kinetic Study Of The Rayleigh-Benard Flows

The present paper employs the direct simulation Monte Carlo (DSMC) method to study the Rayleigh-Benard flows, where the temperature ratio of the upper to lower plate is fixed to 0.1. For a Knudsen number (Kn) of 0.01, as the Rayleigh number (Ra) increases, the flow changes from the thermal conductive state to the convective state at about Ra=1700, and the calculated relation of heat flux through the lower plate versus Ra is in good agreement with classical experimental and theoretical results. For Kn=0.05, the thermal conductive state remains stable, and the increase of Ra cannot trigger thermal instability.

Kinetic studies of gas flows

Our recent studies on kinetic behaviors of gas flows are reviewed in this paper. These flows have a wide range of background, but share a common feature that the flow Knudsen number is larger than 0.01. Thus kinetic approaches such as the direct simulation Monte Carlo method are required for their description. In the past few years, we studied several micro/nano-scale flows by developing novel particle simulation approach, and investigated the flows in low-pressure chambers and at high altitude. In addition, the microscopic behaviors of a couple of classical flow problems were analyzed, which shows the potential for kinetic approaches to reveal the microscopic mechanism of gas flows.

Kinetic Studies of Gas Flows

Our recent studies on kinetic behaviors of gas flows are reviewed in this paper. These flows have a wide range of background, but share a common feature that the flow Knudsen number is larger than 0.01. Thus kinetic approaches such as the direct simulation Monte Carlo method are required for their description. In the past few years, we studied several micro/nano-scale flows by developing novel particle simulation approach, and investigated the flows in low-pressure chambers and at high altitude. In addition, the microscopic behaviors of a couple of classical flow problems were analyzed, which shows the potential for kinetic approaches to reveal the microscopic mechanism of gas flows.

Vaporization of heated materials into discharge plasmas

The vaporization of condensed materials in contact with high-current discharge plasmas is considered. A kinetic numerical method named direct simulation Monte Carlo (DSMC) and analytical kinetic approaches based on the bimodal distribution function approximation are employed. The solution of the kinetic layer problem depends upon the velocity at the outer boundary of the kinetic layer which varies from very small, corresponding to the high-density plasma near the evaporated surface, up to the sound speed, corresponding to evaporation into vacuum. The heavy particles density and temperature at the kinetic and hydrodynamic layer interface were obtained by the analytical method while DSMC calculation makes it possible to obtain the evolution of the particle distribution function within the kinetic layer and the layer thickness.

Rarefied gas dynamics: fundamentals, simulations and micro flows

Strain relaxation and band-gap tunability in ternary InxGa1-xN nanowires

The alloy formation enthalpy and band structure of InGaN nanowires were studied by a combined approach of the valence-force field model, Monte Carlo simulation, and density-functional theory (DFT). For both random and ground-state structures of the coherent InGaN alloy, the nanowire configuration was found to be more favorable for the strain relaxation than the bulk alloy. We proposed an analytical formula for computing the band gap of any InGaN nanowires based on the results from the screened exchange hybrid DFT calculations, which in turn reveals a better band-gap tunability in ternary InGaN nanowires than the bulk alloy.

Analog-digital and digital-analog converters using single-electron and MOS transistors

This paper proposes two kinds of novel single-electron analog-digital conversion (ADC) and digital-analog conversion (DAC) circuits that consist of single-electron transistors (SETs) and metal-oxide-semiconductor (MOS) transistors. The SET/MOS hybrid ADC and DAC circuits possess the merits of the SET circuit and the MOS circuit. We obtain the SPICE macro-modeling code of the SET transistor by studying and fitting the characteristics of the SET with SPICE simulation and Monte Carlo simulation methods. The SPICE macro-modeling code is used for the simulation of the SET/MOS hybrid ADC and DAC circuits. We simulate the performances of the SET/MOS hybrid 3-b ADC and 2-b DAC circuits by using the H-SPICE simulator. The simulation results demonstrate that the hybrid circuits can perform analog-digital and digital-analog data conversion well at room temperature. The hybrid ADC and DAC circuits have advantages as-follows: 1) compared with conventional circuits, the architectures of the circuits are simpler; 2) compared with single electron transistor circuits, the circuits have much larger load capability; 3) the power dissipation of the circuits are lower than uW; 4) the data conversion rate of the circuits can exceed 100 MHz; and 5) the resolution of the ADC and DAC circuits can be increased by the pipeline architectures.

Stress transfer and damage evolution simulations of fiber-reinforced polymer-matrix composites

The stress transfer from broken fibers to unbroken fibers in fiber-reinforced thermosetting polymer-matrix composites and thermoplastic polymer-matrix composites was studied using a detailed finite element model. In order to check the validity of this approach, an epoxy-matrix monolayer composite was used as thermosetting polymer-matrix composite and a polypropylene (PP)-matrix monolayer composite was used as thermoplastic polymer-matrix composite, respectively. It is found that the stress concentrations near the broken fiber element cause damage to the neighboring epoxy matrix prior to the breakage of other fibers, whereas in the case of PP-matrix composites the fibers nearest to the broken fiber break prior to the PP matrix damage, because the PP matrix around the broken fiber element yields. In order to simulate composite damage evolution, a Monte Carlo technique based on a finite element method has been developed in the paper. The finite element code coupled with statistical model of fiber strength specifically written for this problem was used to determine the stress redistribution. Five hundred samples of numerical simulation were carried out to obtain statistical deformation and failure process of composites with fixed fiber volume fraction.

海底金属矿床充水条件及矿坑突水的断层效应研究

The Xinli mine area of Sanshandao mine is adjacent to the Bohai Sea and its main exploitable ore deposit occurs in the undersea rock mass. The mine is the biggest undersea gold mine of China after production. The mine area faces a latent danger of water bursting, even sudden seawater inrush. There is no mature experience in undersea mining in China so far. The vein ore deposit is located in the lower wall of a fault; its possible groundwater sources mainly include bittern, Quaternary pore water and modern seawater. To ensure the safety of undersea mining, to survey the flooding conditions of the ore deposit using proper measures and study the potential seawater inrush pattern are the key technical problems. With the Xinli mine area as a case study, the engineering geological conditions of the Xinli mine area are surveyed in situ, the regional structural pattern and rock mass framework characteristics are found out, the distribution of the structural planes are modeled by a Monte Carlo method and the connectivity coefficients of rock mass structural planes are calculated. The regional hydro-geological conditions are analyzed and the in-situ hydro-geological investigation and sampling are performed in detail, the hydrochemistry and isotopes testing and groundwater dynamic monitoring are conducted, the recharge, runoff, discharge conditions are specified and the sources of flooding are distinguished. Some indices are selected from the testing results to calculate the proportion of each source in some water discharge points and in the whole water discharge of the Xinli mine area. The temporal and spatial variations of each water source of the whole ore deposit flooding are analyzed. According to the special project conditions in the Xinli mine area, the permeability coefficient tensors of the rock mass in Xinli mine area are calculated based on a fracture geometry measurement method, in terms of the connectivity and a few hydraulic testing results, a modified synthetic permeability coefficient are calculated. The hydro-geological conceptual and mathematical model are established，the water yield of mine is predicted using Visual Modflow code. The spreading law of surrounding rock mass deformation and secondary stress are studied by numerical analysis; the intrinsic mechanism of the faults slip caused by the excavation of ore deposit is analyzed. The results show that the development of surrounding rock mass deformation and secondary stress of vein ore deposit in the lower wall of a fault, is different from that in a thick-big ore deposit. The secondary stress caused by the excavation of vein ore deposit in the lower wall of a fault, is mainly distributed in the upper wall of the fault, one surface subsidence center will occur. The influences of fault on the rock mass movement, secondary stress and hydro-geological structures are analyzed; the secondary stress is blocked by the fault and the tensile stress concentration occurs in the rock mass near the fault, the original water blocking structure is destructed and the permeable structure is reconstructed, the primary structural planes begin to expand and newborn fissures occur, so the permeability of the original permeable structure is greatly enhanced, so the water bursting will probably occur. Based on this knowledge, the possible water inrush pattern and position of the Xinli mine area are predicted. Some computer programs are developed using object-oriented design method under the development platform Visual Studio.Net. These programs include a Monte Carlo simulation procedure, a joint diagrammatizing procedure, a structural planes connectivity coefficient calculating procedure, a permeability tensor calculating procedure, a water chemical formula edit and water source fixture conditions calculating procedure. A new computer mapping algorithm of joint iso-density diagram is raised. Based on the powerful spatial data management and icon functions of Geographic Information System, the pit water discharge dynamic monitoring data management information systems are established with ArcView.