Molecular dynamics simulations of the crystal–melt interface mobility in HCP Mg and BCC Fe

The temperature profile around the moving solid–liquid interface during non-equilibrium molecular dynamics (MD) simulations of crystallization and melting is examined for HCP Mg and BCC Fe. An evident spike (valley) is found around the solid–liquid interface during solidification (melting). Considering the effect of a non-uniform temperature distribution, it is found that, if the actual interface temperature is adopted to compute the interface mobility, rather than the thermostat temperature (or the mean temperature of the whole system), the kinetic coefficient is approximately a factor of two larger than previous estimates. Although the magnitude of the kinetic coefficient is larger than the previous estimates, the crystalline anisotropies derived in the current work are consistent with earlier calculations.     

Three-dimensional global analysis of thermal stress and dislocations in a silicon ingot during a unidirectional solidification process with a square crucible

A three-dimensional global model was used to obtain the solution of a thermal field within the entire furnace during a unidirectional solidification process of multicrystalline silicon with a square crucible. Then the thermal stress distribution in the silicon ingot was solved. Based on the solution of thermal stress, relaxation of stress and multiplication of dislocations were performed by using the Haasen–Alexander–Sumino model (HAS model). The influence of crucible constraint on stress levels and dislocations was investigated. It was found that the crucible constraint had significant influence on the thermal stresses and dislocations in the ingot. The results indicated that it is important to reduce the crucible constraint in order to relax thermal stresses and reduce dislocations in a silicon ingot during the solidification process.     

Efficient implementation of the proper outlet flow conditions in blood flow simulations through asymmetric arterial bifurcations

We present an efficient implementation of the proper (in vivo) outlet boundary conditions in detailed, three‐dimensional (3D) and time‐periodic simulations of blood flow through arteries. This is achieved through the intermediate use of an approximate ‘simulant’ model of the outlet pressure/flow relationship corresponding to the full 3D and time‐dependent numerical simulation. This model allows us to efficiently couple the 3D outlet pressure/flow conditions to the equivalent relations due to the downstream arterial network, as obtained from a one‐dimensional approximate model in the form of Fourier frequency impedance coefficients. An adjustable time‐periodic function correction term in the simulant model requires input from the full 3D model that has to run iteratively until convergence. The advantage of the proposed numerical scheme is that it decouples the upstream detailed simulation from the downstream approximate network model offering exceptional versatility. This approach is demonstrated here in a series of detailed 3D simulations of blood flow, performed using the commercial software FLUENT?, through an asymmetric arterial bifurcation. Two cases are considered: first a healthy system patterned after the left main coronary arterial bifurcation, and second a diseased case where an occlusion has developed in one of the daughter vessels, resulting in strengthening the asymmetry of the bifurcation. Rapid convergence of the iterative process was achieved in both cases. Subtle changes occur in the shear patterns of the daughter vessels, whereas the flow distribution is quite different. In the presence of a stenosis additional regions of low shear develop due to inertial effects. Copyright © 2010 John Wiley & Sons, Ltd.     

On the relative impact of subgrid‐scale modelling and conjugate heat transfer in LES of hot jets in cross‐flow over cold plates

This work describes numerical simulations of a hot jet in cross‐flow with applications to anti‐ice systems of aircraft engine nacelles. Numerical results are compared with experimental measurements obtained at ONERA to evaluate the performances of LES in this industrial context. The combination of complex geometries requiring unstructured meshes and high Reynolds number does not allow the resolution of boundary layers so that wall models must be employed. In this framework, the relative influence of subgrid‐scale modelling and conjugate heat transfer in LESs of aerothermal flows is evaluated. After a general overview of the transverse jet simulation results, a LES coupled with a heat transfer solver in the walls is used to show that thermal boundary conditions at the wall have more influence on the results than subgrid scale models. Coupling fluid flow and heat transfer in solids simulations is the only method to specify their respective thermal boundary conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Understanding the interaction of multi-walled carbon nanotubes with mutagenic organic pollutants using computational modeling and biological experiments

