Search

Article Highlights The use of a p-multigrid method significantly reduces the CPU timings for higher values of the spline degree. The Multigrid Reduction in Time method shows both strong and weak scalability up to 2048 cores. Iteration numbers are independent of the number of time steps, mesh width and spline degree.

Abstract Background The preconditioned conjugate gradient (PCG) method is the current method of choice for iterative solving of genetic evaluations. The relative difference between two successive iterates and the relative residual of the system of equations are usually chosen as a termination criterion for the PCG method in animal breeding. However, our initial analyses showed that these two commonly used termination criteria may report that a PCG method applied to a single-step single nucleotide polymorphism best linear unbiased prediction (ssSNPBLUP) is not converged yet, whereas the solutions are accurate enough for practical use. Therefore, the aim of this study was to propose two termination criteria that have been (partly) developed in other fields, but are new in animal breeding, and to compare their behavior to that of the two termination criteria widely used in animal breeding for the PCG method applied to ssSNPBLUP. The convergence patterns of ssSNPBLUP were also compared to the convergence patterns of single-step genomic BLUP (ssGBLUP). Results Building upon previous work, we propose two termination criteria that take the properties of the system of equations into account. These two termination criteria are directly related to the relative error of the iterates with respect to the true solutions. Based on pig and dairy cattle datasets, we show that the preconditioned coefficient matrices of ssSNPBLUP and ssGBLUP have similar properties when using a second-level preconditioner for ssSNPBLUP. Therefore, the PCG method applied to ssSNPBLUP and ssGBLUP converged similarly based on the relative error of the iterates with respect to the true solutions. This similar convergence behavior between ssSNPBLUP and ssGBLUP was observed for both proposed termination criteria. This was, however, not the case for the termination criterion defined as the relative residual when applied to the dairy cattle evaluations. Conclusion Our results showed that the PCG method can converge similarly when applied to ssSNPBLUP and to ssGBLUP. The two proposed termination criteria always depicted these similar convergence behaviors, and we recommend them for comparing convergence properties of different models and for routine evaluations.

The accuracy, stability and computational efficiency of numerical methods on central processing units (CPUs) for the depth-averaged shallow water equations were well covered in the literature. A large number of these methods were already developed and compared. However, on graphics processing units (GPUs), such comparisons are relatively scarce. In this paper, we present the results of comparing two time-integration methods for the shallow water equations on structured grids. An explicit and a semi-implicit time integration method were considered. For the semi-implicit method, the performance of several iterative solvers was compared. The implementation of the semi-implicit method on a GPU in this study was a novel approach for the shallow water equations. This also holds for the repeated red black (RRB) solver that was found to be very efficient on a GPU. Additionally, the results of both methods were compared with several CPU-based software systems for the shallow water flows on structured grids. On a GPU, the simulations were 25 to 75 times faster than on a CPU. Theory predicts an explicit method to be best suited for a GPU due to the higher level of inherent parallelism. It was found that both the explicit and the semi-implicit methods ran efficiently on a GPU. For very shallow applications, the explicit method was preferred because the stability condition on the time step was not very restrictive. However, for deep water applications, we expect the semi-implicit method to be preferred.

Thermal nitric-oxide (NOx) formation in industrial furnaces due to local overheating is a widely known problem. Various industries made significant investments to reduce thermal NOx by varying the operating conditions and designs of the furnace. It is difficult to find the optimal operating conditions that minimize NOx formation in the furnace by trial and error methods. The high temperature in the furnace complicates performing experiments in the furnace. Numerical modeling can provide significant information in such cases. Therefore, the objective of this paper is to obtain a numerical model of the furnace in such a way that the operating conditions can be varied and examined.In this paper, a three-dimensional steady-state finite element model for the anode baking industrial furnace is discussed. The COMSOL Multiphysics software is used for modeling the non-premixed turbulent combustion and the conjugate heat transfer to the insulation lining. The cfMesh software is used for obtaining the mesh. The results show that the simulated temperature agrees well with the measured data from our industrial partner in regions distant from the flames. The analysis shows that by decreasing the fuel mass flow rate and increasing the fuel pipe diameter by 45%, the peak in thermal NOx ppm generated in the furnace decreases by 42%. The model is limited by the use of a single-step chemistry mechanism with an eddy dissipation combustion model and a simplified approach for radiation, such as the P1 approximation model. The model can be further improved by considering a detailed chemistry mechanism model for combustion and a discrete ordinate model for radiation.

