Phase-field simulation of dendritic solidification using a full threaded tree with adaptive meshing

Simulation of the microstructure evolution during solidifi cation is greatly benefi cial to the control of solidifi cation microstructures. A phase-fi eld method based on the full threaded tree(FTT) for the simulation of casting solidifi cation microstructure was proposed in this paper, and the structure of the full threaded tree and the mesh refi nement method was discussed. During dendritic growth in solidifi cation, the mesh for simulation is adaptively refi ned at the liquid-solid interface, and coarsened in other areas. The numerical results of a threedimension dendrite growth indicate that the phase-fi eld method based on FTT is suitable for microstructure simulation. Most importantly, the FTT method can increase the spatial and temporal resolutions beyond the limits imposed by the available hardware compared with the conventional uniform mesh. At the simulation time of 0.03 s in this study, the computer memory used for computation is no more than 10 MB with the FTT method, while it is about 50 MB with the uniform mesh method. In addition, the proposed FTT method is more effi cient in computation time when compared with the uniform mesh method. It would take about 20 h for the uniform mesh method, while only 2 h for the FTT method for computation when the solidifi cation time is 0.17 s in this study.

A three-dimensional robust volume-of-fluid solver based on the adaptive mesh refinement

The present study provides a three-dimensional volume-of-fluid method based on the adaptive mesh refinement technique.The projection method on the adaptive mesh is introduced for solving the incompressible Navier-Stokes equations.The octree structure mesh is employed to solve the flow velocities and the pressure.The developed solver is applied to simulate the deformation of the cubic droplet driven by the surface tension without the effect of the gravity.The numerical results well predict the shape evolution of the droplet.

An Adaptive Nonhydrostatic Atmospheric Dynamical Core Using a Multi-Moment Constrained Finite Volume Method

An adaptive 2 D nonhydrostatic dynamical core is proposed by using the multi-moment constrained finite-volume(MCV) scheme and the Berger-Oliger adaptive mesh refinement(AMR) algorithm. The MCV scheme takes several pointwise values within each computational cell as the predicted variables to build high-order schemes based on single-cell reconstruction. Two types of moments, such as the volume-integrated average(VIA) and point value(PV), are defined as constraint conditions to derive the updating formulations of the unknowns, and the constraint condition on VIA guarantees the rigorous conservation of the proposed model. In this study, the MCV scheme is implemented on a height-based, terrainfollowing grid with variable resolution to solve the nonhydrostatic governing equations of atmospheric dynamics. The AMR grid of Berger-Oliger consists of several groups of blocks with different resolutions, where the MCV model developed on a fixed structured mesh can be used directly. Numerical formulations are designed to implement the coarsefine interpolation and the flux correction for properly exchanging the solution information among different blocks. Widely used benchmark tests are carried out to evaluate the proposed model. The numerical experiments on uniform and AMR grids indicate that the adaptive model has promising potential for improving computational efficiency without losing accuracy.

ADAPTIVE HYBRID CARTESIAN GRID METHOD FOR VORTEX-DOMINATED FLOWS

An efficient compressible Euler equation solver for vortex-dominated flows is presented based on the adaptive hybrid Cartesian mesh and vortex identifying method.For most traditional grid-based Euler solvers,the excessive numerical dissipation is the great obstruction for vortex capturing or tracking problems.A vortex identifying method based on the curl of velocity is used to identify the vortex in flow field.Moreover,a dynamic adaptive mesh refinement(DAMR)process for hybrid Cartesian gird system is employed to track and preserve vortex.To validate the proposed method,a single compressible vortex convection flow is involved to test the accuracy and efficiency of DAMR process.Additionally,the vortex-dominated flow is investigated by the method.The obtained results are shown as a good agreement with the previous published data.

A Numerical Model for Simulating Two-Phase Flow with Adaptive Mesh Refinement

In this study,a numerical model for simulating two-phase flow is developed.The Cartesian grid with Adaptive Mesh Refinement(AMR)is adopted to reduce the computational cost.An explicit projection method is used for the time integration and the Finite Difference Method(FDM)is applied on a staggered grid for the discretization of spatial derivatives.The Volume of Fluid(VOF)method with Piecewise-Linear Interface Calculation(PLIC)is extended to the AMR grid to capture the gas-water interface accurately.A coarse-fine interface treatment method is developed to preserve the flux conservation at the interfaces.Several two-dimensional(2D)and three-dimensional(3D)benchmark cases are carried out for the validation of the model.2D and 3D shear flow tests are conducted to validate the extension of the VOF method to the AMR grid.A 2D linear sloshing case is considered in which the model is proved to have 2nd-order accuracy in space.The efficiency of applying the AMR grid is discussed with a nonlinear sloshing problem.Finally,2D solitary wave past stage and 2D/3D dam break are simulated to demonstrate that the model is able to simulate violent interface problems.

