Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SIMULATION”

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,297 records · Page 72

Assessment of dynamic-screw-pinch-driven, current-scaled MagLIF target implosion performance using 3D magnetohydrodynamic simulations

Analytic studies and two-dimensional “clean” radiation-magnetohydrodynamic (rad-MHD) simulations employing dynamical similarity driver-target scaling prescriptions [Ruiz et al., Phys. Plasmas 30, 032708 (2023)] suggest that Magnetized Liner Inertial Fusion (MagLIF) target implosions can scale to > 10 MJ DT fusion yields when peak drive current is increased beyond 60 MA. We present results from three-dimensional (3D) rad-MHD simulations of similarity-scaled MagLIF target implosions at peak drive currents ranging from 15 to 40 MA. Simulations in this study suggest that magneto-Rayleigh–Taylor instability (MRTI) growth and feedthrough to the fuel region are more severe at higher drive current scales, which reduces the fusion yield compared to prior analytic and 2D clean simulation predictions. In contrast to standard MagLIF, simulations of current-scaled MagLIF target implosions driven by a dynamic screw pinch (DSP) demonstrate reduced MRTI feedthrough and greater fuel magnetization, resulting in improved thermonuclear performance and enhanced performance scaling with peak drive current. DSP drive enables additional scaling of the liner mass to increase liner radius but maintain implosion time, resulting in higher implosion velocities at the expense of increased susceptibility to MRTI. We present a current- and mass-scaled simulated DSP-MagLIF target implosion at the ∼ 40 MA peak current level that produces ignition scale performance, demonstrating a burn-averaged Lawson ignition parameter above unity and DT fusion yield above 1 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Differentiable Simulation of a Liquid Argon Time Projection Chamber

Liquid argon time projection chambers (LArTPCs) are widely used in particle detection for their tracking and calorimetric capabilities. The particle physics community actively builds and improves high-quality simulators for such detectors in order to develop physics analyses in a realistic setting. The ability of these simulators to mimic real, measured data is limited by the modeling of the physical detectors used for data collection. This modeling can be improved by performing dedicated calibration measurements. Conventional approaches calibrate individual detector parameters or processes one at a time. However, the impact of detector processes is entangled, making this a poor description of the underlying physics. We introduce a differentiable simulator that enables a gradient-based optimization, allowing for the first time a simultaneous calibration of all detector parameters. We describe the procedure of making a differentiable simulator, highlighting the challenges of retaining the physics quality of the standard, non-differentiable version while providing meaningful gradient information. We further discuss the advantages and drawbacks of using our differentiable simulator for calibration. Finally, we provide a starting point for extensions to our approach, including applications of the differentiable simulator to physics analysis pipelines.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Data and scripts associated with the manuscript "A Hydrogeophysical Framework to Assess Infiltration During a Simulated Ecosystem-Scale Flooding Experiment”

This package contains the data and scripts used in "A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment" (Adebayo et al., 2023). The Geophysics folder contains two subfolders containing Ground Penetrating Radar (GPR) and Electrical Resistivity Imaging (ERI) datasets acquired between August 23 and August 26, 2021. The GPR subfolder contains raw data (Tempest_Project.GPZ) and processed results using the Sensors and Software EKKO Project software. The ERI subfolder contains two additional subfolders named Static and Timelapse. The Static subfolder contains background resistivity datasets acquired before infiltration, both raw and inverted results using the ResIpy code. The Timelapse subfolder contains repeated resistivity datasets acquired during infiltration on two profiles, along with raw and inverted results using the ResIpy code. The Flow model folder contains ATS model simulation files. The Geophysics-informed subfolder contains simulation datasets for two profiles generated using information from the geophysics data. The National_database-informed subfolder contains simulation datasets for two profiles generated using information from a national soil database. Each subfolder contains initial simulation results, multiple ensembles, and final simulation files. The Sensitivity_analysis folder contains the observation files exported from the flow model, notebooks used in running the code, and sensitivity analysis results showing the sensitivity of different parameters over a 48-hour period. We used mutual information-based global sensitivity analysis for the sensitivity analysis, following Jiang et al., 2022. This study shows a novel approach that integrates geophysical data with hydrological modeling to understand subsurface flow during coastal flooding.

54 ENVIRONMENTAL SCIENCES↗

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗

ADDGALS: Simulated Sky Catalogs for Wide Field Galaxy Surveys

