UPML formulation for truncating conductive media in curvilinear coordinates

This work presents a formulation based on UPML for truncating conductive media by using a local and non-orthogonal coordinate system to solve Maxwell’s equations by the FDTD method. The detailed procedure for obtaining the UPML equations for this case is shown and the complete equation set is provided.     

The weak solution of an antiplane contact problem for electro-viscoelastic materials with long-term memory

We study a mathematical model which describes the antiplane shear deformation of a cylinder in frictionless contact with a rigid foundation. The material is assumed to be electro-viscoelastic with long-term memory, and the friction is modeled with Tresca’s law and the foundation is assumed to be electrically conductive. First we derive the classical variational formulation of the model which is given by a system coupling an evolutionary variational equality for the displacement field with a time-dependent variational equation for the potential field. Then we prove the existence of a unique weak solution to the model. Moreover, the proof is based on arguments of evolution equations and on the Banach fixedpoint theorem.     

Numerical Solution of Newtonian and Non-Newtonian Flows

This paper deals with the numerical solution of Newtonian and non-Newtonian flows. The flows are supposed to be laminar, viscous, incompressible and steady. The model used for non-Newtonian fluids is some variant of power-law. Governing equations in this model are incompressible Navier-Stokes equations. For numerical solution one could use artificial compressibility method with three stage Runge-Kutta finite volume method in cell centered formulation for discretization of space derivatives. Following cases of flows are solwed: flow through a bypass connected to main channel in 2D and 3D and non-Newtonian flow through branching channels in 2D. These results are presented for 2D and 3D case. This problem could have an application in the area of biomedicine. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exterior problem of magnetic field in Darwin model and its numerical solution

In this paper, we consider the solutions of magnetic field in the Darwin model to the Maxwell's equations in 2D unbounded domain. We first deduce the variational formulation and prove the well‐posedness of the weak solution, and then prove the existence and uniqueness of the infinite element solution. Error estimate and the numerical examples are provided. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Differential transformation approach to thermal conductive problems with discontinuous boundary condition

A two-dimensional differential transformation method is employed to reduce the partial differential equations of the non-continuous thermal conductive boundary value problem to a Taylor series in a polynomial form. The partial differential equations are solved by the two-dimensional T-spectra method of differential transformation and by the use of trial initial polynomial conditions. The investigative parameters include the time-step of the differential transformation and the number of sub-domain segments. The numerical simulation results indicate that the proposed approach using the two-dimensional T-spectra method of differential transformation is applicable to the solution of non-continuous thermal conductive boundary value problems.     

Mathematical modelling of compressible viscous fluid flow in a male rotor‐housing gap of screw machines

This contribution is devoted to the mathematical modelling of a compressible viscous fluid flow through a 2‐D model of the male rotor‐housing gap in screw machines. Numerical solution of the nonlinear conservative system of the compressible Navier‐Stokes equations is obtained by means of the cell‐centred finite volume formulation of the explicit two‐step TVD MacCormack scheme proposed by Causon on a structured quadrilateral grid using the own developed numerical code. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stochastic systems in Riemannian manifolds

A formulation of stochastic systems in a Riemannian manifold is given by stochastic differential equations in the tangent bundle of the manifold. Brownian motion is constructed in a compact Riemannian manifold as well as the horizontal lift of this process to the bundle of orthonormal frames. The solution of some stochastic differential equations in the tangent bundle of the manifold is defined by the transformation of the measure for the manifold-valued Brownian motion by a suitable Radon-Nikodym derivative. Real-valued stochastic integrals are defined for this Brownian motion using parallelism along the Brownian paths. A stochastic control problem is formulated and solved for these stochastic systems where a suitable convexity condition is assumed.This research was supported by NSF Grants Nos. GK-32136, ENG-75-06562, and MCS-76-01695.The author wishes to thank D. Gromoll, J. Simons, and J. Thorpe for some helpful conversations on differential geometry.     

粘弹性方程全离散化有限体积元格式及数值模拟

本文利用有限体积元方法研究二维粘弹性方程, 给出一种时间二阶精度的全离散化有限体积元格式, 并给出这种全离散化有限体积元解的误差估计, 最后用数值例子验证数值结果与理论结果是相吻合的. 通过与有限元方法和有限差分方法相比较, 进一步说明了全离散化有限体积元格式是求解二维粘弹性方程数值解的最有效方法之一.     

