H 1 space-time discontinuous finite element method for convection-diffusion equations

An H~1 space-time discontinuous Galerkin (STDG) scheme for convectiondiffusion equations in one spatial dimension is constructed and analyzed. This method is formulated by combining the H~1 Galerkin method and the space-time discontinuous finite element method that is discontinuous in time and continuous in space. The existence and the uniqueness of the approximate solution are proved. The convergence of the scheme is analyzed by using the techniques in the finite difference and finite element methods. An optimal a-priori error estimate in the L~∞ (H~1 ) norm is derived. The numerical exper- iments are presented to verify the theoretical results.     

An extension of the SIMPLE based discontinuous Galerkin solver to unsteady incompressible flows

In this paper, we present a SIMPLE based algorithm in the context of the discontinuous Galerkin method for unsteady incompressible flows. Time discretization is done fully implicit using backward differentiation formulae (BDF) of varying order from 1 to 4. We show that the original equation for the pressure correction can be modified by using an equivalent operator stemming from the symmetric interior penalty (SIP) method leading to a reduced stencil size. To assess the accuracy as well as the stability and the performance of the scheme, three different test cases are carried out: the Taylor vortex flow, the Orr‐Sommerfeld stability problem for plane Poiseuille flow and the flow past a square cylinder. (1) Simulating the Taylor vortex flow, we verify the temporal accuracy for the different BDF schemes. Using the mixed‐order formulation, a spatial convergence study yields convergence rates of k + 1 and k in the L²‐norm for velocity and pressure, respectively. For the equal‐order formulation, we obtain approximately the same convergence rates, while the absolute error is smaller. (2) The stability of our method is examined by simulating the Orr–Sommerfeld stability problem. Using the mixed‐order formulation and adjusting the penalty parameter of the symmetric interior penalty method for the discretization of the viscous part, we can demonstrate the long‐term stability of the algorithm. Using pressure stabilization the equal‐order formulation is stable without changing the penalty parameter. (3) Finally, the results for the flow past a square cylinder show excellent agreement with numerical reference solutions as well as experiments. Copyright © 2015 John Wiley & Sons, Ltd.     

Analytical and numerical solutions of a coupled non-linear system arising in a three-dimensional rotating flow

Hydromagnetic flow between two horizontal plates in a rotating system, where the lower is a stretching sheet and the upper is a porous solid plate (in the presence of a magnetic field), is analyzed. The equations of conservation of mass and momentum are reduced to a system of coupled non-linear ordinary differential equations. These basic non-linear differential equations, for the velocity field (f′,f,g), are solved numerically by using a fourth-order Runge-Kutta integration scheme. The numerical results thus obtained are validated by the analytical results (for small R) obtained by the perturbation technique and presented through graphs. Also, the effects of the non-dimensional parameters R, λ, M² and K² on the velocity field are discussed, and it is shown that for large K², the coriolis force and the magnetic field that act against the pressure gradient cause reverse flow.     

Superconvergence analysis of bi-<Emphasis Type="Italic">k</Emphasis>-degree rectangular elements for two-dimensional time-dependent Schrödinger equation

Superconvergence has been studied for long, and many different numerical methods have been analyzed. This paper is concerned with the problem of superconvergence for a two-dimensional time-dependent linear Schrödinger equation with the finite element method. The error estimate and superconvergence property with order O(h^k+1) in the H¹ norm are given by using the elliptic projection operator in the semi-discrete scheme. The global superconvergence is derived by the interpolation post-processing technique. The superconvergence result with order O(h^k+1 + τ²) in the H¹ norm can be obtained in the Crank-Nicolson fully discrete scheme.     

Variational approach to non-linear boundary value problems for elasto-plastic incompressible bending plate

