Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Quantifying Rates of Erosion Using Historical Aerial Photos - 19478

Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. In this work multi-temporal historical aerial images are analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. Scans of photographic negatives were obtained for multiple photographs and several negatives were scanned with a calibrated photogrammetric scanner by Geomni (formerly Blue Skies Consulting) in Belen, New Mexico. Other negatives were provided by USDA or scanned on a high-resolution, large format media scanner in the United States National Archives Cartographic Research Room, College Park, Maryland (1939). An additional set of aerial photographs was obtained from the USGS EROS Data Center web site. For other years, photographic prints were available for aerial photography. Several prints were scanned with a large-format media scanner at 1200 dpi. Ground control points for selected digital images used in the stereoscopic analyses were identified in ArcMap. Ground control points were collected by matching features which could be identified in the aerial photography as well as the orthophotography or shaded relief imagery from the Lidar dataset. Such features included fence line intersections, road intersections, railroad bridges, and structures, depending on the age of the photography. Many old fence lines and abandoned roads remained visible in the shaded relief imagery while they were completely obscured by tree canopy in modern photography. Ground control elevations were sampled from the Lidar-based DEM or directly from point clouds. The image coordinates and matching three-dimensional projected coordinates of the ground control points were exported to IDL where parameters for a three-dimensional direct linear transformation (DLT) were obtained. Built-in singular value decomposition routines were used to solve an overdetermined set of equations for the eleven parameters of the DLT for a particular historical photograph. A solution to the DLT parameters was obtained for each historical aerial photograph in which features of interest were identified using the stereoscope. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. There was inevitable uncertainty associated with the estimation of points from the historical aerial imagery. This includes estimation of changes in locations of features like gully heads through time. In this work a method to estimate the uncertainty associated with measurements of locations of features using the stereoscopic analysis is developed and applied. For features with locations in space that do not change (e.g., road intersections, etc.) the difference in measurements of locations from the Lidar dataset and that of the estimate from the historical aerial image provides an observation of the measurement error. These points used to estimate the measurement error are referred to as control points. For each of the historical aerial images, a different measurement error distribution is estimated using control points. These measurement errors are combined with the estimated distances from the stereoscopic analyses to produce distributions of gully retreat rates and hillslope erosion rates used along with output from Erosion Working Group (EWG) Landscape Evolution Models (LEMs) to provide statistical probability distributions for input into the PPA Model. Specifically, gully head retreat rate and gully widening rate were characterized using statistical probability distributions. The central limit theorem was then applied to develop rates that represent long-term averages that are used to inform a PPA Model built in GoldSim. The importance of these rates with respect to characterization of future risk can then assessed using sensitivity analyses. This provides an approach to determine if subsequent reduction of the uncertainty associated with the probability distributions of erosion rates will correspond to meaningful reduction in the estimation of future risk from the site. Ultimately this is a general method that utilizes historical aerial imagery to estimate rates of erosion with uncertainty. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Clinical introduction of image lag correction for a cone beam CT system

Purpose: Image lag in the flat-panel detector used for Linac integrated cone beam computed tomography (CBCT) has a degrading effect on CBCT image quality. The most prominent visible artifact is the presence of bright semicircular structure in the transverse view of the scans, known also as radar artifact. Several correction strategies have been proposed, but until now the clinical introduction of such corrections remains unreported. In November 2013, the authors have clinically implemented a previously proposed image lag correction on all of their machines at their main site in Amsterdam. The purpose of this study was to retrospectively evaluate the effect of the correction on the quality of CBCT images and evaluate the required calibration frequency. Methods: Image lag was measured in five clinical CBCT systems (Elekta Synergy 4.6) using an in-house developed beam interrupting device that stops the x-ray beam midway through the data acquisition of an unattenuated beam for calibration. A triple exponential falling edge response was fitted to the measured data and used to correct image lag from projection images with an infinite response. This filter, including an extrapolation for saturated pixels, was incorporated in the authors’ in-house developed clinical CBCT reconstruction software. To investigate the short-term stability of the lag and associated parameters, a series of five image lag measurement over a period of three months was performed. For quantitative analysis, the authors have retrospectively selected ten patients treated in the pelvic region. The apparent contrast was quantified in polar coordinates for scans reconstructed using the parameters obtained from different dates with and without saturation handling. Results: Visually, the radar artifact was minimal in scans reconstructed using image lag correction especially when saturation handling was used. In patient imaging, there was a significant reduction of the apparent contrast from 43 ± 16.7 to 15.5 ± 11.9 HU without the saturation handling and to 9.6 ± 12.1 HU with the saturation handling, depending on the date of the calibration. The image lag correction parameters were stable over a period of 3 months. The computational load was increased by approximately 10%, not endangering the fast in-line reconstruction. Conclusions: The lag correction was successfully implemented clinically and removed most image lag artifacts thus improving the image quality. Image lag correction parameters were stable for 3 months indicating low frequency of calibration requirements.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1991

Robust control of systems with real parameter uncertainty and unmodelled dynamics

