Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Directional Effect”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Direct imaging of electrical switching of antiferromagnetic Néel order in α-Fe 2 O 3 epitaxial films

In this work, we report the direct observation of switching of the Néel vector of antiferromagnetic (AFM) domains in response to electrical pulses in micron-scale Pt/α-Fe 2 O 3 Hall bars using photoemission electron microscopy. Current pulses lead to reversible and repeatable switching with the current direction determining the final state, consistent with Hall effect experiments that probe only the spatially averaged response. Current pulses also produce irreversible changes in domain structure, in and even outside the current path. In both cases only a fraction of the domains switch in response to pulses. Furthermore, the analysis of images taken with different x-ray polarizations shows that the AFM Néel order has an out-of-plane component in equilibrium that is important to consider in analyzing the switching data. These results show that—in addition to effects associated with spin-orbit torques from the Pt layer, which can produce reversible switching—changes in AFM order can be induced by purely thermal effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Enhancement of magnetic domain wall velocity via resonant dissipation of standing wave modes of domain wall structure with perpendicular magnetic anisotropy

The dynamic behaviors of magnetic domain walls have significant implications for developing advanced spintronic devices. In this study, we investigate the intriguing resonance phenomenon within the magnetic domain wall structure, focusing on the dissipation mechanism and its impact on domain wall dynamic motion. By applying a static external magnetic field, we observe a remarkable amplification of domain wall velocity through the resonant excitations of the flexure modes. The resonance of the flexure mode exerts an indirect influence on the domain wall velocity through energy dissipation, leading to a significant enhancement from the predictions of the one-dimensional Walker model. We establish a robust model platform rooted in energy dissipation allows comprehensive understanding the dramatic enhancement of domain wall velocity observed during width-directional flexure mode resonance. By directly quantifying energy dissipation, our approach provides highly effective and accurate estimates of domain wall velocity, surpassing previous methods that focus solely on the instantaneous dynamics of the domain wall structure. These findings offer crucial insights into the fundamental mechanisms governing the dynamics of domain walls and other magnetic systems, such as magnetic skyrmions, bubbles, and magnetic vortices, paving the way for development and optimization of next-generation spintronic devices.

Collective dynamics↗
DOE OSTI2023

Analysis and Compensation of Tx and Rx IQ Imbalances in Uplink IM-OFDMA Systems

Index modulation-based orthogonal division multiple access (IM-OFDMA) has recently been proposed as a potential technique for future wireless communication systems due to its superior spectral efficiency and error performance over conventional multiple access schemes. However, its performance is still under investigation by researchers in light of a variety of scenarios and assumptions. Following this direction, in this Article, the individual and joint effects of the transmitter (Tx) and receiver (Rx) in-phase and quadrature imbalances (IQI) on the bit-error-rate (BER) performance of the uplink IM-OFDMA scheme are investigated by considering whether each user's RF front-end is identical or not. Moreover, to reduce the detrimental effect of IQI, a preamble-based estimation and compensation method is proposed for IM-OFDMA. Closed-form expressions for the average BER of IM-OFDMA are obtained by considering the physical effect of the Tx and Rx IQI. Also, using Monte Carlo simulations, the derived expressions are verified under different system configurations. Furthermore, analytical and simulation results reveal that Tx and Rx IQI cause an error floor in the BER performance of IM-OFDMA schemes. However, the proposed estimation and compensation method not only reduces the impact of IQI but also leads to a better BER performance compared to the case of no IQI, which mainly refers to the frequency diversity caused by the IQI.

42 ENGINEERING↗
DOE OSTI2026

Genetic and metabolic drivers of membrane remodeling in Clostridium thermocellum under alcohol stress

