Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic efficiency”

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.

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

At least 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A global three-dimensional radiation magneto-hydrodynamic simulation of super-eddington accretion disks

We study super-Eddington accretion flows onto black holes using a global three-dimensional radiation magneto-hydrodynamical simulation. We solve the time-dependent radiative transfer equation for the specific intensities to accurately calculate the angular distribution of the emitted radiation. Turbulence generated by the magneto-rotational instability provides self-consistent angular momentum transfer. The simulation reaches inflow equilibrium with an accretion rate ∼220 L {sub Edd}/c {sup 2} and forms a radiation-driven outflow along the rotation axis. The mechanical energy flux carried by the outflow is ∼20% of the radiative energy flux. The total mass flux lost in the outflow is about 29% of the net accretion rate. The radiative luminosity of this flow is ∼10 L {sub Edd}. This yields a radiative efficiency ∼4.5%, which is comparable to the value in a standard thin disk model. In our simulation, vertical advection of radiation caused by magnetic buoyancy transports energy faster than photon diffusion, allowing a significant fraction of the photons to escape from the surface of the disk before being advected into the black hole. We contrast our results with the lower radiative efficiencies inferred in most models, such as the slim disk model, which neglect vertical advection. Our inferred radiative efficiencies also exceed published results from previous global numerical simulations, which did not attribute a significant role to vertical advection. We briefly discuss the implications for the growth of supermassive black holes in the early universe and describe how these results provided a basis for explaining the spectrum and population statistics of ultraluminous X-ray sources.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

COOP 3D ARPA Experiment 109 National Center for Atmospheric Research

Coupled atmospheric and hydrodynamic forecast models were executed on the supercomputing resources of the National Center for Atmospheric Research (NCAR) in Boulder, Colorado and the Ohio Supercomputing Center (OSC)in Columbus, Ohio. respectively. The interoperation of the forecast models on these geographically diverse, high performance Cray platforms required the transfer of large three dimensional data sets at very high information rates. High capacity, terrestrial fiber optic transmission system technologies were integrated with those of an experimental high speed communications satellite in Geosynchronous Earth Orbit (GEO) to test the integration of the two systems. Operation over a spacecraft in GEO orbit required modification of the standard configuration of legacy data communications protocols to facilitate their ability to perform efficiently in the changing environment characteristic of a hybrid network. The success of this performance tuning enabled the use of such an architecture to facilitate high data rate, fiber optic quality data communications between high performance systems not accessible to standard terrestrial fiber transmission systems. Thus obviating the performance degradation often found in contemporary earth/satellite hybrids.

Source record↗
DOE OSTI2019

Performance Analysis and Tank Test Validation of a Hybrid Wave-Current Energy Converter with a Single Power Takeoff

Marine and hydrokinetic (MHK) energy, including ocean waves, tidal current, ocean current and river current, has been recognized as a promising power source due to its full-day availability and high energy potential. At this stage, ocean current energy, tidal energy and ocean wave energy are currently the most competitive sourves among all the categories of MHK. The state of art MHK energy harvesting technology mainly focus on harvesting either ocean wave energy or current energy, but not both. However, a significant amount of ocean waves and tidal/ ocean current coexist in many sites and traditional devices that harvest from a single form of MHK energy, cannot make full use of the coexisting ocean energy. Furthermore, MHK energy harvesting devices need to advance to be cost-effective and competitive with other energy sources. This is difficult to achieve. Ocean wave excitation is irregular, which means that ocean wave height and wave periods are unpredictable and excitation forces on energy harvesting devices can have large variance in amplitude and frequency. Such problems/ restrictions can be possibly addressed by the concept of a hybrid energy converter. In this sense, a hybrid wave-current ocean energy conveter (HWCEC) that simutaneously harvests energy from current and wave with one single power takeoff (PTO) is designed.The wave energy is extracted through relative heaving motion between a floating buoy and a submerged second body, while the current energy is extracted using a marine current turbine (MCT). Energy from both sources are integrated by a hybrid PTO whose concept is based on a mechanical motion rectifier (MMR). In this study, different working modes are investigated together with switching criteria.Simulations were conducted with hydrodynamic coefficients obtained from computational fluid dynamics analysis and boundary element method. Tank tests were conducted for a HWCEC under coexisting wave and current inputs. For comparison, separate baseline tests of a turbine and a two-body point absorber, each acting in isolation, are conducted. Experimental results validate the dynamic modeling and show that a HWCEC can increase the output power with a range between 29-87 percent over either current turbine and wave energy converter acting individually, and it can reduce by up to 70 percent the peak-to-average power ratio compared with the wave energy converter on the tested conditions.Such results demonstrate the potential of the HWCEC as an efficient and cost-effective design.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Functionalized Magnetic Nanoparticles for Technetium Sequestration from Groundwater