During this research period we have made significant progress in the four proposed areas: (1) design of robust controllers via H infinity optimization; (2) design of robust controllers via mixed H2/H infinity optimization; (3) M-delta structure and robust stability analysis for structured uncertainties; and (4) a study on controllability and observability of perturbed plant. It is well known now that the two-Riccati-equation solution to the H infinity control problem can be used to characterize all possible stabilizing optimal or suboptimal H infinity controllers if the optimal H infinity norm or gamma, an upper bound of a suboptimal H infinity norm, is given. In this research, we discovered some useful properties of these H infinity Riccati solutions. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H infinity norm. We also set up a detailed procedure for applying the H infinity theory to robust control systems design. The desire to design controllers with H infinity robustness but H(exp 2) performance has recently resulted in mixed H(exp 2) and H infinity control problem formulation. The mixed H(exp 2)/H infinity problem have drawn the attention of many investigators. However, solution is only available for special cases of this problem. We formulated a relatively realistic control problem with H(exp 2) performance index and H infinity robustness constraint into a more general mixed H(exp 2)/H infinity problem. No optimal solution yet is available for this more general mixed H(exp 2)/H infinity problem. Although the optimal solution for this mixed H(exp 2)/H infinity control has not yet been found, we proposed a design approach which can be used through proper choice of the available design parameters to influence both robustness and performance. For a large class of linear time-invariant systems with real parametric perturbations, the coefficient vector of the characteristic polynomial is a multilinear function of the real parameter vector. Based on this multilinear mapping relationship together with the recent developments for polytopic polynomials and parameter domain partition technique, we proposed an iterative algorithm for coupling the real structured singular value.

Chang, Bor-Chin↗
NASA NTRS1993

Atmospheric moisture fields derived by satellite observations over the tropical Pacific Ocean

Values of precipitable water are retrieved over the tropical and subtropical Pacific Ocean from TOVS infrared and microwave channel brightness temperature and OLR observations by means of stepwise linear regression. The most useful temperature and moisture sensing channels are pre-selected from sensitivity tests of a radiative transfer model. Numerous models are developed and tested against collocated radiosonde observations and Nimbus-7 SMMR precipitable water estimates. For RAOB comparisons, the best estimator used 15 TOVS predictors and captured 71.1 deg percent of the variance (+0.62 g/sq cm standard error) for column precipitable water; for precipitable water of 700-500 mb bulk layer, these values were 71.7 percent and +/- 0.17 g/sq cm. Little skill of estimated precipitable water was obtained for moisture above 500 mb. Regressions were less skillful against SMMR, unless collocation parameters were tightly controlled; SMMR was less acceptable than RAOB's because of observational drift and errors. Generally, the most skillful predictors were boundary layer brightness temperatures of TOVS channels and satellite estimated stability indices. 'Moisture channels' were hardly useful except for estimating middle and upper tropospheric moisture. Additional regression models were constructed testing the sensitivity to different observational and meteorological characteristics. Models which used some in situ observations surface observations or stabilities calculated from RAOB, were the most successful. The best of these explained 87.5 percent of the variance but the regression selected almost no TOVS channels, relying instead on conventional RAOB and surface observations. A set of four regression models were developed, stratifying atmospheric characteristics on the basis of collocated OLR values. These models improved the variance explained by 5.0 percent; the model associated with the highest OLR values (275 W/sq m less than or equal to OLR; that is, no cloud) showed only marginal skill. Precipitable water fields were generated from the best TOVS-only model for seven days in January 1983 and compared with SMMR-estimated fields. OLR fields and ECMWF precipitable water analysis. The TOVS regression model compared favorably to the SMMR analysis in amplitudes and features. It revealed more evolving synoptic signal than the ECMWF analysis. In synoptically active regions, it differed with respect to the OLR analysis, primarily because of actual differences in the vertical distribution of water vapor.

Chung, Hyosang↗
DOE OSTI2022

High-Temperature Linear Receiver Enabled by Multicomponent Aerogels