Clostridium thermocellum is a leading candidate for consolidated bioprocessing of lignocellulosic biomass into biofuels due to its native cellulolytic capabilities. Beyond ethanol, C. thermocellum is being developed as a platform for producing higher-chain alcohols such as isobutanol and n-butanol. However, its physiological adaptations to alcohol stress remain poorly understood. Here, we investigate how C. thermocellum remodels its membrane lipid composition in response to exogenous ethanol, n-butanol, isobutanol, and butyrate. Exposure to linear alcohols such as n-butanol or to organic acids like butyrate increased the proportion of straight-chain fatty acids in the membrane at the expense of branched-chain species, whereas exposure to the branched alcohol isobutanol produced the opposite effect. Isotope tracer experiments demonstrated that C. thermocellum directly incorporates the carbon backbones of exogenous alcohols and acids into fatty acids, providing a mechanistic basis for these contrasting shifts. We show that the bifunctional aldehyde/alcohol dehydrogenase AdhE is essential for the assimilation of exogenous alcohols into fatty acids, acting through its oxidative activity by first oxidizing alcohols to aldehydes and then converting them to acyl-CoA intermediates. Deletion of the pyruvate:ferredoxin oxidoreductase isozyme pfor4 abolished branched-chain fatty acid synthesis, but supplementation with isobutanol restored production, indicating that Pfor4 substitutes for the canonical branched-chain α-keto acid dehydrogenase complex. These findings reveal two distinct routes for branched-chain fatty acid production in C. thermocellum: a Pfor4-dependent pathway from α-keto acid intermediates derived from amino acid synthesis, and an AdhE-dependent salvage pathway that assimilates exogenous branched-chain alcohols.

Acetivibrio thermocellus↗
DOE OSTI2021

Precision measurement of the electron energy-loss function in tritium and deuterium gas for the KATRIN experiment

Abstract The KATRIN experiment is designed for a direct and model-independent determination of the effective electron anti-neutrino mass via a high-precision measurement of the tritium $$\upbeta $$ β -decay endpoint region with a sensitivity on $$m_\nu $$ m ν of 0.2 $$\hbox {eV}/\hbox {c}^2$$ eV / c 2 (90% CL). For this purpose, the $$\upbeta $$ β -electrons from a high-luminosity windowless gaseous tritium source traversing an electrostatic retarding spectrometer are counted to obtain an integral spectrum around the endpoint energy of 18.6 keV. A dominant systematic effect of the response of the experimental setup is the energy loss of $$\upbeta $$ β -electrons from elastic and inelastic scattering off tritium molecules within the source. We determined the energy-loss function in-situ with a pulsed angular-selective and monoenergetic photoelectron source at various tritium-source densities. The data was recorded in integral and differential modes; the latter was achieved by using a novel time-of-flight technique. We developed a semi-empirical parametrization for the energy-loss function for the scattering of 18.6-keV electrons from hydrogen isotopologs. This model was fit to measurement data with a 95% $$\hbox {T}_2$$ T 2 gas mixture at 30 K, as used in the first KATRIN neutrino-mass analyses, as well as a $$\hbox {D}_2$$ D 2 gas mixture of 96% purity used in KATRIN commissioning runs. The achieved precision on the energy-loss function has abated the corresponding uncertainty of $$\sigma (m_\nu ^2)< {{10}^{-2}}{\hbox {eV}^{2}}$$ σ ( m ν 2 ) < 10 - 2 eV 2 [1] in the KATRIN neutrino-mass measurement to a subdominant level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Development of Emergency Vehicle Preemption Strategies on Smart Corridors in a Digital Twin Environment

Emergency Response Vehicles (ERVs), such as firetrucks, ambulances, etc. operate with the purpose of saving lives and mitigating property damage. As such, ERV travel-time reductions may result in significant benefits to the community. A common strategy to improve travel times is Emergency Vehicle Preemption (EVP). EVP seeks to reduce ERV delays by providing the right-of-way to ERVs as they approach an intersection. This study proposes a new Dynamic Preemption strategy that determines the need for preemption prior to the ERV reaching the vicinity of the intersection, utilizing real-time data streams. This paper evaluates the effectiveness of some existing and proposed preempt control strategies using a digital twin testbed consisting of a series of signalized intersections on an urban arterial in Georgia. The best EVP strategy maximizes the improvement in ERV travel time while minimizing the adverse effect of preemption on the traffic in conflicting directions. Therefore, this study evaluates both the positive impact of EVP on the ERV as well as the adverse impact on the cross-street traffic. The study found that the potential exists for significant improvements in ERV travel time with the proposed Dynamic Preemption strategy, with minimal impact to the conflicting traffic. For the simulation corridor there was a 20% reduction in the ERV travel times with the implementation of the Dynamic Preemption strategy, compared to traditional EVP practices.

Roy, Somdut [Georgia Institute of Technology]↗
DOE OSTI2025

RingX: Scalable Parallel Attention for Long-Context Learning on HPC

