Engineering topics

Stroud, D.

Publications and source records attributed to Stroud, D..

10 records

NASA NTRS1991

Bond-orientational order in liquid Si

Bond-orientational order in liquid Si via Monte Carlo simulation in conjuncation with empirical two- and three-body potentials of the form proposed by Stillinger and Weber are studied. Bond-orientational order (BOO) is described in terms of combinations of spherical harmonic functions. Liquid Si is found to have pronounced short-range BOO corresponding to l = 3, as expected for a structure with local tetrahedral order. No long-range BOO is found either in the equilibrium or the supercooled liquid. When the three-body potential is artificially removed, the tetrahedral bond-orientation order disappears and the liquid assumes a close-packed structure.

Wang, Z. Q.↗

NASA NTRS1990

Surface entropy of liquids via a direct Monte Carlo approach - Application to liquid Si

Two methods are presented for a direct Monte Carlo evaluation of the surface entropy S(s) of a liquid interacting by specified, volume-independent potentials. The first method is based on an application of the approach of Ferrenberg and Swendsen (1988, 1989) to Monte Carlo simulations at two different temperatures; it gives much more reliable results for S(s) in liquid Si than previous calculations based on numerical differentiation. The second method expresses the surface entropy directly as a canonical average at fixed temperature.

Wang, Z. Q.↗

NASA NTRS1989

Ginzburg-Landau theory for the solid-liquid interface of bcc elements. II - Application to the classical one-component plasma, the Wigner crystal, and He-4

The previously developed Ginzburg-Landau theory for calculating the crystal-melt interfacial tension of bcc elements to treat the classical one-component plasma (OCP), the charged fermion system, and the Bose crystal. For the OCP, a direct application of the theory of Shih et al. (1987) yields for the surface tension 0.0012(Z-squared e-squared/a-cubed), where Ze is the ionic charge and a is the radius of the ionic sphere. Bose crystal-melt interface is treated by a quantum extension of the classical density-functional theory, using the Feynman formalism to estimate the relevant correlation functions. The theory is applied to the metastable He-4 solid-superfluid interface at T = 0, with a resulting surface tension of 0.085 erg/sq cm, in reasonable agreement with the value extrapolated from the measured surface tension of the bcc solid in the range 1.46-1.76 K. These results suggest that the density-functional approach is a satisfactory mean-field theory for estimating the equilibrium properties of liquid-solid interfaces, given knowledge of the uniform phases.

Zeng, X. C.↗

NASA NTRS1988

Monte Carlo studies of liquid semiconductor surfaces - Si and Ge

The liquid-vapor interface of Si and Ge, and of Si doped with impurities is studied using the empirical Stillinger-Weber two- and three-body potentials. The surface tension of the pure elements is calculated by use of a direct evaluation of the free energy required to create the surface. For Si, both the surface tension and its temperature derivative are in good agreement with experiment. To within numerical accuracy, the free surfaces of both Si and Ge have a monotonically decreasing density, with a 10-90 percent surface width of about 2.2 A in both cases. When large and small model impurities are introduced into pure liquid Si, they are found in the simulations to migrate, respectively, towards the surface and away from the surface. This behavior is consistent with the interpretation that impurities with atoms larger than Si tend to lower the surface tension.

Wang, Z. Q.↗

NASA NTRS1987

Ginzburg-Landau theory for the solid-liquid interface of bcc elements

Consideration is given to a simple order-parameter theory for the interfacial tension of body-centered-cubic solids in which the principal order parameter is the amplitude of the density wave at the smallest nonzero reciprocal-lattice vector of the solid. The parameters included in the theory are fitted to the measured heat of fusion, melting temperature, and solid-liquid density difference, and to the liquid structure factor and its temperature derivative at freezing. Good agreement is found with experiment for Na and Fe and the calculated anisotropy of the surface tension among different crystal faces is of the order of 2 percent. On the basis of various assumptions about the universal behavior of bcc crystals at melting, the formalism predicts that the surface tension is proportional to the heat of fusion per surface atom.

Shih, W. H.↗

NASA NTRS1987

Calculation of the surface tension of liquid metals using a one-component-plasma reference system

The one-component-plasma (OCP) model is used as a reference system instead of the traditional hard-sphere fluid to calculate the liquid-vapor interfacial surface tension of liquid metals within the density functional formalism. The calculated surface tensions of the alkali metals are in excellent agreement with experiment. For the polyvalent metal Al, the result obtained is larger than experimental measurements. It is concluded that the OCP system is not suitable to describe the liquid-vapor phase transition in simple metals which have a nominal plasma parameter larger than the usual freezing value of about 178. The calculated interfacial widths in all cases are narrower than the expected experimental values.

Zeng, X. C.↗

NASA NTRS1985

Thermodynamic properties of liquid Si and Ge

Thermodynamic properties of liquid Si and Ge are calculated using standard variational techniques, with the hard-sphere system as a reference system. Third-order electron-ion pseudopotential contributions are included nonvariationally and are shown to have little effect on the free energy, although a larger one on the derivatives of the free energy. Two different model pseudopotentials are used and give similar results in good agreement with experiment.

Shih, W.-H.↗

NASA NTRS1984

Theoretical study of miscibility and glass-forming trends in mixtures of polystyrene spheres

A theoretical study of glass-forming trends and miscibility in mixtures of polystyrene spheres (polyballs) of different diameters, suspended in an aqueous solution, is presented. The polyballs are assumed to be charged and to interact via a Debye-Hueckel screened Coulomb potential. The Helmholtz free energy is calculated from a variational principle based on the Gibbs-Bogoliubov inequality, in which a mixture of hard spheres of different diameters is chosen as the reference system. It is found that when the charges of the two types of polyballs are sufficiently different, the variationally determined ratio of hard-sphere diameters differs substantially, leading to packing difficulties characteristic of glass formation. The experimentally observed range of glass formation corresponds to a ratio of hard-sphere diameters of 0.8 or less. Calculations of the free energy as a function of concentration indicate that the liquid polyball mixture is stable against the phase separation, even for widely different polyball charges.

Shih, W.-H.↗

NASA NTRS1984

Density-functional calculations of the surface tension of liquid Al and Na

Calculations of the surface tensions of liquid Al and Na are described using the full ionic density functional formalism of Wood and Stroud (1983). Surface tensions are in good agreement with experiment in both cases, with results substantially better for Al than those found previously in the gradient approximation. Preliminary minimization with respect to surface profile leads to an oscillatory profile superimposed on a nearly steplike ionic density disribution; the oscillations have a wavellength of about a hardsphere diameter.

Stroud, D.↗

NASA NTRS1983

Theoretical study of the freezing of polystyrene sphere suspensions

A theoretical study of melting in aqueous suspensions of polystyrene spheres is presented. The Helmholtz free energies of the liquid and solid phase of the suspensions are calculated as functions of sphere number density and electrolyte concentrations. The results tend to show that the freezing curve of such suspensions is that of a conventional classical liquid with repulsive short-range interactions. The interactions can be treated satisfactorily within a slightly modified Debye-Hueckel approximation. The modifications include size correction and a correction for nonlinear screening. The results are confirmed by analogous calculations for the solid phase, and for the line along which liquid and solid free energies are equal.

Shih, W.-H.↗