Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear model”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Relativistic weak lensing from a fully non-linear cosmological density field

In this paper we examine cosmological weak lensing on non-linear scales and show that there are Newtonian and relativistic contributions and that the latter can also be extracted from standard Newtonian simulations. We use the post-Friedmann formalism, a post-Newtonian type framework for cosmology, to derive the full weak-lensing deflection angle valid on non-linear scales for any metric theory of gravity. We show that the only contributing term that is quadratic in the first order deflection is the expected Born correction and lens-lens coupling term. We use this deflection angle to analyse the vector and tensor contributions to the E- and B- mode cosmic shear power spectra. In our approach, once the gravitational theory has been specified, the metric components are related to the matter content in a well-defined manner. Specifying General Relativity, we write down a complete set of equations for a GR+ΛCDM universe for computing all of the possible lensing terms from Newtonian N-body simulations. We illustrate this with the vector potential and show that, in a GR+ΛCDM universe, its contribution to the E-mode is negligible with respect to that of the conventional Newtonian scalar potential, even on non-linear scales. Thus, under the standard assumption that Newtonian N-body simulations give a good approximation of the matter dynamics, we show that the standard ray tracing approach gives a good description for a ΛCDM cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Drag reduction and the Vogel exponent of a flexible beam in transient shear flows

Interactions between a flexible beam and a fluid in a channel are of great relevance to biological hairy surfaces, aquatic vegetation, marine life (e.g., fish gills), and many industrial systems alike. While steady state response of a beam to such flows is fairly well-explored, their behavior in the transient regime is not fully understood. Here, a series of numerical simulations are performed to study the laminar Couette flow of an incompressible viscous fluid past an elastic beam in a two-dimensional channel. The flexible beam is perpendicular to the direction of flow, and its base is fixed to the stationary bottom of the channel. We measure the evolution of the Vogel exponent, drag reduction, and reconfiguration number during the transient and steady-state response of the fluid–structure system for different geometrical and physical properties. Our benchmark shows a good agreement between numerical and experimental observations. Our results show that the system's steady-state response at different bulk-fluid velocities can be reproduced by investigating the shear flow response during the transient regime. We define a new variable that characterizes the evolution of the local velocity profile in the proximity of the free end of the beam and use that to characterize the transient-regime response. The analysis yields insight into the competing effects of elasticity of the beam and non-linear flow response.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Scaling effects in the static and dynamic response of graphite-epoxy beam-columns

Scale model technology represents one method of investigating the behavior of advanced, weight-efficient composite structures under a variety of loading conditions. It is necessary, however, to understand the limitations involved in testing scale model structures before the technique can be fully utilized. These limitations, or scaling effects, are characterized. in the large deflection response and failure of composite beams. Scale model beams were loaded with an eccentric axial compressive load designed to produce large bending deflections and global failure. A dimensional analysis was performed on the composite beam-column loading configuration to determine a model law governing the system response. An experimental program was developed to validate the model law under both static and dynamic loading conditions. Laminate stacking sequences including unidirectional, angle ply, cross ply, and quasi-isotropic were tested to examine a diversity of composite response and failure modes. The model beams were loaded under scaled test conditions until catastrophic failure. A large deflection beam solution was developed to compare with the static experimental results and to analyze beam failure. Also, the finite element code DYCAST (DYnamic Crash Analysis of STructure) was used to model both the static and impulsive beam response. Static test results indicate that the unidirectional and cross ply beam responses scale as predicted by the model law, even under severe deformations. In general, failure modes were consistent between scale models within a laminate family; however, a significant scale effect was observed in strength. The scale effect in strength which was evident in the static tests was also observed in the dynamic tests. Scaling of load and strain time histories between the scale model beams and the prototypes was excellent for the unidirectional beams, but inconsistent results were obtained for the angle ply, cross ply, and quasi-isotropic beams. Results show that valuable information can be obtained from testing on scale model composite structures, especially in the linear elastic response region. However, due to scaling effects in the strength behavior of composite laminates, caution must be used in extrapolating data taken from a scale model test when that test involves failure of the structure.

Jackson, Karen E.↗
DOE OSTI2014

Application of the TITAN Hybrid Deterministic Transport Code to Medical Physics - Paper 118

