Theoretical & Applied Mechanics Letters

Online First
Current Issue
Archive
Most Downloaded

Display Method：

Bedding plane-embedded augmented virtual internal bonds for fracture propagation simulation in shale

Zihan Liu, Zhennan Zhang, Ahmad Ghassemi

Accepted Manuscript

[Abstract] (3) [FullText HTML] (2) [PDF 0KB] (0)

Abstract:
To effectively simulate the fracture propagation in shale, the bedding plane (BP) effect is incorporated into the augmented virtual internal bond (AVIB) constitutive relation through BP tensor. Comparing the BP-embedded AVIB with the theory of transverse isotropy, it is found the approach can represent the anisotropic properties induced by parallel BPs. Through the simulation example, it is found that this method can simulate the stiffness anisotropy of shale and can represent the effect of BPs on hydraulic fracture propagation direction. Compared with the BP-embedded VIB, this method can account for the various Poisson's ratio. It provides a feasible approach to simulate the fracture propagation in shale.

Neuroevolution-enabled adaptation of the Jacobi method for Poisson’s equation with density discontinuities

T.-R. Xiang, X. I. A. Yang, Y.-P. Shi

Accepted Manuscript

[Abstract] (21) [FullText HTML] (12) [PDF 0KB] (0)

Abstract:
Lacking labeled examples of working numerical strategies, adapting an iterative solver to accommodate a numerical issue, e.g., density discontinuities in the pressure Poisson equation, is non-trivial and usually involves a lot of trial and error. Here, we resort to evolutionary neural network. A evolutionary neural network observes the outcome of an action and adapts its strategy accordingly. The process requires no labeled data but only a measure of a network’s performance at a task. Applying neuro-evolution and adapting the Jacobi iterative method for the pressure Poisson equation with density discontinuities, we show that the adapted Jacobi method is able to accommodate density discontinuities.

Error of large-eddy simulation in the wall pressure fluctuation of a turbulent channel flow

Rong Li, Bowen Yang, Zixuan Yang, Shizao Wang, Guowei He

Accepted Manuscript

[Abstract] (178) [FullText HTML] (41) [PDF 0KB] (11)

Abstract:
We analyze the error of large-eddy simulation (LES) in wall pressure fluctuation of a turbulent channel flow. To separate different sources of the error, we conduct both direct numerical simulations (DNS) and LES, and apply an explicit filter on DNS data to obtain filtered DNS (FDNS) data. The error of LES is consequently decomposed into two parts: The first part is the error of FDNS with respect to DNS, which quantifies the influence of the filter operation. The second part is the difference between LES and FDNS induced by the error of LES in velocity field. By comparing the root-mean-square value and the wavenumber-frequency spectrum of the wall pressure fluctuation, it is found that the inaccuracy of the velocity fluctuations is the dominant source that induces the error of LES in the wall pressure fluctuation. The present study provides a basis on future LES studies of the wall pressure fluctuation.

Derivation of FFT numerical bounds of the effective properties of composites and polycristals

Minh-Tan Nguyen, Quy-Dong To, Vincent Monchiet

Corrected proof , doi: 10.1016/j.taml.2021.100236

[Abstract] (76) [FullText HTML] (64) [PDF 3017KB] (2)

Abstract:
In this paper, we provide exact fast Fourier transform (FFT)-based numerical bounds for the elastic properties of composites having arbitrary microstructures. Two bounds, an upper and a lower, are derived by considering usual variational principles based on the strain and the stress potentials. The bounds are computed by solving the Lippmann-Schwinger equation together with the shape coefficients which allow an exact description of the microstructure of the composite. These coefficients are the exact Fourier transform of the characteristic functions of the phases. In this study, the geometry of the microstructure is approximated by polygonals (two-dimensional, 2D objects) and by polyhedrons (three-dimensional, 3D objects) for which exact expressions of the shape coefficients are available. Various applications are presented in the paper showing the relevance of the approach. In the first benchmark example, we consider the case of a composite with fibers. The effective elastic coefficients ares derived and compared, considering the exact shape coefficient of the circular inclusion and its approximation with a polygonal. Next, the homogenized elastic coefficients are derived for a composite reinforced by 2D flower-shaped inclusions and with 3D toroidal-shaped inclusions. Finally, the method is applied to polycristals considering Voronoi tessellations for which the description with polygonals and polyhedrons becomes exact. The comparison with the original FFT method of Moulinec and Suquet is provided in order to show the relevance of these numerical bounds.

