A theory of finite tensile deformation of double-network hydrogels

A continuum damage model was developed to describe the finite tensile deformation of tough double-network (DN) hydrogels synthesized by polymerization of a water-soluble monomer inside a highly crosslinked rigid polyelectrolyte network. Damage evolution in DN hydrogels was characterized by performing loading-unloading tensile tests and oscillatory shear rheometry on DN hydrogels synthesized from 3-sulfopropyl acrylate potassium salt (SAPS) and acrylamide (AAm). The model can explain all the mechanical features of finite tensile deformation of DN hydrogels, including idealized Mullins effect and permanent set observed after unloading, qualitatively and quantitatively. The constitutive equation can describe the finite elasto-plastic tensile behavior of DN hydrogels without resorting to a yield function. It was showed that tensile mechanics of DN hydrogels in the model is controlled by two material parameters which are related to the elastic moduli of first and second networks. In effect, the ratio of these two parameters is a dimensionless number that controls the behavior of material. The model can capture the stable branch of material response during neck propagation where engineering stress becomes constant. Consistent with experimental data, by increasing the elastic modulus of the second network the finite tensile behavior of the DN hydrogel changes from necking to strain hardening.     

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.     

Versatile anisotropic mesh adaptation methodology applied to pure quantity of interest error estimator. Steady,laminar incompressible flow

We introduce a new flexible mesh adaptation approach to efficiently compute a quantity of interest by the finite element method. Efficiently, we mean that the method provides an evaluation of that quantity up to a predetermined accuracy at a lower computational cost than other classical methods. The central pillar of the method is our scalar error estimator based on sensitivities of the quantity of interest to the residuals. These sensitivities result from the computation of a continuous adjoint problem. The mesh adaptation strategy can drive anisotropic mesh adaptation from a general scalar error contribution of each element. The full potential of our error estimator is then reached. The proposed method is validated by evaluating the lift, the drag, and the hydraulic losses on a 2D benchmark case: the flow around a cylinder at a Reynolds number of 20.     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

A framework to perform high-order deconvolution for finite-volume method on simplicial meshes

In this article, a new framework to design high-order approximations in the context of node-centered finite volumes on simplicial meshes is proposed. The major novelty of this method is that it relies on very simple and compact differential operators, which is a critical point to achieve good performances in the High-performance computing context. This method is based on deconvolution between nodal and volume-average values, which can be conducted to any order. The interest of the new method is illustrated through three different applications: mesh-to-mesh interpolation, levelset curvature computation, and numerical scheme for convection. Higher order can also be achieved within the present framework by introducing high-rank tensors. Although these tensors feature much symmetries, their manipulation can quickly become an overwhelming task. For this reason and without loss of generality, the present articles are limited to third-order expansion. This method, although tightly connected to the k-exact schemes theory, does not rely on successive corrections: the high-order property is obtained in a single operation, which makes them more attractive in terms of performances.     

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.     

Synthesis,self-assembly and photo-responsive behavior of AB2 shaped amphiphilic azo block copolymer

In this work, we report the synthesis of AB2 shaped amphiphilic azo block copolymer by macromolecular azo coupling reaction. The product and intermediates were characterized by various methods. The selfassembly in selected solvents and photo-responsive behavior of the copolymer were studied by means of UV–vis spectrophotometry and TEM. Spherical aggregates can be obtained by gradually adding water into the solution of this amphiphilic azo block copolymer. Upon irradiation with polarized UV(365 nm)light, the aggregates can be elongated in the polarized direction.     

Variable spaced particle in mesh-free method to handle wave-floating body interactions

In this work, the motion of a two-dimensional rectangular freely floating body under waves is simulated using Improved Meshless Local Petrov-Galerkin method with Rankine Source function (IMLPG_R) with variable spacing resolutions. The IMLPG_R method is a particle method that solves Navier–Stokes equations using the fractional step method to capture the wave properties. However, many existing particle methods are computationally intensive to model the wave-floating body due to the requirement of fine particles, needing uniform distribution throughout the domain. To improve the computational efficiency and capture the body response properly, variable spaced particle distribution with fine resolution near the floating body and coarse resolution far from the body is implemented. Numerical schemes to handle variable resolutions are reported. An iterative scheme to handle the wave-floating body is implemented in the particle method. Two test cases, one with small wave and another with steep waves, are simulated for uniform particle distribution and the result shows good agreement with literature. Based on this, the performance of the variable spaced particle distribution is tested in coupling with floating body solver. The application of the method for wave impact load from the green water loading of the floating structure is also simulated.