Concentrating solar thermal (CST) technology has significant potential mainly due to its dispatchability and low cost of storage. However, to compete with other sources, including utility-scale solar PV, its final cost (cents/kWh) still needs to be lowered. Cost reduction can be achieved by improving the system level efficiency of the CST plants through the deployment of advanced power cycles, which operate at high temperatures of ~700°C. However, optical, and thermal losses pose a major challenge to the efficiency of such CST systems. The overall aim of this project is to investigate and de-risk a linear solar receiver concept called an Aerogel Insulated Receiver (AIR) that generates high temperatures (up to 700°C) at a low solar concentration ratio (<100) and with a high collection efficiency (optical × receiver). Our prior work has demonstrated the thermal stability1 and optical and heat-insulating properties2,3 of transparent aerogel insulation at a one-inch scale. The focus of this work is on (1) co-optimization of the geometry of the aerogel tile and receiver enclosure to fit a standard parabolic collector (PTC), (2) scale-up of aerogels into 4-inch tiles while preserving key properties, (3) experimental measurement of receiver heat loss (W/m) in a >70-cm test stand and validation of anticipated receiver performance at high temperatures. Regarding (1), appropriate optical and thermal models for a parabolic trough receiver (PTR) are developed and validated. The geometry of the aerogels and the receiver enclosure are co-optimized to maximize the collection efficiency. The model predicts a 54% collection efficiency at 700°C for an AIR design based on flat aerogels. By combining the collection efficiency with the power block efficiency of supercritical CO 2 cycles, we predict >10% improvements in peak plant efficiency relative to existing line-focusing CST systems. The application of curved plasmonic aerogels is predicted to further increase the collection efficiency to 64%. Regarding (2), we demonstrate the successful development of 6-inch-long refractory aerogel tiles with optical, thermal, and stability characteristics consistent with our prior work. This scale-up requires a transition to a larger ALD station and modifying the ALD process variables such as exposure time and the number of precursor doses. Regarding (3), we design and develop an AIR test stand measuring 3 feet in length. Heat loss performance analysis is carried out using the test stand. The results indicate that aerogel insulation can significantly reduce receiver thermal losses at the high operating temperatures required for next-generation PTRs. The experimental results agree with the heat loss performance predicted by our receiver model. Lastly, we conducted preliminary failure mode and effects (FMEA) and techno-economic (TEA) analyses to identify failure mitigation strategies and commercial opportunities, respectively. Overall, this project identifies key opportunities and challenges in deploying aerogel insulating receivers in next-generation line-focusing CST technologies.

14 SOLAR ENERGY↗
DOE OSTI2015

A non-modal analytical method to predict turbulent properties applied to the Hasegawa-Wakatani model

Linear eigenmode analysis often fails to describe turbulence in model systems that have non-normal linear operators and thus nonorthogonal eigenmodes, which can cause fluctuations to transiently grow faster than expected from eigenmode analysis. When combined with energetically conservative nonlinear mode mixing, transient growth can lead to sustained turbulence even in the absence of eigenmode instability. Since linear operators ultimately provide the turbulent fluctuations with energy, it is useful to define a growth rate that takes into account non-modal effects, allowing for prediction of energy injection, transport levels, and possibly even turbulent onset in the subcritical regime. We define such a non-modal growth rate using a relatively simple model of the statistical effect that the nonlinearities have on cross-phases and amplitude ratios of the system state variables. In particular, we model the nonlinearities as delta-function-like, periodic forces that randomize the state variables once every eddy turnover time. Furthermore, we estimate the eddy turnover time to be the inverse of the least stable eigenmode frequency or growth rate, which allows for prediction without nonlinear numerical simulation. We test this procedure on the 2D and 3D Hasegawa-Wakatani model [A. Hasegawa and M. Wakatani, Phys. Rev. Lett. 50, 682 (1983)] and find that the non-modal growth rate is a good predictor of energy injection rates, especially in the strongly non-normal, fully developed turbulence regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Interaction of the body, head, and eyes during walking and turning

Body, head, and eye movements were measured in five subjects during straight walking and while turning corners. The purpose was to determine how well the head and eyes followed the linear trajectory of the body in space and whether head orientation followed changes in the gravito-inertial acceleration vector (GIA). Head and body movements were measured with a video-based motion analysis system and horizontal, vertical, and torsional eye movements with video-oculography. During straight walking, there was lateral body motion at the stride frequency, which was at half the frequency of stepping. The GIA oscillated about the direction of heading, according to the acceleration and deceleration associated with heel strike and toe flexion, and the body yawed in concert with stepping. Despite the linear and rotatory motions of the head and body, the head pointed along the forward motion of the body during straight walking. The head pitch/roll component appeared to compensate for vertical and horizontal acceleration of the head rather than orienting to the tilt of the GIA or anticipating it. When turning corners, subjects walked on a 50-cm radius over two steps or on a 200-cm radius in five to seven steps. Maximum centripetal accelerations in sharp turns were ca.0.4 g, which tilted the GIA ca.21 degrees with regard to the heading. This was anticipated by a roll tilt of the head of up to 8 degrees. The eyes rolled 1-1.5 degrees and moved down into the direction of linear acceleration during the tilts of the GIA. Yaw head deviations moved smoothly through the turn, anticipating the shift in lateral body trajectory by as much as 25 degrees. The trunk did not anticipate the change in trajectory. Thus, in contrast to straight walking, the tilt axes of the head and the GIA tended to align during turns. Gaze was stable in space during the slow phases and jumped forward in saccades along the trajectory, leading it by larger angles when the angular velocity of turning was greater. The anticipatory roll head movements during turning are likely to be utilized to overcome inertial forces that would destabilize balance during turning. The data show that compensatory eye, head, and body movements stabilize gaze during straight walking, while orienting mechanisms direct the eyes, head, and body to tilts of the GIA in space during turning.

Non-NASA Center↗
DOE OSTI2017

Comparison of Virtual Oscillator and Droop Control: Preprint

