Analysis of radial basis function interpolation approach

The radial basis function （RBF） interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical properties in the laboratory on the basis of physical rock datasets, which include the formation factor, viscosity, permeability, and molecular composition. However, this approach does not consider the effect of spatial distribution of the calibration data on the interpolation result. This study proposes a new RBF interpolation approach based on the Freedman＇s RBF interpolation approach, by which the unit basis functions are uniformly populated in the space domain. The inverse results of the two approaches are comparatively analyzed by using our datasets. We determine that although the interpolation effects of the two approaches are equivalent, the new approach is more flexible and beneficial for reducing the number of basis functions when the database is large, resulting in simplification of the interpolation function expression. However, the predicted results of the central data are not sufficiently satisfied when the data clusters are far apart.

HERMITE—BIRKHOFF INTERPOLATION OF SCATTERED DATA BY RADIAL BASIS FUNCTIONS

For Hermite-Birkhoff interpolation of scattered multidumensional data by radial basis function (?),existence and characterization theorems and a variational principle are proved. Examples include (?)(r)=r^b,Duchon's thin-plate splines,Hardy's multiquadrics,and inverse multiquadrics.

Radial Basis Function Interpolation and Galerkin Projection for Direct Trajectory Optimization and Costate Estimation

This work presents a novel approach combining radial basis function(RBF)interpolation with Galerkin projection to efficiently solve general optimal control problems.The goal is to develop a highly flexible solution to optimal control problems,especially nonsmooth problems involving discontinuities,while accounting for trajectory accuracy and computational efficiency simultaneously.The proposed solution,called the RBF-Galerkin method,offers a highly flexible framework for direct transcription by using any interpolant functions from the broad class of global RBFs and any arbitrary discretization points that do not necessarily need to be on a mesh of points.The RBF-Galerkin costate mapping theorem is developed that describes an exact equivalency between the Karush-Kuhn-Tucker(KKT)conditions of the nonlinear programming problem resulted from the RBF-Galerkin method and the discretized form of the first-order necessary conditions of the optimal control problem,if a set of discrete conditions holds.The efficacy of the proposed method along with the accuracy of the RBF-Galerkin costate mapping theorem is confirmed against an analytical solution for a bang-bang optimal control problem.In addition,the proposed approach is compared against both local and global polynomial methods for a robot motion planning problem to verify its accuracy and computational efficiency.

Large Scattered Data Fitting Based on Radial Basis Functions

Solving large radial basis function （RBF） interpolation problem with non-customized methods is computationally expensive and the matrices that occur are typically badly conditioned. In order to avoid these difficulties, we present a fitting based on radial basis functions satisfying side conditions by least squares, although compared with interpolation the method loses some accuracy, it reduces the computational cost largely. Since the fitting accuracy and the non-singularity of coefficient matrix in normal equation are relevant to the uniformity of chosen centers of the fitted RBE we present a choice method of uniform centers. Numerical results confirm the fitting efficiency.

An efficient large-scale mesh deformation method based on MPI/Open MP hybrid parallel radial basis function interpolation

An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation,the present method employs different surface points sets and influence radius for each volume point movement,accompanied with efficient geometry searching strategy.The deformed surface points,also called Control Points (CPs),are stored in each processor.The displacement of spatial points is interpolated by using only 20-50 nearest control points,and the local influence radius is set to 5-20 times the maximum displacement of control points.To shorten the searching time for the nearest control point clouds,an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique.Besides,an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications.Three 3D cases,including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements,are used to test the present mesh deformation method.The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 450 and more than 80% parallel efficiency with 1024 MPI processors.In addition,the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.

3D anisotropic modeling and identification for airborne EM systems based on the spectral-element method

The airborne electromagnetic （AEM） method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element （SE） method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell＇s equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology.

Analysis of a Gyroscope′s Rotor Nonlinear Supported Magnetic Field Based on the B-Spline Wavelet-FEM

A supported framework of a gyroscope＇s rotor is designed and the B-Spline wavelet finite element model of nonlinear supported magnetic field is worked out. A new finite element space is studied in which the scaling function of the B-spline wavelet is considered as the shape function of a tetrahedton. The magnetic field is spited by an artificial absorbing body which used the condition of field radiating, so the solution is unique. The resolution is improved via the varying gradient of the B-spline function under the condition of unchanging gridding. So there are some advantages in dealing with the focus flux and a high varying gradient result from a nonlinear magnetic field. The result is more practical. Plots of flux and in the space is studied via simulating the supported system model. The results of the study are useful in the research of the supported magnetic system for the gyroscope rotor.

The Formulation of Finite Difference RBFWENO Schemes for Hyperbolic Conservation Laws:An Alternative Technique

