Modeling Convection Diffusion with Exponential Upwinding

This paper shows the usefulness of the exponential upwinding technique in convection diffusion computations. In particular, it is demonstrated that, even when convection is dominant, if exponential upwinding is employed in conjunction with either the Jacobi or the Gauss-Seidel iteration process, one can obtain computed solutions that are accurate and free of unphysical

A Fourth-Order Upwinding Embedded Boundary Method(UEBM)for Maxwell’s Equations in Media with Material Interfaces:Part I

In this paper,we present a new fourth-order upwinding embedded boundary method(UEBM)over Cartesian grids,originally proposed in the Journal of Computational Physics[190(2003),pp.159-183.]as a second-order method for treating material interfaces for Maxwell’s equations.In addition to the idea of the UEBM to evolve solutions at interfaces,we utilize the ghost fluid method to construct finite difference approximation of spatial derivatives at Cartesian grid points near the material interfaces.As a result,Runge-Kutta type time discretization can be used for the semidiscretized system to yield an overall fourth-order method,in contrast to the original second-order UEBM based on a Lax-Wendroff type difference.The final scheme allows time step sizes independent of the interface locations.Numerical examples are given to demonstrate the fourth-order accuracy as well as the stability of the method.We tested the scheme for several wave problems with various material interface locations,including electromagnetic scattering of a plane wave incident on a planar boundary and a two-dimensional electromagnetic application with an interface parallel to the y-axis.

Aerodynamic/stealth design of S-duct inlet based on discrete adjoint method

It is a major challenge for the airframe-inlet design of modern combat aircrafts,as the flow and electromagnetic wave propagation in the inlet of stealth aircraft are very complex.In this study,an aerodynamic/stealth optimization design method for an S-duct inlet is proposed.The upwind scheme is introduced to the aerodynamic adjoint equation to resolve the shock wave and flow separation.The multilevel fast multipole algorithm(MLFMA)is utilized for the stealth adjoint equation.A dorsal S-duct inlet of flying wing layout is optimized to improve the aerodynamic and stealth characteristics.Both the aerodynamic and stealth characteristics of the inlet are effectively improved.Finally,the optimization results are analyzed,and it shows that the main contradiction between aerodynamic characteristics and stealth characteristics is the centerline and crosssectional area.The S-duct is smoothed,and the cross-sectional area is increased to improve the aerodynamic characteristics,while it is completely opposite for the stealth design.The radar cross section(RCS)is reduced by phase cancelation for low frequency conditions.The method is suitable for the aerodynamic/stealth design of the aircraft airframe-inlet system.

NUMERICAL SIMULATION ON AERODYNAMIC CHARACTERISTICS OF WINDBLAST PROTECTION AT HIGH SPEED

The computational fluid dynamics （CFD） method is used to investigate the aerodynamic characteristics of the seat/occupant with windblast protection devices. The upwind Osher scheme is used for the spatial discretisation. The detached-eddy simulation （DES） based on the Spalart-Allmaras one-equation turbulence model is ap- plied to the detached viscous flow simulation behind the seat/occupant, with Mach numbers 0.6 and 1.2 at attack angles between --10 and 30°, and at two sideslip angles of 0 and 15°, respectively. The aerodynamic characteristics of seat/occupants with and without windblast protection devices are calculated in cases of the freestream Mach numbers 0. 8 and 1.6, attack angles from 5 to 30°, and three sideslip angles of 0, --20 and --50°, respectively. Results show that simulation results agree well with experimental data. And the occupant is efficiently protected by windblast protection devices.

ADAPTIVE DELAUNAY TRIANGULATION WITH MULTIDIMENSIONAL DISSIPATION SCHEME FOR HIGH-SPEED COMPRESSIBLE FLOW ANALYSIS

