Optimal Dual Schemes for Adaptive Grid Based Hexmeshing

Hexahedral meshes are a ubiquitous domain for the numerical resolution of partial differential equations. Computing a pure hexahedral mesh from an adaptively refined grid is a prominent approach to automatic hexmeshing, and requires the ability to restore the all hex property around the hanging nodes that arise at the interface between cells having different size. The most advanced tools to accomplish this task are based on mesh dualization. These approaches use topological schemes to regularize the valence of inner vertices and edges, such that dualizing the grid yields a pure hexahedral mesh. In this article, we study in detail the dual approach, and propose four main contributions to it: (i) We enumerate all the possible transitions that dual methods must be able to handle, showing that prior schemes do not natively cover all of them; (ii) We show that schemes are internally asymmetric, therefore not only their construction is ambiguous, but different implementative choices lead to hexahedral meshes with different singular structure; (iii) We explore the combinatorial space of dual schemes, selecting the minimum set that covers all the possible configurations and also yields the simplest singular structure in the output hexmesh; (iv) We enlarge the class of adaptive grids that can be transformed into pure hexahedral meshes, relaxing one of the tight topological requirements imposed by previous approaches. Our extensive experiments show that our transition schemes consistently outperform prior art in terms of ability to converge to a valid solution, amount and distribution of singular mesh edges, and element count. Last but not least, we publicly release our code and reveal a conspicuous amount of technical details that were overlooked in previous literature, lowering an entry barrier that was hard to overcome for practitioners in the field.

The pre -breakdown characteristics of weakly ionized media in the high non-uniform electric field using calculations with refined grid cells

The theoretical model of a pre-breakdown ions formation in a liquid dielectric and their flows, caused by high non-uniform electric field is represented. The 3D system of the macroscopic pre-breakdown electrohydrodynamic equations is written. Influence of the electric field on the molecule dissociation rate is taken into account. The system includes the Poison equation for electric field potential, the equation of ions formation and the Navier-Stokes equations with electric force. In addition, results of modelling using refined mesh are presented.

An Adaptive Refined Grid Search Strategy for Assessing Operational Flexibility and Application on Refrigerant Selection

A novel adaptive refined grid search strategy is developed for representative characterization of process feasible region boundaries and accurate estimation of its hypervolume. In particular, a linked list data structure adopted from the field of computer science is used to maintain the grid connectivity information. A uniform perturbation scheme is used to refine the search only near boundaries. The volumetric flexibility index FI_V is calculated directly from a summation of feasible hypercubes in the grid, without the need to apply shape reconstruction techniques. The proposed adaptive grid search strategy can capture complex region shapes with reduced sampling costs and without randomness for better reproducibility. Operational flexibility is optimized traditionally at a process scale. A case study on refrigerant selection is presented to demonstrated that the developed strategy could be combined within a computer-aided molecular design framework for operational flexibility optimization in molecular scale.

Vertical grid refinement for stratocumulus clouds in the radiation scheme of the global climate model ECHAM6.3-HAM2.3-P3

In this study, we implement a vertical grid refinement scheme in the radiation routine of the global aerosol–climate model ECHAM-HAM, aiming to improve the representation of stratocumulus clouds and address the underestimation of their cloud cover. The scheme is based on a reconstruction of the temperature inversion as a physical constraint for the cloud top. On the refined grid, the boundary layer and the free troposphere are separated and the cloud's layer is made thinner. The cloud cover is recalculated either by conserving the cloud volume (SC-VOLUME) or by using the Sundqvist cloud cover routine on the new grid representation (SC-SUND). In global climate simulations, we find that the SC-VOLUME approach is inadequate, as there is a mismatch, in most cases, between the layer of the inversion and the layer of the stratocumulus cloud, which prevents its application and is itself likely caused by an overly low vertical resolution. Additionally, we find that the occurrence frequency of stratocumulus clouds is underestimated in ECHAM-HAM, limiting a priori the potential benefits of a scheme like SC-VOLUME targeting only cloud amount when present. With the SC-SUND approach, the possibility for new clouds to be formed on the refined grid results in a large increase in mean total cloud cover in stratocumulus regions. In both cases, however, the changes exerted in the radiation routine are too weak to produce a significant improvement in the simulated stratocumulus cloud cover. We investigate and discuss the reasons behind this. The grid refinement scheme could be used more effectively for this purpose if implemented directly in the model's cloud microphysics and cloud cover routines, but other possible ways forward are also discussed.

