SEARCH · Search NASA

Results for “Hydrodynamic radius”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2021

The Asymmetric Inner Disk of the Herbig Ae Star HD 163296 in the Eyes of VLTI/MATISSE: Evidence for a Vortex?

Context.A complex environment exists in the inner few astronomical units of planet-forming disks. High-angular-resolution observa-tions play a key role in our understanding of the disk structure and the dynamical processes at work.Aims.In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-massstar HD 163296 from early VLTI/MATISSE observations taken in theL- andN-bands. We put special emphasis on the detection ofpotential disk asymmetries.Methods.We use simple geometric models to fit the interferometric visibilities and closure phases. Our models include a smoothedring, a flat disk with an inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetriesas well. We also perform numerical hydrodynamical simulations of the inner edge of the disk.Results.Our modeling reveals a significant brightness asymmetry in theL-band disk emission. The brightness maximum of the asym-metry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio inthe azimuthal variation is3.5±0.2. Comparing our result on the location of the asymmetry with other interferometric measurements,we confirm that the morphology of ther<0.3au disk region is time-variable. We propose that this asymmetric structure, located in ornear the inner rim of the dusty disk, orbits the star. To find the physical origin of the asymmetry, we tested a hypothesis where a vortexis created by Rossby wave instability, and we find that a unique large-scale vortex may be compatible with our data. The half-lightradius of theL-band-emitting region is0.33±0.01au, the inclination is52◦+5◦−7◦, and the position angle is143◦±3◦. Our models predictthat a non-negligible fraction of theL-band disk emission originates inside the dust sublimation radius forμm-sized grains. Refractorygrains or large (&10μm-sized) grains could be the origin of this emission.N-band observations may also support a lack of smallsilicate grains in the innermost disk (r.0.6au), in agreement with our findings fromL-band data.

J Varga↗

DOE OSTI2025

Dust supply to close binary systems

Binary systems can be born surrounded by circumbinary discs. The gaseous disc around either of the two stellar companions can have its life extended by the supply of mass arriving from the circumbinary disc. The objective of this study is to investigate the gravitational interactions exerted by a compact and eccentric binary system on the circumbinary and circumprimary discs, and the resulting transport of gas and solids between the disc components. We assume that the gas in the system behaves as a fluid, and we model its evolution by means of high-resolution hydrodynamical simulations. Dust grains are modelled as Lagrangian particles that interact with the gas and the stars. Our models indicate that significant fluxes of gas and dust proceed from the circumbinary disc towards the circumprimary disc. For the applied system parameters, grains of certain sizes are segregated outside the tidal gap generated by the stars. Consequently, the size distribution of the transported dust is not continuous, but presents a gap in the millimetre size range. In close binaries, the lifetime of an isolated circumprimary disc is found to be short, ~10 5 years, because of its low mass. However, because of the influx of gas from beyond the tidal gap, the disc around the primary star can survive much longer, ~10 6 years, as long as gas accretion from the circumbinary disc continues. The supply of solids and the extended lifetime of a circumbinary disc also aids in the possible formation of giant planets. Compared to close binary systems without a circumbinary disc, we expect a higher frequency of singleplanet or multiple-planet systems. Additionally, a planetesimal or debris belt can form in the proximity of the truncation radius of the circumprimary disc and/or around the location of the exterior edge of the tidal gap.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

The distribution of [α/Fe] in the Milky Way disc

ABSTRACT Using a sample of red giant stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we infer the conditional distribution $p([\alpha /{\rm Fe}]\, |\, [{\rm Fe}/{\rm H}])$ in the Milky Way disk for the α-elements Mg, O, Si, S, and Ca. In each bin of [Fe/H] and Galactocentric radius R, we model p([α/Fe]) as a sum of two Gaussians, representing ‘low-α’ and ‘high-α’ populations with scale heights $z_1=0.45\, {\rm kpc}$ and $z_2=0.95\, {\rm kpc}$, respectively. By accounting for age-dependent and z-dependent selection effects in APOGEE, we infer the [α/Fe] distributions that would be found for a fair sample of long-lived stars covering all z. Near the Solar circle, this distribution is bimodal at sub-solar [Fe/H], with the low-α and high-α peaks clearly separated by a minimum at intermediate [α/Fe]. In agreement with previous results, we find that the high-α population is more prominent at smaller R, lower [Fe/H], and larger |z|, and that the sequence separation is smaller for Si and Ca than for Mg, O, and S. We find significant intrinsic scatter in [α/Fe] at fixed [Fe/H] for both the low-α and high-α populations, typically ∼0.04-dex. The means, dispersions, and relative amplitudes of this two-Gaussian description, and the dependence of these parameters on R, [Fe/H], and α-element, provide a quantitative target for chemical evolution models and a test for hydrodynamic simulations of disk galaxy formation. We argue that explaining the observed bimodality will probably require one or more sharp transitions in the disk’s gas accretion, star formation, or outflow history in addition to radial mixing of stellar populations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Novel scalings of neutron star properties from analyzing dimensionless Tolman–Oppenheimer–Volkoff equations

