Finite-difference time-domain modeling of curved material interfaces by using boundary condition equations method

To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.

Improved absorbing boundary condition based on linear interpolation for ADI-FDTD method

With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.

Application of inverse method to estimation of boundary conditions during investment casting simulation

Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient （IHTC） between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.

Implicit finite difference method for fractional percolation equation with Dirichlet and fractional boundary conditions

An implicit finite difference method is developed for a one-dimensional frac- tional percolation equation （FPE） with the Dirichlet and fractional boundary conditions. The stability and convergence are discussed for two special cases, i.e., a continued seep- age flow with a monotone percolation coefficient and a seepage flow with the fractional Neumann boundary condition. The accuracy and efficiency of the method are checked with two numerical examples.

A Comparison Study of Deep Galerkin Method and Deep Ritz Method for Elliptic Problems with Different Boundary Conditions

Recent years have witnessed growing interests in solving partial differential equations by deep neural networks,especially in the high-dimensional case.Unlike classical numerical methods,such as finite difference method and finite element method,the enforcement of boundary conditions in deep neural networks is highly nontrivial.One general strategy is to use the penalty method.In the work,we conduct a comparison study for elliptic problems with four different boundary conditions,i.e.,Dirichlet,Neumann,Robin,and periodic boundary conditions,using two representative methods:deep Galerkin method and deep Ritz method.In the former,the PDE residual is minimized in the least-squares sense while the corresponding variational problem is minimized in the latter.Therefore,it is reasonably expected that deep Galerkin method works better for smooth solutions while deep Ritz method works better for low-regularity solutions.However,by a number of examples,we observe that deep Ritz method can outperform deep Galerkin method with a clear dependence of dimensionality even for smooth solutions and deep Galerkin method can also outperform deep Ritz method for low-regularity solutions.Besides,in some cases,when the boundary condition can be implemented in an exact manner,we find that such a strategy not only provides a better approximate solution but also facilitates the training process.

A New Extrapolation Treatment for Boundary Conditions in Lattice Boltzmann Method

The lattice Boltzmann method is widely used in solving hydrodynamics for various systems.The accuracy of the extrapolation method for boundary conditions and Mei’s method for curved boundary is not high enough.We propose a modified extrapolation method for boundary conditions.Numerical simulations conform the higher accuracy and stability compared with the other methods.

Spline Fictitious Boundary Element Alternating Method for Edge Crack Problems with Mixed Boundary Conditions

The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.

ON FICTITIOUS DOMAIN METHOD FOR THE NUMERICAL SOLUTION TO HEAT CONDUCTION EQUATION WITH DERIVATIVE BOUNDARY CONDITIONS

This paper investigates the stability and convergence of some knowndifference schemes for the numerical solution to heat conduction equation withderivative boundary conditions by the fictitious domain method.The discrete vari-ables at the false mesh points are firstly eliminated from the difference schemes andthe local truncation errors are then analyzed in detail.The stability and convergenceof the schemes are proved by energy method.An improvement is proposed to obtainbetter schemes over the original ones.Several numerical examples and comparisonswith other schemes are presented.

On L ∞ Stability and Convergence of Fictitious Domain Method for the Numerical Solution to Parabolic Differential Equation with Derivative Boundary Conditions

This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in L ∞ are proven.

ANALYSIS OF OPEN MICROSTRIP STRUCTURES BY USING DIAKOPTIC METHOD OF LINES COMBINED WITH PERIODIC BOUNDARY CONDITIONS

This paper presents the analysis of open microstrip structures by using diakoptic method of lines (ML) combined with periodic boundary conditions (PBC). The parameters of microstrip patch are obtained from patch current excited by plane wave. Impedance matrix elements are computed by using fast Fourier transform(FFT), and reduced equation is solved by using diakoptic technique. Consequently, the computing time is reduced significantly. The convergence property of simulating open structure by using PBC and the comparison of the computer time between using PBC and usual absorbing boundary condition (ABC) show the validity of the method proposed in this paper. Finally, the resonant frequency of a microstrip patch is computed. The numerical results obtained are in good agreement with those published.

AN IMPROVEMENT FOR LOWER－ORDER PANEL METHOD BASED ON THE DIRICHLET BOUNDARY CONDITION

An improved algorithm for velocity field of general configurations ispresentd for low-order panel method based on the internal Dirichlet boundary condi-tion. A direct calculating method for the velocity distribution by means of a limit pro-cess combining with analytic evaluation of higher-order singular integrals instead of theconventional method of doublet strength gradient is devised in order to avoid the diffi-culty of edge extrapolation of doublet strength. The problem of substantialunderpredictions of the induced drag coefficient obtained from the VSAERO analysisdisappears for the present improved algorithm. Illustrative calculations for several testcases such as swept back wing, swept forward wing and wing-body combination showthat the accuracy of results may be improved and is competitive with high-order panelmethod. In addition, the present direct integral method can be used to evaluate the ve-locity distribution for external flow field correctly, where the method of gradient cannot be used at all.

An extended iterative direct-forcing immersed boundary method in thermo-fluid problems with Dirichlet or Neumann boundary conditions

An iterative direct-forcing immersed boundary method is extended and used to solve convection heat transfer problems.The pressure,momentum source,and heat source at immersed boundary points are calculated simultaneously to achieve the best coupling.Solutions of convection heat transfer problems with both Dirichlet and Neumann boundary conditions are presented.Two approaches for the implementation of Neumann boundary condition,i.e.direct and indirect methods,are introduced and compared in terms of accuracy and computational efficiency.Validation test cases include forced convection on a heated cylinder in an unbounded flow field and mixed convection around a circular body in a lid-driven cavity.Furthermore,the proposed method is applied to study the mixed convection around a heated rotating cylinder in a square enclosure with both iso-heat flux and iso-thermal boundary conditions.Computational results show that the order of accuracy of the indirect method is less than the direct method.However,the indirect method takes less computational time both in terms of the implementation of the boundary condition and the post processing time required to compute the heat transfer variables such as the Nusselt number.It is concluded that the iterative direct-forcing immersed boundary method is a powerful technique for the solution of convection heat transfer problems with stationary/moving boundaries and various boundary conditions.

