Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnitude”

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 1,189 records · Page 66

Airborne Release Fractions from Surrogate Nuclear Waste Fires Containing Lanthanide Nitrates and Depleted Uranium Nitrate in 30% Tributyl Phosphate in Kerosene

Airborne contaminants from fires containing nuclear waste represent significant health hazards and shape the design and operation of nuclear facilities. Much of the data used to formulate DOE-HDBK-3010-94, “Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities,” from the U.S. Department of Energy, were taken over 40 years ago. The objectives of this study were to reproduce experiments from Pacific Northwest Laboratories conducted in June 1973 employing current aerosol measurement methods and instrumentation, develop an enhanced understanding of particulate formation and transport from fires containing nuclear waste, and provide modeling and experimental capabilities for updating current standards and practices in nuclear facilities. A special chamber was designed to conduct small fires containing 25 mL of flammable waste containing lutetium nitrate, ytterbium nitrate, or depleted uranium nitrate. Carbon soot aerosols showed aggregates of primary particles ranging from 20 to 60 nm in diameter. In scanning electron microscopy, ~200-nm spheroidal particles were also observed dispersed among the fractal aggregates. The 200-nm spherical particles were composed of metal phosphates. Airborne release fractions (ARFs) were characterized by leaching filter deposits and quantifying metal concentrations with mass spectrometry. The average mass-based ARF for 238 U experiments was 1.0 × 10 -3 with a standard deviation of 7.5 × 10 -4 . For the original experiments, DOE-HDBK-3010-94 states, “Uranium ARFs range from 2 × 10 -4 to 3 × 10 -3 , an uncertainty of approximately an order of magnitude.” Therefore, current measurements were consistent with DOE-HDBK-3010-94 values. ARF values for lutetium and ytterbium were approximately one to two orders of magnitude lower than 238 U. Furthermore, metal nitrate solubility may have varied with elemental composition and temperature, thereby affecting ARF values for uranium surrogates (Yb and Lu). In addition to ARF data, solution boiling temperatures and evaporation rates can also be deduced from experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparative Evaluation of Spectral Methods for Robust Reactor Noise Estimation

Reactor noise analysis provides a noninvasive means to determine neutron kinetic parameters from stochastic fluctuations in detector signals. However, standard cross-power spectral density (CPSD) analyses can be sensitive to numerical processing choices, which may introduce processing-dependent systematic shifts in estimates of the prompt neutron decay constant (α) and limit reproducibility. This study uses a hybrid multitaper–Welch spectral estimator to analyze subcritical noise measurements from a fast-spectrum critical assembly. The decay constant α was extracted using three frequency-domain methods: the CPSD, the magnitude-squared coherence (MSC), and the generalized magnitude-squared coherence (GMSC). These coherence-based estimators normalize detector auto-spectral structure and are expected to reduce the sensitivity of fitted α values to processing parameters. A Sobol global sensitivity analysis identified which numerical inputs most strongly influence the fitted values of α. All estimators produced a linear dependence of α on inverse count rate, with delayed-critical extrapolations near 1.7 × 10 4 s −1 , in agreement within 8% of MCNP6.3 KOPTS benchmark calculations. Sensitivity results show that while the CPSD depends on both time-bin width and taper selection, the MSC and GMSC are dominated by time-bin width alone, indicating reduced parameter coupling and greater robustness to processing variability. These findings demonstrate the feasibility and practical value of coherence-based spectral estimators for extracting α from reactor noise and support their broader application to multi-detector and irregular datasets in subcritical system characterization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Physics-informed neural networks for heterogeneous poroelastic media

This study presents a novel physics-informed neural network (PINN) framework for modeling poroelasticity in heterogeneous media with material interfaces. The approach introduces a composite neural network (CoNN) where separate neural networks predict displacement and pressure variables for each material. While sharing identical activation functions, these networks are independently trained for all other parameters. To address challenges posed by heterogeneous material interfaces, the CoNN is integrated with the Interface-PINNs (I-PINNs) framework (Sarma et al., Comput. Methods Appl. Mech. Eng. 429: 117135, 2024), allowing different activation functions across material interfaces. Further, this ensures accurate approximation of discontinuous solution fields and gradients. Performance and accuracy of this combined architecture were evaluated against the conventional PINNs approach, a single neural network (SNN) architecture, and the eXtended PINNs (XPINNs) framework through two one-dimensional benchmark examples with discontinuous material properties. The results show that the proposed CoNN with I-PINNs architecture achieves an RMSE that is two orders of magnitude better than the conventional PINNs approach and is at least 40 times faster than the SNN framework. Compared to XPINNs, the proposed method achieves an RMSE at least one order of magnitude better and is 40% faster.

