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.     

Prediction of the drag reduction effect of pulsating pipe flow based on machine learning

Prediction of drag reduction effect caused by pulsating pipe flows is examined using machine learning. First, a large set of flow field data is obtained experimentally by measuring turbulent pipe flows with various pulsation patterns. Consequently, more than 7000 waveforms are applied, obtaining a maximum drag reduction rate and maximum energy saving rate of 38.6% and 31.4%, respectively. The results indicate that the pulsating flow effect can be characterized by the pulsation period and pressure gradient during acceleration and deceleration. Subsequently, two machine learning models are tested to predict the drag reduction rate. The results confirm that the machine learning model developed for predicting the time variation of the flow velocity and differential pressure with respect to the pump voltage can accurately predict the nonlinearity of pressure gradients. Therefore, using this model, the drag reduction effect can be estimated with high accuracy.     

Pattern transition and regulation in a subthalamopallidal network under electromagnetic effect

Although the significant roles of magnetic induction and electromagnetic radiation in the neural system have been widely studied, their influence on Parkinson's disease (PD) has yet to be well explored. By virtue of the magnetic flux variable, this paper studies the transition of firing patterns induced by magnetic induction and the regulation effect of external magnetic radiation on the firing activities of the subthalamopallidal network in basal ganglia. We find: (i) The network reproduces five typical waveforms corresponding to the severity of symptoms: weak cluster, episodic, continuous cluster, episodic, and continuous wave. (ii) Magnetic induction is a double-edged sword for the treatment of PD. Although the increase of magnetic coefficient may lead the physiological firing activity to transfer to pathological firing activity, it also can regulate the pathological intensity firing activity with excessive β-band power transferring to the physiological firing pattern with weak β-band power. (iii) External magnetic radiation could inhibit continuous tremulous firing and β-band power of subthalamic nucleus (STN), which means the severity of symptoms weakened. Especially, the bi-parameter plane of the regulation region shows that a short pulse period of magnetic radiation and a medium level of pulse percentage can well regulate pathological oscillation. This work helps to understand the firing activity of the subthalamopallidal network under electromagnetic effect. It may also provide insights into the mechanisms behind the electromagnetic therapy of PD-related firing activity.     

Viscoelastic properties of ultrathin polycarbonate films by liquid dewetting

A liquid dewetting method for the determination of the viscoelastic properties of ultrathin polymer films has been extended to study thickness effects on the properties of ultrathin polycarbonate (PC) films. PC films with film thicknesses ranging from 4 to 299 nm were placed on glycerol at temperatures from below the macroscopic glass transition temperature (T_g) to above it with the dewetting responses being monitored. It is found that the isothermal creep results for films of the same thickness, but dewetted at different temperatures can be superposed into one master curve, which is consistent with the fact of PC being a thermorheologically simple material. Furthermore, the results show that the T_g of PC thin films is thickness dependent, but the dependence is weaker than the results for freely standing films and similar to literature data for PC films supported on rigid substrates. It was also found that the rubbery plateau region for the PC films stiffens dramatically, but still less than what has been observed for freely standing polycarbonate films. The rubbery stiffening is discussed in terms of a recently reported model that relates macroscopic segmental dynamics with the stiffening. © 2015 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 1559–1566     

Transition-Metal Catalysis of Triene 6π Electrocyclization: The π-Complexation Strategy Realized

Triene 6π electrocyclization, wherein a conjugated triene undergoes a concerted stereospecific cycloisomerization to a cyclohexadiene, is a reaction of great historical and practical significance. In order to circumvent limitations imposed by the normally harsh reaction conditions, chemists have long sought to develop catalytic variants based upon the activating power of metal–alkene coordination. Herein, we demonstrate the first successful implementation of such a strategy by utilizing [(C₅H₅)Ru(NCMe)₃]PF₆ as a precatalyst for the disrotatory 6π electrocyclization of highly substituted trienes that are resistant to thermal cyclization. Mechanistic and computational studies implicate hexahapto transition-metal coordination as responsible for lowering the energetic barrier to ring closure. This work establishes a foundation for the development of new catalysts for stereoselective electrocyclizations.     

Phase transition in the one-dimensional pair-hopping model with unusual one-electron hopping

By modifying the conventional one-electron hopping behavior, we study effects of an occupation-dependent hopping on the ground state of the half-filled one-dimensional pair-hopping model. At weak coupling, the use of bosonization and renormalization-group analysis techniques helps to derive the phase diagram. Such unusual hopping is shown to drive a spin-gap transition and to introduce a new region where the triplet superconducting instability dominates for positively small pair-hopping interaction.     

Fluid‐structure coupling of linear elastic model with compressible flow models

Cavitation erosion is caused in solids exposed to strong pressure waves developing in an adjacent fluid field. The knowledge of the transient distribution of stresses in the solid is important to understand the cause of damaging by comparisons with breaking points of the material. The modeling of this problem requires the coupling of the models for the fluid and the solid. For this purpose, we use a strategy based on the solution of coupled Riemann problems that has been originally developed for the coupling of 2 fluids. This concept is exemplified for the coupling of a linear elastic structure with an ideal gas. The coupling procedure relies on the solution of a nonlinear equation. Existence and uniqueness of the solution is proven. The coupling conditions are validated by means of quasi‐1D problems for which an explicit solution can be determined. For a more realistic scenario, a 2D application is considered where in a compressible single fluid, a hot gas bubble at low pressure collapses in a cold gas at high pressure near an adjacent structure.     

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.