Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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,315 records · Page 73

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

SERPENT 2 and TRIPOLI-4 neutronics benchmarking of the Jules Horowitz reactor

Jules Horowitz Reactor (JHR) is a future center for nuclear material testing in Europe, currently under construction at CEA Cadarache in southern France. It has been modelled extensively with CEA's HORUS3D/N neutron calculation package, that includes a Monte-Carlo neutron transport code TRIPOLI-4. In order to support the future computational analyses on the reactor core physics, VTT's Monte Carlo neutronics code SERPENT 2 has been used to create an alternative model of the JHR, taking advantage of the pre-existing Constructive Solid Geometry (CSG) structures on which the TRIPOLI-4 model was built. This work describes a benchmark, in which two JHR core configurations were analyzed using both of these models. The first configuration is a fresh fuel core at the beginning of life, in which aluminum based dummy-devices are installed in the reactor test positions. The second configuration represents the end-of-cycle equilibrium core with a heterogeneous burnup distribution, in which CALIPSO devices containing cladding material samples are installed at the in-core experimental positions, while the reflector positions are occupied by DEN-REP devices containing UO{sub 2} fuel. Neutron flux spectra and heat deposition are the parameters of interest in this study, due to their high importance in creating the desired experimental irradiation environment. The impact of delayed gamma energy deposition will be studied in detail, since the two models account for it with different methods. This work aims to establish a SERPENT 2 model of JHR, understanding its differences and agreements with TRIPOLI-4, and discussing its potential future applications. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Microwave Analysis with Monte Carlo Methods for ECH Transmission Lines

A new code framework, MORAMC, is presented which model transmission line (TL) systems consisting of overmoded circular waveguide and other components including miter bends and transmission line gaps. The transmission line is modeled as a set of mode converters in series where each component is composed of one or more converters. The parametrization of each mode converter can account for the fabrication tolerances of physically realizable components. These tolerances as well as the precision to which these TL systems can be installed and aligned gives a practical limit to which the uncertainty of the microwave performance of the system can be calculated. Because of this, Monte Carlo methods are a natural fit and are employed to calculate the probability distribution that a given TL can deliver a required power and mode purity. Several examples are given to demonstrate the usefulness of MORAMC in optimizing TL systems.

47 OTHER INSTRUMENTATION↗

Genetic algorithm optimization of nuclear criticality experiment for reduction of intermediate-energy 239 Pu nuclear data uncertainties

Nuclear criticality experiments are conducted to investigate specific nuclear data important for safe handling and storage of fissile materials, reactor design and operation, and the validation of radiation transport codes. Incorrect or uncertain nuclear data can prohibitively impact operational safety limits, reactor licensing, and predictive simulation capability; therefore, integral measurements from criticality experiments are necessary and should be performed frequently. To maximize the impact of the integral measurements, it is important to consider experiment geometry, material selection, and component dimensions. When taking these considerations into account, the experiment design process becomes iterative and very time intensive. This work utilizes a genetic algorithm to efficiently explore potential nuclear criticality experiment designs for the Laboratory Directed Research & Development project PARADIGM (PARallel Approach of Differential and InteGral Measurements) at Los Alamos National Laboratory. In this paper, the building blocks of the genetic algorithm are discussed in detail, the genetic algorithm methodology is verified, and the genetic algorithm is used to produce three candidate experiment models for the final PARADIGM design. The three candidate models produced by the genetic algorithm consist of copper-reflected assemblies containing 14 repeating units of alumina, graphite, boron, and plutonium plates. Furthermore, in addition to the optimization results, final design considerations are also discussed for designs with a height and/or weight very close to or slightly above assembly machine operational limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electronic Excitations in Complex Molecular Environments: Many-Body Green’s Functions Theory in VOTCA-XTP

