• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1125-1134
• /
• 2017

In this paper, we present a new multilevel in space and energy diffusion (MSED) method for solving multigroup diffusion eigenvalue problems. The MSED method can be described as a PI scheme with three additional features: (1) a grey (one-group) diffusion equation used to efficiently converge the fission source and eigenvalue, (2) a space-dependent Wielandt shift technique used to reduce the number of PIs required, and (3) a multigrid-in-space linear solver for the linear solves required by each PI step. In MSED, the convergence of the solution of the multigroup diffusion eigenvalue problem is accelerated by performing work on lower-order equations with only one group and/or coarser spatial grids. Results from several Fourier analyses and a one-dimensional test code are provided to verify the efficiency of the MSED method and to justify the incorporation of the grey diffusion equation and the multigrid linear solver. These results highlight the potential efficiency of the MSED method as a solver for multidimensional multigroup diffusion eigenvalue problems, and they serve as a proof of principle for future work. Our ultimate goal is to implement the MSED method as an efficient solver for the two-dimensional/three-dimensional coarse mesh finite difference diffusion system in the Michigan parallel characteristics transport code. The work in this paper represents a necessary step towards that goal.

• Investigative Magnetic Resonance Imaging
• /
• v.11 no.2
• /
• pp.110-118
• /
• 2007

Purpose: The objective of this work to construct eigenvalue maps that have information of magnitude of three primary diffusion directions using diffusion tensor images. Materials and Methods: To construct eigenvalue maps, we used a 3.0T MRI scanner. We also compared the Moore-Penrose pseudo-inverse matrix method and the SVD (single value decomposition) method to calculate magnitude of three primary diffusion directions. Eigenvalue maps were constructed by calculating of magnitude of three primary diffusion directions. We did investigate the relationship between eigenvalue maps and fractional anisotropy map. Results: Using Diffusion Tensor Images by diffusion tensor imaging sequence, we did construct eigenvalue maps of three primary diffusion directions. Comparison between eigenvalue maps and Fractional Anisotropy map shows what is difference of Fractional Anisotropy value in brain anatomy. Furthermore, through the simulation of variable eigenvalues, we confirmed changes of Fractional Anisotropy values by variable eigenvalues. And Fractional anisotropy was not determined by magnitude of each primary diffusion direction, but it was determined by combination of each primary diffusion direction. Conclusion: By construction of eigenvalue maps, we can confirm what is the reason of fractional anisotropy variation by measurement the magnitude of three primary diffusion directions on lesion of brain white matter, using eigenvalue maps and fractional anisotropy map.

• Kyungpook Mathematical Journal
• /
• v.57 no.2
• /
• pp.211-222
• /
• 2017

In this paper, we consider three inverse eigenvalue problems for a special type of acyclic matrices. The acyclic matrices considered in this paper are described by a graph called a broom on n + m vertices, which is obtained by joining m pendant edges to one of the terminal vertices of a path on n vertices. The problems require the reconstruction of such a matrix from given partial eigen data. The eigen data for the first problem consists of the largest eigenvalue of each of the leading principal submatrices of the required matrix, while for the second problem it consists of an eigenvalue of each of its trailing principal submatrices. The third problem has an eigenvalue and a corresponding eigenvector of the required matrix as the eigen data. The method of solution involves the use of recurrence relations among the leading/trailing principal minors of ${\lambda}I-A$, where A is the required matrix. We derive the necessary and sufficient conditions for the solutions of these problems. The constructive nature of the proofs also provides the algorithms for computing the required entries of the matrix. We also provide some numerical examples to show the applicability of our results.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1990.10a
• /
• pp.22-27
• /
• 1990

Since an eigenvalue problem in structural analysis has been recognized as an important process for the assessment of structural strength, it is usually to be carried out the eigenvalue analysis or buckling analysis of structures when the compression behabiour of the member is dorminant. In general, various variables involved in the eigenvalue problem have also shown their variability. So it is natural to apply the probabilistic analysis into such problem. Since the limit state equation for the eigenvalue analysis or buckling reliability analysis is expressed implicitly in terms of random variables involved, the probabilistic finite element method is combined with the conventional reliability method such as MVFOSM and AFOSM for the determination of probability of failure due to buckling. The accuracy of the results obtained by this method is compared with results from the Monte Carlo simulations. Importance sampling method is specially chosen for overcomming the difficulty in a large simulation number needed for appropriate accurate result. From the results of the case study, it is found that the method developed here has shown good performance for the calculation of probability of buckling failure and could be used for checking the safety of the calculation of probability of buckling failure and could be used for checking the safely of frame structure which might be collapsed by either yielding or buckling.