The TITAN code is a 3-dimensional parallel deterministic transport code that solves the linear Boltzmann equation. The TITAN code has hybrid algorithms that are ideal for medical physics applications: i) SN and characteristics method and ii) SN with a fictitious quadrature set and ray-tracing. These hybrid algorithms have been used in the simulation of computed tomography (CT) and single photon emission computed tomography (SPECT) In this paper, the hybrid algorithms developed for medical physics applications in the TITAN code are discussed. First, the TITAN code's block-oriented approach for utilizing SN and characteristic solvers for CT is presented and results for a simple test case are given. Second, the paper describes a hybrid formulation for SPECT with an SN solver in the phantom, a fictitious quadrature set for directions from the phantom surface passing through the collimator, and a simplified ray-tracing algorithm for particle transport to the detector. Finally, ongoing work is described and concluding remarks are given. The TITAN code's hybrid algorithms allow it to efficiently solve problems with low-scatter regions, which are frequently found in medical physics applications. The TITAN code has been used to solve a simplified CT problem and generate projection images for SPECT simulations. The WCOS technique was developed to improve the TITAN code's collimator representation. For these applications, the TITAN code was shown to require significantly less computation time than MCNP5. Work is currently underway to implement an iterative image reconstruction algorithm in the TITAN code. A script that calls the TITAN code to complete the forward projection step in SPECT has already been used to reconstruction a 2-dimensional phantom. The ability of TITAN to fully model particle transport in the patient in short computation times is expected to be a great advantage to iterative reconstruction. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Detector and Physics Performance at a Muon Collider

A muon collider represents the ideal machine to reach very high center-of-mass energies and luminosities by colliding elementary particles. This is the result of the low level of beamstrahlung and synchrotron radiation compared to linear or circular electron-positron colliders. In contrast with other lepton machines, the design of a detector for a multi-TeV muon collider requires detailed knowledge of the interaction region due to the significant backgrounds created by muon beam decays in the collider ring. The physics reach can be properly evaluated only when the detector performance in such an environment is determined. In this work, the backgrounds generated by muon beams of 750 GeV are characterized and the performance of the tracking system and the calorimeter detector is illustrated. Solutions to minimize the effect of the beam-induced backgrounds are discussed and applied to obtain track and jet reconstruction performance. The μ + μ – → $Hν\bar{ν} → b\bar{b} ν\bar{ν}$ process is fully simulated and reconstructed to demonstrate that physics measurements are possible in this harsh environment. The measurement precision for the Higgs boson coupling to $b\bar{b}$ is evaluated for $\sqrt{s}$=1.5, 3, and 10 TeV and compared to other proposed machines.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1998

Evaluation of risk from space radiation with high-energy heavy ion beams

The most challenging radiation in space consists of fully ionized atomic elements with high energy for which only the few lowest energy ions can be stopped in shielding materials. The health risk from exposure to these ions and their secondary radiations generated in shield materials is poorly understood since there are few human data and a systematic study in relevant animal model systems has not been made. The accuracy of risk prediction is described as the major limiting factor in the management of space radiation risk. The expected impact of systematic studies is examined using the limited available biological data and models. Given the limitations of current predictions, models must be developed that are able to incorporate the required fundamental scientific data into accurate risk estimates. The important radiation components that can be provided for laboratory testing are identified. The use of ground-based accelerator beams to simulate space radiation is explained and quantitative scientific constraints on such facilities are derived. Three facilities, one each in the United States, in Germany and in Japan, currently have the partial capability to satisfy these constraints. A facility has been proposed using the Brookhaven National Laboratory Booster Synchrotron in the United States; in conjuction with other on-site accelerators, it will be able to provide the full range of heavy ion beams and energies required.

NASA Discipline Radiation Health↗
DOE OSTI2015

SU-E-T-652: Quantification of Dosimetric Uncertainty of I-125 COMS Eye Plaque

