Reynolds-stress turbulence model in the KIVA code for engine simulation

A Reynolds-Stress Turbulence Model has been incorporated with success into the KIVA code, a computational fluid dynamics hydrocode for three-dimensional simulation of fluid flow in engines. The newly implemented Reynolds-stress turbulence model greatly improves the robustness of KIVA, which in its original version has only eddy-viscosity turbulence models. Validation of the Reynolds-stress turbulence model is accomplished by conducting pipe-flow and channel-flow simulations, and comparing the computed results with experimental and direct numerical simulation data. Flows in engines of various geometry and operating conditions are calculated using the model, to study the complex flow fields as well as confirm the model’s validity. Results show that the Reynolds-stress turbulence model is able to resolve flow details such as swirl and recirculation bubbles. The model is proven to be an appropriate choice for engine simulations, with consistency and robustness, while requiring relatively low computational effort.

In-situ electrical, mechanical and electrochemical characterizations of one-dimensional nanostructures

One-dimensional nanostructures initiated new aspects to the materials applications due to their superior properties compared to the bulk materials. Properties of nanostructures have been characterized by many techniques and used for various device applications. However, simultaneous correlation between the physical and structural properties of these nanomaterials has not been widely investigated. Therefore, it is necessary to perform in-situ study on the physical and structural properties of nanomaterials to understand their relation. In this work, we will use a unique instrument to perform real time atomic force microscopy (AFM) and scanning tunneling microscopy (STM) of nanomaterials inside a transmission electron microscopy (TEM) system. This AFM/STM-TEM system is used to investigate the mechanical, electrical, and electrochemical properties of boron nitride nanotubes (BNNTs) and Silicon nanorods (SiNRs). BNNTs are one of the subjects of this PhD research due to their comparable, and in some cases superior, properties compared to carbon nanotubes. Therefore, to further develop their applications, it is required to investigate these characteristics in atomic level. In this research, the mechanical properties of multi-walled BNNTs were first studied. Several tests were designed to study and characterize their real-time deformation behavior to the applied force. Observations revealed that BNNTs possess highly flexible structures under applied force. Detailed studies were then conducted to understand the bending mechanism of the BNNTs. Formations of reversible ripples were observed and described in terms of thermodynamic energy of the system. Fracture failure of BNNTs were initiated at the outermost walls and characterized to be brittle. Second, the electrical properties of individual BNNTs were studied. Results showed that the bandgap and electronic properties of BNNTs can be engineered by means of applied strain. It was found that the conductivity, electron concentration and carrier mobility of BNNTs can be tuned as a function of applied stress. Although, BNNTs are considered to be candidate for field emission applications, observations revealed that their properties degrade upon cycles of emissions. Results showed that due to the high emission current density, the temperature of the sample was increased and reached to the decomposition temperature at which the B-N bonds start to break. In addition to BNNTs, we have also performed in-situ study on the electrochemical properties of silicon nanorods (SiNRs). Specifically, lithiation and delithiation of SiNRs were studied by our STM-TEM system. Our observations showed the direct formation of Li22Si5 phases as a result of lithium intercalation. Radial expansion of the anode materials were observed and characterized in terms of size-scale. Later, the formation and growth of the lithium fibers on the surface of the anode materials were observed and studied. Results revealed the formation of lithium islands inside the ionic liquid electrolyte which then grew as Li dendrite toward the cathode material.

Modeling mechanical response of heterogeneous materials

Heterogeneous materials are ubiquitous in nature and as synthetic materials. These materials provide unique combination of desirable mechanical properties emerging from its heterogeneities at different length scales. Future structural and technological applications will require the development of advanced light weight materials with superior strength and toughness. Cost effective design of the advanced high performance synthetic materials by tailoring their microstructure is the challenge facing the materials design community. Prior knowledge of structure-property relationships for these materials is imperative for optimal design. Thus, understanding such relationships for heterogeneous materials is of primary interest. Furthermore, computational burden is becoming critical concern in several areas of heterogeneous materials design. Therefore, computationally efficient and accurate predictive tools are highly essential. In the present study, we mainly focus on mechanical behavior of soft cellular materials and tough biological material such as mussel byssus thread. Cellular materials exhibit microstructural heterogeneity by interconnected network of same material phase. However, mussel byssus thread comprises of two distinct material phases. A robust numerical framework is developed to investigate the micromechanisms behind the macroscopic response of both of these materials. Using this framework, effect of microstuctural parameters has been addressed on the stress state of cellular specimens during split Hopkinson pressure bar test. A voronoi tessellation based algorithm has been developed to simulate the cellular microstructure. Micromechanisms (microinertia, microbuckling and microbending) governing macroscopic behavior of cellular solids are investigated thoroughly with respect to various microstructural and loading parameters. To understand the origin of high toughness of mussel byssus thread, a Genetic Algorithm (GA) based optimization framework has been developed. It is found that two different material phases (collagens) of mussel byssus thread are optimally distributed along the thread. These applications demonstrate that the presence of heterogeneity in the system demands high computational resources for simulation and modeling. Thus, Higher Dimensional Model Representation (HDMR) based surrogate modeling concept has been proposed to reduce computational complexity. The applicability of such methodology has been demonstrated in failure envelope construction and in multiscale finite element techniques. It is observed that surrogate based model can capture the behavior of complex material systems with sufficient accuracy. The computational algorithms presented in this thesis will further pave the way for accurate prediction of macroscopic deformation behavior of various class of advanced materials from their measurable microstructural features at a reasonable computational cost.