Technetium 99 (Tc) is among the most common environmental contaminants at DOE sites and one of the most common risk drivers in low- and high-level waste disposal sites. The majority of Tc is generated from anthropogenic sources, such as nuclear power plants, global weapons, nuclear storage facilities and medical applications. Through these sources, Tc contamination has been unintentionally introduced in to the environment. The most common chemical form of Tc is Tc(VII)O{sub 4}{sup -}. Due to its high solubility and mobility, Tc can enter the food chain and cause adverse health effects to humans. Currently, ion exchange resins and reduction processes are the most common approaches for Tc immobilization. Although these techniques have shown to be effective, they also possess major drawbacks, such as high cost, low adsorption capacity, and complex creation and maintenance. Therefore, development of more efficient and simple technologies for the remediation of Tc-contaminated systems are needed. Functionalized magnetic nanoparticles have been used to remove organic and inorganic contaminants from water resources. These nanoparticles have attracted extensive attention as an adsorbent material due to their large surface area, high efficiency, low-cost, easy functionalization and separation with a magnet. This study seek to develop functionalized magnetic iron oxide nanoparticles for the efficient removal of Tc and other heavy metal contaminants from water resources under ambient conditions. Objectives: Synthesize magnetic iron oxide nanoparticles and functionalize their surface with Cetyltrimethylammonium Bromide (CTAB) and tetraethyl-orthosilicate (TEOS). Characterize the synthesized nanoparticles using scanning electron microscopy (SEM) coupled with energy dispersive X-ray spectroscopy (EDS), Dynamic Light Scattering (DLS) and Zeta PALS. Perform adsorption studies to evaluate their adsorption behavior and capacity for (a) Technetium using Rhenium (ReO{sub 4}{sup -}) as a surrogate and (b) heavy metals, e.g. Cu{sup 2+}. Conclusions: Magnetic iron oxide nanoparticles were successfully functionalized with CTAB and TEOS. The functionalization of the iron oxide nanoparticles affects their surface charge and their hydrodynamic diameter. The addition of CTAB or TEOS decreased the hydrodynamic diameter of the nanoparticles due to repulsive and steric forces. The SEM micrographs show spherical nanoparticles of different sizes. The EDX analysis shows the presence of iron and oxygen from the iron oxide crystalline structure, and the different constituents of the CTAB and TEOS molecules. Proof-of-concept shows the successful adsorption of rhenium (ReO{sub 4}{sup -}) and copper Cu{sup 2+}) onto CTAB-Fe{sub 2}O{sub 3} and TEOS-Fe{sub 2}O{sub 3} respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2012

Magnetic Field Effects on Plasma Plumes