42 ENGINEERING↗

Gravity of two photon decay and its quantum coherence

A linear analytical solution is derived for the gravitational shock wave produced by a particle of mass M that decays into a pair of null particles. The resulting space-time is shown to be unperturbed and isotropic, except for a discontinuous perturbation on a spherical null shell. Formulae are derived for the perturbation as a function of polar angle, as measured by an observer at the origin observing clocks on a sphere at distance R. The effect of the shock is interpreted physically as an instantaneous displacement in time and velocity when the shock passes the clocks. The time displacement is shown to be anisotropic, dominated by a quadrupole harmonic aligned with the particle-decay axis, with a magnitude δτ ~ GM/c 3 , independent of R. The velocity displacement is isotropic. The solution is used to derive the gravitational effect of a quantum state with a superposition of a large number of randomly oriented, statistically isotropic particle decays. This approach is shown to provide a well-controlled approximation to estimate the magnitude of gravitational fluctuations in systems composed of null point particles up to the Planck energy in a causal diamond of duration τ = 2R/c, as well as quantum-gravitational fluctuations of black holes and cosmological horizons. Coherent large-angle quantum distortions of macroscopic geometry from fluctuations up to the Planck scale are shown to grow linearly with the duration, with a variance $\langle$δτ 2 $\rangle$ ~ τtP much larger than that produced in models without causal quantum coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗

Heavy ion irradiation induced failure of gallium nitride high electron mobility transistors: effects of in-situ biasing

Abstract While radiation is known to degrade AlGaN/GaN high-electron-mobility transistors (HEMTs), the question remains on the extent of damage governed by the presence of an electrical field in the device. In this study, we induced displacement damage in HEMTs in both ON and OFF states by irradiating with 2.8 MeV Au 4+ ion to fluence levels ranging from 1.72 × 10 10 to 3.745 × 10 13 ions cm −2 , or 0.001–2 displacement per atom (dpa). Electrical measurement is done in situ , and high-resolution transmission electron microscopy (HRTEM), energy dispersive x-ray (EDX), geometrical phase analysis (GPA), and micro-Raman are performed on the highest fluence of Au 4+ irradiated devices. The selected heavy ion irradiation causes cascade damage in the passivation, AlGaN, and GaN layers and at all associated interfaces. After just 0.1 dpa, the current density in the ON-mode device deteriorates by two orders of magnitude, whereas the OFF-mode device totally ceases to operate. Moreover, six orders of magnitude increase in leakage current and loss of gate control over the 2-dimensional electron gas channel are observed. GPA and Raman analysis reveal strain relaxation after a 2 dpa damage level in devices. Significant defects and intermixing of atoms near AlGaN/GaN interfaces and GaN layer are found from HRTEM and EDX analyses, which can substantially alter device characteristics and result in complete failure.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Structure-property correlations and scaling in the magnetic and magnetocaloric properties of GdCrO 3 particles

The structure, magnetic, and magnetocaloric (MC) properties of orthorhombic nanocrystalline GdCrO 3 with six particle sizes: < d > = 87, 103, 145, 224, 318, and 352 nm are reported. The particle size was tailored by annealing under different temperatures and estimated by scanning electron microscopy. With increase in < d >, Goldschmidt tolerance factor t, orthorhombic strain s, and out-of-plane Cr–O 1 –Cr bond angle first decrease, reaching minimum values for < d > = 224 nm, and then increase for sample with < d > = 318 nm and 352 nm, thus showing a V-shaped variation. Temperature dependence of the magnetization (M) reveals an antiferromagnetic transition at ${T}_{\text{N}}^{\text{Cr}}\sim 168$ K for < d > ≥ 224 nm and ${T}_{\text{N}}^{\text{Cr}}\sim 167$ K for < d > < 224 nm and an essentially d-independent spin-reorientation at T SR = 9 K. M measured at 5 K and 7 T first increases with increase in < d >, reaching maximum value for sample with < d > = 224 nm, and then decreases for samples with < d > = 318 nm and 352 nm, showing an inverted-V variation with < d >. Similar < d >-dependence is observed for the magnetic entropy change (MEC) and relative cooling power (RCP) showing a close relationship between the structural and magnetic properties of GdCrO 3 nanoparticles investigated here. The 224 nm sample with the minimum values of t, s, and Cr–O 1 –Cr bond angle exhibits the maximum value of MEC (–ΔS) = 37.8 J kg –1 K –1 at 5 K under a field variation (ΔH) of 7 T and its large estimated RCP of 623.6 J Kg –1 is comparable with those of typical MC materials. Both (–ΔS) and RCP are shown to scale with the saturation magnetization MS, suggesting that MS is the crucial factor controlling their magnitudes. Assuming (–ΔS) ~ (ΔH) n , the temperature dependence of n for the six samples are determined, n varying between 1.3 at 5 K to n = 2.2 at 130 K in line with its expected magnitudes based on mean-field theory. Furthermore, these results on structure-property correlations and scaling in GdCrO 3 suggest that its MC properties are tunable for potential low-temperature magnetic refrigeration applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

