Search

We show that a universe with a nonminimally coupled scalar field can fit current measurements of the expansion rate of the Universe better than the standard $\Lambda$-Cold Dark Matter model or other minimally coupled dark energy models. While we find a clear improvement in the goodness of fit for this dark energy model with respect to others that have been considered in the recent literature, using various information theoretic criteria, we show that the evidence for it is still inconclusive., Comment: Updated version. Minor changes to the discussion and references added. 5 pages, 2 figures

New constraints on the expansion rate of the Universe seem to favor evolving dark energy in the form of thawing quintessence models, i.e., models for which a canonical, minimally coupled scalar field has, at late times, begun to evolve away from potential energy domination. We scrutinize the evidence for thawing quintessence by exploring what it predicts for the equation of state. We show that, in terms of the usual Chevalier-Polarski-Linder parameters, ($w_0$, $w_a$), thawing quintessence is, in fact, only marginally consistent with a compilation of the current data. Despite this, we embrace the possibility that thawing quintessence is dark energy and find constraints on the microphysics of this scenario. We do so in terms of the effective mass $m^2$ and energy scale $V_0$ of the scalar field potential. We are particularly careful to enforce un-informative, flat priors on these parameters so as to minimize their effect on the final posteriors. While the current data favors a large and negative value of $m^2$, when we compare these models to the standard $\Lambda$CDM model we find that there is scant evidence for thawing quintessence., Comment: Accepted and forthcoming in PRD. Minor changes to discussions. 18 pages, 14 figures

We revisit a class of simple single-field inflation models and demonstrate that they can readily produce a negligible tensor/scalar ratio $r$. Motivated by recent work suggesting the need to introduce higher order operators to stabilise unregulated potentials, as well as by work indicating that such terms can have significant effects on observable predictions, we explicitly construct corrected versions of the quadratic hilltop potential that are motivated by an effective field theory expansion. We employ Markov Chain Monte Carlo (MCMC) methods and optimization techniques to sample viable models and minimize $r$. We find that such potentials can readily lower $r$ values below projected CMB-S4 sensitivity, while still remaining within observable constraints on $n_s$. Furthermore, we find that the minimum $r$ reached for each order of the expansion considered is well-described by a power law $r_{min}(q) \propto q^{-B}$ before asymptoting to a value of $r_{min} \sim 10^{-11}$, where $q$ is the order to which the expansion of $V(\phi)$ is carried out., Comment: Accepted for publication in PRD. Minor adjustmends to the discussion and references added

The paper argues that we ought to conceive of the Dark Energy problem -- the question of how to account for observational data, naturally interpreted as accelerated expansion of the universe -- as a crisis of underdetermined pursuit-worthiness. Not only are the various approaches to the Dark Energy problem evidentially underdetermined; at present, no compelling reasons single out any of them as more likely to be true than the other. More vexingly for working scientists, none of the approaches stands out as uncontroversially preferable over its rivals in terms of its rationally warranted promise, i.e. the reasons to further work on, explore and develop it. We demonstrate this claim by applying a Peircean economic model of pursuit-worthiness in terms of a cognitive cost/benefit estimate -- with the instantiation of theory virtues as key indicators of cognitive gains -- to the four main Dark Energy proposals (the cosmological constant approach, modified gravity, quintessence, and inhomogeneous cosmologies). Our analysis yields that these approaches do not admit of an unambiguous, or uncontroversial, ranking with respect to which ansatz deserves distinguished attention and research efforts. The overall methodological counsel that our analysis underwrites recommends a pragmatic double research strategy forward: to encourage and foster theory pluralism and the search for tests -- with the goal of enhancing the testability of the $\Lambda$CDM model and "testing it to destruction"., Comment: Preprint. Comments welcome!

Regular, automated testing is a foundational principle of modern software development. Numerous widely-used continuous integration systems exist, but they are often not suitable for the unique needs of scientific simulation software. Here we describe the testing infrastructure developed for and used by the Modules for Experiments in Stellar Astrophysics (MESA) project. This system allows the computationally-demanding MESA test suite to be regularly run on a heterogeneous set of computers and aggregates and displays the testing results in a form that allows for the rapid identification and diagnosis of regressions. Regularly collecting comprehensive testing data also enables longitudinal studies of the performance of the software and the properties of the models it generates., Comment: 12 page, 7 figures, Accepted to ApJS

There is compelling evidence that the Universe is undergoing a late phase of accelerated expansion. One of the simplest explanations for this behaviour is the presence of dark energy. A plethora of microphysical models for dark energy have been proposed. The hope is that, with the ever increasing precision of cosmological surveys, it will be possible to precisely pin down the model. We show that this is unlikely and that, at best, we will have a phenomenological description for the microphysics of dark energy. Furthermore, we argue that the current phenomenological prescriptions are ill-equipped for shedding light on the fundamental theory of dark energy., Comment: Updated the figures to have higher resolution. 14 pages, 8 figures

