Sharif Digital Repository

Simulation and controlling nonlinear process is one of the challenging problems in academic and industrial research. Saving data for keeping history of process or arising number of calculation for optimization in control design, is forcing to decrease number of equation or order of system. One of the methods to decrease the order of distillation process is using the travelling wave theory. This solution shows god accuracy for simulation the behavior of distillation column. On the other hand finite difference method as a solver for distillation equation takes some time to converge to steady state condition and the numbers of calculation in this method are high. Finite difference method is... 

A finite difference technique is developed to predict the second stage creep displacement rates and stress analysis of a short fiber metal matrix composite subjecting to a constant axial load. The exponential law is adopted to describe the matrix creep behavior. Also, a method for prediction of interfacial debonding at fiber/matrix interface is developed using a stress based method. The obtained results could greatly help to better understand the flow pattern of matrix material and the load transfer mechanism between fiber and matrix. The stress components and strain rates are also validated by the available FEM and experimental results  

In this thesis, localization of a postulated noise from limited neutron detectors sparsely distributed throughout the core of a typical VVER-1000 reactor is investigated. For this purpose, developing a 2-D neutron noise simulator for hexagonal geometries based on the 2-group diffusion approximation, the reactor dynamic transfer function is calculated. The box-scheme finite difference method is first developed for hexagonal geometries, to be used for spatial discretisation of both 2-D 2-group static and noise diffusion equations. Using the discretised static equations, a 2-D 2-group static simulator (HEXDIF-2) is developed which its results are benchmarked against the well-known CITATION... 

It is the goal of this thesis to develop an analytical model for the hydro-mechanical deep drawing (HDD) process of an axisymmetric sheet metal with the fixed gap method to evaluate the effects of some assumptions such as: proportional loading, plane stress, and constant thickness conditions. The effect of parameters on the HDD process is also studied. The main model is developed with considering the normal stress and part thickness change, non-proportional loading, bending and unbending effects at the top of the cup wall. The interrelationships between geometrical and mechanical variables are obtained in the finite difference form based on the incremental strain theory, thereby being solved... 

The main goal of this project is analysis, design and implementation of scanning near field optical system for detection of biological species. The activities fall in two main category. Theoretical and experimental. In theoretical part, after studying different models describing near field interaction, we have developed software for computationally analysis of nonlinear interaction of light with nanostructures, considering third order nonlinear susceptibility and dispersion behavior of permittivity for metallic nanostructures. The software implements three dimensional finite difference time domain (FDTD) method for analysis of interaction of electromagnetic wave with matter. In developed... 

In this study, fluid-solid interaction (FSI) is simulated computationally by using a high-order accurate numerical method. The two-dimensional incompressible viscous flows are considered in the fluid domain. The primary problem with solutions of the incompressible Navier–Stokes equations is the difficulty of coupling changes in the velocity field with changes in the pressure field while satisfying the continuity equation. Herein, the artificial compressibility method is used to overcome this difficulty. Preconditioning is implemented to reduce the stiffness of the system of equations to increase the convergence rate of the solution. Using preconditioning, physical solutions even at low... 

This thesis is about modeling of liquid-vapor two phase flows through nozzle. Liquid-vapor two phase flow is very applicable in industries such as boilers, expansion valves, refrigeration and sudden failures in pipelines. Two models are used for modeling of fluid flow through nozzle, transient and steady state. In Both model the solution field is considered as saturated liquid and vapor and the mass, momentum and energy equations as well as equation of state is used for describing of fluid flow properties. Homogeneous equilibrium model is used and for complete modeling of fluid flow heat transfer and friction force is also considered.
Numerical solutions are used for solving of the... 

Cancer invasion of tissue is a key aspect of the growth and spread of cancer and is crucial in the process of metastatic spread i.e. the growth of secondary cancers. Invasion consists in cancer cells secreting various matrix degrading enzymes (MDEs) which destroy the surronding tissue or extracellular matrix (ECM). Through a combination of proliferation and migration, the cancer cells then actively spread locally into the surrounding tissue. Thus processes occuring at the level of individual cells eventually give rise to processes occuring at the tissue level. In this thesis we introduce a new type of multiscale model describing the process of cancer invasion of tissue.Our multiscale model... 

In the present study, the numerical simulation of incompressible laminar and turbulent flows using a high-order finite difference lattice Boltzmann method is presented. To handle curved geometries with non uniform grids, the incompressible form of lattice Boltzmann equation is transformed into the generalized curvilinear coordinates and the spatial derivatives of the resulting equation are discretized using the fifth-order WENO scheme. The advantage of using the WENO-LBM developed is that it needs less number of grid points and remains stable even at high Reynolds number flows. For the temporal term, the fourth-order explicit Rung-Kutta scheme is adopted for laminar flow calculations and... 

With the transportation progress and appearance of transit industry, necessity of modern vehicles redounded to appearance of heavier trucks with higher speeds, moving on the roads. In addition, with the progress in mechanical engineering and automobile industry, appearance of such trucks is growing increasingly. So, bridges as one of the structures that civil engineers design and construct, nowadays are subjected to higher moving dynamic loads in comparison with the past. As a result, lots of investigation and researches in the universities and institutions all over the world are being conducting on the effects of the vehicle speeds on the dynamic stresses of the bridges. Bridge codes... 