Many-body Green’s functions theory within the GW approximation and the Bethe-Salpeter Equation (BSE) is implemented in the open-source VOTCA-XTP software, aiming at the calculation of electronically excited states in complex molecular environments. Based on Gaussian-type atomic orbitals and making use of resolution of identity techniques, the code is designed specifically for nonperiodic systems. Application to a small molecule reference set successfully validates the methodology and its implementation for a variety of excitation types covering an energy range from 2 to 8 eV in single molecules. Further, embedding each GW-BSE calculation into an atomistically resolved surrounding, typically obtained from Molecular Dynamics, accounts for effects originating from local fields and polarization. Using aqueous DNA as a prototypical system, different levels of electrostatic coupling between the regions in this GW-BSE/MM setup are demonstrated. Particular attention is paid to charge-transfer (CT) excitations in adenine base pairs. It is found that their energy is extremely sensitive to the specific environment and to polarization effects. The calculated redshift of the CT excitation energy compared to a nucelobase dimer treated in vacuum is of the order of 1 eV, which matches expectations from experimental data. Predicted lowest CT energies are below that of a single nucleobase excitation, indicating the possibility of an initial (fast) decay of such an UV excited state into a binucleobase CT exciton. Furthermore, the results show that VOTCA-XTP’s GW-BSE/MM is a powerful tool to study a wide range of types of electronic excitations in complex molecular environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transitioning GlideinWMS, a multi domain distributed workload manager, from GSI proxies to tokens and other granular credentials

GlideinWMS is a distributed workload manager that has been used in production for many years to provision resources for experiments like CERN’s CMS, many Neutrino experiments, and the OSG. Its security model was based mainly on GSI (Grid Security Infrastructure), using X.509 certificate proxies and VOMS (Virtual Organization Membership Service) extensions. Even when other credentials, like SSH keys, were used to authenticate with resources, proxies were also added all the time, to establish the identity of the requestor and the associated memberships or privileges. This single credential was used for everything and was, often implicitly, forwarded wherever needed. The addition of identity and access tokens and the phase-out of GSI forced us to reconsider the security model of GlideinWMS, to handle multiple credentials which can differ in type, technology, and functionality. Both identity tokens and access tokens are supported. GSI proxies even if no more mandatory, are still used, together with various JWT (JSON Web Token) based tokens and other certificates. The functionality of the credentials, defined by issuer, audience, and scope, also differ: a credential can allow access to a computing resource, or can protect the GlideinWMS framework from tampering, or can grant read or write access to storage, can provide an identity for accounting or auditing, or can provide a combination of any the formers. Furthermore, the tools in use do not include automatic forwarding and renewal of the new credentials so credential lifetime and renewal requirements became part of the discussion as well. In this paper, we will present how GlideinWMS was able to change its design and code to respond to all these changes.

97 MATHEMATICS AND COMPUTING↗

A new compact solid-state neutral particle analyser at ASDEX Upgrade: Setup and physics modeling

At ASDEX Upgrade (AUG), a new compact solid-state detector has been installed to measure the energy spectrum of fast neutrals based on the principle described by Shinohara et al. [Rev. Sci. Instrum. 75, 3640 (2004)]. The diagnostic relies on the usual charge exchange of supra-thermal fast-ions with neutrals in the plasma. Therefore, the measured energy spectra directly correspond to those of confined fast-ions with a pitch angle defined by the line of sight of the detector. Experiments in AUG showed the good signal to noise characteristics of the detector. It is energy calibrated and can measure energies of 40-200 keV with count rates of up to 140 kcps. The detector has an active view on one of the heating beams. The heating beam increases the neutral density locally; thereby, information about the central fast-ion velocity distribution is obtained. The measured fluxes are modeled with a newly developed module for the 3D Monte Carlo code F90FIDASIM [Geiger et al., Plasma Phys. Controlled Fusion 53, 65010 (2011)]. The modeling allows to distinguish between the active (beam) and passive contributions to the signal. Thereby, the birth profile of the measured fast neutrals can be reconstructed. This model reproduces the measured energy spectra with good accuracy when the passive contribution is taken into account.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Notification of Class 1 Permit Modification to Attachment B, Part A Application for the Los Alamos National Laboratory Hazardous Waste Facility Permit, EPA ID # NM0890010515

