Accuracy analysis of immersed boundary method using method of manufactured solutions

The immersed boundary method is an effective technique for modeling and simulating fluid-structure interactions especially in the area of biomechanics.This paper analyzes the accuracy of the immersed boundary method.The procedure contains two parts,i.e.,the code verification and the accuracy analysis.The code verification provides the confidence that the code used is free of mistakes,and the accuracy analysis gives the order of accuracy of the immersed boundary method.The method of manufactured solutions is taken as a means for both parts.In the first part,the numerical code employs a second-order discretization scheme,i.e.,it has second-order accuracy in theory.It matches the calculated order of accuracy obtained in the numerical calculation for all variables.This means that the code contains no mistake,which is a premise of the subsequent work.The second part introduces a jump in the manufactured solution for the pressure and adds the corresponding singular forcing terms in the momentum equations.By analyzing the discretization errors,the accuracy of the immersed boundary method is proven to be first order even though the discretization scheme is second order.It has been found that the coarser mesh may not be sensitive enough to capture the influence of the immersed boundary,and the refinement on the Lagrangian markers barely has any effect on the numerical calculation.

Verification of an Explicitly Coupled Thermal-Phase Field-Mechanical Electromagnetic (TPME) Framework by the Method of Manufactured Solutions

An explicitly coupled two-dimensional (2D) multiphysics finite element method (FEM) framework comprised of thermal, phase field, mechanical and electromagnetic (TPME) equations was developed to simulate the conversion of solid kerogen in oil shale to liquid oil through in-situ pyrolysis by radio frequency heating. Radio frequency heating as a method of in-situ pyrolysis represents a tenable enhanced oil recovery method, whereby an applied electrical potential difference across a target oil shale formation is converted to thermal energy, heating the oil shale and causing it to liquify to become liquid oil. A number of in-situ pyrolysis methods are reviewed but the focus of this work is on the verification of the TPME numerical framework to model radio frequency heating as a potential dielectric heating process for enhanced oil recovery. Very few studies exist which describe production from oil shale;furthermore, there are none that specifically address the verification of numerical models describing radio frequency heating. As a result, the Method of Manufactured Solutions (MMS) was used as an analytical verification method of the developed numerical code. Results show that the multiphysics finite element framework was adequately modeled enabling the simulation of kerogen conversion to oil as a part of the analysis of a TPME numerical model.

Effect of regularized delta function on accuracy of immersed boundary method

The immersed boundary method is an effective technique for modeling and simulating fluid-structure interactions especially in the area of biomechanics. The effect of the regularized delta function on the accuracy is an important subject in the property study. A method of manufactured solutions is used in the research. The computational code is first verified to be mistake-free by using smooth manufactured solutions. Then, a jump in the manufactured solution for pressure is introduced to study the accuracy of the immersed boundary method. Four kinds of regularized delta functions are used to test the effect on the accuracy analysis. By analyzing the discretization errors, the accuracy of the immersed boundary method is proved to be first-order. The results show that the regularized delta function cannot improve the accuracy, but it can change the discretization errors in the entire computational domain.

Accuracy analysis of gradient reconstruction on isotropic unstructured meshes and its effects on inviscid flow simulation

The accuracy of gradient reconstruction methods on unstructured meshes is analyzed both mathematically and numerically.Mathematical derivations reveal that,for gradient reconstruction based on the Green-Gauss theorem(the GG methods),if the summation of first-and-lower-order terms does not counterbalance in the discretized integral process,which rarely occurs,second-order accurate approximation of face midpoint value is necessary to produce at least first-order accurate gradient.However,gradient reconstruction based on the least-squares approach(the LSQ methods)is at least first-order on arbitrary unstructured grids.Verifications are performed on typical isotropic grid stencils by analyzing the relationship between the discretization error of gradient reconstruction and the discretization error of the face midpoint value approximation of a given analytic function.Meanwhile,the numerical accuracy of gradient reconstruction methods is examined with grid convergence study on typical isotropic grids.Results verify the phenomenon of accuracy degradation for the GG methods when the face midpoint value condition is not satisfied.The LSQ methods are proved to be at least first-order on all tested isotropic grids.To study gradient accuracy effects on inviscid flow simulation,solution errors are quantified using the Method of Manufactured Solutions(MMS)which was validated before adoption by comparing with an exact solution case,i.e.,the 2-dimensional(2D)inviscid isentropic vortex.Numerical results demonstrate that the order of accuracy(OOA)of gradient reconstruction is crucial in determining the OOA of numerical solutions.Solution accuracy deteriorates seriously if gradient reconstruction does not reach first-order.