Abstract Background The preconditioned conjugate gradient (PCG) method is an iterative solver of linear equations systems commonly used in animal breeding. However, the PCG method has been shown to encounter convergence issues when applied to single-step single nucleotide polymorphism BLUP (ssSNPBLUP) models. Recently, we proposed a deflated PCG (DPCG) method for solving ssSNPBLUP efficiently. The DPCG method introduces a second-level preconditioner that annihilates the effect of the largest unfavourable eigenvalues of the ssSNPBLUP preconditioned coefficient matrix on the convergence of the iterative solver. While it solves the convergence issues of ssSNPBLUP, the DPCG method requires substantial additional computations, in comparison to the PCG method. Accordingly, the aim of this study was to develop a second-level preconditioner that decreases the largest eigenvalues of the ssSNPBLUP preconditioned coefficient matrix at a lower cost than the DPCG method, in addition to comparing its performance to the (D)PCG methods applied to two different ssSNPBLUP models. Results Based on the properties of the ssSNPBLUP preconditioned coefficient matrix, we proposed a second-level diagonal preconditioner that decreases the largest eigenvalues of the ssSNPBLUP preconditioned coefficient matrix under some conditions. This proposed second-level preconditioner is easy to implement in current software and does not result in additional computing costs as it can be combined with the commonly used (block-)diagonal preconditioner. Tested on two different datasets and with two different ssSNPBLUP models, the second-level diagonal preconditioner led to a decrease of the largest eigenvalues and the condition number of the preconditioned coefficient matrices. It resulted in an improvement of the convergence pattern of the iterative solver. For the largest dataset, the convergence of the PCG method with the proposed second-level diagonal preconditioner was slower than the DPCG method, but it performed better than the DPCG method in terms of total computing time. Conclusions The proposed second-level diagonal preconditioner can improve the convergence of the (D)PCG methods applied to two ssSNPBLUP models. Based on our results, the PCG method combined with the proposed second-level diagonal preconditioner seems to be more efficient than the DPCG method in solving ssSNPBLUP. However, the optimal combination of ssSNPBLUP and solver will most likely be situation-dependent.

Abstract Background The single-step single nucleotide polymorphism best linear unbiased prediction (ssSNPBLUP) method, such as single-step genomic BLUP (ssGBLUP), simultaneously analyses phenotypic, pedigree, and genomic information of genotyped and non-genotyped animals. In contrast to ssGBLUP, SNP effects are fitted explicitly as random effects in the ssSNPBLUP model. Similarly, principal components associated with the genomic information can be fitted explicitly as random effects in a single-step principal component BLUP (ssPCBLUP) model to remove noise in genomic information. Single-step genomic BLUP is solved efficiently by using the preconditioned conjugate gradient (PCG) method. Unfortunately, convergence issues have been reported when solving ssSNPBLUP by using PCG. Poor convergence may be linked with poor spectral condition numbers of the preconditioned coefficient matrices of ssSNPBLUP. These condition numbers, and thus convergence, could be improved through the deflated PCG (DPCG) method, which is a two-level PCG method for ill-conditioned linear systems. Therefore, the first aim of this study was to compare the properties of the preconditioned coefficient matrices of ssGBLUP and ssSNPBLUP, and to document convergence patterns that are obtained with the PCG method. The second aim was to implement and test the efficiency of a DPCG method for solving ssSNPBLUP and ssPCBLUP. Results For two dairy cattle datasets, the smallest eigenvalues obtained for ssSNPBLUP (ssPCBLUP) and ssGBLUP, both solved with the PCG method, were similar. However, the largest eigenvalues obtained for ssSNPBLUP and ssPCBLUP were larger than those for ssGBLUP, which resulted in larger condition numbers and in slow convergence for both systems solved by the PCG method. Different implementations of the DPCG method led to smaller condition numbers, and faster convergence for ssSNPBLUP and for ssPCBLUP, by deflating the largest unfavourable eigenvalues. Conclusions Poor convergence of ssSNPBLUP and ssPCBLUP when solved by the PCG method are related to larger eigenvalues and larger condition numbers in comparison to ssGBLUP. These convergence issues were solved with a DPCG method that annihilates the effect of the largest unfavourable eigenvalues of the preconditioned coefficient matrix of ssSNPBLUP and of ssPCBLUP on the convergence of the PCG method. It resulted in a convergence pattern, at least, similar to that of ssGBLUP.