Non-linear boundary value problems for inelastic isotropic homogeneous incompressible bending plate, within the range of J₂-deformation theory, are considered. An existence of the weak solution of the non-linear problem with clamped boundary condition is obtained in H²(Ω) by using monotone operator theory and Browder-Minty theorem. For linearization of the non-linear problem a monotone iteration scheme is constructed. It is shown that the sequence of potentials obtained from the sequence of approximate solutions (i.e. iterations), is a monotone decreasing one. Convergence of the iteration process in H²-norm is proved by using the convexity argument. Numerical solutions, based on finite-difference scheme, are given for linear bending problems with rigid clamped as well as simply supported boundary conditions. Further numerical examples are presented to illustrate the convergence of approximate solutions and monotonicity of the potentials as applied to the non-linear problems.     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

Effective thermal conductivity of foam concretes: Homogenization schemes vs experimental data and FEM simulations

This paper focuses on the prediction of thermal conductivity of foam concretes. Various analytical forms for their effective thermal conductivity according to their porosity p (air voids volume) have been first derived based on five well-known Mean-Field Homogenization (MFH) schemes. These predictions were found to be very close for low porosities but move away gradually with increasing porosity. Thus, in order to determine the best homogenization scheme predicting the effective thermal conductivity of foam concretes, MFH predictions were confronted with experimental data obtained on nineteen foam concretes and also with numerical results obtained from 3d Finite Element Method (FEM) simulations conducted on an idealized foam concrete. These comparisons have shown that the normalized effective thermal conductivity of foam concrete is closely framed by the power law (1- p)^3/2 given by the Differential scheme and by the hyperbolic law (1- p)/(1+ p/2) given by the Mori-Tanaka scheme.     

Measurements of non-equilibrium and equilibrium temperature behind a strong shock wave in simulated martian atmosphere

Non-equilibrium radiation measurements behind strong shock wave for simulated Martian atmosphere are presented in this paper. The shock wave is established in a hydrogen oxygen combustion driven shock tube. Timeresolved spectra of the Δv = 0 sequence of the B ²Σ⁺ → X ²Σ⁺ electronic transition of CN have been observed through optical emission spectroscopy (OES). A new method, which is based on fitting high resolution spectrum for rotational and vibrational temperatures measurement, is proposed to diagnose temperature distribution behind the shock wave. It is estimated that the current scheme has the maximum deviation less than 8% (1σ) for vibrational temperature measurement through detailed analysis of the influence of the uncertainties of spectroscopic constants and spectral resolution. Radiation structure of the shock layer, including induction, relaxation and equilibrium process, and corresponding rotational and vibrational temperatures are obtained through time gating OES diagnostics with sub-microsecond temporal resolution. The present extensive results will strongly benefit the reaction rate estimation and computational fluid dynamics (CFD) code validation in high enthalpy Mars reentry chemistry.     

New nonconforming finite element method for solving transient Naiver-Stokes equations

For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral finite element spaces. The semi- and full-discrete schemes of the stabilized method are studied based on the pressure projection and a variational multi-scale method. It has some attractive features： avoiding higher-order derivatives and edge-based data structures, adding a discrete velocity term only on the fine scale, being effective for high Reynolds number fluid flows, and avoiding increased computation cost. For the full-discrete scheme, it has second-order estimations of time and is unconditionally stable. The presented numerical results agree well with the theoretical results.     

An adaptive viscoelastic stress splitting scheme and its applications: AVSS/SI and AVSS/SUPG

We report an adaptive viscoelastic stress splitting (AVSS) scheme, which was found to be extremely robust in the numerical simulation of viscoelastic flow involving steep stress boundary layers. The scheme is different from the elastic viscous split stress (EVSS) in that the local Newtonian component is allowed to depend adaptively on the magnitude of the local elastic stress. Two decoupled versions of the scheme were implemented for the Upper Convected Maxwell (UCM) model: the streamline integration (AVSS/SI), and the streamline upwind Petrov-Galerkin (AVSS/SUPG) methods of integrating the stress. The implementations were benchmarked against the known analytic Poiseuille solution, and no upper limiting Weissenberg number was found (the computation was stopped at Weissenberg number of O(10⁴)). The flow past a sphere in a tube was solved next, giving convergent results up to a Weissenberg number of 3.2 with the AVSS/SI and 1.55 with the AVSS/SUPG (both were decoupled schemes; without the adaptive scheme, the limiting Weissenberg number for the decoupled streamline integration method was about 0.3). These results show that the decoupled AVSS is more stable than the corresponding EVSS, and the SI is more robust than SUPG in solving the constitutive equation of hyperbolic type. In addition, we found a very long stress wake behind the sphere, and a weak vortex in the rear stagnation region at a Weissenberg number above W_i of about 1.6, which does not seem to increase in size or strength with increasing W_i.     

