New Results on Global Rank Axioms of Poset Matroids

An excellent introduction to the topic of poset matroids is due to M. Barnabei, G. Nicoletti and L. Pezzoli. On the basis of their work, we have obtained the global rank axioms for poset matroids.In this paper, we study the special integral function f and obtain a new class of poset matroids from the old ones, and then we generalize this result according to the properties of f. Almost all of these results can be regarded as the application of global rank axioms for poset matroids. The main results in our paper have, indeed, investigated the restriction of the basis of the poset matroid, and we give them the corresponding geometric interpretation.     

Adaptive mesh refinement for high‐resolution finite element schemes

New a posteriori error indicators based on edgewise slope‐limiting are presented. The L₂‐norm is employed to measure the error of the solution gradient in both global and element sense. A second‐order Newton–Cotes formula is utilized in order to decompose the local gradient error from a ??₁ finite element solution into a sum of edge contributions. The slope values at edge midpoints are interpolated from the two adjacent vertices. Traditional techniques to recover (superconvergent) nodal gradient values from consistent finite element slopes are reviewed. The deficiencies of standard smoothing procedures—L₂‐projection and the Zienkiewicz–Zhu patch recovery—as applied to nonsmooth solutions are illustrated for simple academic configurations. The recovered gradient values are corrected by applying a slope limiter edge‐by‐edge so as to satisfy geometric constraints. The direct computation of slopes at edge midpoints by means of limited averaging of adjacent gradient values is proposed as an inexpensive alternative. Numerical tests for various solution profiles in one and two space dimensions are presented to demonstrate the potential of this postprocessing procedure as an error indicator. Finally, it is used to perform adaptive mesh refinement for compressible inviscid flow simulations. Copyright © 2006 John Wiley & Sons, Ltd.     

On monotonicity and boundedness properties of linear multistep methods

In this paper an analysis is provided of nonlinear monotonicity and boundedness properties for linear multistep methods. Instead of strict monotonicity for arbitrary starting values we shall focus on generalized monotonicity or boundedness with Runge-Kutta starting procedures. This allows many multistep methods of practical interest to be included in the theory. In a related manner, we also consider contractivity and stability in arbitrary norms.

     

A comparative study of the performance of high resolution advection schemes in the context of data assimilation

High resolution advection schemes have been developed and studied to model propagation of flows involving sharp fronts and shocks. So far the impact of these schemes in the framework of inverse problem solution has been studied only in the context of linear models. A detailed study of the impact of various slope limiters and the piecewise parabolic method (PPM) on data assimilation is the subject of this work, using the nonlinear viscous Burgers equation in 1?D. Also provided are results obtained in 2?D using a global shallow water equations model. The results obtained in this work may point out to suitability of these advection schemes for data assimilation in more complex higher dimensional models. Copyright © 2005 John Wiley & Sons, Ltd.     

Isomerism in Ag and non-yrast levels in Pd and Rh, studied in β decay

The β decay of ⁹⁶Ag (Z=47,N=49) was investigated by measuring positrons, X rays as well as β-delayed protons and γ rays. The γ radiation was studied by means of germanium detectors and a NaI total-absorption spectrometer. Two β-decaying isomers in ⁹⁶Ag were established with half-lives of 4.40(6) and 6.9(6) s and tentative spin–parity assignments of (8⁺) and (2⁺), respectively. For both isomers, the intensities of β transitions to low-lying levels of ⁹⁶Pd (Z=46,N=50) and β-delayed proton decays to levels in ⁹⁵Rh (Z=45,N=50) were measured. Several new ⁹⁶Pd levels were firmly established. The level energies, their γ decays and the Gamow–Teller decay of ⁹⁶Ag are compared to shell-model predictions. A new low-lying level in ⁹⁵Rh was found at 680 keV excitation energy. Through a comparison with low-lying states of N=50 isotones, this level is interpreted as the first excited 7/2⁺ state built on the proton 9/2⁺ ground state. The assignments of further excited states in ⁹⁵Rh are discussed.     

