A possible counterexample to well posedness of entropy solutions and to Godunov scheme convergence

A particular case of initial data for the two-dimensional Euler equations is studied numerically. The results show that the Godunov method does not always converge to the physical solution, at least not on feasible grids. Moreover, they suggest that entropy solutions (in the weak entropy inequality sense) are not well posed.

     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Geometric Characterizations for Subgroups of PU（1, n; C）

In this paper, we discuss non-elementary subgroups and discontinuous subgroups of PU(1,n; C), and give their geometric characterizations.     

Error estimation,adaptivity and data transfer in enriched plasticity continua to analysis of shear band localization

In this paper, an adaptive FE analysis is presented based on error estimation, adaptive mesh refinement and data transfer for enriched plasticity continua in the modelling of strain localization. As the classical continuum models suffer from pathological mesh-dependence in the strain softening models, the governing equations are regularized by adding rotational degrees-of-freedom to the conventional degrees-of-freedom. Adaptive strategy using element elongation is applied to compute the distribution of required element size using the estimated error distribution. Once a new mesh is generated, state variables and history-dependent variables are mapped from the old finite element mesh to the new one. In order to transfer the history-dependent variables from the old to new mesh, the values of internal variables available at Gauss point are first projected at nodes of old mesh, then the values of the old nodes are transferred to the nodes of new mesh and finally, the values at Gauss points of new elements are determined with respect to nodal values of the new mesh. Finally, the efficiency of the proposed model and computational algorithms is demonstrated by several numerical examples.     

A free boundary problem related to thermal insulation

We study a free boundary problem arising from the theory of thermal insulation. The outstanding feature of this set optimization problem is that the boundary of the set being optimized is not a level surface of a harmonic function, but rather a hypersurface along which a harmonic function satisfies a Robin condition. We show that minimal sets exist, satisfy uniform density estimates, and, under some geometric conditions, have “locally flat” boundaries.     

Kinematics of layered reinforced-concrete planar beam finite elements with embedded transversal cracking

In this work crack formation and development is addressed and implemented in a planar layered reinforced-concrete beam element. The crack initiation and growth is described using the strength criterion in conjunction with exact kinematics of the interlayer connection. In this way a novel embedded-discontinuity beam finite element is derived in which the tensile stresses in concrete at the crack position reaching the tensile strength will trigger a crack to open. Since the element is multi-layered, in this way the crack is allowed to propagate through the depth of the beam. The cracked layer(s) will involve discontinuity in the cross-sectional rotation equal to the crack-profile angle, as well as a discontinuity in the position vector of the layer’s reference line. A bond–slip relationship is superimposed onto this model in a kinematically consistent manner with reinforcement being treated as an additional layer of zero thickness with its own material parameters and a constitutive law implemented in the multi-layered beam element.     

Three-dimensional interfacial fracture analysis of a one-dimensional hexagonal quasicrystal coating

In this paper, the three-dimensional (3D) interfacial fracture is analyzed in a one-dimensional (1D) hexagonal quasicrystal (QC) coating structure under mechanical loading. A planar interface crack with arbitrary shape is studied by a displacement discontinuity method. Fundamental solutions of interfacial concentrated displacement discontinuities are obtained by the Hankel transform technique, and the corresponding boundary integral-differential equations are constructed with the superposition principle. Green’s functions of constant interfacial displacement discontinuities within a rectangular element are derived, and a boundary element method is proposed for numerical simulation. The singularity of stresses near the crack front is investigated, and the stress intensity factors (SIFs) as well as energy release rates (ERRs) are determined. Finally, relevant influencing factors on the fracture behavior are discussed.     

Third‐order WENO scheme with a new smoothness indicator

In this article, we have devised a new reference smoothness indicator for third‐order weighted essentially non‐oscillatory (WENO) scheme to achieve desired order of convergence at critical points. In the context of the weighted essentially non‐oscillatory scheme, reference smoothness indicator is constructed in such a way that it satisfies the sufficient condition on the weights for the third‐order convergence. The goal is to construct a reference smoothness indicator such that the resulted scheme have to achieve the required order of accuracy even if the first two derivatives vanish but not the third derivative. The construction of such reference smoothness indicator is not possible through a linear combination of local smoothness indicators only. We have proposed a reference smoothness indicator to be of the fourth order of accuracy on three‐point stencil that contains the linear combination of the first derivative information of the local and global stencils. The performance enhancement of the WENO scheme through this reference smoothness indicator is verified through the standard numerical experiments. Numerical results indicate that the new scheme provides better results in comparison with the earlier third‐order WENO schemes like WENO‐JS and WENO‐Z. Copyright © 2017 John Wiley & Sons, Ltd.     

An efficient low‐dissipative WENO filter scheme

This paper presents a simple and efficient procedure developed for tracing discontinuities in flow fields. Numerical experiments are carried out to test the new sensor coupled with the associated nonlinear WENO dissipation filter developed to suppress the numerical dissipation. The tests show that, for a problem containing shocks and vortices, the implementation of the new sensor and an optimized WENO scheme can obtain a stable and highly resolved solution. The numerical experiments demonstrated that the new filter scheme performs efficiently both in parallel and serial running for one‐dimensional inviscid flow problems. Direct numerical simulation of a Mach 5 turbulent boundary layer over a flat plate was carried out to demonstrate the applicability of the scheme to the DNS practices. Copyright © 2011 John Wiley & Sons, Ltd.     

孔洞裂纹问题的一种数值方法

研究孔洞与裂纹的相互作用问题,通过把适于单一裂纹的Bueckner原理扩充到含有多孔洞多裂纹的一般体系,将原问题分解为承受远处载荷不含裂纹不含孔洞的均匀问题,和在远处不承受载荷但在裂纹面上和孔洞表面上承受面力的多孔洞多裂纹问题.于是,以应力强度因子作为参量的问题可以通过考虑后者来解决,而利用笔者提出的杂交位移不连续法,这种多孔洞多裂纹问题是容易数值求解的.算例说明该数值方法对分析平面弹性介质中孔洞与裂纹的相互作用既简单又有效.