Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored.

An augmented method for free boundary problems with moving contact lines

An augmented immersed interface method (IIM) is proposed for simulating one-phase moving contact line problems in which a liquid drop spreads or recoils on a solid substrate. While the present two-dimensional mathematical model is a free boundary problem, in our new numerical method, the fluid domain enclosed by the free boundary is embedded into a rectangular one so that the problem can be solved by a regular Cartesian grid method. We introduce an augmented variable along the free boundary so that the stress balancing boundary condition is satisfied. A hybrid time discretization is used in the projection method for better stability. The resultant Helmholtz/Poisson equations with interfaces then are solved by the IIM in an efficient way. Several numerical tests including an accuracy check, and the spreading and recoiling processes of a liquid drop are presented in detail. (C) 2010 Elsevier Ltd. All rights reserved.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Direct Numerical Simulation of Shock/Turbulent Boundary Layer Interaction in a Supersonic Compression Ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9. The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition. Both the mean wall pressure and the velocity profiles agree with those of the experimental data, which validates the simulation. The turbulent kinetic energy budget in the separation region is analyzed. Results show that the turbulent production term increases fast in the separation region, while the turbulent dissipation term reaches its peak in the near-wall region. The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation. Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble, the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance, but rather the instability of separation bubble.

Effects of Magnetohydrodynamic Interaction-Zone Position on Shock-Wave/Boundary-Layer Interaction

A third-order weighted essentially nonoscillatory and non-free-parameter difference scheme magnetohydrodynamic solver has been established to investigate the mechanisms of magnetohydrodynamics controlling separation induced by an oblique shock wave impinging on a flat plate. The effects of magnetohydrodynamic interaction-zone location on the separation point, reattachment point, separation-bubble size, and boundary-layer velocity profiles are analyzed. The results show that there exists a best location for the magnetohydrodynamic zone to be applied, where the separation point is delayed the farthest, and the separation bubble is decreased up to about 50% in size compared to the case without magnetohydrodynamic control, which demonstrated the promising of magnetohydrodynamics suppressing the separation induced by shock-wave/boundary-layer interactions.

Boundary Characterization Experiment Series Overview

Ocean acoustic propagation and reverberation in continental shelf regions is often controlled by the seabed and sea surface boundaries. A series of three multi-national and multi-disciplinary experiments was conducted between 2000-2002 to identify and measure key ocean boundary characteristics. The frequency range of interest was nominally 500-5000 Hz with the main focus on the seabed, which is generally considered as the boundary of greatest importance and least understood. Two of the experiments were conducted in the Mediterranean in the Strait of Sicily and one experiment in the North Atlantic with sites on the outer New Jersey Shelf (STRATAFORM area) and on the Scotian Shelf. Measurements included seabed reflection, seabed, surface, and biologic scattering, propagation, reverberation, and ambient noise along with supporting oceanographic, geologic, and geophysical data. This paper is primarily intended to provide an overview of the experiments and the strategies that linked the various measurements together, with detailed experiment results contained in various papers in this volume and other sources

extending the boundary of spreadsheet programming: lessons learned from chinese governmental projects

ICSE

Effects of Pr dopant on grain boundary and electrical properties of Ce5.2Sm0.8MoO15-delta

Material formulated as Ce5.2Sm0.8-xPrxMo15-(delta) (x=0.08) was prepared by adding small amounts of Pr dopant in oxide Ce5.2SM0.8-xPrxMoO15-delta. Structural and electrical properties were investigated by means of X-ray diffraction (XRD), Raman spectroscopy, X-ray photoelectron spectroscopy (XPS), field emission scanning electron microscopy (FE-SEM) and AC impedance spectroscopy. The effect of small amounts of Pr on microstructure and electrical conductivity was discussed. It was showed that the material doped with Pr has a lot of dents and small openings, which provide channels for oxygen ions, resulting in lower grain boundary and total conductivity activation energy. Thus the corresponding grain boundary conductivity and total conductivity of the material were improved notably. The grain boundary conductivity of the material doped with Pr is 6.79 X 10(-3) S center dot cm(-1) at 500 degrees C, which is twice as large as that without Pr (5.61 X 10(-5) S center dot cm(-1)).

Grain boundary conductivity of high purity neodymium-doped ceria nanosystem with and without the doping of molybdenum oxide

Solid solutions of Ce1-xNdxO2-x/2 (0.05 <= x <= 0.2) and (Ce1-xNdx)(0.95)MO0.05O2-delta (0.05 <= x <= 0.2) have been synthesized by a modified sol-gel method. Both materials have very low content of SiO2 (similar to 27 ppm). Their structures and ionic conductivities were characterized by X-ray diffraction (XRD), field-emission scanning electron microscopy (FESEM) and electrochemical impedance spectroscopy (M). The XRD patterns indicate that these materials are single phases with a cubic fluorite structure. The powders calcined at 300 degrees C with a crystal size of 5.7 nm have good sinterability, and the relative density could reach above 96% after being sintered at 1450 degrees C. With the addition Of MoO3, the sintering temperature could be decreased to 1250 degrees C. Impedance spectroscopy measurement in the temperature range of 250-800 degrees C indicates that a sharp increase of conductivity is observed when a small amount of Nd2O3 is added into ceria, of which Ce0.85Nd0.15O1.925 (15NDC) shows the highest conductivity. With the addition of a small amount Of MoO3, the grain boundary conductivity of 15NDC at 600 degrees C increases from 2.56 S m(-1) to 5.62 S m(-1).

Boundary effect of relief structure on crystallization of diblock copolymer in thin films

In this Letter, crystal growth of a symmetric crystalline-amorphous diblock copolymer, poly(styrene-b-epsilon-caprolactone) (PS-b-PCL), in thin films was investigated by atomic force microscopy (AFM), Relief structures of holes and islands were formed during annealing the film at the molten state, and the in situ observation of subsequent crystal growth at room temperature indicated that the crystals were preferred to occur at the edge of holes or islands and grew into the interior area. It was concluded that the stretched PCL blocks at the edge of relief structures, caused by material transportation or deformation of the interface, could act as nucleation agents during polymer crystallization. The crystal growth rate of individual lamellae varied both from lamellae to lamellae and in time, but the area occupied by crystals increased constantly with time. At 22 degreesC, the growth rate was 1.2 x 10(-2) mum(2)/min with the scan size 2 x 2 mum(2).

Growth rates of different polymorphs from interspherulitic boundary profiles

Quantitative data on the crystallization kinetics of polymorphic polymers can be derived from the investigation of gross spherulitic morphology formed in isothermal conditions. Depending on distance between centers, and the time lag between their formation and relative growth rates, various types of boundary lines can be generated by the impinging of two spherical bodies whose radii increase linearly with time, In polymorphic polymers, different types of spherulites often develop simultaneously at different rates from sporadic or predetermined nuclei. In same cases, the so-called growth transformation, in which a nucleus of the fast growing specie is formed at the tip of an advancing lamella of the slower crystal form, provides an alternative mode of nucleation, It is shown that if only one event of growth transformation takes place at the front of a slow growing body, the fast growing spherulite swallows the parent one and the resultant shape of interspherulitic boundary is described by two symmetrical logarithmic spirals whose parameters can be extracted from micrographs taken at the end of crystallization. These concepts are applied to determine the radial growth rate of gamma form spherulites of polypivalolactone in a wide range of temperatures through analysis of the alpha/gamma interspherulitic profiles formed in isothermal conditions and direct measurement of the growth rate of the alpha counterparts at the same temperature.