This document contains a notification for a Class 1 Permit Modification to modify the Los Alamos National Laboratory (LANL) Hazardous Waste Facility Permit (Permit) issued to the Department of Energy and the Los Alamos National Security, LLC, collectively known as the Permittees, in November 2010. All changes have been made to Attachment B, Part A Application, and a full, updated form, Revision 10.0, is enclosed. The modification updates the facility contact and address information, and the process codes and design capacity section for Technical Area (TA) 14; TA-16; TA-39; TA-54 West; TA-54, Area G; and TA-55 within the Attachment B, Part A Application. The process codes and design capacity section for TA-14, TA-16, and TA-39 was updated to incorporate a statement that the Permittees have notified of the intent to close units at these locations. The process design capacity for TA-54, Area G has been updated from 3,664,150 gallons to 4,346,590 gallons to account for the TA-54, Area G Pad 11 container storage unit that was added when the Permit was issued in 2010 but at the time did not previously store hazardous waste. The number of units at TA-54, Area G was also updated to 9 units to include this new unit. The process design capacity for TA-55 has also been updated from 178,500 gallons to 207,600 gallons to incorporate the addition of the TA-55-185 container storage unit that was also added when the Permit was issued. The number of storage units at TA-55 was updated to 7 to include the unit. The process design capacity at TA-54 West was changed to 11,660 gallons to correct a typographical error in previous versions of the form.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Application Results on Early Exascale Hardware

This Exascale Computing Project (ECP) milestone report summarizes the status of 27 of the 31 ECP Applications Development (AD) subprojects at the end of FY21. In November and December of 2021, a comprehensive assessment of AD projects was conducted by the ECP leadership along with external subject matter experts (SMEs). (NNSA application projects are reviewed separately using the ASC milestone process.) The AD review committee—consisting of the AD lead, AD deputy, Level 3 (L3), and at least one external project SME—was tasked with evaluating each project’s progress relative to ECP project goals specified in the FY21 timeline. Key areas of focus were code maturity and performance on pre-exascale systems, an in-depth analysis of final key performance parameter (KPP) verification contracts, and future R&D priorities in the final year of ECP and beyond. As such, this report contains not only an accurate snapshot of each subproject’s current status but also represents a broad account of successes and challenges in porting large scientific applications to DOE’s next-generation high-performance computing architectures – the Frontier and Aurora systems.

97 MATHEMATICS AND COMPUTING↗

ATOMIC Simulations and Experimental Data for CaCO3 Mixtures

