Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched field processing”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

THE CIRCUMGALACTIC MEDIUM OF MASSIVE GALAXIES AT z {approx} 3: A TEST FOR STELLAR FEEDBACK, GALACTIC OUTFLOWS, AND COLD STREAMS

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift {approx}3, using the Eris suite of cosmological hydrodynamic ''zoom-in'' simulations. The reference run adopts a blastwave scheme for supernova feedback that produces large-scale galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. The effect of the local UV radiation field is added in post-processing. The CGM (defined as all gas at R > 0.2 R{sub vir} = 10 kpc, where R{sub vir} is the virial radius) contains multiple phases having a wide range of physical conditions, with more than half of its heavy elements locked in a warm-hot component at T > 10{sup 5} K. Synthetic spectra, generated by drawing sightlines through the CGM, produce interstellar absorption-line strengths of Ly{alpha}, C II, C IV, Si II, and Si IV as a function of the galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high redshift by Steidel et al. The covering factor of absorbing material declines less rapidly with impact parameter for Ly{alpha} and C IV compared to C II, Si IV, and Si II, with Ly{alpha} remaining strong (W{sub Ly{alpha}} > 300 mA) to {approx}> 5 R{sub vir} = 250 kpc. Only about one third of all the gas within R{sub vir} is outflowing. The fraction of sightlines within one virial radius that intercept optically thick, N{sub H{sub I}}>10{sup 17.2} cm{sup -2} material is 27%, in agreement with recent observations by Rudie et al. Such optically thick absorption is shown to trace inflowing ''cold'' streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar by previous episodes of star formation in the main host and in nearby dwarfs, are the origin of strong (N{sub C{sub II}}>10{sup 13} cm{sup -2}) C II absorption with a covering factor of 22% within R{sub vir} and 10% within 2 R{sub vir}. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large O VI halo, with a typical column density N{sub O{sub VI}} {approx}> 10{sup 14} cm{sup -2} and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. Our zoom-in simulations of this single system appear then to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-T-434: Development of a Fan-Beam Optical Scanner Using CMOS Array for Small Field Dosimetry

Purpose: To design and construct a second generation optical computed tomography (OCT) system using a fan-beam with a CMOS array detector for the 3D dosimetry with polymer gel and radiochromic solid dosimeters. The system was specifically designed for the small field dosimetry. Methods: The optical scanner used a fan-beam laser, which was produced from a collimated red laser beam (λ=620 nm) with a 15-degree laser-line generating lens. The fan-beam was sent through an index-matching bath which holds the sample stage and a sample. The emerging laser light was detected with a 2.54 cm-long CMOS array detector (512 elements). The sample stage rotated through the full 360 degree projection angles at 0.9-degree increments. Each projection was normalized to the unirradiated sample at the projection angle to correct for imperfections in the dosimeter. A larger sample could be scanned by using a motorized mirror and linearly translating the CMOS detector. The height of the sample stage was varied for a full 3D scanning. The image acquisition and motor motion was controlled by a computer. The 3D image reconstruction was accomplished by a fan-beam reconstruction algorithm. All the software was developed inhouse with MATLAB. Results: The scanner was used on both PRESAGE and PAGAT gel dosimeters. Irreconcilable refraction errors were seen with PAGAT because the fan beam laser line refracted away from the detector when the field was highly varying in 3D. With PRESAGE, this type of error was not seen. Conclusion: We could acquire tomographic images of dose distributions by the new OCT system with both polymer gel and radiochromic solid dosimeters. Preliminary results showed that the system was more suited for radiochromic solid dosimeters since the radiochromic dosimeters exhibited minimal refraction and scattering errors. We are currently working on improving the image quality by thorough characterization of the OCT system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

The Principle of Maximum Conformality