Multiple spatial and temporal scales method for numerical simulation of non-classical heat conduction problems: one dimensional case

A multiple spatial and temporal scales method is studied to simulate numerically the phenomenon of non-Fourier heat conduction in periodic heterogeneous materials. The model developed is based on the higher-order homogenization theory with multiple spatial and temporal scales in one dimensional case. The amplified spatial scale and the reduced temporal scale are introduced respectively to account for the fluctuations of non-Fourier heat conduction due to material heterogeneity and non-local effect of the homogenized solution. By separating the governing equations into various scales, the different orders of homogenized non-Fourier heat conduction equations are obtained. The reduced time dependence is thus eliminated and the fourth-order governing differential equations are derived. To avoid the necessity of C¹ continuous finite element implementation, a C⁰ continuous mixed finite element approximation scheme is put forward. Numerical results are shown to demonstrate the efficiency and validity of the proposed method.     

Adaptive fuzzy H ∞ tracking design of SISO uncertain nonlinear fractional order time-delay systems

In this paper, a fuzzy logic controller equipped with training algorithms is developed such that the H ^?? tracking performance should be satisfied for a model-free nonlinear fractional order time delay system which is infinite dimensional in nature and time delay is a source of instability. In order to deal with the linguistic uncertainties caused from delay terms, the adaptive time delay fuzzy logic system is constructed to approximate the unknown time delay system functions. By incorporating Lyapunov stability criterion with H ^?? tracking design technique, the free parameters of the adaptive fuzzy controller can be tuned on line by output feedback control law and adaptive law. Moreover, the tracking error and external disturbance can be attenuated to arbitrary desired level. The numerical results show the effectiveness of the proposed adaptive H ^?? tracking scheme.     

Effective elastic properties of porous materials: Homogenization schemes vs experimental data

In this paper, we focus on the prediction of elastic moduli of isotropic porous materials made of a solid matrix having a Poisson's ratio vm of 0.2. We derive simple analytical formulae for these effective moduli based on well-known Mean-Field Eshelby-based Homogenization schemes. For each scheme, we find that the normalized bulk, shear and Young's moduli are given by the same form depending only on the porosity p. The various predictions are then confronted with experimental results for the Young's modulus of expanded polystyrene (EPS) concrete. The latter can be seen as an idealized porous material since it is made of a bulk cement matrix, with Poisson's ratio 0.2, containing spherical mono dispersed EPS beads. The Differential method predictions are found to give a very good agreement with experimental results. Thus, we conclude that when vm=0.2, the normalized effective bulk, shear and Young's modulus of isotropic porous materials can be well predicted by the simple form (1 − p)² for a large range of porosity p ranging between 0 and 0.56.     

Similarity solution of mixed convection over a horizontal moving plate

Investigation to the mixed convective heat and mass transfer over a horizontal plate has been carried out. By applying transformation group theory to analysis of the governing equations of continuity, momentum, energy and diffusion, we show the existence of similarity solution for the problem provided that the temperature and concentration at the wall are proportional to x ^4/(7-5n) and that the moving speed of the plate is proportional to x ^(3-n)/(7-5n), and further obtain a similarity representation of the problem. The similarity equations have been solved numerically by a fourth-order Runge–Kutta scheme. The numerical results obtained for Pr=0.72 and various values of the parameters Sc, K ₁, K ₂ and K ₃ reveals the influence of the parameters on the flow, heat and mass transfer behavior.     

Stream function finite element solution of Stokes flow with penalties

