Variation of added mass and its application to the calculation of amplitude response for a circular cylinder

In the present study, analyzed are the variation of added mass for a circular cylinder in the lock-in ( synchronization) range of vortex-induced vibration (VIV) and the relationship between added mass and natural frequency. A theoretical minimum value of the added mass coefficient for a circular cylinder at lock-in is given. Developed are semi-empirical formulas for the added mass of a circular cylinder at lock-in as a function of flow speed and mass ratio. A comparison between experiments and numerical simulations shows that the semi-empirical formulas describing the variation of the added mass for a circular cylinder at lock-in are better than the ideal added mass. In addition, computation models such as the wake oscillator model using the present formulas can predict the amplitude response of a circular cylinder at lock-in more accurately than those using the ideal added mass.

Measurement and Calculation of Vibrational Temperature behind Strong Shock Waves

In the present paper the measured values of vibrational temperature behind strong shock waves are compared with theoretical ones. The histories of vibrational temperature behind strong shock waves in a shock tube were measured using two monochromators. The test gas was pure nitrogen at 100-300Pa, and the speeds of shock waves were 5.0-6.0km/s. The electronic temperature of N-2(+) was also approximately determined from experiment and compared with the experimental vibrational temperature. The results show that the presented calculational method is effective, and the electronic energy of N2+ is excited much faster than its vibrational energy. One Langmuir probe was used to determine the effective time of region 2. The influence of viscosity in the shock tube is also analyzed.

Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

Part I. Approximate Hartree-Fock wavefunctions one-electron properties and electronic structure of the water molecule. Part II. Perturbation-variational calculation of the nuclear spin-spin isotope coupling constant in HD

Part I

Several approximate Hartree-Fock SCF wavefunctions for the ground electronic state of the water molecule have been obtained using an increasing number of multicenter s, p, and d Slater-type atomic orbitals as basis sets. The predicted charge distribution has been extensively tested at each stage by calculating the electric dipole moment, molecular quadrupole moment, diamagnetic shielding, Hellmann-Feynman forces, and electric field gradients at both the hydrogen and the oxygen nuclei. It was found that a carefully optimized minimal basis set suffices to describe the electronic charge distribution adequately except in the vicinity of the oxygen nucleus. Our calculations indicate, for example, that the correct prediction of the field gradient at this nucleus requires a more flexible linear combination of p-orbitals centered on this nucleus than that in the minimal basis set. Theoretical values for the molecular octopole moment components are also reported.

Part II

The perturbation-variational theory of R. M. Pitzer for nuclear spin-spin coupling constants is applied to the HD molecule. The zero-order molecular orbital is described in terms of a single 1s Slater-type basis function centered on each nucleus. The first-order molecular orbital is expressed in terms of these two functions plus one singular basis function each of the types e^-r/r and e^-r ln r centered on one of the nuclei. The new kinds of molecular integrals were evaluated to high accuracy using numerical and analytical means. The value of the HD spin-spin coupling constant calculated with this near-minimal set of basis functions is J_HD = +96.6 cps. This represents an improvement over the previous calculated value of +120 cps obtained without using the logarithmic basis function but is still considerably off in magnitude compared with the experimental measurement of J_HD = +43 0 ± 0.5 cps.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.

Theoretical investigation on the one- and two-photon absorption properties for a series of symmetrical charge transfer fluorene-phenylene and fluorene-thiophene derivatives

One- and two-photon absorption properties of a series of fluorene derivatives with symmetrical charge transfer D-IT-D and A-IT-A structural motifs have been theoretically investigated with ZINDO/S method. The optimized structures and the characterization of frontier molecular orbitals were obtained by using AMI calculations. Two-photon absorption properties of molecules have been studied using three-state model. The calculation results have shown that fluorene-thiophene derivatives exhibit larger two-photon absorption cross-section as compared with other studied molecules. To illustrate the results, the crucial effects of thiophene ring on fluorenethiophene derivatives and the net charge changes on the pi-conjugated bridges are analyzed theoretically. (c) 2006 Elsevier B.V. All rights reserved.

Theoretical analysis of optical properties of dielectric coatings dependence on substrate subsurface defects

A model for refractive index of stratified dielectric substrate was put forward according to theories of inhomogeneous coatings. The substrate was divided into surface layer, subsurface layer and bulk layer along the normal direction of its surface. Both the surface layer (separated into N-1 sublayers of uniform thickness) and subsurface layer (separated into N-2 sublayers of uniform thickness), whose refractive indices have different statistical distributions, are equivalent to inhomogeneous coatings, respectively. And theoretical deduction was carried Out by employing characteristic matrix method of optical coatings. An example of mathematical calculation for optical properties of dielectric coatings had been presented. The computing results indicate that substrate subsurface defects can bring about additional bulk scattering and change propagation characteristic in thin film and Substrate. Therefore, reflectance, reflective phase shift and phase difference of an assembly of coatings and substrate deviate from ideal conditions. The model will provide some beneficial theory directions for improving optical properties of dielectric coatings via substrate surface modification. (c) 2005 Elsevier B.V. All rights reserved.