The paper investigates the historical and contemporary pursuit-worthiness of cosmic inflation-the rationale for working on it (rather than necessarily the evidential support for claims to its approximate truth): what reasons existed, and exist, that warrant inflation's status as the mainstream paradigm studied, explored, and further developed by the majority of the cosmology community? We'll show that inflation exemplifies various salient theory virtues: explanatory depth, unifying/integrative power, fertility and positive heuristics, the promotion of understanding, and the prospect (and passing) of novel benchmark tests. This, we'll argue, constitutes inflation's auspicious promise. It marks inflation as preferable over both the inflation-less Hot Big Bang Model, as well as rivals to inflation: inflation, we maintain, rightly deserved, and continues to deserve, the concerted research efforts it has enjoyed., Comment: Forthcoming in British Journal for the Philosophy of Science. Updated to fix minor typos

We show that the dynamical common core of the recently-discovered non-relativistic geometric trinity of gravity is Maxwell gravitation. Moreover, we explain why no analogous distinct dynamical common core exists in the case of the better-known relativistic geometric trinity of gravity., Comment: Forthcoming in Philosophy of Physics

I apply Dawid's Meta-Empirical Assessment (MEA) methodology to the theory of cosmological inflation. I argue that applying this methodology does not currently offer a compelling case for ascribing non-empirical confirmation to cosmological inflation. In particular, I argue that despite displaying strong instances of Unexpected Explanatory Coherence (UEA), it is premature to evaluate the theory on the basis of the No Alternatives Argument (NAA). More significantly though, I argue that the theory of cosmological inflation fails to sustain a convincing Meta-Inductive Argument (MIA) because the empirical evidence and theoretical successes that it seeks to draw meta-empirical support from do not warrant a meta-inductive inference to inflation. I conclude by assessing how future developments could pave the way towards crafting a more compelling case for the non-empirical confirmation of cosmological inflation., Comment: Accepted in Studies in History and Philosophy of Science

This work presents a novel methodology for analysis and control of nonlinear fluid systems using neural networks. The approach is demonstrated on four different study cases being the Lorenz system, a modified version of the Kuramoto-Sivashinsky equation, a streamwise-periodic 2D channel flow, and a confined cylinder flow. Neural networks are trained as models to capture the complex system dynamics and estimate equilibrium points through a Newton method, enabled by backpropagation. These neural network surrogate models (NNSMs) are leveraged to train a second neural network, which is designed to act as a stabilizing closed-loop controller. The training process employs a recurrent approach, whereby the NNSM and the neural network controller (NNC) are chained in closed loop along a finite time horizon. By cycling through phases of combined random open-loop actuation and closed-loop control, an iterative training process is introduced to overcome the lack of data near equilibrium points. This approach improves the accuracy of the models in the most critical region for achieving stabilization. Through the use of L1 regularization within loss functions, the NNSMs can also guide optimal sensor placement, reducing the number of sensors from an initial candidate set. The datasets produced during the iterative training process are also leveraged for conducting a linear stability analysis through a modified dynamic mode decomposition approach. The results demonstrate the effectiveness of computationally inexpensive neural networks in modeling, controlling, and enabling stability analysis of nonlinear systems, providing insights into the system behaviour and offering potential for stabilization of complex fluid systems., Comment: 30 pages, 22 figures, under consideration for publication

Modal decomposition techniques are important tools for the analysis of unsteady flows and, in order to provide meaningful insights with respect to coherent structures and their characteristic frequencies, the modes must possess a robust spatial support. In this context, although widely used, methods based on singular value decomposition (SVD) may produce modes that are difficult to interpret when applied to problems dominated by intermittent and transient events. Fortunately, specific modal decomposition techniques have been recently developed to analyze such problems. However, a proper comparison between existing methods is still lacking from the literature. Therefore, this work compares two recent methods: the fast adaptive multivariate empirical mode decomposition (FA-MVEMD) and the multi-resolution dynamic mode decomposition (mrDMD). These techniques are employed here for the study of flow databases involving transient and intermittent dynamics. Specifically, the investigated problems include an SD7003 airfoil subjected to deep dynamic stall conditions, and a steady NACA0012 airfoil operating at a transitional Reynolds number. In the former case, the methods are employed to investigate the onset and evolution of the dynamic stall vortex (DSV), while in the latter case, intermittent vortex pairing is analyzed. We show that the combination of a multidimensional EMD with the Hilbert transform provides modes with superior spatial support when compared to the mrDMD, also allowing the characterization of instantaneous frequencies of coherent structures. Moreover, the EMD also condenses a larger amount of information within a single intrinsic mode function (IMF).

The geometric trinity of gravity comprises three distinct formulations of general relativity: (i) the standard formulation describing gravity in terms of spacetime curvature, (ii) the teleparallel equivalent of general relativity describing gravity in terms of spacetime torsion, and (iii) the symmetric teleparallel equivalent of general relativity (STEGR) describing gravity in terms of spacetime non-metricity. In this article, we complete a geometric trinity of non-relativistic gravity, by (a) taking the non-relativistic limit of STEGR to determine its non-relativistic analogue, and (b) demonstrating that this non-metric theory is equivalent to the Newton--Cartan theory and its teleparallel equivalent, i.e., the curvature and the torsion based non-relativistic theories that are both geometrised versions of classical Newtonian gravity., Comment: Accepted version to appear in General Relativity and Gravitation. Updated with minor changes to the discussion

