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

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.

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.

A Variational Multiscale Method for Particle Dispersion Modeling in the Atmosphere

A LES model is proposed to predict the dispersion of particles in the atmosphere in the context of Chemical,Biological,Radiological and Nuclear(CBRN)applications.The code relies on the Finite Element Method(FEM)for both the fluid and the dispersed solid phases.Starting from the Navier-Stokes equations and a general description of the FEM strategy,the Streamline Upwind Petrov-Galerkin(SUPG)method is formulated putting some emphasis on the related assembly matrix and stabilization coefficients.Then,the Variational Multiscale Method(VMS)is presented together with a detailed illustration of its algorithm and hierarchy of computational steps.It is demonstrated that the VMS can be considered as a more general version of the SUPG method.The final part of the work is used to assess the reliability of the implemented predictor/multicorrector solution strategy.

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格式。结果表明，该方法可以取得高清晰度的流场结构。

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.

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.

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.

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.

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.

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.

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 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.

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.

UPWIND SPLITTING SCHEME FOR CONVECTION-DIFFUSION EQUATIONS

WT5,5”BX] A new class of numerical schemes is proposed to solve convection diffusion equations by combining the upwind technique and the method of operator splitting. For every time step, the multi dimensional approximation is performed in several independent directions alternatively, while the upwind technique is applied to treat the convection term in every individual direction. This scheme possesses maximum principle. Stability and convergence are analysed by energy method.[WT5,5”HZ]

LONG-TIME CONVERGENCE OF GENERALIZED DIFFERENCE METHOD FOR NAVIER-STOKES EQUATIONS

In this paper, we first provide a generalized difference method for the 2-dimensional Navier-Stokes equations by combing the ideas of staggered scheme m and generalized upwind scheme in space, and by backward Euler time-stepping. Then we apply the abstract framework of to prove its long-time convergence. Finally, a numerical example for solving driven cavity flows is given.

Developing shock-capturing difference methods

A new shock-capturing method is proposed which is based on upwind schemes and flux-vector splittings. Firstly, original upwind schemes are projected along characteristic directions. Secondly, the amplitudes of the characteristic decompositions are carefully controlled by limiters to prevent non-physical oscillations. Lastly, the schemes are converted into conservative forms, and the oscillation-free shock-capturing schemes are acquired. Two explicit upwind schemes （2nd-order and 3rd-order） and three compact upwind schemes （3rd-order, 5th-order and 7th-order） are modified by the method for hyperbolic systems and the modified schemes are checked on several one-dimensional and two-dimensional test cases. Some numerical solutions of the schemes are compared with those of a WENO scheme and a MP scheme as well as a compact-WENO scheme. The results show that the method with high order accuracy and high resolutions can capture shock waves smoothly.

Inlet flow disturbance effects on direct numerical simulation of incompressible round jet at Reynolds number 2500

This letter reports inlet flow disturbance effects on direct numerical simulation of incompressible round jet at Reynolds number 2500.The simulation employs an accurate projection method in which a sixth order biased upwind difference scheme is used for spatial discretization of nonlinear convective terms,with a fourth order central difference scheme used in the discretization of the divergence of intermediate velocity.Carefully identifying reveals that the inlet flow disturbance has some influences on the distribution pattern of mean factor of swirling strength intermittency.With the increase of inlet disturbance magnitude jet core cone slightly shortens,observable differences occur in the centerline velocity and its fluctuations,despite the negligible impacts on the least square fitted centerline velocity decay constant(B_u)and distribution parameter(K_u)for velocity profile in self-similar region.