Direct numerical simulation of the interaction of isotropic turbulence with a shock wave using shock-fitting

The interaction of three-dimensional isotropic turbulence with a plane shock at Mach numbers of $M=2.0$ and $M=3.0$ is investigated via direct numerical simulation. The numerical scheme is based on a characteristic-type formulation of the Navier–Stokes equations and uses fifth-order upwind schemes in space, a fourth order Runge Kutta scheme in time and a shock-fitting as inlet condition. The isotropic turbulence was generated in a separate computation based on a prescribed energy spectrum. This turbulent flow is considered as frozen, and is convected through the shock with a prescribed average shock speed. An FFT interpolation is used to obtain the upstream values at the instantaneous shock location. Turbulence enhancement is observed, and the evolution of velocity fluctuations as well as turbulence microscales are in good agreement with the behaviour observed using shock-capturing. To cite this article: J. Sesterhenn et al., C. R. Mecanique 333 (2005).     

Non-compact Groups,Coherent States,Relativistic Wave Equations and the Harmonic Oscillator

Relativistic geometrical action for a quantum particle in the superspace is analyzed from theoretical group point of view. To this end an alternative technique of quantization outlined by the authors in a previous work and that is based in the correct interpretation of the square root Hamiltonian, is used. The obtained spectrum of physical states and the Fock construction consist of Squeezed States which correspond to the representations with the lowest weights and with four possible (non-trivial) fractional representations for the group decomposition of the spin structure. From the theory of semigroups the analytical representation of the radical operator in the superspace is constructed, the conserved currents are computed and a new relativistic wave equation is proposed and explicitly solved for the time-dependent case. The relation with the Relativistic Schr?dinger equation and the Time-dependent Harmonic Oscillator is analyzed and discussed.      

Simultaneous densitometric determination of anthelmintic drug albendazole and its metabolite albendazole sulfoxide by HPTLC in human plasma and pharmaceutical formulations

A new, simple, accurate and precise high‐performance thin‐layer chromatographic method has been developed and validated for simultaneous determination of an anthelmintic drug, albendazole, and its active metabolite albendazole, sulfoxide. Planar chromatographic separation was performed on aluminum‐backed layer of silica gel 60G F₂₅₄ using a mixture of toluene–acetonitrile–glacial acetic acid (7.0:2.9:0.1, v /v/v) as the mobile phase. For quantitation, the separated spots were scanned densitometrically at 225 nm. The retention factors (R _f) obtained under the established conditions were 0.76 ± 0.01 and 0.50 ± 0.01 and the regression plots were linear (r ² ≥ 0.9997) in the concentration ranges 50–350 and 100–700 ng/band for albendazole and albendazole sulfoxide, respectively. The method was validated for linearity, specificity, accuracy (recovery) and precision, repeatability, stability and robustness. The limit of detection and limit of quantitation found were 9.84 and 29.81 ng/band for albendazole and 21.60 and 65.45 ng/band for albendazole sulfoxide, respectively. For plasma samples, solid‐phase extraction of analytes yielded mean extraction recoveries of 87.59 and 87.13% for albendazole and albendazole sulfoxide, respectively. The method was successfully applied for the analysis of albendazole in pharmaceutical formulations with accuracy ≥99.32%.     

A continuum model accounting for defect and mass densities in solids with inelastic material behaviour

In this paper the analysis of structures with inelastic material behaviour is considered taking into account the evolution of defects and changes in mass density. The underlying kinematical concept of an oriented continuum is general enough to describe the micro- and macrobehaviour of material bodies appropriately. Based on the logical and consistent variational arguments for a Lagrangian functional the dynamic balance laws, boundary and transversality conditions, all related to the evolution of defect density and mass changes, are derived for macro- and microstresses of deformational as well as of configurational type. The adopted procedure, which formally leaves the balance laws unaltered, leads to the additional balance law for changes in defect density and additional boundary conditions for the changes in mass and defect densities. Driving forces or affinities, associated with the evolution of defect and mass densities, and a generalization of the J-integral representing the thermodynamic forces on defects are obtained. A nonlocal constitutive model accounting for changes in the defect density is presented.     

AN ADAPTIVE UNSTRUCTURED TRI-TREE ITERATIVE SOLVER FOR MIXED FINITE ELEMENT FORMULATION OF THE STOKES EQUATIONS