Here, we will discuss our numerical studies of plasma jets and loops, of basic interest for plasma propulsion and plasma astrophysics. Space plasma propulsion systems require strong guiding magnetic fields known as magnetic nozzles to control plasma flow and produce thrust. Propulsion methods currently being developed that require magnetic nozzles include the VAriable Specific Impulse Magnetoplasma Rocket (VASIMR) [1] and magnetoplasmadynamic thrusters. Magnetic nozzles are functionally similar to de Laval nozzles, but are inherently more complex due to electromagnetic field interactions. The two crucial physical phenomenon are thrust production and plasma detachment. Thrust production encompasses the energy conversion within the nozzle and momentum transfer to a spacecraft. Plasma detachment through magnetic reconnection addresses the problem of the fluid separating efficiently from the magnetic field lines to produce maximum thrust. Plasma jets similar to those of VASIMR will be studied with particular interest in dual jet configurations, which begin as a plasma loops between two nozzles. This research strives to fulfill a need for computational study of these systems and should culminate with a greater understanding of the crucial physics of magnetic nozzles with dual jet plasma thrusters, as well as astrophysics problems such as magnetic reconnection and dynamics of coronal loops.[2] To study this problem a novel, hybrid kinetic theory and single fluid magnetohydrodynamic (MHD) solver known as the Magneto-Gas Kinetic Method is used.[3] The solver is comprised of a "hydrodynamic" portion based on the Gas Kinetic Method and a "magnetic" portion that accounts for the electromagnetic behaviour of the fluid through source terms based on the resistive MHD equations. This method is being further developed to include additional physics such as the Hall effect. Here, we will discuss the current level of code development, as well as numerical simulation results

Ebersohn, F.↗
DOE OSTI2013

THE BIGGEST EXPLOSIONS IN THE UNIVERSE. II

One of the leading contenders for the origin of supermassive black holes (SMBHs) at z ∼> 7 is catastrophic baryon collapse in atomically cooled halos at z ∼ 15. In this scenario, a few protogalaxies form in the presence of strong Lyman-Werner UV backgrounds that quench H{sub 2} formation in their constituent halos, preventing them from forming stars or blowing heavy elements into the intergalactic medium prior to formation. At masses of 10{sup 8} M{sub ☉} and virial temperatures of 10{sup 4} K, gas in these halos rapidly cools by H lines, in some cases forming 10{sup 4}-10{sup 6} M{sub ☉} Population III stars and, a short time later, the seeds of SMBHs. Instead of collapsing directly to black holes (BHs), some of these stars died in the most energetic thermonuclear explosions in the universe. We have modeled the explosions of such stars in the dense cores of line-cooled protogalaxies in the presence of cosmological flows. In stark contrast to the explosions in diffuse regions in previous simulations, these supernovae briefly engulf the protogalaxy, but then collapse back into its dark matter potential. Fallback drives turbulence that efficiently distributes metals throughout the interior of the halo and fuels the rapid growth of nascent BHs at its center. The accompanying starburst and X-ray emission from these line-cooled galaxies easily distinguish them from more slowly evolving neighbors and might reveal the birthplaces of SMBHs on the sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Multiscale method based on coupled lattice‐Boltzmann and Langevin‐dynamics for direct simulation of nanoscale particle/polymer suspensions in complex flows

Summary A hybrid computational method coupling the lattice‐Boltzmann (LB) method and a Langevin‐dynamics (LD) method is developed to simulate nanoscale particle and polymer (NPP) suspensions in the presence of both thermal fluctuation and long‐range many‐body hydrodynamic interactions (HIs). Brownian motion of the NPP is explicitly captured by a stochastic forcing term in the LD method. The LD method is two‐way coupled to the nonfluctuating LB fluid through a discrete LB forcing source distribution to capture the long‐range HI. To ensure intrinsically linear scalability with respect to the number of particles, a Eulerian‐host algorithm for short‐distance particle neighbor search and interaction is developed and embedded to LB‐LD framework. The validity and accuracy of the LB‐LD approach are demonstrated through several sample problems. The simulation results show good agreements with theory and experiment. The LB‐LD approach can be favorably incorporated into complex multiscale computational frameworks for efficiently simulating multiscale multicomponent particulate suspension systems such as complex blood suspensions.

Liu, Zixiang↗
DOE OSTI2022

Collision of localized shocks in AdS 5 as a series expansion in transverse gradients

