Some problems in boundary element programming of axisymmetric elasticity

Axisymmetric problems in elasticity can be reduced to two dimensional ones, but they are a little more complicated than plane problems. Therefore, some special problems will be encountered in the boundary element programming of axisymmetric elasticity. In this paper, the methods to treat these problems and some remarks are given according to our experience in programming. Numerical examples are presented for the checking of these treatments.     

Perfectly matched layer boundary condition for two-dimensional Euler equations in generalized coordinate system

In the past, perfectly matched layer (PML) equations have been constructed in Cartesian and spherical coordinates. In this article, the focus is on the development of a PML absorbing technique for treating numerical boundaries, especially those with unbounded domains, in a generalized coordinate system for a flow in an arbitrary direction. The PML equations for two-dimensional Euler equations are developed in split form through a space–time transformation involving a complex variable transformation with the application of a pseudo-mean-flow in the PML domain. A numerical solver is developed using conventional numerical schemes without employing any form of filtering or artificial dissipation to solve the governing PML equations for two-dimensional Euler equations in a generalized coordinate system. Physical domains of arbitrary shapes are considered and numerical simulations are carried out to validate and demonstrate the effectiveness of the PML as an absorbing boundary condition in generalized coordinates.     

A note on characteristic decomposition for two-dimensional Euler system in van der Waals fluids

This paper is devoted to the generalization of a well-known result on reducible equations by Courant and Friedrichs [7] and a motivational result on compressible Euler system within the context of ideal gases by Li et al. [10]. The characteristic decomposition technique has been used to prove that any hyperbolic state, adjacent to a constant state, is simple for a pseudo-steady isentropic irrotational flow, modeled by Euler equations, in van der Waals fluids. Furthermore, this result is extended to full Euler system in self-similar coordinates provided the pseudo-flow characteristics are extending into a constant state.     

Geometric non-linear analysis of space frame with rotation greater than 90°, with Euler angles and quasi-fixed local axes system

The geometric non-linear analysis of space frames with definition of the local axes through rotations similar to Euler angles, presents difficulties when the rotations reach 90° in the fixed local axes. To surmount these difficulties this work presents a new technique using a system coined quasi-fixed local axes. In this technique, the direction cosines of these axes of the actual increment of loading are defined through the direction cosines of the moving local axes of the last iteration of the previous increment and of the quasi-fixed local axes of the previous increment (or fixed, in case of the first increment). The technique of definition of the quasi-fixed local axes is used with updated Lagrangian reference. The numerical examples presented herein show the great numerical stability of the technique.     

An Euler system source term that develops prototype Z‐pinch implosions intended for the evaluation of shock‐hydro methods

In this paper, a phenomenological model for a magnetic drive source term for the momentum and total energy equations of the Euler system is described. This body force term is designed to produce a Z‐pinch like implosion that can be used in the development and evaluation of shock‐hydrodynamics algorithms that are intended to be used in Z‐pinch simulations. The model uses a J × B Lorentz force, motivated by a 0‐D analysis of a thin shell (or liner implosion), as a source term in the equations and allows for arbitrary current drives to be simulated. An extension that would include the multi‐physics aspects of a proposed combined radiation hydrodynamics (rad‐hydro) capability is also discussed. The specific class of prototype problems that are developed is intended to illustrate aspects of liner implosions into a near vacuum and with idealized pre‐fill plasma effects. In this work, a high‐resolution flux‐corrected‐transport method implemented on structured overlapping meshes is used to demonstrate the application of such a model to these idealized shock‐hydrodynamic studies. The presented results include an asymptotic solution based on a limiting‐case thin‐shell analytical approximation in both (x, y) and (r, z). Additionally, a set of more realistic implosion problems that include density profiles approximating plasma pre‐fill and a set of perturbed liner geometries that excite a hydro‐magnetic like Rayleigh–Taylor instability in the implosion dynamics are demonstrated. Finally, as a demonstration of including and evaluating multiphysics effects in the Euler system, a simple radiation model is included and self‐convergence results for two types of (r, z) implosions are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

亏损振系广义状态向量灵敏度的移频算法

由于亏损系统的状态向量系线性相关,故无法实现与非亏损系统类似的解耦及控制,这为亏损系统灵敏度分析带来很大困难。针对这个问题本文摒弃一般状态向量系,首先引入广义状态向量及其伴随向量理论,利用它们具有的良好规范正交关系,实现了灵敏度控制方程解耦,其次使用移频技术,实现了亏损系统状态矩阵的广义状态向量的灵敏度分析。该方法可用于求解亏损系统各频段广义特征向量的灵敏度,算法结构紧凑,易于理解和实施。数值算例结果表明,该算法是可行和有效的。     

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

A fractional step scheme with modified characteristic finite differences running in a parallel arithmetic is presented to simulate a nonlinear percolation system of multilayer dynamics of fluids in a porous medium with moving boundary values. With the help of theoretical techniques including the change of regions, piecewise threefold quadratic interpolation, calculus of variations, multiplicative commutation rule of difference operators, multiplicative commutation rule of difference operators, decomposition of high order difference operators, induction hypothesis, and prior estimates, an optimal order in l ² norm is displayed to complete the convergence analysis of the numerical algorithm. Some numerical results arising in the actual simulation of migration-accumulation of oil resources by this method are listed in the last section.