Adaptive Delaunay triangulation is combined with the cell-centered upwinding algorithm to analyze inviscid high-speed compressible flow problems. The multidimensional dissipation scheme was developed and included in the upwinding algorithm for unstructured triangular meshes to improve the computed shock wave resolution. The solution accuracy is further improved by coupling an error estimation procedure to a remeshing algorithm that generates small elements in regions with large change of solution gradients, and at the same time, larger elements in other regions. The proposed scheme is further extended to achieve higher-order spatial and temporal solution accuracy. Efficiency of the combined procedure is evaluated by analyzing supersonic shocks and shock propagation behaviors for both the steady and unsteady high-speed compressible flows.

APPLICATIONS OF HIGHER-ORDER UPWIND SCHEME FOR PROJECTILES IN TRANSONIC FLOW

针对弹丸跨音速粘性绕流流场进行了数值仿真研究，并与实验进行了对比。数值仿真的控制方程为雷诺平均的Navier-Stokes方程，湍流模型采用Baldwin-Lomax模型，应用了四阶的MUSCL TVD格式。结果表明，该方法可以取得高清晰度的流场结构。

PARALLELIZED UPWIND FLUX SPLITTING SCHEME FOR SUPERSONIC REACTING FLOWS ON UNSTRUCTURED HYBRID MESHES

A parallelized upwind flux splitting scheme for supersonic reacting flows on hybrid meshes is presented. The complexity of super/hyper-sonic combustion flows makes it necessary to establish solvers with higher resolution and efficiency for multi-component Euler/N-S equations. Hence, a spatial second-order van Leer type flux vector splitting scheme is established by introducing auxiliary points in interpolation, and a domain decomposition method used on unstructured hybrid meshes for obtaining high calculating efficiency. The numerical scheme with five-stage Runge-Kutta time step method is implemented to the simulation of combustion flows, including the supersonic hydrogen/air combustion and the normal injection of hydrogen into reacting flows. Satisfying results are obtained compared with limited references.

NUMERICAL METHOD OF MIXED FINITE VOLUME-MODIFIED UPWIND FRACTIONAL STEP DIFFERENCE FOR THREE-DIMENSIONAL SEMICONDUCTOR DEVICE TRANSIENT BEHAVIOR PROBLEMS

Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.

GROUP VELOCITY CONTROL SCHEME WITH LOW DISSIPATION

In order to prevent smearing the discontinuity, a modified term is added to the third order Upwind Compact Difference scheme to lower the dissipation error. Moreover, the dispersion error is controled to hold back the non physical oscillation by means of the group velocity control. The scheme is used to simulate the interactions of shock density stratified interface and the disturbed interface developing to vortex rollers. Numerical results are satisfactory.

UNSTEADY/STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.

PREDICTING THE CONSEQUENCES OF SEAWATER INTRUSION AND PROTECTION PROJECTS

The simulation of this process and the effects of protection projects lays the foundation of its effective control and defence. The mathematical model of the problem and upwind splitting alternating direction method were presented. Using this method, the numerical simulation of seawater intrusion in Laizhou Bay Area of Shandong Provivce was finished. The numerical results turned out to be identical with the real measurements, so the prediction of the consequences of protection projectects is reasonable.

Uniform Convergence Analysis of Finite Difference Scheme for Singularly Perturbed Delay Differential Equation on an Adaptively Generated Grid

Adaptive grid methods are established as valuable computational technique in approximating effectively the solutions of problems with boundary or interior layers. In this paper,we present the analysis of an upwind scheme for singularly perturbed differential-difference equation on a grid which is formed by equidistributing arc-length monitor function.It is shown that the discrete solution obtained converges uniformly with respect to the perturbation parameter.Numerical experiments illustrate in practice the result of convergence proved theoretically.

Application of Mixed Differential Quadrature Method for Solving the Coupled Two-Dimensional Incompressible Navier-Stokes Equation and Heat Equation

The traditional differential quadrature method was improved by using theupwind difference scheme for the convective terms to solve the coupled two-dimensionalincompressible Navier-stokes equations and heat equation. The new method was compared with theconventional differential quadrature method in the aspects of convergence and accuracy. The resultsshow that the new method is more accurate, and has better convergence than the conventionaldifferential quadrature method for numerically computing the steady-state solution.