Canonically, `classic' tests of general relativity (GR) include perihelion precession, the bending of light around stars, and gravitational redshift; `modern' tests have to do with, inter alia, relativistic time delay, equivalence principle tests, gravitational lensing, strong field gravity, and gravitational waves. The orthodoxy is that both classic and modern tests of GR afford experimental confirmation of that theory in particular. In this article, we question this orthodoxy, by showing there are classes of both relativistic theories (with spatiotemporal geometrical properties different from those of GR) and non-relativistic theories (in which the lightcones of a relativistic spacetime are `widened') which would also pass such tests. Thus, (a) issues of underdetermination in the context of GR loom much larger than one might have thought, and (b) given this, one has to think more carefully about what exactly such tests in fact are testing., Comment: This is a revised version of the previous draft

The paper re-examines the principal methodological questions, arising in the debate over the cosmological standard model's postulate of Dark Matter vs. rivalling proposals that modify standard (Newtonian and general-relativistic) gravitational theory, the so-called Modified Newtonian Dynamics (MOND) and its subsequent extensions. What to make of such seemingly radical challenges of cosmological orthodoxy? In the first part of our paper, we assess MONDian theories through the lens of key ideas of major 20th century philosophers of science (Popper, Kuhn, Lakatos, and Laudan), thereby rectifying widespread misconceptions and misapplications of these ideas common in the pertinent MOND-related literature. None of these classical methodological frameworks, which render precise and systematise the more intuitive judgements prevalent in the scientific community, yields a favourable verdict on MOND and its successors -- contrary to claims in the MOND-related literature by some of these theories' advocates; the respective theory appraisals are largely damning. Drawing on these insights, the paper's second part zooms in on the most common complaint about MONDian theories, their ad-hocness. We demonstrate how the recent coherentist model of ad-hocness captures, and fleshes out, the underlying -- but too often insufficiently articulated -- hunches underlying this critique. MONDian theories indeed come out as severely ad hoc: they do not cohere well with either theoretical or empirical-factual background knowledge. In fact, as our complementary comparison with the cosmological standard model's Dark Matter postulate shows, with respect to ad-hocness, MONDian theories fare worse than the cosmological standard model., Comment: forthcoming in Studies in History and Philosophy of Science

A standard line in the contemporary philosophical literature has it that physical theories are equivalent only when they agree on their empirical content, where this empirical content is often understood as being encoded in the equations of motion of those theories. In this article, we question whether it is indeed the case that the empirical content of a theory is exhausted by its equations of motion, showing that (for example) considerations of boundary conditions play a key role in the empirical equivalence (or otherwise) of theories. Having argued for this, we show that philosophical claims made by Knox (2011) that general relativity is equivalent to teleparallel gravity, and by Weatherall (2016) that electromagnetism in the Faraday tensor formalism is equivalent to electromagnetism in the vector potential formalism, can both be called into question. We then show that properly considering the role of boundary conditions in theory structure can potentially restore these claims of equivalence and close with some remarks on the pragmatics of adjudications on theory identity., Comment: Updated to match published version

In this work, we propose a framework whereby flow imaging data is leveraged to extract relevant information from flowfield visualizations. To this end, a vision transformer (ViT) model is developed to predict the unsteady pressure distribution over an airfoil under dynamic stall from images of the flowfield. The network is capable of identifying relevant flow features present in the images and associate them to the airfoil response. Results demonstrate that the model is effective in interpolating and extrapolating between flow regimes and for different airfoil motions, meaning that ViT-based models may offer a promising alternative for sensors in experimental campaigns and for building robust surrogate models of complex unsteady flows. In addition, we uniquely treat the image semantic segmentation as an image-to-image translation task that infers semantic labels of structures from the input images in a supervised way. Given an input image of the velocity field, the resulting convolutional neural network (CNN) generates synthetic images of any corresponding fluid property of interest. In particular, we convert the velocity field data into pressure in order to subsequently estimate the pressure distribution over the airfoil in a robust manner.

Current design constraints have encouraged the studies of aeroacoustic fields around compressible jet flows. The present work addresses the numerical study of unsteady turbulent jet flows as a preparation for future aeroacoustic analyses of main engine rocket plumes. An in-house large eddy simulation tool is used in order to reproduce high fidelity results of compressible jet flows. The large eddy simulation formulation is written using a second order numerical scheme for a finite difference spatial discretization. Numerical simulations of perfectly expanded jets are performed and the results are compared to the literature. Dynamic mode decompositions (DMD) of the jet flow, using large size three-dimensional snapshots, are performed. Three variables are analyzed, namely, the velocity magnitude, the vorticity magnitude and the divergence of velocity. In particular, two frequencies are identified and they are linked to flow structures observed in experiments performed by other authors in the literature. The spatial shapes of the corresponding dynamic modes are also discussed., Comment: AIAA aviation paper 2017; AIAA 2017-0497. 55th AIAA Aerospace Sciences Meeting. January 2017. arXiv admin note: substantial text overlap with arXiv:2301.00651; text overlap with arXiv:2212.12365, arXiv:2212.12353, arXiv:2301.00817