一类三维时空二阶精度高分辨率MmB差分格式

直接把Nessyahu和Tadmor^[1,2]的思想推广到三维非线性双曲型守恒律情形，以交错形式Lax—Friedrichs格式为基本模块，使用二阶分片线性逼近代替一阶分片常数逼近，减少了Lax—Friedrichs格式的过多数值粘性，通过对混合导数离散形式的适当处理，构造了一类不须解Riemann问题、具有时空二阶精度高分辨率的MmB差分格式。这些差分格式很容易推广到向量系统中去。最后，一些数值模拟计算结果也证明了这些差分格式的有效性。     

Nonlinear Dynamics of Shells: Theory,Finite Element Formulation,and Integration Schemes

The paper is concerned with a dynamical formulation of a recently established shell theory capable to catch finite deformations and falls within the class of geometrically exact shell theories. A basic aspect is the design of time integration schemes which preserve specific features of the continuous system such as conservation of momentum, angular momentum, and energy when the applied forces allow to. The integration method differs from the one recently proposed by Simo and Tarnow in being applicable without modifications to shell formulations with linear as well as nonlinear configuration spaces and in being independent of the nonlinearities involved in the strain-displacement relations. A finite element formulation is presented and various examples of nonlinear shell dynamics including large overall and chaotic motions are considered.     

Hybrid schemes to compute contact discontinuities in Euler equations with any EOSSchémas hybrides pour la simulation des discontinuités de contact des équations d'Euler avec une loi d'état quelconque

We propose here some explicit hybrid schemes which enable accurate computation of Euler equations with arbitrary (analytic or tabulated) equation of state (EOS). The method is valid for the exact Godunov scheme and some approximate Godunov schemes. To cite this article: T. Gallouët et al., C. R. Mecanique 330 (2002) 445–450.     

Transient shock wave flows in tubes with a sudden change in cross section

This paper describes propagation of shock waves within circular cross-section shock tubes with a sudden area change in cross section. A dispersion-controlled scheme was used to solve the Euler equations assuming axisymmetric flows. For experimental visualizations an aspheric cylindrical test section was designed to keep collimated incident light rays parallel once they were reflected or refracted on the inner and outer surfaces of the test section. For effective comparisons with experimental results, equivalent numerical interferograms were constructed to demonstrate effectiveness of the numerical method and verify the observed shock-wave phenomena. The numerical method was used to calculate three further cases with variations of the initial shock-wave Mach number and the flow geometry to clarify the role of these parameters. Complex transient shock-wave phenomena, such as shock-wave reflection, shock/vortex interaction and shock-wave focusing were observed in these cases, and interpreted with shock wave theory. In addition, the research clearly shows that combination of CFD with experiments is effective to highlight physical phenomena in axisymmetric flows. Received 15 June 1996 / Accepted 20 December 1996     

A fourth order accurate finite volume scheme for numerical simulations of turbulent Rayleigh–Bénard convection in cylindrical containers

A finite volume scheme, which is based on fourth order accurate central differences in spatial directions and on a hybrid explicit/semi-implicit time stepping scheme, was developed to solve the incompressible Navier–Stokes and energy equations on cylindrical staggered grids. This includes a new fourth order accurate discretization of the velocity and temperature fields at the singularity of the cylindrical coordinate system and a new stability condition [J. Appl. Numer. Anal. Comput. Math. 1 (2004) 315–326]. The method was applied in direct numerical simulations of turbulent Rayleigh–Bénard convection for different Rayleigh numbers $\mathit{Ra}={10}^{\gamma }$, $\gamma =5\text{,}\dots \text{,}8$, in wide cylinders with the aspect ratios $a\equiv H/R=0.2$ and $a=0.4$ (where R denotes the radius and H – the height of the cylinder). To cite this article: O. Shishkina, C. Wagner, C. R. Mecanique 333 (2005).