Phase-field modeling of dendritic growth under forced flow based on adaptive finite element method

A mathematical model combined projection algorithm with phase-field method was applied. The adaptive finite element method was adopted to solve the model based on the non-uniform grid, and the behavior of dendritic growth was simulated from undercooled nickel melt under the forced flow. The simulation results show that the asymmetry behavior of the dendritic growth is caused by the forced flow. When the flow velocity is less than the critical value, the asymmetry of dendrite is little influenced by the forced flow. Once the flow velocity reaches or exceeds the critical value, the controlling factor of dendrite growth gradually changes from thermal diffusion to convection. With the increase of the flow velocity, the deflection angle towards upstream direction of the primary dendrite stem becomes larger. The effect of the dendrite growth on the flow field of the melt is apparent. With the increase of the dendrite size, the vortex is present in the downstream regions, and the vortex region is gradually enlarged. Dendrite tips appear to remelt. In addition, the adaptive finite element method can reduce CPU running time by one order of magnitude compared with uniform grid method, and the speed-up ratio is proportional to the size of computational domain.

Adaptive Finite Element Modeling of Marine Controlled-Source Electromagnetic Fields in Two-Dimensional General Anisotropic Media

In this paper, we extend the scope of numerical simulations of marine controlled-source electromagnetic (CSEM) fields in a particular case of anisotropy (dipping anisotropy) to the general case of anisotropy by using an adaptive finite element approach. In comparison to a dipping anisotropy case, the first order spatial derivatives of the strike-parallel components arise in the partial differential equations for generally anisotropic media, which cause a non-symmetric linear system of equations for finite element modeling. The adaptive finite element method is employed to obtain numerical solutions on a sequence of refined unstructured triangular meshes, which allows for arbitrary model geometries including bathymetry and dipping layers. Numerical results of a 2D anisotropic model show both anisotropy strike and dipping angles have great influence on the marine CSEM responses.

Adaptive Finite Element Method Based on Optimal Error Estimates for Linear Elliptic Problems on Nonconvex Polygonal Domains

The subject of this work is to propose adaptive finite element methods based on an optimal maximum norm error control estimate.Using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions,the proposed procedures are analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that they generate the correct type of refinement and lead to the desired control under consideration.

3-D Marine CSEM Modeling in General Anisotropic Media by Using an Adaptive Finite Element Approach Based on the Vector-Scalar Potential

We present three-dimensional(3-D)modeling method of marine controlled-source electromagnetic(CSEM)fields in general anisotropic media using an adaptive finite element approach based on the vector-scalar potential.The modeling is based on the governing Helmholtz equations in the vector-scalar potential system.Unstructured tetrahedral grids are employed,which can exactly simulate the terrain relief and complex electrical structures.Moreover,based on the gradient recovery technology,the adaptive finite element approach is used to drive the mesh refinement,and make the finite element solutions converge gradually to the exact solutions.The primary-secondary field approach is used to improve the numerical accuracy of CSEM fields near the source point,where the primary field is calculated by using the quasi-analytical formula.The accuracy of this approach is verified by a one-dimensional model.Two 3-D models are used to demonstrate the effectiveness of the adaptive mesh refinement and the influences of dipping anisotropy layer on the marine CSEM responses for both inline and broadside geometries.The complex synthetic model is simulated to show the capability and flexibility of the approach for geometrically complex situations.

ADAPTIVE FINITE ELEMENT METHOD BASED ON OPTIMAL ERROR ESTIMATES FOR LINEAR ELLIPTIC PROBLEMS ON CONCAVE CORNER DOMAINS

An adaptive finite element procedure designed for specific computational goals is presented,using mesh refinement strategies based on optimal or nearly optimal a priori error estimates for the finite element method and using estimators of the local regularity of the unknown exact solution derived from computed approximate solutions.The proposed procedure is analyzed in detail for a non-trivial class of corner problems and shown to be efficient in the sense that the method can generate the correct type of refinements and lead to the desired control under consideration.

Quasi-optimal complexity of adaptive finite element method for linear elasticity problems in two dimensions

This paper introduces an adaptive finite element method （AFEM） using the newest vertex bisection and marking exclusively according to the error estimator without special treatment of oscillation. By the combination of the global lower bound and the localized upper bound of the posteriori error estimator, perturbation of oscillation, and cardinality of the marked element set, it is proved that the AFEM is quasi-optimal for linear elasticity problems in two dimensions, and this conclusion is verified by the numerical examples.

