Numerical modelling of devices and advanced transmission schemes embracing Raman effect

This thesis presents a theoretical investigation on applications of Raman effect in optical fibre communication as well as the design and optimisation of various Raman based devices and transmission schemes. The techniques used are mainly based on numerical modelling. The results presented in this thesis are divided into three main parts. First, novel designs of Raman fibre lasers (RFLs) based on Phosphosilicate core fibre are analysed and optimised for efficiency by using a discrete power balance model. The designs include a two stage RFL based on Phosphosilicate core fibre for telecommunication applications, a composite RFL for the 1.6 μm spectral window, and a multiple output wavelength RFL aimed to be used as a compact pump source for fiat gain Raman amplifiers. The use of Phosphosilicate core fibre is proven to effectively reduce the design complexity and hence leads to a better efficiency, stability and potentially lower cost. Second, the generalised Raman amplified gain model approach based on the power balance analysis and direct numerical simulation is developed. The approach can be used to effectively simulate optical transmission systems with distributed Raman amplification. Last, the potential employment of a hybrid amplification scheme, which is a combination between a distributed Raman amplifier and Erbium doped amplifier, is investigated by using the generalised Raman amplified gain model. The analysis focuses on the use of the scheme to upgrade a standard fibre network to 40 Gb/s system.

Multiphase modelling of the characteristics of close coupled gas atomization

Not withstanding the high demand of metal powder for automotive and High Tech applications, there are still many unclear aspects of the production process. Only recentlyhas supercomputer performance made possible numerical investigation of such phenomena. This thesis focuses on the modelling aspects of primary and secondary atomization. Initially two-dimensional analysis is carried out to investigate the influence of flow parameters (reservoir pressure and gas temperature principally) and nozzle geometry on final powder yielding. Among the different types, close coupled atomizers have the best performance in terms of cost and narrow size distribution. An isentropic contoured nozzle is introduced to minimize the gas flow losses through shock cells: the results demonstrate that it outperformed the standard converging-diverging slit nozzle. Furthermore the utilization of hot gas gave a promising outcome: the powder size distribution is narrowed and the gas consumption reduced. In the second part of the thesis, the interaction of liquid metal and high speed gas near the feeding tube exit was studied. Both axisymmetric andnon-axisymmetric geometries were simulated using a 3D approach. The filming mechanism was detected only for very small metal flow rates (typically obtained in laboratory scale atomizers). When the melt flow increased, the liquid core overtook the adverse gas flow and entered in the high speed wake directly: in this case the disruption isdriven by sinusoidal surface waves. The process is characterized by fluctuating values of liquid volumes entering the domain that are monitored only as a time average rate: it is far from industrial robustness and capability concept. The non-axisymmetric geometry promoted the splitting of the initial stream into four cores, smaller in diameter and easier to atomize. Finally a new atomization design based on the lesson learned from previous cases simulation is presented.

CFD and bubble column reactors:simulation and experiment

The simulation of two-phase flow in bubble columns using commercially available software fromFluent Incorporated is presented here. Data from a bubble column with a ratio of height to thecolumn diameter of 5 : 1 are compared with simulations and experimental results for time-averaged velocity and Reynolds stress proles are used to validate transient, two-dimensional simulations.The models are based on multiphase biological reactors with applications in the food industry. An example case of the mass transfer of oxygen through the liquid phase is also presented.

Molecular Dynamics implementation of BN2D or 'Mercedes Benz' water model

Two-dimensional 'Mercedes Benz' (MB) or BN2D water model (Naim, 1971) is implemented in Molecular Dynamics. It is known that the MB model can capture abnormal properties of real water (high heat capacity, minima of pressure and isothermal compressibility, negative thermal expansion coefficient) (Silverstein et al., 1998). In this work formulas for calculating the thermodynamic, structural and dynamic properties in microcanonical (NVE) and isothermal-isobaric (NPT) ensembles for the model from Molecular Dynamics simulation are derived and verified against known Monte Carlo results. The convergence of the thermodynamic properties and the system's numerical stability are investigated. The results qualitatively reproduce the peculiarities of real water making the model a visually convenient tool that also requires less computational resources, thus allowing simulations of large (hydrodynamic scale) molecular systems. We provide the open source code written in C/C++ for the BN2D water model implementation using Molecular Dynamics.

Next generation high density three-dimensional magnetic recording systems

