In this paper, we discuss the stability and error estimates of the fully discrete schemes for linear conservation laws, which consists of an arbitrary Lagrangian–Eulerian discontinuous Galerkin method in space and explicit total variation diminishing Runge–Kutta (TVD-RK) methods up to third order accuracy in time. The scaling arguments and the standard energy analysis are the key techniques used in our work. We present a rigorous proof to obtain stability for the three fully discrete schemes under suitable CFL conditions. With the help of the reference cell, the error equations are easy to establish and we derive the quasi-optimal error estimates in space and optimal convergence rates in time. For the Euler-forward scheme with piecewise constant elements, the second order TVD-RK method with piecewise linear elements and the third order TVD-RK scheme with polynomials of any order, the usual CFL condition is required, while for other cases, stronger time step restrictions are needed for the results to hold true. More precisely, the Euler-forward scheme needs τ ≤ ρh 2 and the second order TVD-RK scheme needs $ \tau \le \rho {h}^{\frac{4}{3}}$ for higher order polynomials in space, where τ and h are the time and maximum space step, respectively, and ρ is a positive constant independent of τ and h .
more »
« less
Novel conservative methods for adaptive force softening in collisionless and multispecies N -body simulations
Modelling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, and planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multiphysics simulations where species with different ‘softening laws’ may interact. We therefore derive a generalized form of the energy-and-momentum conserving gravitational equations of motion, applicable to arbitrary rules used to determine the force softening, together with consistent associated time-step criteria, interaction terms between species with different softening laws, and arbitrary maximum/minimum softenings. We also derive new methods to maintain better accuracy and conservation when symmetrizing forces between particles. We review and extend previously discussed adaptive softening schemes based on the local neighbour particle density, and present several new schemes for scaling the softening with properties of the gravitational field, i.e. the potential or acceleration or tidal tensor. We show that the ‘tidal softening’ scheme not only represents a physically motivated, translation and Galilean invariant and equivalence-principle respecting (and therefore conservative) method but also imposes negligible time-step or other computational penalties, ensuring that pairwise two-body scattering is small compared to smooth background forces and can resolve outstanding challenges in properly capturing tidal disruption of substructures (minimizing artificial destruction) while also avoiding excessive N-body heating. We make all of this public in the GIZMO code.
more » « less- NSF-PAR ID:
- 10464515
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 525
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 5951-5977
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Lifting line theory describes the cumulative effect of shed vorticity from finite span lifting surfaces. In this work, the theory is reformulated to improve the accuracy of the actuator line model (ALM). This model is a computational tool used to represent lifting surfaces, such as wind-turbine blades in computational fluid dynamics. In ALM, blade segments are represented by means of a Gaussian body force distribution with a prescribed kernel size. Prior analysis has shown that a representation of the blade using an optimal kernel width $\unicode[STIX]{x1D716}^{opt}$ of approximately one quarter of the chord size results in accurate predictions of the velocity field and loads along the blades. Also, simulations have shown that use of the optimal kernel size yields accurate representation of the tip-vortex size and the associated downwash resulting in accurate predictions of the tip losses. In this work, we address the issue of how to represent the effects of finite span wings and tip vortices when using Gaussian body forces with a kernel size larger than the optimal value. This question is relevant in the context of coarse-scale large-eddy simulations that cannot afford the fine resolutions required to resolve the optimal kernel size. For this purpose, we present a filtered lifting line theory for a Gaussian force distribution. Based on the streamwise component of the vorticity transport equation, we develop an analytical model for the induced velocity resulting from the spanwise changes in lift force for an arbitrary kernel scale. The results are used to derive a subfilter-scale velocity model that is used to correct the velocity along the blade when using kernel sizes larger than $\unicode[STIX]{x1D716}^{opt}$ . Tests are performed in large-eddy simulation of flow over fixed wings with constant and elliptic chord distributions using various kernel sizes. Results show that by using the proposed subfilter velocity model, kernel-size independent predictions of lift coefficient and total lift forces agree with those obtained with the optimal kernel size.more » « less
-
more » « less
Abstract The evolution of complex dentitions in mammals was a major innovation that facilitated the expansion into new dietary niches, which imposed selection for tight form–function relationships. Teeth allow mammals to ingest and process food items by applying forces produced by a third-class lever system composed by the jaw adductors, the cranium, and the mandible. Physical laws determine changes in jaw adductor (biting) forces at different bite point locations along the mandible (outlever), thus, individual teeth are expected to experience different mechanical regimes during feeding. If the mammal dentition exhibits functional adaptations to mandible feeding biomechanics, then teeth are expected to have evolved to develop mechanically advantageous sizes, shapes, and positions. Here, we present bats as a model system to test this hypothesis and, more generally, for integrative studies of mammal dental diversity. We combine a field-collected dataset of bite forces along the tooth row with data on dental and mandible morphology across 30 bat species. We (1) describe, for the first time, bite force trends along the tooth row of bats; (2) use phylogenetic comparative methods to investigate relationships among bite force patterns, tooth, and mandible morphology; and (3) hypothesize how these biting mechanics patterns may relate to the developmental processes controlling tooth formation. We find that bite force variation along the tooth row is consistent with predictions from lever mechanics models, with most species having the greatest bite force at the first lower molar. The cross-sectional shape of the mandible body is strongly associated with the position of maximum bite force along the tooth row, likely reflecting mandibular adaptations to varying stress patterns among species. Further, dental dietary adaptations seem to be related to bite force variation along molariform teeth, with insectivorous species exhibiting greater bite force more anteriorly, narrower teeth and mandibles, and frugivores/omnivores showing greater bite force more posteriorly, wider teeth and mandibles. As these craniodental traits are linked through development, dietary specialization appears to have shaped intrinsic mechanisms controlling traits relevant to feeding performance.
-
more » « less
Abstract During behavior, the work done by actuators on the body can be resisted by the body's inertia, elastic forces, gravity, or viscosity. The dominant forces that resist actuation have major consequences on the control of that behavior. In the literature, features and actuation of locomotion, for example, have been successfully predicted by nondimensional numbers (e.g. Froude number and Reynolds number) that generally express the ratio between two of these forces (gravitational, inertial, elastic, and viscous). However, animals of different sizes or motions at different speeds may not share the same dominant forces within a behavior, making ratios of just two of these forces less useful. Thus, for a broad comparison of behavior across many orders of magnitude of limb length and cycle period, a dimensionless number that includes gravitational, inertial, elastic, and viscous forces is needed. This study proposes a nondimensional number that relates these four forces: the phase shift (ϕ) between the displacement of the limb and the actuator force that moves it. Using allometric scaling laws, ϕ for terrestrial walking is expressed as a function of the limb length and the cycle period at which the limb steps. Scale-dependent values of ϕ are used to explain and predict the electromyographic (EMG) patterns employed by different animals as they walk.
-
more » « less
Abstract We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure formation simulations of the universe. The new framework couples the HACC gravitational
N -body solver with a modern smoothed-particle hydrodynamics (SPH) approach called conservative reproducing kernel SPH (CRKSPH). CRKSPH utilizes smoothing functions that exactly interpolate linear fields while manifestly preserving conservation laws (momentum, mass, and energy). The CRKSPH method has been incorporated to accurately model baryonic effects in cosmology simulations—an important addition targeting the generation of precise synthetic sky predictions for upcoming observational surveys. CRK-HACC inherits the codesign strategies of the HACC solver and is built to run on modern GPU-accelerated supercomputers. In this work, we summarize the primary solver components and present a number of standard validation tests to demonstrate code accuracy, including idealized hydrodynamic and cosmological setups, as well as self-similarity measurements.