We introduce a computational framework to more efficiently calculate the collision of localized shocks in five dimensional asymptotically Anti-de Sitter space. We expand the Einstein equations in transverse gradients and find that our numerical results agree well with exact solutions already at first order in the expansion. Moreover, the Einstein equations at first order in transverse gradients can be decoupled into two sets of differential equations. Here, the bulk fields of one of these sets has only a negligible contribution to boundary observables, such that the computation on each time slice can be simplified to the solution of several planar shockwave equations plus four further differential equations for each transverse plane ‘pixel’. At the cost of errors of ≲ 10% at the hydrodynamization time and for low to mid rapidities, useful numerical solutions can be sped up by roughly one order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Hohlraum x-ray preheat asymmetry measurement at the ICF capsule via Mo ball fluorescence imaging

In inertial confinement fusion, penetrating asymmetric hohlraum preheat radiation (>1.8 keV, which includes high temperature coronal M-band emission from laser spots) can lead to asymmetric ablation front and ablator–fuel interface hydrodynamic instability growth in the imploding capsule. First experiments to infer the preheat asymmetries at the capsule were performed on the National Ignition Facility for high density carbon (HDC) capsules in low density fill (0.3 mg/cc 4He) Au hohlraums by time resolved imaging of 2.3 keV fluorescence emission of a smaller Mo sphere placed inside the capsule. Measured Mo emission is pole hot (P2 > 0) since M-band is generated mainly by the outer laser beams as their irradiance at the hohlraum wall is 5× higher than for the inner beams. P2 has a large swing vs time, giving insight into the laser heated hohlraum dynamics. P4 asymmetry is small at the sphere due to efficient geometric smoothing of hohlraum P4 asymmetries at large hohlraum-to-capsule radii ratios. The asymmetry at the HDC capsule is inferred from the Mo emission asymmetry accounting for the Mo/HDC radius difference and HDC capsule opacity.

Dewald, E. L. (ORCID:0000000340546967)↗
DOE OSTI2015

A physical approach to the identification of high-Z mergers: Morphological classification in the stellar mass domain

At z ≳ 1, the distinction between merging and “normal” star-forming galaxies based on single band morphology is often hampered by the presence of large clumps which result in a disturbed, merger-like appearance even in rotationally supported disks. In this paper we discuss how a classification based on canonical, non-parametric structural indices measured on resolved stellar mass maps, rather than on single-band images, reduces the misclassification of clumpy but not merging galaxies. We calibrate the mass-based selection of mergers using the MIRAGE hydrodynamical numerical simulations of isolated and merging galaxies which span a stellar mass range of 10{sup 9.8}–10{sup 10.6} M{sub ⊙} and merger ratios between 1:1–1:6.3. These simulations are processed to reproduce the typical depth and spatial resolution of observed Hubble Ultra Deep Field (HUDF) data. We test our approach on a sample of real z≃2 galaxies with kinematic classification into disks or mergers and on ∼100 galaxies in the HUDF field with photometric/spectroscopic redshift between 1.5 ≤ z ≤ 3 and M > 10{sup 9.4} M{sub ⊙}. We find that a combination of the asymmetry A{sub MASS} and M{sub 20, MASS} indices measured on the stellar mass maps can efficiently identify real (major) mergers with ≲20% contamination from clumpy disks in the merger sample. This mass-based classification cannot be reproduced in star-forming galaxies by H-band measurements alone, which instead result in a contamination from clumpy galaxies which can be as high as 50%. Moreover, we find that the mass-based classification always results in a lower contamination from clumpy galaxies than an H-band classification, regardless of the depth of the imaging used (e.g., CANDELS versus HUDF).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Exploring Flood Predictability in Taiwan through Coupled Atmospheric–Hydrological and High-Performance Hydrodynamic Models

Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.

Coupled models↗
DOE OSTI2019

Harnessing billions of tasks for a scalable portable hydrodynamic simulation of the merger of two stars

