COMPACT FINITE DIFFERENCE-FOURIER SPECTRAL METHOD FOR THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

High accuracy compact finite difference methods and their applications

Numerical simulation of complex flow fields with multi-scale structures is one of the most important and challenging branches of computational fluid dynamics. From linear analysis and numerical experiments it has been discovered that the higher-order accurate method can give reliable and efficient computational results, as well as better resolution of the complex flow fields with multi-scale structures. Compact finite difference schemes, which feature higher-order accuracy and spectral-like resolution with smaller stencils and easier application of boundary conditions, has attracted more and more interest and attention.

Two Energy-Preserving Compact Finite Difference Schemes for the Nonlinear Fourth-Order Wave Equation

In this paper,two fourth-order compact finite difference schemes are derived to solve the nonlinear fourth-order wave equation which can be viewed as a generalized model from the nonlinear beam equation.Differing from the existing compact finite difference schemes which preserve the total energy in a recursive sense,the new schemes are proved to per-fectly preserve the total energy in the discrete sense.By using the standard energy method and the cut-off function technique,the optimal error estimates of the numerical solutions are established,and the convergence rates are of O(h^(4)+τ^(2))with mesh-size h and time-step τ.In order to improve the computational efficiency,an iterative algorithm is proposed as the outer solver and the double sweep method for pentadiagonal linear algebraic equations is introduced as the inner solver to solve the nonlinear difference schemes at each time step.The convergence of the iterative algorithm is also rigorously analyzed.Several numerical results are carried out to test the error estimates and conservative properties.

A Compact Finite Volume Scheme for the Multi-Term Time Fractional Sub-Diffusion Equation

In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obtain the compact finite volume scheme have high order accuracy. We use a compact operator to deal with spatial direction;then we can get the compact finite volume scheme. It is proved that the finite volume scheme is unconditionally stable and convergent in L∞-norm. The convergence order is O(τ2-α + h4). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.

On Numerical Modelling of Industrial Powder Compaction Processes for Large Deformation of Endochronic Plasticity at Finite Strains

Compaction processes are one the most important par ts of powder forming technology. The main applications are focused on pieces for a utomotive, aeronautic, electric and electronic industries. The main goals of the compaction processes are to obtain a compact with the geometrical requirements, without cracks, and with a uniform distribution of density. Design of such proc esses consist, essentially, in determine the sequence and relative displacements of die and punches in order to achieve such goals. A.B. Khoei presented a gener al framework for the finite element simulation of powder forming processes based on the following aspects; a large displacement formulation, centred on a total and updated Lagrangian formulation; an adaptive finite element strategy based on error estimates and automatic remeshing techniques; a cap model based on a hard ening rule in modelling of the highly non-linear behaviour of material; and the use of an efficient contact algorithm in the context of an interface element fo rmulation. In these references, the non-linear behaviour of powder was adequately desc ribed by the cap plasticity model. However, it suffers from a serious deficiency when the stress-point reaches a yield surface. In the flow theory of plasticit y, the transition from an elastic state to an elasto-plastic state appears more or less abruptly. For powder material it is very difficult to define the locati on of yield surface, because there is no distinct transition from elastic to ela stic-plastic behaviour. Results of experimental test on some hard met al powder show that the plastic effects were begun immediately upon loading. In such mater ials the domain of the yield surface would collapse to a point, so making the di rection of plastic increment indeterminate, because all directions are normal to a point. Thus, the classical plasticity theory cannot deal with such materials and an advanced constitutive theory is necessary. In the present paper, the constitutive equations of powder materials will be discussed via an endochronic theory of plasticity. This theory provides a unifi ed point of view to describe the elastic-plastic behaviour of material since it places no requirement for a yield surface and a ’loading function’ to disting uish between loading an unloading. Endochronic theory of plasticity has been app lied to a number of metallic materials, concrete and sand, but to the knowledge of authors, no numerical scheme of the model has been applied to powder material . In the present paper, a new approach is developed based on an endochronic rate independent, density-dependent plasticity model for describing the isothermal deformation behavior of metal powder at low homologous temperature. Although the concept of yield surface has not been explicitly assumed in endochronic theory, it is shown that the cone-cap plasticity yield surface (Fig.1), which is the m ost commonly used plasticity models for describing the behavior of powder materi al can be easily derived as a special case of the proposed endochronic theory. Fig.1 Trace of cone-cap yield function on the meridian pl ane for different relative density As large deformation is observed in powder compaction process, a hypoelastic-pl astic formulation is developed in the context of finite deformation plasticity. Constitutive equations are stated in unrotated frame of reference that greatly s implifies endochronic constitutive relation in finite plasticity. Constitutive e quations of the endochronic theory and their numerical integration are establish ed and procedures for determining material parameters of the model are demonstra ted. Finally, the numerical schemes are examined for efficiency in the model ling of a tip shaped component, as shown in Fig.2. Fig.2 A shaped tip component. a) Geometry, boundary conditio n and finite element mesh; b) density distribution at final stage of