The attention mechanism has become foundational for remarkable AI breakthroughs since the introduction of the Transformer, driving the demand for increasingly longer context to power frontier models such as large-scale reasoning language models and high-resolution image/video generators. However, its quadratic computational and memory complexities present substantial challenges. Current state-of-the-art parallel attention methods, such as ring attention, are widely adopted for long-context training but utilize a point-to-point communication strategy that fails to fully exploit the capabilities of modern HPC network architectures. In this work, we propose ringX, a scalable family of parallel attention methods optimized explicitly for HPC systems. By enhancing workload partitioning, refining communication patterns, and improving load balancing, ringX achieves up to 3.4 × speedup compared to conventional ring attention on the Frontier supercomputer. Optimized for both bi-directional and causal attention mechanisms, ringX demonstrates its effectiveness through training benchmarks of a Vision Transformer (ViT) on a climate dataset and a Generative Pre-Trained Transformer (GPT) model, Llama3 8B. Our method attains an end-to-end training speedup of approximately 1.5 × in both scenarios. To our knowledge, the achieved 38% model FLOPs utilization (MFU) for training Llama3 8B with a 1M-token sequence length on 4,096 GPUs represents one of the highest training efficiencies reported for long-context learning on HPC systems. Our code implementation is available at https://github.com/jqyin/ringX-attention.

Yin, Junqi [ORNL] (ORCID:0000000338435520)↗
DOE OSTI2020

Improved recovery of cell-cycle gene expression in Saccharomyces cerevisiae from regulatory interactions in multiple omics data

Gene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task. Results We employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression in Saccharomyces cerevisiae . Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators ( STE12-TEC1 and RAP1-HAP1-MSN4 ). Conclusion Our findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Ultrafast optical modulation of vibrational strong coupling in ReCl(CO) 3 (2,2-bipyridine)

Polaritons – hybrid light-matter states formed from the strong coupling of a bright molecular transition with a confined photonic mode – may offer new opportunities for optical control of molecular behavior. Vibrational strong coupling (VSC) has been reported to impact ground-state chemical reactivity, but its influence on electronic excited-state dynamics remains unexplored. Here, we take a first step towards excited-state VSC by demonstrating optical modulation of the ReCl(CO) 3 (bpy), (bpy = 2,2-bipyridine) complex under VSC using femtosecond ultraviolet (UV)-pump/infrared (IR)-probe spectroscopy. We establish ground-state VSC of ReCl(CO) 3 (bpy) in a microfluidic Fabry-Pérot cavity equipped with indium tin oxide (ITO)-coated mirrors. ITO is effectively dichroic as it is reflective in the IR and transmissive in the UV-visible and therefore minimizes optical interference. Excitation with UV pump light drives ReCl(CO) 3 (bpy) into a manifold of electronic excited states that subsequently undergo non-radiative relaxation dynamics. We probe the transient response of the strongly-coupled system in the mid-IR, observing both Rabi contraction and cavity-filtered excited-state absorption signatures. We reconstruct the intrinsic response of intracavity molecules from the transient cavity transmission spectra to enable quantitative comparison with extracavity control experiments. We report no changes in the excited-state dynamics of ReCl(CO) 3 (bpy) under ground-state VSC. However, we do observe significant amplification of transient vibrational signals due to classical cavity-enhanced optical effects. This effort lays the groundwork to pursue direct excited-state VSC aimed at modulating photochemical reactivity.

cavity-enhanced spectroscopy↗
DOE OSTI2020

Modeling of Atom Interferometer Accelerometer

This report presents the theoretical effort to model and simulate the atom-interferometer accelerometer operating in a highly mobile environment. Multitudes of non-idealities may occur in such a rapidly-changing environment with a large acceleration whose amplitude and direction both change quickly. We studied the undesired effect of high mobility in the atom-interferometer accelerator in a detailed model and a simulator. The undesired effects include the atom cloud's movement during Raman pulses, the Doppler effect due to the relative movement between the atom-cloud and the supporting platform, the finite atom cloud temperature, and the lateral movement of the atom cloud. We present the relevant feed-forward mitigation strategies for each identified non-ideality to neutralize the impact and obtain accurate acceleration measurements.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

I.1 Heavy-Duty Diesel Combustion (Sandia National Laboratories)