EMC3-EIRENE simulation of first wall recycling fluxes in W7-X with relation to H-alpha measurements

In the Wendelstein 7-X stellarator, the main locations of particle sources are expected to be the carbon divertors, baffles and graphite heat shield first wall. In this paper, the heat shield is implemented in EMC3-EIRENE to understand the expected areas and magnitudes of the recycling flux to this component. It is found that in the simulation the heat shield is not a significant source of recycling neutrals. The areas of simulated recycling flux are shown to correlate well with footprints of plasma-wetting seen in post-experimental campaign in-vessel inspection photos. EMC3-EIRENE reconstruction of line-integrated H-alpha measurements at the heat shield indicate that the majority of emission does not come from local recycling neutrals. Rather, the H-alpha signals at the heat shield are dominated by ionization of neutrals which have leaked from the divertor/baffle region into the midplane. The magnitude of the H-alpha line emission from the synthetic reconstruction is consistent with the experiment, indicating that a large overestimation of heat shield recycling would occur if these measurements were assumed to be from local recycling sources. In the future, it may be possible to obtain some information of local recycling from the heat shield since it was found that the majority of the recycling flux occurs on two well-localized areas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Recent progress on neoclassical impurity transport in stellarators with implications for a stellarator reactor

Accumulation of impurities in the core of the plasma is a potential problem for any fusion reactor, and is thought to be especially problematic for stellarators. We use the collisional transport code SFINCS to calculate, in a reactor stellarator scenario, the effects of flux-surface electrostatic potential variation on tungsten transport. The potential variation is found to have a large effect on the transport, reversing the sign of the transport coefficients related to the radial electric field and the bulk ion temperature gradient, and also increasing the magnitude of the former by two orders of magnitude. We optimize the electrostatic potential variation to minimize impurity accumulation, resulting in a 20%–50% reduction of the inward impurity flux. Furthermore, we investigate whether the electrostatic potential variation caused by fast particles due to ion-cyclotron resonance heating can be used to reduce the impurity accumulation, but find the effect to be small in the cases considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D simulations of inertial confinement fusion implosions part 2: systematic flow anomalies and impact of low modes on performances in OMEGA experiments

We present the first 3D radiation-hydrodynamic simulations of directly driven inertial confinement fusion implosions with an inline package for polarized crossed-beam energy transfer, which were used to assess the impact of the current distributed polarization rotators (DPRs) on OMEGA as well as other known sources of asymmetry. Applied to OMEGA implosions, the simulations predict bang times with no need for ad hoc multipliers, as well as yields—if you separately account for the impacts of imprint and fuel age. The magnitude of the flow is well reproduced when the low mode sources are large, whereas the modeling of the stalk is thought to be required to match the flow magnitude in the remaining cases. For the cases explored in more detail, polarized cross-beam energy transfer (CBET)—the only known systematic drive asymmetry, brought the results closest to the measured flow vectors. The remaining discrepancies are shown to be stemming from the limited knowledge of the laser pointing modes. For typical current levels of beam mispointing, power imbalance, target offset, and asymmetry caused by polarized CBET, low modes degrade the yield by more than 40%. The current strategy of attempting to compensate the mode-1 asymmetry with a preimposed target offset recovers only about one-third of the losses caused by the low modes due to the dynamic nature of the multiple asymmetries and the presence of low modes other than l = 1. Furthermore, addressing the root causes of the drive asymmetries is apt to be more beneficial. To that end, one possible solution to the specific issue of polarized CBET (10 µm DPRs) is shown to work well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spectral transition of multiscale turbulence in the tokamak pedestal