Current design constraints have encouraged the studies of aeroacoustics fields around compressible jet flows. The present work addresses the numerical study of subgrid scale modeling for unsteady turbulent jet flows as a preliminary step for future aeroacoustic analyses of main engine rocket plumes. An in-house large eddy simulation (LES) tool is developed in order to reproduce high fidelity results of compressible jet flows. In the present study, perfectly expanded jets are considered because the authors want to emphasize the effects of the jet mixing phenomena. The large eddy simulation formulation is written using the finite difference approach, with an explicit time integration and using a second order spatial discretization. The energy equation is carefully discretized in order to model the energy equation of the filtered Navier-Stokes formulation. The classical Smagorinsky model, the dynamic Smagorinsky model and the Vreman models are the chosen subgrid scale closures for the present work. Numerical simulations of perfectly expanded jets are performed and compared with the literature in order to validate and compare the performance of each subgrid closure in the solver., Comment: conference paper - AIAA Aviation 2016. arXiv admin note: substantial text overlap with arXiv:2212.12353, arXiv:2212.12365

An in-house large eddy simulation tool is developed in order to reproduce high fidelity results of compressible jet flows. The large eddy simulation formulation is written using the finite difference approach, with an explicit time integration and using a second order spatial discretization. The energy equation is carefully discretized in order to model the energy equation of the filtered Navier-Stokes formulation. Such nu\-me\-ri\-cal studies are very expensive and demand high performance computing. Message passage interface protocols are implemented into the code in order to perform parallel computations. The present work addresses the computational performance of the solver running on up to 400 processors in parallel. Different mesh configurations, whose size varies from approximately 5.9 million points to approximately 1.0 billion points, are evaluate in the current paper. Speedup and efficiency curves are evaluated in order to assess the strong scalability of the solver., Comment: AIAA aviation conference 2017. arXiv admin note: substantial text overlap with arXiv:2301.00651, arXiv:2212.12365, arXiv:2212.12353

Current design constraints have encouraged the studies of aeroacoustics fields around compressible jet flows. The present work addresses the numerical study of unsteady turbulent jet flows for aeroacoustic analyses of main engine rocket plumes. A novel large eddy simulation (LES) tool is developed in order to reproduce high fidelity results of compressible jet flows which could be used for aeroacoustic studies with the Ffowcs Williams and Hawkings approach. The numerical solver is an upgrade of an existing Reynolds-averaged Navier-Stokes solver previously developed in the group. The original framework is rewritten in a modern fashion and intensive parallel computation capabilities have been added to the code. The LES formulation is written using the finite difference approach. The energy equation is carefully discretized in order to model the energy equation of the filtered Navier-Stokes formulation. The classical Smagorinsky model is the chosen subgrid scale closure for the present work. Numerical simulations of perfectly expanded jets are performed and compared with the literature in order to validate the new solver. Moreover, speedup and the computational performance of the code are evaluated and discussed. Flow results are used for an initial evaluation of the noise radiated from the rocket plume., Comment: 25 pages conference paper

The present work is concerned with a study of large eddy simulations (LES) of unsteady turbulent jet flows. In particular, the present analysis is focused on the effects of the subgrid-scale modeling used when a second-order spatial discretization methodology is employed for the numerical simulations. The present effort addresses perfectly expanded supersonic jets, because the authors want to emphasize the effects of the jet mixing phenomena. The LES formulation is discretized using the finite difference approach, after the equations are rewritten in a generalized coordinate system. Both space and time discretizations are second-order accurate and an explicit time march is adopted. Special care is dedicated to the discretization of the energy equation to appropriately model the set of filtered equations appearing in the LES formulation. The classical Smagorinsky, the dynamic Smagorinsky and the Vreman models are the subgrid-scale closures selected for the present work. The computational results are compared to data in the literature to validate the present simulation tool. Results indicate that the characteristics of numerical discretization can be as important as the effects of the subgrid-scale models for such low-order spatial discretization schemes. A detailed analysis is presented for the performance of each subgrid closure in the numerical context here considered., Comment: 44 pages article

We develop and apply a multi-dimensional account of explanatory depth towards a comparative analysis of inflationary and bouncing paradigms in primordial cosmology. Our analysis builds on earlier work due to Azhar and Loeb (2021) that establishes initial conditions fine-tuning as a dimension of explanatory depth relevant to debates in contemporary cosmology. We propose dynamical fine-tuning and autonomy as two further dimensions of depth in the context of problems with instability and trans-Planckian modes that afflict bouncing and inflationary approaches respectively. In the context of the latter issue, we argue that the recently formulated trans-Planckian censorship conjecture leads to a trade-off for inflationary models between dynamical fine-tuning and autonomy. We conclude with the suggestion that explanatory preference with regard to the different dimensions of depth is best understood in terms of differing attitudes towards heuristics for future model building., Comment: fixed typos in accepted version. Forthcoming in The British Journal for the Philosophy of Science (2023)