Supplyment

Accepted Manuscript

[Abstract] (77) [FullText HTML] (34) [PDF 2408KB] (0)

Abstract:
$v.abstractInfoEn

Analytical and Numerical Studies for Seiches in a Closed Basin with Various Geometric Shape of Longitudinal Section

I. Magdalena, Nadhira karima, H. Q. Rif'atin

Accepted Manuscript

[Abstract] (169) [FullText HTML] (41) [PDF 3043KB] (1)

Abstract:
Seiches are long-period standing waves with a unique period called a natural resonant period, during which the phenomenon of resonance occurs. The occurrence of resonance in coastal areas can cause destruction to surrounding natural and man-made structures. By determining the resonant period of a given basin, we can pinpoint the conditions that allow waves to achieve resonance. In this study, a mathematical model is developed from the shallow water equations to examine seiches and resonances in various types of closed basin. The developed model is solved analytically using the separation of variables method to determine the seiches' fundamental resonant periods. Comparisons between the analytical solutions and experimental measurements for resonant periods are also provided. It is shown that the analytical resonant period confirms the experimental data for closed basin of various geometric profiles. Using a finite volume method on a staggered grid, the model is solved numerically to simulate the wave profile when resonance phenomenon occurs. Through those numerical simulations, we also obtain the fundamental resonant period for each basin which agrees with the derived analytical period.

Computing mean fields with known Reynolds stresses at steady state

Xianwen Guo, Zhenhua Xia, Shiyi Chen

Accepted Manuscript

[Abstract] (170) [FullText HTML] (41) [PDF 3354KB] (4)

Abstract:
With the rising of modern data science, data-driven turbulence modeling with the aid of machine learning algorithms is becoming a new promising field. Many approaches are able to achieve better Reynolds stress prediction, with much lower modeling error (

\begin{document}$\epsilon_{\rm{M}}$\end{document}

), than traditional RANS models, but they still suffer from numerical error and stability issues when the mean velocity fields are estimated by solving Reynolds-averaged Navier-Stokes (RANS) equations with the predicted Reynolds stresses. This fact illustrates that the error of solving the RANS equations (

\begin{document}$\epsilon_{\rm{P}}$\end{document}

) is also very important for a RANS simulation. In the present work, the error

\begin{document}$\epsilon_{\rm{P}}$\end{document}

is studied separately by using the Reynolds stresses obtained from direct numerical simulation (DNS)/highly resolved large-eddy simulation to minimize the modeling error

\begin{document}$\epsilon_{\rm{M}}$\end{document}

, and the sources of

\begin{document}$\epsilon_{\rm{P}}$\end{document}

are derived mathematically. For the implementations with known Reynolds stresses solely, we suggest to run an auxiliary RANS simulation to make a first guess on

\begin{document}$\nu_{\rm{t}}^*$\end{document}

and

\begin{document}$S_{ij}^0$\end{document}

. With around 10 iterations, the error of the streamwise velocity component could be reduced by about one-order of magnitude in flow over periodic hills. The present work is not to develop a new RANS model, but to clarify the facts that obtaining mean field with known Reynolds stresses is nontrivial and that the nonlinear part of the Reynolds stresses is very important in flow problems with separations. The proposed approach to reduce

\begin{document}$\epsilon_{\rm{P}}$\end{document}

may be very useful for the a posteriori applications of the data-driven turbulence models.

Effect of interfacial stiffness on the stretchability of metal/elastomer bilayers under in-plane biaxial tension

Zheng Jia, Teng Li

Accepted Manuscript

[Abstract] (255) [FullText HTML] (38) [PDF 0KB] (1)