Adaptive Moving Mesh Central-Upwind Schemes for Hyperbolic System of PDEs:Applications to Compressible Euler Equations and Granular Hydrodynamics

We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.

Adaptive Finite Element Method for Steady Convection-Diffusion Equation

This paper examines the numerical solution of the convection-diffusion equation in 2-D. The solution of this equation possesses singularities in the form of boundary or interior layers due to non-smooth boundary conditions. To overcome such singularities arising from these critical regions, the adaptive finite element method is employed. This scheme is based on the streamline diffusion method combined with Neumann-type posteriori estimator. The effectiveness of this approach is illustrated by different examples with several numerical experiments.

Cavitation Evolution Around a Twist Hydrofoil by Large Eddy Simulation(LES)with Mesh Adaption

The cavitating flow around a Delft Twist-11 hydrofoil is simulated using the large eddy simulation approach.The volume-of-fluid method incorporated with the Schnerr-Sauer cavitation model is utilized to track the water-vapor interface.Adaptive mesh refinement(AMR)is also applied to improve the simulation accuracy automatically.Two refinement levels are conducted to verify the dominance of AMR in predicting cavitating flows.Results show that cavitation features,including the U-type structure of shedding clouds,are consistent with experimental observations.Even a coarse mesh can precisely capture the phase field without increasing the total cell number significantly using mesh adaption.The predicted shedding frequency agrees fairly well with the experimental data under refinement level 2.This study illustrates that AMR is a promising approach to achieve accurate simulations for multiscale cavitating flows within limited computational costs.Finally,the force element method is currently adopted to investigate the lift and drag fluctuations during the evolution of cavitation structure.The mechanisms of lift and drag fluctuations due to cavitation and the interaction between vorticity forces and cavitation are explicitly revealed.

Evolution of vortices in the wake of an ARJ21 airplane: Application of the lift-drag model

Wake separation is crucial to aircraft landing safety and is an important factor in airport operational efficiency.The near-ground evolution characteristics of wake vortices form the foundation of the wake separation system design.In this study,we analysed the near-ground evolution of vortices in the wake of a domestic aircraft ARJ21 initialised by the lift-drag model using large eddy simulations based on an adaptive mesh.Evolution of wake vortices formed by the main wing,flap and horizontal tail was discussed in detail.The horizontal tail vortices are the weakest and dissipate rapidly,whereas the flap vortices are the strongest and induce the tip vortex to merge with them.The horizontal tail and flap of an ARJ21 do not significantly influence the circulation evolution,height change and movement trajectory of the wake vortices.The far-field evolution of wake vortices can therefore be analysed using the conventional wake vortex model.

A study of a supersonic capsule/rigid disk-gap-band parachute system using large-eddy simulation

The aerodynamic performances and flow features of the capsule/rigid disk-gap-band(DGB)parachute system from the Mach number 1.8 to 2.2 are studied.We use the adaptive mesh refinement(AMR),the hybrid tuned center-difference and weighted essentially non-oscillatory(TCD-WENO)scheme,and the large-eddy simulation(LES)with the stretched-vortex subgrid model.The simulations reproduce complex interaction of the flow structures,including turbulent wakes and bow shocks,as well as bow shocks and expansion waves.The results show that the calculated aerodynamic drag coefficient agrees well with the previou simulation.Both the aerodynamic drag coefficient and the aerodynamic drag oscillation of the parachute system decrease with the increase of the initial Mach number of the fluid.It is found that the position and angle of the bow shock ahead of the canopy change as the Mach number increases,which makes the flow inside the canopy and the turbulent wake behind the canopy more complex and unstable.

Numerical Solution of Quasilinear Singularly Perturbed Problems by the Principle of Equidistribution

In this paper, the numerical solution and its error analysis of quasilinear singular perturbation two-point boundary value problems based on the principle of equidistribution are given. On the non-uniform grid of the uniformly distributed arc-length monitor function, the solution of the simple upwind scheme is obtained. It is proved that the adaptive simple upwind scheme based on the principle of equidistribution has uniform convergence for small perturbation parameters. Numerical experiments are carried out and the error analysis are confirmed.