Large eddy simulations are performed to study tonal noise generation by a NACA0012 airfoil at an angle of attack ${\alpha = 3}$ deg. and freestream Mach number of ${M_{\infty} = 0.3}$. Different Reynolds numbers are analyzed spanning ${0.5 \times 10^5 \le Re \le 4 \times 10^5}$. Results show that the flow patterns responsible for noise generation appear from different laminar separation bubbles, including one observed over the airfoil suction side and another near the trailing edge, on the pressure side. For lower Reynolds numbers, intermittent vortex dynamics on the suction side results in either coherent structures or turbulent packets advected towards the trailing edge. Such flow dynamics also affects the separation bubbles on the pressure side, which become intermittent. Despite the irregular occurrence of laminar-turbulent transition, the noise spectrum depicts a main tone with multiple equidistant secondary tones. Increasing the Reynolds number leads to a permanent turbulent regime on the suction side that reduces the coherence level causing only small scale turbulent eddies to be observed. Furthermore, the laminar separation bubble on the suction side almost vanishes while that on the pressure side becomes more pronounced and permanent. As a consequence, the dominant noise generation mechanism becomes the vortex shedding along the wake., Comment: 31 pages, 19 figures

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement., Comment: 50 pages, 29 figures; Accepted to ApJS

The physics of shock-boundary layer interactions (SBLIs) in a supersonic turbine cascade at Mach 2.0 and Reynolds number 395,000, based on the axial chord, is investigated through a wall-resolved large eddy simulation. Special attention is given to the characterization of the low-frequency dynamics of the separation bubbles using flow visualization, spectral analysis, space-time cross correlations, and flow modal decomposition. The mean flowfield shows different shock structures formed on both sides of the airfoil. On the suction side, an oblique shock impinges on the turbulent boundary layer, whereas a Mach reflection interacts with the pressure side boundary layer. The interactions taking place in the present turbine cascade show similarities and discrepancies with respect to more canonical cases. For example, the characteristic frequencies of the shock/bubble motions are comparable to those described in the literature of canonical cases. However, the suction side bubble leads to compression waves that do not coalesce into a separation shock, and a thin bubble forms on the pressure side despite the strong normal shock from the Mach reflection. Instantaneous flow visualizations illustrate elongated streamwise structures on the incoming boundary layers and their interactions with the shocks and separation bubbles. The space-time cross-correlations reveal that the near-wall streaks drive the motion of the suction side separation bubble, which in turn promotes oscillations of the reattachment shock and shear layer flapping. Organized motions in the SBLIs and their corresponding characteristic frequencies and spatial support are identified using proper orthogonal decomposition., Comment: 47 pages, 25 figures. Published in the Physical Review Fluids (https://journals.aps.org/prfluids/abstract/10.1103/PhysRevFluids.7.094602)

It is now well-known that Newton-Cartan theory is the correct geometrical setting for modelling the quantum Hall effect. In addition, in recent years edge modes for the Newton-Cartan quantum Hall effect have been derived. However, the existence of these edge modes has, as of yet, been derived using only orthodox methodologies involving the breaking of gauge-invariance; it would be preferable to derive the existence of such edge modes in a gauge-invariant manner. In this article, we employ recent work by Donnelly and Freidel in order to accomplish exactly this task. Our results agree with known physics, but afford greater conceptual insight into the existence of these edge modes: in particular, they connect them to subtle aspects of Newton-Cartan geometry and pave the way for further applications of Newton-Cartan theory in condensed matter physics., Comment: 23 pages, forthcoming in Foundations of Physics

This work is motivated by CFD simulations from a realistic landing gear performed modeling only the half bottom of the aircraft fuselage [15]. Hence, in this previous analysis, the aeroacoustic predictions had to employ incomplete FWH surfaces. In the present study, a discussion on the closed surface requirement is provided for the frequency domain FWH formulation, although we believe that the approach presented here can also be applied for the time-domain methodology. We show that all sources computed on the closed permeable surface may play an important role if the acoustic prediction is sought at observer positions distributed along the entire circular arc. However, results obtained for 2D and 3D model problems show that, if wisely designed, a permeable FWH surface composed only by a finite patch, i.e., an incomplete surface, can be employed to accurately predict the noise at specific observer positions. In this case, the surface must be placed between the line of sight from the sources to the observers, for example, in flyover or sideline locations. In the current applications, results from finite surfaces agree with the expected values as long as the source magnitudes decay on elements distant from the true incident source, i.e., the airframe. Hence, in order for this approach to be valid, the sources cannot be truncated. The planar setup is beneficial in airframe noise applications since it avoids contamination by quadrupole sources crossing the boundaries and also may simplify the surface design., Comment: 12 pages, 9 figures

Extremum seeking control (ESC) and its slope seeking generalization are applied in a high-fidelity flow simulation framework for reduction of acoustic noise generated by a NACA0012 airfoil. Two Reynolds numbers are studied for which different noise generation mechanisms are excited. For a low Reynolds number flow, the scattering of vortex shedding at the airfoil trailing edge produces tonal noise while, for a moderate Reynolds number case, boundary layer instabilities scatter at the trailing edge leading to noise emission at multiple tones superimposed on a broadband hump. Different control setups are investigated and they are configured to either find an optimal steady actuator intensity or an optimal position for a blowing/suction device. Implementation details are discussed regarding the control modules and design of digital filters., Comment: 39 pages, 32 figures

A large eddy simulation is performed to study secondary tones generated by a NACA0012 airfoil at angle of attack of $\alpha = 3^{\circ}$ with freestream Mach number of $M_{\infty} = 0.3$ and Reynolds number of $Re = 5 \times 10^4$. Laminar separation bubbles are observed over the suction side and near the trailing edge, on the pressure side. Flow visualization and spectral analysis are employed to investigate vortex shedding aft of the suction side separation bubble. Vortex interaction results in merging or bursting such that coherent structures or turbulent packets are advected towards the trailing edge leading to different levels of noise emission. Despite the intermittent occurrence of laminar-turbulent transition, the noise spectrum depicts a main tone with multiple equidistant secondary tones. To understand the role of flow instabilities on the tones, the linearized Navier-Stokes equations are examined in its operator form through bi-global stability and resolvent analyses, and by time evolution of disturbances using a matrix-free method. These linear global analyses reveal amplification of disturbances over the suction side separation bubble. Non-normality of the linear operator leads to further transient amplification due to modal interaction among eigenvectors. Two-point, one time autocovariance calculations of pressure along the spanwise direction elucidate aspects of the acoustic feedback loop mechanism in the non-linear solutions. This feedback process is self-sustained by acoustic waves radiated from the trailing edge, which reach the most sensitive flow location between the leading edge and the separation bubble, as identified by the resolvent analysis. Leading edge disturbances arising from secondary diffraction and phase interference among the most unstable frequencies computed in the eigenspectrum are also shown to have an important role in the feedback loop., Comment: 31 pages

A data-driven closure modeling based on proper orthogonal decomposition (POD) temporal modes is used to obtain stable and accurate reduced order models (ROMs) of unsteady compressible flows. Model reduction is obtained via Galerkin and Petrov-Galerkin projection of the non-conservative compressible Navier-Stokes equations. The latter approach is implemented using the least-squares Petrov-Galerkin (LSPG) technique and the present methodology allows pre-computation of both Galerkin and LSPG coefficients. Closure is performed by adding linear and non-linear coefficients to the original ROMs and minimizing the error with respect to the POD temporal modes. In order to further reduce the computational cost of the ROMs, an accelerated greedy missing point estimation (MPE) hyper-reduction method is employed. A canonical compressible cylinder flow is first analyzed and serves as a benchmark. The second problem studied consists of the turbulent flow over a plunging airfoil undergoing deep dynamic stall. For the first case, linear and non-linear closure coefficients are both low in intrusiveness, capable of providing results in excellent agreement with the full order model. Regularization of calibrated models is also straightforward for this case. On the other hand, the dynamic stall flow is significantly more challenging, specially when only linear coefficients are used. Results show that non-linear calibration coefficients outperform their linear counterparts when a POD basis with fewer modes is used in the reconstruction. However, determining a correct level of regularization is more complicated with non-linear coefficients. Hyper-reduced models show good results when combined with non-linear calibration and an appropriate sized POD basis.

An analysis of calibration for reduced-order models (ROMs) is presented in this work. The Galerkin and least-squares Petrov-Galerkin (LSPG) methods are tested on compressible flows involving a disparity of temporal scales. A novel calibration strategy is proposed for the LSPG method and two test cases are analyzed. The first consists of a subsonic airfoil flow where boundary layer instabilities are responsible for trailing-edge noise generation and the second comprises a supersonic airfoil flow with a transient period where a detached shock wave propagates upstream at the same time that shock-vortex interaction occurs at the trailing edge. Results show that calibration produces stable and long-time accurate for both cases. In order to reduce the computational costs of the LSPG models, an accelerated greedy missing point estimation (MPE) algorithm is employed for hyper-reduction. For the first case investigated, LSPG solutions obtained with hyper-reduction show good comparison with those obtained by the full order model. However, for the supersonic case the transient features of the flow need to be properly captured by the sampled points. Otherwise, the dynamics of the moving shock wave are not fully recovered. The impact of different time-marching schemes is also assessed and, differently than reported in literature, Galerkin models are shown to be more accurate than those computed by LSPG when the non-conservative form of the Navier-Stokes equations are solved. For the supersonic case, the Galerkin and LSPG models (without hyper-reduction) capture the overall dynamics of the detached and oblique shock waves along the airfoil. However, when shock-vortex interaction occurs at the trailing-edge, the Galerkin ROM is able to capture the high-frequency fluctuations from vortex shedding while the LSPG presents a more dissipative solution, not being able to recover the flow dynamics.

Acoustics loads are rocket design constraints which push researches and engineers to invest efforts in the aeroacoustics phenomena which is present on launch vehicles. Therefore, an in-house computational fluid dynamics tool is developed in order to reproduce high-fidelity results of supersonic jet flows for aeroacoustic analogy applications. The solver is written using the large eddy simulation formulation that is discretized using a finite-difference approach and an explicit time integration. Numerical simulations of supersonic jet flows are very expensive and demand efficient high-performance computing. Therefore, non-blocking message passage interface protocols and parallel input/output features are implemented into the code in order to perform simulations which demand up to one billion degrees of freedom. The present work evaluates the parallel efficiency of the solver when running on a supercomputer with a maximum theoretical peak of 127.4 TFLOPS. Speedup curves are generated using nine different workloads. Moreover, the validation results of a realistic flow condition are also presented in the current work., Comment: Journal article from Journal of the Brazilian Society of Mechanical Sciences and Engineering (41, Article number: 547 (2019))

New regulations are imposing noise emissions limitations for the aviation industry which are pushing researchers and engineers to invest efforts in studying the aeroacoustics phenomena. Following this trend, an in-house computational fluid dynamics tool is build to reproduce high fidelity results of supersonic jet flows for aeroacoustic analogy applications. The solver is written using the large eddy simulation formulation that is discretized using a finite difference approach and an explicit time integration. Numerical simulations of supersonic jet flows are very expensive and demand efficient high-performance computing. Therefore, non-blocking message passage interface protocols and parallel Input/Output features are implemented into the code in order to perform simulations which demand up to one billion grid points. The present work addresses the evaluation of code improvements along with the computational performance of the solver running on a computer with maximum theoretical peak of 2.727 PFlops. Different mesh configurations, whose size varies from a few hundred thousand to approximately one billion grid points, are evaluated in the present paper. Calculations are performed using different workloads in order to assess the strong and weak scalability of the parallel computational tool. Moreover, validation results of a realistic flow condition are also presented in the current work., Comment: 13 pages journal article. arXiv admin note: text overlap with arXiv:2003.08746

Cosmological models with a dynamical dark energy field typically lead to a modified propagation of gravitational waves via an effectively time-varying gravitational coupling $G(t)$. The local variation of this coupling between the time of emission and detection can be probed with standard sirens. Here we discuss the role that Lunar Laser Ranging (LLR) and binary pulsar constraints play in the prospects of constraining $G(t)$ with standard sirens. In particular, we argue that LLR constrains the matter-matter gravitational coupling $G_N(t)$, whereas binary pulsars and standard sirens constrain the quadratic kinetic gravity self-interaction $G_{gw}(t)$. Generically, these two couplings could be different in alternative cosmological models, in which case LLR constraints are irrelevant for standard sirens. We use the Hulse-Taylor pulsar data and show that observations are highly insensitive to time variations of $G_{gw}(t)$ yet highly sensitive to $G_N(t)$. We thus conclude that future gravitational waves data will become the best probe to test $G_{gw}(t)$, and will hence provide novel constraints on dynamical dark energy models., Comment: Updated to published version

Direct numerical simulations are carried out to investigate the flow features responsible for secondary tones arising in trailing-edge noise at moderate Reynolds numbers. Simulations are performed for a NACA 0012 airfoil at freestream Mach numbers 0.1, 0.2 and 0.3 for angle of incidence 0 deg. and for Mach number 0.3 at 3 deg. angle of incidence. The Reynolds number based on the airfoil chord is fixed at $Re_c=10^5$. Flow configurations are investigated where noise generation arises from the scattering of boundary layer instabilities at the trailing edge. Results show that noise emission has a main tone with equidistant secondary tones, as discussed in literature. An interesting feature of the present flows at zero incidence is shown; despite the geometric symmetry, the flows become non-symmetric with a separation bubble only on one side of the airfoil. A separation bubble is also observed for the non-zero incidence flow. For both angles of incidence analyzed, it is shown that low-frequency motion of the separation bubbles induce a frequency modulation of the flow instabilities developed along the airfoil boundary layer. When the airfoil is at 0 deg. angle of attack an intense amplitude modulation is also observed in the flow quantities, resulting in a complex vortex interaction mechanism at the trailing edge. Both amplitude and frequency modulations directly affect the velocity and pressure fluctuations that are scattered at the trailing edge, what leads to secondary tones in the acoustic radiation.

We consider cosmological models with a dynamical dark energy field, and study the presence of three types of commonly found instabilities, namely ghost (when fields have negative kinetic energy), gradient (negative momentum squared) and tachyon (negative mass squared). In particular, we study the linear scalar perturbations of theories with two interacting scalar fields as a proxy for a dark energy and matter fields, and explicitly show how canonical transformations relate these three types of instabilities with each other. We generically show that low-energy ghosts are equivalent to tachyonic instabilities, and that high-energy ghosts are equivalent to gradient instabilities. Via examples we make evident the fact that whenever one of these fields exhibits an instability then the entire physical system becomes unstable, with an unbounded Hamiltonian. Finally, we discuss the role of interactions between the two fields, and show that whereas most of the time interactions will not determine whether an instability is present or not, they may affect the timescale of the instability. We also find exceptional cases in which the two fields are ghosts and hence the physical system is seemingly unstable, but the presence of interactions actually lead to stable solutions. These results are very important for assessing the viability of dark energy models that may exhibit ghost, gradient or tachyonic modes., Comment: 11 pages. Updated to match published version

High-fidelity simulations are performed to study active flow control techniques for alleviating deep dynamic stall of a SD7003 airfoil in plunging motion. The flow Reynolds number is $Re=60{,}000$ and the freestream Mach number is $M=0.1$. Numerical simulations are performed with a finite difference based solver that incorporates high-order compact schemes for differentiation, interpolation and filtering on a staggered grid. A mesh convergence study is conducted and results show good agreement with available data in terms of aerodynamic coefficients. Different spanwise arrangements of actuators are implemented to simulate blowing and suction at the airfoil leading edge. We observe that, for a specific frequency range of actuation, mean drag and drag fluctuations are substantially reduced while mean lift is maintained almost unaffected, especially for a 2D actuator setup. For this frequency range, 2D flow actuation disrupts the formation of the dynamic stall vortex, what leads to drag reduction due to a pressure increase along the airfoil suction side, towards the trailing edge region. At the same time, pressure is reduced on the suction side near the leading edge, increasing lift and further reducing drag.

We present a numerical methodology for construction of reduced order models, ROMs, of fluid flows through the combination of flow modal decomposition and regression analysis. Spectral proper orthogonal decomposition, SPOD, is applied to reduce the dimensionality of the model and, at the same time, filter the POD temporal modes. The regression step is performed by a deep feedforward neural network, DNN, and the current framework is implemented in a context similar to the sparse identification of non-linear dynamics algorithm, SINDy. A discussion on the optimization of the DNN hyperparameters is provided for obtaining the best ROMs and an assessment of these models is presented for a canonical nonlinear oscillator and the compressible flow past a cylinder. Then, the method is tested on the reconstruction of a turbulent flow computed by a large eddy simulation of a plunging airfoil under dynamic stall. The reduced order model is able to capture the dynamics of the leading edge stall vortex and the subsequent trailing edge vortex. For the cases analyzed, the numerical framework allows the prediction of the flowfield beyond the training window using larger time increments than those employed by the full order model. We also demonstrate the robustness of the current ROMs constructed via deep feedforward neural networks through a comparison with sparse regression. The DNN approach is able to learn transient features of the flow and presents more accurate and stable long-term predictions compared to sparse regression.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development., Comment: 57 pages, 57 figures; Accepted to ApJS

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction., Comment: 7 pages, 5 figures; Accepted for publication in ApJ Letters

After an optical peak, a classical or recurrent nova settles into a brief (days to years) period of quasi-stable thermonuclear burning in a compact configuration nearly at the white dwarf (WD) radius. During this time, the underlying WD becomes visible as a strong emitter of supersoft X-rays. Observations during this phase have revealed oscillations in the X-ray emission with periods on the order of tens of seconds. A proposed explanation for the source of these oscillations are internal gravity waves excited by nuclear reactions at the base of the hydrogen-burning layer. In this work, we present the first models exhibiting unstable surface $g$-modes with periods similar to oscillation periods found in galactic novae. However, when comparing mode periods of our models to the observed oscillations of several novae, we find that the modes which are excited have periods shorter than that observed., Comment: 13 pages, 7 figures. Accepted to ApJ

We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models from the pre main-sequence to core O-depletion at solar and subsolar metallicities for a total of 2000 Monte Carlo stellar models. For each stellar model, we independently and simultaneously sample 665 thermonuclear reaction rates and use them in a $\texttt{MESA}$ in situ reaction network that follows 127 isotopes from $^{1}$H to $^{64}$Zn. With this framework we survey the core mass, burning lifetime, composition, and structural properties at five different evolutionary epochs. At each epoch we measure the probability distribution function of the variations of each property and calculate Spearman Rank-Order Correlation coefficients for each sampled reaction rate to identify which reaction rate has the largest impact on the variations on each property. We find that uncertainties in $^{14}$N$(p,\gamma)^{15}$O, triple-$\alpha$, $^{12}$C$(\alpha,\gamma)^{16}$O, $^{12}$C($^{12}$C,$p$)$^{23}$Na, $^{12}$C($^{16}$O,$p$)$^{27}$Al, $^{16}$O($^{16}$O,$n$)$^{31}$S, $^{16}$O($^{16}$O,$p$)$^{31}$P, and $^{16}$O($^{16}$O,$\alpha$)$^{28}$Si reaction rates dominate the variations of the properties surveyed. We find that variations induced by uncertainties in nuclear reaction rates grow with each passing phase of evolution, and at core H-, He-depletion are of comparable magnitude to the variations induced by choices of mass resolution and network resolution. However, at core C-, Ne-, and O-depletion, the reaction rate uncertainties can dominate the variation causing uncertainty in various properties of the stellar model in the evolution towards iron core-collapse., Comment: Accepted for publication in ApJS, 33 pages, 21 figures