A key problem in making precise perturbative QCD predictions is the uncertainty in determining the renormalization scale of the running coupling {alpha}{sub s}({mu}{sup 2}). It is common practice to guess a physical scale {mu} = Q which is of order of a typical momentum transfer Q in the process, and then vary the scale over a range Q/2 and 2Q. This procedure is clearly problematic since the resulting fixed-order pQCD prediction will depend on the renormalization scheme, and it can even predict negative QCD cross sections at next-to-leading-order. Other heuristic methods to set the renormalization scale, such as the 'principle of minimal sensitivity', give unphysical results for jet physics, sum physics into the running coupling not associated with renormalization, and violate the transitivity property of the renormalization group. Such scale-setting methods also give incorrect results when applied to Abelian QED. Note that the factorization scale in QCD is introduced to match nonperturbative and perturbative aspects of the parton distributions in hadrons; it is present even in conformal theory and thus is a completely separate issue from renormalization scale setting. The PMC provides a consistent method for determining the renormalization scale in pQCD. The PMC scale-fixed prediction is independent of the choice of renormalization scheme, a key requirement of renormalization group invariance. The results avoid renormalon resummation and agree with QED scale-setting in the Abelian limit. The PMC global scale can be derived efficiently at NLO from basic properties of the PQCD cross section. The elimination of the renormalization scheme ambiguity using the PMC will not only increases the precision of QCD tests, but it will also increase the sensitivity of colliders to new physics beyond the Standard Model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Dilepton radiation by vector mesons and off-shell partons in the hot and dense medium

According to the dynamical quasiparticle model (DQPM) - matched to reproduce lattice QCD results in thermodynamic limit, - the constituents of the strongly interacting quark-gluon plasma (sQGP) are massive and off-shell quasi-particles (quarks and gluons) with broad spectral functions. In order to address the electromagnetic radiation of the sQGP, we derive off-shell cross sections of qq-bar{yields}{gamma}*, qq-bar{yields}{gamma}*+g and qg{yields}{gamma}*q(q-barg{yields}{gamma}*q-bar) reactions taking into account the effective propagators for quarks and gluons from the DQPM. Dilepton production in In+In collisions at 158 AGeV is studied by implementing these processes into the parton-hadron-string dynamics (PHSD) transport approach. The microscopic PHSD transport approach describes the full evolution of the heavy-ion collision: from the dynamics of quasi-particles in the sQGP phase (when the local energy density is above {approx}1 GeV/fm{sup 3}) through hadronization and to the following hadron interactions and off-shell propagation after the hadronization. A comparison to the data of the NA60 Collaboration shows that the low mass dilepton spectra are well described by including a collisional broadening of vector mesons, while the spectra in the intermediate mass range are dominated by off-shell quark-antiquark annihilation, quark Bremsstrahlung and gluon-Compton scattering in the nonperturbative QGP. In particular, the observed softening of the m{sub T} spectra at intermediate masses (1 GeV {<=} M {<=} 3 GeV) is approximately reproduced.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Blind identification of full-field vibration modes of output-only structures from uniformly-sampled, possibly temporally-aliased (sub-Nyquist), video measurements

Enhancing the spatial and temporal resolution of vibration measurements and modal analysis could significantly benefit dynamic modelling, analysis, and health monitoring of structures. For example, spatially high-density mode shapes are critical for accurate vibration-based damage localization. In experimental or operational modal analysis, higher (frequency) modes, which may be outside the frequency range of the measurement, contain local structural features that can improve damage localization as well as the construction and updating of the modal-based dynamic model of the structure. In general, the resolution of vibration measurements can be increased by enhanced hardware. Traditional vibration measurement sensors such as accelerometers have high-frequency sampling capacity; however, they are discrete point-wise sensors only providing sparse, low spatial sensing resolution measurements, while dense deployment to achieve high spatial resolution is expensive and results in the mass-loading effect and modification of structure's surface. Non-contact measurement methods such as scanning laser vibrometers provide high spatial and temporal resolution sensing capacity; however, they make measurements sequentially that requires considerable acquisition time. As an alternative non-contact method, digital video cameras are relatively low-cost, agile, and provide high spatial resolution, simultaneous, measurements. Combined with vision based algorithms (e.g., image correlation or template matching, optical flow, etc.), video camera based measurements have been successfully used for experimental and operational vibration measurement and subsequent modal analysis. However, the sampling frequency of most affordable digital cameras is limited to 30–60 Hz, while high-speed cameras for higher frequency vibration measurements are extremely costly. This work develops a computational algorithm capable of performing vibration measurement at a uniform sampling frequency lower than what is required by the Shannon-Nyquist sampling theorem for output-only modal analysis. In particular, the spatio-temporal uncoupling property of the modal expansion of structural vibration responses enables a direct modal decoupling of the temporally-aliased vibration measurements by existing output-only modal analysis methods, yielding (full-field) mode shapes estimation directly. Then the signal aliasing properties in modal analysis is exploited to estimate the modal frequencies and damping ratios. Furthermore, the proposed method is validated by laboratory experiments where output-only modal identification is conducted on temporally-aliased acceleration responses and particularly the temporally-aliased video measurements of bench-scale structures, including a three-story building structure and a cantilever beam.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Heliostats with Adjustable Shape for High Concentration throughout the Day