This data consists of simulations and experimental measurements of laser-induced breakdown spectroscopy (LIBS). The simulations are produced by ATOMIC, a general purpose plasma modeling and kinetics code that has been designed to compute emission (or absorption) spectra from plasmas [1] and are used to develop a statistical characterization of matrix effects. Our overall suite of simulations includes contains several sets of simulations: training and validation sets of simulations for three and four element mixtures of calcium, carbon, oxygen, and nitrogen (included to account for atmosphere) along with simulations of the individual elements. The 4-element simulations include the mixture of all four elements mentioned and for each of the four individual elements. The 3-element simulations include output for the mixture of calcium, carbon, oxygen and for these three individual elements. The training data were produced using a 600-run design, shown in Figure 1, that varies input parameters temperature (T), electron density (Ne), and proportion of the elements calcium, carbon, oxygen, and nitrogen (Ca; C; O; N) for the 4 element output. The 3-element output includes all parameters except for the proportion of nitrogen. The element proportions (all the variables but T and Ne) sum to one and are unused in the single-element simulations. The validation data was produced with a 80-run design shown in Figure 2. The training and validation simulation outputs for the 4-element simulations for the mixture and for the single element calcium are shown as sample simulations in Figures 3 and 4 respectively. The simulations produce spectra over a range of 190nm - 950nm that roughly mimics the range collected by the SciAps Z-300 LIBS instrument that was used for the experimental data. The measured spectra for a CaCO3 (which may include contribution from Earth's atmosphere) in the experiment is shown in in Figure 5. All files are kept in directories whose names indicate the elemental composition (CaCO3, Ca, C, O, or N), number of elements (3 or 4), and purpose (training, which is not labeled in the file name, or validation) with file names numbered to indicate the line in the design files used to produce the simulation. The designs are provided as text files with names indicating their purpose. The experimental data is provided as a CSV file. [1] J Colgan, EJ Judge, DP Kilcrease, and JE Barefield II. Ab-initio modeling of an iron laser-induced plasma: Comparison between theoretical and experimental atomic emission spectra. Spectrochimica Acta Part B: Atomic Spectroscopy, 97:65{73}, 2014.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Viscous/inviscid interaction procedure for high-amplitude oscillating airfoils

A coupling procedure between an inviscid potential theoretical panel method and an unsteady interactive boundary layer calculation procedure has been developed for calculating unsteady airloads on oscillating airfoils in the high-incidence regime. A coordinate transformation into a time-dependent surface coordinate has been incorporated into the boundary layer code to avoid a quasi-steady treatment of the viscous flow equations in the leading edge region. The upper and lower sides of the airfoil are calculated with the complete set of unsteady equations. The interactive region of the boundary layer calculation can be straightforwardly extended over the complete upper surface of the airfoil while accounting for interaction areas in the region of strong adverse pressure gradients with incipient separation. Comparisons of sample calculations with experimental data show reasonably good agreement even for unsteady drag coefficients.

Geissler, W.↗

Reliability analysis of laminated CMC components through shell subelement techniques

An updated version of the integrated design program Composite Ceramics Analysis and Reliability Evaluation of Structures (C/CARES) was developed for the reliability evaluation of ceramic matrix composites (CMC) laminated shell components. The algorithm is now split into two modules: a finite-element data interface program and a reliability evaluation algorithm. More flexibility is achieved, allowing for easy implementation with various finite-element programs. The interface program creates a neutral data base which is then read by the reliability module. This neutral data base concept allows easy data transfer between different computer systems. The new interface program from the finite-element code Matrix Automated Reduction and Coupling (MARC) also includes the option of using hybrid laminates (a combination of plies of different materials or different layups) and allows for variations in temperature fields throughout the component. In the current version of C/CARES, a subelement technique was implemented, enabling stress gradients within an element to be taken into account. The noninteractive reliability function is now evaluated at each Gaussian integration point instead of using averaging techniques. As a result of the increased number of stress evaluation points, considerable improvements in the accuracy of reliability analyses were realized.

Starlinger, Alois↗

Analysis of a Bank Control Guidance for Aerocapture at Uranus

We apply the Fully Numerical Predictor-corrector Aerocapture Guidance to capture a vehicle into orbit for a mission at Uranus. Using the Genesis flight mechanics simulation, we analyze both undispersed and dispersed trajectories in order to tune the guidance parameters. We then assess the performance of the guidance using Monte Carlo analyses. Properly accounting for the oblateness of Uranus within the guidance proves to be critical. We identify high-to-low density gradients as the cause of large target orbit misses. Intentionally targeting a steeper flight-path angle at entry interface mitigates the risk of exiting the atmosphere on a hyperbolic orbit. Finally, We deliver the guidance code and tuned parameters for integration into the Program to Optimize Simulated Trajectories II.

Daniel A Matz↗

Simulation of divertor heat flux width on EAST by BOUT++ transport code

The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance of the AMOEBA Water Model in the Vicinity of QM Solutes: A Diagnosis Using Energy Decomposition Analysis

The importance of incorporating solvent polarization effects into the modeling of solvation processes has been well-recognized, and therefore a new generation of hybrid quantum mechanics/molecular mechanics (QM/MM) approaches that accounts for this effect is desirable. We present a fully self-consistent, mutually polarizable QM/MM scheme using the AMOEBA force field, in which the total energy of the system is variationally minimized with respect to both the QM electronic density and the MM induced dipoles. This QM/AMOEBA model is implemented through the Q-Chem/LibEFP code interface and then applied to the evaluation of solute-solvent interaction energies for various systems ranging from the water dimer to neutral and ionic solutes (NH 3 , NH 4 + , CN - ) surrounded by increasing numbers of water molecules (up to 100). In order to analyze the resulting interaction energies, we also utilize an energy decomposition analysis (EDA) scheme which identifies contributions from permanent electrostatics, polarization, and van der Waals (vdW) interaction for the interaction between the QM solute and the solvent molecules described by AMOEBA. This facilitates a component-wise comparison against full QM calculations where the corresponding energy components are obtained via a modified version of the absolutely localized molecular orbitals (ALMO)-EDA. The results show that the present QM/AMOEBA model can yield reasonable solute-solvent interaction energies for neutral and cationic species, while further scrutiny reveals that this accuracy highly relies on the delicate balance between insufficiently favorable permanent electrostatics and softened vdW interaction. For anionic solutes where the charge penetration effect becomes more pronounced, the QM/MM interface turns out to be unbalanced. These results are consistent with and further elucidate our findings in a previous study using a slightly different QM/AMOEBA model (Dziedzic et al. J. Chem. Phys. 2016, 145, 124106). The implications of these results for further refinement of this model are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Finite element analysis of two disk rotor system

A finite element model of simple horizontal rotor system is developed for evaluating its dynamic behaviour. The model is based on Timoshenko beam element and accounts for the effect of gyroscopic couple and other rotational forces. Present rotor system consists of single shaft which is supported by bearings at both ends and two disks are mounted at different locations. The natural frequencies, mode shapes and orbits of rotating system for a specific range of rotation speed are obtained by developing a MATLAB code for solving the finite element equations of rotary system. Consequently, Campbell diagram is plotted for finding a relationship between natural whirl frequencies and rotation of the rotor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SOLPS-ITER validation with TCV L-mode discharges

This work presents a quantitative test of SOLPS-ITER simulations against tokamak a configuration variable (TCV) L-mode experiments. These simulations account for drifts, currents, kinetic neutrals, and carbon impurities providing the most complete edge transport simulations for TCV to date. The comparison is performed on nominally identical discharges carried out to assess the effectiveness of TCV’s divertor baffles in the framework of the European Plasma Exhaust program and employs numerous edge diagnostics providing a detailed code-experiment benchmark for TCV. The simulations show a qualitative consistency, but the quantitative differences remain, which are assessed herein. It is found that, for a given separatrix density, the simulations most notably yield a colder, and denser, divertor state with a higher divertor neutral pressure than measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The alanine detector in BNCT dosimetry: Dose response in thermal and epithermal neutron fields

Purpose: The response of alanine solid state dosimeters to ionizing radiation strongly depends on particle type and energy. Due to nuclear interactions, neutron fields usually also consist of secondary particles such as photons and protons of diverse energies. Various experiments have been carried out in three different neutron beams to explore the alanine dose response behavior and to validate model predictions. Additionally, application in medical neutron fields for boron neutron capture therapy is discussed. Methods: Alanine detectors have been irradiated in the thermal neutron field of the research reactor TRIGA Mainz, Germany, in five experimental conditions, generating different secondary particle spectra. Further irradiations have been made in the epithermal neutron beams at the research reactors FiR 1 in Helsinki, Finland, and Tsing Hua open pool reactor in HsinChu, Taiwan ROC. Readout has been performed with electron spin resonance spectrometry with reference to an absorbed dose standard in a {sup 60}Co gamma ray beam. Absorbed doses and dose components have been calculated using the Monte Carlo codes FLUKA and MCNP. The relative effectiveness (RE), linking absorbed dose and detector response, has been calculated using the Hansen and Olsen alanine response model. Results: The measured dose response of the alanine detector in the different experiments has been evaluated and compared to model predictions. Therefore, a relative effectiveness has been calculated for each dose component, accounting for its dependence on particle type and energy. Agreement within 5% between model and measurement has been achieved for most irradiated detectors. Significant differences have been observed in response behavior between thermal and epithermal neutron fields, especially regarding dose composition and depth dose curves. The calculated dose components could be verified with the experimental results in the different primary and secondary particle fields. Conclusions: The alanine detector can be used without difficulty in neutron fields. The response has been understood with the model used which includes the relative effectiveness. Results and the corresponding discussion lead to the conclusion that application in neutron fields for medical purpose is limited by its sensitivity but that it is a useful tool as supplement to other detectors and verification of neutron source descriptions.

07 ISOTOPE AND RADIATION SOURCES↗