An Adaptive Finite Element Method Based on Optimal Error Estimates for Linear Elliptic Problems

The subject of the work is to propose a series of papers about adaptive finite element methods based on optimal error control estimate. This paper is the third part in a series of papers on adaptive finite element methods based on optimal error estimates for linear elliptic problems on the concave corner domains. In the preceding two papers (part 1:Adaptive finite element method based on optimal error estimate for linear elliptic problems on concave corner domain; part 2:Adaptive finite element method based on optimal error estimate for linear elliptic problems on nonconvex polygonal domains), we presented adaptive finite element methods based on the energy norm and the maximum norm. In this paper, an important result is presented and analyzed. The algorithm for error control in the energy norm and maximum norm in part 1 and part 2 in this series of papers is based on this result.

A Goal-Oriented Adaptive Finite Element Method for 3D Resistivity Modeling Using Dual-Error Weighting Approach

A goal-oriented adaptive finite element（FE） method for solving 3D direct current（DC） resistivity modeling problem is presented. The model domain is subdivided into unstructured tetrahedral elements that allow for efficient local mesh refinement and flexible description of complex models. The elements that affect the solution at each receiver location are adaptively refined according to a goal-oriented posteriori error estimator using dual-error weighting approach. The FE method with adapting mesh can easily handle such structures at almost any level of complexity. The method is demonstrated on two synthetic resistivity models with analytical solutions and available results from integral equation method, so the errors can be quantified. The applicability of the numerical method is illustrated on a resistivity model with a topographic ridge. Numerical examples show that this method is flexible and accurate for geometrically complex situations.

Adaptive hp finite element method for fluorescence molecular tomography with simplified spherical harmonics approximation

Recently,the simplified spherical harmonics equations(SP)model has at tracted much att entionin modeling the light propagation in small tissue ggeometriesat visible and near-infrared wave-leng ths.In this paper,we report an eficient numerical method for fluorescence moleeular tom-ography(FMT)that combines the advantage of SP model and adaptive hp finite elementmethod(hp-FEM).For purposes of comparison,hp-FEM and h-FEM are,respectively applied tothe reconstruction pro cess with diffusion approximation and SPs model.Simulation experiments on a 3D digital mouse atlas and physical experiments on a phantom are designed to evaluate thereconstruction methods in terms of the location and the reconstructed fluorescent yield.Theexperimental results demonstrate that hp-FEM with SPy model,yield more accurate results thanh-FEM with difusion approximation model does.The phantom experiments show the potentialand feasibility of the proposed approach in FMT applications.

Adaptive Finite Element Method Based on Optimal Error Estimtes for Linear Elliptic Problems on Concave Corner Domains (continuation)

This paper is the third part in a series of papers on adaptive finite elementmethods based on optimal error estimates for linear elliptic problems on the concavecorner domains. In this paper, a result is obtained. The algorithms for error controlboth in the energy norm and in the maximum norm presented in part 1 and part 2 ofthis series arc based on this result.

A Multi-Mesh Adaptive Finite Element Approximation to Phase Field Models

In this work,we propose an efficient multi-mesh adaptive finite element method for simulating the dendritic growth in two-and three-dimensions.The governing equations used are the phase field model,where the regularity behaviors of the relevant dependent variables,namely the thermal field function and the phase field function,can be very different.To enhance the computational efficiency,we approximate these variables on different h-adaptive meshes.The coupled terms in the system are calculated based on the implementation of the multi-mesh h-adaptive algorithm proposed by Li(J.Sci.Comput.,pp.321-341,24(2005)).It is illustrated numerically that the multi-mesh technique is useful in solving phase field models and can save storage and the CPU time significantly.

An Adaptive Finite Element PML Method for the Acoustic Scattering Problems in Layered Media

The paper concerns the numerical solution for the acoustic scattering problems in a two-layer medium.The perfectly matched layer(PML)technique is adopted to truncate the unbounded physical domain into a bounded computational domain.An a posteriori error estimate based adaptive finite element method is developed to solve the scattering problem.Numerical experiments are included to demonstrate the efficiency of the proposed method.

WEAK-DUALITY BASED ADAPTIVE FINITE ELEMENT METHODS FOR PDE-CONSTRAINED OPTIMIZATION WITH POINTWISE GRADIENT STATE-CONSTRAINTS