Virtual oscillator control and droop control are two techniques that can be used to ensure synchronization and power sharing of parallel inverters in islanded operation. VOC relies on the implementation of non-linear Van der Pol oscillator equations in the control system of the inverter, acting upon the time-domain instantaneous inverter current and terminal voltage. On the other hand, DC explicitly computes active and reactive power produced by the inverter and relies on limited bandwidth low-pass filters. Even though both methods can be engineered to produce the same steady-state characteristics, their dynamic performances are significantly different. This paper presents analytical and experimental results that aim to compare both methods. It is shown that VOC is inherently faster and enables minimizing the circulating currents. The results are verified using three 120V, 1kW inverters.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Photoelectron Spectroscopy and Theoretical Studies of Anion-pi Interactions: Binding Strength and Anion Specificity

Proposed in theory and confirmed to exist, anion–π interactions have been recognized as new and important non-covalent binding forces. Despite extensive theoretical studies, numerous crystal structural identifications, and a plethora of solution phase investigations, intrinsic anion–π interaction strengths that are free from complications of condensed phases’ environments, have not been directly measured in the gas phase. Herein we present a joint photoelectron spectroscopic and theoretical study on this subject, in which tetraoxacalix[2]arene[2]triazine 1, an electron-deficient and cavity self-tunable macrocyclic was used as a charge-neutral molecular host to probe its interactions with a series of anions with distinctly different shapes and charge states (spherical halides Cl⁻, Br⁻, I⁻, linear thiocyanate SCN⁻, trigonal planar nitrate NO₃⁻, pyramidic iodate IO₃⁻, and tetrahedral sulfate SO₄²⁻). The binding energies of the resultant gaseous 1:1 complexes (1•Cl⁻,1•Br⁻, 1•I⁻, 1•SCN⁻, 1•NO₃⁻, 1•IO₃⁻ and 1•SO₄²⁻) were directly measured experimentally, exhibiting substantial non-covalent interactions with pronounced anion specific effects. The binding strengths of Cl⁻, NO₃⁻, IO₃⁻ with 1 are found to be strongest among all singly charged anions, amounting to ca. 30 kcal/mol, but only about 40% of that between 1 and SO₄²⁻. Quantum chemical calculations reveal that all anions reside in the center of the cavity of 1 with anion–π binding motif in the complexes’ optimized structures, where 1 is seen to be able to self-regulate its cavity structure to accommodate anions of different geometries and three-dimensional shapes. Electron density surface and natural bond orbital charge distribution analysis further support anion–π binding formation. The calculated binding energies of the anions and 1 nicely reproduce the experimentally estimated electron binding energy increase. This work illustrates that size-selective photoelectron spectroscopy combined with theoretical calculations represent a powerful technique to probe intrinsic anion–π interactions and has potential to provide quantitative guest-host molecular binding strengths and unravel fundamental insights in specific anion recognitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Investigations into the Ternary NaF-KF-UF4 Salt System – Phase A

A knowledge gap exists in the data and understanding of fresh fuel salt and irradiated multicomponent fuel salt systems thermophysical properties. Quantifying these properties is necessary for the design and construction of test reactors, as well as the licensing of future commercial molten-salt reactors. To facilitate thermal property determination on a proposed fuel salt composition for Seaborg Technologies, several samples containing depleted uranium tetrafluoride (UF4), sodium fluoride (NaF), and potassium fluoride (KF) were blended, and a melt temperature analysis was performed. From the melting temperature analysis, it was determined that sample Seaborg-7, a ternary salt composition of 26.4UF4-24.7KF-48.9NaF (mol%), was very near a ternary eutectic point. Therefore, thermal properties such as melting temperature, salt stability, density, heat capacity, thermal diffusivity, and viscosity were experimentally determined on the Seaborg-7 salt. These measurements document the baseline properties of fresh fuel salt as a function of temperature, where future experiments on irradiated fuel salt will provide a holistic perspective on the change of thermophysical properties during reactor operations. Several precision instruments were used to collect property data, and instrument calibrations and data collection were performed and documented in a standardized and reproducible manner with meticulous detail. This process ensured that the measurement procedures and resulting data can readily be duplicated elsewhere. The Seaborg-7 salt was shown to be stable at temperatures up to 900°C, as no mass change was observed upon repeated heating and cooling. The peak melting temperature was determined to be 547°C (557°C endset). The enthalpy of fusion (??H?_fus^o) was determined to be 167.5 ± 2.7 J/g while the enthalpy of crystallization (??H?_c^o) was determined to be -147.8 ± 13.3 J/g. In addition to the eutectic melting peak, upon heating, several pre eutectic peaks were observed, occurring at 470°C (onset) and 499°C (peak). Specific heat capacity measurements showed a slightly increasing trend with respect to temperature in the solid phase, while the liquid-specific heat capacity showed a somewhat flat trend with an average value of 106.1 ± 1.24 J/mol·K between 600 to 800°C. Three independent trials using the Seaborg-7 salt determined the density to be ?(T) = 4.908 – 0.000363·T(°C), validated between 32 to 200°C, and ?(T) = 4.808 – 0.00113·T(°C), validated between ~575 to 850°C. Thermal diffusivity was determined for the liquid state and is represented by the linear equation y = 0.1581 + 0.000207·T(°C) between 550 to 850°C. The viscosity was determined from 600 to 800°C and is represented by the exponential fit equation, ? (mPa·s) = 736.58e^(-0.006·T(°C)). This report documents the conclusion of fuel salt thermophysical property measurements for the Seaborg SPP, Phase A project.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗
DOE OSTI2016