The Tolman–Oppenheimer–Volkoff (TOV) equations govern the radial evolution of pressure and energy density in static neutron stars (NSs) in hydrodynamical equilibrium. Using the reduced pressure and energy density with respect to the NS central energy density, the original TOV equations can be recast into dimensionless forms. While the traditionally used integral approach for solving the original TOV equations require an input nuclear Equation of State (EOS), the dimensionless TOV equations can be anatomized by using the reduced pressure and energy density as polynomials of the reduced radial coordinate without using any input nuclear EOS. It has been shown in several of our recent works that interesting and novel perspectives about NS core EOS can be extracted directly from NS observables by using the latter approach. Our approach is based on intrinsic and perturbative analyses of the dimensionless (IPAD) TOV equations (IPAD-TOV). In this review article, we first discuss the length and energy density scales of NSs as well as the dimensionless TOV equations for scaled variables and their perturbative solutions near NS cores. We then review several new insights into NS physics gained from solving perturbatively the scaled TOV equations. Whenever appropriate, comparisons with the traditional approach from solving the original TOV equations will be made. In particular, we first show that the nonlinearity of the TOV equations basically excludes a linear EOS for dense matter in NS cores. We then show that perturbative analyses of the scaled TOV equations enable us to reveal novel scalings of the NS mass, radius and the compactness with certain combinations of the NS central pressure and energy density. Thus, observational data on either mass, radius or compactness can be used to constrain directly the core EOS of NS matter independent of the still very uncertain nuclear EOS models. As examples, the EOS of the densest visible matter in our Universe before the most massive neutron stars collapse into black holes (BHs) as well as the central EOS of a canonical or a 2.1 solar mass NS are extracted without using any nuclear EOS model. In addition, we show that causality in NSs sets an upper bound of about 0.374 for the ratio of pressure over energy density and correspondingly a lower limit for trace anomaly in supra-dense matter. We also demonstrate that the strong-field gravity plays a fundamental role in extruding a peak in the density/radius profile of the speed of sound squared (SSS) in massive NS cores independent of the nuclear EOS. Finally, some future perspectives of NS research using the new approach reviewed here by solving perturbatively the dimensionless TOV equations are outlined.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Assessment of dynamic-screw-pinch-driven, current-scaled MagLIF target implosion performance using 3D magnetohydrodynamic simulations

Analytic studies and two-dimensional “clean” radiation-magnetohydrodynamic (rad-MHD) simulations employing dynamical similarity driver-target scaling prescriptions [Ruiz et al., Phys. Plasmas 30, 032708 (2023)] suggest that Magnetized Liner Inertial Fusion (MagLIF) target implosions can scale to > 10 MJ DT fusion yields when peak drive current is increased beyond 60 MA. We present results from three-dimensional (3D) rad-MHD simulations of similarity-scaled MagLIF target implosions at peak drive currents ranging from 15 to 40 MA. Simulations in this study suggest that magneto-Rayleigh–Taylor instability (MRTI) growth and feedthrough to the fuel region are more severe at higher drive current scales, which reduces the fusion yield compared to prior analytic and 2D clean simulation predictions. In contrast to standard MagLIF, simulations of current-scaled MagLIF target implosions driven by a dynamic screw pinch (DSP) demonstrate reduced MRTI feedthrough and greater fuel magnetization, resulting in improved thermonuclear performance and enhanced performance scaling with peak drive current. DSP drive enables additional scaling of the liner mass to increase liner radius but maintain implosion time, resulting in higher implosion velocities at the expense of increased susceptibility to MRTI. We present a current- and mass-scaled simulated DSP-MagLIF target implosion at the ∼ 40 MA peak current level that produces ignition scale performance, demonstrating a burn-averaged Lawson ignition parameter above unity and DT fusion yield above 1 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Understanding ICF hohlraums using NIF gated laser-entrance-hole images