A finite element method for solution of the stream function formulation of Stokes flow is developed. The method involves complete cubic non-conforming (C⁰) triangular Hermite elements. This element fails the patch test. To correct the element and produce a convergent method we employ a penalty method to weakly enforce the desired continuity constraint on the normal derivative across the inter-element boundaries. Successful use of the method is demonstrated to require reduced integration of the inter-element penalty with a 1-point Gauss rule. Error estimates relate the optimal choice of penalty parameter to mesh size and are corroborated by numerical convergence studies. The need for reduced integration is interpreted using rank relations for an associated hybrid method.     

Multivalued penalty method for evolution variational inequalities with <Emphasis Type="Italic">λ</Emphasis><Subscript>0</Subscript>-pseudomonotone multivalued maps

We consider evolution variational inequalities with λ ₀-pseudomonotone maps. The main properties of these maps are investigated. By using the finite-difference method, we prove the property of strong solvability for the class of evolution variational inequalities with λ ₀-pseudomonotone maps. Using the penalty method for multivalued maps, we show the existence of weak solutions of evolution variational inequalities on closed convex sets. The class of multivalued penalty operators is constructed. We also consider a model example to illustrate this theory.     

Synchronization of PWL function-based 2D and 3D multi-scroll chaotic systems

We study the synchronization of a piecewise linear function-based chaotic system. That system generates multiple scrolls in multiple directions (two- and three-directions) on phase space. In this scenario, the design of a controller based on Generalized Hamiltonian forms is possible. As function of control signals, we propose a master?Cslave synchronization scheme using 2 ⁿ ?1 combinations to drive a nonlinear state observer. Associated with this, the piecewise linear functions of the slave are directly controlled by the state-variables of the master system. We computed the synchronization error for each combinations. Besides, the circuit synthesis based on operational amplifiers validates our synchronization scheme by means of SPICE simulations. We observed that the synchronization error at circuit level depends on the number of the control signals used. Our numerical and SPICE simulation results are in agreement showing the usefulness of the proposed approach.     

Mixed convection boundary layer flow over a horizontal circular cylinder with Newtonian heating

The steady mixed convection boundary layer flow over a horizontal circular cylinder, generated by Newtonian heating in which the heat transfer from the surface is proportional to the local surface temperature, is considered in this study. The governing boundary layer equations are first transformed into a system of non-dimensional equations via the non-dimensional variables, and then into non-similar equations before they are solved numerically using a numerical scheme known as the Keller-box method. Numerical solutions are obtained for the skin friction coefficient Re ^1/2 C _f and the local wall temperature θ _w (x) as well as the velocity and temperature profiles with two parameters, namely the mixed convection parameter λ and the Prandtl number Pr.     

Finite difference analysis of unsteady natural convection flow along an inclined plate with variable surface temperature and mass diffusion

An analysis is performed to study the transient laminar natural convection flows along an inclined semi-infinite flat plate in which the wall temperatureT _w ^′ and species concentration on the wallC _w ^′ vary as the power of the axial co-ordinate in the formT _w ^′ (x)=T _∞ ^′ +ax ⁿ andC _w ^′ =C _∞ ^′ +bx ^m respectively. The dimensionless governing equations considered here are unsteady, two-dimensional, coupled and non-linear integro-differential equations. A finite difference scheme of Crank-Nicolson type is employed to solve the problem. The velocity, temperature, concentration, skin friction, Nusselt number and Sherwood number are studied in detail for various sets of values of parameters. Correlation equations are also established for Nusselt number and Sherwood number in terms of parameters.     

An automated vortex detection scheme using the wavelet transform of the<Emphasis Type="Italic"> d</Emphasis><Subscript>2</Subscript> field

A vortex detection scheme using the wavelet transform of the discriminant of the eigenvalues of the velocity gradient matrix (d ₂) is presented in this paper. The use of d ₂ field results in better eduction capability over the previously used test fields such as enstrophy since it automatically distinguishes shear layers from vortices. Level sets are used to refine the shape of the vortex without causing a huge computational penalty. Further, the algorithm is easily automated to aid batch processing. The detection scheme was applied to swirl flow fields, successfully estimating the vortex location, the core radius and the vortex shape.