Development of a CO 2 Chemical Sensor for Downhole CO 2 Monitoring in Carbon Sequestration

Geologic storage of carbon dioxide (CO 2 ) has been proposed as a viable means for reducing anthropogenic CO 2 emissions. The means for geological sequestration of CO 2 is injection of supercritical CO 2 underground, which requires the CO 2 to remain either supercritical, or in solution in the water/brine present in the underground formation. However, there are aspects of geologic sequestration that need further study, particularly in regards to safety. To date, none of the geologic sequestration locations have been tested for storage integrity under the changing stress conditions that apply to the sequestration of very large amounts of CO 2 . Establishing environmental safety and addressing public concerns require widespread monitoring of the process in the deep subsurface. In addition, studies of subsurface carbon sequestration such as flow simulations, models of underground reactions and transports require a comprehensive monitoring process to accurately characterize and understand the storage process. Real-time information about underground CO 2 movement and concentration change is highly helpful for: (1) better understanding the uncertainties present in CO 2 geologic storage; (2) improvement of simulation models; and (3) evaluation of the feasibility of geologic CO 2 storage. Current methods to monitor underground CO 2 storage include seismic, geoelectric, isotope and tracer methods, and fluid sampling analysis. However, these methods commonly resulted low resolution, high cost, and the inability to monitor continuously over the long time scales of the CO 2 storage process. A preferred way of monitoring in-situ underground CO 2 migration is to continuous measure CO 2 concentration change in brine during the carbon storage process. An approach to obtain the real time information on CO 2 concentration change in formation solution is highly demanded in carbon storage to understand the CO 2 migration subsurface and to answer the public safety problem. The objective of the study is to develop a downhole CO 2 sensor that can in-situ, continuously monitor CO 2 concentration change in deep saline. The sensor is a Severinghaus-type CO 2 sensor with small size, which renders it can be embedded in monitoring well casing or integrated with pressure/temperature transducers, enabling the development of “smart” wells. The studies included: (1) prepare and characterize metal-oxide electrodes. Test the electrodes response to pH change. Investigate different ions and brine concentration effects on the electrode’s performance. Study the stability of the electrode in brine solution; (2) fabricate a downhole CO 2 sensor with the metal-oxide electrodes prepared in the laboratory. Test the performance of the CO 2 sensor in brine solutions. Study high pressure effects on the performance of the sensor; (3) design and conduct CO 2 /brine coreflooding experiments with the CO2 sensor. Monitor CO 2 movement along the core and test the performance of the sensor in coreflooding tests. Develop a data acquisition system that can digitize the sensor’s output voltage. Our completed research has resulted in deep understanding of downhole CO 2 sensor development and CO 2 monitoring in CO 2 storage process. The developed downhole CO 2 sensor included a metal-oxide electrode, a gas-permeable membrane, a porous steel cup, and a bicarbonate-based internal electrolyte solution. Iridium oxide-based electrode was prepared and used for preparation the CO 2 sensor. The prepared iridium oxide-based electrode displayed a linearly response to pH change. Different factors such as different ions and ions concentration, temperature, and pressure effects on the electrode performance on pH response were investigated. The results indicated that the electrode exhibited a good performance even in high salt concentration of produced water. To improve the electrode performance under high pressure, IrO 2 nanoparticles with the particle size in the range of 1-2 nm were prepared and electrodeposited on stainless steel substrate by cyclic voltammetry. It was observed that the thin film of iridium oxide was formed on the substrate surface and such iridium oxide-based electrode displayed excellent performance under high pressure for longer term. A downhole CO 2 sensor with the iridium oxide-based electrode was prepared. The working principle of the CO 2 sensor is based on the measurement of the pH change of the internal electrolyte solution caused by the hydrolysis of CO 2 and then determination of the CO 2 concentration in water. The prepared downhole CO 2 sensor had the size of diameter of 0.7 in. and length of 1.5 in. The sensor was tested under the pressures of 500 psi, 2,000 psi, and 3,000 psi. A linear correlation was observed between the sensor potential change and dissolved CO 2 concentration in water. The response time of the CO 2 sensor was in the range of 60-100 minutes. Further tests indicated that the CO 2 sensor exhibited good reproducibility under high pressure. A CO 2 /brine coreflooding system was constructed to simulate the real-world CO 2 storage process. The prepared downhole CO 2 sensor was loaded in the system to monitor CO 2 movement during CO 2 /brine coreflooding test. The results indicated that the sensor could detect CO 2 movement in the tests. Further studies showed that the sensor could be recovered by brine flooding after CO 2 /brine flushed the core. The results of the coreflooding tests demonstrated that the sensor had potential application for CO 2 monitoring in carbon sequestration. A data acquisition system for the downhoe CO 2 sensor was developed and coded. The system converted the sensor output signal into digital data and transported the data from downhole to wellhead surface. The data acquisition system was tested and evaluated in the laboratory with the prepared sensor for data collection.

36 MATERIALS SCIENCE↗
DOE OSTI2015

