Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scenario”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Benzyl trichloroacetimidates as derivatizing agents for phosphonic acids related to nerve agents by EI-GC-MS during OPCW proficiency test scenarios

Abstract The use of benzyl trichloroacetimidates for the benzylation of phosphonic acid nerve agent markers under neutral, basic, and slightly acidic conditions is presented. The benzyl-derived phosphonic acids were detected and analyzed by Electron Ionization Gas Chromatography–Mass Spectrometry (EI-GC–MS). The phosphonic acids used in this work included ethyl-, cyclohexyl- and pinacolyl methylphosphonic acid, first pass hydrolysis products from the nerve agents ethyl N -2-diisopropylaminoethyl methylphosphonothiolate (VX), cyclosarin (GF) and soman (GD) respectively. Optimization of reaction parameters for the benzylation included reaction time and solvent, temperature and the effect of the absence or presence of catalytic acid. The optimized conditions for the derivatization of the phosphonic acids specifically for their benzylation, included neutral as well as catalytic acid (< 5 mol%) and benzyl 2,2,2-trichloroacetimidate in excess coupled to heating the mixture to 60 °C in acetonitrile for 4 h. While the neutral conditions for the method proved to be efficient for the preparation of the p -methoxybenzyl esters of the phosphonic acids, the acid-catalyzed process appeared to provide much lower yields of the products relative to its benzyl counterpart. The method’s efficiency was tested in the successful derivatization and identification of pinacolyl methylphosphonic acid (PMPA) as its benzyl ester when present at a concentration of ~ 5 μg/g in a soil matrix featured in the Organisation for the Prohibition of Chemical Weapons (OPCW) 44th proficiency test (PT). Additionally, the protocol was used in the detection and identification of PMPA when spiked at ~ 10 μg/mL concentration in a fatty acid-rich liquid matrix featured during the 38 th OPCW-PT. The benzyl derivative of PMPA was partially corroborated with the instrument's internal NIST spectral library and the OPCW central analytical database (OCAD v.21_2019) but unambiguously identified through comparison with a synthesized authentic standard. The method’s MDL (LOD) values for the benzyl and the p -methoxybenzyl pinacolyl methylphosphonic acids were determined to be 35 and 63 ng/mL respectively, while the method’s Limit of Quantitation (LOQ) was determined to be 104 and 189 ng/mL respectively in the OPCW-PT soil matrix evaluated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Context. Analysis of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE) data suggests the existence of a clear distinction between two sequences of disc stars in the [α/Fe] versus [Fe/H] abundance ratio space, known as the high- and low-α sequence, respectively. This dichotomy also emerges from an analysis of the vertical distribution of the [α/Fe] abundance ratio. Aims. We aim to test whether the revised two-infall chemical evolution models designed to reproduce the low- and high-α sequences in the [α/Fe] versus [Fe/H] ratios in the solar neighbourhood are also capable of predicting the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Methods. Along with the chemical composition of the simple stellar populations born at different Galactic times predicted by our reference chemical evolution models in the solar vicinity, we provide their maximum vertical height above the Galactic plane |zmax| computed assuming the relation between the vertical action and stellar age in APOGEE thin-disc stars. Result. The vertical distribution of the [Mg/Fe] abundance ratio predicted by the reference chemical evolution models is in agreement with that observed when combining the APOGEE DR16 data (chemical abundances) with the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-α sequence stars). Including the high-α disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z|< 2 kpc. However, our model predicts an overly flat (almost constant) growth of the maximum vertical height |zmax| quantity as a function of [Mg/Fe] for high-α objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are that: (i) the APOGEE sample with |z|< 2 kpc is more likely than ours to be contaminated by halo stars, causing the median values to be kinematically hotter, and (ii) external perturbations – such as minor mergers – that the Milky Way experienced in the past could have heated up the disc, and the heating of the orbits cannot be modeled by only scattering processes. Assuming a disc dissection based on chemistry for APOGEE-DR16 stars (|z|< 2 kpc), the observed |zmax| distributions for high-α and low-α sequences are in good agreement with our model predictions if we consider the errors in the vertical action estimates in the calculation. Moreover, a better agreement between predicted and observed stellar distributions at different Galactic vertical heights is achieved if asteroseismic ages are included as a constraint in the best-fit model calculations. Conclusions. The signature of a delayed gas infall episode, which gives rise to a hiatus in the star formation history of the Galaxy, are imprinted both in the [Mg/Fe] versus [Fe/H] relation and in vertical distribution of [Mg/Fe] abundances in the solar vicinity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Reduced-model scrape-off layer turbulence (nSOLT) simulations comparing three fueling scenarios