Adaptive finite element methods for optimization problems for second order linear el- liptic partial differential equations subject to pointwise constraints on the l2-norm of the gradient of the state are considered. In a weak duality setting, i.e. without assuming a constraint qualification such as the existence of a Slater point, residual based a posteriori error estimators are derived. To overcome the lack in constraint qualification on the continuous level, the weak Fenchel dual is utilized. Several numerical tests illustrate the performance of the proposed error estimators.Mathematics subject classification： 65N30, 90C46, 65N50, 49K20, 49N15, 65K10.

Adaptive mixed least squares Galerkin/Petrov finite element method for stationary conduction convection problems

An adaptive mixed least squares Galerkin/Petrov finite element method （FEM） is developed for stationary conduction convection problems. The mixed least squares Galerkin/Petrov FEM is consistent and stable for any combination of discrete velocity and pressure spaces without requiring the Babuska-Brezzi stability condition. Using the general theory of Verfiirth, the posteriori error estimates of the residual type are derived. Finally, numerical tests are presented to illustrate the effectiveness of the method.

Multi-mesh Adaptive Finite Element Algorithms for Constrained Optimal Control Problems Governed By Semi-Linear Parabolic Equations

In this paper, we derive a posteriori error estimators for the constrained optimal control problems governed by semi-linear parabolic equations under some assumptions. Then we use them to construct reliable and efficient multi-mesh adaptive finite element algorithms for the optimal control problems. Some numerical experiments are presented to illustrate the theoretical results.

Strategies for h-Adaptive Refinement for a Finite Element Treatment of Harmonic Oscillator Schrodinger Eigenproblem

A Schrodinger eigenvalue problem is solved for the 219 quantum simple harmonic oscillator using a finite element discretization of real space within which elements are adaptively spatially refined. We compare two competing methods of adaptively discretizing the real-space grid on which computations are performed without modifying the standard polynomial basis-set traditionally used in finite element interpolations; namely, （i） an application of the Kelly error estimator, and （ii） a refinement based on the local potential level. When the performance of these methods are compared to standard uniform global refinement, we find that they significantly improve the total time spent in the eigensolver.

ERROR REDUCTION IN ADAPTIVE FINITE ELEMENT APPROXIMATIONS OF ELLIPTIC OBSTACLE PROBLEMS

We consider an adaptive finite element method （AFEM） for obstacle problems associated with linear second order elliptic boundary value problems and prove a reduction in the energy norm of the discretization error which leads to R-linear convergence. This result is shown to hold up to a consistency error due to the extension of the discrete multipliers （point functionals） to H^-1 and a possible mismatch between the continuous and discrete coincidence and noncoincidence sets. The AFEM is based on a residual-type error estimator consisting of element and edge residuals. The a posteriori error analysis reveals that the significant difference to the unconstrained case lies in the fact that these residuals only have to be taken into account within the discrete noncoincidence set. The proof of the error reduction property uses the reliability and the discrete local efficiency of the estimator as well as a perturbed Galerkin orthogonality. Numerical results are given illustrating the performance of the AFEM.

Adaptive Finite Element Approximations for a Class of Nonlinear Eigenvalue Problems in Quantum Physics

In this paper,we study an adaptive finite element method for a class of nonlinear eigenvalue problems resulting from quantum physics that may have a nonconvex energy functional.We prove the convergence of adaptive finite element approximations and present several numerical examples of micro-structure of matter calculations that support our theory.

Retrieving Topological Information of Implicitly Represented Diffuse Interfaces with Adaptive Finite Element Discretization

We consider the finite element based computation of topological quantities of implicitly represented surfaces within a diffuse interface framework.Utilizing an adaptive finite element implementation with effective gradient recovery techniques,we discuss how the Euler number can be accurately computed directly from the numerically solved phase field functions or order parameters.Numerical examples and applications to the topological analysis of point clouds are also presented.

A HIGH ORDER ADAPTIVE FINITE ELEMENT METHOD FOR SOLVING NONLINEAR HYPERBOLIC CONSERVATION LAWS

In this note, we apply the h-adaptive streamline diffusion finite element method with a small mesh-dependent artificial viscosity to solve nonlinear hyperbolic partial differential equations, with the objective of achieving high order accuracy and mesh efficiency. We compute the numerical solution to a steady state Burgers equation and the solution to a converging-diverging nozzle problem. The computational results verify that, by suitably choosing the artificial viscosity coefficient and applying the adaptive strategy based on a posterior error estimate by Johnson et al., an order of N^-3/2 accuracy can be obtained when continuous piecewise linear elements are used, where N is the number of elements.