Abstract We present a method for creating simulated galaxy catalogs with realistic galaxy luminosities, broadband colors, and projected clustering over large cosmic volumes. The technique, denoted Addgals (Adding Density Dependent GAlaxies to Lightcone Simulations), uses an empirical approach to place galaxies within lightcone outputs of cosmological simulations. It can be applied to significantly lower-resolution simulations than those required for commonly used methods such as halo occupation distributions, subhalo abundance matching, and semi-analytic models, while still accurately reproducing projected galaxy clustering statistics down to scales of r ∼ 100 h −1 kpc . We show that Addgals catalogs reproduce several statistical properties of the galaxy distribution as measured by the Sloan Digital Sky Survey (SDSS) main galaxy sample, including galaxy number densities, observed magnitude and color distributions, as well as luminosity- and color-dependent clustering. We also compare to cluster–galaxy cross correlations, where we find significant discrepancies with measurements from SDSS that are likely linked to artificial subhalo disruption in the simulations. Applications of this model to simulations of deep wide-area photometric surveys, including modeling weak-lensing statistics, photometric redshifts, and galaxy cluster finding, are presented in DeRose et al., and an application to a full cosmology analysis of Dark Energy Survey (DES) Year 3 like data is presented in DeRose et al. We plan to publicly release a 10,313 square degree catalog constructed using Addgals with magnitudes appropriate for several existing and planned surveys, including SDSS, DES, VISTA, Wide-field Infrared Survey Explorer, and Rubin Observatory’s Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗

Subsetted model output for aerosol-related variables from a LASSO-CACTI WRF simulation

Aerosol-related model output from a Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m through the use of four nested domains, labeled d1 through d4. The simulations are of convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. This version of the output is a collection of variables subsetted from the raw WRF output to reduce the file overhead for users not needing the full raw dataset. The subset files are in netCDF format with the height-based variables on the raw, model eta levels. The aerosol properties are consistent with the Thompson aerosol-aware cloud microphysics scheme where the simulated aerosol are proxies for representative “water friendly” and “ice friendly” aerosol.

54 ENVIRONMENTAL SCIENCES↗

Subsetted model output for static variables from a LASSO-CACTI WRF simulation

Static model output from a Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m through the use of four nested domains, labeled d1 through d4. The simulations are of convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. This version of the output is a collection of variables subsetted from the raw WRF output to reduce the file overhead for users not needing the full raw dataset. The subset files are in netCDF format with the height-based variables on the raw, model eta and soil levels. The variables in this file are ones that do not change during a simulation, such as the terrain height and map-scale factors.

54 ENVIRONMENTAL SCIENCES↗

Mission simulation as an approach to develop requirements for automation in Advanced Life Support Systems

This paper examines mission simulation as an approach to develop requirements for automation and robotics for Advanced Life Support Systems (ALSS). The focus is on requirements and applications for command and control, control and monitoring, situation assessment and response, diagnosis and recovery, adaptive planning and scheduling, and other automation applications in addition to mechanized equipment and robotics applications to reduce the excessive human labor requirements to operate and maintain an ALSS. Based on principles of systems engineering, an approach is proposed to assess requirements for automation and robotics using mission simulation tools. First, the story of a simulated mission is defined in terms of processes with attendant types of resources needed, including options for use of automation and robotic systems. Next, systems dynamics models are used in simulation to reveal the implications for selected resource allocation schemes in terms of resources required to complete operational tasks. The simulations not only help establish ALSS design criteria, but also may offer guidance to ALSS research efforts by identifying gaps in knowledge about procedures and/or biophysical processes. Simulations of a planned one-year mission with 4 crewmembers in a Human Rated Test Facility are presented as an approach to evaluation of mission feasibility and definition of automation and robotics requirements.

NASA Discipline Number 61-10↗

DSMC Shock Simulation of Saturn Entry Probe Conditions

