Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Compact model”

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 271 records · Page 15

Comparison of Fission Product Release Predictions using PARFUME and BISON with Results from the AGR-3/4 Irradiation Experiment

The PARFUME (PARticle Fuel ModEl) fuel performance modeling code and the BISON nuclear fuel performance application built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library were used to predict the fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the third and fourth irradiation experiment of the Advanced Gas Reactor (AGR-3/4) Fuel Development and Qualification Program. The fuel performance modeling codes PARFUME and BISON modeled the AGR-3/4 irradiation experiment using the fuel compact time-averaged volume averaged (TAVA) daily temperatures for a total irradiation duration of 369.1 effective full power days (EFPD) to predict the release fraction of the fission product silver (Ag-110m) from a representative TRISO-coated fuel particle from AGR-3/4 compacts. Post-irradiation examination (PIE) measurements provided data on the release of these fission products in the compacts outside of the silicon carbide (SIC) layer. The PARFUME and BISON results were then compared to the silver release measured from compact gamma scanning. The results showed good agreement between PARFUME and BISON but both codes under-predicted the silver release fraction for all the compacts. In addition, BISON was used to model and predict the fission product concentration radial profile outside of the compacts in capsules’ inner and outer rings. These rings were either comprised of matrix and/or structural graphite. To obtain the concentration profiles of silver, cesium, and strontium, a sorption isotherm model was developed in BISON to capture the effects of fission product transport across the gaps between the concentric rings. The general shape of the concentration radial profiles as calculated by BISON were similar in the inner ring (IR) but varied in the outer ring (OR) depending on the fission product of interest or capsule temperature. Using this methodology and model, BISON now has the capability to aid in developing new fission product diffusion coefficients for matrix or structural graphite materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multiscale Porous High-temperature Heat Exchanger Using Ceramic Co-extrusion

In this project, our MIT, Purdue, and GE team aims to design, model, fabricate, and test a novel high temperature, compact, and durable ceramic heat exchanger to be operated under high temperature and pressure conditions for aerospace applications. Our approach is grounded in introducing multiscale porosity, i.e., centimeter-scale channels embedded with micrometer-scale channels, into the ceramic heat exchanger to significantly improve its heat transfer performance and mechanical strength while maintaining minimal pressure losses. We first developed high-fidelity thermal-fluid-mechanical model capable of precisely capturing the heat transfer rate, temperature profile, pressure drop, and mechanical stress throughout the entire heat exchanger design. Guided by our model, we identified the optimal design parameters for the SiC heat exchanger body and manifolds. Then, we established a completed fabrication procedure to create multiscale features in the ceramic heat exchanger, including co-extrusion, lamination, burnout, and sintering. We fabricated multiple heat exchanger bodies consisting of 6 × 6 and 3 × 3 centimeter-scale channels where each individual centimeter-scale channel comprises 625 crack-free microchannels with 90 μm × 90 μm opening. Owing to the multiscale features, our fabricated heat exchanger bodies exhibited desirable mechanical strength with 156 MPa flexural strength under 1300 Celsius degree. Despite the demonstrated highly tailorable microscopic features and superior mechanical strength, we identified delamination due to the complex interaction among ceramic, polymer, and gas species can be a critical challenge to create fully defect-free heat exchanger, which requires further fundamental investigations in future works. To test the heat exchanger performance, we constructed a high-temperature and high-pressure experimental apparatus that can be safely operated under 400 Celsius degree and 4 bar. With insights gained from mechanistic modeling, material development, and detailed characterization, a cost model was finally developed to understand the market potential of the developed technology, where a cost of $43,000 Celsius degree/kW was envisioned. This project developed a transformative approach to high-performance heat exchanger design. The thermal-fluid-mechanical design approach developed in this project can serve as a generic tool to guide the design of various heat-exchangers operated under high-temperature and high-pressure conditions. The material fabrication approach established in this project can be a useful guide for ceramic processing at extreme length scales.

42 ENGINEERING↗

HTS Cable Conductor for Compact Fusion Tokamak Solenoids

