Modified characteristic finite difference fractional step method for moving boundary value problem of percolation coupled system

For the coupled system with moving boundary values of multilayer dynamicsof fluids in porous media,a characteristic finite difference fractional step scheme appli-cable to the parallel arithmetic is put forward.Some techniques,such as the change ofregions,the calculus of variations,the piecewise threefold quadratic interpolation,themultiplicative commutation rule of difference operators,the decomposition of high orderdifference operators,and the prior estimates,are adopted.The optimal order estimatesin the l2norm are derived to determine the error in the approximate solution.This nu-merical method has been successfully used to simulate the flow of migration-accumulationof the multilayer percolation coupled system.Some numerical results are well illustratedin this paper.     

NUMERICAL METHOD AND APPLICATION FOR THREEDIMENSIONAL NONLINEAR SYSTEM OF DYNAMICS OF FLUIDS IN POROUS MEDIA

For the system of multilayer dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources.     

NUMERICAL METHOD FOR THREE-DIMENSIONAL NONLINEAR CONVECTION-DOMINATED PROBLEM OF DYNAMICS OF FLUIDS IN POROUS MEDIA

For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fractional steps techniques are needed to convert a multi-dimensional problem into a series of successive one-dimensional problems. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators, and the theory of prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources and predicting the consequences of seawater intrusion and protection projects.     

Characteristic fractional step finite difference method for nonlinear section coupled system

For the section coupled system of multilayer dynamics of fluids in porous media, a parallel scheme modified by the characteristic finite difference fractional steps is proposed for a complete point set consisting of coarse and fine partitions. Some tech- niques, such as calculus of variations, energy method, twofold-quadratic interpolation of product type, multiplicative commutation law of difference operators, decomposition of high order difference operators, and prior estimates, are used in theoretical analysis. Optimal order estimates in 12 norm are derived to show accuracy of the second order approximation solutions. These methods have been used to simulate the problems of migration-accumulation of oil resources.     

Discrete and non-local elastica

In this paper, the buckling and post-buckling behavior of an elastic lattice system referred to as the discrete elastica problem is investigated using an equivalent non-local continuum approach. The geometrically exact post-buckling analysis of the elastic chain, also called Hencky system, is first numerically solved using the shooting method. This discrete physical model is also mathematically equivalent to a finite difference formulation of the continuum elastica. Starting from the exact difference equations of the discrete problem, a continualization method is applied for approximating the difference operators by differential ones, in order to better characterize the discrete system by an enriched continuous one. It is shown that the new continuum associated with the discrete system exactly fits the discrete elastica post-buckling problem, where the non-locality is of Eringen׳s type (also called stress gradient non-local model). An asymptotic expansion is performed for both the discrete and the non-local continuum models, in order to approximate the post-buckling branches of the discrete system. Some numerical investigations show the efficiency of the non-local approach, especially for capturing the scale effects inherent to the cell size of the lattice model.     

The Time-Scaled Trapezoidal Integration Rule for Discrete Fractional Order Controllers

In this paper, the time-scaled trapezoidal integration rule for discretizing fractional order controllers is discussed. This interesting proposal is used to interpret discrete fractional order control (FOC) systems as control with scaled sampling time. Based on this time-scaled version of trapezoidal integration rule, discrete FOC can be regarded as some kind of control strategy, in which strong control action is applied to the latest sampled inputs by using scaled sampling time. Namely, there are two time scalers for FOC systems: a normal time scale for ordinary feedback and a scaled one for fractional order controllers. A new realization method is also proposed for discrete fractional order controllers, which is based on the time-scaled trapezoidal integration rule. Finally, a one mass position 1/s ^k control system, realized by the proposed method, is introduced to verify discrete FOC systems and their robustness against saturation non-linearity.     

Testing the accuracy of the overlap criterion

Here we investigate the accuracy of the overlap criterion when applied to a simple near-integrable model in both its 2D and 3D versions. To this end, we consider, respectively, two and three quartic oscillators as the unperturbed system, and couple the degrees of freedom by a cubic, non-integrable perturbation. For both systems we compute the unperturbed resonances up to order O(ε²), and model each resonance by means of the pendulum approximation in order to estimate the theoretical critical value of the perturbation parameter for a global transition to chaos. We perform several surface of sections for the bi-dimensional case to derive an empirical value to be compared to our theoretical estimation. Although both values are of the same order of magnitude, there is a significant difference between them. For the 3D case a numerical estimate is attained that we observe matches quite well the critical value resulting from theoretical means. This confirms once again that calculating resonances up to O(ε²) suffices in order the overlap criterion to work out.     