The transition in the turbulence spectrum from ion-scale dominated regimes to multiscale transport regimes that couple ion and electron scales is studied with gyrokinetic simulations of turbulent transport. The simulations are based on DIII-D high-confinement mode (H-mode) plasma parameters in the tokamak pedestal. The transition is initiated by varying the ion temperature gradient. To our knowledge, no full multiscale simulations of pedestal-like transport have been done previously. The experimental parameters lie in a bifurcation region between the two regimes. At long wavelengths, a complex, ion-direction hybrid mode is the dominant linearly unstable drift wave, while an electron temperature gradient-driven mode is unstable at short wavelengths. In the transition from the multiscale branch to the ion-scale branch, the magnitude of the ion-scale poloidal wavenumber spectrum of the nonlinear turbulent energy flux increases and the magnitude of the high-wavenumber spectrum decreases. Here, the decrease in the electron-scale transport is due to nonlinear mixing with ion-scale fluctuations and the ion-scale-driven zonal flows. A shift in the total energy associated with the fluctuating electrostatic potential intensity from dominantly drift kinetic energy in the multiscale regime to dominantly potential intensity in the ion-scale regime is well-correlated with the trend in the total energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementation of Pfirsch–Schlüter parallel flow in x-ray imaging crystal spectrometer inversion analysis

The x-ray imaging crystal spectrometer (XICS) tomographic inversion code for Wendelstein 7-X (W7-X) has been modified to consider the effects of parallel flows and has been applied to analyze measurements taken during recent experimental campaigns. Previous analysis neglected the effects of parallel flows due to the primarily perpendicular geometry of the sightlines and the small magnitude predicted by neoclassical theory. To reconsider these effects, the incompressibility condition for plasma flows is used to calculate the parallel Pfirsch–Schlüter flow component for the equilibrium configuration. By incorporating this condition along with the geometry of the sightlines—i.e. the fractional contributions of perpendicular and parallel flows—, an updated expression for the measured flow is used for the profile inversion. Application of this modified inversion code to data from W7-X shows that the magnitude of the radial electric field and the flux surface averaged perpendicular flow are reduced by approximately a factor of 2 and brought into better agreement with neoclassical predictions and charge exchange recombination spectroscopy measurements.

Pfirsch–Schlüter flows↗

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗

Strong surface magnetic field generation in relativistic short pulse laser–plasma interaction with an applied seed magnetic field

Abstract While plasma often behaves diamagnetically, we demonstrate that the laser irradiation of a thin opaque target with an embedded target-transverse seed magnetic field B seed can trigger the generation of an order-of-magnitude stronger magnetic field with opposite sign at the target surface. Strong surface field generation occurs when the laser pulse is relativistically intense and results from the currents associated with the cyclotron rotation of laser-heated electrons transiting through the target and the compensating current of cold electrons. We derive a predictive scaling for this surface field generation, B gen ∼ −2 πB seed Δ x / λ 0 (in the large spot size limit), where Δ x is the target thickness and λ 0 is the laser wavelength, and conduct 1D and 2D particle-in-cell simulations to confirm its applicability over a wide range of conditions. We additionally demonstrate that both the seed and surface-generated magnetic fields can have a strong impact on application-relevant plasma dynamics, for example substantially altering the overall expansion and ion acceleration from a μ m-thick laser-irradiated target with a kilotesla-level seed magnetic field.

Weichman, K. (ORCID:0000000234877922)↗

Spontaneous skyrmion conformal lattice and transverse motion during dc and ac compression

Abstract We use atomistic-based simulations to investigate the behavior of ferromagnetic skyrmions being continuously compressed against a rigid wall under dc and ac drives. The compressed skyrmions can be annihilated close to the wall and form a conformal crystal with both a size and a density gradient, making it distinct from conformal crystals observed previously for superconducting vortices and colloidal particles. For both dc and ac driving, the skyrmions can move transverse to the compression direction due to a combination of density and size gradients. Forces in the compression direction are converted by the Magnus force into transverse motion. Under ac driving, the amount of skyrmion annihilation is reduced and we find a skyrmion Magnus ratchet pump. We also observe shear banding in which skyrmions near the wall move up to twice as fast as skyrmions further from the wall. When we vary the magnitude of the applied drive, we find a critical current above which the skyrmions are completely annihilated during a time scale that depends on the magnitude of the drive. By varying the magnetic parameters, we find that the transverse motion is strongly dependent on the skyrmion size. Smaller skyrmions are more rigid, which interferes with the size gradient and destroys the transverse motion. We also confirm the role of the size gradient by comparing our atomistic simulations with a particle-based model, where we find that the transverse motion is only transient. Our results are relevant for applications where skyrmions encounter repulsive magnetic walls, domain walls, or interfaces.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Milky Way white dwarfs as sub-GeV to multi-TeV dark matter detectors