Significant progress has been made recently in the US fusion community to develop a strategic plan to enable engineering design and construction of a fusion pilot plant (FPP). Princeton Plasma Physics Laboratory (PPPL) is working on developing high current density HTS conductors for next fusion experiments. Partnering with the US industry, we are evaluating feasibility and affordability of Conductor on Round Core (CORC ® ) developed by Advanced Conductor Technologies (ACT) for the next compact fusion tokamak facility. High current density achieved by a CORC ® cable based on a four-layer model coil recently tested at National High Magnetic Field Laboratory (NHMFL) motivated its consideration for low cost, reduced size fusion magnet application. This is of interest to PPPL because of its scalability to tokamak central solenoid (CS) coils in terms of required flux swings for plasma startup operations. Furthermore, partnering with ACT, we designed and built a two-layer model coil solenoid directly wound with CORC to demonstrate its applicability for compact solenoids such as that used in the national spherical torus experiment (NSTX), NSTX-upgrade (NSTX-U) and the US sustained high-power density test facility (SHPD). The ~160 mm diameter solenoid wound by a two-layer CORC is being tested in early 2022 under electromagnetic cyclic loading at NHMFL in a unique 160 mm bore, 14 T background field magnet facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling gas breakthrough and flow phenomena through engineered barrier systems using a discrete fracture approach

Compacted clays are being considered to build engineered barrier systems (EBS) intended for the safe isolation of high-level nuclear waste (HLW) and spent nuclear fuel (SNF). The corrosion of the metallic canister containing the HLW/SNF will lead to the generation and buildup of the gas pressure in the more internal part of the clay buffer. This phenomenon would eventually trigger the formation and propagation of fractures in the clay barrier, jeopardizing its safety functions. Here, in this work, we propose to use the fragmentation technique (MFT) to model evolving fractures in clays triggered by gas pressurization. The MFT has been successfully used to model the formation of fractures in concrete, drying cracks in soil, hydraulic and thermo-fractures in rocks. In this work, we extend the MFT to deal with multiphase fluid flow in deformable porous media, and we upgraded a fully coupled computer finite element code using the extended technique. The proposed approach is first verified against analytical solutions and is then applied to model gas breakthrough experiments in clays. A very satisfactory performance of the method is observed in all the analyses, showing the potential of the MFT to tackle multiphase flow problems in deformable porous media with evolving discontinuities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Aspects of the Domain-Wall Standard Model

Abstract In the “Domain-Wall Standard Model,” all the Standard Model (SM) fields are localized in certain domains in a non-compact 5D space-time. While the SM is realized as a 4D effective theory at low energies, the model involves Kaluza–Klein (KK) modes of the SM particles. In this paper we introduce two simple solvable examples which lead to domain-wall configurations for the SM particles and their KK-modes. Based on the examples, we address a variety of phenomenologies of the Domain-Wall SM, such as the KK-mode gauge boson phenomenology at the Large Hadron Collider (LHC), the effect of the KK-mode SM fermions on Higgs boson phenomenology, and the KK-mode fermion search at the LHC with its decay into a corresponding SM fermion and a Nambu–Goldstone boson associated with a spontaneous breaking of the translational invariance in the 5th dimension. We also propose a simple unified picture of localizing all the SM fields.

Physics↗

Additive Manufactured Compact Microwave Absorbers

A high-performance, compact microwave absorber was created using Fused Deposition Modeling (FDM) 3D printing. Both a narrowband and a broadband absorber were created. The narrowband absorber was designed at 4.9 GHz, mid-band in WR-187 waveguide. The broadband absorber tried to achieve the best attenuation across the entire 3.95 to 5.85 GHz band. Two types of carbon loaded polylactic acid (PLA) plastic and one type of unloaded PLA were 3D printed with variable percentages of air to achieve different values of effective dielectric constant and loss tangent. The absorber comprised five or six rectangular pieces of these plastic materials. The thickness and fill factor values for each piece were optimized to minimize reflection through fast analytic modeling in MATLAB®. The results were then verified by HFSS® simulation as well. The stack progressed from the lowest loss and lowest dielectric constant to the highest at the shorting end. The final narrowband load had simulated return loss of 87 dB at 4.9 GHz with an analytic solution in MATLAB. The measured return loss of the 3D printed attenuator was 73 dB at 4.929 GHz. The total length of the absorber was 2.44 inches. A commercial absorber for WR-187 with return loss of 40 dB has length of 13 inches. The experiment proves that an effective and compact microwave absorber can be created using 3D printing.