Our motivation is to develop more efficient heliostats that can provide commercially viable solar thermal power at temperatures > 800°C. Such high temperatures will enable high-temperature industrial processes, as well as electrical generation after sunset with high efficiency. The importance of this research is that such heliostats have the potential to substantially expand the global use of solar energy, by adding solar thermal power as a major component. Thermal solar currently accounts for only 1% of all solar power (with PV being the rest), with heliostat fields providing just 0.25%. Our goals have been 1) to demonstrate a technical improvement for heliostats that can enable fields of them to more efficiently power receivers and reactors, and 2) to show a path to low-cost mass production. Our solution uses new opto-mechanical technology to correct a fundamental deficiency of present heliostats that limits their concentration, namely that they have fixed shape. Most of today’s heliostat research does not address this, but is directed simply toward cost reduction in an effort to make heliostats commercially viable. We are motivated to explore also improving heliostat efficiency, which can be done by continually changing their shape to maximize the concentration of sunlight throughout the day. This is not a new concept, but it has never been implemented in a practical, cost-effective way that approaches the theoretical limit to concentration while also improving mechanical performance; this is our goal. Our major accomplishments have been: 1) We have realized the planned design, construction, and test of a prototype heliostat that achieves the required shape changes in an 8 m2 single-piece glass mirror. The mirror is attached to a steel support frame that is automatically mechanically twisted by the heliostat drives that orient the mirror to direct sunlight to the tower-mounted receiver. Closed-loop tracking is done using a new beamsplitter camera that exploits the target-oriented mount configuration. Field tests of the heliostat show that the light is reflected through the day to always form a disc image of the sun, as needed to obtain the highest concentration. 2) We have developed the design for a field of 431 heliostats to deliver annual average of 1 MW of thermal power at 3,000 sun concentration, matched to a high-temperature ≥ 1000°C chemical reactor. 3) We have also developed, beyond the original stated goals of the project, a new concept for closed-loop tracking and shape-sensing for all the heliostats in the above field, using just 6 cameras around the concentrated reactor focus. Our research adds to the understanding of solar thermal energy by its demonstration of the technical effectiveness of a higher performing heliostat, and by its concept for a new powerful method for real-time tracking and shape sensing in the field, as described above. We have studied the economic feasibility of fields of our twisting heliostats to provide high- temperature heat at a price competitive with that of burning gas, to satisfy the DOE’s studied zero-emissions scenario, where the gas price has to include the cost of carbon capture. The project has the potential to greatly benefit the public if it helps limit global warming by 1) reducing carbon emission from industrial heating, which is currently a major contributor to the 40-billion-ton annual increase in atmospheric CO 2 . 2) Ultimately, the technology could prove to be the least expensive method to power direct air capture of CO 2 on the very large scale needed to remove the 1 trillion-ton excess of CO 2 already in the atmosphere.

14 SOLAR ENERGY↗
DOE OSTI2010

NEAR-INFRARED SPECTROSCOPY OF THE EXTRASOLAR PLANET HR 8799 b