Purpose: To compare dosimetrically three plan calculation systems (Plato, Varian Brachytherapy, and in-house-made Excel) available for I-125 COMS eye plaque treatment with measurement. Methods: All systems assume homogeneous media and calculations are based on a three-dimensional Cartesian coordinates, Plato and Brachytherapy Planning are based on AAPM TG-43 and the in-house Excel program only on inverse square corrections. Doses at specific depths were measured with EBT3 Gafchromic film from a fully loaded and a partially loaded 16 mm plaque (13 and 8 seeds respectively, I-125, model 6711 GE, Oncura). Measurements took place in a water tank, utilizing solid water blocks and a 3D-printed plaque holder. Taking advantage that gafchromic film has low energy dependence, a dose step wedge was delivered with 6 MV photon beam from a Varian 2100 EX linac for calibration. The gray-scale to dose in cGy was obtained with an Epson Expression 10000 XL scanner in the green channel. Treatment plans were generated for doses of 2200 cGy to a depth of 7 mm, and measurements were taken on a sagittal plane. Results: The calculated dose at the prescription point was 2242, 2344, and 2211 cGy with Excel, Brachyvision and Plato respectively for a fully loaded plaque, for the partially loaded plaque the doses were 2266, 2477, and 2193 cGy respectively. At 5 mm depth the doses for Brachyvision and Plato were comparable (3399 and 3267 cGy respectively), however, the measured dose in film was 3180 cGy which was lower by as much as 6.4% in the fully loaded plaque and 7.6% in the partially loaded plaque. Conclusion: Careful methodology and calibration are essential when measuring doses at specific depth due to the sensitivity and rapid dose fall off of I-125.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2000

Novel Integration Radial and Axial Magnetic Bearing

Typically, fully active magnetically suspended systems require one axial and two radial magnetic bearings. Combining radial and axial functions into a single device allows for more compact and elegant packaging. Furthermore, in the case of high-speed devices such as energy storage flywheels, it is beneficial to minimize shaft length to keep rotor mode frequencies as high as possible. Attempts have been made to combine radial and axial functionality, but with certain drawbacks. One approach requires magnetic control flux to flow through a bias magnet reducing control effectiveness, thus resulting in increased resistive losses. This approach also requires axial force producing magnetic flux to flow in a direction into the rotor laminate that is undesirable for minimizing eddy-current losses resulting in rotational losses. Another approach applies a conical rotor shape to what otherwise would be a radial heteropolar magnetic bearing configuration. However, positional non-linear effects are introduced with this scheme and the same windings are used for bias, radial, and axial control adding complexity to the controller and electronics. For this approach, the amount of axial capability must be limited. It would be desirable for an integrated radial and axial magnetic bearing to have the following characteristics; separate inputs for radial and axial control for electronics and control simplicity, all magnetic control fluxes should only flow through their respective air gaps and should not flow through any bias magnets for minimal resistive losses, be of a homopolar design to minimize rotational losses, position related non-linear effects should be minimized, and dependent upon the design parameters, be able to achieve any radial/axial force or power ratio as desired. The integrated radial and axial magnetic bearing described in this paper exhibits all these characteristics. Magnetic circuit design, design equations, and magnetic field modeling results will be presented.

Blumenstock, Kenneth↗
DOE OSTI2025

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Rapid adaptation to Coriolis force perturbations of arm trajectory

1. Forward reaching movements made during body rotation generate tangential Coriolis forces that are proportional to the cross product of the angular velocity of rotation and the linear velocity of the arm. Coriolis forces are inertial forces that do not involve mechanical contact. Virtually no constant centrifugal forces will be present in the background when motion of the arm generates transient Coriolis forces if the radius of body rotation is small. 2. We measured the trajectories of arm movements made in darkness to a visual target that was extinguished as movement began. The reaching movements were made prerotation, during rotation at 10 rpm in a fully enclosed rotating room, and postrotation. During testing the subject was seated at the center of the room and pointed radially. Neither visual nor tactile feedback about movement accuracy was present. 3. In experiment 1, subjects reached at a fast or slow rate and their hands made contact with a horizontal surface at the end of the reach. Their initial perrotary movements were highly significantly deviated relative to prerotation in both trajectories and end-points in the direction of the transient Coriolis forces that had been generated during the reaches. Despite the absence of visual and tactile feedback about reaching accuracy, all subjects rapidly regained straight movement trajectories and accurate endpoints. Postrotation, transient errors of opposite sign were present for both trajectories and endpoints. 4. In a second experiment the conditions were identical except that subjects pointed just above the location of the extinguished target so that no surface contact was involved. All subjects showed significant initial perrotation deviations of trajectories and endpoints in the direction of the transient Coriolis forces. With repeated reaches the trajectories, as viewed from above, again became straight, but there was only partial restoration of endpoint accuracy, so that subjects reached in a straight line to the wrong place. Aftereffects of opposite sign were transiently present in the postrotary movements. 5. These observations fail to support current equilibrium point models, both alpha and lambda, of movement control. Such theories would not predict endpoint errors under our experimental conditions, in which the Coriolis force is absent at the beginning and end of a movement. Our results indicate that detailed aspects of movement trajectory are being continuously monitored on the basis of proprioceptive feedback in relation to motor commands. Adaptive compensations can be initiated after one perturbation despite the absence of either visual or tactile feedback about movement trajectory and endpoint error. Moreover, movement trajectory and end-point can be remapped independently.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗
NASA NTRS2013