The newly available ns-gated laser-entrance-hole (LEH) imager on the National Ignition Facility provides routine, non-perturbative measurements of the x-ray emission from laser-heated plasmas inside the hohlraum as viewed at 19° to the hohlraum axis through one of its LEHs. Multiple images are acquired for a series of times and filter-selected x-ray energy bands within a single shot. The images provide time dependent data on phenomena including the effective radius of the LEH, the length of the gold-plasma “bubble” evolving off the interior wall surface heated by the outer beams, the evolving radius of the x-ray heated hohlraum wall, and the radius of the ablation front of the fusion capsule. Here, these measurements are explained and illustrated with sample data. These techniques are then applied to understand hohlraum behavior as a function of gas fill. For hohlraums with helium gas fill densities of 0.15 to 0.30 mg/cm 3 , synthetic images computed from simulations agree well with experimental gated LEH images when an inhibited heat transport model [Jones et al., Phys. Plasmas 24, 056312 (2017)] is used. This model can be adjusted to reproduce the expansion rate of the laser-heated plasma bubble in such a way as to improve agreement with the images. At the higher 0.6 mg/cc gas fill, the experimental images show more pronounced 3D features, resulting in slightly less good agreement with the 2D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗

NASA NTRS1994

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗

NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Boundary layer transition↗

NASA NTRS1999

Hydrodynamic Stability and Magnetic Reconnection in Disks and Stars

The purpose of this grant is to study parametric instability. The simplest example of parametric instability is a harmonic oscillator with a periodic modulation of the spring constant. If the modulation frequency is close to twice the natural frequency of the oscillator, the amplitude of oscillation tends to grow exponentially. The growth rate is proportional to the strength of the modulation, but it also depends upon the closeness to resonance of the two frequencies, and upon natural damping rate or "Q" of the oscillator. Parametric instabilities are very common in physics. A familiar example is a jogger's ponytail--normally a very strongly damped pendulum, it can be destabilized by the variation in effective gravity during the jogger's stride. Observation confirms that the period of the pendulum is half that of the jogger's vertical motion. In astrophysics, parametric instability may occur by external tidal forcing, or by interaction among eigenmodes. In the latter case, an energetic eigenmode may destabilize modes of half its frequency, provided some weak nonlinearity exists to couple them. Under a previous Astrophysical Theory grant (NAGW-2419), the PI discovered a parametric instability of tidally forced disks such as the accretion disks in cataclysmic variables and X ray binaries [2]. The destabilized modes are tightly-wound, incompressible, three-dimensional waves analogous to g-modes and r-modes in stars. Later work has confirmed our analysis [4]. It was hoped that these modes might provide a source of turbulence and angular momentum transport in accretion disks. However, a follow-up investigation of this instability by local numerical simulations, although confirming the analytically estimated growth rates, found negligible angular momentum flux [3]. Other work, partly supported by the ATP, now strongly indicates that the transport mechanism in such disks is magnetohydrodynamic turbulence [6]. Nevertheless, the parametric mechanism may truncate the outer edges of disks in close binaries [2], and it may be important in disks of very low ionization such as protostellar disks, or even cataclysmic-variable disks in quiescence where the MHD mechanism may be ineffective [5]. All analyses up to 1996 were done in a local approximation where the orbital frequency, shear rate, and tidal field were treated as constants. The locally computed growth rate turns out to depend strongly on radius, and it was unclear how to average these local rates to obtain the correct global rate. This is a critical issue for accretion disks in close binaries, because the local growth rate is comparable to the orbital frequency towards the outer edge of the disk but decreases rapidly inwards. Paper #1 examined this issue in a simplified global model where the destabilizing terms vary with position. We found that the global growth rate is essentially equal to the maximum local rate, provided that the latter is smoothed over a radial range equal to the distance that the destabilized wave propagates at its group speed in one growth time. Thus, in an accretion disk, waves would grow rapidly in the outer parts but would propagate both inwards and outwards at a maximum group speed of order the disk thickness divided by the orbital period.