36 MATERIALS SCIENCE↗

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗

Updating and Refining BISON TRISO Assessment

Several advanced reactor concepts consider tristructural isotropic (TRISO) fuel in their designs. To support the advanced reactor industry, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program within the U.S. Department of Energy (DOE) is developing modeling and simulation capabilities within the BISON fuel performance code. Assessment of these capabilities against available experi- mental data, such as the Advanced Gas Reactor (AGR) program, is essential in instilling confidence in the code predictions. Associated uncertainties in model inadequacy, experimental measurements, and manufacturing are also of importance. This report discusses the updating and standardization of the existing assessment suite in BISON for TRISO fuel performance analysis. The cases considered the AGR-1, AGR-2, and AGR-3/4 experimental programs. The report also includes a preliminary investigation into the sources of model inadequacy associated with silver (Ag) release fractions using the AGR-2 Compact 2-1-2 as a case study. The results indicate that uncertainties in the operating temperatures supplied to the model and the form of those temperature conditions (spatially averaged versus spatially resolved) both affect Ag release predictions. This suggests that previously observed model inadequacy can be attributed to both physical models and the operating conditions supplied to them. These findings will be leveraged to guide development of assessment cases for additional experiments and refinement of physical models for behaviors that potentially influence Ag release, such as palladium penetration.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy

A significant number of hypervelocity stars with velocities between 1500 − 2500 km s −1 have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new AREPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M ⊙ remains bound, moving at a velocity of ∼2800 km s −1 . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.

Astronomy and AstroPhysics↗

SEC ‐ SAXS / MC Ensemble Structural Studies of the Microtubule Binding Protein Cdt1 Show Monomeric, Folded‐Over Conformations

ABSTRACT Cdt1 is a mixed folded protein critical for DNA replication licensing and it also has a “moonlighting” role at the kinetochore via direct binding to microtubules and the Ndc80 complex. However, it is unknown how the structure and conformations of Cdt1 could allow it to participate in these multiple, unique sets of protein complexes. While robust methods exist to study entirely folded or unfolded proteins, structure–function studies of combined, mixed folded/disordered proteins remain challenging. In this work, we employ orthogonal biophysical and computational techniques to provide structural characterization of mitosis‐competent human Cdt1. Thermal stability analyses shows that both folded winged helix domains1 are unstable. CD and NMR show that the N‐terminal and linker regions are intrinsically disordered. DLS shows that Cdt1 is monomeric and polydisperse, while SEC‐MALS confirms that it is monomeric at high concentrations, but without any apparent inter‐molecular self‐association. SEC‐SAXS enabled computational modeling of the protein structures. Using the program SASSIE, we performed rigid body Monte Carlo simulations to generate a conformational ensemble of structures. We observe that neither fully extended nor extremely compact Cdt1 conformations are consistent with SAXS. The best‐fit models have the N‐terminal and linker disordered regions extended into the solution and the two folded domains close to each other in apparent “folded over” conformations. We hypothesize the best‐fit Cdt1 conformations could be consistent with a function as a scaffold protein that may be sterically blocked without binding partners. Our study also provides a template for combining experimental and computational techniques to study mixed‐folded proteins.

Cell Biology↗

Hybrid Star Models in the Light of New Multimessenger Data

Abstract Recent astrophysical mass inferences of compact stars HESS J1731-347 and PSR J0952-0607, with extremely small and large masses respectively, as well as the measurement of the neutron skin of Ca in the CREX experiment challenge and constrain the models of dense matter. We examine the concept of hybrid stars—objects containing quark cores surrounded by nucleonic envelopes—as models that account for these new data along with other inferences. We employ a family of 81 nucleonic equations of state (EOSs) with variable skewness and slope of symmetry energy at saturation density and a constant speed-of-sound EOS for quark matter. For each nucleonic EOS, a family of hybrid EOSs is generated by varying the transition density, the energy jump, and the speed of sound. These models are tested against the data from GW170817 and J1731-347, which favor low-density soft EOS and J0592-0607 and J0740+6620, which require high-density stiff EOS. The addition of J0592-0607's mass measurement to the constraints has no significant impact on the parameter space of the admissible EOS, but allows us to explore the potential effect of pulsars more massive than J0740+6620, if such exists. We then examine the occurrence of twin configurations and quantify the ranges of masses and radii that they can possess. It is shown that including J1731-347 data favors EOSs that predict low-mass twins with M ≲ 1.3 M ⊙ that can be realized if the deconfinement transition density is low. If combined with large speed of sound in quark matter such models allow for maximum masses of hybrid stars in 2.0–2.6 M ⊙ .