NASA Tech Briefs, April 2013

Topics covered include: Fully Integrated, Miniature, High-Frequency Flow Probe Utilizing MEMS Leadless SOI Technology; Nanoscale Surface Plasmonics Sensor With Nanofluidic Control; Advanced Dispersed Fringe Sensing Algorithm for Coarse Phasing Segmented Mirror Telescopes; Neural Network Back-Propagation Algorithm for Sensing Hypergols; Bulk Moisture and Salinity Sensor; Change-Based Satellite Monitoring Using Broad Coverage and Targetable Sensing; Circularly Polarized Microwave Antenna Element with Very Low Off-Axis Cross-Polarization; Ultra-Low Heat-Leak, High-Temperature Superconducting Current Leads for Space Applications; Flash Cracking Reactor for Waste Plastic Processing; An Automated Safe-to-Mate (ASTM) Tester; Wireless Chalcogenide Nanoionic-Based Radio-Frequency Switch; Compute Element and Interface Box for the Hazard Detection System; DOT Transmit Module; Composite Aerogel Multifoil Protective Shielding; Li-Ion Electrolytes with Improved Safety and Tolerance to High-Voltage Systems; Polymer-Reinforced, Non-Brittle, Lightweight Cryogenic Insulation; Controlled, Site-Specific Functionalization of Carbon Nanotubes with Diazonium Salts; Regenerable Sorbent for CO2 Removal; Sprayable Aerogel Bead Compositions With High Shear Flow Resistance and High Thermal Insulation Value; Lexan Linear Shaped Charge Holder with Magnets and Backing Plate; Robotic Ankle for Omnidirectional Rock Anchors; Wind, Wave, and Tidal Energy Without Power Conditioning; An Active Heater Control Concept to Meet IXO Type Mirror Module Thermal-Structural Distortion Requirement; Waterless Clothes-Cleaning Machine; Integrated Electrical Wire Insulation Repair System; LVGEMS Time-of-Flight Mass Spectrometry on Satellites; Surface Inspection Tool for Optical Detection of Surface Defects; Per-Pixel, Dual-Counter Scheme for Optical Communications; Certification-Based Process Analysis; Surface Navigation Using Optimized Waypoints and Particle Swarm Optimization; Smart-Divert Powered Descent Guidance to Avoid the Backshell Landing Dispersion Ellipse; Estimating Foreign-Object-Debris Density from Photogrammetry Data; Adaptive Sampling of Spatiotemporal Phenomena with Optimization Criteria; Building a 2.5D Digital Elevation Model From 2D Imagery; Eyes on the Earth 3D; Target Trailing With Safe Navigation for Maritime Autonomous Surface Vehicles; Adams-Based Rover Terramechanics and Mobility Simulator - ARTEMIS; ISTP CDF Skeleton Editor; Uplink Summary Generator (ULSGEN) Version 1.0; Robotics On-Board Trainer (ROBoT); Software Engineering Tools for Scientific Models; Automatic Data Filter Customization Using a Genetic Algorithm; Tracker Toolkit; Towards Efficient Scientific Data Management Using Cloud Storage; On a Formal Tool for Reasoning About Flight Software Cost Analysis; A Nanostructured Composites Thermal Switch Controls Internal and External Short Circuit in Lithium Ion Batteries; Spacecraft Crew Cabin Condensation Control; and Functional Near-Infrared Spectroscopy Signals Measure Neuronal Activity in the Cortex.

