SEARCH · Engineering Papers
Results for “distributed average tracking”
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.
Secondary-electron distribution for heavy ions.
Space averaged and extended track pictures of delta ray correction to transverse energy distribution associated with heavy ion secondary electrons
Radial Distribution of Electron Spectra from High-Energy Ions
The average track model describes the response of physical and biological systems using radial dose distribution as the key physical descriptor. We report on an extension of this model to describe the average distribution of electron spectra as a function of radial distance from an ion. We present calculations of these spectra for ions of identical linear energy transfer (LET), but dissimilar charge and velocity to evaluate the differences in electron spectra from these ions. To illustrate the usefulness of the radial electron spectra for describing effects that are not described by electron dose, we consider the evaluation of the indirect events in microdosimetric distributions for ions. We show that folding our average electron spectra model with experimentally determined frequency distributions for photons or electrons provides a good representation of radial event spectra from high-energy ions in 0.5-2 micrometer sites.
Predictions of Solar Cycle 24: How are We Doing?
Predictions of solar activity are an essential part of our Space Weather forecast capability. Users are requiring usable predictions of an upcoming solar cycle to be delivered several years before solar minimum. A set of predictions of the amplitude of Solar Cycle 24 accumulated in 2008 ranged from zero to unprecedented levels of solar activity. The predictions formed an almost normal distribution, centered on the average amplitude of all preceding solar cycles. The average of the current compilation of 105 predictions of the annual-average sunspot number is 106 +/- 31, slightly lower than earlier compilations but still with a wide distribution. Solar Cycle 24 is on track to have a below-average amplitude, peaking at an annual sunspot number of about 80. Our need for solar activity predictions and our desire for those predictions to be made ever earlier in the preceding solar cycle will be discussed. Solar Cycle 24 has been a below-average sunspot cycle. There were peaks in the daily and monthly averaged sunspot number in the Northern Hemisphere in 2011 and in the Southern Hemisphere in 2014. With the rapid increase in solar data and capability of numerical models of the solar convection zone we are developing the ability to forecast the level of the next sunspot cycle. But predictions based only on the statistics of the sunspot number are not adequate for predicting the next solar maximum. I will describe how we did in predicting the amplitude of Solar Cycle 24 and describe how solar polar field predictions could be made more accurate in the future.
Applications of amorphous track models in radiation biology
The average or amorphous track model uses the response of a system to gamma-rays and the radial distribution of dose about an ion's path to describe survival and other cellular endpoints from proton, heavy ion, and neutron irradiation. This model has been used for over 30 years to successfully fit many radiobiology data sets. We review several extensions of this approach that address objections to the original model, and consider applications of interest in radiobiology and space radiation risk assessment. In the light of present views of important cellular targets, the role of target size as manifested through the relative contributions from ion-kill (intra-track) and gamma-kill (inter-track) remains a critical question in understanding the success of the amorphous track model. Several variations of the amorphous model are discussed, including ones that consider the radial distribution of event-sizes rather than average electron dose, damage clusters rather than multiple targets, and a role for repair or damage processing.
Simulating self-powered neutron detector responses to infer burnup-induced power distribution perturbations in next-generation light water reactors
Understanding how 3D power distribution will be monitored throughout reactor core volumetric space in next-generation nuclear power reactors is crucial to the design, deployment, and licensing of these reactors. Although numerous techniques exist for 3D power distribution monitoring based on the response of both in situ and ex situ sensors currently implemented or proposed for use in the US reactor fleet, crucial details about these techniques are often unclear. The publicly available documentation does not include information such as how well these techniques are characterized and optimized in their implementations and the levels of uncertainty in the inferred 3D power distribution. The work described herein investigated a recently developed 3D power distribution inferencing method as applied to two next-generation reactor simulations: (1) the NuScale small modular reactor design and (2) the Westinghouse AP1000 design, both of which contain in-core strings of vanadium self-powered neutron detectors (SPNDs). This investigation considered a range of SPND string sensor densities, as well as a range of 3D power distribution axial segment sizes. In this work, SPND response simulation is informed by neutron flux calculations in representative homogenized cores. For the different sensor densities and power distribution axial segment sizes in these simulations, the average solution error, solver iterations, and run time were tracked to parameterize the sensor-core configuration.
Fission track astrology of three Apollo 14 gas-rich breccias
The three Apollo 14 breccias 14301, 14313, and 14318 all show fission xenon due to the decay of Pu-244. To investigate possible in situ production of the fission gas, an analysis was made of the U-distribution in these three breccias. The major amount of the U lies in glass clasts and in matrix material and no more than 25% occurs in distinct high-U minerals. The U-distribution of each breccia is discussed in detail. Whitlockite grains in breccias 14301 and 14318 found with the U-mapping were etched and analyzed for fission tracks. The excess track densities are much smaller than indicated by the Xe-excess. Because of a preirradiation history documented by very high track densities in feldspar grains, however, it is impossible to attribute the excess tracks to the decay of Pu-244. A modified track method has been developed for measuring average U-concentrations in samples containing a heterogeneous distribution of U in the form of small high-U minerals. The method is briefly discussed, and results for the rocks 14301, 14313, 14318, 68815, 15595, and the soil 64421 are given.
Characterization of Operational Performance of Ka-Band Links in Deep Space Network
This paper presents a characterization of performance of 32-GHz Ka-band link in an operational environment. The data come from tracking of Kepler spacecraft by the NASA Deep Space Network (DSN) in the past few years. Ka-band link offers a significant signal-to-noise advantage compared to the more common X- or S-band; however, it is subject to more signal fluctuation caused by the weather. Kepler is the first mission supported by the DSN that solely rely on the higher deep space Ka-band communications link to return its high-rate science data. A well characterization of the operational performance of Kepler would benefit future Ka-band missions, especially for those operating with smaller link margin. The study examines how weather conditions at the DSN facilities (e.g., winds, clouds, rains) affect the received signal, particularly on telemetry data. It addresses questions such as how often the weather affects the link and how much degradation the link could suffer. Among the 22 Ka-band passes in 2012, heavy clouds affected one pass and two passes were impacted by high winds. The adverse weather caused at time as much as 3-dB instantaneous change in the signal to noise ratio (SNR). Estimates of degradation to the SNR as a function of wind speeds are captured based on observations. This study also quantifies the probability distribution of the variation of received signal power. With such information, future missions can better plan their link design. An optimal link design aims for just having enough margin to realize the data return with a targeted probability, with neither having too much margin that it reduces the downlink data rate nor insufficient reserve that causes frequent data outages. The study tries to provide some answers to the questions such as how much does SNR vary from one tracking pass to the next, and what is the cumulative distribution of signal fluctuation. Regarding the signal variation over many tracking passes, we found that the averaged symbol SNR (SSNR) for each pass changed by as much as 4 dB. Some variations were due to geometry such as the changing distance between spacecraft and Earth and different antenna pointing elevations. Other variations seem to be random in nature, reflecting the randomness of operational environments. Some passes were found to be quite stable, with a standard deviation of symbol SNR around 0.25 dB while others had greater variation, up to 1.4 dB. Overall, the cumulative distribution of the symbol SNR reflects that 50% of the fluctuations was less than 0.6 dB, 90% of fluctuation was within 1.6 dB and 95% within 2.2 dB. Some of the challenges in data processing, validation and modeling were captured in this paper. We observed some inconsistent measurements where operational data significantly deviate from a normal expectation. These inconsistent measurements made it hard to develop an accurate and consistent performance model, such as the degradation of signal SNR as a function of high wind. The difficulty was compounded by the fact that there were only a few Ka-band passes observed in high winds.
Track Structure Model for Radial Distributions of Electron Spectra and Event Spectra from High-Energy Ions
An analytic method is described for evaluating the average radial electron spectrum and the radial and total frequency-event spectrum for high-energy ions. For high-energy ions, indirect events make important contributions to frequency-event spectra. The method used for evaluating indirect events is to fold the radial electron spectrum with measured frequency-event spectrum for photons or electrons. The contribution from direct events is treated using a spatially restricted linear energy transfer (LET). We find that high-energy heavy ions have a significantly reduced frequency-averaged final energy (yF) compared to LET, while relativistic protons have a significantly increased yF and dose-averaged lineal energy (yD) for typical site sizes used in tissue equivalent proportional counters. Such differences represent important factors in evaluating event spectra with laboratory beams, in space- flight, or in atmospheric radiation studies and in validation of radiation transport codes. The inadequacy of LET as descriptor because of deviations in values of physical quantities, such as track width, secondary electron spectrum, and yD for ions of identical LET is also discussed.
Lidar Altimeter Measurements of Canopy Structure: Methods and Validation for Closed Canopy, Broadleaf Forests
Lidar altimeter observations of vegetated landscapes provide a time-resolved measure of laser pulse backscatter energy from canopy surfaces and the underlying ground. Airborne lidar altimeter data was acquired using the Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) for a successional sequence of four, closed-canopy, deciduous forest stands in eastern Maryland. The four stands were selected so as to include a range of canopy structures of importance to forest ecosystem function, including variation in the height and roughness of the outer-most canopy surface and the vertical organization of canopy stories and gaps. The character of the SLICER backscatter signal is described and a method is developed that accounts for occlusion of the laser energy by canopy surfaces, transforming the backscatter signal to a canopy height profile (CHP) that quantitatively represents the relative vertical distribution of canopy surface area. The transformation applies an increased weighting to the backscatter amplitude as a function of closure through the canopy and assumes a horizontally random distribution of the canopy components. SLICER CHPs, averaged over areas of overlap where lidar ground tracks intersect, are shown to be highly reproducible. CHP transects across the four stands reveal spatial variations in vegetation, at the scale of the individual 10 m diameter laser footprints, within and between stands. Averaged SLICER CHPs are compared to analogous height profile results derived from ground-based sightings to plant intercepts measured on plots within the four stands. Tbe plots were located on the segments of the lidar ground tracks from which averaged SLICER CHPs were derived, and the ground observations were acquired within two weeks of the SLICER data acquisition to minimize temporal change. The differences in canopy structure between the four stands is similarly described by the SLICER and ground-based CHP results, however a Chi-square test of similarity documents differences that are statistically significant. The differences are discussed in terms of measurement properties that define the smoothness of the resulting CHPs and Lidar Altimeter Measurements of Canopy Structure - Harding et al. canopy properties that may vertically bias the CHP representations of canopy structure. The statistical differences are most likely due to the more noisy character of the ground-based CHPs, especially high in the canopy where ground-based sightings are rare resulting in an underestimate of canopy surface area and height, and to departures from the assumption of horizontal randomness which bias the CHPs toward the observer (upward for SLICER and downward for ground-based CHPs). The results demonstrate that the SLICER observations reliably provide a measure of canopy structure that reveals ecologically interesting structural variations such as those characterizing a successional sequence of closed-canopy, broadleaf forest stands.
Second Cancers After Fractionated Radiotherapy: Stochastic Population Dynamics Effects
When ionizing radiation is used in cancer therapy it can induce second cancers in nearby organs. Mainly due to longer patient survival times, these second cancers have become of increasing concern. Estimating the risk of solid second cancers involves modeling: because of long latency times, available data is usually for older, obsolescent treatment regimens. Moreover, modeling second cancers gives unique insights into human carcinogenesis, since the therapy involves administering well characterized doses of a well studied carcinogen, followed by long-term monitoring. In addition to putative radiation initiation that produces pre-malignant cells, inactivation (i.e. cell killing), and subsequent cell repopulation by proliferation can be important at the doses relevant to second cancer situations. A recent initiation/inactivation/proliferation (IIP) model characterized quantitatively the observed occurrence of second breast and lung cancers, using a deterministic cell population dynamics approach. To analyze ifradiation-initiated pre-malignant clones become extinct before full repopulation can occur, we here give a stochastic version of this I I model. Combining Monte Carlo simulations with standard solutions for time-inhomogeneous birth-death equations, we show that repeated cycles of inactivation and repopulation, as occur during fractionated radiation therapy, can lead to distributions of pre-malignant cells per patient with variance >> mean, even when pre-malignant clones are Poisson-distributed. Thus fewer patients would be affected, but with a higher probability, than a deterministic model, tracking average pre-malignant cell numbers, would predict. Our results are applied to data on breast cancers after radiotherapy for Hodgkin disease. The stochastic IIP analysis, unlike the deterministic one, indicates: a) initiated, pre-malignant cells can have a growth advantage during repopulation, not just during the longer tumor latency period that follows; b) weekend treatment gaps during radiotherapy, apart from decreasing the probability of eradicating the primary cancer, substantially increase the risk of later second cancers.
Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas
Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.
Distributed Optimal Power Management for Battery Energy Storage Systems: A Novel Accelerated Tracking ADMM Approach
Optimal power management (OPM) is critical for large-scale battery energy storage systems. Today’s methods often require formidable computational effort due to the design based on centralized numerical optimization. Thus, this paper investigates computationally distributed OPM where the agents based on the cells communicate over a network to cooperatively solve the OPM problem. We propose an accelerated tracking alternating direction method of multipliers (ADMM) algorithm to solve the distributed OPM. The proposed algorithm embeds dynamic average consensus and Nesterov’s acceleration technique in the ADMM algorithm. Not only is the proposed algorithm fully distributed without a need for fusion or aggregating nodes, but it also accelerates convergence. The paper formulates the OPM in a model predictive control framework where it seeks to regulate the charging/discharging power of each battery cell to minimize the total power losses and promote balanced use of the constituent cells while complying with the safety constraints. The paper provides ample simulation results to demonstrate the effectiveness and advantages of the proposed distributed OPM in terms of computation and convergence.
An airfoil-based synthetic actuator disk model for wind turbine aerodynamic and structural analysis
Here, this study introduces an airfoil-based refinement technique to enhance the Actuator Disk Model (ADM) for improved wind turbine aerodynamic load prediction and structural simulation in conjunction with Large Eddy Simulations of the wind flow. While ADM offers higher computational efficiency than the more detailed but resource-intensive Actuator Line Model (ALM), it traditionally lacks the resolution needed to capture the localized blade forces accurately. To address this limitation, we introduce a refinement technique that uses airfoil-specific data and employs interpolation-based grid point refinement, achieving ALM-comparable accuracy while preserving ADM's efficiency. Unlike conventional ADM that provides only rotor-disk averaged forces, our synthetic method tracks transient aerodynamic load variations over multiple blade revolutions, allowing us to calculate the distributions of maximum and minimum loads during typical cycles. Applied to the NREL 5 MW reference turbine, our enhanced ADM accurately predicts key aerodynamic parameters (angle of attack, axial velocity, lift, drag, axial and tangential forces along the blades) as well as structural responses (blade tip deflection, maximum stress, and stress concentration). Our results show that the tip deflection ranges from 2.33m (3.69 % of blade length) to 4.28m (6.79 %), with maximum stress concentration occurring near the blade root. This research demonstrates that a refined synthetic ADM approach can serve as a computationally efficient alternative for both aerodynamic analysis and structural simulation of wind turbine blades subjected to realistic wind fields.
A Neural Network Parametrization of Volumetric Cloud Fraction Profiles Using Satellite Observations and MERRA-2 Reanalysis Meteorological Data
Clouds play a crucial role in regulating the hydrologic cycle and Earth's radiative energy budget, yet they are often poorly represented in global climate models (GCMs). This study applies deep machine learning techniques to develop a physical parameterization of volumetric cloud fraction (VCF), the fraction of a 3-D grid volume occupied by clouds using satellite lidar-radar measurements. The neural network (NN) captures the complicated relationships between observed VCF profiles and collocated meteorological variables from MERRA-2 reanalysis data. Our results show that the NN model, particularly a sequence-to-sequence long short-term memory (LSTM) network with a sixfactor loss function, effectively learns the underlying cloud physical processes. The NN model outperforms MERRA-2 reanalysis in representing low-level clouds in tropical and subtropical regions and low- and middle-level clouds over midlatitude storm-track regions, and also improves VCF histograms. These improvements are reflected in the vertical distributions of zonally, meridionally, and globally averaged VCFs, geographic distributions of low-, middle-, and high-level clouds, and seasonal variations in monthly-mean VCF. Furthermore, the NN predictions effectively capture the El Niño-Southern Oscillation (ENSO) effects and other interannual variations. The NN parameterization is further evaluated through a sensitivity analysis, in which a single predictor is perturbed at a time. This reveals that relative humidity (RH) is the dominant factor influencing variations in globally averaged VCF at low and middle altitudes, followed by temperature. At higher altitudes, temperature becomes the primary driver of VCF through its effect on RH. Changes in wind components had minimal impact on globally averaged VCF.
Estimation of spatially restricted LET using track structure models
The spatial distribution of energy deposition is an important determinant in the formation of biologically significant lesions. It has been widely realized that Linear Energy Transfer (LET) being an average quantity is not sufficient to describe the situation at a submicroscopic scale. To remedy this to some extent 'energy-cut-off' values are sometimes used but since they are related to secondary electron energy and only indirectly to their range they are also not adequate although they may be easily calculated. 'Range-restricted LET' appears to be better but its determination is usually quite involved. Xapsos (1992) suggested a semi-empirical approximation based on a modified Bethe-formula which contains a number of assumption which are difficult to verify. A simpler and easier way is to use existing beam-models which describe energy deposition around an ion's path. They all agree that the energy density (i. e., energy deposited per unit mass) decreases with the inverse square of the distance from the track center. This simple dependence can be used to determine the fraction of total LET which is deposited in a cylinder of a given radius. As an example our own beam model. Energy density depends on distance x (measured in m) from the track center according to the presented formula.
Advances in Measuring Antarctic Sea-Ice Thickness and Ice-Sheet Elevations with ICESat Laser Altimetry
NASA's Ice, Cloud and Land Elevation Satellite (ICESat) has been measuring elevations of the Antarctic ice sheet and sea-ice freeboard elevations with unprecedented accuracy. Since February 20,2003, data has been acquired during three periods of laser operation varying from 36 to 54 days, which is less than the continuous operation of 3 to 5 years planned for the mission. The primary purpose of ICESat is to measure time-series of ice-sheet elevation changes for determination of the present-day mass balance of the ice sheets, study of associations between observed ice changes and polar climate, and estimation of the present and future contributions of the ice sheets to global sea level rise. ICESat data will continue to be acquired for approximately 33 days periods at 3 to 6 month intervals with the second of ICESat's three lasers, and eventually with the third laser. The laser footprints are about 70 m on the surface and are spaced at 172 m along-track. The on-board GPS receiver enables radial orbit determinations to an accuracy better than 5 cm. The orbital altitude is around 600 km at an inclination of 94 degrees with a 8-day repeat pattern for the calibration and validation period, followed by a 91 -day repeat period for the rest of the mission. The expected range precision of single footprint measurements was 10 cm, but the actual range precision of the data has been shown to be much better at 2 to 3 cm. The star-tracking attitude-determination system should enable footprints to be located to 6 m horizontally when attitude calibrations are completed. With the present attitude calibration, the elevation accuracy over the ice sheets ranges from about 30 cm over the low-slope areas to about 80 cm over areas with slopes of 1 to 2 degrees, which is much better than radar altimetry. After the first period of data collection, the spacecraft attitude was controlled to point the laser beam to within 50 m of reference surface tracks over the ice sheets. Detection of ice elevation changes over select areas of the ice sheet is demonstrated with using both crossover analysis and precise-repeat track analysis. Sea ice freeboard-height distributions over the Antarctic sea pack are derived over distances of 50 km and converted into maps of average freeboard thickness and sea-ice thickness.
Self-Consistent Magnetosphere-Ionosphere Coupling and Associated Plasma Energization Processes
Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.