Simple and robust h-adaptive shock-capturing method for flux reconstruction framework

In this paper,a simple and robust shock-capturing method is developed for the Flux Reconstruction(FR)framework by combining the Adaptive Mesh Refinement(AMR)technique with the positivity-preserving property.The adaptive technique avoids the use of redundant meshes in smooth regions,while the positivity-preserving property makes the solver capable of providing numerical solutions with physical meaning.The compatibility of these two significant features relies on a novel limiter designed for mesh refinements.It ensures the positivity of solutions on all newly created cells.Therefore,the proposed method is completely positivity-preserving and thus highly robust.It performs well in solving challenging problems on highly refined meshes and allows the transition of cells at different levels to be completed within a very short distance.The performance of the proposed method is examined in various numerical experiments.When solving Euler equations,the technique of Local Artificial Diffusivity(LAD)is additionally coupled to damp oscillations.More importantly,when solving Navier-Stokes equations,the proposed method requires no auxiliaries and can provide satisfying numerical solutions directly.The implementation of the method becomes rather simple.

A Study on CFL Conditions for the DG Solution of Conservation Laws on Adaptive Moving Meshes

The selection of time step plays a crucial role in improving stability and efficiency in the Discontinuous Galerkin(DG)solution of hyperbolic conservation laws on adaptive moving meshes that typically employs explicit stepping.A commonly used selection of time step is a direct extension based on Courant-Friedrichs-Levy(CFL)conditions established for fixed and uniform meshes.In this work,we provide a mathematical justification for those time step selection strategies used in practical adaptive DG computations.A stability analysis is presented for a moving mesh DG method for linear scalar conservation laws.Based on the analysis,a new selection strategy of the time step is proposed,which takes into consideration the coupling of theα-function(that is related to the eigenvalues of the Jacobian matrix of the flux and the mesh movement velocity)and the heights of the mesh elements.The analysis also suggests several stable combinations of the choices of theα-function in the numerical scheme and in the time step selection.Numerical results obtained with a moving mesh DG method for Burgers’and Euler equations are presented.For comparison purpose,numerical results obtained with an error-based time step-size selection strategy are also given。

LSFEM Implementation of MHD Numerical Solver

Many problems in physics are inherently of multi-scale nature. The issues of MHD turbulence or magnetic reconnection, namely in the hot and sparse, almost collision-less astrophysical plasmas, can stand as clear examples. The Finite Element Method (FEM) with adaptive gridding appears to be the appropriate numerical implementation for handling the broad range of scales contained in such high Lundquist-number MHD problems. In spite the FEM is now routinely used in engineering practice in solid-state and fluid dynamics, its usage for MHD simulations has recently only begun and only few implementations exist so far. In this paper we present our MHD solver based on the Least-Square FEM (LSFEM) formulation. We describe the transformation of the MHD equations into form required for finding the LSFEM functional and some practical issues in implementation of the method. The algorithm was tested on selected problems of ideal (non-resistive) and resistive MHD. The tests show the usability of LSFEM for solving MHD equations.

Large eddy simulation of cavitating flows with dynamic adaptive mesh refinement using OpenFOAM

Cavitating flows are dominated by large gradients of physical properties and quantities containing complicated interfacial structures and lots of multi-scale eddies that need to be accurately characterized using a high-resolution mesh.The present work,within OpenFOAM,proposes an effective modeling framework using the large eddy simulation(LES)approach along with the volume of fluid(VOF)method to simulate the two-phase flow system and applies the Schnerr-Sauer model to calculate the mass-transfer rate between water and vapor.The adaptive mesh refinement(AMR)which is a powerful tool for allocating high-resolution grids only to the region of the greatest concern is adopted for improving the solution of interfacial structures.The effect of grid size is firstly investigated and the time-averaged quantities are verified against the experimental data,and then simulations of cavitating flows are successfully achieved to precisely characterize the features of cavitation with automatically and dynamically refining the mesh.As the refinement only takes place in the interfacial region,high-precision simulations can be achieved with limited computational resources,and the method shows promising prospects for modeling of the multi-scale,time-critical and computationally intensive cavitating flows.

Multiscale Basis Functions for Singular Perturbation on Adaptively Graded Meshes