We propose a highly scalable demonstration of a portable asynchronous many-task programming model and runtime system applied to a grid-based adaptive mesh refinement hydrodynamic simulation of a double white dwarf merger with 14 levels of refinement that spans 17 orders of magnitude in astrophysical densities. The code uses the portable C++ parallel programming model that is embodied in the HPX library and being incorporated into the ISO C++ standard. The model reflects a significant shift from existing bulk synchronous parallel programming models under consideration for exascale systems. Through the use of the Futurization technique, seemingly sequential code is transformed into wait-free asynchronous tasks. We demonstrate the potential of our model by showing results from strong scaling runs on National Energy Research Scientific Computing Center’s Cori system (658,784 Intel Knight’s Landing cores) that achieve a parallel efficiency of 96.8% using billions of asynchronous tasks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

ATMOSPHERIC CIRCULATION OF HOT JUPITERS: DAYSIDE–NIGHTSIDE TEMPERATURE DIFFERENCES

The full-phase infrared light curves of low-eccentricity hot Jupiters show a trend of increasing dayside-to-nightside brightness temperature difference with increasing equilibrium temperature. Here, we present a three-dimensional model that explains this relationship, in order to provide insight into the processes that control heat redistribution in tidally locked planetary atmospheres. This three-dimensional model combines predictive analytic theory for the atmospheric circulation and dayside–nightside temperature differences over a range of equilibrium temperatures, atmospheric compositions, and potential frictional drag strengths with numerical solutions of the circulation that verify this analytic theory. The theory shows that the longitudinal propagation of waves mediates dayside–nightside temperature differences in hot Jupiter atmospheres, analogous to the wave adjustment mechanism that regulates the thermal structure in Earth’s tropics. These waves can be damped in hot Jupiter atmospheres by either radiative cooling or potential frictional drag. This frictional drag would likely be caused by Lorentz forces in a partially ionized atmosphere threaded by a background magnetic field, and would increase in strength with increasing temperature. Additionally, the amplitude of radiative heating and cooling increases with increasing temperature, and hence both radiative heating/cooling and frictional drag damp waves more efficiently with increasing equilibrium temperature. Radiative heating and cooling play the largest role in controlling dayside–nightside temperature differences in both our analytic theory and numerical simulations, with frictional drag only being important if it is stronger than the Coriolis force. As a result, dayside–nightside temperature differences in hot Jupiter atmospheres increase with increasing stellar irradiation and decrease with increasing pressure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Applications of tribology to determine attrition by wear of particulate solids in CFB systems

In recent years, much attention has been focused on the development of novel technologies for carbon capture and chemicals production that utilize a circulating fluidized bed (CFB) configuration; examples include chemical looping combustion and circulation of temperature swing adsorbents in a CFB configuration for CO 2 capture. A major uncertainty in determining the economic feasibility of these technologies is the required solids makeup rate, which, among other factors, is due to impact and wear attrition at various locations, including standpipes, cyclones, and the gas jets in fluid beds. While correlations have been developed that estimate the attrition rates at these areas, these correlations are dependent on constants that are uncertain without extensive experiment in the corresponding unit operation. Thus, it is difficult to determine the attrition rate a priori without performing extensive experiments on the materials or scaling up entirely. In this work, the authors outline a methodology for predictive attrition based on fundamental material properties from fields of tribology—specifically, the study of wear—to the knowledge of forces and sliding distances determined from hydrodynamic models to develop basic attrition models for novel CFB systems. The equations are derived for the standpipe and cyclone, which are common components found in CFBs, and the cyclone equation is compared to experimental data of attrition in the literature. The cyclone equation derived in this work results in an abrasion rate based on (1) material properties such as particle density and hardness, (2) inlet velocity, and (3) cyclone geometry. According to this equation, increasing the diameter of the cyclone and the solids inlet velocity tends to increase the rate of abrasion of the catalyst, while decreasing the hardness increases the abrasion rate. The functionality of the increasing attrition rate with velocity increase implies that increasing the efficiency of the cyclone may also increase the attrition rate via abrasion. With modifications to the severity coefficient term to include the solids loading, the cyclone equation derived in this work fits data from Reppenhagen and Werther with a coefficient of determination (R2) of 92%.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

A comprehensive review of nonlinear oscillators in hydrokinetic energy harnessing using flow-induced vibrations