To solve conservation laws,efficient schemes such as essentially nonoscillatory(ENO)and weighted ENO(WENO)have been introduced to control the Gibbs oscillations.Based on radial basis functions(RBFs)with the classical WENO-JS weights,a new type of WENO schemes has been proposed to solve conservation laws[J.Guo et al.,J.Sci.Comput.,70(2017),pp.551–575].The purpose of this paper is to introduce a new formulation of conservative finite difference RBFWENO schemes to solve conservation laws.Unlike the usual method for reconstructing the flux functions,the flux function is generated directly with the conservative variables.Comparing with Guo and Jung(2017),the main advantage of this framework is that arbitrary monotone fluxes can be employed,while in Guo and Jung(2017)only smooth flux splitting can be used to reconstruct flux functions.Several 1D and 2D benchmark problems are prepared to demonstrate the good performance of the new scheme.

Numerical simulations of vortex-induced vibrations of a flexible riser with different aspect ratiosin uniform and shear currents

This paper aimed at describing numerical simulations of vortex-induced vibrations（VIVs） of a long flexible riser with different length-to-diameter ratio（aspect ratio） in uniform and shear currents. Three aspect ratios were simulated： L/D= 500, 750 and 1 000. The simulation was carried out by the in-house computational fluid dynamics（CFD） solver viv-FOAM-SJTU developed by the authors, which was coupled with the strip method and developed on the OpenFOAM platform. Moreover, the radial basis function（RBF） dynamic grid technique is applied to the viv-FOAM-SJTU solver to simulate the VIV in both in-line（IL） and cross-flow（CF） directions of flexible riser with high aspect ratio. The validation of the benchmark case has been completed. With the same parameters, the aspect ratio shows a significant influence on VIV of a long flexible riser. The increase of aspect ratio exerted a strong effect on the IL equilibrium position of the riser while producing little effect on the curvature of riser. With the aspect ratio rose from 500 to 1 000, the maximum IL mean displacement increased from 3 times the diameter to 8 times the diameter. On the other hand, the vibration mode of the riser would increase with the increase of aspect ratio. When the aspect ratio was 500, the CF vibration was shown as a standing wave with a 3-（rd） order single mode. When the aspect ratio was 1 000, the modal weights of the 5-（th） and 6-（th） modes are high, serving as the dominant modes. The effect of the flow profile on the oscillating mode becomes more and more apparent when the aspect ratio is high, and the dominant mode of riser in shear flow is usually higher than that in uniform flow. When the aspect ratio was 750, the CF oscillations in both uniform flow and shear flow showed multi-mode vibration of the 4-（th） and 5-（th） mode. While, the dominant mode in uniform flow is the 4-（th） order, and the dominant mode in shear flow is the 5-（th） order.

RBFs-MSA Hybrid Method for Mesh Deformation

Simulating unsteady flow phenomena involving moving boundaries is a challenging task,one key requirement of which is a reliable and fast algorithm to deform the computational mesh.Radial basis functions（RBFs） interpolation is a very simple and robust method to deform the mesh.However,the number of operations and the requirement of memory storage will be increased rapidly as the number of grid nodes increases,which limits the application of RBFs to three-dimensional（3D） moving mesh.Moving submesh approach（MSA） is an efficient method,but its robustness depends on the method used to deform the background mesh.A hybrid method which combines the benefits of MSA and RBFs interpolation,which is called RBFs-MSA,has been presented.This hybrid method is proved to be robust and efficient via several numerical examples.From the aspect of the quality of deforming meshes,this hybrid method is comparable with the RBFs interpolation;from the aspect of computing efficiency,one test case shows that RBFs-MSA is about two orders of magnitude faster than RBFs interpolation.For these benefits of RBFs-MSA,the new method is suitable for unsteady flow simulation which refers to boundaries movement.

Numerical studies of static aeroelastic effects on grid fin aerodynamic performances

The grid fin is an unconventional control surface used on missiles and rockets. Although aerodynamics of grid fin has been studied by many researchers, few considers the aeroelastic effects.In this paper, the static aeroelastic simulations are performed by the coupled viscous computational fluid dynamics with structural flexibility method in transonic and supersonic regimes. The developed coupling strategy including fluid–structure interpolation and volume mesh motion schemes is based on radial basis functions. Results are presented for a vertical and a horizontal grid fin mounted on a body. Horizontal fin results show that the deformed fin is swept backward and the axial force is increased. The deformations also induce the movement of center of pressure, causing the reduction and reversal in hinge moment for the transonic flow and the supersonic flow,respectively. For the vertical fin, the local effective incidences are increased due to the deformations so that the deformed normal force is greater than the original one. At high angles of attack, both the deformed and original normal forces experience a sudden reduction due to the interference of leeward separated vortices on the fin. Additionally, the increment in axial force is shown to correlate strongly with the increment in the square of normal force.