SU-E-T-675: Remote Dosimetry with a Novel PRESAGE Formulation

Purpose: 3D-gel dosimetry provides high-resolution treatment validation; however, scanners aren’t widely available. In remote dosimetry, dosimeters are shipped out from a central base institution to a remote site for irradiation, then shipped back for scanning and analysis, affording a convenient service for treatment validation to institutions lacking the necessary equipment and resources. Previous works demonstrated the high-resolution performance and temporal stability of PRESAGE. Here the newest formulation is investigated for remote dosimetry use. Methods: A new formulation of PRESAGE was created with the aim of improved color stability post irradiation. Dose sensitivity was determined by irradiating cuvettes on a Varian Linac (6MV) from 0–15Gy and measuring change in optical density at 633nm. Sensitivity readings were tracked over time in a temperature control study to determine long-term stability. A large volume study was performed to evaluate the accuracy for remote dosimetry. A 1kg dosimeter was pre-scanned, irradiated on-site with an 8Gy 4field box treatment, post-scanned and shipped to Princess Margaret Hospital for remote reading on an identical scanner. Results: Dose sensitivities ranged from 0.0194–0.0295 ΔOD/(Gy*cm)—similar to previous formulations. Post-irradiated cuvettes stored at 10°C retained 100% initial sensitivity over 5 days and 98.6% over 10 weeks while cuvettes stored at room temperature fell to 95.8% after 5 days and 37.4% after 10 weeks. The immediate and 5-day scans of the 4field box dosimeter data was reconstructed, registered to the corresponding eclipse dose-distribution, and compared with analytical tools in CERR. Immediate and 5-day scans looked visually similar. Line profiles revealed close agreement aside from a slight elevation in dose at the edge in the 5-day readout. Conclusion: The remote dosimetry formulation exhibits excellent temporal stability in small volumes. While immediate and 5-day readout scans of large volume dosimeters show promising agreement, further development is required to reduce an apparent time dependent edge elevation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Final Report for Project: Impacts of stratification and non-equilibrium winds and waves on hub-height winds

This project used a combination of turbulence-resolving large-eddy simulations, single-column modeling (where turbulence is parameterized), and currently available observations to improve, assess, and develop a parameterization of the impact of non-equilibrium wave states and stratification on the buoy-observed winds to establish reliable wind data at the turbine hub-height level. Analysis of turbulence-resolving simulations and observations illuminates the non-linear coupling between the atmosphere and the undulating sea surface. This analysis guides modification of existing boundary layer parameterizations to include wave influences for upward extrapolation of surface-based observations through the turbine layer. Our surface roughness modifications account for the interaction between stratification and the effects of swell’s amplitude and wavelength as well as swell’s relative motion with respect to the mean wind direction. The single-column version of the open source Weather and Research Forecasting (WRF) model (Skamarock et al., 2008) serves as our platform to test our proposed planetary boundary layer parameterization modifications that account for wave effects on marine atmospheric boundary layer flows. WRF has been widely adopted for wind resource analysis and forecasting. The single column version is particularly suitable to development, analysis, and testing of new boundary layer parameterizations. We utilize WRF’s single-column version to verify and validate our proposed modifications to the Mellor-Yamada-Nakanishi-Niino (MYNN) boundary layer parameterization (Nakanishi and Niino, 2004). We explore the implications of our modifications for two-way coupling between WRF and wave models (e.g.,Wavewatch III). The newly implemented parameterization accounting for marine atmospheric boundary layer-wave coupling is then tested in three-dimensional WRF simulations at grid sizes near 1 km. These simulations identify the behavior of simulated winds at the wind plant scale. Overall project conclusions include; In the presence of fast-moving swell (significant wave height Hs = 6.4 m, and phase speed cp = 18 ms -1 ), the atmospheric boundary layer grows more rapidly when waves propagate opposite to the winds compared to when winds and waves are aligned. Pressure drag increases by nearly a factor of 2 relative to the turbulent stress for the extreme case where waves propagate at 180° compared to the pressure gradient forcing. Net wind speed reduces by nearly 15% at hub-height for the 180°-case compared to the 0°-case, and turbulence intensities increase by nearly a factor of 2. These impacts diminish with decreasing wave age; Stratification increases hub height wind speeds and increases the vertical shear of the mean wind across the rotor plane. Fortuitously, this stability-induced enhanced shear does not influence turbulence intensity at hub height, but does increase (decrease) turbulence intensity below (above) hub height. Increased stability also increases the wave-induced pressure stress by ~ 10%; Off the East Coast of the United States during Coupled Boundary Layers Air-Sea Transfer - Low Wind (CBLAST-Low), cases with short fetch include thin stable boundary layers with depths of only a few tens of meters. In the coastal zone, the relationship between the mean wind and the surface fiction velocity (u*(V )) is significantly related to wind direction for weak winds but is not systematically related to the air sea difference of virtual potential temperature, δθv; since waves generally propagate from the south at the Air-Sea Interaction Tower (ASIT) tower, these results suggest that under weak wind conditions waves likely influence surface stress more than stratification does; and Winds and waves are frequently misaligned in the coastal zone. Stability conditions persist for long duration. Over a four year period, the Forschungsplattformen in Nord- und Ostsee Nr. 1 (FINO1) tower (a site with long fetch) primarily experienced weakly-unstable conditions, while stability at the ASIT tower (with a larger influence of offshore winds) experiences a mix of both unstable and stable conditions, where the summer months are predominantly stable. Wind-wave misalignment likely explains the large scatter in observed non-dimensional surface roughness under swell-dominated conditions. Andreas et al.’s (2012) relationship between u* and the 10-m wind speed under predicts the increased u* produced by wave-induced pressure drag produced by misaligned winds and waves. Incorporating wave-state (speed and direction) influences in parameterizations improves predictive skill. In a broad sense, these results suggest that one needs information on winds, temperature, and wave state to upscale buoy measurements to hub-height and across the rotor plane. Our parameterization of wave-state influences on surface drag has been submitted for inclusion in the next publicly available release. In combination, our project elucidates the impacts of two important physical processes (non-equilibrium wind/waves and stratification) on the atmosphere within which offshore turbines operate. This knowledge should help guide and inform manufacturers making critical decisions surrounding design criteria of future turbines to be deployed in the coastal zone. Reductions in annually averaged hub height wind speed error using our new wave-state-aware surface layer parameterization are relatively modest. However since wind turbine power production depends on the wind speed cubed, the error in estimated power production is close to 5%; which is significant and can substantially impact wind resource assessment and decision making with regards to the viability of particular location for a wind plant location. For a single 30-hour forecast, significant reductions in wind speed prediction errors can yield substantially improved wind power forecast skill, thereby mitigating costs and/or increasing revenue through improved; forecasting for maintenance operations and planning; day-ahead forecasting for power trading and resource allocation; and short-term forecasting for dispatch and grid balancing.