The 2D scrape-off-layer turbulence code (nSOLT) includes 1D Boltzmann neutral–plasma interactions, a model of divertor recycling (introduced here), and a fixed source of plasma concentrated at the core-side boundary. Three fueling methods are considered herein: (1) neutral injection in the far-SOL is accomplished by specifying the density of Franck–Condon distributed neutrals streaming in from the boundary. (2) Divertor recycling is modeled by injecting a fraction of the particle parallel flux in the scrape-off layer (SOL) back into the edge region as a source of plasma. (3) A constant source fuels the edge plasma from the core-side boundary to model pellet injection. For machine parameters (B, R m , and L // ) illustrative of the MAST-U device, and for a deuterium plasma, turbulent equilibria are obtained that share the same plasma fueling rate for each of the three fueling methods, with only one of the sources on in each case. In the presence of self-consistent turbulence, quasi-steady plasma and neutral (deuterium) profiles, fueling efficiencies, SOL transparencies, and heat flux widths are compared. Characteristics of the turbulent fluctuations, including skewness, cross-phases, and power spectra, are described. The calculated fueling efficiencies, SOL transparencies to neutral penetration, and many of the turbulent properties are remarkably similar for all three fueling methods despite significant differences in the plasma profiles. The nonlinear states of the three cases are dominated by separatrix-spanning vortex cells that control particle and heat losses into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Empirical probability and machine learning analysis of m , n = 2, 1 tearing mode onset parameter dependence in DIII-D H-mode scenarios

m, n = 2, 1 tearing mode onset empirical probability and machine learning analyses of a multiscenario DIII-D database of over 14 000 H-mode discharges show that the normalized plasma beta, the rotation profile, and the magnetic equilibrium shape have the strongest impact on the 2,1 tearing mode stability, in qualitative agreement with neoclassical tearing modes (m and n are the poloidal and toroidal mode numbers, respectively). In addition, 2,1 tearing modes are most likely to destabilize when n > 1 tearing modes are already present in the core plasma. The covariance matrix of tearing sensitive plasma parameters takes a nearly block-diagonal form, with the blocks incorporating thermodynamic, current and safety factor profile, separatrix shape, and plasma flow parameters, respectively. This suggests a number of paths to improved stability at fixed pressure and edge safety factor primarily by preserving a minimum of 1 kHz differential rotation, increasing the minimum safety factor above unity, using upper single null magnetic configuration, and reducing the core impurity radiation. In addition, lower triangularity, lower elongation, and lower pedestal pressure may also help to improve stability. Furthermore, the electron and ion temperature, collisionality, resistivity, internal inductance, and the parallel current gradient appear to only weakly correlate with the 2,1 tearing mode onsets in this database.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Oxide dissolution and oxygen diffusion scenarios in niobium and implications on the Bean–Livingston barrier in superconducting cavities

We generalize a native Nb 2 O 5 dissolution model [G. Ciovati, Appl. Phys. Lett. 89, 022507 (2006)] to sequential overlayer dissolutions, multilayer dissolution, and realistic temperature profiles, which may be applicable to other materials. The model is applied to secondary ion mass spectrometry depth profile measurements for varying temperature profiles and two-step oxide dissolution in Nb and found to agree well. In the context of the Meissner screening response due to impurity profiles on the length scale of the London penetration depth, the shallow diffusion of O impurities results in a substantial decrease in the peak supercurrent density near the surface. In this framework, oxide dissolution and oxygen diffusion can account for a rise in peak supportable magnetic field in SRF cavities with baking time and a suppression after the optimal baking time is reached, in good agreement with peak-field baking temperatures and times as well as recent quench field measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Structural changes across thermodynamic maxima in supercooled liquid tellurium: a water-like scenario

Liquid polymorphism is an intriguing phenomenon which has been found in a few single-component systems, the most famous being water. By supercooling liquid Te to more than 130 K below its melting point and performing simultaneous small-angle and wide-angle X-ray scattering measurements, we observe clear maxima in its thermodynamic response functions around 615 K, suggesting the possible existence of liquid polymorphism. A close look at the underlying structural evolution shows the development of intermediate-range order upon cooling, most strongly around the thermodynamic maxima, which we attribute to bond-orientational ordering. The striking similarities between our results and those of water, despite the lack of hydrogen-bonding and tetrahedrality in tellurium, indicate that water-like anomalies may be a general phenomenon among liquid systems with competing bond- and density-ordering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of Postulated Accident Scenarios for the Transformational Challenge Reactor

The Transformational Challenge Reactor (TCR) is a helium-cooled, yttrium-hydride-moderated reactor that was designed for the U.S. Department of Energy Office of Nuclear Energy. A key objective of the TCR was to employ advanced manufacturing techniques in a nuclear system and demonstrate their potential for revolutionizing the nuclear reactor design process. One purpose of the present work is to demonstrate the safety of the TCR under postulated accidents. Based on RELAP5-3D and COMSOL analyses, the TCR remained below all current safety limits and far below the expected failure limits for the core materials. Another purpose of this work is to provide useful insights and recommendations regarding the application of RELAP5-3D to gas-cooled or other advanced reactors. A novel approach was implemented for simultaneously modeling conduction and radiation in RELAP5-3D, which was found to provide reasonable predictions of radial core, vessel, and ex-vessel heat transfer during postulated events. A multicode approach was also applied, in which high-fidelity COMSOL calculations were used to tune the radial heat transfer parameters in RELAP5-3D. The tuned RELAP5-3D model demonstrated comparable peak temperature predictions as COMSOL, despite a coarse treatment of the core in RELAP5-3D consisting of only two lumped heat structures. This high-fidelity tuning approach enabled enhanced accuracy as well as minimal complexity within the RELAP5-3D model, even for complex fuel geometric designs as in the TCR. Finally, investigations were made into the potential for flow reversal during a pressurized loss-of-forced-flow event in the TCR. The TCR is designed with downward helium flow through the core during normal operation. The RELAP5-3D model predicted that this downward flow would persist, without flow reversal, up to several days after the circulator trip. This was attributed to natural circulation hysteresis effects as have been noted in similar thermofluidic systems. Although flow stagnation and eventual reversal did not lead to unsafe TCR conditions, interesting spatial effects were observed which may have safety relevance for other reactor system designs and coolant types that are designed for downward core flow during normal operation, warranting closer investigation of the flow reversal phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Best-case scenarios for neutrino capture experiments

A direct discovery of the cosmic neutrino background would bring to a closure the searches for relic left-over radiation predicted by the Hot Big Bang cosmology. Recently, the KATRIN experiment put a limit on the local relic neutrino overdensity with respect to the cosmological predicted average value at η ≲ 10 11 [ Phys. Rev. Lett. 129 (2022) 011806]. In this work, we first examine to what extent such values of η are conceivable. We show that even under cavalier assumptions, a cosmic origin of η ≳ 10 4 seems out of reach (with the caveat of forming bound objects under a new force,) but find that a hypothetical local source of low-energy neutrinos could achieve η ~ 10 11 . Second, when such values are considered, we point out that the experimental signature in KATRIN and other neutrino-capture experiments changes, contrary to what has hitherto been assumed. Our results are model-independent and maximally accommodating as they only assume the Pauli exclusion principle. As intermittent physics target in the quest for C ν B detection, we identify an experimental sensitivity to η ~ 10 4 for which conceivable sources exist; to resolve the effect of a degenerate Fermi gas for such overdensity an energy resolution of 10 meV is required.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Color halo scenario of charmonium-like hybrids

The internal structures of J PC =1 -- ,(o,1,2) -+ charmonium-like hybrids are investigated under lattice QCD in the quenched approximation. We define the Bethe-Salpeter wave function (Φn(r)) in the Coulomb gauge as the matrix element of a spatially extended hybrid-like operator ($\bar{c}cg$) between the vacuum and n-th state for J PC , with r being the spatial separation between a localized $\bar{c}c$ component and the chromomagnetic strength tensor. These wave functions exhibit some similarities for states with the aforementioned different quantum numbers, and their r-behaviors (no node for the ground states and one node for the first excited states) imply that r can be a meaningful dynamical variable for these states. Additionally, the mass splittings of the ground states and first excited states of charmonium-like hybrids in these channels are obtained for the first time to be approximately 1.2-1.4 GeV. These results do not support the flux-tube description of heavy-quarkonium-like hybrids in the Born-Oppenheimer approximation. In contrast, a charmonium-like hybrid can be viewed as a "color halo" charmonium for which a relatively localized color octet $\bar{c}c$ is surrounded by gluonic degrees of freedom, which can readily decay into a charmonium state along with one or more light hadrons. The color halo picture is compatible with the decay properties of Y (42600) and suggests LHCb and BelleII to search for (0,1,2) -+ charmonium-like hybrids in Xc0,1,2N and J/ψω(φ) final states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development of an integrated core–edge scenario using the super H-mode

An optimized pedestal regime called the super-H (SH) mode is leveraged to couple a fusion relevant core plasma with a high density scrape-off layer appropriate for realistic reactor power exhaust solutions. Recent DIII-D experiments have expanded the operating space of the SH regime using advanced control algorithms and investigated optimization of impurity seeding, deuterium gas puffing, and 3D magnetic perturbations. Simultaneous real-time control of the pedestal density and radiated power with in-vessel coils and nitrogen seeding enable optimal coupled divertor and pedestal conditions. Four case studies are analysed with varied levels of radiated power in the divertor volume ranging from 0 (no seeding) to 8.5 MW radiated from carbon and nitrogen emission. Plasmas with a 4.5 MW radiated power target establish a radiative mantle, leading to divertor temperatures of ~16 eV while maintaining SH-mode, and with only marginal impact on the pedestal and core performance. Increased levels of N 2 seeding with a 7.5 MW radiated power target facilitate detachment onset and divertor temperatures <5 eV, with no degradation in stored energy and the operational point remaining inside the SH-mode channel for $ > 2.5{\tau }_{\text{E}}$. Finally, a 8.5 MW radiated power target leads to partial detachment, which is so far associated with the loss of access to SH-mode pedestal conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Studies of the outer-off-midplane lower hybrid wave launch scenario for plasma start-up on the TST-2 spherical tokamak

Abstract Establishment of an efficient central solenoid (CS) free tokamak plasma start-up method may lead to an economical fusion reactor. CS-free start-up using lower hybrid (LH) waves has been studied on the TST-2 spherical tokamak. Plasma current of about a quarter of CS-driven discharges has been obtained fully non-inductively using the outer-midplane and top LH launchers. Recently, an outer-off-midplane LH launcher was developed to achieve higher plasma current by optimizing for core absorption and minimal fast electron losses. Using the (outer-)off-midplane launcher, fully non-inductive plasma current start-up up to about 8 kA was achieved. Coupled ray-tracing and Fokker–Planck simulation was performed on equilibria reconstructed with an extended MHD model. It was found that the experimentally observed plasma current was in reasonable agreement with the numerical simulation. The simulation predicted appreciable orbit losses for the off-midplane launcher driven discharge at the present parameters, which was consistent with the experimentally observed x-ray radiation characteristics. The simulation showed that the current density was saturated for the present off-midplane launcher discharges and higher density and higher LH power was necessary to achieve higher plasma current.

Physics↗
DOE OSTI2025

Operational space for lower hybrid heating scenarios in the full tungsten environment of WEST

In tungsten—W—Environment in Steady-state Tokamak (WEST), the lower hybrid current drive (LHCD) system is key for achieving long pulse operation by providing most of the non-inductive plasma current, as well as a crucial source of electron heating. Therefore, determining the operational space for its application is fundamental. In the present study, the LHCD operational space is deeply analyzed for 0.5 MA pulses. This space is bounded by three limits: (i) the ratio of the LHCD power over density must be above a threshold to compensate tungsten radiation with enough core heating, (ii) the line-averaged density must be high enough to allow good coupling of the hybrid wave with the plasma, and (iii) fast electron ripple losses must be below a limit to avoid reaching a thermal threshold on plasma-facing components. If the tungsten radiation peak or burn-through phase is not safely overcome, a maximum electron temperature of 1.5 keV is obtained, confinement is degraded, and magnetohydrodynamic activity is frequently triggered, potentially causing a disruption. From experimental measurements and interpretative simulations, we highlight the main mechanisms that prevent the plasma from heating up during LHCD power ramp-up. Three parameters play a major role: plasma density, tungsten concentration and LHCD power deposition. A strategy to overcome this limitation is found: a precise density ramp-up performed simultaneously with the increase in LHCD power. Additionally, we show that boronization greatly facilitates the burn-through of tungsten by lowering its content during the heating phase. Finally, taking into account the three constraints given above, the LHCD operational space is determined at power ramp-up and during constant heating phases.

lower hybrid heating and current drive↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
DOE OSTI2023

Computational Fluid Dynamics Simulations to Assess Spatial Variability and Optimal Ventilation Scenarios for Biological Laboratory Exposures

A significant amount of uncertainty exists regarding potential human exposure to laboratory biomaterials and organisms in Biosafety Level 2 (BSL-2) research laboratories. Computational fluid dynamics (CFD) modeling is proposed as a way to better understand potential impacts of different combinations of biomaterials, laboratory manipulations, and exposure routes on risks to laboratory workers. Here, in this study, we use CFD models to simulate airborne concentrations of contaminants in an actual BSL-2 laboratory under different configurations. Results show that ventilation configuration, sampling location, and contaminant source location can significantly impact airborne concentrations and exposures. Depending on the source location and airflow patterns, the transient and time-integrated concentrations varied by several orders of magnitude. Contaminant plumes from sources located near a return vent (or exhaust like a fume hood or ventilated biosafety cabinet) are likely to be more contained than sources that are further from the exhaust. Having a direct flow between the source and the exhaust (through-flow condition) may reduce potential exposures to individuals outside the air flow path. Designing a BSL-2 room with ventilation and airflow patterns that maximize through-flow conditions to the return/exhaust vents and minimize dispersion and mixing throughout the room is, therefore, recommended. CFD simulations can also be used to assist in characterizing the impacts of supply and return vent locations, room layout, and source locations on spatial and temporal contaminant concentrations. In addition, proper placement of particle sensors can also be informed by CFD simulations to provide additional characterization and monitoring of potential exposures in BSL-2 facilities.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2023

Convergence of small scale Ly α structure at high- z under different reionization scenarios

ABSTRACT The Ly α forest (LAF) at z > 5 probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-z, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25 − 10 h−1 Mpc and 5 − 80 h−1 kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at z = 5.5. In observable Ly α power, convergence within 10 per cent is achieved for k < 0.1 s km–1, but larger k shows deviation of up to 20 per cent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early (z = 9) and later (z = 6) reionizations. We conclude that at z ∼ 5.5, resolutions of 10 kpc are necessary for convergence of LAF power at k < 0.1 s km–1, while higher-k modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.

79 ASTRONOMY AND ASTROPHYSICS↗