More Like this

1. Frictional Damping from Biomimetic Scales

   https://doi.org/10.1038/s41598-019-50944-0  

   Ali, Hessein ; Ebrahimi, Hossein ; Ghosh, Ranajay ( December 2019 , Scientific Reports)

   null (Ed.)

   Abstract Stiff scales adorn the exterior surfaces of fishes, snakes, and many reptiles. They provide protection from external piercing attacks and control over global deformation behavior to aid locomotion, slithering, and swimming across a wide range of environmental condition. In this report, we investigate the dynamic behavior of biomimetic scale substrates for further understanding the origins of the nonlinearity that involve various aspect of scales interaction, sliding kinematics, interfacial friction, and their combination. Particularly, we study the vibrational characteristics through an analytical model and numerical investigations for the case of a simply supported scale covered beam. Our results reveal for the first time that biomimetic scale beams exhibit viscous damping behavior even when only Coulomb friction is postulated for free vibrations. We anticipate and quantify the anisotropy in the damping behavior with respect to curvature. We also find that unlike static pure bending where friction increases bending stiffness, a corresponding increase in natural frequency for the dynamic case does not arise for simply supported beam. Since both scale geometry, distribution and interfacial properties can be easily tailored, our study indicates a biomimetic strategy to design exceptional synthetic materials with tailorable damping behavior. 

   more » « less
2. Inverse design of mesoscopic models for compressible flow using the Chapman-Enskog analysis

   https://doi.org/10.1186/s42774-020-00059-2  

   Chen, Tao ; Wang, Lian-Ping ; Lai, Jun ; Chen, Shiyi ( February 2021 , Advances in Aerodynamics)

   In this paper, based on simplified Boltzmann equation, we explore the inverse-design of mesoscopic models for compressible flow using the Chapman-Enskog analysis. Starting from the single-relaxation-time Boltzmann equation with an additional source term, two model Boltzmann equations for two reduced distribution functions are obtained, each then also having an additional undetermined source term. Under this general framework and using Navier-Stokes-Fourier (NSF) equations as constraints, the structures of the distribution functions are obtained by the leading-order Chapman-Enskog analysis. Next, five basic constraints for the design of the two source terms are obtained in order to recover the NSF system in the continuum limit. These constraints allow for adjustable bulk-to-shear viscosity ratio, Prandtl number as well as a thermal energy source. The specific forms of the two source terms can be determined through proper physical considerations and numerical implementation requirements. By employing the truncated Hermite expansion, one design for the two source terms is proposed. Moreover, three well-known mesoscopic models in the literature are shown to be compatible with these five constraints. In addition, the consistent implementation of boundary conditions is also explored by using the Chapman-Enskog expansion at the NSF order. Finally, based on the higher-order Chapman-Enskog expansion of the distribution functions, we derive the complete analytical expressions for the viscous stress tensor and the heat flux. Some underlying physics can be further explored using the DNS simulation data based on the proposed model.

    

   more » « less
3. The Scaling of Vortical Electron Acceleration in Thin-current Magnetic Reconnection and Its Implications in Solar Flares

   https://doi.org/10.3847/1538-4357/ad09e6  

   Crawford, C. ; Che, H. ; Benz, A. O. ( January 2024 , The Astrophysical Journal)

   To investigate how magnetic reconnection (MR) accelerates electrons to a power-law energy spectrum in solar flares, we explore the scaling of a kinetic model proposed by Che & Zank (CZ) and compare it to observations. Focusing on thin current sheet MR particle-in-cell (PIC) simulations, we analyze the impact of domain size on the evolution of the electron Kelvin–Helmholtz instability (EKHI). We find that the duration of the growth stage of the EKHI ($t_G\sim {\mathrm{\Omega }}_e^{-1}$) is short and remains nearly unchanged because the electron gyrofrequency Ω_eis independent of domain size. The quasi-steady stage of the EKHI (t_MR) dominates the electron acceleration process and scales linearly with the size of the simulations asL/v_A0, wherev_A0is the Alfvén speed. We use the analytical results obtained by CZ to calculate the continuous temporal evolution of the electron energy spectra from PIC simulations and linearly scale them to solar flare observational scales. For the first time, an electron acceleration model predicts the sharp two-stage transition observed in typical soft–hard–harder electron energy spectra, implying that the electron acceleration model must be efficient with an acceleration timescale that is a small fraction of the duration of solar flares. Our results suggest that we can use PIC MR simulations to investigate the observational electron energy spectral evolution of solar flares if the ratiot_MR/t_Gis sufficiently small, i.e., ≲10%.

    

   more » « less