17 WIND ENERGY↗
DOE OSTI2014

On the feasibility of comprehensive high-resolution 3D remote dosimetry

Purpose: This study investigates the feasibility of remote high-resolution 3D dosimetry with the PRESAGE®/Optical-CT system. In remote dosimetry, dosimeters are shipped out from a central base institution to a remote institution for irradiation, then shipped back to the base institution for subsequent readout and analysis. Methods: Two nominally identical optical-CT scanners for 3D dosimetry were constructed and placed at the base (Duke University) and remote (Radiological Physics Center) institutions. Two formulations of PRESAGE® (SS1, SS2) radiochromic dosimeters were investigated. Higher sensitivity was expected in SS1, which had higher initiator content (0.25% bromotrichloromethane), while greater temporal stability was expected in SS2. Four unirradiated PRESAGE® dosimeters (two per formulation, cylindrical dimensions 11 cm diameter, 8.5–9.5 cm length) were imaged at the base institution, then shipped to the remote institution for planning and irradiation. Each dosimeter was irradiated with the same simple treatment plan: an isocentric 3-field “cross” arrangement of 4 × 4 cm open 6 MV beams configured as parallel opposed laterals with an anterior beam. This simple plan was amenable to accurate and repeatable setup, as well as accurate dose modeling by a commissioned treatment planning system (Pinnacle). After irradiation and subsequent (within 1 h) optical-CT readout at the remote institution, the dosimeters were shipped back to the base institution for remote dosimetry readout 3 days postirradiation. Measured on-site and remote relative 3D dose distributions were registered to the Pinnacle dose calculation, which served as the reference distribution for 3D gamma calculations with passing criteria of 5%/2 mm, 3%/3 mm, and 3%/2 mm with a 10% dose threshold. Gamma passing rates, dose profiles, and color-maps were all used to assess and compare the performance of both PRESAGE® formulations for remote dosimetry. Results: The best agreements between the Pinnacle plan and dosimeter readout were observed in PRESAGE® formulation SS2. Under 3%/3 mm 3D gamma passing criteria, passing rates were 91.5% ± 3.6% (SS1) and 97.4% ± 2.2% (SS2) for immediate on-site dosimetry, 96.7% ± 2.4% (SS1) and 97.6% ± 0.6% (SS2) for remote dosimetry. These passing rates are well within TG119 recommendations (88%–90% passing). Under the more stringent criteria of 3%/2 mm, there is a pronounced difference [8.0 percentage points (pp)] between SS1 formulation passing rates for immediate and remote dosimetry while the SS2 formulation maintains both higher passing rates and consistency between immediate and remote results (differences ≤ 1.2 pp) at all metrics. Both PRESAGE® formulations under study maintained high linearity of dose response (R{sup 2} > 0.996) for 1–8 Gy over 14 days with response slope consistency within 4.9% (SS1) and 6.6% (SS2), and a relative dose distribution that remained stable over time was demonstrated in the SS2 dosimeters. Conclusions: Remote 3D dosimetry was shown to be feasible with a PRESAGE® dosimeter formulation (SS2) that exhibited relative temporal stability and high accuracy when read off-site 3 days postirradiation. Characterization of the SS2 dose response demonstrated linearity (R{sup 2} > 0.998) over 14 days and suggests accurate readout over longer periods of time would be possible. This result provides a foundation for future investigations using remote dosimetry to study the accuracy of advanced radiation treatments. Further work is planned to characterize dosimeter reproducibility and dose response over longer periods of time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Stable electrospray signal on a microfabricated glass chip with three-dimensional open edge and tiered depth geometries