We present 2.12-2.23 {mu}m high contrast integral field spectroscopy of the extrasolar planet HR 8799 b. Our observations were obtained with OSIRIS on the Keck II telescope and sample the 2.2 {mu}m CH{sub 4} feature, which is useful for spectral classification and as a temperature diagnostic for ultracool objects. The spectrum of HR 8799 b is relatively featureless, with little or no methane absorption, and does not exhibit the strong CH{sub 4} seen in T dwarfs of similar absolute magnitudes. The spectrum is consistent with field objects from early-L to T4 (3{sigma} confidence level), with a best-fitting type of T2. A similar analysis of the published 1-4 {mu}m photometry shows the infrared spectral energy distribution (SED) matches L5-L8 field dwarfs, especially the reddest known objects which are believed to be young and/or very dusty. Overall, we find that HR 8799 b has a spectral type consistent with L5-T2, although its SED is atypical compared to most field objects. We fit the 2.2 {mu}m spectrum and the infrared SED using the Hubeny and Burrows, Burrows et al., and Ames-Dusty model atmosphere grids, which incorporate non-equilibrium chemistry, non-solar metallicities, and clear and cloudy variants. No models agree with all of the data, but those with intermediate clouds produce significantly better fits. The largest discrepancy occurs in the J band, which is highly suppressed in HR 8799 b. Models with high eddy diffusion coefficients and high metallicities are somewhat preferred over those with equilibrium chemistry and solar metallicity. The best-fitting effective temperatures range from 1300 to 1700 K with radii between {approx}0.3 and 0.5 R{sub Jup}. These values are inconsistent with evolutionary model-derived values of 800-900 K and 1.1-1.3 R{sub Jup} based on the luminosity of HR 8799 b and the age of HR 8799, a discrepancy that probably results from imperfect atmospheric models or the limited range of physical parameters covered by the models. The low temperature inferred from evolutionary models indicates that HR 8799 b is {approx}400 K cooler than field L/T transition objects, providing further evidence that the L/T transition is gravity-dependent. With an unusually dusty photosphere, an exceptionally low luminosity for its spectral type, and hints of extreme secondary physical parameters, HR 8799 b appears to be unlike any class of field brown dwarf currently known.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SU-F-T-436: A Method to Evaluate Dosimetric Properties of SFGRT in Eclipse TPS

Purpose: The objective was to develop a method for dose distribution calculation of spatially-fractionated-GRID-radiotherapy (SFGRT) in Eclipse treatment-planning-system (TPS). Methods: Patient treatment-plans with SFGRT for bulky tumors were generated in Varian Eclipse version11. A virtual structure based on the GRID pattern was created and registered to a patient CT image dataset. The virtual GRID structure was positioned on the iso-center level together with matching beam geometries to simulate a commercially available GRID block made of brass. This method overcame the difficulty in treatment-planning and dose-calculation due to the lack o-the option to insert a GRID block add-on in Eclipse TPS. The patient treatment-planning displayed GRID effects on the target, critical structures, and dose distribution. The dose calculations were compared to the measurement results in phantom. Results: The GRID block structure was created to follow the beam divergence to the patient CT images. The inserted virtual GRID block made it possible to calculate the dose distributions and profiles at various depths in Eclipse. The virtual GRID block was added as an option to TPS. The 3D representation of the isodose distribution of the spatially-fractionated beam was generated in axial, coronal, and sagittal planes. Physics of GRID can be different from that for fields shaped by regular blocks because the charge-particle-equilibrium cannot be guaranteed for small field openings. Output factor (OF) measurement was required to calculate the MU to deliver the prescribed dose. The calculated OF based on the virtual GRID agreed well with the measured OF in phantom. Conclusion: The method to create the virtual GRID block has been proposed for the first time in Eclipse TPS. The dosedistributions, in-plane and cross-plane profiles in PTV can be displayed in 3D-space. The calculated OF’s based on the virtual GRID model compare well to the measured OF’s for SFGRT clinical use.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗
DOE OSTI2013

Combining prior day contours to improve automated prostate segmentation