This work is the first work using patterned soft underlayers in multilevel three-dimensional vertical magnetic data storage systems. The motivation stems from an exponentially growing information stockpile, and a corresponding need for more efficient storage devices with higher density. The world information stockpile currently exceeds 150EB (ExaByte=1x1018Bytes); most of which is in analog form. Among the storage technologies (semiconductor, optical and magnetic), magnetic hard disk drives are posed to occupy a big role in personal, network as well as corporate storage. However; this mode suffers from a limit known as the Superparamagnetic limit; which limits achievable areal density due to fundamental quantum mechanical stability requirements. There are many viable techniques considered to defer superparamagnetism into the 100's of Gbit/in2 such as: patterned media, Heat-Assisted Magnetic Recording (HAMR), Self Organized Magnetic Arrays (SOMA), antiferromagnetically coupled structures (AFC), and perpendicular magnetic recording. Nonetheless, these techniques utilize a single magnetic layer; and can thusly be viewed as two-dimensional in nature. In this work a novel three-dimensional vertical magnetic recording approach is proposed. This approach utilizes the entire thickness of a magnetic multilayer structure to store information; with potential areal density well into the Tbit/in2 regime. ^ There are several possible implementations for 3D magnetic recording; each presenting its own set of requirements, merits and challenges. The issues and considerations pertaining to the development of such systems will be examined, and analyzed using empirical and numerical analysis techniques. Two novel key approaches are proposed and developed: (1) Patterned soft underlayer (SUL) which allows for enhanced recording of thicker media, (2) A combinatorial approach for 3D media development that facilitates concurrent investigation of various film parameters on a predefined performance metric. A case study is presented using combinatorial overcoats of Tantalum and Zirconium Oxides for corrosion protection in magnetic media. ^ Feasibility of 3D recording is demonstrated, and an emphasis on 3D media development is emphasized as a key prerequisite. Patterned SUL shows significant enhancement over conventional "un-patterned" SUL, and shows that geometry can be used as a design tool to achieve favorable field distribution where magnetic storage and magnetic phenomena are involved. ^

Eulerian-Lagrangian two phase debris flow model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier-Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.

Numerical heat transfer analysis of carbon-based foams for use in thermal protection system

The applicability of carbon-based foams as an insulating or active cooling material in thermal protection systems (TPSs) of space vehicles is considered using a computer modeling. This study focuses on numerical investigation of the performance of carbon foams for use in TPSs of space vehicles. Two kinds of carbon foams are considered in this study. For active cooling, the carbon foam that has a thermal conductivity of 100 W/m-k is used and for the insulation, the carbon foam having a thermal conductivity of 0.225 W/m-k is used. A 3D geometry is employed to simulate coolant flow and heat transfer through carbon foam model. Gambit has been used to model the 3D geometry and the numerical simulation is carried out in FLUENT. Numerical results from this thesis suggests that the use of CFOAM and HTC carbon foams in TPS's may effectively protect the aluminum structure of the space shuttle during reentry of the space vehicle.

Eulerian-Lagrangian Two Phase Debris Flow Model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier–Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.

Modelagem numérica de uma estrutura de contenção de estacas espaçadas atirantadas em areia

Retaining walls design involves factors such as plastification, loading and unloading, pre-stressing, excessive displacements and earth and water thrust. Furthermore, the interaction between the retained soil and the structure is rather complex and hard to predict. Despite the advances in numerical simulation techniques and monitoring of forces and displacements with field instrumentation, design projects are still based on classical methods, whose simplifying assumptions may overestimate structural elements of the retaining wall. This dissertation involves a three-dimensional numerical study on the behavior of a retaining wall using the finite element method (FEM). The retaining wall structure is a contiguous bored pile wall with tie-back anchors. The numerical results were compared to data obtained from field instrumentation. The influence of the position of one or two layers of anchors and the effects of the construction of a slab bounded at the top of the retaining wall was evaluated. Furthermore, this study aimed at investigating the phenomenon of arching in the soil behind the wall. Arching was evaluated by analyzing the effects of pile spacing on horizontal stresses and displacements. Parametric analysis with one layers of anchors showed that the smallest horizontal displacements of the structure were achieved for between 0.3 and 0.5 times the excavation depth. Parametric analyses with two anchor layers showed that the smallest horizontal displacements were achieve for anchors positioned in depths of 0.4H and 0.7H. The construction of a slab at the top of the retaining wall decreased the horizontal displacements by 0.14% times the excavation depth as compared to analyses without the slab. With regard to the arching , analyzes showed an optimal range of spacing between the faces of the piles between 0.4 and 0.6 times the diameter of the pile

Numerical Investigation on the Heat Transfer Enhancement Using Micro/Nano Phase-Change Particulate Flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re=600-900 and phase change material particle concentrations ¡Ü0.25 , as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (0.5-0.7) slurry flow. By using the two newly-defined parameters, named effectiveness factor and performance index, respectively, it is found that there exists an optimal relation between the channel design parameters, particle volume fraction, Reynolds number, and the wall heat flux. The influence of the particle volume fraction, particle size, and the particle viscosity, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed.

Modeling, Analysis and Performance Comparison of Two Direct Sampling DCSK Receivers under Frequency Nonselective Fading Channels

Two direct sampling correlator-type receivers for differential chaos shift keying (DCSK) communication systems under frequency non-selective fading channels are proposed. These receivers operate based on the same hardware platform with different architectures. In the first scheme, namely sum-delay-sum (SDS) receiver, the sum of all samples in a chip period is correlated with its delayed version. The correlation value obtained in each bit period is then compared with a fixed threshold to decide the binary value of recovered bit at the output. On the other hand, the second scheme, namely delay-sum-sum (DSS) receiver, calculates the correlation value of all samples with its delayed version in a chip period. The sum of correlation values in each bit period is then compared with the threshold to recover the data. The conventional DCSK transmitter, frequency non-selective Rayleigh fading channel, and two proposed receivers are mathematically modelled in discrete-time domain. The authors evaluated the bit error rate performance of the receivers by means of both theoretical analysis and numerical simulation. The performance comparison shows that the two proposed receivers can perform well under the studied channel, where the performances get better when the number of paths increases and the DSS receiver outperforms the SDS one.

