An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （...An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.展开更多
To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any compl...To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any complicated boundary treatment or artificial diffusion, this method treated curved boundaries directly by utilizing the inherent merit of unstructured mesh in fitting curves. The PLIC (Piecewise Linear Interface Calculation) method was adopted to obtain a second-order accurate linearized reconstruction approximation and the MLER (Modified Lagrangian-Eulerian Re-map) method was introduced to advect fluid volumes on unstructured mesh. Moreover, an analytical relation for the interface’s line constant vs. the volume clipped by the interface was developed so as to improve the method’s efficiency. To validate this method, a comprehensive series of large straining advection tests were performed. Numerical results provide convincing evidences for the method’s high volume conservative accuracy and second-order shape error convergence rate. Also, a dramatic improvement on computational accuracy over its unstructured triangular mesh counterpart is checked.展开更多
In this paper, the accuracy of Chang＇s unstructured space-time conservation element and solution element （CE/SE） scheme is analysed for the first time. Based on a redefinition of conservation elements and solution ...In this paper, the accuracy of Chang＇s unstructured space-time conservation element and solution element （CE/SE） scheme is analysed for the first time. Based on a redefinition of conservation elements and solution elements, an improved two-dimensional （2D） unstructured CE/SE scheme with an adjustable parameter β is proposed to accurately capture shock waves. The new scheme can be applied to any type of grid without special treatnmnt. Compared with Chang＇s original parameter a, larger/5 dose not cost extra computational resources. Numerical tests reveal that the new scheme is not only clear in physical concept, compact and highly accurate but also more capable of capturing shock waves than the popular fifth-order accurate weighted essentially non-oscillatory scheme.展开更多
The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can ov...The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.展开更多
文摘An upwind scheme based on the unstructured mesh is developed to solve ideal 2-D magnetohydrodynamics （MHD） equations. The inviscid fluxes are approximated by using the modified advection upstream splitting method （AUSM） scheme, and a 5-stage explicit Runge-Kutta scheme is adopted in the time integration. To avoid the influence of the magnetic field divergence created during the simulation, the hyperbolic divergence cleaning method is introduced. The shock-capturing properties of the method are verified by solving the MHD shock-tube problem. Then the 2-D nozzle flow with the magnetic field is numerically simulated on the unstructured mesh. Computational results demonstrate the effects of the magnetic field and agree well with those from references.
基金the National Natural Science Foundation ofChina under Grant No. 50779043, 50779045
文摘To simulate two-dimensional free-surface flows with complex boundaries directly and accurately, a novel VOF (Volume-of-fluid) method based on unstructured quadrilateral mesh is presented. Without introducing any complicated boundary treatment or artificial diffusion, this method treated curved boundaries directly by utilizing the inherent merit of unstructured mesh in fitting curves. The PLIC (Piecewise Linear Interface Calculation) method was adopted to obtain a second-order accurate linearized reconstruction approximation and the MLER (Modified Lagrangian-Eulerian Re-map) method was introduced to advect fluid volumes on unstructured mesh. Moreover, an analytical relation for the interface’s line constant vs. the volume clipped by the interface was developed so as to improve the method’s efficiency. To validate this method, a comprehensive series of large straining advection tests were performed. Numerical results provide convincing evidences for the method’s high volume conservative accuracy and second-order shape error convergence rate. Also, a dramatic improvement on computational accuracy over its unstructured triangular mesh counterpart is checked.
基金supported by the National Natural Science Foundation of China (Grant Nos. 10732010 and 10972010)
文摘In this paper, the accuracy of Chang＇s unstructured space-time conservation element and solution element （CE/SE） scheme is analysed for the first time. Based on a redefinition of conservation elements and solution elements, an improved two-dimensional （2D） unstructured CE/SE scheme with an adjustable parameter β is proposed to accurately capture shock waves. The new scheme can be applied to any type of grid without special treatnmnt. Compared with Chang＇s original parameter a, larger/5 dose not cost extra computational resources. Numerical tests reveal that the new scheme is not only clear in physical concept, compact and highly accurate but also more capable of capturing shock waves than the popular fifth-order accurate weighted essentially non-oscillatory scheme.
文摘The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.