Regulatory drivers and market demands for lower pollutant emissions, lower carbon dioxide emissions, and lower fuel consumption motivate the development of cleaner and more fuel-efficient engine operating strategies. Most current production heavy-duty diesel engines use a combination of both in-cylinder and exhaust emissions-control strategies to achieve these goals. The emissions and efficiency performance of in-cylinder strategies depend strongly on flow and mixing processes that can be influenced by using multiple fuel injections. Past work performed under this project showed that adding a second injection can reduce soot to levels below what would have been produced by an unchanged first injection, thereby increasing load while decreasing soot and potentially reducing brake specific fuel consumption. Information characterizing the important in-cylinder processes with multiple injections has been gleaned from ensemble-averaged planar laser-induced incandescence (PLII) imaging visualizing the soot cloud and planar induced fluorescence (PLIF) of OH characterizing the soot oxidation regions. PLII showed a consistent disruption of the first injection soot cloud by the second injection. In conjunction with OH-PLIF, differences in soot oxidation patterns for multiple injections compared to single injections were observed. This understanding was further enhanced in FY20, when high-speed imaging resolving the above-mentioned effects in a single cycle were combined with direct numerical simulations investigating the multiple-injection ignition process on the microscopic level of turbulence and chemistry interaction. In FY21, these findings in conjunction with findings from other researchers published in the scientific literature were composed into a preliminary multiple-injection conceptual model of fuel-mixing, injection and ignition processes. Remaining key research questions were also highlighted. In addition, wall heat flux was investigated experimentally and with numerical simulations to understand the potential of multiple injections to reduce the engine heat losses and further enhance the efficiency.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Improved Understanding of Zircaloy Hydrogen Pickup Mechanisms in BWRs. NSUF Report (Final Report)

Using an instrumented autoclave operating at 288°C and 1100 psi, electrical impedance spectroscopy (EIS) measurements were taken from three neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). The three high priority fuel channel samples were tested at 288°C using a simulated PWR water with the AC frequency range extending down through 10 -4 Hz, with corresponding DC measurements to validate the resistance measurements. EIS measurements were also taken at room temperature using a simulated PWR water from four neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). and four water rod samples (Zircaloy-2). EIS measurements were also performed on unirradiated Zircaloy-2 tubes that were previously oxidized in autoclaves at various temperatures for different time durations to obtain oxide layers of different thicknesses. These experiments proved to be more challenging than anticipated because of the small size of the samples, which required development of fixtures to allow valid testing of the materials, some of which are radioactive, under pressure and temperature. Equivalent circuits were used to model the data. Some samples resulted in Nyquist plots where the arcs were too close to each other resulting in the inability to separate the different oxide morphological thicknesses. Microstructural characterization in the form of optical microscopy, scanning electron microscopy and (scanning) transmission electron microscopy was performed to document the microstructural features of the oxidized Zr-2 and Zr-4 samples (oxide and metal/oxide interface) and to interpret the EIS results. At the start of this work, it was assumed that despite the laminar cracks and porosity visible in both optical and SEM, the oxide scale is dense near the metal, effectively preventing the water from contacting the metal directly. While it is generally accepted that the large laminar cracks pores allow the water to percolate into the oxide film, it is unknown just how deep it can go, and how interconnected is the porosity, which will dictate how much the porosity is connected electrically to the water within the oxide scale. The EIS results demonstrated that there is in fact a distribution of oxide thicknesses over the various samples, and within each material condition, that are communicating to the liquid. This is derived from the fact that Nyquist plots have an arc radius whose center point is often below the Z’ axis, implying that instead of a single oxide thickness, these oxides have a distribution of thicknesses. This distribution of thicknesses produces arcs that are not clearly separable in most of the samples.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Adaptive Angular Quadrature Scheme for a backwards-in-time Method of Characteristics Solution to the Radiative Transfer Equation [Slides]

Radiative transfer/radiation transport are important problems to solve in astrophysics and high energy density physics. Various methods exist to solve radiation transport, such as Monte Carlo (MC), Discrete Ordinates (S N ), Method of Characteristics (MOC), and the spherical harmonics (P N ) method. Method of Characteristics requires “launching” of rays in discrete directions. Unresolved details of angular mesh create ray effects and can miss sources in the domain. Ray effects can lead to unphysical “stepping” in solution and incorrect energy deposition. Adaptive quadrature schemes can be used to detect and mitigate these effects. The Method of Characteristics (MOC) is a common method for solving hyperbolic PDEs in radiation transport and supersonic flow problems. Generally in MOC for radiation transport, virtual particles are tracked from birth to the end of a timestep. This requires interpolation to go from final location to cell averaged or corner values of angular intensity. Backwards-in-Time (BIT) particle tracking avoids this by prescribing the final position of the virtual particle at the cell nodes/corners. Angular intensities are computed at time k + 1 by launching ray back to previous timestep(s), or t = 0. Scheme allows solution to be computed as the characteristic ray is traced backwards in time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Estimating Leaf Area Index in Row Crops Using Wheel-Based and Airborne Discrete Return Light Detection and Ranging Data