Purpose: To improve the accuracy of automatically segmented prostate, rectum, and bladder contours required for online adaptive therapy. The contouring accuracy on the current image guidance [image guided radiation therapy (IGRT)] scan is improved by combining contours from earlier IGRT scans via the simultaneous truth and performance level estimation (STAPLE) algorithm. Methods: Six IGRT prostate patients treated with daily kilo-voltage (kV) cone-beam CT (CBCT) had their original plan CT and nine CBCTs contoured by the same physician. Three types of automated contours were produced for analysis. (1) Plan: By deformably registering the plan CT to each CBCT and then using the resulting deformation field to morph the plan contours to match the CBCT anatomy. (2) Previous: The contour set drawn by the physician on the previous day CBCT is similarly deformed to match the current CBCT anatomy. (3) STAPLE: The contours drawn by the physician, on each prior CBCT and the plan CT, are deformed to match the CBCT anatomy to produce multiple contour sets. These sets are combined using the STAPLE algorithm into one optimal set. Results: Compared to plan and previous, STAPLE improved the average Dice's coefficient (DC) with the original physician drawn CBCT contours to a DC as follows: Bladder: 0.81 {+-} 0.13, 0.91 {+-} 0.06, and 0.92 {+-} 0.06; Prostate: 0.75 {+-} 0.08, 0.82 {+-} 0.05, and 0.84 {+-} 0.05; and Rectum: 0.79 {+-} 0.06, 0.81 {+-} 0.06, and 0.85 {+-} 0.04, respectively. The STAPLE results are within intraobserver consistency, determined by the physician blindly recontouring a subset of CBCTs. Comparing plans recalculated using the physician and STAPLE contours showed an average disagreement less than 1% for prostate D98 and mean dose, and 5% and 3% for bladder and rectum mean dose, respectively. One scan takes an average of 19 s to contour. Using five scans plus STAPLE takes less than 110 s on a 288 core graphics processor unit. Conclusions: Combining the plan and all prior days via the STAPLE algorithm to produce treatment day contours is superior to the current standard of deforming only the plan contours to the daily CBCT. STAPLE also improves the precision, with a substantial decrease in standard deviation, a key for adaptive therapy. Geometrically and dosimetrically accurate contours can be automatically generated with STAPLE on prostate region kV CBCT in a time scale suitable for online adaptive therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Impact of cosmic variance on the galaxy-halo connection for Lyα emitters

In this paper we study the impact of cosmic variance and observational uncertainties in constraining the mass and occupation fraction, f{sub occ}, of dark matter (DM) halos hosting Lyα-emitting galaxies (LAEs) at high redshift. To this end, we construct mock catalogs from an N-body simulation to match the typical size of observed fields at z = 3.1 (∼1 deg{sup 2}). In our model a DM halo with mass in the range M{sub min}<M{sub h}<M{sub max} can only host one detectable LAE at most. We proceed to explore the parameter space determined by M{sub min}, M{sub max}, and f{sub occ} with a Markov Chain Monte Carlo algorithm using the angular correlation function and the LAEs’ number density as observational constraints. We find that the preferred minimum and maximum masses in our model span a wide range 10{sup 10.0}h{sup −1}M{sub ⊙}⩽M{sub min}⩽10{sup 11.1}h{sup −1}M{sub ⊙}, 10{sup 11.0}h{sup −1}M{sub ⊙}⩽M{sub max}⩽10{sup 13.0}h{sup −1}M{sub ⊙}, followed by a wide range in the occupation fraction 0.02⩽f{sub occ}⩽0.30. As a consequence, the median mass, M {sub 50}, of all the consistent models has a large uncertainty M{sub 50}=3.16{sub −2.37}{sup +9.34}×10{sup 10} h{sup −1}M{sub ⊙}. However, we find that the same individual models have a relatively tight 1σ scatter around the median mass ΔM{sub 1σ}=0.55{sub −0.31}{sup +0.11} dex. We are also able to show that f{sub occ} is uniquely determined by M{sub min}, regardless of M{sub max}. We argue that upcoming large surveys covering at least 25 deg{sup 2} should be able to put tighter constraints on M{sub min} and f{sub occ} through the LAE number density distribution width constructed over several fields of ∼1 deg{sup 2}.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2001

Strategies for Engaging NASA Earth Scientists in K-12 Education and Public Outreach