• Structural Engineering and Mechanics
• /
• v.52 no.5
• /
• pp.1051-1067
• /
• 2014

The structures of concern in this study are subject to two types of forces: dead loads from the acceleration imposed on the structures as well as the installed operation machines and the additional adjustable forces. We wish to determine the critical values of the adjustable forces when buckling of the structures occurs. The mathematical statement of such a problem gives rise to a constrained eigenvalue problem (CEVP) in which the dominant eigenvalue is subject to an equality constraint. A numerical algorithm for solving the CEVP is proposed in which an iterative method is employed to identify an interval embracing the target eigenvalue. The algorithm is applied to four engineering application examples finding the critical loads of a fixed-free beam subject to its own body force, two plane structures and one wide-flange beam using shell elements when acceleration force is present. The accuracy is demonstrated using the first example whose classical solution exists. The significance of the equality constraint in the EVP is shown by comparing the solutions without the constraint on the eigenvalue. Effectiveness and accuracy of the numerical algorithm are presented.

• Journal of Mechanical Science and Technology
• /
• v.17 no.3
• /
• pp.370-379
• /
• 2003

A study of free vibration of rectangular Mindlin plates is presented. The analysis is based on the Chebyshev pseudospectral method, which uses test functions that satisfy the boundary conditions as basis functions. The result shows that rapid convergence and accuracy as well as the conceptual simplicity are achieved when the pseudospectral method is applied to the solution of eigenvalue problems. Numerical examples of rectangular Mindlin plates with clamped and simply supported boundary conditions are provided for various aspect ratios and thickness-to length ratios.

• Water Engineering Research
• /
• v.5 no.4
• /
• pp.177-183
• /
• 2004

This paper introduces an eigenvalue method of transforming the hyperbolic partial differential equations of a particular unsteady friction water hammer model into characteristic form. This method is based on the solution of the corresponding one-dimensional Riemann problem that transforms hyperbolic quasi-linear equations into ordinary differential equations along the characteristic directions, which in this case arises as the eigenvalues of the system. A mathematical justification and generalization of the eigenvalues method is provided and this approach is compared to the traditional characteristic method.

• Communications of the Korean Mathematical Society
• /
• v.26 no.4
• /
• pp.591-601
• /
• 2011

In this paper, we deal with the discrete p-Laplacian operators with a potential term having the smallest nonnegative eigenvalue. Such operators are classified as its smallest eigenvalue is positive or zero. We discuss differences between them such as an existence of solutions of p-Laplacian equations on networks and properties of the energy functional. Also, we give some examples of Poisson equations which suggest a difference between linear types and nonlinear types. Finally, we study characteristics of the set of a potential those involving operator has the smallest positive eigenvalue.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.5-16
• /
• 1999

We introduce a new projector for accelerating convergence of a symmetric eigenvalue problem Ax = x, and devise a power/Lanczos hybrid algorithm. Acceleration can be achieved by removing the hard-to-annihilate nonsolution eigencomponents corresponding to the widespread eigenvalues with modulus close to 1, by estimating them accurately using the Lanczos method. However, the additional Lanczos results can be obtained without expensive matrix-vector multiplications but a very small amount of extra work, by utilizing simple power-Lanczos interconversion algorithms suggested. Numerical experiments are given at the end.

• Journal of Mechanical Science and Technology
• /
• v.17 no.10
• /
• pp.1458-1465
• /
• 2003

An eigenvalue analysis of the circular Mindlin plates with free boundary conditions is presented. The analysis is based on the Chebyshev-Fourier pseudospectral method. Even though the eigenvalues of lower vibration modes tend to convergence more slowly than those of higher vibration modes, the eigenvalues converge for sufficiently fine pseudospectral grid resolutions. The eigenvalues of the axisymmetric modes are computed separately. Numerical results are provided for different grid resolutions and for different thickness-to-radius ratios.