Minimum-dissipation model for large-eddy simulation using symmetry-preserving discretization in OpenFOAM

The minimum-dissipation model is applied to channel flow up to $Re_\tau = 2000$, flow past a circular cylinder at $Re=3900$, and flow over periodic hills at $Re=10595$. Numerical simulations were performed in OpenFOAM which is based on the finite volume methods. We used both symmetry-preserving and standard second-order accurate discretization methods in OpenFOAM on structured meshes. The results are compared to DNS and experimental data. The results of channel flow demonstrate a static QR model performs equally well as the dynamic models while reducing the computational cost. The model constant of $C=0.024$ gives the most accurate prediction, and the contribution of the sub-grid model decreases with the increase of the mesh resolution and becomes very small (less than 0.2 molecular viscosity) if a fine mesh is used. Furthermore, the QR model is able to predict the mean and rms velocity accurately up to $Re_\tau = 2000$ without a wall damping function. The symmetry-preserving discretization outperforms the standard OpenFOAM discretization at $Re_\tau=1000$. The results for the flow over a cylinder show that the mean velocity, drag coefficient, and lift coefficient are in good agreement with the experimental data and the central difference schemes conjugated with the QR model predict better results. The various comparisons carried out for flows over periodic hills demonstrate the need to use central difference schemes in OpenFOAM in combination with the minimum dissipation model. The best model constant is again $C=0.024$. The single wind turbine simulation shows that the QR model is capable of predicting accurate results in complex rotating scenarios.Comment: 6 pages; 12 figures; ETMM14 conference paper. arXiv admin note: substantial text overlap with arXiv:2309.0441

Too hot to handle: heat stroke in multiple system atrophy

Item does not contain fulltex

Magnetic-field-induced propulsion of jellyfish-inspired soft robotic swimmers

The multifaceted appearance of soft robots in the form of swimmers, catheters, surgical devices, and drug-carrier vehicles in biomedical and microfluidic applications is ubiquitous today. Jellyfish-inspired soft robotic swimmers (jellyfishbots) have been fabricated and experimentally characterized by several researchers that reported their swimming kinematics and multimodal locomotion. However, the underlying physical mechanisms that govern magnetic-field-induced propulsion are not yet fully understood. Here, we use a robust and efficient computational framework to study the jellyfishbot swimming kinematics and the induced flow field dynamics through numerical simulation. We consider a two-dimensional model jellyfishbot that has flexible lappets, which are symmetric about the jellyfishbot center. These lappets exhibit flexural deformation when subjected to external magnetic fields to displace the surrounding fluid, thereby generating the thrust required for propulsion. We perform a parametric sweep to explore the jellyfishbot kinematic performance for different system parameters—structural, fluidic, and magnetic. In jellyfishbots, the soft magnetic composite elastomeric lappets exhibit temporal and spatial asymmetries when subjected to unsteady external magnetic fields. The average speed is observed to be dependent on both these asymmetries, quantified by the glide magnitude and the net area swept by the lappet tips per swimming cycle, respectively. We observe that a judicious choice of the applied magnetic field and remnant magnetization profile in the jellyfishbot lappets enhances both these asymmetries. Furthermore, the dependence of the jellyfishbot swimming speed upon the net area swept (spatial asymmetry) is twice as high as the dependence of speed on the glide ratio (temporal asymmetry). Finally, functional relationships between the swimming speed and different kinematic parameters and nondimensional numbers are developed. Our results provide guidelines for the design of improved jellyfish-inspired magnetic soft robotic swimmers

Low-Dissipation Simulation Methods and Models for Turbulent Subsonic Flow

The simulation of turbulent flows by means of computational fluid dynamics is highly challenging. The costs of an accurate direct numerical simulation (DNS) are usually too high, and engineers typically resort to cheaper coarse-grained models of the flow, such as large-eddy simulation (LES). To be suitable for the computation of turbulence, methods should not numerically dissipate the turbulent flow structures. Therefore, energy-conserving discretizations are investigated, which do not dissipate energy and are inherently stable because the discrete convective terms cannot spuriously generate kinetic energy. They have been known for incompressible flow, but the development of such methods for compressible flow is more recent. This paper will focus on the latter: LES and DNS for turbulent subsonic flow. A new theoretical framework for the analysis of energy conservation in compressible flow is proposed, in a mathematical notation of square-root variables, inner products, and differential operator symmetries. As a result, the discrete equations exactly conserve not only the primary variables (mass, momentum and energy), but also the convective terms preserve (secondary) discrete kinetic and internal energy. Numerical experiments confirm that simulations are stable without the addition of artificial dissipation. Next, minimum-dissipation eddy-viscosity models are reviewed, which try to minimize the dissipation needed for preventing sub-grid scales from polluting the numerical solution. A new model suitable for anisotropic grids is proposed: the anisotropic minimum-dissipation model. This model appropriately switches off for laminar and transitional flow, and is consistent with the exact sub-filter tensor on anisotropic grids. The methods and models are first assessed on several academic test cases: channel flow, homogeneous decaying turbulence and the temporal mixing layer. As a practical application, accurate simulations of the transitional flow over a delta wing have been performed

When Does Eddy Viscosity Damp Subfilter Scales Sufficiently?

