Sharif Digital Repository

Obtaining the set of algebraic equations that directly correspond to a physical phenomenon has been viable in the recent direct discrete method (DDM). Although this method may find its roots in physical and geometrical considerations, there are still some degrees of freedom that one may suspect optimize-able. Here we have used the information embedded in the corresponding adjoint equation to form a local functional, which in turn by its minimization, yield suitable dual mesh positioning  

In this article we introduce the cylindrical construction, as an edge-replacement procedure admitting twists on both ends of the hyperedges, generalizing the concepts of lifts and Pultr templates at the same time. We prove a tensor-hom duality for this construction and we show that not only a large number of well-known graph constructions are cylindrical but also the construction and its dual give rise to some new graph constructions, applications and results. To show the applicability of the main duality we introduce generalized Grötzsch, generalized Petersen-like and Coxeter-like graphs and we prove some coloring properties of these graphs  

A morphing-body optimization method is introduced to accelerate adjoint-based shape optimization techniques. The optimization process solves the flow and adjoint equations around a continuously deforming body whose shape is controlled by the cost function. Effect of various parameters on the efficiency of the scheme is studied. It is found that, for the best performance of the algorithm, the morphing rate of the airfoil should be restricted, since larger rates foster oscillations and lower values are not computationally feasible. Moreover, the iterative procedure in the adjoint solver should be adapted to the iteration scheme in the flow solver and to the morphing rate  

To quantify the core barrel motion (CBM) in a pressurized water reactor, a scaling factor can be calculated to convert the Root Mean Square (RMS) value of the ex-core signals (%) to the core barrel motion amplitude (mil) (Thompson et al., 1980). In the current paper, a scaling factor is calculated using the direct and adjoint methods for a typical VVER-1000 reactor. The scaling factor is calculated using the perturbed parameters that result from CBM. The results show that the calculated scaling factors are not the same in one and two-dimensional modeling, and strongly depend on the ex-core detector location. The linearity assumption of relative detector response versus the small displacement... 

To solve neutron transport equation in multigroup approach, in addition to weighting function and number of energy groups, proper selection of the group boundaries have high importance for the accuracy of the calculations. In the current paper, the bilinear combination of forward and adjoint neutron spectra is used for the optimization of 69 energy group structure of WIMSD5 lattice physics code. To remedy the energy self-shielding effect, homogeneous adjoint and forward BN equations on an ultrafine energy group structure have been solved to obtain the ultrafine forward and adjoint spectra. The coarse group intervals are selected to have equal values of bilinear function in each... 

Accurate calculation of kinetic parameters is of utmost importance in the safety analysis of a nuclear reactor. In the current paper, two approaches are investigated to evaluate these parameters in energy phase space. In the first approach, these parameters are derived from an energy-continuous form of the forward and adjoint transport equations and then integrals with respect to the energy variable are replaced by weighted summations over the energy groups, while in the second approach these parameters are extracted from the multi-group forward equation and its associate adjoint equation in which their multigroup constants are weighted by forward spectrum. The difference of weighting... 

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... 

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... 

Abstract In this paper, the solution of multi-group neutron/adjoint equation using Finite Element Method (FEM) for hexagonal and rectangular reactor cores is reported. The spatial discretization of the neutron diffusion equation is performed based on two different Finite Element Methods (FEMs) using unstructured triangular elements generated by Gambit software. Calculations are performed using Galerkin and Generalized Least Squares FEMs; based on which results are compared. Using the power iteration method for the neutron and adjoint calculations, the neutron and adjoint flux distributions with the corresponding eigenvalues are obtained. The results are then validated against the valid... 

In the present study, the neutron noise, i.e. the stationary fluctuation of the neutron flux around its mean value is calculated based on the 2G, 3D neutron diffusion theory. To this end, the static neutron calculation is performed at the first stage. The spatial discretization of the neutron diffusion equation is performed based on linear approximation of Galerkin Finite Element Method (GFEM) using unstructured tetrahedron elements. Using power iteration method, neutron flux and corresponding eigen-value are obtained. The results are then benchmarked against the valid results for VVER-1000 (3D) benchmark problem. In the second stage, the neutron noise equation is solved using GFEM and... 