FE simulation and experimental tests of high-strength structural bolts under tension

Experimental tests have been completed for high-strength 8.8 bolts for studying their mechanical performance subjected to tensile loading. As observed from these tests, failure of structural bolts has been identified as in one of two ways: threads stripping and necking of the threaded portion of the bolt shank, which is possibly due to the degree of fit between internal and external threads. Following the experimental work, a numerical approach has been developed for demonstration of the tensile performance with proper consideration of tolerance class between bolts and nuts. The degree of fit between internal and external threads has been identified as a critical factor affecting failure mechanisms of high-strength structural bolts in tension, which is caused by the machining process. In addition, different constitutive material laws have been taken into account in the numerical simulation, demonstrating the entire failure mechanism for structural bolts with different tolerance classes in their threads. It is also observed that the bolt capacities are closely associated with their failure mechanisms.

Extracting bridge frequencies from dynamic responses of two passing vehicles

A conventional way to identify bridge frequencies is utilizing vibration data measured directly from the bridge. A drawback with this approach is that the deployment and maintenance of the vibration sensors are generally costly and time-consuming. One of the solutions is in a drive-by approach utilizing vehicle vibrations while the vehicle passes over the bridge. In this approach, however, the vehicle vibration includes the effect of road surface roughness, which makes it difficult to extract the bridge modal properties. This study aims to examine subtracting signals of two trailers towed by a vehicle to reduce the effect of road surface roughness. A simplified vehicle-bridge interaction model is used in the numerical simulation; the vehicle - trailer and bridge system are modeled as a coupled model. In addition, a laboratory experiment is carried out to verify results of the simulation and examine feasibility of the damage detection by the drive-by method.

Investigation of the Impact of Dual-Lane Axle Spacing on Lateral Load Distribution, 2016

The spacing of adjacent wheel lines of dual-lane loads induces different lateral live load distributions on bridges, which cannot be determined using the current American Association of State Highway and Transportation Officials (AASHTO) Load and Resistance Factor Design (LRFD) or Load Factor Design (LFD) equations for vehicles with standard axle configurations. Current Iowa law requires dual-lane loads to meet a five-foot requirement, the adequacy of which needs to be verified. To improve the state policy and AASHTO code specifications, it is necessary to understand the actual effects of wheel-line spacing on lateral load distribution. The main objective of this research was to investigate the impact of the wheel-line spacing of dual-lane loads on the lateral load distribution on bridges. To achieve this objective, a numerical evaluation using two-dimensional linear elastic finite element (FE) models was performed. For simulation purposes, 20 prestressed-concrete bridges, 20 steel bridges, and 20 slab bridges were randomly sampled from the Iowa bridge database. Based on the FE results, the load distribution factors (LDFs) of the concrete and steel bridges and the equivalent lengths of the slab bridges were derived. To investigate the variations of LDFs, a total of 22 types of single-axle four-wheel-line dual-lane loads were taken into account with configurations consisting of combinations of various interior and exterior wheel-line spacing. The corresponding moment and shear LDFs and equivalent widths were also derived using the AASHTO equations and the adequacy of the Iowa DOT five-foot requirement was evaluated. Finally, the axle weight limits per lane for different dual-lane load types were further calculated and recommended to complement the current Iowa Department of Transportation (DOT) policy and AASHTO code specifications.

Dynamic Simulation Methods for Evaluating Motor Vehicle and Roadway Design and Resolving Policy Issues, 1990

This study had three objectives: (1) to develop a comprehensive truck simulation that executes rapidly, has a modular program construction to allow variation of vehicle characteristics, and is able to realistically predict vehicle motion and the tire-road surface interaction forces; (2) to develop a model of doweled portland cement concrete pavement that can be used to determine slab deflection and stress at predetermined nodes, and that allows for the variation of traditional thickness design factors; and (3) to implement these two models on a work station with suitable menu driven modules so that both existing and proposed pavements can be evaluated with respect to design life, given specific characteristics of the heavy vehicles that will be using the facility. This report summarizes the work that has been performed during the first year of the study. Briefly, the following has been accomplished: A two dimensional model of a typical 3-S2 tractor-trailer combination was created. A finite element structural analysis program, ANSYS, was used to model the pavement. Computer runs have been performed varying the parameters defining both vehicle and road elements. The resulting time specific displacements for each node are plotted, and the displacement basin is generated for defined vehicles. Relative damage to the pavement can then be estimated. A damage function resulting from load replications must be assumed that will be reflected by further pavement deterioration. Comparison with actual damage on Interstate 80 will eventually allow verification of these procedures.