Regarding the issue that significant amount of energy consumption in the world is dedicated to heating and cooling of the buildings, by using traditional methods of heating and cooling, the environment is facing serious problems like green house gases. There were various techniques for decreasing the amount of contaminants stem from this process. Heat-exchanger energy piles are one of the most common methods that will result in economic usage of energy resources. Assessing the long-term behavior of the energy piles requires comprehensive understanding of the complex interaction between soil and pile subjected to mechanical and thermal loadings. Several numerical and analytical methods have... 

Because of seismic flexibility, aesthetics, economic advantage and rapid construction of reinforced soil walls, application of these walls has raised tremendously since 1970’s. The early reinforced soil retaining walls used metallic srips to reinforce the soil with precast concrete panels as facing, however, since 1990’s using of geosynthetics to reinforce these walls with modular concrete block facing have replaced the old fashion. In this study, a finite difference code was developed in FLAC software. After verification of the numerical modeling procedure, a wall with 3.6m height reinforced by geosynthetics was simulated. In the course of numerical simulations the stage of counstruction,... 

Liquefaction is an important phenomenon in geotechnical engineering which can cause severe damages to structures. Liquefaction-induced lateral spreading is defined as the lateral displacement in mild slopes or level grounds ending in free faces (such as quay walls) triggered by liquefaction in subsurface soil layers. During recent years, extensive studies have been conducted around the world documenting liquefaction induced lateral spreading and its effects on deep foundations. In the present study, a series of shaking table experiments which were previously conducted at Sharif University of Technology are numerically simulated using the three dimensional finite difference based program,... 

In the present study, the numerical simulation of the panel flutter in compressible inviscid flow is performed by the compact finite difference method. For this purpose, the 2D compressible Euler equations written in the arbitrary Lagrange-Eulerian form are considered and the resulting system of equations in the generalized curvilinear coordinates is solved by the fourth-order compact finite-difference method. An appropriate nonlinear filter is applied for the shock capturing and for the solution to be stable. The governing equation for the panel is also numerically solved by using the fourth-order compact finite difference method. The time integration in the flow domain is made by the... 

The numerical solution of the parabolized Navier-Stokes (PNS) equations for accurate computation of hypersonic axisymmetric flowfield with real gas effects is obtained by using the fourth-order compact finite-difference method. The PNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock fitting procedure is utilized in the compact PNS scheme to obtain accurate solutions in the vicinity of the shock. To stabilize the numerical solution, numerical dissipation term and filters are used. The main advantage of the present formulation is that the basic flow variables... 

In the present study, the preconditioned incompressible Navier‐Stokes equations with the artificial compressibility (AC) method formulated in the generalized curvilinear coordinates are numerically solved by using a high‐order compact finite‐difference scheme for accurately and efficiently computing the incompressible flows. A fourth‐order compact finite‐difference scheme is utilized to discretize the spatial derivative terms of the resulting system of equations and the time integration is carried out based on the dual time‐stepping method. The capability of the proposed solution methodology for computing the steady and unsteady incompressible viscous flows in a wide range of Reynolds... 

In this research, a simple yet efficient two-dimensional model of proton exchange membrane fuel cell (PEMFC) is developed using stream function-vorticity formulation. The model accounts for the fluid flow and masstransport processes in cathode gas channel and gas diffuser layer. The governing equations are decoupled and solved usingtwo efficient approaches, Alternating Direction Implicit (ADI) and False Transient Method (FTM) in a modified geometry to simulate the fully-developed condition consistently. An in-house code was developed to perform a parametric study of fuel cell performance. The results show the variations of important parameters along the cathode channel such as velocity... 

Innovations in laser, microwaves, and similar technologies have significantly advanced thermal treatments for diseases or even injuries concerning skin tissue. For a thorough understanding in the underlying mechanisms of bioheat transfer behavior of skin,a1D unsteady non-dimensional hyperbolic model of heat transfer through this tissue with metabolic heat generation which is subject to specific boundary conditions, is solved numerically using the finite difference method. A thermal shock is generated at the base of the tissue, which moves forward with a finite speed. A Numerical solution for a simple one-layer skin tissue is obtained. Then, the effects of various parameters, time step,... 

In this work, a high-order accurate gas kinetic scheme based on the compact finite-difference Boltzmann method (CFDBM) is developed and applied for simulating the compressible rarefied gas flows. Here, the Shakhov model of the Boltzmann equation is considered and the spatial derivative term in the resulting equation is discretized by using the fourth-order compact finite-difference method and the time integration is performed by using the third-order TVD Runge-Kutta method. A filtering procedure with three discontinuity-detecting sensors is applied and examined for the stabilization of the solution method especially for the problems involving the discontinuity regions such as the shock. The... 

Neutron noise equations, which are obtained by assuming small perturbations of macroscopic cross sections around a steady-state neutron field and by subsequently taking the Fourier transform in the frequency domain, have been usually solved by analytical techniques or by resorting to diffusion theory, but in this thesis, in order to increase of accuracy of neutron noise calculation, has been used transport approximation for neutron noise calculation and the Monte Carlo method has been used to solve transport equation of the neutron noise in the frequency domain. Since the transport equation of the neutron noise is a complex equation, a new Monte Carlo technique for treating complex-valued...