Investigation of nanofluid flow in a vertical channel considering polynomial boundary conditions by Akbari-Ganji’s method

In this research,a vertical channel containing a laminar and fully developed nanofluid flow is investigated.The channel surface’s boundary conditions for temperature and volume fraction functions are considered qth-order polynomials.The equations related to this problem have been extracted and then solved by the AGM and validated through the Runge-Kutta numerical method and another similar study.In the study,the effect of parameters,including Grashof number,Brownian motion parameter,etc.,on the motion,velocity,temperature,and volume fraction of nanofluids have been analyzed.The results demonstrate that increasing the Gr number by 100%will increase the velocity profile function by 78%and decrease the temperature and fraction profiles by 20.87%and 120.75%.Moreover,rising the Brownian motion parameter in five different sizes(0.1,0.2,0.3,0.4,and 0.5)causes lesser velocity,about 24.3%at first and 4.35%at the last level,and a maximum 52.86%increase for temperature and a 24.32%rise for ψ occurs when N b rises from 0.1 to 0.2.For all N_(t) values,at least 55.44%,18.69%,for F(η),andΩ(η),and 20.23%rise for ψ(η)function is observed.Furthermore,enlarging the N r parameter from 0.25 to 0.1 leads F(η)to rise by 199.7%,fluid dimensionless temperature,and dimensional volume fraction to decrease by 18%and 92.3%.In the end,a greater value of q means a more powerful energy source,amplifying all velocity,temperature,and volume fraction functions.The main novelty of this research is the combined convection qth-order polynomials boundary condition applied to the channel walls.Moreover,The AMG semi-analytical method is used as a novel method to solve the governing equations.

A simple and robust method for essential boundary condition treatments in mesh-less methods

To improve the treatment efficiency of essential boundary condition in mesh-less methods, a simple and robust method is proposed in this paper. Rising weight of nodes in the construction of trail function, specified for essential boundary condition, can make the trail function pass through these nodes. And then, the trail function can satisfy the essential boundary condition previously by setting diagonal element to 1 or multiplying diagonal element by a big number in FEM. The MLS method is adopted to validate this method, and it is proved that this method is eostless and robust in most of mesh-less methods.

The Algebraic Immersed Interface and Boundary Method for Elliptic Equations with Jump Conditions

A new simple fictitious domain method, the algebraic immersed interface and boundary (AIIB) method, is presented for elliptic equations with immersed interface conditions. This method allows jump conditions on immersed interfaces to be discretized with a good accuracy on a compact stencil. Auxiliary unknowns are created at existing grid locations to increase the degrees of freedom of the initial problem. These auxiliary unknowns allow imposing various constraints to the system on interfaces of complex shapes. For instance, the method is able to deal with immersed interfaces for elliptic equations with jump conditions on the solution or discontinuous coefficients with a second order of spatial accuracy. As the AIIB method acts on an algebraic level and only changes the problem matrix, no particular attention to the initial discretization is required. The method can be easily implemented in any structured grid code and can deal with immersed boundary problems too. Several validation problems are presented to demonstrate the interest and accuracy of the method.

Free Vibration Analysis of Rectangular Plate with Cutouts under Elastic Boundary Conditions in Independent Coordinate Coupling Method

Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.

Highly Efficient Method for Solving Parabolic PDE with Nonlocal Boundary Conditions

In this work, a highly efficient algorithm is developed for solving the parabolic partial differential equation (PDE) with the nonlocal condition. For this purpose, we employ orthogonal Chelyshkov polynomials as the basis. The convergence analysis of the proposed scheme is derived. Numerical experiments are carried out to explain the efficiency and precision of the proposed scheme. Furthermore, the reliability of the scheme is verified by comparisons with assured existing methods.

A Split-Step Predictor-Corrector Method for Space-Fractional Reaction-Diffusion Equations with Nonhomogeneous Boundary Conditions

A split-step second-order predictor-corrector method for space-fractional reaction-diffusion equations with nonhomogeneous boundary conditions is presented and analyzed for the stability and convergence.The matrix transfer technique is used for spatial discretization of the problem.The method is shown to be unconditionally stable and second-order convergent.Numerical experiments are performed to confirm the stability and secondorder convergence of the method.The split-step predictor-corrector method is also compared with an IMEX predictor-corrector method which is found to incur oscillatory behavior for some time steps.Our method is seen to produce reliable and oscillatioresults for any time step when implemented on numerical examples with nonsmooth initial data.We also present a priori reliability constraint for the IMEX predictor-corrector method to avoid unwanted oscillations and show its validity numerically.

Application of the double absorbing boundary condition in seismic modeling

We apply the newly proposed double absorbing boundary condition（DABC）（Hagstrom et al., 2014） to solve the boundary reflection problem in seismic finite-difference（FD） modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer（PML）, the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.

ON SOLUTION TO THE NAVIER-STOKES EQUATIONS WITH NAVIER SLIP BOUNDARY CONDITION FOR THREE DIMENSIONAL INCOMPRESSIBLE FLUID

In this article, we prove the existence and uniqueness of solutions of the NavierStokes equations with Navier slip boundary condition for incompressible fluid in a bounded domain of R^3. The results are established by the Galerkin approximation method and improved the existing results.