4. Generalized Fluid Models of the Braginskii Type

   https://doi.org/10.3847/1538-4365/ac5044  

   Hunana, P. ; Passot, T. ; Khomenko, E. ; Martínez-Gómez, D. ; Collados, M. ; Tenerani, A. ; Zank, G. P. ; Maneva, Y. ; Goldstein, M. L. ; Webb, G. M. ( June 2022 , The Astrophysical Journal Supplement Series)

   Abstract Several generalizations of the well-known fluid model of Braginskii (1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multispecies plasmas with arbitrary masses and temperatures, where all of the fluid moments are described by their evolution equations. The 21-moment model contains two “heat flux vectors” (third- and fifth-order moments) and two “viscosity tensors” (second- and fourth-order moments). The Braginskii model is then obtained as a particular case of a one ion–electron plasma with similar temperatures, with decoupled heat fluxes and viscosity tensors expressed in a quasistatic approximation. We provide all of the numerical values of the Braginskii model in a fully analytic form (together with the fourth- and fifth-order moments). For multispecies plasmas, the model makes the calculation of the transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasistatic approximation that makes some Braginskii coefficients divergent in a weakly collisional regime. Importantly, we show that the heat fluxes and viscosity tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the fourth-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad. 

   more » « less
5. Stochastic theory and direct numerical simulations of the relative motion of high-inertia particle pairs in isotropic turbulence

   https://doi.org/10.1017/jfm.2016.859  

   Dhariwal, Rohit ; Rani, Sarma L. ; Koch, Donald L. ( February 2017 , Journal of Fluid Mechanics)

   The relative velocities and positions of monodisperse high-inertia particle pairs in isotropic turbulence are studied using direct numerical simulations (DNS), as well as Langevin simulations (LS) based on a probability density function (PDF) kinetic model for pair relative motion. In a prior study (Rani et al. , J. Fluid Mech. , vol. 756, 2014, pp. 870–902), the authors developed a stochastic theory that involved deriving closures in the limit of high Stokes number for the diffusivity tensor in the PDF equation for monodisperse particle pairs. The diffusivity contained the time integral of the Eulerian two-time correlation of fluid relative velocities seen by pairs that are nearly stationary. The two-time correlation was analytically resolved through the approximation that the temporal change in the fluid relative velocities seen by a pair occurs principally due to the advection of smaller eddies past the pair by large-scale eddies. Accordingly, two diffusivity expressions were obtained based on whether the pair centre of mass remained fixed during flow time scales, or moved in response to integral-scale eddies. In the current study, a quantitative analysis of the (Rani et al. 2014) stochastic theory is performed through a comparison of the pair statistics obtained using LS with those from DNS. LS consist of evolving the Langevin equations for pair separation and relative velocity, which is statistically equivalent to solving the classical Fokker–Planck form of the pair PDF equation. Langevin simulations of particle-pair dispersion were performed using three closure forms of the diffusivity – i.e. the one containing the time integral of the Eulerian two-time correlation of the seen fluid relative velocities and the two analytical diffusivity expressions. In the first closure form, the two-time correlation was computed using DNS of forced isotropic turbulence laden with stationary particles. The two analytical closure forms have the advantage that they can be evaluated using a model for the turbulence energy spectrum that closely matched the DNS spectrum. The three diffusivities are analysed to quantify the effects of the approximations made in deriving them. Pair relative-motion statistics obtained from the three sets of Langevin simulations are compared with the results from the DNS of (moving) particle-laden forced isotropic turbulence for $St_{\unicode[STIX]{x1D702}}=10,20,40,80$ and $Re_{\unicode[STIX]{x1D706}}=76,131$ . Here, $St_{\unicode[STIX]{x1D702}}$ is the particle Stokes number based on the Kolmogorov time scale and $Re_{\unicode[STIX]{x1D706}}$  is the Taylor micro-scale Reynolds number. Statistics such as the radial distribution function (RDF), the variance and kurtosis of particle-pair relative velocities and the particle collision kernel were computed using both Langevin and DNS runs, and compared. The RDFs from the stochastic runs were in good agreement with those from the DNS. Also computed were the PDFs $\unicode[STIX]{x1D6FA}(U|r)$ and $\unicode[STIX]{x1D6FA}(U_{r}|r)$ of relative velocity $U$ and of the radial component of relative velocity $U_{r}$ respectively, both PDFs conditioned on separation $r$ . The first closure form, involving the Eulerian two-time correlation of fluid relative velocities, showed the best agreement with the DNS results for the PDFs. 

   more » « less