A new volume‐preserving and continuous interface reconstruction method for 2D multi‐material flow

A new two‐dimensional interface reconstruction method that ensures continuity of the interface and preserves volume fractions is presented here. It is made of two steps, first, the minimization of a cost functional based on volume fractions least square errors by using dynamic programming, a fast and efficient scheme well known in image processing, and then a local correction phase. In each cell, the interface is made of two line segments joining two edges. This new interface reconstruction method, called Dynamic Programming Interface Reconstruction has been coupled with various advection schemes, among them the Lagrange + remap scheme. With a reasonable computational cost, it has been observed in various test cases that Dynamic Programming Interface Reconstruction is more accurate and less diffusive compared with other existing classical reconstruction methods. Copyright © 2017 John Wiley & Sons, Ltd.     

膨胀球面运动时的形状问题

算出了膨胀球面运动时的形状,从而说明光的波前为什么不因为长度收缩而成为椭球面.     

Implicit large eddy simulation using the high‐order correction procedure via reconstruction scheme

We investigate implicit large eddy simulation of the Taylor–Green vortex, Comte‐Bellot–Corrsin experiment, turbulent channel flow and transitional and turbulent flow over an SD7003 airfoil using the high‐order unstructured correction procedure via reconstruction (CPR) scheme, also known as the flux reconstruction scheme. We employ P1 (second‐order) to P5 (sixth‐order) spatial discretizations. Results show that the CPR scheme can accurately predict turbulent flows without the addition of a sub‐grid scale model. Numerical dissipation, concentrated at the smallest resolved scales, is found to filter high‐frequency content from the solution. In addition, the high‐order schemes are found to be more accurate than the low‐order schemes on a per degree of freedom basis for the canonical test cases we consider. These results motivate the further investigation and use of the CPR scheme for simulating turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Hybrid approach for two dimensional damage localization using piezoelectric smart aggregates

In the paper, a novel approach for damage localization in reinforced concrete plates, based on the computational analysis of piezoelectric smart aggregates, has been presented. The hybrid approach for damage localization is based on two criteria: wave propagation energy and time of flight. The comprehensive numerical analysis using standard and explicit finite element method has been conducted. In addition, the proposed algorithm of the hybrid method has been coded in MATLAB. The approach has been verified numerically using different square reinforced concrete plate models, considering different number, position and size of damage, as well as different number and position of the piezoelectric smart aggregates. Obtained results confirm the successful application of the novel approach to the damage localization.     

On the relationship between the dynamic behavior and nanoscale staggered structure of the bone

Bone, a typical load-bearing biological material, composed of ordinary base materials such as organic protein and inorganic mineral arranged in a hierarchical architecture, exhibits extraordinary mechanical properties. Up to now, most of previous studies focused on its mechanical properties under static loading. However, failure of the bone occurs often under dynamic loading. An interesting question is: Are the structural sizes and layouts of the bone related or even adapted to the functionalities demanded by its dynamic performance? In the present work, systematic finite element analysis was performed on the dynamic response of nanoscale bone structures under dynamic loading. It was found that for a fixed mineral volume fraction and unit cell area, there exists a nanoscale staggered structure at some specific feature size and layout which exhibits the fastest attenuation of stress waves. Remarkably, these specific feature sizes and layouts are in excellent agreement with those experimentally observed in the bone at the same scale, indicating that the structural size and layout of the bone at the nanoscale are evolutionarily adapted to its dynamic behavior. The present work points out the importance of dynamic effect on the biological evolution of load-bearing biological materials.     

Dynamic bulk and shear moduli due to grain-scale local fluid flow in fluid-saturated cracked poroelastic rocks: Theoretical model

Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori–Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.     

Simulation of elastic wave propagation using cellular automata and peridynamics,and comparison with experiments

Peridynamics is a non-local continuum mechanics formulation that can handle spatial discontinuities as the governing equations are integro-differential equations which do not involve gradients such as strains and deformation rates. This paper employs bond-based peridynamics. Cellular Automata is a local computational method which, in its rectangular variant on interior domains, is mathematically equivalent to the central difference finite difference method. However, cellular automata does not require the derivation of the governing partial differential equations and provides for common boundary conditions based on physical reasoning. Both methodologies are used to solve a half-space subjected to a normal load, known as Lamb’s Problem. The results are compared with theoretical solution from classical elasticity and experimental results. This paper is used to validate our implementation of these methods.     

Finite element modeling of micropolar-based phononic crystals

The performance of a Cosserat/micropolar solid as a numerical vehicle to represent dispersive media is explored. The study is conducted using the finite element method with emphasis on Hermiticity, positive definiteness, principle of virtual work and Bloch–Floquet boundary conditions. The periodic boundary conditions are given for both translational and rotational degrees of freedom and for the associated force- and couple-traction vectors. Results in terms of band structures for different material cells and mechanical parameters are provided.     

An efficient front‐tracking method for simulation of multi‐density bubbles

In this article, a masked bubble strategy is proposed using the front‐tracking method when simulation of multi‐density bubbles to reduce remarkably the computational cost from both the RAM usage and the number of computations at each time step comparing with the regular method. In the masked bubble strategy, instead of using full domain to update the properties at each time step, each bubble is considered as enclosed in the smallest box required to compute the properties based on the Peskin's function, which needs at least two full mesh sizes from both sides of the interface of each bubble in any directions. To show the performance of the masked bubble strategy in the front‐tracking method, we study the multi‐density bubbles motion in a curved duct flow induced by a pressure gradient in the absence of gravity. To solve the density Poisson equation, the parallel direct solver scheme is tested. The comparison of numerical simulations at the same conditions indicates that the parallel direct solver scheme under the masked bubble strategy considerably reduces the computational time and RAM usage relative to the regular full‐domain method, providing using simulations on finer grid resolutions. Copyright © 2016 John Wiley & Sons, Ltd.