Abstract:
Flexible electronic devices are often subjected to large and repeated deformation, so that their functional components such as metal interconnects need to sustain strains up to tens of percent, which is far beyond the intrinsic deformability of metal materials (~1%). To meet the stringent requirements of flexible electronics, metal/elastomer bilayers, a stretchable structure that consists of a metal film adhered to a stretchable elastomer substrate, have been developed to improve the stretch capability of metal interconnects. Previous studies have predicted that the metal/elastomer bilayers are much more stretchable than freestanding metal films. However, these investigations usually assume perfect bonding between the metal and elastomer layers. In this work, the effect of the metal/elastomer interface with a finite interfacial stiffness on the stretchability of bilayer structures is analyzed. The results show that the assumption of perfect interface (with infinite interfacial stiffness) may lead to an overestimation of the stretchability of bilayer structures. It is also demonstrated that increased adhesion between the metal and elastomer layers can enhance the stretchability of the metal layer.

Experimental study of the effect of slotted blades on the Savonius wind turbine performance

Dominicus Danardono Dwi Prija Tjahjana, Zainal Arifin, Suyitno Suyitno, Wibawa Endra Juwana, Aditya Rio Prabowo, Catur Harsito

Accepted Manuscript

[Abstract] (166) [FullText HTML] (32) [PDF 3791KB] (6)

Abstract:
This study investigates the effect of Reynolds number on the performance of Savonius wind turbine with slotted blades. The turbine performance investigation was based on the torque coefficient (Ct), power coefficient (Cp), and tip speed ratio (TSR). The experiment used two number of blade configuration, blade overlap ratio of 10%, 12.5% and 20%, slotted position of 15%, 20%, 25% and 35%, and also slotted gap width of 3 mm, 5 mm, 7 mm, and 9 mm. The wind speed carried out in this experiment are 5.94 m/s, 6.46 m/s, 6.99 m/s, and 7.27 m/s, which are generated from the fan blowers as a wind source. The Savonius turbine with 10% overlap ratio shows the best performance. The highest Cp obtained is 0.138 by the variation of a 3 mm gap with Re of 1.44 × 10⁴ and 0.526 tip speed ratio (TSR).

Fundamental kinematics laws of interstitial fluid flows on vascular walls

Yajun Yin, Hongyi Li, Gang Peng, Xiaobin Yu, Yiya Kong

Accepted Manuscript

[Abstract] (89) [FullText HTML] (39) [PDF 2521KB] (4)

Abstract:
In the previous studies, the phenomenon that the interstitial fluid (ISF) can flow along tunica adventitia of the arteries and veins in both human and animal bodies was reported. On the basis of these studies, this paper aims to: (i) summarize the basic properties of the ISF flows in the walls of arteries and veins, (ii) combine the basic properties with axiomaticism and abstract the axiom for ISF flows, and (iii) propose three fundamental laws of the ISF flow, (i.e., the existence law, the homotropic law and the reverse law). The three laws provide solid theoretical basement for exploring the kinematic patterns of interstitial fluid flow in the cardiovascular system.

Analytical solutions for inflation of pre-stretched elastomeric circular membranes under uniform pressure

Jianghong Yuan, Xianlin Liu, Haibin Xia, Yin Huang

Accepted Manuscript

[Abstract] (175) [FullText HTML] (34) [PDF 2851KB] (11)

Abstract:
Elastomeric membranes are frequently used in several emerging fields such as soft robotics and flexible electronics. For convenience of the structural design, it is very attractive to find simple analytical solutions to well describe their elastic deformations in response to external loadings. However, both the material/geometrical nonlinearity and the deformation inhomogeneity due to boundary constraints make it much challenging to get an exact analytical solution. In this paper, we focus on the inflation of a pre-stretched elastomeric circular membrane under uniform pressure, and derive an approximate analytical solution of the pressure-volume curve based upon a reasonable assumption on the shape of the inflated membrane. Such an explicit expression enables us to quantitatively design the material and geometrical parameters of the pre-stretched membrane to generate a target pressure-volume curve with prescribed peak point and initial slope. This work would be of help in the simplified mechanical design of structures involving elastomeric membranes.

Why neural networks apply to scientific computing?

Shaoqiang Tang, Yang Yang

Accepted Manuscript

[Abstract] (92) [FullText HTML] (37) [PDF 3294KB] (2)

Abstract:
In recent years, neural networks have become an increasingly powerful tool in scientific computing. The universal approximation theorem asserts that a neural network may be constructed to approximate any given continuous function at desired accuracy. The backpropagation algorithm further allows efficient optimization of the parameters in training a neural network. Powered by GPU's, effective computations for scientific and engineering problems are thereby enabled. In addition, we show that finite element shape functions may also be approximated by neural networks.