This paper compares and assesses several numerical methods that solve the steady-state power flow problem on integrated transmission-distribution networks. The integrated network model consists of a balanced transmission and an unbalanced distribution network. It is important to analyze these integrated electrical power systems due to the changes related to the energy transition. We classified the existing integration methods as unified and splitting methods. These methods can be applied to homogeneous (complete three-phase) and hybrid (single-phase/three-phase) network models, which results in four approaches in total. These approaches were compared on their accuracy and numerical performance—CPU time and number of iterations—to demonstrate their applicability on large-scale electricity networks. Furthermore, their sensitivity towards the amount of distributed generation and the addition of multiple distribution feeders was investigated. The methods were assessed by running power flow simulations using the Newton–Raphson method on several integrated power systems up to 25,000 unknowns. The assessment showed that unified methods applied to hybrid networks performed the best on these test cases. The splitting methods are advantageous when complete network data sharing between system operators is not allowed. The use of high-performance techniques for larger test cases containing multiple distribution networks will make the difference in speed less significant.

The emissions from the industrial furnaces impact the environment. Among the various factories, those having anode baking furnaces are working on reducing the pollutant emissions. The aerodynamics in the furnace influences the emissions due to the high dependence of combustion and radiation phenomena on the mixing characteristics. Therefore, this paper aims to establish the numerical simulation results for the three-dimensional turbulent flow in a single section of an anode baking furnace with a high rate of fuel injection. The stabilized non-linear finite element approach on the Reynolds-averaged Navier-Stokes (RANS) equation is used with COMSOLMultiphysics. The turbulent viscosity ratio is highly sensitive to the mesh for the standard k-ϵ model. The requirements of the Cartesian and refined mesh near the jet development region is explained. The comparison of meshes generated by two meshing tools namely cfMesh and COMSOL Multiphysics default Mesher is carried out. The high numerical diffusion in the flow models due to the coarser mesh leads to convergence but deficit the precision in the results. This paper shows that the mesh generated by cfMesh with flow aligned refinement combined with the non-linear finite element solver in COMSOL Multiphysics proves to provide accurate results of turbulent quantities.

In this paper, we propose a fast linear power flow method using a constant impedance load model to simulate both the entire Low Voltage (LV) and Medium Voltage (MV) networks in a single simulation. Accuracy and efficiency of this linear approach are validated by comparing it with the Newton power flow algorithm and a commercial network design tool Vision on various distribution networks including real network data. Results show that our method can be as accurate as classical Nonlinear Power Flow (NPF) methods using a constant power load model and additionally, it is much faster than NPF computations. In our research, it is shown that voltage problems can be identified more efficiently when MV and LV are integrally evaluated. Moreover, Numerical Analysis (NA) techniques are applied to the Large Linear Power Flow (LLPF) problem with 27 million nonzeros in order to improve the computation time by studying the properties of the linear system. Finally, the original computation times of LLPF problems with real and complex components are reduced by 2.8 times and 5.7 times, respectively.