Sharif Digital Repository

In this paper, two new approximation functions are introduced. These new techniques, which are referred herein as the multi-triangles method (MTM) and weighted multi-triangles method (WMTM) are applied for the approximation of unknowns and their derivatives at the points of interest. The approximations are performed in terms of the unknowns corresponding to the field nodes which are the vertices of the region surrounding the desired point and determined by Delaunay triangulations. The capability and accuracy of the proposed approximation functions are compared with the other approximating techniques in the meshless finite volume (MFV) frame work for some benchmark problems. Numerical... 

In pebble bed reactors, the core is filled with thousands of graphite and fuel pebbles. Fuel pebbles in these reactors consist of TRISO particles, which are embedded in a graphite matrix stochastically. The reactor core is also stochastically filled with pebbles. These two stochastic geometries comprise the so-called double heterogeneous nature of this type of reactor. In this paper, a pebble bed reactor, the HTR-10, is used to demonstrate a treatment of this double heterogeneity using the MCNP5 Monte Carlo code and MATLAB programming. In this technique, TRISO particles are modeled in a pebble using the expanded FILL and LATTICE features of MCNP5. MATLAB is used to generate random numbers... 

An improvement to the Direct Discrete Method (DDM), also known as the Cell Method, has been discussed. The improvement is based on a duality theorem between the primal and dual complexes. Also, the analog counterpart of the Integral operator has been derived in this paper. The multi-group neutron diffusion is then derived, directly in a discrete algebraic form, according to this procedure. A numerical example has shown that this method would yield a high order of convergence (approximately 4.6) if its parameters are adjusted suitably. Finally, the method is applied to the 2D IAEA benchmark problem, and has shown to yield accurate solutions with a reasonably low number of unknowns  

Creating hydraulically induced fractures in oil/gas reservoirs is one of the methods for Enhanced Oil Recovery (EOR) that has been applied extensively in petroleum industry in recent years. Despite its popularity, the design process of Hydraulic Fracture treatment is mostly empirical based on the previous experiences gained in the oil-rich formation. The reason lies in the complexities involved in the Hydraulic Fracture process including interacting effects of fluid(s) flowand solid deformations, injection of non-Newtonian fluids in the porous media, leak-off of the injected fluid into the formation, complex geometry of the induced fracture in the intact or naturally fractured rock,... 

A robust multiple model adaptive control strategy using fuzzy fusion (RMMAC-FF) is presented in this paper. The main idea in multi-model controllers is to identify the best model of the system at any instant of time and apply the appropriate control input to it. RMMAC-FF, integrates a fuzzy robust controller, with the fuzzy multiple model adaptive estimation and a fuzzy switching to come up with a new strong methodology to control complex systems. Simulation results of the RMMAC-FF on a two-cart system, used as a benchmark problem, verify the theory and confirm the effectiveness of the proposed controller  

In this paper, a new interaction prediction approach for hierarchical control of non-linear large-scale systems is presented. The proposed approach uses a new gradient-type coordination scheme which is robust with respect to the parameters' variation, and also has a good convergence rate. In classical coordination strategies, which can be divided into the gradient-type and substitution-type approaches, it is not possible to improve the robustness and the convergence rate at the same time, since by increasing one the other decreases. The proposed approach has the main advantages of the gradient-type algorithms in being independent of the parameter's variation and also the initial guess of the... 

In this paper, two adaptive refinement techniques are presented for enhancing the capability of the new kind of finite volume method called meshless finite volume (MFV) method in which the moving least squares (MLS) approximation technique is employed. The proposed algorithms perform by inserting new nodes in large error regions identified using the Zienkiewicz-Zhu (Z-Z) error estimator and the T-Belytschko (TB) stress recovery scheme. In the first method referred as CV-based adaptive refinement, the new nodes are inserted at the vertices of control volumes with higher errors. The second method, referred as GA-based adaptive refinement, contains two schemes where an adaptive refinement... 

An improved incompressible smoothed particle hydrodynamics (ISPH) method is presented, which employs first-order consistent discretization schemes both for the first-order and second-order spatial derivatives. A recently introduced wall boundary condition is implemented in the context of ISPH method, which does not rely on using dummy particles and, as a result, can be applied more efficiently and with less computational complexity. To assess the accuracy and computational efficiency of this improved ISPH method, a number of two-dimensional incompressible laminar internal flow benchmark problems are solved and the results are compared with available analytical solutions and numerical data.... 

Hierarchical reinforcement learning facilitates learning in large and complex domains by exploiting subtasks and creating hierarchical structures using these subtasks. Subtasks are usually defined through finding subgoals of the problem. Providing mechanisms for autonomous subgoal discovery and skill acquisition is a challenging issue in reinforcement learning. Among the proposed algorithms, a few of them are successful both in performance and also efficiency in terms of the running time of the algorithm. In this paper, we study four methods for subgoal discovery which are based on graph partitioning. The idea behind the methods proposed in this paper is that if we partition the transition... 

In the present study, the development of a neutron noise simulator, DYN-ACNEM, using the Average Current Nodal Expansion Method (ACNEM) in 2-G, 2-D hexagonal geometries is reported. In first stage, the static neutron calculation is performed. The neutron/adjoint flux distribution and corresponding eigen-values are calculated using the algorithm developed based on power iteration method by considering the coarse meshes. The results of the static calculation are validated against the well-known IAEA-2D benchmark problem. In the second stage, the dynamic calculation is performed in the frequency domain in which the dimension of the variable space of the noise equations is lower than the time... 

