Three dimensional numerical simulation of thermohaline currents of the Persian Gulf

Observational data and a three dimensional numerical model (POM) are used to investigate the Persian Gulf outflow structure and its spreading pathway into the Oman Sea. The model is based on orthogonal curvilinear coordinate system in horizontal and train following coordinate (sigma coordinate) system in vertical. In the simulation, the horizontal diffusivity coefficients are calculated form Smogorinsky diffusivity formula and the eddy vertical diffusivities are obtained from a second turbulence closure model (namely Mellor-Yamada level 2.5 model of turbulence). The modeling area includes the east of the Persian Gulf, the Oman Sea and a part of the north-east of the Indian Ocean. In the model, the horizontal grid spacing was assumed to be about 3.5 km and the number of vertical levels was set to 32. The simulations show that the mean salinity of the PG outflow does not change substantially during the year and is about 39 psu, while its temperature exhibits seasonal variations. These lead to variations in outflow density in a way that is has its maximum density in late winter (March) and its minimum in mid-summer (August). At the entrance to the Oman Sea, the PG outflow turns to the right due to Coriolis Effect and falls down on the continental slope until it gains its equilibrium depth. The highest density of the outflow during March causes it to sink more into the deeper depths in contrast to that of August which the density is the lowest one. Hence, the neutral buoyancy depths of the outflow are about 500 m and 250 m for March and August respectively. Then, the outflow spreads in its equilibrium depths in the Oman Sea in vicinity of western and southern boundaries until it approach the Ras al Hamra Cape where the water depth suddenly begins to increase. Therefore, during March, the outflow that is deeper and wider relative to August, is more affected by the steep slope topography and as a result of vortex stretching mechanism and conservation of potential vorticity it separates from the lateral boundaries and finally forms an anti-cyclonic eddy in the Oman Sea. But during August the outflow moves as before in vicinity of lateral boundaries. In addition, the interaction of the PG outflow with tide in the Strait of Hormuz leads to intermittency in outflow movement into the Oman Sea and it could be the major reason for generations of Peddy (Peddies) in the Oman Sea.

Two-dimensional numerical simulation of the channel electron in an In0.52Al0.48As/In0.53Ga0.48As HEMT

A two-dimensional quantum model based on the solution of Schrodinger and Poisson equations is first presented for In0.52Al0.48As/In0.53Ga0.47As/InP HEMT. According to the model, the two-dimensional distributions of electron density and transverse electric field in the channel of InAlAs/InGaAs HEMT are discussed.

Numerical simulation of railway track supporting system using finite-infinite and thin layer elements under impulsive loads

The present study deals with two dimensional, numerical simulation of railway track supporting system subjected to dynamic excitation force. Under plane strain condition, the coupled finite-infinite elements to represent the near and far field stress distribution and thin layer interface element was employed to model the interfacial behavior between sleepers and ballast. To account for the relative debonding, slipping and crushing that could take place in the contact area between the sleepers and ballast, modified Mohr-Coulomb criterion was adopted. Furthermore an attempt has been made to consider the elasto-plastic material non-linearity of the railway track supporting media by employing different constitutive models to represent steel, concrete and supporting materials. Based on the proposed physical and constitutive modeling a code has been developed for dynamic loads. The applicability of the developed F.E code has been demonstrated by analyzing a real railway supporting structure.

Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.

Numerical simulation of a new two-dimensional variable-order fractional percolation equation in non-homogeneous porous media

Percolation flow problems are discussed in many research fields, such as seepage hydraulics, groundwater hydraulics, groundwater dynamics and fluid dynamics in porous media. Many physical processes appear to exhibit fractional-order behavior that may vary with time, or space, or space and time. The theory of pseudodifferential operators and equations has been used to deal with this situation. In this paper we use a fractional Darcys law with variable order Riemann-Liouville fractional derivatives, this leads to a new variable-order fractional percolation equation. In this paper, a new two-dimensional variable-order fractional percolation equation is considered. A new implicit numerical method and an alternating direct method for the two-dimensional variable-order fractional model is proposed. Consistency, stability and convergence of the implicit finite difference method are established. Finally, some numerical examples are given. The numerical results demonstrate the effectiveness of the methods. This technique can be used to simulate a three-dimensional variable-order fractional percolation equation.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