This paper presents the microfabrication and performance of a three-dimensional electrospray ionization (ESI) emitter tip made from glass, which achieves stable current signals important for chemical analysis. Our novel microfabrication process and custom-built signal conditioning hardware provides the advantage of providing accurate features and steady signals. The fabrication process relies on standard microfabrication techniques (i.e., deposition, photolithography, and wet etching). This fabrication method involves the novel application of two layers of positive and negative photoresists in addition to Parafilm® wax tape. Open edge and tiered depth details were successfully created from a multilayer planar mask. This is a benefit for integrated miniaturized and microfluidic systems that often require micro features for their functionality but relatively large millimeter size features for their physical periphery. We demonstrate the fundamental performance of electrospray with this microfluidic chip. The emitter tip was fixed on a linear axis stage with high resolution (10 μm) to finely control the tip distance from a metal counter electrode plate. Further, a custom printed circuit board system was built to safely control four voltages applied to the microchip ports from a single high voltage power supply. To readily form the electrospray, non-aqueous solvents were used for their low viscosity and a constant voltage of +2.7 kV was applied to the sheath electrospray microchannel. The liquid being sprayed was 80/20 (v/v) methanol/acetonitrile with 0.1% acetic acid in the sheath microchannel and with ammonium acetate (10–40 mM) in its remaining microchannels. The electrospray signal was measured in response to varying the distance (1.4 to 1.6mm) between the electrospray emitter tip and a metal counter electrode plate in addition to the varying concentration of the background electrolyte, ammonium acetate. Stable and repeatable electrospray signal showed linear relationships with emitter tip distance and concentration (r 2 ≥ 0.95).

42 ENGINEERING↗
DOE OSTI2016

Investigation of the phase relations in the U-Al-Ge ternary system: Influence of the Al/Ge substitution on the properties of the intermediate phases

The phase relations within the U-Al-Ge ternary system were studied for two isothermal sections, at 673 K for the whole Gibbs triangle and at 1173 K for the concentration range 25–100 at% U. The identification of the phases, their composition ranges and stability were determined by x-ray powder diffraction, scanning electron microscopy coupled to energy dispersive spectroscopy and differential thermal analysis. The tie-lines and the solubility domains were determined for the U-Ge and U-Al binaries, the UAl{sub 3}-UGe{sub 3} solid-solution and for the unique ternary intermediate phase U{sub 3}Al{sub 2−x}Ge{sub 3+x}. The experimental isopleth section of the pseudo-binary UAl{sub 3}-UGe{sub 3} reveals an isomorphous solid solution based on the Cu{sub 3}Au-type below the solidus. The U{sub 3}Al{sub 2−x}Ge{sub 3+x} solid solution extends for −0.1≤x≤1.35 and −0.2≤x≤1.5 at 673 K and 1173 K respectively. It crystallizes in the I-centered tetragonal symmetry. The reciprocal lattice of several compositions of the U{sub 3}Al{sub 2−x}Ge{sub 3+x} solid solution was examined by electron diffraction at room temperature, revealing the presence of a c-glide plane. Their crystal structure was refined by single crystal x-ray diffraction suggesting an isomorphous solid solution best described with the non-centrosymmetric space group I4cm in the paramagnetic domain. The magnetic measurements confirm the ferromagnetic ordering of the solid solution U{sub 3}Al{sub 2−x}Ge{sub 3+x} with an increase of Tc with the Al content. The thermal variation of the specific heat bear out the magnetic transitions with some delocalized character of the uranium 5f electrons. - Graphical abstract: The phase relations within the U-Al-Ge ternary system were experimentally assessed for two isothermal sections, at 673 K for the whole Gibbs triangle and at 1173 K for the concentration range 25–100 at% U. A complete UAl{sub 3}-UGe{sub 3} solid-solution based on the Cu{sub 3}Au-type forms below the solidus. A unique ternary phase showing a large homogeneity domain, U{sub 3}Al{sub 2−x}Ge{sub 3+x} for −0.1≤x≤1.35 and −0.2≤x≤1.5 at 673 K and 1173 K respectively has been evidenced. It is best described with the non-centrosymmetric space group I4cm above room temperature. A linear increase of the ferromagnetic ordering is observed with the Al content. - Highlights: • Isothermal sections of the U-Al-Ge system were investigated for 673 K and 1173 K. • An isomorphous solid-solution UAl{sub 3}-UGe{sub 3} forms for the whole composition range. • U{sub 3}Al{sub 2−x}Ge{sub 3+x} the unique ternary phase to form exists for a large homogeneity domain. • U{sub 3}Al{sub 2−x}Ge{sub 3+x} is best described in I4cm space group above room temperature. • The ferromagnetic transition of U{sub 3}Al{sub 2−x}Ge{sub 3+x} linearly increases with the Al content.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