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.     

Oracle Inequalities for Efromovich–Pinsker Blockwise Estimates

Oracle inequality is a relatively new statistical tool for the analysis of nonparametric adaptive estimates. Oracle is a good pseudo-estimate that is based on both data and an underlying estimated curve. An oracle inequality shows how well an adaptive estimator mimics the oracle for a particular underlying curve. The most advanced oracle inequalities have been recently obtained by Cavalier and Tsybakov (2001) for Stein type blockwise estimates used in filtering a signal from a stationary white Gaussian process. The authors also conjecture that a similar result can be obtained for Efromovich–Pinsker (EP) type blockwise estimators where their approach, based on Stein's formula for risk calculation, does not work. This article proves the conjecture and extends it upon more general models which include not stationary and dependent processes. Other possible extensions, a discussion of practical implications and a numerical study are also presented.     

A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

Complex genetic evolution of artificial self‐replicators in cellular automata

It is widely believed that evolutionary dynamics of artificial self‐replicators realized in cellular automata (CA) are limited in diversity and adaptation. Contrary to this view, we show that complex genetic evolution may occur within simple CA. The evolving self‐replicating loops (“evoloops”) we investigate exhibit significant diversity in macro‐scale morphologies and mutational biases, undergoing nontrivial genetic adaptation by maximizing colony density and enhancing sustainability against other species. Nonmutable subsequences enable genetic operations that alter fitness differentials and promote long‐term evolutionary exploration. These results demonstrate a unique example of genetic evolution hierarchically emerging from local interactions between elements much smaller than individual replicators. © 2004 Wiley Periodicals, Inc. Complexity 10: 33–39, 2004     

Development of xylose-fermenting yeasts for ethanol production at high acetic acid concentrations

Mutants resistant to comparatively high levels of acetic acid were isolated from the xylose-fermenting yeastsCandida shehatae andPichia stipitis by adapting these cultures to increasing concentrations of acetic acid grown in shake-flask cultures. These mutants were tested for their ability to ferment xylose in presence of high acetic acid concentrations, in acid hydrolysates of wood, and in hardwood spent sulfite liquor, and compared with their wild-type counterparts and between themselves. TheP. stipitis mutant exhibited faster fermentation times, better tolerance to acid hydrolysates, and tolerance to lower pH.     

Mass shifting roles of negative kernel mass in density estimation

f(x) is a univariate density in C ⁴ with bounded support. For any n and sufficiently small kernel bandwidths, the symmetric appendage of any negative mass, –U, to any smooth unimodal symmetric kernel of order p=2 shifts expected estimator mass from regions where f(x)>0 to regions where f(x)<0. For large n, the mean automatic kernel adaptation induced by –U is analyzed in the simplest MISE reduction scenario: The symmetric appendage of –U to the uniform kernel K(x, X) over MISE-optimal bandwidths reduces MISE by shifting K(x, X) mass asymmetrically across the observation X in the direction of decreasing |f(x)|.     

Adaptation of Cu/ZnO/Al_2O_3 to Temperature Change in Methanol Synthesis from CO_2 Hydrogenation

The induction behavior in CO2 hydrogenation was studied by varying the reaction temperature to investigate the adaptation of the Cu/ZnO/Al2O3 catalyst to the temperature change. The results indicated that a used catalyst had a tendency to keep the last running state in new reaction conditions for MeOH formation, and that this tendency was related to the difference in Cu/Cun+ ratio caused by CO2 and CO produced at different reaction temperatures. However, the reverse water-gas shift reaction (RWGS) induced at four temperatures was completely different from that of methanol synthesis. It implied that the two so-called competitive reactions in CO2+H2, RWGS and methanol synthesis, have different active centers.     

A moving mesh algorithm for aero‐thermo‐mechanical modelling in turbomachinery

This paper describes the development of a mesh deformation method used for aero‐thermo‐mechanical coupling of turbo‐engine components. The method is based on the nonlinear solution of an elastic medium analogy, solved using finite element discretisation and modified to let the boundary nodes be free to slide over the deflected surfaces. This sliding technique relies on a B‐spline reconstruction of the moving boundary and increases the robustness of the method in situations where the boundary deflection field presents significant gradients or large relative motion between two distinct boundaries. The performance of the method is illustrated with the application to an interstage cavity of a turbine assembly, subjected to the deformations computed by a coupled thermo‐mechanical analysis of the engine component. Copyright © 2011 John Wiley & Sons, Ltd.