Numerical analysis of gas emission rule from a goaf of tailing roadway

Targeting the problem of large amounts of gas emission from the goaf of the No.14201 working face in the Shaqu coal mine of Huajin Coking Coal Co. Ltd., we used a negative exponential function to describe the attenuation process of gas emission in goaf (the stable source) based on the principle of field flow. Equations of two-component flow (gas and air) and seep- age-diffusion in a heterogeneous goaf flow field are solved by means of numerical simulation and fluid mechanics principles of air movement and gas distribution during gas emission from goaf. The results indicate that the air diversion volume has a negative, exponential relation with the volume of gas emitted from goaf to the working face and is clearly inversely related to gas concentra- tion. We calculated the minimum amount of air diversion and distributed air volume in the tailing roadway required for safe pro- duction.

UPWIND LOCAL DIFFERENTIAL QUADRATURE METHOD FOR SOLVING COUPLED VISCOUS FLOW AND HEAT TRANSFER EQUATIONS

The differential quadrature method (DQM) has been applied successfully to solve numerically many problems in the fluid mechanics. But it is only limited to the flow problems in regular regions. At the same time, here is no upwind mechanism to deal with the convective property of the fluid flow in traditional DQ method. A local differential quadrature method owning upwind mechanism (ULDQM) was given to solve the coupled problem of incompressible viscous flow and heat transfer in an irregular region. For the problem of flow past a contraction channel whose boundary does not parallel to coordinate direction, the satisfactory numerical solutions were obtained by using ULDQM with a few grid points. The numerical results show that the ULDQM possesses advantages including well convergence, less computational workload and storage as compared with the low-order finite difference method.

NUMERICAL SCHEMES WITH HIGH ORDER OF ACCURACY FOR THE COMPUTATION OF SHOCK WAVES

High order accurate scheme is highly desirable for Slow computation with shocks. After analysis has been made for the reason of the generation of non-physical oscillations around the shock in numerical computations, a third-order, upwind biased, shock capturing scheme was proposed. Also, a new shock fitting method, called pseudo shock fitting method, was suggested, which in principle can be with any order of accuracy. Test cases for one dimensional flows show that the new method is very satisfactory.

Streamline upwind finite element method using 6-node triangular element with adaptive remeshing technique for convective-diffusion problems

A streamline upwind finite element method using 6-node triangular element is presented. The method is applied to the convection term of the governing transport equation directly along local streamlines. Several convective-diffusion examples are used to evaluate efficiency of the method. Results show that the method is monotonic and does not produce any oscillation. In addition, an adaptive meshing technique is combined with the method to further increase accuracy of the solution, and at the same time, to minimize computational time and computer memory requirement.

Streamline upwind finite element method for conjugate heat transfer problems

This paper presents a combined finite element method for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow. The streamline upwind finite element method is used for the analysis of thermal viscous flow in the fluid region, whereas the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the proposed method is to consistently couple heat transfer along the fluid-solid interface. Three test cases, i.e. conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the efficiency of the present method.

Third-Order Upwind Schemes for Convection Equations

Aim To construct a third order upwind scheme for convection equation. Methods Upwind Lagrange interpolation was used. Results and Conclusion The schemes L p stability for p∈ is proved. Numerical examples show that performance of the third order upwind scheme is better than that of most second order schemes.

SUPG finite element method based on penalty function for lid-driven cavity flow up to Re = 27500

A streamline upwind/Petrov-Galerkin （SUPG） finite element method based on a penalty function is pro- posed for steady incompressible Navier-Stokes equations. The SUPG stabilization technique is employed for the for- mulation of momentum equations. Using the penalty function method, the continuity equation is simplified and the pres- sure of the momentum equations is eliminated. The lid-driven cavity flow problem is solved using the present model. It is shown that steady flow simulations are computable up to Re = 27500, and the present results agree well with previous solutions. Tabulated results for the properties of the primary vortex are also provided for benchmarking purposes.