An efficient and accurate fully discrete finite element method for unsteady incompressible Oldroyd fluids with large time step

This paper proposes a second‐order accuracy in time fully discrete finite element method for the Oldroyd fluids of order one. This new approach is based on a finite element approximation for the space discretization, the Crank–Nicolson/Adams–Bashforth scheme for the time discretization and the trapezoid rule for the integral term discretization. It reduces the nonlinear equations to almost unconditionally stable and convergent systems of linear equations that can be solved efficiently and accurately. Here, the numerical simulations for L², H¹ error estimates of the velocity and L² error estimates of the pressure at different values of viscoelastic viscosities α, different values of relaxation time λ₁, different values of null viscosity coefficient μ₀ are shown. In addition, two benchmark problems of Oldroyd fluids with different solvent viscosity μ and different relaxation time λ₁ are simulated. All numerical results perfectly match with the theoretical analysis and show that the developed approach gives a high accuracy to simulate the Oldroyd fluids under a large time step. Furthermore, the difference and the connection between the Newton fluids and the viscoelastic Oldroyd fluids are displayed. Copyright © 2015 John Wiley & Sons, Ltd.     

Anticlastic behavior of flat plates

The purpose of this paper is to show that the readings from strain gages can be used effectively to compute small transverse deflections in a rectangular plate and, further, show that the theory developed by Lamb for the rectangular-plate problem agrees with experiment. A numerical procedure is developed, based on the trapezoidal rule, which determines the transverse deflections from the readings of strain gages mounted to the top and bottom surface of a rectangular plate subjected to large longitudinal curvatures. It is shown using the strain-gage technique that experiment agrees with Lamb's theory forb ² /Rt ratios up to 50.     

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.     

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.     

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.     

On the symmetry of deformable crystals