Astronomy & Astrophysics↗

Accelerated Fuel Qualification of Fast Modular Reactor Fuel in a Thermal Reactor: Modeling and Simulation Paired with Irradiation Testing

The accelerated fuel qualification (AFQ) methodology is applied by simulating accelerated fuel tests of the General Atomics Electromagnetic Systems’ fuel system for its 44-MW(electric) gas-cooled, fast-spectrum fast modular reactor (FMR). This fuel is comprised of UO 2 pellets in SiGA® cladding, a silicon carbide ceramic matrix composite. Fast reactors, like the FMR, offer many benefits, including high fuel utilization and flexibility, but may require a lengthy material design process if tests are performed using fast neutron irradiation alone. A thermal neutron irradiation can instead be used to rapidly test how well key components of the current material models extend to high burnup. Thermal neutrons produce a different radial power distribution within the pin than fast neutrons. However, the temperature and burnup values for the two neutron types are comparable, and the differences between the simulated fuel responses are relatively small, demonstrating the weak sensitivity of the physics-based fuel model calculations on the neutron type and the irradiation rate. Furthermore, the deformation of the SiGA cladding saturates after about 1 displacement per atom for both neutron spectra. In an accelerated fuel test, the irradiation time required to reach the target fuel burnup can be reduced by a factor of 3 by using a small rodlet with a 45% smaller pellet diameter while maintaining the same linear power. Therefore, the time for data collection up to high burnup can be significantly reduced while maintaining the same temperature profile, which largely determines the material response. Tests of fuel rodlets of standard and compact size will be carried out in the Idaho National Laboratory’s Advanced Test Reactor (ATR), including full size and compact rodlets with varying gap sizes. By applying physics-based mechanistic modeling and simulation in accordance with the AFQ methodology, this type of compact rodlet testing in a thermal test reactor captures the necessary phenomena to test fuel material models up to high burnup and to simulate the expected impact of fast neutron radiation on the fuel in FMR operations. Finally, this approach to testing fast reactor fuels in existing thermal test reactors, paired with advanced physics-based mechanistic modeling and simulation, is expected to be applicable to a range of advanced fuels and will decrease the overall fuel qualification timeframe from decades to years.

Advanced test reactor (ATR)↗

Universal bounds on CFT Distance Conjecture

For any unitary conformal field theory in two dimensions with the central charge c, we prove that, if there is a nontrivial primary operator whose conformal dimension ∆ vanishes in some limit on the conformal manifold, the Zamolodchikov distance t to the limit is infinite, the approach to this limit is exponential ∆ = exp(−αt + O(1)), and the decay rate obeys the universal bounds c−1/2 ≤ α ≤ 1. In the limit, we also find that an infinite tower of primary operators emerges without a gap above the vacuum and that the conformal field theory becomes locally a tensor product of a sigma-model in the large radius limit and a compact theory. As a corollary, we establish a part of the Distance Conjecture about gravitational theories in three-dimensional anti-de Sitter space. In particular, our bounds on α indicate that the emergence of exponentially light states is inevitable as the moduli field corresponding to t rolls beyond the Planck scale along the steepest path and that this phenomenon can begin already at the curvature scale of the bulk geometry. We also comment on implications of our bounds for gravity in asymptotically flat spacetime by taking the flat space limit and compare with the Sharpened Distance Conjecture.

AdS-CFT Correspondence↗

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗

Preliminary Reporting of Thermophysical Property Measurements for the Ghareb Formation

