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.     

Stress Distribution in an Infinite Elastic Body with a Locally Curved Fiber in a Geometrically Nonlinear Statement

Within the framework of a piecewise homogeneous body model, with the use of three-dimensional geometrically nonlinear exact equations of elasticity theory, a method for determining the stress—strain state in unidirectional fibrous composites with locally curved fibers is developed for the case where the interaction between the fibers is neglected. All the investigations are carried out for an infinite elastic body containing a single locally curved fiber. Numerical results illustrating the effect of geometrical nonlinearity on the distribution of the self-balanced normal and shear stresses acting on the interface and arising as a result of local curving of the fiber are presented.__________Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 41, No. 4, pp. 433–448, July–August, 2005.     

Numerical analysis of turbulent structure in compound meandering open channel by algebraic Reynolds stress model

The turbulent flow in a compound meandering channel with a rectangular cross section is one of the most complicated turbulent flows, because the flow behaviour is influenced by several kinds of forces, including centrifugal forces, pressure‐driven forces and shear stresses generated by momentum transfer between the main channel and the flood plain. Numerical analysis has been performed for the fully developed turbulent flow in a compound meandering open‐channel flow using an algebraic Reynolds stress model. The boundary‐fitted coordinate system is introduced as a method for coordinate transformation in order to set the boundary conditions along the complicated shape of the meandering open channel. The turbulence model consists of transport equations for turbulent energy and dissipation, in conjunction with an algebraic stress model based on the Reynolds stress transport equations. With reference to the pressure–strain term, we have made use of a modified pressure–strain term. The boundary condition of the fluctuating vertical velocity is set to zero not only for the free surface, but also for computational grid points next to the free surface, because experimental results have shown that the fluctuating vertical velocity approaches zero near the free surface. In order to examine the validity of the present numerical method and the turbulent model, the calculated results are compared with experimental data measured by laser Doppler anemometer. In addition, the compound meandering open channel is clarified somewhat based on the calculated results. As a result of the analysis, the present algebraic Reynolds stress model is shown to be able to reasonably predict the turbulent flow in a compound meandering open channel. Copyright © 2005 John Wiley & Sons, Ltd.     

Raman properties embeddedof GaSb nanoparticlesin SiO2 films

The Raman shifts of nanocrystalline GaSb excited by an Ar^ ion laser at wavelengths 514.5, 496.5, 488.0, 476.5,and 457.9nm are studied by an SPEX-1403 laser Raman spectrometer respectively, and they are explained by phonon confinement, tensile stress, resonant Raman scattering and quantum size effects. The Stokes and anti-Stokes Raman spectra of GaSb nanocrystals strongly support the Raman feature of GaSb nanocrystals. The calculated optical spectra compare well with experimental data on Raman scattering GaSb nanocrystals.     

On a free boundary problem for the Reynolds equation derived from the Stokes system with Tresca boundary conditions

The asymptotic behaviour of a Stokes flow with Tresca free boundary friction conditions when one dimension of the fluid domain tends to zero is studied. A specific Reynolds equation associated with variational inequalities is obtained and uniqueness is proved.     

Closed-form solutions for a rectangular layer with central crack under anti-plane shear stress

We consider the problem of determining the stress distributionin a finite rectangular elastic layer containing a Griffithcrack which is opened by internal shear stress acting alongthe length of the crack. The mode III crack is assumed to belocated in the middle plane of the rectangular layer. The followingtwo problems are considered: (A) the central crack is perpendicularto the two fixed lateral surfaces and parallel to the othertwo stress-free surfaces; (B) all the lateral surfaces of therectangular layer are clamped and the central crack is parallelto the two lateral surfaces. By using Fourier transformations,we reduce the solution of each problem to the solution of dualintegral equations with sine kernels and a weight function whichare solved exactly. Finally, we derive closed-form expressionsfor the stress intensity factor at the tip of the crack andthe numerical values for the stress intensity factor at theedges of the cracks are presented in the form of tables.     

Au(111)上烷烃硫醇SAMs表面应力的理论研究

应用已建立的关于金属表面吸附层中表面应力的统计热力学理论 ,计算了Au(111)上烷烃硫醇SAMs的表面应力及其与烷烃硫醇链长、吸附覆盖度的定量关系 .计算结果与实验相符 ,较好地解释了Berger等人的实验结果 ,特别是解决了在表面应力符号性质上理论与实验的矛盾 .在表面吸附层应力的多种物理起源中 ,通过底物的分子间作用力有着决定性的贡献 ,揭示了分子的吸附能间接地起着重要作用 .这与阴离子化学吸附体系Cl-/Au(111)的有关研究结果相同 .     

The study of the influence of uniaxial stress on impurity complexes in silicon

The influence of external uniaxial stress on the different indium-donor complexes in silicon has been studied using the perturbed γ –γ angular correlation (PAC) method. Such effect of an applied stress is detected by means of the probe atoms situated at different complexes in the sample. The current results showed that the responses of the probes in an extrinsic silicon samples are found to be dissimilar for the same value of stress. Such change in the local environments of the probe atoms could be associated with the various strain field created by the implantations of varied size of impurities. The phosphorous implantation in silicon has even lead to the complete absence of observable effect of the applied stress suggesting significant lose of the elasticity of the sample.     

A numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. Copyright © 2004 John Wiley & Sons, Ltd.