Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn implicit kinetic scheme for multiscale heat transfer problem accounting for phonon dispersion and polarization
74
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
An efficient implicit kinetic scheme is developed to solve the stationary phonon Boltzmann transport equation (BTE) based on the non-gray model including the phonon dispersion and polarization. Due to the wide range of the dispersed phonon mean free paths, the phonon transport under the non-gray model is essentially multiscale, and has to be solved differently and appropriately for varied phonon frequencies and branches. The proposed implicit kinetic scheme is composed of a microscopic iteration and a macroscopic iteration. The microscopic iteration is capable of automatically adapting with varied phonon mean free path of each phonon frequency and branch through solving the phonon BTE. The energy transfer of all phonons is gathered together by the microscopic iteration to evaluate the heat flux. The temperature field is predicted through a macroscopic heat transfer equation according to the heat flux, and the equilibrium state in the phonon BTE is also updated. The combination of the phonon BTE solver and the macroscopic equation makes the present method very efficient in a wide length scale. Three numerical tests, including the cross-plane, in-plane and nano-porous heat transfer in silicon, validate that the present scheme can handle with the phonon dispersion and polarization correctly and predict the multiscale heat transfer phenomena efficiently in a wide range. The present method could be tens of times faster than the typical implicit DOM and keeps the same amount of the memory requirements as the Fourier solver for multiscale heat transfer problem.
rate research
Read More
A novel hybrid computational method based on the discrete-velocity (DV) approximation, including the lattice-Boltzmann (LB) technique, is proposed. Numerical schemes for the kinetic equations are used in regions of rarefied flows, and LB schemes are employed in continuum flow zones. The schemes are written under the finite-volume (FV) formulation to achieve the flexibility of local mesh refinement. The truncated Hermite polynomial expansion is used for matching of DV and LB solutions. Special attention is paid to preserving conservation properties in the coupling algorithm. The test results obtained for the Couette flow of a rarefied gas are in excellent agreement with the benchmark solutions, mostly thanks to mesh refinement (both in the physical and velocity spaces) in the Knudsen layer.
Magnetization dynamics in magnetic materials is modeled by the Landau-Lifshitz-Gilbert (LLG) equation. In the LLG equation, the length of magnetization is conserved and the system energy is dissipative. Implicit and semi-implicit schemes have been used in micromagnetics simulations due to their unconditional numerical stability. In more details, implicit schemes preserve the properties of the LLG equation, but solve a nonlinear system of equations per time step. In contrast, semi-implicit schemes only solve a linear system of equations, while additional operations are needed to preserve the length of magnetization. It still remains unclear which one shall be used if both implicit and semi-implicit schemes are available. In this work, using the implicit Crank-Nicolson (ICN) scheme as a benchmark, we propose to make this implicit scheme semi-implicit. It can be proved that both schemes are second-order accurate in space and time. For the unique solvability of nonlinear systems of equations in the ICN scheme, we require that the temporal step size scales quadratically with the spatial mesh size. It is numerically verified that the convergence of the nonlinear solver becomes slower for larger temporal step size and multiple magnetization profiles are obtained for different initial guesses. The linear systems of equations in the semi-implicit CN (SICN) scheme are unconditionally uniquely solvable, and the condition that the temporal step size scales linearly with the spatial mesh size is needed in the convergence of the SICN scheme. In terms of numerical efficiency, the SICN scheme achieves the same accuracy as the ICN scheme with less computational time. Based on these results, we conclude that a semi-implicit scheme is superior to its implicit analog both theoretically and numerically, and we recommend the semi-implicit scheme in micromagnetics simulations if both methods are available.
In this paper we present and implement the Palindromic Discontinuous Galerkin (PDG) method in dimensions higher than one. The method has already been exposed and tested in [4] in the one-dimensional context. The PDG method is a general implicit high order method for approximating systems of conservation laws. It relies on a kinetic interpretation of the conservation laws containing stiff relaxation terms. The kinetic system is approximated with an asymptotic-preserving high order DG method. We describe the parallel implementation of the method, based on the StarPU runtime library. Then we apply it on preliminary test cases.
In this paper, authors focus effort on improving the conventional discrete velocity method (DVM) into a multiscale scheme in finite volume framework for gas flow in all flow regimes. Unlike the typical multiscale kinetic methods unified gas-kinetic scheme (UGKS) and discrete unified gas-kinetic scheme (DUGKS), which concentrate on the evolution of the distribution function at the cell interface, in the present scheme the flux for macroscopic variables is split into the equilibrium part and the nonequilibrium part, and the nonequilibrium flux is calculated by integrating the discrete distribution function at the cell center, which overcomes the excess numerical dissipation of the conventional DVM in the continuum flow regime. Afterwards, the macroscopic variables are finally updated by simply integrating the discrete distribution function at the cell center, or by a blend of the increments based on the macroscopic and the microscopic systems, and the multiscale property is achieved. Several test cases, involving unsteady, steady, high speed, low speed gas flows in all flow regimes, have been performed, demonstrating the good performance of the multiscale DVM from free molecule to continuum Navier-Stokes solutions and the multiscale property of the scheme is proved.
The Reynolds-Averaged Navier-Stokes equations and the Large-Eddy Simulation equations can be coupled using a transition function to switch from a set of equations applied in some areas of a domain to the other set in the other part of the domain. Following this idea, different time integration schemes can be coupled. In this context, we developed a hybrid time integration scheme that spatially couples the explicit scheme of Heun and the implicit scheme of Crank and Nicolson using a dedicated transition function. This scheme is linearly stable and second-order accurate. In this paper, an extension of this hybrid scheme is introduced to deal with a temporal adaptive procedure. The idea is to treat the time integration procedure with unstructured grids as it is performed with Cartesian grids with local mesh refinement. Depending on its characteristic size, each mesh cell is assigned a rank. And for two cells from two consecutive ranks, the ratio of the associated time steps for time marching the solutions is $2$. As a consequence, the cells with the lowest rank iterate more than the other ones to reach the same physical time. In a finite-volume context, a key ingredient is to keep the conservation property for the interfaces that separate two cells of different ranks. After introducing the different schemes, the paper recalls briefly the coupling procedure, and details the extension to the temporal adaptive procedure. The new time integration scheme is validated with the propagation of 1D wave packet, the Sods tube, and the transport of a bi-dimensional vortex in an uniform flow.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
