A Revisit to CMFD Schemes: Fourier Analysis and Enhancement

The coarse-mesh finite difference (CMFD) scheme is a very effective nonlinear diffusion acceleration method for neutron transport calculations. CMFD can become unstable and fail to converge when the computational cell optical thickness is relatively large in k-eigenvalue problems or diffusive fixed-source problems. Some variants and fixups have been developed to enhance the stability of CMFD, including the partial current-based CMFD (pCMFD), optimally diffusive CMFD (odCMFD), and linear prolongation-based CMFD (lpCMFD). Linearized Fourier analysis has proven to be a very reliable and accurate tool to investigate the convergence rate and stability of such coupled high-order transport/low-order diffusion iterative schemes. It is shown in this paper that the use of different transport solvers in Fourier analysis may have some potential implications on the development of stabilizing techniques, which is exemplified by the odCMFD scheme. A modification to the artificial diffusion coefficients of odCMFD is proposed to improve its stability. In addition, two explicit expressions are presented to calculate local optimal successive overrelaxation (SOR) factors for lpCMFD to further enhance its acceleration performance for fixed-source problems and k-eigenvalue problems, respectively.

Three-Dimensional Cartesian Grid Method for the Simulations of Flows with Shock Waves in the Domains with Varying Boundaries

This paper is devoted to the development and quantitative evaluation of the properties of the numerical algorithm of the Cartesian grid method for three-dimensional (3D) simulation of shock-wave propagation in areas of varying shape. The detailed description of the algorithm is presented. The algorithm is relatively simple to implement and does not require solving the problem of determination of the shape of the body’s boundary intersection with regular computational cell. The accuracy of the algorithm is demonstrated by comparing the simulated and experimental data in the problems of the interaction of a shock wave (SW) with a nonmoving sphere and a moving particle.

Damage and Failure in a Statistical Crack Model

In this paper I will take a close look at a statistical crack model (SCM) as is used in engineering computer codes to simulate fracture at high strain rates. My general goal is to understand the macroscopic behavior effected by the microphysical processes incorporated into an SCM. More specifically, I will assess the importance of including local interactions between cracks into the growth laws of an SCM. My strategy will be to construct a numerical laboratory that represents a single computational cell containing a realization of a statistical distribution of cracks. The cracks will evolve by the microphysical models of the SCM, leading to quantifiable damage and failure of the computational cell. I will use the numerical data generated by randomly generated ensembles of the fracture process to establish scaling laws that will modify and simplify the implementation of the SCM into large scale engineering codes.

Computational cell cycle analysis of single cell RNA-Seq data

The variation in gene expression profiles of cells captured in different phases of the cell cycle can interfere with cell type identification and functional analysis of single cell RNA-Seq (scRNA-Seq) data. In this paper, we introduce SC1CC (SC1 Cell Cycle analysis tool), a computational approach for clustering and ordering single cell transcriptional profiles according to their progression along cell cycle phases. We also introduce a new robust metric, Gene Smoothness Score (GSS) for assessing the cell cycle based order of the cells. SC1CC is available as part of the SC1 web-based scRNA-Seq analysis pipeline, publicly accessible at https://sc1.engr.uconn.edu/.

Hydraulic Transient Analysis of Sewer Pipe Systems Using a Non-Oscillatory Two-Component Pressure Approach

On the basis of the two-component pressure approach, we developed a numerical model to capture mixed transient flows in close conduit systems. To achieve this goal, an innovative Godunov finite-volume numerical scheme is proposed to suppress the spurious numerical oscillations occurring during rapid pipe pressurization. To dissipate the spurious numerical oscillations, we admit artificial numerical viscosity to the numerical scheme through applying a proposed Harten, Lax, and van Leer (HLL) Riemann solver for calculating the numerical fluxes at the computational cell interfaces. The proposed solver controls the magnitude of the numerical viscosity through adjusting the left and right wave velocities. A wave velocity calculator is proposed to optimally distribute the numerical viscosity over several computational cells around the computational cell in which the pressurization front is located. The proposed solver admits significant artificial numerical viscosity when the pipe pressurization is imminent and automatically reduces it in other places; in this way the numerical diffusion and data smearing is minimized. The validity of the proposed model is justified by the aid of several test cases in which the numerical results are compared with both experimental data and the results obtained from analytical methods. The results reveal that the proposed model succeeds in completely removing the spurious numerical oscillations, even when the pipe acoustic speed is over 1000 m/s. The numerical results also show that the model can successfully capture occurrence of negative pressures during the course of transient flow.