Source record↗
DOE OSTI2025

Across the Scales of the Nucleus: Understanding Short Range Correlations from Medium Modification to Probe Independence

The atomic nucleus presents an intricate system due to the non-linear forces described by Quantum Chromodynamics (QCD) that govern its structure. The range of scales involved is remarkable; the most massive nuclei weigh approximately five orders of magnitude more than the quarks that compose them. The nucleus can be analyzed at various levels, from quarks to hadrons to the nucleus as a whole. Short-Range Correlations (SRCs) within the nucleus play a significant role that spans these diverse scales. At the most fundamental level, SRCs influence the interaction between nucleons. The nucleon-nucleon (NN) interaction, arising from QCD, is crucial in determining nuclear properties. SRCs serve as valuable probes for measuring this NN interaction, as the nucleons within SRCs become effectively decoupled from the rest of the nucleus. Multiple experimental techniques, including electron scattering, have been employed to investigate the NN interaction through SRCs. However, our first project demonstrates that inclusive measurements alone are inadequate to constrain this interaction fully. Moving to the scale of the nucleus, SRCs contribute to the high-momentum tail of the nuclear spectral function. While the low-momentum region is characterized by nucleons exhibiting bulk properties, nucleons begin to pair into SRCs at higher momenta. Our research aims to bridge the understanding between the mean-field portion of the nucleus and its high-momentum SRC components. Additionally, SRCs affect the quark structure of protons, as evidenced by the EMC effect, which indicates that quarks behave differently when protons are embedded within a nucleus—an effect referred to as medium modification. This thesis explores the correlation between SRCs and medium modification across various experimental setups. Finally, we seek to establish an interpretation of the nuclear ground state. Accomplishing this requires demonstrating that our SRC observables are independent of the probe’s scale and scheme. The concluding project of this thesis illustrates how we utilize triple coincidence quasi-elastic scattering across a range of (Q2) values to develop a model-dependent framework for understanding SRC distributions within the nucleus’s ground-state wavefunction.