Large eddy simulation (LES) seeks to predict the dynamics of spatially filtered turbulent flows. The very essence is that the LES-solution contains only scales of size ≥Δ, where Δ denotes some user-chosen length scale. This property enables us to perform a LES when it is not feasible to compute the full, turbulent solution of the Navier-Stokes equations. Therefore, in case the large eddy simulation is based on an eddy viscosity model we determine the eddy viscosity such that any scales of size <Δ are dynamically insignificant. In this paper, we address the following two questions: how much eddy diffusion is needed to (a) balance the production of scales of size smaller than Δ; and (b) damp any disturbances having a scale of size smaller than Δ initially. From this we deduce that the eddy viscosity νe has to depend on the invariants q = ½tr(S^2) and r =−⅓tr(S^3) of the (filtered) strain rate tensor S. The simplest model is then given by νe = 3/2(Δ/π)^2|r|/q. This model is successfully tested for a turbulent channel flow (Reτ = 590).

Mixed modeling for large-eddy simulation: The single-layer and two-layer minimum-dissipation-Bardina models

Predicting the behavior of turbulent flows using large-eddy simulation requires modeling of the subgrid-scale stress tensor. This tensor can be approximated using mixed models, which combine the dissipative nature of functional models with the capability of structural models to approximate out-of-equilibrium effects. We propose a mathematical basis to mix (functional) eddy-viscosity models with the (structural) Bardina model. By taking an anisotropic minimum-dissipation (AMD) model for the eddy viscosity, we obtain the (single-layer) AMD-Bardina model. In order to also obtain a physics-conforming model for wall-bounded flows, we further develop this mixed model into a two-layer approach: the near-wall region is parameterized with the AMD-Bardina model, whereas the outer region is computed with the Bardina model. The single-layer and two-layer AMD-Bardina models are tested in turbulent channel flows at various Reynolds numbers, and improved predictions are obtained when the mixed models are applied in comparison to the computations with the AMD and Bardina models alone. The results obtained with the two-layer AMD-Bardina model are particularly remarkable: both first- and second-order statistics are extremely well predicted and even the inflection of the mean velocity in the channel center is captured. Hence, a very promising model is obtained for large-eddy simulations of wall-bounded turbulent flows at moderate and high Reynolds numbers.Comment: 29 pages, 14 figures, 3 tables; revised, accepted manuscrip

Intrinsic Detectivity Limits of Organic Near-Infrared Photodetectors

Organic photodetectors (OPDs) with a performance comparable to that of conventional inorganic ones have recently been demonstrated for the visible regime. However, near-infrared photodetection has proven to be challenging and, to date, the true potential of organic semiconductors in this spectral range (800–2500&nbsp;nm) remains largely unexplored. In this work, it is shown that the main factor limiting the specific detectivity (D*) is non-radiative recombination, which is also known to be the main contributor to open-circuit voltage losses. The relation between open-circuit voltage, dark current, and noise current is demonstrated using four bulk-heterojunction devices based on narrow-gap donor polymers. Their maximum achievable D* is calculated alongside a large set of devices to demonstrate an intrinsic upper limit of D* as a function of the optical gap. It is concluded that OPDs have the potential to be a useful technology up to 2000&nbsp;nm, given that high external quantum efficiencies can be maintained at these low photon energies

A symmetry-preserving second-order time-accurate PISO-based method

A new conservative symmetry-preserving second-order time-accurate PISO-based pressure-velocity coupling for solving the incompressible Navier-Stokes equations on unstructured collocated grids is presented in this paper. This new method for implicit time stepping is an extension of the conservative symmetry-preserving incremental-pressure projection method for explicit time stepping and unstructured collocated meshes of Trias et al. [35]. In order to assess and compare both methods, we have implemented them within one unified solver in the open source code OpenFOAM where we use a Butcher array to prescribe the Runge-Kutta method. Thus, by changing the entries of the Butcher array, explicit and diagonally implicit Runge-Kutta schemes can be combined into one solver. We assess the energy conservation properties of the implemented discretisation methods and the temporal consistency of the selected Runge-Kutta schemes using Taylor-Green vortex and lid-driven cavity flow test cases. Finally, we use a more complex turbulent channel flow test case in order to further assess the performance of the presented new conservative symmetry-preserving incremental-pressure PISO-based method. Although both implemented methods are based on a symmetry-preserving discretisation, we show they still produce a small amount of numerical dissipation when the total pressure is directly solved from a Poisson equation. When an incremental-pressure approach is used, where a pressure correction is solved from a Poisson equation, both methods are effectively fully-conservative. For high-fidelity simulations of incompressible turbulent flows, it is highly desirable to use fully-conservative methods. For such simulations, the presented numerical methods are therefore expected to have large added value, since they pave the way for the execution of truly energy-conservative high-fidelity simulations in complex geometries. Furthermore, both methods are implemented in OpenFOAM, which is widely used within the CFD community, so that a large part of this community can benefit from the developed and implemented numerical methods

The Herschel Virgo Cluster Survey: IV. Resolved dust analysis of spiral galaxies

We present a resolved dust analysis of three of the largest angular size spiral galaxies, NGC 4501 and NGC 4567/8, in the Herschel Virgo Cluster Survey (HeViCS) Science Demonstration field. Herschel has unprecedented spatial resolution at far-infrared wavelengths and with the PACS and SPIRE instruments samples both sides of the peak in the far infrared spectral energy distribution (SED).We present maps of dust temperature, dust mass, and gas-to-dust ratio, produced by fitting modified black bodies to the SED for each pixel. We find that the distribution of dust temperature in both systems is in the range ~19 - 22 K and peaks away from the centres of the galaxies. The distribution of dust mass in both systems is symmetrical and exhibits a single peak coincident with the galaxy centres. This Letter provides a first insight into the future analysis possible with a large sample of resolved galaxies to be observed by Herschel.Comment: Letter accepted for publication in A&A (Herschel special issue