Finite element modeling of consolidation of composite laminates

Advanced fiber reinforced polymer composites have been increasingly applied to various structural components. One of the important processes to fabricate high performance laminated composites is an autoclave assisted prepreg lay-up. Since the quality of laminated composites is largely affected by the cure cycle, selection of an appropriate cure cycle for each application is important and must be optimized. Thus, some fundamental model of the consolidation and cure processes is necessary for selecting suitable parameters for a specific application. This article is concerned with the ＂flow-compaction＂ model during the autoclave processing of composite materials. By using a weighted residual method, two-dimensional finite element formulation for the consolidation process of thick thermosetting composites is presented and the corresponding finite element code is developed. Numerical examples, including comparison of the present numerical results with one-dimensional and twodimensional analytical solutions, are given to illustrate the accuracy and effectiveness of the proposed finite element formulation. In addition, a consolidation simulation of AS4/3501-6 graphite/epoxy laminate is carded out and compared with the experimental results available in the literature.

Explicit High-Order Method to Solve Coupled Nonlinear Schrödinger Equations

Models of the coupled nonlinear Schrödinger equations submit various critical physical phenomena with a typical equation for optical fibres with linear refraction. In this article, we will presuppose the Compact Finite Difference method with Runge-Kutta of order 4 (explicit) method, which is sixth-order and fourth-order in space and time respectively, to solve coupled nonlinear Schrödinger equations. Many methods used to solve coupled nonlinear Schrödinger equations are second order in time and need to use extra-technique to rise up to fourth-order as Richardson Extrapolation technique. The scheme obtained is immediately fourth-order in one step. This approach is a conditionally stable method. The conserved quantities and the exact single soliton solution indicate the competence and accuracy of the article’s suggestion schemes. Furthermore, the article discusses the two solitons interaction dynamics.

UNIFORM ERROR BOUNDS OF A CONSERVATIVE COMPACT FINITE DIFFERENCE METHOD FOR THE QUANTUM ZAKHAROV SYSTEM IN THE SUBSONIC LIMIT REGIME

In this paper,we consider a uniformly accurate compact finite difference method to solve the quantum Zakharov system(QZS)with a dimensionless parameter 0<ε≤1,which is inversely proportional to the acoustic speed.In the subsonic limit regime,i.e.,when 0<ε?1,the solution of QZS propagates rapidly oscillatory initial layers in time,and this brings significant difficulties in devising numerical algorithm and establishing their error estimates,especially as 0<ε?1.The solvability,the mass and energy conservation laws of the scheme are also discussed.Based on the cut-off technique and energy method,we rigorously analyze two independent error estimates for the well-prepared and ill-prepared initial data,respectively,which are uniform in both time and space forε∈(0,1]and optimal at the fourth order in space.Numerical results are reported to verify the error behavior.

Fast High Order Algorithm for Three-Dimensional Helmholtz Equation Involving Impedance Boundary Condition with Large Wave Numbers

Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundary value problems. Finite difference method is widely applied to solving these problems due to its ease of use. However, when the wave number is large, the pollution effects are still a major difficulty in obtaining accurate numerical solutions. We develop a fast algorithm for solving three-dimensional Helmholtz boundary problems with large wave numbers. The boundary of computational domain is discrete based on high-order compact difference scheme. Using the properties of the tensor product and the discrete Fourier sine transform method, the original problem is solved by splitting it into independent small tridiagonal subsystems. Numerical examples with impedance boundary conditions are used to verify the feasibility and accuracy of the proposed algorithm. Results demonstrate that the algorithm has a fourth- order convergence in and -norms, and costs less CPU calculation time and random access memory.

FOURTH-ORDER COMPACT SCHEMES FOR HELMHOLTZ EQUATIONS WITH PIECEWISE WAVE NUMBERS IN THE POLAR COORDINATES

In this paper, fourth-order compact finite difference schemes are proposed for solving Helmholtz equation with piecewise wave numbers in polar coordinates with axis-symmetric and in some cases that the solution depends both of independent variables. The idea of the immersed interface method is applied to deal with the discontinuities in the wave number and certain derivatives of the solution. Numerical experiments are included to confirm the accuracy and efficiency of the proposed method.

Aerodynamic Noise Propagation Simulation using Immersed Boundary Method and Finite Volume Optimized Prefactored Compact Scheme