We apply the multiscale basis functions for the singularly perturbed reaction-diffusion problem on adaptively graded meshes,which can provide a good balance between the numerical accuracy and computational cost.The multiscale space is built through standard finite element basis functions enriched with multiscale basis functions.The multiscale basis functions have abilities to capture originally perturbed information in the local problem,as a result our method is capable of reducing the boundary layer errors remarkably on graded meshes,where the layer-adapted meshes are generated by a given parameter.Through numerical experiments we demonstrate that the multiscale method can acquire second order convergence in the L^(2)norm and first order convergence in the energy norm on graded meshes,which is independent ofε.In contrast with the conventional methods,our method is much more accurate and effective.

Adaptive Conservative Cell Average Spectral Element Methods for Transient Wigner Equation in Quantum Transport

A new adaptive cell average spectral element method(SEM)is proposed to solve the time-dependent Wigner equation for transport in quantum devices.The proposed cell average SEM allows adaptive non-uniform meshes in phase spaces to reduce the high-dimensional computational cost of Wigner functions while preserving exactly the mass conservation for the numerical solutions.The key feature of the proposed method is an analytical relation between the cell averages of the Wigner function in the k-space(local electron density for finite range velocity)and the point values of the distribution,resulting in fast transforms between the local electron density and local fluxes of the discretized Wigner equation via the fast sine and cosine transforms.Numerical results with the proposed method are provided to demonstrate its high accuracy,conservation,convergence and a reduction of the cost using adaptive meshes.

Adaptive mesh refinement immersed boundary method for simulations of laminar flows past a moving thin elastic structure

One of the critical issues in numerical simulation of fluid-structure interaction problems is inaccuracy of the solutions,especially for flows past a stationary thin elastic structure where large deformations occur.High resolution is required to capture the flow characteristics near the fluid-structure interface to enhance accuracy of the solutions within proximity of the thin deformable body.Hence,in this work,an algorithm is developed to simulate fluid-structure interactions of moving deformable structures with very thin thicknesses.In this algorithm,adaptive mesh refinement(AMR)is integrated with immersed boundary finite element method(IBFEM)with two-stage pressure-velocity corrections.Despite successive interpolation of the flow field by IBM,the governing equations were solved using a fixed structured mesh,which significantly reduces the computational time associated with mesh reconstruction.The cut-cell IBM is used to predict the body forces while FEM is used to predict deformation of the thin elastic structure in order to integrate the motions of the fluid and solid at the interface.AMR is used to discretize the governing equations and obtain solutions that efficiently capture the thin boundary layer at the fluid-solid interface.The AMR-IBFEM algorithm is first verified by comparing the drag coefficient,lift coefficient,and Strouhal number for a benchmark case(laminar flow past a circular cylinder at Re=100)and the results showed good agreement with those of other researchers.The algorithm is then used to simulate 2-D laminar flows past stationary and moving thin structures positioned perpendicular to the freestream direction.The results also showed good agreement with those obtained from the arbitrary Lagrangian-Eulerian(ALE)algorithm for elastic thin boundaries.It is concluded that the AMR-IBFEM algorithm is capable of predicting the characteristics of laminar flow past an elastic structure with acceptable accuracy(error of-0.02%)with only-1%of the computational time for simulations with full mesh refinement.

AnAdaptive Moving Mesh Method for Two-Dimensional Relativistic Hydrodynamics

This paper extends the adaptive moving mesh method developed by Tang and Tang[36]to two-dimensional(2D)relativistic hydrodynamic(RHD)equations.The algorithm consists of two“independent”parts:the time evolution of the RHD equations and the(static)mesh iteration redistribution.In the first part,the RHD equations are discretized by using a high resolution finite volume scheme on the fixed but nonuniform meshes without the full characteristic decomposition of the governing equations.The second part is an iterative procedure.In each iteration,the mesh points are first redistributed,and then the cell averages of the conservative variables are remapped onto the new mesh in a conservative way.Several numerical examples are given to demonstrate the accuracy and effectiveness of the proposed method.

An Adaptive Mesh Refinement Strategy for Immersed Boundary/Interface Methods

An adaptive mesh refinement strategy is proposed in this paper for the Immersed Boundary and Immersed Interface methods for two-dimensional elliptic interface problems involving singular sources.The interface is represented by the zero level set of a Lipschitz functionϕ(x,y).Our adaptive mesh refinement is done within a small tube of|ϕ(x,y)|δwith finer Cartesian meshes.The discrete linear system of equations is solved by a multigrid solver.The AMR methods could obtain solutions with accuracy that is similar to those on a uniform fine grid by distributing the mesh more economically,therefore,reduce the size of the linear system of the equations.Numerical examples presented show the efficiency of the grid refinement strategy.