A Geometry-based framework for modeling the complexity of origami folding

Samuel Schulman, Xin Ning

Accepted Manuscript

[Abstract] (149) [FullText HTML] (34) [PDF 2737KB] (1)

Abstract:
This paper presents a quantitative framework to analyze the complexity of folding origami structures from flat membranes. Extensive efforts have realized intricate origami patterns with desired functions such as mechanical properties, packaging efficiency, and deployment behavior. However, the complexity associated with the manufacturing or folding of origami patterns has not been explored. Understanding how difficult origami structures are to make, and how much time they require to form is crucial information to determining the practical feasibility of origami designs and future applications such as robotic origami assembly in space. In this work, we develop this origami complexity metric by modeling the geometric properties and crease formation of the origami structure, from which it outputs crease and pattern complexity values and a prediction of the time to complete the pattern assembly, based on the characteristics of the operator. The framework is experimentally validated by fabricating various Miura-ori origami paper models.

Complexity analysis of the Portevin-Le Chatelier in an Al alloy at different temperatures

Jianfei Xu, Guoxiong Chen, Shihua Fu

Corrected proof , doi: 10.1016/j.taml.2021.100232

[Abstract] (129) [FullText HTML] (53) [PDF 2611KB] (16)

Abstract:
The complexity of the Portevin-Le Chatelier (PLC) effect in an Al alloy at different temperatures was analyzed by modified multiscale entropy. The results show that three evolutions of entropy with scale factor, i.e., near zero, monotonically increasing and peak-shape, were observed corresponding to the smooth curves, type-A serrations and type-B/-C serrations, respectively. The scale factor at the peak was one-third of the average serration period. The sample entropy increased initially and then decreased with temperature, which was opposite to the critical strain. It is also suggested that the type-A serrations corresponded to self-organized criticality and the type-B/-C serrations corresponded to chaos through the evolutions of entropy with scale factor.

Analytical solutions for sediment concentration in waves based on linear diffusivity

Yiqin Xie, Jifu Zhou, Xu Wang, Yanrong Kuai, Yongjun Lu

Corrected proof

[Abstract] (198) [FullText HTML] (53) [PDF 2742KB] (2)

Abstract:
Two kinds of analytical solutions are derived through Laplace transform for the equation that governs wave-induced suspended sediment concentration with linear sediment diffusivity under two kinds of bottom boundary conditions, namely the reference concentration (Dirichlet) and pickup function (Numann), based on a variable transformation that is worked out to transform the governing equation into a modified Bessel equation. The ability of the two analytical solutions to describe the profiles of suspended sediment concentration is discussed by comparing with different experimental data. And it is demonstrated that the two analytical solutions can well describe the process of wave-induced suspended sediment concentration, including the amplitude and phase and vertical profile of sediment concentration. Furthermore, the solution with boundary condition of pickup function provides better results than that of reference concentration in terms of the phase-dependent variation of concentration.

A model for universal spatial variations of temperature fluctuations in turbulent Rayleigh-Bénard convection

Xiaozhou He, Eberhard Bodenschatz, Guenter Ahlers

Corrected proof

[Abstract] (196) [FullText HTML] (53) [PDF 2444KB] (7)

Abstract:
We propose a theoretical model for spatial variations of the temperature variance

\begin{document}$\sigma^2(z,r)$\end{document}

(z is the distance from the sample bottom and r the radial coordinate) in turbulent Rayleigh-Bénard convection (RBC). Adapting the “attached-eddy” model of shear flow to the plumes of RBC, we derived an equation for

\begin{document}$\sigma^2$\end{document}

which is based on the universal scaling of the normalized RBC temperature spectra. This equation includes both logarithmic and power-law dependences on

\begin{document}$z/\lambda_{th}$\end{document}

, where

\begin{document}$\lambda_{th}$\end{document}

is the thermal boundary layer thickness. The equation parameters depend on r and the Prandtl number Pr, but have only an extremely weak dependence on the Rayleigh number Ra Thus our model provides a near-universal equation for the temperature variance profile in turbulent RBC.