Oxygen transport in a three-dimensional microvascular network incorporated with early tumour growth and preexisting vessel cooption: Numerical simulation study

We propose a dynamic mathematical model of tissue oxygen transport by a preexisting three-dimensional microvascular network which provides nutrients for an in situ cancer at the very early stage of primary microtumour growth. The expanding tumour consumes oxygen during its invasion to the surrounding tissues and cooption of host vessels. The preexisting vessel cooption, remodelling and collapse are modelled by the changes of haemodynamic conditions due to the growing tumour. A detailed computational model of oxygen transport in tumour tissue is developed by considering (a) the time-varying oxygen advection diffusion equation within the microvessel segments, (b) the oxygen flux across the vessel walls, and (c) the oxygen diffusion and consumption with in the tumour and surrounding healthy tissue. The results show the oxygen concentration distribution at different time points of early tumour growth. In addition, the influence of preexisting vessel density on the oxygen transport has been discussed. The proposed model not only provides a quantitative approach for investigating the interactions between tumour growth and oxygen delivery, but also is extendable to model other molecules or chemotherapeutic drug transport in the future study.

Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.

Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Numerical simulation of coherent structure in two-dimensional compressible mixing layers

The coherent structure in two-dimensional mixing layers is simulated numerically with the compressible Navier-Stokes equations. The Navier-Stokes equations are discretized with high-order accurate upwind compact schemes. The process of development of flow structure is presented: loss of stability, development of Kelvin-Helmholtz instability, rolling up and pairing. The time and space development of the plane mixing layer and influence of the compressibility are investigated.

Influence of liquid bridge volume on the onset of oscillation in floating-zone convection .2. Numerical-simulation

The onset of oscillation in the floating zone convection driven by the gradient of surface tension was studied numerically for an unsteady and two-dimensional model, and studies were concentrated on the influence of liquid bridge volume on the onset of oscillation in comparison with the experimental results in the Paper I. The numerical results agree with the experimental ones presented in the previous paper, in which the distributions of critical applied temperature difference depending on the volume of liquid bridge and a gap range of liquid volume in marginal stability curve were obtained.

Three-dimensional supercell simulation of novel semiconductor nanostructures

In this thesis we investigate atomic scale imperfections and fluctuations in the quantum transport properties of novel semiconductor nanostructures. For this purpose, we have developed a numerically efficient supercell model of quantum transport capable of representing potential variations in three dimensions. This flexibility allows us to examine new quantum device structures made possible through state-of-the-art semiconductor fabrication techniques such as molecular beam epitaxy and nanolithography. These structures, with characteristic dimensions on the order of a few nanometers, hold promise for much smaller, faster and more efficient devices than those in present operation, yet they are highly sensitive to structural and compositional variations such as defect impurities, interface roughness and alloy disorder. If these quantum structures are to serve as components of reliable, mass-produced devices, these issues must be addressed.

In Chapter 1 we discuss some of the important issues in resonant tunneling devices and mention some of thier applications. In Chapters 2 and 3, we describe our supercell model of quantum transport and an efficient numerical implementation. In the remaining chapters, we present applications.

In Chapter 4, we examine transport in single and double barrier tunneling structures with neutral impurities. We find that an isolated attractive impurity in a single barrier can produce a transmission resonance whose position and strength are sensitive to the location of the impurity within the barrier. Multiple impurities can lead to a complex resonance structure that fluctuates widely with impurity configuration. In addition, impurity resonances can give rise to negative differential resistance. In Chapter 5, we study interface roughness and alloy disorder in double barrier structures. We find that interface roughness and alloy disorder can shift and broaden the n = 1 transmission resonance and give rise to new resonance peaks, especially in the presence of clusters comparable in size to the electron deBroglie wavelength. In Chapter 6 we examine the effects of interface roughness and impurities on transmission in a quantum dot electron waveguide. We find that variation in the configuration and stoichiometry of the interface roughness leads to substantial fluctuations in the transmission properties. These fluctuations are reduced by an attractive impurity placed near the center of the dot.