A comprehensive review of hydrokinetic energy converters based on alternating lift technology (ALT) is provided. Emphasis is on nonlinear oscillators based on Flow Induced Vibration (FIV) or Oscillation (FIO). Due to strong coupling in Fluid-Structure Interaction (FSI), and in order to maximize the hydrokinetic harnessed energy, design of nonlinear oscillators and analysis by model tests or computational fluid dynamics dominates this area. Research confirmed that the nonlinear oscillator can harvest energy from a stochastic excitation modeled by a generic wide spectrum, and overcome the most severe oscillator limitations: specifically, the need for continuous frequency tuning due to the narrow bandwidth response, and low efficiency outside the narrow bandwidth oscillator response. This review covers the following aspects of nonlinear oscillators in ALT converters: (1) Geometric changes in oscillator cross-section; e.g., circular, square, rectangular, or trilateral shapes. (2) Passive turbulence control of FIV/FIO. (3) Position based nonlinear stiffness. (4) Multi-cylinder synergistic FIV/FIO. (5) Mechanically linked oscillators. (6) Velocity-based, nonlinear, adaptive harnessing damping.

42 ENGINEERING↗
DOE OSTI2015

Triggering collapse of the presolar dense cloud core and injecting short-lived radioisotopes with a shock wave. IV. effects of rotational axis orientation

Both astronomical observations of the interaction of Type II supernova remnants (SNRs) with dense interstellar clouds as well as cosmochemical studies of the abundances of daughter products of short-lived radioisotopes (SLRIs) formed by supernova nucleosynthesis support the hypothesis that the Solar System's SLRIs may have been derived from a supernova. This paper continues a series devoted to examining whether or not such a shock wave could have triggered the dynamical collapse of a dense, presolar cloud core and simultaneously injected sufficient abundances of SLRIs to explain the cosmochemical evidence. Here, we examine the effects of shock waves striking clouds whose spin axes are oriented perpendicular, rather than parallel, to the direction of propagation of the shock front. The models start with 2.2 M{sub ⊙} cloud cores and shock speeds of 20 or 40 km s{sup −1}. Central protostars and protoplanetary disks form in all models, although with their disk spin axes aligned somewhat randomly. The disks derive most of their angular momentum not from the initial cloud rotation, but from the Rayleigh-Taylor fingers that also inject shock wave SLRIs. Injection efficiencies, f{sub i}, the fraction of the incident shock wave material injected into the collapsing cloud core, are ∼0.04–0.1 in these models, similar to when the rotation axis is parallel to the shock propagation direction. Evidently, altering the rotation axis orientation has only a minor effect on the outcome, strengthening the case for this scenario as an explanation for the Solar System's SLRIs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Active galactic nucleus feedback in an isolated elliptical galaxy: The effect of strong radiative feedback in the kinetic mode

Based on two-dimensional high-resolution hydrodynamic numerical simulation, we study the mechanical and radiative feedback effects from the central active galactic nucleus (AGN) on the cosmological evolution of an isolated elliptical galaxy. The inner boundary of the simulation domain is carefully chosen so that the fiducial Bondi radius is resolved and the accretion rate of the black hole is determined self-consistently. It is well known that when the accretion rates are high and low, the central AGNs will be in cold and hot accretion modes, which correspond to the radiative and kinetic feedback modes, respectively. The emitted spectrum from the hot accretion flows is harder than that from the cold accretion flows, which could result in a higher Compton temperature accompanied by a more efficient radiative heating, according to previous theoretical works. Such a difference of the Compton temperature between the two feedback modes, the focus of this study, has been neglected in previous works. Significant differences in the kinetic feedback mode are found as a result of the stronger Compton heating. More importantly, if we constrain models to correctly predict black hole growth and AGN duty cycle after cosmological evolution, we find that the favored model parameters are constrained: mechanical feedback efficiency diminishes with decreasing luminosity (the maximum efficiency being ≅ 10{sup –3.5}), and X-ray Compton temperature increases with decreasing luminosity, although models with fixed mechanical efficiency and Compton temperature can be found that are satisfactory as well. We conclude that radiative feedback in the kinetic mode is much more important than previously thought.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