Here, we show that Milky Way white dwarfs are excellent targets for dark matter (DM) detection. Using Fermi and H.E.S.S. Galactic center gamma-ray data, we investigate sensitivity to DM annihilating within white dwarfs into long-lived or boosted mediators and producing detectable gamma rays. Depending on the Galactic DM distribution, we set new constraints on the spin-independent scattering cross section down to 10 -45 -10 -41 cm 2 in the sub-GeV DM mass range, which is multiple orders of magnitude stronger than existing limits. For a generalized NFW DM profile, we find that our white dwarf constraints exceed spin-independent direct detection limits across most of the sub-GeV to multi-TeV DM mass range, achieving sensitivities as low as about 10 -46 cm 2 . In addition, we improve earlier versions of the DM capture calculation in white dwarfs, by including the low-temperature distribution of nuclei when the white dwarf approaches crystallization. This yields smaller capture rates than previously calculated by a factor of a few up to two orders of magnitude, depending on white dwarf size and the astrophysical system.

79 ASTRONOMY AND ASTROPHYSICS↗

Color Transformations of Photometric Measurements of Galaxies in Optical and Near-infrared Wide-field Imaging Surveys*

Abstract Over the past 2 decades, wide-field photometric surveys in optical and infrared domains reached a nearly all-sky coverage thanks to numerous observational facilities operating in both hemispheres. However, subtle differences among exact realizations of Johnson and SDSS photometric systems require one to convert photometric measurements into the same system prior to analysis of composite data sets originating from multiple surveys. It turns out that the published photometric transformations lead to substantial biases when applied to integrated photometry of galaxies from the corresponding catalogs. Here we present photometric transformations based on piece-wise linear approximations of integrated photometry of galaxies in the optical surveys SDSS, DECaLS, BASS, MzLS, DES, DELVE, KiDS, VST ATLAS, and the near-infrared surveys UKIDSS, UHS, VHS, and VIKING. We validate our transformations by constructing k -corrected color–magnitude diagrams of non-active galaxies and measuring the position and tightness of the “red sequence.” We also provide transformations for aperture magnitudes and show how they are affected by the image quality difference among the surveys. We present the implementation of the derived transformations in python and idl and also a web-based color transformation calculator for galaxies. By comparing DECaLS and DES, we identified systematic issues in DECaLS photometry for extended galaxies, which we attribute to the photometric software package used by DECaLS. As an application of our method, we compiled two multi-wavelength photometric catalogs for over 200,000 low- and intermediate-redshift galaxies originating from CfA FAST and Hectospec spectral archives.

Astronomy & Astrophysics↗

The Electromagnetic Characteristics of the Tianlai Cylindrical Pathfinder Array

A great challenge for 21 cm intensity mapping experiments is the strong foreground radiation which is orders of magnitude brighter than the 21 cm signal. Removal of the foreground takes advantage of the fact that its frequency spectrum is smooth while the redshifted 21 cm signal spectrum is stochastic. However, a complication is the non-smoothness of the instrument response. This paper describes the electromagnetic simulation of the Tianlai cylinder array, a pathfinder for 21 cm intensity mapping experiments. Due to the vast scales involved, a direct simulation requires a large amount of computing resources. We have made the simulation practical by using a combination of methods: first simulate a single feed, then an array of feed units, finally with the feed array and a cylindrical reflector together, obtain the response for a single cylinder. We studied its radiation pattern, bandpass response and the effects of mutual coupling between feed units, and compared the results with observation. Many features seen in the measurement result are reproduced well in the simulation, especially the oscillatory features which are associated with the standing waves on the reflector. The mutual coupling between feed units is quantified with S-parameters, which decrease as the distance between the two feeds increases. Based on the simulated S-parameters, we estimate the correlated noise which has been seen in the visibility data, and the results show very good agreement with the data in both magnitude and frequency structures. These results provide useful insights on the problem of 21 cm signal extraction for real instruments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