Denniston, A. W. Denniston [University of Tel Aviv↗
DOE OSTI2017

Reconstructing Solar Particle Event Spectra from Absorbed Dose Measurements

The space radiation environment is unlike anything encountered on Earth. It is dominated by heavy charged particles consisting of highly charged nuclei with atomic numbers ranging from 1 to 94 with a broad spectrum of kinetic energies. The three primary sources of radiation in space are trapped belts, galactic cosmic rays (GCRs), and solar particle events (SPEs). The trapped belts are only of concern near bodies with significant magnetic fields and GCR flux spectra are relatively well-predicted using models that account for solar modulation of the local interstellar spectrum. SPEs, which consist primarily of protons, are difficult to predict in terms of occurrence, duration, intensity, and energy spectrum. SPE intensity and energy spectrum are critical for understanding radiation exposure-related health effects for astronauts in space. The National Aeronautics and Space Administration (NASA) relies on data from the Geostationary Operational Environmental Satellite series to characterize intensities and energy spectra of SPEs near Earth. However, these data are not applicable for spacecraft in other regions of the solar system. To understand the radiation environment for astronauts embarking on travel beyond Earth orbit, NASA is developing the Hybrid Electronic Radiation Assessor (HERA), a distributed dosimetry system that employs a solid-state silicon sensor coupled with the Timepix readout chip. The pixelated nature of the detector permits detailed analyses of charged particle tracks which can then be used to characterize the fluence spectrum at the location of the detector. Despite the unprecedented level of detail regarding the radiation environment provided by these detectors, knowledge of the fluence spectrum at all locations within the vehicle, specifically inside of astronauts, is desired to fully characterize the associated risk of harmful health effects. Radiation transport of the external radiation environment, considering the shielding provided by the vehicle, other objects, and even astronauts, is performed to obtain this information. Thus, there is clearly a need to define the external radiation environment using measurements performed within the vehicle. The number of detectors that may be flown as a part of HERA is limited by logistics, including data, power, and mass restrictions. The current configuration for crewed flight consists of six detectors at various locations with in the Orion Multi-Purpose Crew Vehicle (MPCV). It follows that an external proton fluence spectrum can be characterized with at most six energy bins with traditional matrix inversion techniques. The purpose of the present study is to determine whether it is possible to reconstruct SPE spectra with greater resolution from absorbed dose measurements and knowledge of the response of individual detectors using linear programming.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2015

SU-D-213-01: Transparent Photon Detector For The Online Monitoring Of IMRT Beams

Purpose: An innovative Transparent Detector for Radiotherapy (TraDeRa) has been developed. The detector aims at real-time monitoring of modulated beam ahead of the patient during delivery sessions, with a field cover up to 40×40 cm {sup 2}. Methods: TraDeRa consists in a pixelated matrix of ionization chambers with a patented electrodes design. An in-house designed specific integrated circuit allows to extract the signal and provides a real-time map of beam intensity and shape, at the linac pulse-scale. The measurements under irradiation are made with a 6 MV clinical X-Ray beam. Dose calculations are performed with the Monte Carlo code PENELOPE, modeling the full accelerator head and the TraDeRa detector. Results: A 2 % attenuation of the beam was measured in the presence of TraDeRa and the PENELOPE dosimetric study showed no significant modification of the photon beam properties. TraDeRa detects error leaf position as small as 1 mm compared to a reference field, for both static and modulated fields. In addition, measurements are accurate over a large dynamic range from low intensity signals, as inter-leaves leaks, to very high intensities as obtained on the medical line of the European Synchrotron Radiation Facility. The detector is fully operational for conventional and high dose rate beams as FFF modes (up to 2400 MU/min). Conclusion: The current version of TraDeRa shows promising results for IMRT quality assurance (QA), allowing pulse-scale monitoring of the beam and high sensitivity for errors detection. The attenuation is small enough not to hinder the irradiation while keeping the beam upstream of the patient under constant control. A final prototype under development will include 1600 independent electrodes, half of them with a high resolution centered on the beam axis. This compact detector provides an independent set of measurements for a better QA. Funding support : This work was supported by the LABEX PRIMES (ANR-11-LABX-0063) of Universite de Lyon, within the program “Investissements d’Avenir” (ANR-11-IDEX-0007) operated by the French National Research Agency (ANR) and within the project “INSPIRA” operated by the OSEO institution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Cosmic bubble and domain wall instabilities I: parametric amplification of linear fluctuations

This is the first paper in a series where we study collisions of nucleated bubbles taking into account the effects of small initial (quantum) fluctuations in a fully 3+1-dimensional setting. In this paper, we consider the evolution of linear fluctuations around highly symmetric though inhomogeneous backgrounds. In particular, we demonstrate that a large degree of asymmetry develops over time from tiny initial fluctuations superposed upon planar and SO(2,1) symmetric backgrounds. These fluctuations are inevitable consequences of zero-point vacuum oscillations, so excluding them by enforcing a high degree of spatial symmetry is inconsistent in a quantum treatment. To simplify the analysis we consider the limit of two colliding planar walls, with mode functions for the fluctuations characterized by the wavenumber transverse to the collision direction and a longitudinal shape along the collision direction x, which we solve for. In the linear regime, the fluctuations obey a linear wave equation with a time- and space-dependent mass m{sub eff}(x,t). In situations where the walls collide multiple times, m{sub eff} oscillates in time. We use Floquet theory to study the evolution of the fluctuations and generalize the calculations familiar from the preheating literature to the case with many coupled degrees of freedom. The inhomogeneous case has bands of unstable transverse wavenumbers k{sub ⊥} whose corresponding mode functions grow exponentially. By examining the detailed spatial structure of the mode functions in x, we identify both broad and narrow parametric resonance generalizations of the homogeneous m{sub eff}(t) case of preheating. The unstable k{sub ⊥} modes are longitudinally localized, yet can be described as quasiparticles in the Bogoliubov sense. We define an effective occupation number and show they are created in bursts for the case of well-defined collisions in the background. The transverse-longitudinal coupling accompanying nonlinearity radically breaks this localized particle description, with nonseparable 3D modes arising that will be studied in subsequent papers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