Based on the immersed boundary method （IBM） and the finite volume optimized pre-factored compact （FVOPC） scheme, a numerical simulation of noise propagation inside and outside the casing of a cross flow fan is estab- lished. The unsteady linearized Euler equations are solved to directly simulate the aero-acoustic field. In order to validate the FVOPC scheme, a simulation case： one dimensional linear wave propagation problem is carried out using FVOPC scheme, DRP scheme and HOC scheme. The result of FVOPC is in good agreement with the ana- lytic solution and it is better than the results of DRP and HOC schemes, the FVOPC is less dispersion and dissi- pation than DRP and HOC schemes. Then, numerical simulation of noise propagation problems is performed. The noise field of 36 compact rotating noise sources is obtained with the rotating velocity of 1000r/min. The PML absorbing boundary condition is applied to the sound far field boundary condition for depressing the numerical reflection. Wall boundary condition is applied to the casing. The results show that there are reflections on the casing wall and sound wave interference in the field. The FVOPC with the IBM is suitable for noise propagation problems under the complex geometries for depressing the dispersion and dissipation, and also keeping the high order precision.

Numerical Simulation of Time-Harmonic Waves in Inhomogeneous Media using Compact High Order Schemes

In many problems,one wishes to solve the Helmholtz equation with variable coefficients within the Laplacian-like term and use a high order accurate method(e.g.,fourth order accurate)to alleviate the points-per-wavelength constraint by reducing the dispersion errors.The variation of coefficients in the equation may be due to an inhomogeneous medium and/or non-Cartesian coordinates.This renders existing fourth order finite difference methods inapplicable.We develop a new compact scheme that is provably fourth order accurate even for these problems.We present numerical results that corroborate the fourth order convergence rate for several model problems.

A Consistent Fourth-Order Compact Finite Difference Scheme for Solving Vorticity-Stream Function Form of Incompressible Navier-Stokes Equations

The inconsistent accuracy and truncation error in the treatment of boundary usually leads to performance defects,such as decreased accuracy and even numerical instability,of the entire computational method,especially for higher order methods.In this work,we construct a consistent fourth-order compact finite difference scheme for solving two-dimensional incompressible Navier-Stokes(N-S)equations.In the pro-posed method,the main truncation error term of the boundary scheme is kept the same as that of the interior compact finite difference scheme.With such a feature,the nu-merical stability and accuracy of the entire computation can be maintained the same as the interior compact finite difference scheme.Numerical examples show the effec-tiveness and accuracy of the present consistent compact high order scheme in L^(∞).Its application to two dimensional lid-driven cavity flow problem further exhibits that un-der the same condition,the computed solution with the present scheme is much close to the benchmark in comparison to those from the 4^(th)order explicit scheme.The compact finite difference method equipped with the present consistent boundary technique im-proves much the stability of the whole computation and shows its potential application to incompressible flow of high Reynolds number.

Optimal Error Estimates of Compact Finite Difference Discretizations for the Schrodinger-Poisson System

We study compact finite difference methods for the Schrodinger-Poisson equation in a bounded domain and establish their optimal error estimates under proper regularity assumptions on wave functionψand external potential V(x).The CrankNicolson compact finite difference method and the semi-implicit compact finite difference method are both of order O(h^(4)+τ^(2))in discrete l^(2),H^(1) and l^(∞) norms with mesh size h and time step τ.For the errors of compact finite difference approximation to the second derivative and Poisson potential are nonlocal,thus besides the standard energy method and mathematical induction method,the key technique in analysis is to estimate the nonlocal approximation errors in discrete l^(∞) and H^(1) norm by discrete maximum principle of elliptic equation and properties of some related matrix.Also some useful inequalities are established in this paper.Finally,extensive numerical results are reported to support our error estimates of the numerical methods.

A COMPACT FOURTH-ORDER FINITE DIFFERENCE SCHEME FOR THE IMPROVED BOUSSINESQ EQUATION WITH DAMPING TERMS

In this paper, a compact finite difference method is presented for solving the initial boundary value problems for the improved Boussinesq equation with damping terms. The fourth-order equation can be transformed into a first-order ordinary differential system, and then, the classical Pad4 approximation is used to discretize spatial derivative in the non- linear partial differential equations. The resulting coefficient matrix for the semi-discrete scheme is tri-diagonal and can be solved efficiently. In order to maintain the same order of convergence, the classical fourth-order Runge-Kutta method is the preferred method for explicit time integration. Soliton-type solutions are used to evaluate the accuracy of the method, and various numerical experiments are designed to test the different effects of the damping terms.

A Compact Eulerian Interface-Capturing Algorithm for Compressible Multimaterial Elastic-Plastic Flows with Mie-Gr uneisen Equation of State