This work describes the direct simulation Monte Carlo (DSMC) investigation of Saturn entry probe scenarios and the influence of non-equilibrium phenomena on Saturn entry conditions. The DSMC simulations coincide with rarefied hypersonic shock tube experiments of a hydrogen-helium mixture performed in the Electric Arc Shock Tube (EAST) at NASA Ames Research Center. The DSMC simulations are post-processed through the NEQAIR line-by-line radiation code to compare directly to the experimental results. Improved collision cross-sections, inelastic collision parameters, and reaction rates are determined for a high temperature DSMC simulation of a 7-species H2-He mixture and an electronic excitation model is implemented in the DSMC code. Simulation results for 27.8 and 27.4 kms shock waves are obtained at 0.2 and 0.1 Torr respectively and compared to measured spectra in the VUV, UV, visible, and IR ranges. These results confirm the persistence of non-equilibrium for several centimeters behind the shock and the diffusion of atomic hydrogen upstream of the shock wave. Although the magnitude of the radiance did not match experiments and an ionization inductance period was not observed in the simulations, the discrepancies indicated where improvements are needed in the DSMC and NEQAIR models.

Modeling↗

Between-Simulator Comparison of Motion Filter Order and Break Frequency Interaction Effects on Pilot Control Behavior

This paper investigates the interaction effects of motion filter order and break frequency on pilots' manual control behavior and control performance using two simulators. Eighteen pilots performed the experiment in the Vertical Motion Simulator (VMS) at NASA Ames Research Center and twenty pilots in the SIMONA Research Simulator at Delft University of Technology. The experiment used a full-factorial design with three motion filter orders (first-, second-, and third-order) and two filter break frequencies (0.5 and 2.0 rad/s), in addition to reference no-motion and full-motion conditions. Key task variables, such as the quality of the motion and visual cues and the characteristics of the sidestick, were matched across both simulators. Overall, the expected effects of filter order and break frequency variations were found, with both increasing order and increasing break frequency causing pilots to use less motion feedback in their control strategy, resulting in lower tracking performance. Furthermore, across the wide range of filter orders tested in the experiment, the existing Sinacori-Schroeder motion fidelity criterion was found to be a good predictor of the interaction effects of both filter settings on pilot control behavior. For the same motion condition, there was a consistent offset in the results between simulators, due to the more high-gain control strategy adopted by a number of the VMS pilots. Still, the observed relative trends in pilot control behavior and performance between motion conditions were equivalent in both simulators and thus accurately replicated.

manual control↗

The Earth Model Column Collaboratory (EMC2) v1.1: An Open-Source Ground-Based Lidar and Radar Instrument Simulator and Subcolumn Generator for Large-Scale Models

Climate models are essential for our comprehensive understanding of Earth's atmosphere and can provide critical insights on future changes decades ahead. Because of these critical roles, today's climate models are continuously being developed and evaluated using constraining observations and measurements obtained by satellites, airborne, and ground-based instruments. Instrument simulators can provide a bridge between the measured or retrieved quantities and their sampling in models and field observations while considering instrument sensitivity limitations. Here we present the Earth Model Column Collaboratory (EMC2), an open-source ground-based lidar and radar instrument simulator and subcolumn generator, specifically designed for large-scale models, in particular climate models, but also applicable to high-resolution model output. EMC2 provides a flexible framework enabling direct comparison of model output with ground-based observations, including generation of subcolumns that may statistically represent finer model spatial resolutions. In addition, EMC2 emulates ground-based (and air- or space-borne) measurements while remaining faithful to large-scale models' physical assumptions implemented in their cloud or radiation schemes. The simulator uses either single particle or bulk particle size distribution lookup tables, depending on the selected scheme approach, to perform the forward calculations. To facilitate model evaluation, EMC2 also includes three hydrometeor classification methods, namely, radar- and sounding-based cloud and precipitation detection and classification, lidar-based phase classification, and a Cloud Feedback Model Intercomparison Project Observational Simulator Package (COSP) lidar simulator emulator. The software is written in Python, is easy to use, and can be straightforwardly customized for different models, radars, and lidars. Following the description of the logic, functionality, features, and software structure of EMC2, we present a case study of highly supercooled mixed-phase cloud based on measurements from the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE). We compare observations with the application of EMC2 to outputs from four configurations of the NASA Goddard Institute for Space Studies (GISS) climate model (ModelE3) in single-column model (SCM) mode and from a large-eddy simulation (LES) model. We show that two of the four ModelE3 configurations can form and maintain highly supercooled precipitating cloud for several hours, consistent with observations and LES. While our focus is on one of these ModelE3 configurations, which performed slightly better in this case study, both of these configurations and the LES results post-processed with EMC2 generally provide reasonable agreement with observed lidar and radar variables. As briefly demonstrated here, EMC2 can provide a lightweight and flexible framework for comparing the results of both large-scale and high-resolution models directly with observations, with relatively little overhead and multiple options for achieving consistency with model microphysical or radiation scheme physics.