Engagement of the Earth Science research community in formal education at the kindergarten through high school level and in various aspects of informal education and in professional development of practitioners in related fields has been and continues to be a challenge. A range of approaches is being used and new ones are constantly being tried. Fundamental to our strategies is an understanding of the priorities, skills, academic experiences, motivation, rewards and work experiences of most scientists. It is within this context that efforts to engage a scientist in education efforts are attempted. A key strategy is to limit our requests to activities where the scientist's contribution of time and expertise can have the most impact. Don't waste the scientist's time! Time is one of their most prized resources, it is extremely valuable to you, and to them, we treat their time like a treasured resource. The clearer a scientist's role, their unique contribution and the finite nature of their effort, the more likely they are to participate. It is critical that commitments made to scientists are kept. If they want and can do more, great! Don't expect or assume more will be forthcoming. Another approach that we use is to create periodic venues that, among other things, serve to identify individuals who have an interest or inclination to con , tribute to education efforts. Once identified we strive to determine their interests so that we can make the best match between their interests and the needs of the education program or efforts. In this way, we try to make the best use of their time while engaging them in efforts which will be personally rewarding, and will further the overall education objectives. In addition, we try to make it easier for scientists to participate by providing focused training, such as development of their interviewing skills, and exposure to key concepts, knowledge and skills which are well known among educators but are not common knowledge among scientists. Another strategy with which we are still struggling is how to create and provide career meaningful rewards for individuals who demonstrate excellence in education equivalent to excellence in science. We do not yet have a yardstick to measure excellence in education nor is there a consensus among scientific peers that these two can be equivalent. None-the-less, methods to identify excellence in education, such as the peer review process, are being tried with some success. Use of solicitation and selection of educational efforts via a peer review process that is the same as for scientific research to identify unique, interesting and creative ideas has been somewhat effective. Furthermore, the application of the same peer review process to the output of an educational effort is used to validate and ensure the quality of this output. An example will be used to illustrate the application of some of these strategies to the development of a high school Earth and Space System Science curriculum created in partnership with a local county school system (Anne Arundel County, MD).

Meeson, Blanche W.↗
NASA NTRSDate not supplied

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗
DOE OSTI2019

Cross-scale modeling of storm surge, tide, and inundation in Mid-Atlantic Bight and New York City during Hurricane Sandy, 2012

Driven by high-resolution NAM (North America Mesoscale) and ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric fields, a 3-D unstructured grid SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model) was applied to simulate total water level in the Mid-Atlantic Bight during Hurricane Sandy in 2012. The simulated storm surge, tide, significant wave height, and wave peak period results compared favorably with NOAA observations along US East Coast, Long Island Sound, and New York Harbor. The maximum total water level at The Battery in the New York Harbor was accurately simulated with an absolute error of less than 0.08 m (out of 2.9 m peak surge) and a timing difference less than 10 min, which includes the effects on the order of 0.1–0.3 m (5–10%) from the coastal setup induced by the surface waves. The scenario comparisons of (1) “NAM” versus “ECMWF”, (2) “2-D” versus “3-D”, and (3) “with” versus “without” wind wave model were examined. The 3-D barotropic model forced by ECMWF including the effects of wind waves performs the best, attributed to wave-induced radiation stress and reduction of bottom stress when the 3-D version is used. Simultaneously, the ELCIRC-Sub model, using a 5-m topography/bathymetry sub-grid and a regular 200-m resolution finite volume computational model grid, was developed to simulate the street-level inundation in New York City. The momentum and mass fluxes calculated by the coarser base grid model were effectively coupled with the sub-grid so that running a full-blown high-resolution base model is not required. The ELCIRC-Sub model uses MPI (Message Passing Interface) parallel computing to enlarge the coverage of the land surface and an efficient non-linear solver to improve the accuracy of the wetting-and-drying processes. The temporal comparisons of modeled water level with NOAA's tidal stations and USGS' rapid-deployed gauges showed overall performance with an average error of 0.1 m. Particularly, ELCIRC-Sub captured the profile of the highest peak surge (3.9 m), which rose 2 m within 3 hour at Kings Point, NY. The spatial comparisons of the modeled peak water level at 80 surveyed locations showed an average error less than 0.13 m. The modeled maximum inundation extent also matched well with an 80% hit rate against FEMA's Hurricane Sandy maximum flooding extent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

SU-E-J-26: A Novel Technique for Markerless Self-Sorted 4D-CBCT Using Patient Motion Modeling: A Feasibility Study

