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.

Combined adaptive meshing technique and characteristic-based split algorithm for viscous incompressible flow analysis

A combined characteristic-based split algorithm and all adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented. Tile method uses the three-node triangular element with equal-order interpolation functions for all variables of tile velocity components and pressure. The main advantage of the combined nlethod is that it inlproves the sohltion accuracy by coupling an error estinla- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change ill sohmtion gradients, mid at the same time, larger elements in the other regions. The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder, for their transient and steady-state flow behaviors.

DEVELOPMENT AND APPLICATION OF ADAPTIVE MESH TECHNIQUE IN THREE DIMENSIONAL NUMERICAL SIMULATION OF WELDING PROCESS

The adaptive mesh mesh technique is developed and applied in three dimensional numerical simulation of welding process on the base of the commercial software. Special user subroutine is worked out to accom- plish this function.This technique can make the dense mesh moving simultaneously with the heat source while the other area of the structure with much coarser mesh, greatly reducing the number of nodes and elements in the analysis.Temperature field,displacement and stress distributions during welding pro- cess me analyzed by FEM method with adaptive mesh and the analysis is also conducted with normal FEM method. The temperature field,displacement and stress distributions obtained with both methods are shown in contrast. The results show that the temperature fields and the displacement distributions of simulation on adaptive mesh correspond wery well with that of without adaptive mesh. Though the stress distributions have some difference,but the trends of the stress distribution are corresponding.The com- parison of the computing time of the two meshes indicates that the adaptive that the adaptive mesh can greatly reduce the calculation time when used for welding process.

Two-Dimensional Euler Adaptive Mesh Method on Detonation

The adaptive mesh refinement （AMR） method is applied in the 2-D Euler multi-component elasticplastic hydrodynamics code （MEPH2Y）. It is applied on detonation. Firstly, the AMR method is described, including a conservative spatial interpolation, the time integration methodology with the adapitve time increment and an adaptive computational region method. The advantage of AMR technique is exhibited by numerical examples, including the 1-D C-J detonation and the 2-D implosion ignited from a single point. Results show that AMR can promote the computational efficiency, keeping the accuracy in interesting regions.

High-order schemes for predicting computational aeroacoustic propagation with adaptive mesh refinement

High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical issues associated with high-order schemes on an adaptively refined mesh, such as stability and accuracy are addressed. Several CAA benchmark problems are used to demonstrate the feasibility and efficiency of the approach.

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.

METHOD FOR ADAPTIVE MESH GENERATION BASED ON GEOMETRICAL FEATURES OF 3D SOLID

In order to provide a guidance to specify the element size dynamically during adaptive finite element mesh generation, adaptive criteria are firstly defined according to the relationships between the geometrical features and the elements of 3D solid. Various modes based on different datum geometrical elements, such as vertex, curve, surface, and so on, are then designed for generating local refined mesh. With the guidance of the defmed criteria, different modes are automatically selected to apply on the appropriate datum objects to program the element size in the local special areas. As a result, the control information of element size is successfully programmed covering the entire domain based on the geometrical features of 3D solid. A new algorithm based on Delatmay triangulation is then developed for generating 3D adaptive finite element mesh, in which the element size is dynamically specified to catch the geometrical features and suitable tetrahedron facets are selected to locate interior nodes continuously. As a result, adaptive mesh with good-quality elements is generated. Examples show that the proposed method can be successfully applied to adaptive finite element mesh automatic generation based on the geometrical features of 3D solid.

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.

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.

Coupled flood and sediment transport modelling with adaptive mesh refinement

Coupled flood and sediment transport modelling in large-scale domains has for long been hindered by the high computational cost.Adaptive mesh refinement is one of the viable ways to solving this problem without degrading the accuracy.This goal can be accomplished through mesh adaptation,e.g.,mesh coarsening and refining based on the dynamic regime of the flow and sediment transport along with bed evolution.However,previous studies in this regard have been limited to cases either without involving sediment transport or featuring flow-sediment-bed decoupling and the assumption of sediment transport capacity,which are not generally justified.Here,a coupled hydrodynamic and non-capacity sediment transport model is developed on adaptive non-uniform rectangular mesh.The proposed model is validated against experimental tests and numerical results based on the fixed meshes.It is demonstrated that the proposed model can properly capture shock waves,resolve the wetting/drying transition and reproduce morphological evolution.Compared with models based on the fixed meshes,the proposed model features great advantage in computational efficiency and holds promise for wide applications.

High-order compact finite volume methods on unstructured grids with adaptive mesh refinement for solving inviscid and viscous flows

In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.

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.

DISCONTINUITY-CAPTURING FINITE ELEMENT COMPUTATION OF UNSTEADY FLOW WITH ADAPTIVE UNSTRUCTURED MESH

A discontinuity-capturing scheme of finite element method(FEM)is proposed.The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unsteady flows,which exhibits the capability of capturing the shock waves and/or thin shear layers accurately in an unsteady viscous flow at high Reynolds number. In particular,a new testing variable,i.e.,the disturbed kinetic energy E,is suggested and used in the adaptive mesh computation,which is universally applicable to the capturing of both shock waves and shear layers in the inviscid flow and viscous flow at high Reynolds number.Based on several calculated examples,this approach has been proved to be effective and efficient for the calculations of compressible and incompressible flows.

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.

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.

3D Simulation of Flow with Free Surface Based on Adaptive Octree Mesh System

The technique of adaptive tree mesh is an effective way to reduce computational cost through automatic adjustment of cell size according to necessity. In the present study, the 2D numerical N-S solver based on the adaptive quadtree mesh system was extended to a 3D one, in which a spatially adaptive oetree mesh system and multiple particle level set method were adopted for the convenience to deal with the air-water-structure multiple-medium coexisting domain. The stretching process of a dumbbell was simulated and the results indicate that the meshes are well adaptable to the free surface. The collapsing process of water column impinging a circle cylinder was simulated and from the results, it can be seen that the processes of fluid splitting and merging are properly simulated. The interaction of second-order Stokes waves with a square cylinder was simulated and the obtained drag force is consistent with the result by the Morison＇s wave force formula with the coefficient values of the stable drag component and the inertial force component bein~ set as 2.54.

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.

A Computational Framework for Parachute Inflation Based on Immersed Boundary/Finite Element Approach

A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.

COMBINED DELAUNAY TRIANGULATION AND ADAPTIVE FINITE ELEMENT METHOD FOR CRACK GROWTH ANALYSIS

The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around crack tips and large elements in the other regions. The resulting stress intensity factors and simulated crack propagation behavior are used to evaluate the effectiveness of the procedure. Three sample problems of a center cracked plate, a single edge cracked plate and a compact tension specimen, are simulated and their results assessed.

Dynamic Adaptive Finite Element Analysis of Acoustic Wave Propagation Due to Underwater Explosion for Fluid-structure Interaction Problems

In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.