Earth Model Column Collaboratory↗

On-orbit/On-surface Servicing, Assembly, and Manufacturing (OSAM) Architecture Simulation System (OASiS)

On-orbit/on-surface servicing, assembly, and manufacturing (OSAM) –recently renamed“in-space servicing, assembly, and manufacturing”(ISAM) –will aid in the development of vital research-enabling structuresin space, such as the In-space Assembled Telescope, the Lunar Gateway, and the Lunar Safe Haven. As these structures become larger and more complex, developers will need new methods for verifying and validating their OSAM mission architectures at full-scale prior to launch, without paying for extensive facility upgrades or additional hardware prototypes. Under the OSAM Architecture Simulation System (OASiS) project, modular OSAM modeling and simulation capabilities were developed to help rapidly prototype and evaluate OSAM technologies and operational concepts.Theywerethenusedto demonstrate and evaluate data from the simulation of a simpleOSAM operation against ananaloghardware setup. Using collected test data, the teamevaluatedthe accuracy of the simulation environment against the behavior of real hardwareandgainedinsight on how to further improve OSAM modeling and simulation. With thesenewcapabilities and many more to come, developers will be able to test complex OSAM operations on multiple scales, iterate on technology designs faster, achieve higher operational accuracy, and ultimately reduce mission development costs. A comprehensive, high-fidelity simulation environment will ensure the successful servicing, assembly, and manufacturing of structures that help us study and explore our universe.

isam↗

Validation of (not-historical) large-event near-fault ground-motion simulations for use in civil engineering applications

Ground-motion simulations generated from physics-based wave propagation models are gaining increasing interest in the engineering community for their potential to inform the performance-based design and assessment of infrastructure residing in active seismic areas. A key prerequisite before the ground-motion simulations can be used with confidence for application in engineering domains is their comprehensive and rigorous investigation and validation. This article provides a four-step methodology and acceptance criteria to assess the reliability of simulated ground motions of not historical events, which includes (1) the selection of a population of real records consistent with the simulated scenarios, (2) the comparison of the distribution of Intensity Measures (IMs) from the simulated records, real records, and Ground-Motion Prediction Equations (GMPEs), (3) the comparison of the distribution of simple proxies for building response, and (4) the comparison of the distribution of Engineering Demand Parameters (EDPs) for a realistic model of a structure. Specific focus is laid on near-field ground motions (<10km) from large earthquakes (M w 7), for which the database of real records for potential use in engineering applications is severely limited. The methodology is demonstrated through comparison of (2490) near-field synthetic records with 5 Hz resolution generated from the Pitarka et al (2019) kinematic rupture model with a population of (38) pulse-like near-field real records from multiple events and, when applicable, with NGA-W2 GMPEs. Lastly, the proposed procedure provides an effective method for informing and advancing the science needed to generate realistic ground-motion simulations, and for building confidence in their use in engineering domains.

42 ENGINEERING↗

Ingrained: An Automated Framework for Fusing Atomic-Scale Image Simulations into Experiments

To fully leverage the power of image simulation to corroborate and explain patterns and structures in atomic resolution microscopy, an initial correspondence between the simulation and experimental image must be established at the outset of further high accuracy simulations or calculations. Furthermore, if simulation is to be used in context of highly automated processes or high-throughput optimization, the process of finding this correspondence itself must be automated. In this work, "ingrained," an open-source automation framework which solves for this correspondence and fuses atomic resolution image simulations into the experimental images to which they correspond, is introduced. Here, the overall "ingrained" workflow, focusing on its application to interface structure approximations, and the development of an experimentally rationalized forward model for scanning tunneling microscopy simulation are described.

36 MATERIALS SCIENCE↗

Realistic large eddy and dispersion simulation experiments during project sagebrush phase 1

Typical dispersion simulations use mesoscale models at grid spacing greater than 1 km to simulate pollutant transport. Advances in computing power have allowed Large Eddy Simulations (LES), (grid spacing less than 100 m) to be run with boundary conditions in the grayscale regime (grid size of 1 km–100 m). We used the WRF mesoscale/LES model and the HYSPLIT particle dispersion model to simulate SF 6 transport from Project Sagebrush Phase 1, using a grayscale-aware scheme with 4 nests of 3 km, 1 km, 200 m, 67 m. The WRF/LES simulation was in satisfactory agreement with meteorological observations but the HYSPLIT SF6 simulations indicated insufficient sub-grid turbulence kinetic energy, especially in the lowest 100 m. This resulted in over predictions of the surface SF6 concentration and insufficient vertical transport. Best agreement with the Sagebrush observations was obtained with the WRF/LES large-scale and sub-grid turbulence, enhanced with a HYSPLIT parameterized turbulence profile.