Effects of silicon content and cooling rate on mechanical properties of heavy section ductile cast iron

The effects of Si and cooling rate are investigated for their effect on the mechanical properties and microstructure. Three alloys were chosen with varying C and Si contents and an attempt to keep the remainder of the elements present constant. Within each heat, three test blocks were poured. Two blocks had chills – one with a fluid flowing through it to cool it (active chill) and one without the fluid (passive) – and the third block did not have a chill. Cooling curves were gathered and analyzed. The mechanical properties of the castings were correlated to the microstructure, cooling rate and Si content of each block. It was found that an increase in Si content increased the yield stress, tensile strength and hardness but decreased the impact toughness, elongation and Young’s modulus. The fast cooling rates produced by the chills caused a high nodule count in the castings along with a fine ferrite grain size and a high degree of nodularity. The fine microstructures, in turn, increased the strength and ductile to brittle transition temperature (DBTT) of the castings. The fast cooling rate was not adequate to overcome the dramatic increase in DBTT that is caused by the addition of Si.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

IN-SITU ELECTRICAL, MECHANICAL AND ELECTROCHEMICAL CHARACTERIZATIONS OF ONE-DIMENSIONAL NANOSTRUCTURES

One-dimensional nanostructures initiated new aspects to the materials applications due to their superior properties compared to the bulk materials. Properties of nanostructures have been characterized by many techniques and used for various device applications. However, simultaneous correlation between the physical and structural properties of these nanomaterials has not been widely investigated. Therefore, it is necessary to perform in-situ study on the physical and structural properties of nanomaterials to understand their relation. In this work, we will use a unique instrument to perform real time atomic force microscopy (AFM) and scanning tunneling microscopy (STM) of nanomaterials inside a transmission electron microscopy (TEM) system. This AFM/STM-TEM system is used to investigate the mechanical, electrical, and electrochemical properties of boron nitride nanotubes (BNNTs) and Silicon nanorods (SiNRs). BNNTs are one of the subjects of this PhD research due to their comparable, and in some cases superior, properties compared to carbon nanotubes. Therefore, to further develop their applications, it is required to investigate these characteristics in atomic level. In this research, the mechanical properties of multi-walled BNNTs were first studied. Several tests were designed to study and characterize their real-time deformation behavior to the applied force. Observations revealed that BNNTs possess highly flexible structures under applied force. Detailed studies were then conducted to understand the bending mechanism of the BNNTs. Formations of reversible ripples were observed and described in terms of thermodynamic energy of the system. Fracture failure of BNNTs were initiated at the outermost walls and characterized to be brittle. Second, the electrical properties of individual BNNTs were studied. Results showed that the bandgap and electronic properties of BNNTs can be engineered by means of applied strain. It was found that the conductivity, electron concentration and carrier mobility of BNNTs can be tuned as a function of applied stress. Although, BNNTs are considered to be candidate for field emission applications, observations revealed that their properties degrade upon cycles of emissions. Results showed that due to the high emission current density, the temperature of the sample was increased and reached to the decomposition temperature at which the B-N bonds start to break. In addition to BNNTs, we have also performed in-situ study on the electrochemical properties of silicon nanorods (SiNRs). Specifically, lithiation and delithiation of SiNRs were studied by our STM-TEM system. Our observations showed the direct formation of Li22Si5 phases as a result of lithium intercalation. Radial expansion of the anode materials were observed and characterized in terms of size-scale. Later, the formation and growth of the lithium fibers on the surface of the anode materials were observed and studied. Results revealed the formation of lithium islands inside the ionic liquid electrolyte which then grew as Li dendrite toward the cathode material.