Accurate knowledge of thermophysical properties of rock is vital to develop meaningful models of high level nuclear waste emplacement scenarios. The Israel Atomic Energy Commission is considering storing high level nuclear waste in the Ghareb formation, a porous kerogen bearing chalk. Sandia is supporting this effort with an evolving lab- based geomechanics testing program. We have completed measurements of thermal properties up to 275C and room temperature hydrostatic compaction measurements. We report thermal conductivity, thermal diffusivity, specific heat, and mass loss from our thermal measurements, and we report bulk moduli and porosity loss from our compaction measurements. These values are crucial for the numerical models to simulate heat transfer and formation compressibility around a heat generating repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparison of Fission Product Release Predictions using PARFUME with Results from the AGR-2 Irradiation Experiment

The PARFUME (PARticle FUel ModEl) code was used to predict fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the second irradiation experiment (AGR-2) of the Advanced Gas Reactor Fuel Development and Qualification Program. The PARFUME model for the Advanced Gas reactor (AGR)GR-2 experiment used the fuel compact volume average temperature for each of the 560 days of irradiation to calculate the release of fission products of silver, cesium, and strontium from a representative particle for a select number of AGR -2 compacts. In addition, safety tests were performed on 15 compacts ranging from 1500 to 1800°C to determine fission product release at temperatures that bound reactor accident conditions. PARFUME was used to calculate the fission product release of silver, cesium, strontium, and krypton during these safety tests. Post-irradiation examination (PIE) measurements provided the data on release of fission products from the fuel compacts and fuel particles, and retention of fission products in the compacts outside the silicon carbide (SiC) layer. The predicted fraction release from PARFUME was then compared to PIE measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling and analysis of the tritium fuel cycle for ARC- and STEP-class D-T fusion power plants

Abstract The limited tritium resources available for the first fusion power plants (FPPs) make fuel self-sufficiency and tritium inventory minimization leading issues in FPP design. This work builds on the model proposed by Abdou et al (2020 Nucl. Fusion 61 013001), which analyzed the fuel cycle (FC) of a DEMOnstration nuclear FPP-class FPP with a time-dependent system-level model. Here, we use a modified version of their model to analyze the FC of an Affordable, Robust, Compact (ARC)-class tokamak and two versions of a Spherical Tokamak for Energy Production (STEP)-class tokamak. The ARC-class tokamak breeds tritium in a 2LiF + BeF 2 liquid immersion blanket, while the STEP-class tokamak breeds tritium utilizing either a liquid-lithium blanket design or an encapsulated breeding blanket. A time-dependent system-level model is developed in Matlab Simulink ® to simulate the evolution of tritium flows and tritium inventories in the FC. The main goals of this work are to assess tritium self-sufficiency of the ARC- and STEP-class designs and to determine quantitative design requirements that can be used to analyze the adequacy of a proposed FC system. These design requirements are aimed at achieving a low tritium inventory doubling time ( t d ) and a low start-up inventory ( I s t a r t u p ) while keeping the required tritium breeding ratio (TBR r ) as low as possible. We also consider how improvements in FC technology and POs affect TBR r and I s t a r t u p . The model results show that TBR r for ARC- and STEP-class FPPs should be achievable if the tritium burn efficiency (TBE) reaches 0.5%–1% (TBR r < 1.2). This assumes significant, but attainable, improvements over current abilities. However, the model results indicate that an FPP must achieve ambitious performance targets, including FPP availability > 70%, tritium processing time < 4 h, and the implementation of direct internal recycling (DIR). If future research yields major improvements to achievable TBE, it may be possible to achieve tritium self-sufficiency while operating at lower availability and without implementing DIR.

Physics↗

Core formation in high-z massive haloes: heating by post-compaction satellites and response to AGN outflows

ABSTRACT Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass ≥1012M⊙ at z ∼ 2. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing active galactic nucleus (AGN)-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses >1011.3M⊙, which undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion on to haloes ≥1012M⊙, these satellites heat up the cusps in half a virial time at z ∼ 2. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGNs indeed become effective in haloes ≥1012M⊙, where the black hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot circumgalactic medium. For simulations to reproduce the dynamical friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black hole accretion and clumpy response to AGN.

Dekel, Avishai↗