Purpose: To develop an automatic markerless 4D-CBCT projection sorting technique by using a patient respiratory motion model extracted from the planning 4D-CT images. Methods: Each phase of onboard 4D-CBCT is considered as a deformation of one phase of the prior planning 4D-CT. The deformation field map (DFM) is represented as a linear combination of three major deformation patterns extracted from the planning 4D-CT using principle component analysis (PCA). The coefficients of the PCA deformation patterns are solved by matching the digitally reconstructed radiograph (DRR) of the deformed volume to the onboard projection acquired. The PCA coefficients are solved for each single projection, and are used for phase sorting. Projections at the peaks of the Z direction coefficient are sorted as phase 1 and other projections are assigned into 10 phase bins by dividing phases equally between peaks. The 4D digital extended-cardiac-torso (XCAT) phantom was used to evaluate the proposed technique. Three scenarios were simulated, with different tumor motion amplitude (3cm to 2cm), tumor spatial shift (8mm SI), and tumor body motion phase shift (2 phases) from prior to on-board images. Projections were simulated over 180 degree scan-angle for the 4D-XCAT. The percentage of accurately binned projections across entire dataset was calculated to represent the phase sorting accuracy. Results: With a changed tumor motion amplitude from 3cm to 2cm, markerless phase sorting accuracy was 100%. With a tumor phase shift of 2 phases w.r.t. body motion, the phase sorting accuracy was 100%. With a tumor spatial shift of 8mm in SI direction, phase sorting accuracy was 86.1%. Conclusion: The XCAT phantom simulation results demonstrated that it is feasible to use prior knowledge and motion modeling technique to achieve markerless 4D-CBCT phase sorting. National Institutes of Health Grant No. R01-CA184173 Varian Medical System.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2010

Efficiency Improvement of Nitride-Based Solid State Light Emitting Materials -- CRADA Final Report

The development of In{sub x}Ga{sub 1-x} N/GaN thin film growth by Molecular Beam Epitaxy has opened a new route towards energy efficient solid-state lighting. Blue and green LED's became available that can be used to match the whole color spectrum of visible light with the potential to match the eye response curve. Moreover, the efficiency of such devices largely exceeds that of incandescent light sources (tungsten filaments) and even competes favorably with lighting by fluorescent lamps. It is, however, also seen in Figure 1 that it is essential to improve on the luminous performance of green LED's in order to mimic the eye response curve. This lack of sufficiently efficient green LED's relates to particularities of the In{sub x}Ga{sub 1-x}N materials system. This ternary alloy system is polar and large strain is generated during a lattice mismatched thin film growth because of the significantly different lattice parameters between GaN and InN and common substrates such as sapphire. Moreover, it is challenging to incorporate indium into GaN at typical growth temperatures because a miscibility gap exists that can be modified by strain effects. As a result a large parameter space needs exploration to optimize the growth of In{sub x}Ga{sub 1-x}N and to date it is unclear what the detailed physical processes are that affect device efficiencies. In particular, an inhomogeneous distribution indium in GaN modifies the device performance in an unpredictable manner. As a result technology is pushed forward on a trial and error basis in particular in Asian countries such as Japan and Korea, which dominate the market and it is desirable to strengthen the competitiveness of the US industry. This CRADA was initiated to help Lumileds Lighting/USA boosting the performance of their green LED's. The tasks address the distribution of the indium atoms in the active area of their blue and green LED's and its relation to internal and external quantum efficiencies. Procedures to measure the indium distribution with near atomic resolution were developed and applied to test samples and devices that were provided by Lumilids. Further, the optical performance of the device materials was probed by photoluminescence, electroluminescence and time resolved optical measurements. Overall, the programs objective is to provide a physical basis for the development of a simulation program that helps making predictions to improve the growth processes such that the device efficiency can be increased to about 20%. Our study addresses all proposed aspects successfully. Carrier localization, lifetime and recombination as well as the strain-induced generation of electric fields were characterized and modeled. Band gap parameters and their relation to the indium distribution were characterized and modeled. Electron microscopy was developed as a unique tool to measure the formation of indium clusters on a nanometer length scale and it was demonstrated that strain induced atom column displacements can reliably be determined in any materials system with a precision that approaches 2 pm. The relation between the local indium composition x and the strain induced lattice constant c(x) in fully strained In{sub x}Ga{sub 1-x}N quantum wells was found to be: c(x) = 0.5185 + {alpha}x with {alpha} = 0.111 nm. It was concluded that the local indium concentration in the final product can be modulated by growth procedures in a predictable manner to favorably affect external quantum efficiencies that approached target values and that internal quantum efficiencies exceeded them.

36 MATERIALS SCIENCE↗
DOE OSTI2010

MID-INFRARED VARIABILITY FROM THE SPITZER DEEP WIDE-FIELD SURVEY