Laser-induced porous graphene on Polyimide/PDMS composites and its kirigami-inspired strain sensor

Hao Wang, Zifeng Zhao, Panpan Liu, Xiaogang Guo

Accepted Manuscript

[Abstract] (189) [FullText HTML] (42) [PDF 2897KB] (3)

Abstract:
The laser-induced porous graphene (LIG) prepared in a straightforward fabrication method is presented, and its applications in stretchable strain sensors to detect the applied strain are also explored. The LIG formed on the polyimide/polydimethylsiloxane (PI/PDMS) composite exhibits a naturally high stretchability (over 30%), bypassing the transfer printing process compared to the one prepared by laser scribing on PI films. The PI/PDMS composite with LIG shows tunable mechanical and electronic performances with different PI particle concentrations in PDMS. The good cyclic stability and almost linear response of the prepared LIG's resistance with respect to tensile strain provide its access to wearable electronics. To improve the PDMS/PI composite stretchability, we designed and optimized a kirigami-inspired strain sensor with LIG on the top surface, dramatically increasing the maximum strain value that in linear response to applied strain from 3% to 79%.

Bifurcation mechanism of interfacial electrohydrodynamic gravity-capillary waves near the minimum phase speed under a horizontal electric field

Gexing Xu, Zhi Lin

Corrected proof , doi: 10.1016/j.taml.2021.100224

[Abstract] (127) [FullText HTML] (69) [PDF 2928KB] (5)

Abstract:
We investigate the evolution of interfacial gravity-capillary waves propagating along the interface between two dielectric fluids under the action of a horizontal electric field. There is a uniform background flow in each layer, and the relative motion tends to induce Kelvin-Helmholtz (KH) instability. The combined effects of gravity, surface tension and electrically induced forces are all taken into account. Under the short-wave assumption, the expansion and truncation method of Dirichlet-Neumann (DN) operators is applied to derive a reduced dynamical model. When KH instability is suppressed linearly by a considerably large electric field, our numerical results reveal that in certain regions of parameter space, nonlinear symmetric traveling wave solutions can be found near the minimum phase speed. Additionally, the detailed bifurcation structures are presented together with typical wave profiles.

Meandering motions within the viscous sublayer: Supplementary materials

S. Santosh Kumar, Xinyi Huang, Xiang Yang, Jiarong Hong

Corrected proof

[Abstract] (190) [FullText HTML] (54) [PDF 2904KB] (3)

Abstract:
$v.abstractInfoEn

Design of crashworthy attenuator structures as a part of vehicle safety against impact: Application of waste aluminum can-based material

Laksmana Widi Prasetya, Aditya Rio Prabowo, Ubaidillah Ubaidillah, Iwan Istanto, Nur Azmah Binti Nordin

Accepted Manuscript

[Abstract] (196) [FullText HTML] (45) [PDF 3535KB] (6)

Abstract:
The impact attenuator is an essential system in both race cars and urban vehicles. The structure of an impact attenuator serves as a safety barrier between the impacted surface and the driver in an accident. Attenuator materials tend to have a high price; thus, alternative materials were explored in the current work, i.e., used cans from food and beverage containers. The study deployed a nonlinear finite element algorithm to calculate a series of impacts on the attenuator structures. The thickness of the cans and velocity of the impact were considered as the main parameters. Analysis results concluded that the attenuator's average energy was 16,000 J for a can thickness of 1 mm. This value is more than two times the 0.5 mm thick used cans. The attenuator's new design was then matched with an attenuator regulation, and the results surpassed the standard value of 7350 J.

Meandering motions within the viscous sublayer

S. Santosh Kumar, Xinyi Huang, Xiang Yang, Jiarong Hong

Accepted Manuscript

[Abstract] (274) [FullText HTML] (54) [PDF 4109KB] (1)

Abstract:
Impact of viscous sublayer scale roughness elements on large scale flows have not been fully understood and require high resolution 3D flow measurements to unravel. However, existing approaches fail to provide sufficient resolution for such measurements to fully resolve the sublayer. In this study, we use digital Fresnel reflection holography to capture 3D flows within the viscous sublayer at sub-viscous resolutions. The measurement highlights the presence of novel flow structures at the scale of the sublayer, with strong spanwise meandering motions, of 2-3 viscous wall units, indicating a highly unsteady and accelerating flow within. The probability distribution of accelerations shows a stretched exponential shapes characteristic of highly intermittent turbulence seen under isotropic flows. The presence of flow structures even at the scale of the sublayer, i.e., below y⁺~5, points to the effectiveness of roughness elements in modulating the large scale flow.