Leaf area index (LAI) is an important variable for characterizing plant canopy in crop models. It is traditionally defined as the total one-sided leaf area per unit ground area and is estimated by both direct and indirect methods. This paper explores the effectiveness of using light detection and ranging (LiDAR) data to estimate LAI for sorghum and maize with different treatments at multiple times during the growing season from both a wheeled vehicle and Unmanned Aerial Vehicles. Linear and nonlinear regression models are investigated for prediction utilizing statistical and plant structure-based features extracted from the LiDAR point cloud data with ground reference obtained from an in-field plant canopy analyzer (indirect method). Results based on the value of the coefficient of determination ( R 2 ) and root mean squared error for predictive models ranged from ∼0.4 in the early season to ∼0.6 for sorghum and ∼0.5 to 0.80 for maize from 40 Days after Sowing to harvest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Wind Turbine Gust Load Alleviation with Active Flow Control

Large wind turbine rotors are becoming more common in utility-scale wind power, especially for offshore wind plants. However, the trend toward large rotors can be limited by their ability to manage dynamic and extreme loads. To provide a safety margin for the rotor design and avoid catastrophic events such as tower strikes, extreme loads need to be controlled. The objective of this study is to develop and evaluate a feedback control system to alleviate extreme loads and reduce blade deflections under gust events using active flow control devices. We also propose a modification in the turbine controller to achieve further reduction in extreme loads. The extreme load reductions are evaluated under gust wind conditions with direction changes according to the IEC standard. The effects of the gust alleviation controller on turbine performance and fatigue loads are investigated as well. With the deployment of the gust alleviation controller and modified turbine baseline controller, the extreme loads and deflections reduce by up to 23%. The energy captured by the turbine is not affected by the proposed gust alleviation controller. The fatigue loads of various wind turbine components are either reduced or remain unchanged.

17 WIND ENERGY↗
DOE OSTI2022

Curled-Skewed Wakes behind Yawed Wind Turbines Subject to Veered Inflow

This work presents a new engineering analytical model that predicts the effect of both the turbine yaw misalignment and the inflow wind veer on the wake flow distribution downwind of a wind turbine. To consider the veered inflow, two methods were examined. In the first method, the curled shape of the wake due to the yaw offset is initially modelled. The wake shape is then laterally skewed at each height due to the wind veer based on the assumption that the turbine wake is transported downstream by the incoming flow. The second method is a more realistic approach that accounts for the effect of wind veer on the wind velocity direction and the yaw angle seen by the wind turbine. This models the wake region in a local coordinate system defined based on the wind direction at each height. A coordinate transformation is then performed to represent the wake flow distribution in the global coordinate system attached to the ground. The results show that while the two methods provide similar outputs for small variations in the wind direction across the rotor, the difference becomes more evident with an increase in wind veer. High-fidelity simulations for a turbine subject to a neutral atmospheric boundary layer were employed to validate model predictions for different operating conditions.

17 WIND ENERGY↗
DOE OSTI2021

Constraining Mornings and Evenings on Distant Worlds: A new Semianalytical Approach and Prospects with Transmission Spectroscopy

The technique of transmission spectroscopy—the variation of a planetary radius with wavelength due to opacity sources in the planet’s terminator region—has been one of the most successful in the characterization of exoplanet atmospheres to date, providing key insights into the composition and structure of these distant worlds. A common assumption made when using this technique, however, is that the variations are the same in the entire terminator region. In reality, the morning and evening terminators might have distinct temperature, pressure, and thus compositional profiles due to the inherent 3D nature of the planet, which would, in turn, give rise to different spectra on each side of it. Constraining those might be fundamental for our understanding of not only the weather patterns in these distant worlds but also the planetary formation signatures that might only be possible to extract once these features are well understood. Motivated by this physical picture, in this work, we perform a detailed study of the observational prospects of detecting this effect. We present an open-source semianalytical framework with which this information can be extracted directly from transit light curves and perform a detailed study of the prospects of detecting the effect with current missions, such as TESS, and upcoming ones, such as JWST. Our results show that these missions show great promise for the detection of this effect. Transmission spectroscopy studies with JWST, in particular, could provide spectra of each of the limbs, allowing us to convey 3D information previously accessible only via phase curves.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