Goodman, Jeremy↗

DOE OSTI2020

Gravitational-wave signal of a core-collapse supernova explosion of a 15M ⊙ star

In this paper we report on the gravitational-wave signal computed in the context of a three-dimensional simulation of a core-collapse supernova explosion of a 15 M ⊙ star. The simulation was performed with our neutrino hydrodynamics code chimera . We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (less than 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the standing accretion shock instability (SASI), and high-frequency emission (greater than 600 Hz) emanates from the proto–neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto–neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical–acoustic origin, to it. We compare the results of our three-dimensional model analysis with those obtained from the model’s two-dimensional counterpart and find a significant reduction in the strain amplitudes in the former case, as well as significant reductions in all related quantities. Our dominant proto–neutron star gravitational wave emission is not well characterized by emission from surface g modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto–neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. This simple analysis indicates that the spectrum of gravitational wave emission between approximately 20 Hz and approximately 1 kHz, stemming from neutrino-driven convection, the SASI, accretion onto the proto–neutron star, and proto–neutron star convection, will be accessible for a Galactic event.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1996

Multigrid Methods for EHL Problems

In many bearings and contacts, forces are transmitted through thin continuous fluid films which separate two contacting elements. Objects in contact are normally subjected to friction and wear which can be reduced effectively by using lubricants. If the lubricant film is sufficiently thin to prevent the opposing solids from coming into contact and carries the entire load, then we have hydrodynamic lubrication, where the lubricant film is determined by the motion and geometry of the solids. However, for loaded contacts of low geometrical conformity, such as gears, rolling contact bearings and cams, this is not the case due to high pressures and this is referred to as Elasto-Hydrodynamic Lubrication (EHL) In EHL, elastic deformation of the contacting elements and the increase in fluid viscosity with pressure are very significant and cannot be ignored. Since the deformation results in changing the geometry of the lubricating film, which in turn determines the pressure distribution, an EHL mathematical model must simultaneously satisfy the complex elasticity (integral) and the Reynolds lubrication (differential) equations. The nonlinear and coupled nature of the two equations makes numerical calculations computationally intensive. This is especially true for highly loaded problems found in practice. One novel feature of these problems is that the solution may exhibit sharp pressure spikes in the outlet region. To this date both finite element and finite difference methods have been used to solve EHL problems with perhaps greater emphasis on the use of the finite difference approach. In both cases, a major computational difficulty is ensuring convergence of the nonlinear equations solver to a steady state solution. Two successful methods for achieving this are direct iteration and multigrid methods. Direct iteration methods (e.g Gauss Seidel) have long been used in conjunction with finite difference discretizations on regular meshes. Perhaps one of the best examples of the application of such methods is the recent Effective Influence Method of Dowson and Wang. Multigrid methods have also been used with great success by Venner and Venner and Lubrecht with a good summary being given by Venner. As both these finite difference discretization based approaches appear to provide an efficient way of solving EHL problems, it is important to understand their relative merits. This paper is a first attempt at providing such an understanding in the context of EHL point contact problem, (contact of two spheres), in which the contact zone is a point and an ellipse or circle for unloaded and loaded dry contacts respectively. Since the film thickness and the contact width are generally small compared to the local radius of curvature of the two surfaces, the reduced geometry of the surfaces in the contact area can be accurately approximated to the contact between a paraboloid and a flat surface. The layout of the remainder of this paper is as follows. In section 2 we introduce the form of the equations to be solved. The Effective Influence Newton Method is described in Section 3 while Section 4 describes the Multigrid method to be used. Sections 5 and 6 describe the test problems to be used in the comparison between the two methods and compare the performance of the two methods. Section 7 concludes the paper with an argument of the two methods and suggests some future research directions.

Nurgat, Elyas↗

DOE OSTI2021

The antesonic condition for the explosion of core-collapse supernovae – II. Rotation and turbulence

In the problem of steady free fall on to a standing shockwave around a central mass, the ‘antesonic’ condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\le 3/16$, where cT is the isothermal sound speed in the subsonic post-shock flow, and vesc is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and the shock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models of accretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a ~5 percent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by ~20 to 40 percent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multidimensional simulations explode more easily than their one-dimensional counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2021

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

boundary layer transition↗

NASA NTRS2021

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Boundary Layer Transition↗