Theoretical analysis of influence of substrate microdefects on polarized light properties of dielectric coatings

A model for refractive index of stratified dielectric substrate was presented according to inhomogeneous coatings theories. The substrate was divided into surface layer, subsurface layer and bulk layer along the normal direction of its surface. The former two layers were equivalent to inhomogeneous coatings. Theoretical deduction was executed by employing the characteristic matrix method of optical coatings, and one mathematical calculation example was presented. The results indicate that reflectance, reflective phase shift and phase difference of polarized light deviate from ideal conditions. It shows that substrate microdefects can induce volume scattering and change propagation characteristic of light both in coatings and substrate. (c) 2005 Elsevier GmbH. All rights reserved.

Calculation and measurement of the magnetostriction of multilayered composites

A new theoretical model that predicts the magnetostriction of multilayered composites has been developed. The model takes into account the shear stress between the composite layers and consequently predicts a non-uniform strain along their thickness. The model has been experimentally validated by producing composites formed from three materials with different magnetostrains and mechanical properties, and controlled layer thicknesses in the order of micrometers. Deformations of several ppm, up to 7.5% of the saturation magnetostrain were measured between the edge and the centre of such composites. © 2006 Elsevier B.V. All rights reserved.

Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.

Theoretical investigation on the phase stability of Nd2Fe14B and site preference of V, Cr, Mn, Zr and Nb

The stability of the excellent permanent magnetic compound Nd2Fe14B and substitution of Fe in the compound by V, Cr, Mn, Zr and Nb are investigated by using interatomic pair potentials which are converted from lattice-inversion method. Calculation shows that the substitution always makes the cell volume larger, and the increase of the volume is almost linear with substituent concentration. The calculated cohesive energy shows that the preferential order of substitution of Fe is Nb, V, Cr, Mn, Zr. Nevertheless, all the five substituting elements should most preferentially replace Fe in the j(2)' site, which has the greatest space among all six Fe sites. (C) 2005 Elsevier B.V. All rights reserved.

Theoretical analysis of gate voltage-controlled subband states in an AlxGa1-xN/GaN heterostructure

The subband structure and inter-subband transition as a function of gate voltage are determined by solving the Schrodinger and Poisson equations self-consistently in an AlxGa1-xN/GaN heterostructure. Different aluminum mole fraction and thickness of AlxGa1-xN barrier are considered. Calculation results show that energy difference between the first and second subband covers a wide range (from several tens to hundreds milli-electron volt) by applying different gate voltage, which corresponds to the midinfrared and long-wave infrared wavelength scope. Furthermore, such a modulation on the subband transition energy is much more pronounced for the structure with thin barrier. When the applied positive gate voltage is increased, the triangle well formed at the interface turns to be deeper and narrower, which enhances the confinement for electrons. As a result, the overlap between electron wave function at two subbands increases, and thus the optical intersubband transition also enhances its intensity. This tendency is in good agreement with the available data in the literature. (c) 2005 Elsevier B.V. All rights reserved.

Theoretical analysis about the influence of channel layer thickness on the 2D electron gas and its distribution in InP-based high-electron-mobility transistors

The principle of high-electron-mobility transistor (HEMT) and the property of two-dimensional electron gas (2DEG) have been analyzed theoretically. The concentration and distribution of 2DEG in various channel layers are calculated by numerical method. Variation of 2DEG concentration in different subband of the quantum well is discussed in detail. Calculated results show that sheet electron concentration of 2DEG in the channel is affected slightly by the thickness of the channel. But the proportion of electrons inhabited in different subbands can be affected by the thickness of the channel. When the size of channel lies between 20-25 nm, the number of electrons occupying the second subband reaches the maximum. This result can be used in parameter design of materials and devices.

Theoretical analysis of characteristics of GaxIn1-xNyAs1-y/GaAs quantum well lasers with different intermediate layers

Based on the band anticrossing model, the effects of the strain-compensated layer and the strain-mediated layer on the band structure, gain and differential gain of GaInNAs/GaAs quantum well lasers have been investigated. The results show that the GaNAs barrier has a disadvantage in increasing the density of states in the conduction band. Meanwhile, the multilayer quantum wells need higher transparency carrier density than the GaInNAs/GaAs single quantum well with the same wavelength. However, they help to suppress the degradation of the differential gain. The calculation also shows that from the viewpoint of band structure, the strain-compensated structure and the strain-mediated structure have similar features.

Theoretical analysis of the bandgap for the intermixed GaInP/AlGaInP quantum wells

Compositional distribution of the quantum well and barrier after quantum well intermixing for GaInP/AlGaInP system was theoretically analyzed on the basis of atom diffusion law. With the compositional distribution result, the valence subband structure of the intermixed quantum well was calculated on the basis of 6x6 Luttinger-Kohn Hamiltonian, including spin-orbit splitting effects. TO get more accurate results in the calculation, a full 6-band problem was solved without axial approximation, which had been widely used in the Luttinger-Kohn model to simplify the computational efforts, since there was a strong warping in the GaInP valence band. At last, the bandgap energy of the intermixed quantum well was obtained and the calculation result is of much importance in the analysis of quantum well intermixing experiments.