An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.     

Formation of low-temperature deformation-induced segregations of nickel in Fe–Ni-based austenitic alloys

ABSTRACT

The authors present the results of an investigation in Fe–Ni-Cr austenitic alloys of the low-temperature deformation-induced segregations of nickel that form in the micro regions being (i) located close to grain- and subgrain boundaries and (ii) characteristic of the concentration and magnetic inhomogeneities indicated by the appearance of a dark diffraction contrast at the electron diffraction patterns taken from these regions typical (at the same time) of an enhanced value of Curie temperature. The observed effects were connected with the micro distortions caused by the local change of lattice parameter because of an increase in nickel concentration, as well as in the result of a magnetostriction dilatation. Using methods of the X-ray energy dispersive spectroscopy (XEDS) and atomic-probe body-section radiography (tomography – APT) has made it possible to determine the borders of those regions of austenite that were characteristic of an enhanced concentration of nickel in the fields of the localisation of a deformation-induced segregation of nickel in the vicinity of grain (subgrain) boundaries of austenitic alloys of the types Fe–13Cr–30Ni and Fe–37Ni–3Ti.     

ZnO nanowire lines and bundles: Template-deformation-guided alignment for patterned field-electron emitters

One-dimensional ZnO materials have been promising for field-emission (FE) application, but how to facially control the alignment of ZnO emitters is still a great challenge especially for patterned display application. Here, we report the fabrication of novel ZnO nanowire (NW) line and bundle arrays for patterned field-electron emitters. The effects of PS template size and heating time on the resulted ZnO nanoarrays were systematically studied. The deformation degree of PS templates was controlled and hence utilized to adjust the alignment of electrochemically deposited ZnO arrays. It was found that the length of NW lines and the density of NW bundles can effectively tuned by the PS template heating time. The optimal FE performance with turn-on electric field as low as of 4.4 V μm⁻¹ and the field-enhancement factor as high as of 1450 were achieved through decreasing the screening effect among the patterned field-electron emitters.     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

Quadrilateral surface meshes without self-intersecting dual cycles for hexahedral mesh generation

Several promising approaches for hexahedral mesh generation work as follows: Given a prescribed quadrilateral surface mesh they first build the combinatorial dual of the hexahedral mesh. This dual mesh is converted into the primal hexahedral mesh, and finally embedded and smoothed into the given domain. Two such approaches, the modified whisker weaving algorithm by Folwell and Mitchell, as well as a method proposed by the author, rely on an iterative elimination of certain dual cycles in the surface mesh. An intuitive interpretation of the latter method is that cycle eliminations correspond to complete sheets of hexahedra in the volume mesh.

Although these methods can be shown to work in principle, the quality of the generated meshes heavily relies on the dual cycle structure of the given surface mesh. In particular, it seems that difficulties in the hexahedral meshing process and poor mesh qualities are often due to self-intersecting dual cycles. Unfortunately, all previous work on quadrilateral surface mesh generation has focused on quality issues of the surface mesh alone but has disregarded its suitability for a high-quality extension to a three-dimensional mesh.

In this paper, we develop a new method to generate quadrilateral surface meshes without self-intersecting dual cycles. This method reuses previous b-matching problem formulations of the quadrilateral mesh refinement problem. The key insight is that the b-matching solution can be decomposed into a collection of simple cycles and paths of multiplicity two, and that these cycles and paths can be consistently embedded into the dual surface mesh.

A second tool uses recursive splitting of components into simpler subcomponents by insertion of internal two-manifolds. We show that such a two-manifold can be meshed with quadrilaterals such that the induced dual cycle structure of each subcomponent is free of self-intersections if the original component satisfies this property. Experiments show that we can achieve hexahedral meshes with a good quality.     

双奇核136La的集体带结构

通过重离子核反应与在束γ谱的实验技术, 对A=130缺中子核区的双奇核¹³⁶La的高自旋态进行了研究, 所用核反应为¹³⁰Te(¹¹B,5n). 实验结果扩展了¹³⁶La的能级纲图, 包括3个集体转动带, 最高自旋态达20h. 对于\uppi h_{11/2}\otimes \upnu h_{11/2}$~带, 观测到了旋称反转与集体回弯现象. 通过系统学比较, 对旋称反转特性进行了讨论. 由推转壳模型的计算指出, 此集体回弯是由一对中子的角动量顺排引起的. 另外两个集体带为具有~$\gamma\approx -60^\circ$~的扁椭形变带, 其可能的组态为: $\uppi h_{11/2}\otimes \upnu g_{7/2}h_{11/2}^2$~与~$\uppi g_{7/2}\otimes\upnu g_{7/2}^2 d_{5/2} h_{11/2}^2$.     

奇异摄动罗宾问题在Bakhvalov-Shishkin网格上的一致收敛有限差分法

考虑一个奇异摄动罗宾问题在Bakhvalov-Shishkin网格上的迎风差分策略,得到在改进的Shishkin网格上迎风策略是关于ε一致的一阶L∞模收敛的.数值实验证实了此理论结果,显示估计是稳健的.     

Stability and Nonlinear Deformation of Cylindrical Grids

The general stability of single-layer cylindrical grids is studied in linear and nonlinear formulations. Dependence of the general buckling load on the geometry and stiffness parameters of a grid is established in an analytical form. Such grids are numerically analyzed for stability. It is established that the general buckling load is much less than the local buckling load. Typical general buckling modes are found. It is shown that such grids are weakly sensitive to imperfections     

Numerical Prediction of the Strain State of Thermoset Matrix/Mineral Filler Composite Plates

The numerical prediction of the fields of inelastic strains (the linear invariant of the tensor of inelastic strains) in thermoset polyester/marble filler composite plates is discussed. A uniformly distributed load is applied to the plates, which lie on a steel base. The strain fields are predicted by means of the boundary element method by using creep test data for the composites and the polyester matrix itself. Identical creep tests were performed for two ages of the materials (1 month and 13 years), which allowed evaluating the aging effect. The study is carried out in two stages. At the first stage, the application of the generalized Maxwell-Gurevich equation to the thermoset matrix/mineral filler composite is demonstrated. The model parameters determined from the experimental creep data is used for the second stage, where the state of inelastic strains in the plates is predicted by applying the boundary element method. The influence of composite formulation (filler content) and physical aging of the polyester matrix on the state of inelastic strains in the plates is shown.__________Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 41, No. 2, pp. 145–156, March–April, 2005.