A C0-WEAK GALERKIN FINITE ELEMENT METHOD FOR THE TWO-DIMENSIONAL NAVIER-STOKES EQUATIONS IN STREAM-FUNCTION FORMULATION

We propose and analyze a $C^0$-weak Galerkin (WG) finite element method for the numerical solution of the Navier-Stokes equations governing 2D stationary incompressible flows. Using a stream-function formulation, the system of Navier-Stokes equations is reduced to a single fourth-order nonlinear partial differential equation and the incompressibility constraint is automatically satisfied. The proposed method uses continuous piecewise-polynomial approximations of degree $k+2$ for the stream-function $\psi$ and discontinuous piecewise-polynomial approximations of degree $k+1$ for the trace of $\nabla\psi$ on the interelement boundaries. The existence of a discrete solution is proved by means of a topological degree argument, while the uniqueness is obtained under a data smallness condition. An optimal error estimate is obtained in $L^2$-norm, $H^1$-norm and broken $H^2$-norm. Numerical tests are presented to demonstrate the theoretical results.     

Modeling and analysis of the unilateral contact of a piezoelectric body with a conductive support

We consider a mathematical model which describes the quasistatic process of contact between a piezoelectric body and an electrically conductive support, the so-called foundation. We model the material's behavior with a nonlinear electro-viscoelastic constitutive law; the contact is frictionless and is described with the Signorini condition and a regularized electrical conductivity condition. We derive a variational formulation for the problem and then we prove the existence of a unique weak solution to the model. The proof is based on arguments of nonlinear equations with multivalued maximal monotone operators and fixed point. Then we introduce a fully discrete scheme, based on the finite element method to approximate the spatial variable and the backward Euler scheme to discretize the time derivatives. We treat the unilateral contact conditions by using an augmented Lagrangian approach. We implement this scheme in a numerical code then we present numerical simulations in the study of two-dimensional test problems, together with various comments and interpretations.     

基于拟Legendre多项式求解一类分数阶微分方程

本文基于移位的Legendre多项式构造一类新的正交拟Legendre多项式求解一类分数阶微分方程.用阶数随所求未知函数的微分的阶数而变化的拟Legendre多项式逼近未知函数;利用分数阶积分的性质推导拟Legendre多项式的积分算子阵,结合算子矩阵的思想和Tau方法,将问题转化为求解代数方程组的问题.最后,给出数值算例证明该方法的有效性.     

Fractional-order Legendre functions for solving fractional-order differential equations

In this article, a general formulation for the fractional-order Legendre functions (FLFs) is constructed to obtain the solution of the fractional-order differential equations. Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. Therefore, an efficient and reliable technique for the solution of them is too important. For the concept of fractional derivative we will adopt Caputo’s definition by using Riemann–Liouville fractional integral operator. Our main aim is to generalize the new orthogonal functions based on Legendre polynomials to the fractional calculus. Also a general formulation for FLFs fractional derivatives and product operational matrices is driven. These matrices together with the Tau method are then utilized to reduce the solution of this problem to the solution of a system of algebraic equations. The method is applied to solve linear and nonlinear fractional differential equations. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

The use of Chebyshev Polynomials in the space-time least-squares spectral element method

Chebyshev polynomials of the first kind are employed in a space-time least-squares spectral element formulation applied to linear and nonlinear hyperbolic scalar equations. No stabilization techniques are required to render a stable, high order accurate scheme. In parts of the domain where the underlying exact solution is smooth, the scheme exhibits exponential convergence with polynomial enrichment, whereas in parts of the domain where the underlying exact solution contains discontinuities the solution displays a Gibbs-like behavior. An edge detection method is employed to determine the position of the discontinuity. Piecewise reconstruction of the numerical solution retrieves a monotone solution. Numerical results will be given in which the capabilities of the space-time formulation to capture discontinuities will be demonstrated.     

Analytical solution for the fully coupled static response of variable stiffness composite beams