In the present paper, development of the three-dimensional (3D) computational code based on Galerkin finite element method (GFEM) for solving the multigroup forward/adjoint diffusion equation in both rectangular and hexagonal geometries is reported. Linear approximation of shape functions in the GFEM with unstructured tetrahedron elements is used in the calculation. Both criticality and fixed source calculations may be performed using the developed GFEM-3D computational code. An acceptable level of accuracy at a low computational cost is the main advantage of applying the unstructured tetrahedron elements. The unstructured tetrahedron elements generated with Gambit software are used in the... 

The purpose of the present study is to develop the neutron diffusion solver and neutron noise simulator based on the Boundary Element Method (BEM). The 2-D, 2-G neutron/adjoint diffusion equation and corresponding neutron/adjoint noise equation were solved using the mentioned method. The developed neutron static and noise simulator based on the finite element method gives accurate results when the more number of the elements is used. The motivation of the present research is to use the boundary element method to reduce the computational cost. The boundary element method attempts to use the given boundary conditions to fit boundary values into the integral equation, rather than values... 

Classic identifiability analysis of flow barriers in incompressible single-phase flow reveals that it is not possible to identify the location and permeability of low-permeability barriers from production data (wellbore pressures and rates), and that only averaged reservoir properties in between wells can be identified. We extend the classic analysis by including compressibility effects. We use two approaches: a twin experiment with synthetic production data for use with a time-domain parameter-estimation technique, and a transfer-function formalism in the form of bilaterally coupled four-ports allowing for an analysis in the frequency domain. We investigate the identifiability, from noisy... 

The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS, DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time. These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation. Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters... 

In this paper, we establish the existence and uniqueness of the mild solution for stochastic Volterra equation with a non-self-adjoint operator. The specific Volterra equation that we consider is a generalization of the fractional differential equation. To obtain the mild solution for the case of multiplicative problem, the resolvent property of the linear perturbation of a sectorial operator will be considered. © 2023 World Scientific Publishing Company  

Neutron noise induced by propagating disturbances in VVER-type reactor core is addressed in this paper. The spatial discretization of the governing equations is based on the box-scheme finite difference method for triangular-z geometry. Using the derived equations, a 3-D 2-group neutron noise simulator (called TRIDYN-3) is developed for hexagonal-structured reactor core, by which the discrete form of both the forward and adjoint reactor dynamic transfer functions (in the frequency domain) can be calculated. In addition, both types of noise sources, namely point-like and traveling perturbations, can be modeled by TRIDYN-3. The results are then benchmarked in different cases. Considering the... 

In this paper, localization of a noise source 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 boxscheme 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. The dynamic state is assumed in the frequency domain which leads to discarding of the time disrcetisation. The developed 2-D 2- group neutron noise simulator calculates... 

In this paper, the development of a neutron noise simulator for hexagonal-structured reactor cores using both the forward and the adjoint methods is reported. The spatial discretisation of both 2-D 2-group static and dynamic equations is based on a developed box-scheme finite difference method for hexagonal mesh boxes. Using the power iteration method for the static calculations, the 2-group neutron flux and its adjoint with the corresponding eigenvalues are obtained by the developed static simulator. The results are then benchmarked against the well-known CITATION computer code. The dynamic calculations are performed in the frequency domain which leads to discarding of the time... 

A shock control bump (SCB) is a flow control method which uses local small deformations in a flexible wing surface to considerably reduce the strength of shock waves and the resulting wave drag in transonic flows. Most of the reported research is devoted to optimization in a single flow condition. Here, we have used a multi-point adjoint optimization scheme to optimize shape and location of the SCB. Practically, this introduces transonic airfoils equipped with the SCB which are simultaneously optimized for different off-design transonic flight conditions. Here, we use this optimization algorithm to enhance and optimize the performance of SCBs in two benchmark airfoils, i.e., RAE-2822 and...