We use the multi-epoch, mid-infrared Spitzer Deep Wide-Field Survey to investigate the variability of objects in 8.1 deg{sup 2} of the NOAO Deep Wide Field Survey Booetes field. We perform a Difference Image Analysis of the four available epochs between 2004 and 2008, focusing on the deeper 3.6 and 4.5 {mu}m bands. Out of 474, 179 analyzed sources, 1.1% meet our standard variability selection criteria that the two light curves are strongly correlated (r>0.8) and that their joint variance ({sigma}{sub 12}) exceeds that for all sources with the same magnitude by 2{sigma}. We then examine the mid-IR colors of the variable sources and match them with X-ray sources from the XBooetes survey, radio catalogs, 24 {mu}m selected active galactic nucleus (AGN) candidates, and spectroscopically identified AGNs from the AGN and Galaxy Evolution Survey (AGES). Based on their mid-IR colors, most of the variable sources are AGNs (76%), with smaller contributions from stars (11%), galaxies (6%), and unclassified objects, although most of the stellar, galaxy, and unclassified sources are false positives. For our standard selection criteria, 11%-12% of the mid-IR counterparts to X-ray sources, 24 {mu}m AGN candidates, and spectroscopically identified AGNs show variability. The exact fractions depend on both the search depth and the selection criteria. For example, 12% of the 1131 known z>1 AGNs in the field and 14%-17% of the known AGNs with well-measured fluxes in all four Infrared Array Camera bands meet our standard selection criteria. The mid-IR AGN variability can be well described by a single power-law structure function with an index of {gamma} {approx} 0.5 at both 3.6 and 4.5 {mu}m, and an amplitude of S {sub 0} {approx_equal} 0.1 mag on rest-frame timescales of 2 yr. The variability amplitude is higher for shorter rest-frame wavelengths and lower luminosities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Validation of deformable image registration algorithms on CT images of ex vivo porcine bladders with fiducial markers

Purpose: The spatial accuracy of deformable image registration (DIR) is important in the implementation of image guided adaptive radiotherapy techniques for cancer in the pelvic region. Validation of algorithms is best performed on phantoms with fiducial markers undergoing controlled large deformations. Excised porcine bladders, exhibiting similar filling and voiding behavior as human bladders, provide such an environment. The aim of this study was to determine the spatial accuracy of different DIR algorithms on CT images ofex vivo porcine bladders with radiopaque fiducial markers applied to the outer surface, for a range of bladder volumes, using various accuracy metrics. Methods: Five excised porcine bladders with a grid of 30–40 radiopaque fiducial markers attached to the outer wall were suspended inside a water-filled phantom. The bladder was filled with a controlled amount of water with added contrast medium for a range of filling volumes (100–400 ml in steps of 50 ml) using a luer lock syringe, and CT scans were acquired at each filling volume. DIR was performed for each data set, with the 100 ml bladder as the reference image. Six intensity-based algorithms (optical flow or demons-based) implemented in theMATLAB platform DIRART, a b-spline algorithm implemented in the commercial software package VelocityAI, and a structure-based algorithm (Symmetric Thin Plate Spline Robust Point Matching) were validated, using adequate parameter settings according to values previously published. The resulting deformation vector field from each registration was applied to the contoured bladder structures and to the marker coordinates for spatial error calculation. The quality of the algorithms was assessed by comparing the different error metrics across the different algorithms, and by comparing the effect of deformation magnitude (bladder volume difference) per algorithm, using the Independent Samples Kruskal-Wallis test. Results: The authors found good structure accuracy without dependency on bladder volume difference for all but one algorithm, and with the best result for the structure-based algorithm. Spatial accuracy as assessed from marker errors was disappointing for all algorithms, especially for large volume differences, implying that the deformations described by the registration did not represent anatomically correct deformations. The structure-based algorithm performed the best in terms of marker error for the large volume difference (100–400 ml). In general, for the small volume difference (100–150 ml) the algorithms performed relatively similarly. The structure-based algorithm exhibited the best balance in performance between small and large volume differences, and among the intensity-based algorithms, the algorithm implemented in VelocityAI exhibited the best balance. Conclusions: Validation of multiple DIR algorithms on a novel physiological bladder phantom revealed that the structure accuracy was good for most algorithms, but that the spatial accuracy as assessed from markers was low for all algorithms, especially for large deformations. Hence, many of the available algorithms exhibit sufficient accuracy for contour propagation purposes, but possibly not for accurate dose accumulation.

60 APPLIED LIFE SCIENCES↗