This work describes a detailed mathematical procedure in relation to a novel third‐order WENO scheme for the inviscid term of a system of nonlinear equations in the generalized grid system. The scheme developed minimizes the linear and nonlinear sources of dissipation error associated with the classical fifth‐order WENO scheme. The former is minimized by optimizing the resolving efficiency of the scheme whereas the latter is minimized by fixing the accuracy at the second‐order critical point via redefining the nonlinear weights. Moreover, the spectral property of second‐order viscous derivative, approximated by the single and double applications of the standard fourth‐order central finite difference scheme, is presented. The two‐dimensional Euler and Navier–Stokes equations in the generalized grids are mainly pursued. For the robustness in terms of capturing discontinuous and smooth structures, particularly two problems, which are difficult to handle in Cartesian grids, are chosen for discussion. The first one deals with a supersonic shock hitting the circular cylinder and generating all the possible flow inconsistencies. The other one deals with a subsonic flow over a circular cylinder at the incompressible limit. The numerical results are found to be in good agreement with the experimental data.
- NSF-PAR ID:
- 10184460
- Date Published:
- Journal Name:
- Nonlinear Dynamics
- ISSN:
- 0924-090X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
The strong asymptotic stabilization of 3D hyperbolic dynamics is achieved by a damped 2D elastic structure. The model is a Neumann wave-type equation with low regularity coupling conditions given in terms of a nonlinear von Karman plate. This problem is motivated by the elimination of aeroelastic instability (sustained oscillations of bridges, airfoils, etc.) in engineering applications. Empirical observations indicate that the subsonic wave-plate system converges to equilibria. Classical approaches which decouple the plate and wave dynamics have fallen short. Here, we operate on the model as it appears in the engineering literature with no regularization and achieve stabilization by microlocalizing the Neumann boundary data for the wave equation (given through the plate dynamics). We observe a compensation by the plate dynamics precisely where the regularity of the 3D Neumann wave is compromised (in the characteristic sector).more » « less
-
At Cal Poly Humboldt, undergraduate researchers and faculty have constructed a torsion-pendulum experiment that seeks to measure gravitational interactions below test mass separations of 100 microns. The aim of this experiment is to look for deviations in the weak equivalence principle (WEP) and inverse-square law (ISL). The scale at which this experiment operates is within an untested range at the submillimeter scale. This apparatus’s torsion pendulum consists of equal masses with differing materials arranged as a composition dipole. The twist of this configuration is measured as an attractor mass oscillates in a parallel-plate arrangement nearby. The oscillation creates a time-dependent torque on the pendulum which can be studied for deviations in the WEP and ISL. At present, an active leveling scheme has been implemented to mediate the apparatus’s long-term tilt variations. This scheme has been optimized through the use of a power supply and proportional-integral-derivative (PID) loop that mitigates the variations in tilt by applying a voltage to a resistor attached to one of the apparatus legs. The applied voltage causes thermal expansion of the apparatus leg support structure, thus correcting and modulating the tilt of this experiment.more » « less
-
more » « less
The preBötzinger Complex (preBötC) encodes inspiratory time as rhythmic bursts of activity underlying each breath. Spike synchronization throughout a sparsely connected preBötC microcircuit initiates bursts that ultimately drive the inspiratory motor patterns. Using minimal microcircuit models to explore burst initiation dynamics, we examined the variability in probability and latency to burst following exogenous stimulation of a small subset of neurons, mimicking experiments. Among various physiologically plausible graphs of 1000 excitatory neurons constructed using experimentally determined synaptic and connectivity parameters, directed Erdős-Rényi graphs with a broad (lognormal) distribution of synaptic weights best captured the experimentally observed dynamics. preBötC synchronization leading to bursts was regulated by the efferent connectivity of spiking neurons that are optimally tuned to amplify modest preinspiratory activity through input convergence. Using graph-theoretic and machine learning-based analyses, we found that input convergence of efferent connectivity at the next-nearest neighbor order was a strong predictor of incipient synchronization. Our analyses revealed a crucial role of synaptic heterogeneity in imparting exceptionally robust yet flexible preBötC attractor dynamics. Given the pervasiveness of lognormally distributed synaptic strengths throughout the nervous system, we postulate that these mechanisms represent a ubiquitous template for temporal processing and decision-making computational motifs.
SIGNIFICANCE STATEMENT Mammalian breathing is robust, virtually continuous throughout life, yet is inherently labile: to adapt to rapid metabolic shifts (e.g., fleeing a predator or chasing prey); for airway reflexes; and to enable nonventilatory behaviors (e.g., vocalization, breathholding, laughing). Canonical theoretical frameworks—based on pacemakers and intrinsic bursting—cannot account for the observed robustness and flexibility of the preBötzinger Complex rhythm. Experiments reveal that network synchronization is the key to initiate inspiratory bursts in each breathing cycle. We investigated preBötC synchronization dynamics using network models constructed with experimentally determined neuronal and synaptic parameters. We discovered that a fat-tailed (non-Gaussian) synaptic weight distribution—a manifestation of synaptic heterogeneity—augments neuronal synchronization and attractor dynamics in this vital rhythmogenic network, contributing to its extraordinary reliability and responsiveness. -
more » « less
Abstract Subduction zone accretionary prisms are commonly modeled as elastic structures where permanent deformation is accommodated by faulting and folding of otherwise elastic materials, yet accretionary prisms may exhibit other deformation styles over relatively short time scales. In this study, we use 6.5‐year (2014–2021) Sentinel‐1 interferometric synthetic aperture radar (InSAR) time‐series of post‐seismic deformation in the Makran accretionary prism of southeast Pakistan to characterize non‐linear viscoelastic deformation within an active accretionary prism on short timescales (months to years). We constructed a series of 3‐D finite‐element models of the Makran subduction zone, including an accretionary prism, and constrained the elastic thickness of the upper wedge and the flow‐law parameters (power‐law exponent, activation enthalpy, and pre‐exponential constant) of the lower wedge through forward model fits to the InSAR time‐series. Our results show that the prism is elastically thin (8–12 km) and the non‐linear viscoelastic relaxation of the deep portions of the prism alone can sufficiently explain the post‐seismic surface deformation. Our best fitting flow‐law parameters (
n = 3.76 ± 0.39,Q = 82.2 ± 37.73 kJ mol−1, andA = 10−3.36±4.69) are consistent with triggering of low temperature dislocation creep within fluid‐saturated siliciclastic rocks. We believe that the fluids necessary for this weakening originate from sedimentary underplating and/or the presence the hydrocarbons. The presence of power‐law rheology within the lower wedge impacts the estimated plate coupling and the stress state in the subduction system, with respect to the conventional elastic wedge model, and hence should to be considered in future earthquake cycle models.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s11071-020-05805-1