CFD–DEM Simulation of Sand-Retention Mechanisms in Slurry Flow

The primary motivation of this paper is to investigate the sand-retention mechanisms that occur at the opening of sand filters. Various retention mechanisms under various conditions are explored that have a particulate flow with a low concentration of sand particles (called slurry flow) such as particle shape, size, and concentration. The computational fluid dynamic (CFD)–discrete element method (DEM) model is applied to predict the retention mechanisms under steady flow conditions of the well-bore. By using coupled CFD–DEM (CFD to model the fluid flow, and DEM to model the particle flow), the physics involved in the retention mechanisms is studied. The coarse grid unresolved and the smoothed unresolved (refined grid unresolved) coupling approaches implemented in STAR-CCM+ (SIEMENS PLM) are used to transfer data between the fluid and solid phases and calculate the forces. The filter slots under investigation have different geometries: straight, keystone, wire-wrapped screen (WWS) and seamed slot and the particles are considered with different shapes and different aspect ratios and size distributions. The flow regime is laminar in all simulations conducted. The CFD–DEM model is validated from the perspectives of particle–fluid, particle–particle, and particle–wall interactions. Verification of the CFD–DEM model is conducted by mesh sensitivity analysis to investigate the coupling resolution between the CFD and DEM. By simulation of numerous slurry flow scenarios, three retention mechanisms including surface deposition, size exclusion, and sequential arching of particles are observed. However, the concentration of particles is too diluted to result in multiparticle arch formation. In the simulations, various conditions are tested to give us an insight into the parameters and conditions that could affect the occurrence of the retention mechanisms. As an example, the importance of the gravity force and interaction forces on retention mechanisms are confirmed at the microscale in comparison with others forces involved in retention mechanisms such as the drag force, lift force, cohesive force, buoyancy force, and virtual mass force.

Vertical grid refinement for stratocumulus clouds in the radiation scheme of a global climate model

In this study, we implement a vertical grid refinement scheme in the radiation routine of the global aerosol-climate model ECHAM-HAM, aiming to improve the representation of stratocumulus clouds and address the underestimation of their cloud cover. The scheme is based on a reconstruction of the temperature inversion as a physical constraint for the cloud top. On the refined grid, the boundary layer and the free troposphere are separated and the cloud's layer is made thinner. The cloud cover is re-calculated either by conserving the cloud volume (SC-VOLUME) or by using the Sundqvist cloud cover routine on the new grid representation (SC-SUND). In global climate simulations, we find that the SC-VOLUME approach is inadequate, as in most cases there is a mismatch between the layer of the inversion and of the stratocumulus cloud, which prevents its application and is itself likely caused by too-low vertical resolution. With the SC-SUND approach, the possibility for new clouds to be formed on the refined grid results in a large increase in mean total cloud cover in stratocumulus regions. In both cases, however, the changes exerted in the radiation routine are too weak to produce a significant improvement of the simulated stratocumulus cloud cover. The grid refinement scheme could be used more effectively for this purpose if implemented directly in the model's cloud microphysics and cloud cover routines.

A multi-block adaptive solving technique based on lattice Boltzmann method

In this paper, a CFD parallel adaptive algorithm is self-developed by combining the multi-block Lattice Boltzmann Method (LBM) with Adaptive Mesh Refinement (AMR). The mesh refinement criterion of this algorithm is based on the density, velocity and vortices of the flow field. The refined grid boundary is obtained by extending outward half a ghost cell from the coarse grid boundary, which makes the adaptive mesh more compact and the boundary treatment more convenient. Two numerical examples of the backward step flow separation and the unsteady flow around circular cylinder demonstrate the vortex structure of the cold flow field accurately and specifically.