We prove existence, uniqueness, and regularity properties for a solution u of the Bellman-Dirichlet equation of dynamic programming: (1) $$\left\{ \begin{gathered} \max {\text{ }}\{ L^i u + f^i = 0{\text{ in }}\Omega \hfill \\ i{\text{ = 1,2 }} \hfill \\ u{\text{ = 0 on }}\partial \Omega , \hfill \\ \end{gathered} \right.$$ where L ¹ and L ² are two second order, uniformly elliptic operators. The method of proof is to rephrase (1) as a variational inequality for the operator K=L ²(L ¹)^?1 in L ²(Ω) and to invoke known existence theorems. For sufficiently nice f ¹ and f ² we prove in addition that u is in H ³(Ω)?C ^2,α(Ω) (for some 0<α<1) and hence is a classical solution of (1).     

On a formulation for anisotropic elastoplasticity at finite strains invariant with respect to the intermediate configuration

The paper is concerned with a formulation of anisotropic finite strain inelasticity based on the multiplicative decomposition of the deformation gradient F=F_eF_p. A major feature of the theory is its invariance with respect to rotations superimposed on the inelastic part of the deformation gradient. The paper motivates and shows how such an invariance can be achieved. At the heart of the formulation is the mixed-variant transformation of the structural tensor, defined as the tensor product of the privileged directions of the material as given in a reference configuration, under the action of F_p. Issues related to the plastic material spin are discussed in detail. It is shown that, in contrast to the isotropic case, any flow function formulated purely in terms of stress quantities, necessarily exhibits a non-vanishing plastic material spin. The possible construction of spin-free rates is discussed as well, where it is shown that the flow rule must then depend not only on the stress but on the strain as well.     

A compact finite difference method on staggered grid for Navier–Stokes flows

Compact finite difference methods feature high‐order accuracy with smaller stencils and easier application of boundary conditions, and have been employed as an alternative to spectral methods in direct numerical simulation and large eddy simulation of turbulence. The underpinning idea of the method is to cancel lower‐order errors by treating spatial Taylor expansions implicitly. Recently, some attention has been paid to conservative compact finite volume methods on staggered grid, but there is a concern about the order of accuracy after replacing cell surface integrals by average values calculated at centres of cell surfaces. Here we introduce a high‐order compact finite difference method on staggered grid, without taking integration by parts. The method is implemented and assessed for an incompressible shear‐driven cavity flow at Re = 10³, a temporally periodic flow at Re = 10⁴, and a spatially periodic flow at Re = 10⁴. The results demonstrate the success of the method. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermodynamics applied to gradient theories involving the accumulated plastic strain: The theories of Aifantis and Fleck and Hutchinson and their generalization

We discuss the physical nature of flow rules for rate-independent (gradient) plasticity laid down by Aifantis and by Fleck and Hutchinson. As a central result we show that:

•

the flow rule of Fleck and Hutchinson is incompatible with thermodynamics unless its nonlocal term is dropped. If the underlying theory is augmented by a general defect energy dependent on γ^p and ∇γ^p, then compatibility with thermodynamics requires that its flow rule reduce to that of Aifantis.

We establish this result (and others) within a general framework obtained by combining a virtual-power principle of Fleck and Hutchinson with the first two laws of thermodynamics—balance of energy and the Clausius-Duhem inequality—under isothermal conditions.     

Combined convection and wall conduction effects in laminar pipe flow: numerical predictions and experimental validation under uniform wall heating

Laminar combined convection in horizontal circular ducts is investigated both numerically and experimentally, under uniform wall heating. A series of experiments for the heating of water in a long horizontal copper tube are simulated numerically in order to assess the reliability of the theoretical results. Peripheral and axial wall conduction effects, inherently present in the experiments, are accounted for in the numerical model. The cross validation of experimental and numerical data allows significant conclusions to be reached on conjugate conduction and convection with buoyancy effects in horizontal duct flows. Buoyancy is considered for values of the modified Rayleigh number,Ra _qo , up to 5·10⁶; the forced convection contribution is considered for two values of the entry Reynolds number,Re _o=500 and 1000.     

Poisson–Nernst–Planck Systems for Ion Flow with Density Functional Theory for Hard-Sphere Potential: I–V Relations and Critical Potentials. Part II: Numerics

We consider a one-dimensional steady-state Poisson–Nernst–Planck type model for ionic flow through membrane channels. Improving the classical Poisson–Nernst–Planck models where ion species are treated as point charges, this model includes ionic interaction due to finite sizes of ion species modeled by hard sphere potential from the Density Functional Theory. The resulting problem is a singularly perturbed boundary value problem of an integro-differential system. We examine the problem and investigate the ion size effect on the current–voltage (I–V) relations numerically, focusing on the case where two oppositely charged ion species are involved and only the hard sphere components of the excess chemical potentials are included. Two numerical tasks are conducted. The first one is a numerical approach of solving the boundary value problem and obtaining I–V curves. This is accomplished through a numerical implementation of the analytical strategy introduced by Ji and Liu in [Poisson–Nernst–Planck systems for ion flow with density functional theory for hard-sphere potential: I–V relations and critical potentials. Part I: Analysis, J. Dyn. Differ. Equ. (to appear)]. The second task is to numerically detect two critical potential values V _c and V ^c .The existence of these two critical values is first realized for a relatively simple setting and analytical approximations of V _c and V ^c are obtained in the above mentioned reference. We propose an algorithm for numerical detection of V _c and V ^c without using any analytical formulas but based on the defining properties and numerical I–V curves directly. For the setting in the above mentioned reference, our numerical values for V _c and V ^c agree well with the analytical predictions. For a setting including a nonzero permanent charge in which case no analytic formula for the I–V relation is available now, our algorithms can still be applied to find V _c and V ^c numerically.     

Strictly stable high order difference approximations for computational aeroacoustics

High order finite difference approximations with improved accuracy and stability properties have been developed for computational aeroacoustics (CAA). One of our new difference operators corresponds to Tam and Webb's DRP scheme in the interior, but is modified near the boundaries to be strictly stable. A unified formulation of the nonlinear and linearized Euler equations is used, which can be extended to the Navier–Stokes equations. The approach has been verified for 1D, 2D and axisymmetric test problems. We have simulated the sound propagation from a rocket launch before lift-off. To cite this article: B. Müller, S. Johansson, C. R. Mecanique 333 (2005).     

Conditional Relative Acceleration Statistics and Relative Dispersion Modelling

We have used direct numerical simulation results for the Eulerian velocity difference probability density function and the mean acceleration difference conditioned on the velocity difference, to explore some of the assumptions underlying the formulation of Lagrangian stochastic models for relative dispersion. We focussed on the ability of the models to quantitatively represent Richardson’s t ³-law and in particular the value of Richardson’s constant. As a result of intermittency, with decreasing spatial separation and with increasing Reynolds number these Eulerian quantities become more extreme and the model predictions for Richardson’s constant also become more extreme (larger). This is in contrast with recent numerical simulations showing that Richardson’s constant depends only weakly on Reynolds number. We conclude that, at least in the present Lagrangian stochastic modelling framework, in two-particle models (and presumably in multi-particle models) intermittency must be included explicitly in the dissipation rate as well as in the relative velocity statistics.