Research on refined reconstruction method of airfoil pressure based on compressed sensing

Xuan Zhao, Lin Du, Xuhao Peng, Zichen Deng, Weiwei Zhang

Corrected proof , doi: 10.1016/j.taml.2021.A004

[Abstract] (202) [FullText HTML] (69) [PDF 3132KB] (8)

Abstract:
The placement of pressure taps on the surface of the wind tunnel test model is an important means to obtain the surface pressure distribution. However, limited by space location and experimental cost, it is difficult to arrange enough pressure measuring taps on the surface of complex models to obtain complete pressure distribution information, thus it is impossible to obtain accurate lift and moment characteristics through integration. The paper proposes a refined reconstruction method of airfoil surface pressure based on compressed sensing, which can reconstruct the pressure distribution with high precision with less pressure measurement data. Tests on typical airfoil subsonic flow around flow show that the accuracy of lift and moment after the pressure integration reconstructed by 4-8 measuring points can meet the requirements of the national military standard. The algorithm is robust to noise, and provides a new idea for obtaining accurate force data from sparse surface pressure tests in engineering.

Select All

Display Method: |

Special Issue

Multiscale mechanics

G.W. He, G.D. Jin

2021, 11(1): 1 -1. doi: 10.1016/j.taml.2021.100238

[Abstract] (121) [FullText HTML] (93) [PDF 2196KB] (33)

Abstract:
$v.abstractInfoEn

Article

A study of inner-outer interactions in turbulent channel flows by interactive POD

Hongping Wang, Qi Gao

2021, 11(1): 2 -15. doi: 10.1016/j.taml.2021.100222

[Abstract] (193) [FullText HTML] (92) [PDF 3932KB] (24)

Abstract:
The amplitude and frequency modulation of near-wall flow structures by the large-scale motions in outer regions is studied in turbulent channel flows. The proper orthogonal decomposition (POD) method is applied to investigate the interactions between the near-wall motions and the large-scale flow modes of the outer regions based on two datasets from direct numerical simulation of turbulent channel flows at Reynolds numbers of 550 and 1000. The fluctuations in the fields

\begin{document}$u^+$\end{document}

\begin{document}$v^+$\end{document}

\begin{document}$w^+$\end{document}

and Reynolds shear stress

\begin{document}$-(uv)^+$\end{document}

are studied to understand the mechanism of amplitude and frequency modulation of the near-wall structures by the outer large-scale motions. The amplitude modulation coefficient of the Reynolds shear stress is larger than that of the velocity components. The frequency modulation effect has an opposite influence in the spanwise direction compared to the streamwise direction. The streamwise characteristic frequency increases with increasing large-scale velocity. However, the spanwise characteristic frequency exhibits a decreasing trend with increasing large-scale velocity in the near-wall region.

Letter

A study on equivalence of nonlinear energy dissipation between first-order computational homogenization (FOCH) and reduced-order homogenization (ROH) methods

Jiajia Yue, Zifeng Yuan

2021, 11(1): 16 -21. doi: 10.1016/j.taml.2021.100225

[Abstract] (115) [FullText HTML] (82) [PDF 2931KB] (5)

Abstract:
Nowadays, studies on the mechanism of macro-scopic nonlinear behavior of materials by accumulation of micro-scopic degradation are attracting more attention from researchers. Among numerous approaches, multiscale methods have been proved as powerful and practical approaches in predicting macro-scopic material status by averaging and homogenizing physical information from associated micro-scopic material behavior. Usually in mechanical problem, the stress, consistent material modulus, and possible material state variables are quantities in interest through the upscaling process. However, the energy-related quantities are not studied much. Some initiative work has been done in the early year including but not limited to the Hill–Mandel condition in multiscale framework, which gives that the macro-scopic elastic strain energy density can be computed by volumetric averaging of that in the micro-scale. However, in the nonlinear analysis, the energy dissipation is an important quantity to measure the degradation status. In this manuscript, two typical multiscale methods, the first-order computational homogenization (FOCH) and reduced-order homogenization (ROH), are adopted to numerically analyze a fiber-reinforced composite material with capability in material nonlinearity. With numerical experiments, it can be shown that energy dissipation is the same for both approaches.