In this study, a jet propulsion JP combustor is studied numerically to investigate the combustor wall temperature influences on the soot characteristics emitted at its exhaust. There are a number of ways to control the combustor wall temperatures benefiting from different wall cooling technologies. Irrespective of using different high technology cooling systems, it is important to recognize how the wall temperature can affect the soot emission from one specific JP engine. Before examining the main combustor, it is important to assess the accuracy of the computational fluids dynamics (CFD) tool via solving a benchmark problem. In this regard, the predicted flame structure for the benchmark... 

Synthesis of the optimum conventional (with non-sharp separations) distillation column sequence (DCS) is a challenging problem, in the field of chemical process design and optimization, due to its huge search space and combinatorial nature. In this paper, a novel procedure for the synthesis of optimum Conventional Distillation Column Sequence is proposed. The proposed method is based on evolutionary algorithms. The main criterion used to screen alternative DCS's is the Total Annual Cost (TAC). In order to estimate the TAC of each DCS alternative all columns that exist in the DCS are designed using short-cut methods. The performance of the proposed method and other alternatives are compared... 

In this paper, a new interaction prediction approach is presented for optimal control of nonlinear large-scale systems. The proposed approach uses a new gradient-type coordination scheme which has a larger convergence region with respect to the parameters' variation, and also has a good convergence rate. In this approach, the coordination vector is updated using the gradient of coordination error. This type of coordination considerably reduces the number of iterations. The robustness and the convergence rate of the proposed approach against the best classical interaction prediction approaches are shown through simulations of a benchmark problem. © 2008 IEEE  

This paper numerically studies the interaction of a flowing granular material with an entrainable granular bed, while materials are mixed at the interface of two materials. The rheological behavior of this granular mixture is characterized by a generalized viscoplastic model that includes local volume fraction of materials as well as their physical properties, i.e. size, density, and friction angle. Additionally, the effect of the dynamics of entrained bed-type particles on the rheology of the granular mixture is considered. The governing equations of the flow are discretized using the Incompressible Smoothed Particle Hydrodynamics (SPH) method in which mixing of particles can be... 

Various methods for solving the forward/adjoint equation in hexagonal and rectangular geometries are known in the literatures. In this paper, the solution of multigroup forward/adjoint equation using Finite Element Method (FEM) for hexagonal and rectangular reactor cores is reported. The spatial discretization of equations is based on Galerkin FEM (GFEM) using unstructured triangle elements. Calculations are performed for both linear and quadratic approximations of the shape function; based on which results are compared. Using power iteration method for the forward and adjoint calculations, the forward and adjoint fluxes with the corresponding eigenvalues are obtained. The results are then... 

In the present study, the neutron noise, i.e. The stationary fluctuation of the neutron flux around its mean value, is calculated in 2-group forward and adjoint diffusion theory for both hexagonal and rectangular reactor cores. To this end, the static neutron calculation is performed at the first stage. The spatial discretization of equations is based on linear approximation of Galerkin Finite Element Method (GFEM) using unstructured triangle elements. Using power iteration method, forward and adjoint fluxes with the corresponding eigenvalues are obtained. The results are then benchmarked against the valid results for BIBLIS-2D and IAEA-2D benchmark problems and DONJON computer code. The... 

A modified weakly compressible smoothed particle hydrodynamics (WCSPH) is presented, which utilizes consistent discretization schemes for spatial derivatives in the flow equations. Here, each SPH particle is considered as a computational point that represents a specific part of the fluid. To overcome non-physical oscillations that usually arise in standard WCSPH, we modified the mass conservation equation by using a numerical filter. This modification is based on the difference between two discretization schemes used for the term ∇{dot operator}∇Pρ. Furthermore, a new implementation of wall boundary condition in SPH is introduced. This condition is imposed on the pressure of wall boundary... 

A new hybrid approach for dynamic optimization problems with continuous search spaces is presented. The proposed approach hybridizes efficient features of the particle swarm optimization in tracking dynamic changes with a new evolutionary procedure. In the proposed dynamic hybrid PSO (DHPSO) algorithm, the swarm size is varied in a self-regulatory manner. Inspired from the microbial life, the particles can reproduce infants and the old ones die. The infants are especially reproduced by high potential particles and located near the local optimum points, using the quadratic interpolation method. The algorithm is adapted to perform in continuous search spaces, utilizing continuous movement of... 

We present an optimal solution procedure for minimizing total weighted resource tardiness penalty costs in the resource-constrained project scheduling problem. In this problem, we assume the constrained renewable resources are limited to very expensive equipments and machines that are used in other projects and are not available in all periods of time of a project. In other words, for each resource, there is a dictated ready date as well as a due date such that no resource can be available before its ready date but the resources are permitted to be used after their due dates by paying penalty cost depending on the resource type. We also assume that only one unit of each resource type is... 

Weakly Compressible Smoothed Particle Hydrodynamics (WCSPH) can lead to non-physical oscillations in the pressure and density fields when simulating incompressible flow problems. This in turn may result in tensile instability and sometimes divergence. In this paper, it is shown that this difficulty originates from the specific form of spatial discretization used for the pressure term when solving the mass conservation equation. After describing the pressure-velocity decoupling problem associated with the so-called colocated grid methods, a modified approach is presented that overcomes this problem using a different discretization scheme for the second derivative of pressure. The modified...