Source author record

Chunlei Liang

Chunlei Liang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR physics.comp-ph physics.flu-dyn astro-ph.IM math.NA

Catalog footprint

What is connected

6works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Simulations of Global Solar Convection with a Fully Compressible CHORUS++ Code

Fluid-dynamics simulations of global solar convection are a critically important tool for assessing the dynamics of the solar interior. However, simulation studies with a fully compressible hydrodynamics code are not yet common. The CHORUS++ code solves robustly and efficiently the fully compressible hydrodynamics equations using a compact local spectral method and semi-unstructured grid system. Using the CHORUS++ code, we simulate the solar interior plasma flows from 0.7 to 0.99 of the solar radius using the values of the solar total luminosity and the sidereal rotation rate. In this paper, we analyze the simulated global flow structures before the statistically stable state to assess the compressibility of the plasma flows obtained with the fully compressible hydrodynamic code. The divergence of mass flux with the compressible model is overall small in the examined state before reaching the equilibrium state of the convection, which implies that the differences between the fully compressible flows and those obtained with the anelastic, incompressible, or linear-equation models are small in the simulated inner part of the convection zone. Although a fully-relaxed stationary convection has not been achieved yet in the examined state, the model qualitatively reproduces the solar-type differential rotation. Simulations for longer periods are needed to achieve the system relaxation state. This work, assessing the early phase of the hydrodynamic evolution of solar convection, is our first step toward a better understanding of the nature of solar convection and the dynamo processes.

preprint2022arXiv

High-Order Large-Eddy Simulations of a Wind Turbine in Ducted and Open-Rotor Configurations

High-order large-eddy simulations are performed to study the performance and flow field of a ducted wind turbine operating at different tip speed ratios. To evaluate the effects of the duct, simulations with the same tip speed ratios are also performed on the corresponding open-rotor turbine. It is found that the ducted turbine consistently obtains higher power outputs than the open-rotor counterpart, and the duct itself enhances flow turbulence and blade trailing-edge vortices but weakens tip and hub vortices. Flow bifurcation is observed at the largest tip speed ratio and is identified to be caused by blade blockage effects. Comparative simulations are also performed on both turbines under different yaw angles. It is noticed that the ducted configuration is insensitive to small yaw angles and maintains higher power outputs than the open-rotor configuration at all yaw angles. Moreover, it is observed that the wakes of both configurations recover more quickly as the yaw angle increases.

preprint2020arXiv

An open boundary condition for high-order solutions of magnetohydrodynamics on unstructured grids

In this paper a characteristics-based open boundary condition (CBC) is proposed for the magnetohydrodynamic (MHD) system of equations. The algorithm is carefully designed and implemented in the context of a high-order flux reconstruction (FR) scheme under the Generalized Lagrange Multiplier (GLM)-MHD system of equations. It is implemented by adding the contribution of the characteristic equation directly to the corrected flux term in the FR scheme dispensing with solving time-dependent characteristic equations along boundary faces. The CBC method is shown to be more accurate and robust than commonly used zero normal derivative (ZND) and approximate Riemann solver boundary conditions (ARBC) in solving 1D, 2D, and 3D test problems. The CBC method is successfully applied to simulate challenging problems of magnetic reconnection for which other options failed to get stable results over long-period time integration.

preprint2015arXiv

A Compressible High-Order Unstructured Spectral Difference Code for Stratified Convection in Rotating Spherical Shells

We present a novel and powerful Compressible High-ORder Unstructured Spectral-difference (CHORUS) code for simulating thermal convection and related fluid dynamics in the interiors of stars and planets. The computational geometries are treated as rotating spherical shells filled with stratified gas. The hydrodynamic equations are discretized by a robust and efficient high-order Spectral Difference Method (SDM) on unstructured meshes. The computational stencil of the spectral difference method is compact and advantageous for parallel processing. CHORUS demonstrates excellent parallel performance for all test cases reported in this paper, scaling up to 12,000 cores on the Yellowstone High-Performance Computing cluster at NCAR. The code is verified by defining two benchmark cases for global convection in Jupiter and the Sun. CHORUS results are compared with results from the ASH code and good agreement is found. The CHORUS code creates new opportunities for simulating such varied phenomena as multi-scale solar convection, core convection, and convection in rapidly-rotating, oblate stars.

preprint2015arXiv

A Simple, Efficient, High-order Accurate Sliding-Mesh Interface Approach to the Spectral Difference Method on Coupled Rotating and Stationary Domains

This paper presents a simple, efficient, and high-order accurate sliding-mesh interface approach to the spectral difference (SD) method. We demonstrate the approach by solving the two-dimensional compressible Navier-Stokes equations on quadrilateral grids. This approach is an extension of the straight mortar method originally designed for stationary domains by Kopriva, it employs curved dynamic mortars on sliding-mesh interfaces to couple rotating and stationary domains. On the nonconforming sliding-mesh interfaces, the related variables are first projected from cell faces to mortars to compute common fluxes, and then the common fluxes are projected back from the mortars to the cell faces to ensure conservation. To verify the spatial order of accuracy of the sliding-mesh spectral difference (SSD) method, both inviscid and viscous flow cases are tested. It is shown that the SSD method preserves the high-order accuracy of the SD method. Meanwhile, the SSD method is found to be very efficient in terms of computational cost. This novel sliding-mesh interface method is very suitable for parallel processing with domain decomposition. It can be applied to a wide range of problems, such as the aerodynamics of rotorcraft, wind turbines, and oscillating wing wind power generators, etc.

preprint2014arXiv

Turbulent convection in the Sun: modeling in unstructured meshes

We adopted an unstructured hydrodynamical solver CharLES to the problem of global convection in the Sun. With the aim to investigate the properties of solar turbulent convection and reproduce differential rotation pattern. We performed simulations in two spherical shells, with 1.3 and 10 million cells. In the first, coarse mesh, the solution does not reproduce realistic convection, and is dominated by numerical effects. In the second mesh, thermal conduction leads to cooling of bottom layers, that could not be compensated by solar irradiance. More simulations in the 10M cells mesh should be performed to investigate the influence of transport coefficients and numerical effects. Our estimate of the code performance suggests, that realistic simulations in even finer grids could be performed for reasonable computational cost.