Stable heat jet approach for temperature control of Fermi-Pasta-Ulam beta chain

Baiyili Liu, Qian Zhang, Shaoqiang Tang

2021, 11(1): 22 -27. doi: 10.1016/j.taml.2021.100226

[Abstract] (103) [FullText HTML] (85) [PDF 2869KB] (5)

Abstract:
In this paper, we propose a stable heat jet approach for accurate temperature control of the nonlinear Fermi-Pasta-Ulam beta (FPU-

\begin{document}$ \bf{\beta}$\end{document}

) chain. First, we design a stable nonlinear boundary condition, with coefficients determined by a machine learning technique. Its stability can be proved rigorously. Based on this stable boundary condition, we derive a two-way boundary condition complying with phonon heat source, and construct stable heat jet approach. Numerical tests illustrate the stability of the boundary condition and the effectiveness in eliminating boundary reflections. Furthermore, we extend the boundary condition formulation with more atoms, and train the coefficients to eliminate extreme short waves by machine learning technique. Under this extended boundary condition, the heat jet approach is effective for high temperature, and may be adopted for multiscale computation of atomic motion at finite temperature.

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

Linfeng Jiang, Cheng Wang, Shuang Liu, Chao Sun, Enrico Calzavarini

2021, 11(1): 28 -34. doi: 10.1016/j.taml.2021.100227

[Abstract] (154) [FullText HTML] (125) [PDF 3015KB] (15)

Abstract:
We successfully perform the three-dimensional tracking in a turbulent fluid flow of small axisymmetrical particles that are neutrally-buoyant and bottom-heavy, i.e., they have a non-homogeneous mass distribution along their symmetry axis. We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate. The experiment is complemented by a series of simulations based on realistic Navier–Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravitational torque stability number and by the particle aspect ratio. We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential orientation and tumbling rate of the particles. We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.

Wall-resolved large-eddy simulation of turbulent channel flows with rough walls

Shilong Li, Xiaolei Yang, Guodong Jin, Guowei He

2021, 11(1): 35 -41. doi: 10.1016/j.taml.2021.100228

[Abstract] (154) [FullText HTML] (103) [PDF 16306KB] (19)

Abstract:
Turbulent flows over rough surfaces widely exist in nature and industry. Investigating its mechanism is of theoretical and practical significance. In this work we simulate the turbulent channel flow with rough walls using large-eddy simulation with rough elements resolved using the curvilinear immersed boundary method and compare the results obtained in this work with those in the paper by Yuan and Piomelli (J. Fluid Mech. 2014), where the volume of fluid method was employed for modeling rough elements. The mean streamwise velocity profiles predicted by the two methods agree well with each other. Differences in Reynolds stresses and dispersive stresses are observed, which are attributed to the different approaches in dealing with the complex geometry of the rough surface.

Constrained large-eddy simulation of turbulent flow over inhomogeneous rough surfaces

Wen Zhang, Minping Wan, Zhenhua Xia, Jianchun Wang, Xiyun Lu, Shiyi Chen

2021, 11(1): 42 -47. doi: 10.1016/j.taml.2021.100229

[Abstract] (124) [FullText HTML] (95) [PDF 3112KB] (12)

Abstract:
In this work we extend the method of the constrained large-eddy simulation (CLES) to simulate the turbulent flow over inhomogeneous rough walls. In the original concept of CLES, the subgrid-scale (SGS) stress is constrained so that the mean part and the fluctuation part of the SGS stress can be modelled separately to improve the accuracy of the simulation result. Here in the simulation of the rough-wall flows, we propose to interpret the extra stress terms in the CLES formulation as the roughness-induced stress so that the roughness inhomogeneity can be incorporated by modifying the formulation of the constrained SGS stress. This is examined with the simulations of the channel flow with the spanwise alternating high/low roughness strips. Then the CLES method is employed to investigate the temporal response of the turbulence to the change of the wall condition from rough to smooth. We demonstrate that the temporal development of the internal boundary layer is just similar to that in a spatial rough-to-smooth transition process, and the spanwise roughness inhomogeneity has little impact on the transition process.