An iterative adaptive equation solver for solving the implicit Stokes equations simultaneously with tri-tree grid generation is developed. The tri-tree grid generator builds a hierarchical grid structure which is mapped to a finite element grid at each hierarchical level. For each hierarchical finite element grid the Stokes equations are solved. The approximate solution at each level is projected onto the next finer grid and used as a start vector for the iterative equation solver at the finer level. When the finest grid is reached, the equation solver is iterated until a tolerated solution is reached. In order to reduce the overall work, the element matrices are integrated analytically beforehand. The efficiency and behaviour of the present adaptive method are compared with those of the previously developed iterative equation solver which is preconditioned by incomplete LU factorization with coupled node fill-in. The efficiency of the incomplete coupled node fill-in preconditioner is shown to be largely dependent on the global node numbering. The preconditioner is therefore tested for the natural node ordering of the tri-tree grid generator and for different ways of sorting the nodes.     

DYNAMICS ANALYSIS OF STOCHASTIC SPATIAL FLEXIBLE MULTIBODY SYSTEM 1)

轻质、高精度的柔性多体系统被广泛应用于实际工程领域中.由于实际设计公差、制造误差及环境温度等多种不确定因素的存在,使得柔性多体系统的结构参数(物理参数和几何参数)表现出随机性.具有随机结构参数的动力学模型能够客观地反映出真实系统的动力学行为,且结构参数的不确定性对空间柔性多体系统动力学响应的影响是不容忽视的.针对具有多个随机参数的空间柔性多体系统,提出了一种基于广义alpha算法的非侵入式随机柔性多体系统动力学计算方法.采用绝对节点坐标公式(absolute node coordinate formulation, ANCF)来描述柔性体, 推导建立多体系统动力学模型.利用混沌多项式展开(polynomial chaos expansion, PCE)法构建系统随机动力学方程的代理模型,然后将随机响应面法(stochastic response surface method, SRSM)嵌入广义-alpha方法中,分别采用改进抽样的回归方法(regression method of improved sampling, RMIS)和单项求容积法则(Monte Carlo simulation, MCR)来确定样本点.将数值计算结果与蒙特卡洛模拟(Monte Carlo simulation, MCS)结果进行对比, 验证了所提算法的有效性.在相同的定积分精度的条件下,根据单项求容积法则确定的样本点的计算结果稳定性更强, 且其计算效率更高.     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.     

Global structural behaviour of thin and moderately thick monoclinic spherical shells using a mixed shear deformation model

Summary The static and dynamic responses of anisotropic spherical shells under a uniformly distributed transverse load are investigated. Analytical solutions using the mixed variational formulation are presented for spherical shells subjected to various boundary conditions. Numerical results of a refined mixed first-order shear deformation theory for natural frequencies, critical buckling, center deflections and stresses are compared with those obtained using the classical shell theory. A variety of simply-supported and clamped boundary conditions are considered and comparisons with the existing literature are made. The sample numerical results presented herein for global structural behaviour of monoclinic spherical shells should serve as references for future comparisons.     

Implicit and explicit integration in the solution of the absolute nodal coordinate differential/algebraic equations

This investigation is concerned with the use of an implicit integration method with adjustable numerical damping properties in the simulation of flexible multibody systems. The flexible bodies in the system are modeled using the finite element absolute nodal coordinate formulation (ANCF), which can be used in the simulation of large deformations and rotations of flexible bodies. This formulation, when used with the general continuum mechanics theory, leads to displacement modes, such as Poisson modes, that couple the cross section deformations, and bending and extension of structural elements such as beams. While these modes can be significant in the case of large deformations, and they have no significant effect on the CPU time for very flexible bodies; in the case of thin and stiff structures, the ANCF coupled deformation modes can be associated with very high frequencies that can be a source of numerical problems when explicit integration methods are used. The implicit integration method used in this investigation is the Hilber–Hughes–Taylor method applied in the context of Index 3 differential-algebraic equations (HHT-I3). The results obtained using this integration method are compared with the results obtained using an explicit Adams-predictor-corrector method, which has no adjustable numerical damping. Numerical examples that include bodies with different degrees of flexibility are solved in order to examine the performance of the HHT-I3 implicit integration method when the finite element absolute nodal coordinate formulation is used. The results obtained in this study show that for very flexible structures there is no significant difference in accuracy and CPU time between the solutions obtained using the implicit and explicit integrators. As the stiffness increases, the effect of some ANCF coupled deformation modes becomes more significant, leading to a stiff system of equations. The resulting high frequencies are filtered out when the HHT-I3 integrator is used due to its numerical damping properties. The results of this study also show that the CPU time associated with the HHT-I3 integrator does not change significantly when the stiffness of the bodies increases, while in the case of the explicit Adams method the CPU time increases exponentially. The fundamental differences between the solution procedures used with the implicit and explicit integrations are also discussed in this paper.     

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.