An analytical solution is presented for the 3D static response of variable stiffness non-uniform composite beams. Based on Euler-Bernoulli theory, a set of governing differential equations are obtained, in which four degrees of freedom are fully coupled. For the variable stiffness beam, the governing field equations have variable coefficients reflecting the stiffness variation along the beam. Using the direct integration technique, the general analytical solution is derived in the integral form and the closed-form expressions of the obtained solutions are presented employing a series expansion approximation. The series expansion representation enables the proposed approach to be applicable for variable stiffness composite beams with arbitrary span-wise variation of properties. As an alternative solution, the Chebyshev collocation method is applied to the proposed formulation to verify the results obtained from the analytical solution. A number of variable stiffness composite beams made by fibre steering with various boundary conditions and stacking sequences are considered as the test cases. The static response are presented based on the analytical solution and Chebyshev collocation method and excellent agreement is observed for all test cases. The proposed model presents a reliable and efficient approach for capturing the complicated behaviour of variable stiffness non-uniform composite beams.     

A perturbation-boundary integral equation method for transient dynamics of nonlinear elastic materials

This study presents a formulation for solving the transient dynamics of nonlinear elastic materials. By using a perturbation expansion to linearize the basic equations and applying the Laplace transform to the subsequent perturbation equations, the boundary value problem of the transformed equations is further reduced to various boundary integral equations. After discretization of the integral equations, these are solved numerically, completing the solution in the Laplace transform space. Performing a numerical inversion of the Laplace transform yields the solution of the problem in the time domain.     

A boundary element procedure for transient wave propagations in two-dimensional isotropic elastic media

An integral formulation of the elastodynamic equations is presented and discretized to develope a numerical solution procedure. Constant space and linear time dependent interpolation functions are implemented. The expression obtained, the boundary element equations can be solved using a time-stepping scheme. An example is given for comparison with known solutions.     

Strong solutions to the equations of electrically conductive magnetic fluids

We study the equations of flow of an electrically conductive magnetic fluid, when the fluid is subjected to the action of an external applied magnetic field. The system is formed by the incompressible Navier–Stokes equations, the magnetization relaxation equation of Bloch type and the magnetic induction equation. The system takes into account the Kelvin and Lorentz force densities. We prove the local-in-time existence of the unique strong solution to the system equipped with initial and boundary conditions. We also establish a blow-up criterion for the local strong solution.     

A modified Chan-Vese model and its theoretical proof

We present a modified Chan-Vese functional and give its theoretical proof. By using the geometric heat flow method to all the Euler-Lagrange equations, a system of evolution equations in level set formulation is derived. We study the existence of solution to this system by Schauder fixed point theorem and the implicit function theorem in Banach space. This variational formulation can detect interior and exterior boundaries of desired object(s) in color images.     

Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects

An analysis solution method (ASM) is proposed for analyzing arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystal (QC) media. The extended displacement discontinuity (EDD) boundary integral equations governing three-dimensional (3D) crack problems are transferred to simplified integral-differential forms by introducing some complex quantities. The proposed ASM is based on the analogy between these EDD boundary equations for 3D planar cracks problems of 2D hexagonal QCs and those in isotropic thermoelastic materials. Mixed model crack problems under combined normal, tangential and thermal loadings are considered in 2D hexagonal QC media. By virtue of ASM, the solutions to 3D planar crack problems under various types of loadings for 2D hexagonal QCs are formulated through comparison to the corresponding solutions of isotropic thermoelastic materials which have been studied intensively and extensively. As an application, analytical solutions of a penny-shaped crack subjected uniform distributed combined loadings are obtained. Especially, the analytical solutions to a penny-shaped crack subjected to the anti-symmetric uniform thermal loading are first derived for 2D hexagonal QCs. Numerical solutions obtained by EDD boundary element method provide a way to verify the validity of the presented formulation. The influences of phonon-phason coupling effect on fracture parameters of 2D hexagonal QCs are assessed.     

Iterative solvers and stabilisation for mixed electrostatic and magnetostatic formulations

Mixed electrostatic and magnetostatic finite element formulations are considered. Solution methods for the resulting indefinite algebraic systems are investigated. Methods developed for the mixed formulations of the Stokes equations are modified in order to apply to the Maxwell equations: an efficient block preconditioner is proposed and a stabilised formulation is described. The different methods are applied to 2D and 3D examples.