A multi-scale algorithm for dislocation creep at elevated temperatures

Lichao Yuan, Yujie Wei

2021, 11(1): 48 -53. doi: 10.1016/j.taml.2021.100230

[Abstract] (122) [FullText HTML] (72) [PDF 2840KB] (10)

Abstract:
Dislocation creep at elevated temperatures plays an important role for plastic deformation in crystalline metals. When using traditional discrete dislocation dynamics (DDD) to capture this process, we often need to update the forces on N dislocations involving ~N² interactions. In this letter, we introduce a multi-scale algorithm to speed up the calculations by dividing a sample of interest into sub-domain grids: dislocations within a characteristic area interact following the conventional way, but their interaction with dislocations in other grids are simplified by lumping all dislocations in another grid as a super one. Such a multi-scale algorithm lowers the computational load to ~N^1.5. We employed this algorithm to model dislocation creep in Al-Mg alloy. The simulation leads to a power-law creep rate in consistent with experimental observations. The stress exponent of the power-law creep is a resultant of dislocations climb for ~5 and viscous dislocations glide for ~3.

Numerical investigation on motion of an ellipsoidal particle inside confined microcavity flow

Tianrou Miao, Zuoli Xiao

2021, 11(1): 54 -59. doi: 10.1016/j.taml.2021.100234

[Abstract] (89) [FullText HTML] (66) [PDF 3045KB] (5)

Abstract:
The behaviors of a neutrally buoyant ellipsoidal particle in vortical flow confined by a microcavity are numerically studied using the Lattice-Boltzmann method. For specific initial position, an isolated ellipsoid may develop a stable limit cycle orbit inside microcavity due to the interaction between particle and the carrier flow. It is observed that ellipsoidal particles of different shapes exhibit two different stable rotational modes depending on the initial orientation and lateral position. A prolate spheroid tends to enter a tumbling mode whereas an oblate spheroid is apt to achieve a rolling mode. The evolution of rotational velocities along the stable orbit is also analyzed for particles of defferent shapes.

Physics-informed deep learning for digital materials

Zhizhou Zhang, Grace X. Gu

2021, 11(1): 60 -66. doi: 10.1016/j.taml.2021.100220

[Abstract] (284) [FullText HTML] (111) [PDF 2873KB] (15)

Abstract:
In this work, a physics-informed neural network (PINN) designed specifically for analyzing digital materials is introduced. This proposed machine learning (ML) model can be trained free of ground truth data by adopting the minimum energy criteria as its loss function. Results show that our energy-based PINN reaches similar accuracy as supervised ML models. Adding a hinge loss on the Jacobian can constrain the model to avoid erroneous deformation gradient caused by the nonlinear logarithmic strain. Lastly, we discuss how the strain energy of each material element at each numerical integration point can be calculated parallelly on a GPU. The algorithm is tested on different mesh densities to evaluate its computational efficiency which scales linearly with respect to the number of nodes in the system. This work provides a foundation for encoding physical behaviors of digital materials directly into neural networks, enabling label-free learning for the design of next-generation composites.

Experimental investigation of sound absorption in a composite absorber

Nansha Gao, Hong Hou

2021, 11(1): 67 -72. doi: 10.1016/j.taml.2021.100221

[Abstract] (129) [FullText HTML] (72) [PDF 2785KB] (4)

Abstract:
A composite absorber made of a polyurethane sponge and multi-layer micro-perforated plates is presented in this study. Results from an acoustic impedance tube test show that the polyurethane sponge can exhibits higher low-frequency sound absorption in front of the micro-perforated plate, while sound absorption at medium and high-frequencies remains low. The physical mechanism behind this is that the micro-perforated plate increases the denpth cavity. If the polyurethane sponge is placed behind the micro-perforated plate, the amplitude of the original absorption peak will remain constant, but the absorption peaks will shift to lower frequencies. The reason for this phenomenon is that porous materials with low flow resistance can be approximately equivalent to fluid, which not only does not affect the resonance absorption coefficient of micro-perforated plate, but also makes the peaks move to low frequency. This study has the potential applications in the sound absorption design of composite structure.