Carbon nanotubes (CNTs) are very promising materials to remove pollutants from the environment. To develop safe, efficient technologies, it is necessary to understand the mechanisms of interaction between CNTs and pollutants. This requires innovative, interdisciplinary approaches. Detailed chemical analysis of the CNTs along with computational modeling can provide important information about the mechanisms of interaction. If biological experiments are included in these studies, useful complementary information is obtained. To exemplify the use of this approach, we present a case study in which detailed calculations and the Salmonella mutagenicity assay were applied to elucidate how multi-walled CNTs interact with 1-nitropyrene, an important mutagenic pollutant.     

The influence of lattice structure on the dynamic performance of metal hollow sphere agglomerates

This paper addresses the dynamic performances of partially bonded or unbonded metallic hollow sphere structure (MHSS). Under the condition that the densities of MHSS are the same, three morphologies, that is, simple cubic (SC), body-centered cubic (BCC) and face-centered cubic (FCC) are investigated to see the topological effect of lattice structure on the dynamic performance of MHS agglomerates.     

Simulation of oscillatory flow in an aortic bifurcation using FVM and FEM: A comparative study of implementation strategies

Cardiovascular diseases are one of the major causes of long‐term morbidity and mortality in human beings. The nearly epidemic increase in prevalence of such diseases poses a serious threat to public health and calls for efficient methods of diagnosis and treatment. Non‐invasive diagnostic procedures such as MRI are often used in this context; however, these are limited in terms of spatial and temporal resolution and do not provide information on time‐dependent pressures and wall shear stresses—key quantities considered to be partially responsible for the formation and development of related pathologies. The present study is concerned with the numerical simulation of oscillatory flow through the abdominal aortic bifurcation. Computational fluid dynamics simulation of oscillatory flow in a branched geometry at high Reynolds numbers poses considerable challenges. The present study reports a detailed comparison of simulations performed with a finite volume and a finite element method, two approaches with significant differences in their discretization strategy, treatment of boundary conditions and other numerical aspects. Both solvers were parallelized, using loop parallelization of the BiCGStab linear solver for the finite volume and domain decomposition based on the Schur complement method for the finite element technique. The experience gained with these two approaches for the solution of flow in a bifurcation forms the focus of this study. Although similar results were obtained for both methods, the computation time required for convergence was found to be significantly smaller for the finite element approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Liquid-crystal phase equilibria of Lennard-Jones chains

Liquid-crystal phase equilibria of Lennard-Jones chain fluids and the solubility of a Lennard-Jones gas in the coexisting phases are calculated from Monte Carlo simulations. Direct phase equilibria calculations are performed using an expanded formulation of the Gibbs ensemble. Monomer densities, order parameters, and equilibrium pressures are reported for the coexisting isotropic and nematic phases of: (1) linear Lennard-Jones chains, (2) a partially-flexible Lennard-Jones chain, and (3) a binary mixture of linear Lennard-Jones chains. The effect of chain length is determined by calculating the isotropic-nematic coexistence of linear Lennard-Jones chain fluids made of 8, 10, and 12 segments (8-, 10-, 12-mer). The effect of molecular flexibility on the isotropic-nematic equilibrium is studied for a Lennard-Jones 10-mer chain fluid with one freely-jointed segment at the end of the chain. An isotropic-nematic phase split and fractionation are reported for a binary mixture of linear 7-mer and 12-mer chains. Simulation results are compared with theoretical results as obtained from a recently developed analytical equation of state based on perturbation theory. Excellent agreement between theory and simulations is observed. The solubility of a monomer Lennard-Jones gas in the coexisting isotropic and nematic phases is estimated using the Widom test-particle insertion method. A linear relationship between solubility difference and density difference at isotropic-nematic coexistence is observed. It is shown that gas solubility is independent of the nematic ordering of the fluid, at constant temperature and density conditions.     

机载平台式惯导系统导电环故障机理分析及解决措施

导电环是平台式惯导系统惯性平台的重要组成部分,主要承担惯性平台台体内各种信号与系统的传输和交联,导电环发生故障直接会导致惯导系统失效。针对导电环三种故障模式:接触电阻变大、绝缘强度降低以及烧坏,本文采用微动磨损和电接触理论,结合作者的实际工作经验,对故障机理进行了分析,提出了相应的改进措施。这些措施均在生产中得到落实。落实措施后的故障统计表明,采取的改进措施有效,导电环的故障率与原来相比降低了38%。