This paper presents an Eulerian diffuse-interface method using a high-order compact difference scheme for simulating elastic-plastic flows with the Mie–Gruneisen(MG)equation of state(EoS).For simulations of multimaterial problems,numerical errors were generated in the material discontinuities owing to inconsistent treatment of the convective terms.Based on the normal-stress-based mechanical equilibrium assumption for elastic-plastic solids,we introduce an improved form of the consistent localized artificial diffusivity(LAD)method to ensure an oscillation-free interface for velocity and normal stress.The proposed algorithm uses a hyperelastic model.A mixture type of the model system was formed by combining the conservation equations for the basic conserved variables,an equation of a unified deviatoric tensor describing solid deformation,and an additional set of equations for solving the material quantities in the MG EoS.Several one-and two-dimensional problems with various discontinuities,including the elastic-plastic Richtmyer–Meshkov instability,were considered for testing the proposed method.

A STUDY ON NUMERICAL METHOD OF NAVIER-STOKES EQUATION AND NON-LINEAR EVOLUTION OF THE COHERENT STRUCTURES IN A LAMINAR BOUNDARY LAYER

A new method for direct numerical simulation of incompressible Navier-Stokes equations is studied in the paper. The compact finite difference and the non-linear terms upwind compact finite difference schemes on non-uniform meshes in x and y directions are developed respectively. With the Fourier spectral expansion in the spanwise direction, three-dimensional N-S equation are converted to a system of two-dimensional explicit-implicit The treatment of equations. The third-order mixed scheme is employed the three-dimensional for time integration. non-reflecting outflow boundary conditions is presented, which is important for the numerical simulations of the problem of transition in boundary layers, jets, and mixing layer. The numerical results indicate that high accuracy, stabilization and efficiency are achieved by the proposed numerical method. In addition, a theory model for the coherent structure in a laminar boundary layer is also proposed, based on which the numerical method is implemented to the non-linear evolution of coherent structure. It is found that the numerical results of the distribution of Reynolds stress, the formation of high shear layer, and the event of ejection and sweeping, match well with the observed characteristics of the coherent structures in a turbulence boundary layer.

INTERACTIVE STUDY BETWEEN IDENTICAL COHERENT STRUCTURES IN THE WALL REGION OF A TURBULENT BOUNDARY LAYER

A theoretical model for identica l coherent structures in the wall region of a turbulent boundary layer was propo sed, using the idea of general resonant triad of the hydrodynamic stability. The evolution of the structures in the wall region of a turbulent boundary layer wa s studied by combining the compact finite differences of high numerical accuracy and the Fourier spectral hybrid method for solving the three dimensional Navier -Stokes equations. In this method, the third order mixed explicit-implicit sch eme was applied for the time integration. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, an d the sixth-order center compact schemes for the derivatives in spectral space were introduced, respectively. The fourth-order compact schemes satisfied by th e velocities and pressure in spectral space was derived. As an application, the method was implemented to the wall region of a turbulent boundary to study the e volution of identical coherent structures. It is found that the numerical result s are satisfactory.

A FAST HIGH ORDER METHOD FOR ELECTROMAGNETIC SCATTERING BY LARGE OPEN CAVITIES

In this paper, the electromagnetic scattering from a rectangular large open cavity embedded in an infinite ground plane is studied. By introducing a nonlocal artificial boundary condition, the scattering problem from the open cavity is reduced to a bounded domain problem. A compact fourth order finite difference scheme is then proposed to discrete the cavity scattering model in the rectangular domain, and a special treatment is enforced to approximate the boundary condition, which makes truncation errors reach （.9（h4） in the whole computational domain. A fast algorithm, exploiting the discrete Fourier transfor- mation in the horizontal and a Gaussian elimination in the vertical direction, is employed, which reduces the discrete system to a much smaller interface system. An effective pre- conditioner is presented for the BICGstab iterative solver to solve this interface system. Numerical results demonstrate the remarkable accuracy and efficiency of the proposed method. In particular, it can be used to solve the cavity model for the large wave number up to 600π.

A HIGH-ORDER PAD SCHEME FOR KORTEWEG-DE VRIES EQUATIONS

A high-order finite difference Pade scheme also called compact scheme for solving Korteweg-de Vries （KdV） equations, which preserve energy and mass conservations, was developed in this paper. This structure-preserving algorithm has been widely applied in these years for its advantage of maintaining the inherited properties. For spatial discretization, the authors obtained an implicit compact scheme by which spatial derivative terms may be approximated through combining a few knots. By some numerical examples including propagation of single soliton and interaction of two solitons, the scheme is proved to be effective.