54 ENVIRONMENTAL SCIENCES↗

GeoMicro3D: A novel numerical model for simulating the reaction process and microstructure formation of alkali-activated slag

Highlights: • The GeoMicro3D model was innovatively developed for simulating the reaction process and microstructure formation of AAS. • In GeoMicro3D, the dissolution of slag and reactions of aqueous ions are fundamentally described in the local lattice cell. • A novel approach is proposed to speed up the nucleation simulations. • GeoMicro3D was implemented and verified with the relevant experimental data and thermodynamic calculation results. For the first time, this study developed a novel model, named GeoMicro3D, to simulate the reaction process and microstructure formation of alkali-activated slag. The GeoMicro3D model consists of four modules that are designed to simulate, respectively: (i) the initial spatial distribution of real-shape slag particles in alkaline activator, (ii) the dissolution of slag and diffusion of ions via the transition state theory and lattice Boltzmann method, respectively, (iii) the spatial distribution of reaction products using a nucleation probability theory, and (iv) the chemical reactions with thermodynamic modelling. Afterwards the GeoMicro3D model was implemented and verified. The simulation results were discussed and compared with the relevant experimental data and thermodynamic calculation results using GEMS. A good agreement was found in the comparisons, showing the strong simulation capability of GeoMicro3D.

36 MATERIALS SCIENCE↗

Hermes-3: Multi-component plasma simulations with BOUT++

A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. Hermes-3 is built on plasma simulation framework BOUT++, consolidating earlier SD1D and Hermes models into a single code that can be configured at run-time to solve plasma models in 1D, 2D or 3D, either for transport (steady-state) or turbulent (time-evolving) problems, with an arbitrary number of ion and neutral species. Here, we describe the improved numerical algorithms and software design that have been implemented in Hermes-3. To demonstrate the capabilities of this tool, applications relevant to the boundary of tokamak plasmas are presented: 1D simulations of diveror plasmas evolving equations for all charge states of neon and deuterium; 2D transport simulations of tokamak equilibria in single-null X-point geometry with plasma ion and neutral atom species; and simulations of the time-dependent propagation of plasma filaments (blobs). Hermes-3 is publicly available on Github under the GPL-3 open source license. The repository includes documentation and a suite of unit, integrated and convergence tests.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surrogate-driven design optimization with uncertainty constraints in Monte Carlo simulations

In multi-objective design tasks, the computational cost increases rapidly when high-fidelity simulations are used to evaluate objective functions. Surrogate models help mitigate this cost by approximating the simulation output, simplifying the design process. However, under high uncertainty, surrogate models trained on noisy data can produce inaccurate predictions, as their performance depends heavily on the quality of training data. This study investigates the impact of data uncertainty on two multi-objective design problems modelled using Monte Carlo transport simulations: a neutron moderator and an ion-to-neutron converter. For each, a grid search was performed using five different tally uncertainty levels to generate training data for neural network surrogate models. These models were then optimized using NSGA-III. The recovered Pareto-fronts were analyzed across uncertainty levels: in the moderator problem, normalized hypervolume dropped from 0.886 at 1.0% uncertainty to 0.748 at 10% uncertainty, while in the converter problem it remained near 0.50 for all cases. Average simulation times were also compared to evaluate the trade-off between accuracy and computational cost. Results show that the influence of simulation uncertainty is strongly problem-dependent. In the neutron moderator case, higher uncertainties led to exaggerated objective sensitivities and distorted Pareto-fronts, reducing normalized hypervolume. In contrast, the ion-to-neutron converter task was less affected—low-fidelity simulations produced results similar to those from high-fidelity data. These findings suggest that a fixed-fidelity approach is not optimal. Surrogate models can recover the Pareto-front under noisy conditions, and multi-fidelity studies help identify suitable uncertainty levels for each problem to balance efficiency and accuracy.

07 ISOTOPE AND RADIATION SOURCES↗