Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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,009 records · Page 56

Understanding the origin of Tabula Rasa-induced defects in n-type Cz c-Si: The case of nitrogen atmosphere

Phosphorus-doped Czochralski-grown silicon (Cz-Si) has been gaining market share in the large-scale manufacturing of high-efficiency silicon (Si)-based photovoltaic (PV) devices thanks to higher carrier lifetimes than their boron-doped counterpart. However, the fabrication of n-type Cz-Si based solar cells often requires process steps with much higher temperatures and longer times than p-type Silicon. Defect interaction with the high temperatures during such processes tend to be detrimental to the n-type Cz-Si carrier lifetime, therefore limiting the final device efficiency. Short thermal anneals before cell processing, known as Tabula Rasa (TR), have been proposed to mitigate the thermally induced lifetime degradation during n-type Cz-Si solar cell fabrication. This work thoroughly investigates the defects responsible for the lifetime degradation after TR in a N 2 atmosphere treatment. We use temperature-injection-dependent lifetime spectroscopy and the thickness variation method to decouple the effects of TR treatment in the bulk and the surface of the n-type Cz-Si wafers. Using the defect parameter contour mapping (DPCM), we identify the defect energy level (E t ) and the capture cross-section ratio (k) of the most likely process-induced defect, which aligns with previously proposed Si vacancy-associated defects. The DPCM reveals that these vacancy-associated defects have a shallow energy level E t - E v ~0.13 eV and very efficient electron capture cross section k~600. Unexpectedly, the bulk degradation due to vacancy defects in the volume of the wafer, is accompanied by a significant increase in the surface recombination as well. Through evaluating the surface recombination velocity temperature- and injection dependence, we show that after TR, at room temperature and for an injection level of 10 15 cm -3 , in a wafer passivated with a-Si:H(i) the surface recombination dominates the overall lifetime response. Here we hypothesize that the near surface vacancy-associated bulk defects play a role in lowering the electron diffusion current into the a-Si:H(i) from the c-Si(n) reducing the field-effect passivation.

14 SOLAR ENERGY↗

Kepler Mission: A Wide-FOV Photometer Designed to Determine the Frequency of Earth-Size and Larger Planets Around Solar-like stars

The first step in discovering the extent of life in our galaxy is to determine the number of terrestrial planets in the habitable zone (HZ). The Kepler Mission is designed around a 0.95 m aperture Schmidt-type telescope with an array of 42 CCDs designed to continuously monitor the brightness of 100,000 solar-like stars to detect the transits of Earth-size and larger planets. The photometer is scheduled to be launched into heliocentric orbit in 2007. Measurements of the depth and repetition time of transits provide the size of the planet relative to the star and its orbital period. When combined with ground-based spectroscopy of these stars to fix the stellar parameters, the true planet radius and orbit scale, hence the position relative to the HZ are determined. These spectra are also used to discover the relationships between the characteristics of planets and the stars they orbit. In particular, the association of planet size and occurrence frequency with stellar mass and metallicity will be investigated. At the end of the four year mission, hundreds of terrestrial planets should be discovered in and near the HZ of their stars if such planets are common. A null result would imply that terrestrial planets in the HZ occur in less than 1% of the stars and that life might be quite rare. Based on the results of the current doppler-velocity discoveries, detection of a thousand giant planets is expected. Information on their albedos and densities of those giants showing transits will be obtained.

Borucki, William↗

Impact Damage Tolerance of Composite Laminates with Through-The-Thickness Stitches

The ability of through-the-thickness stitches to contain damage during a low-velocity impact event and increase the residual strength of stitched panels was investigated in this research. The impact response, spread of interlaminar delaminations, dent depth, surface damage, and static residual strength after impact were studied for carbon-epoxy specimens fabricated from HTS40 TohoTenax standard modulus fibers, stitched together with Vectran 1200 denier thread and infused with API 1078 resin after through-the-thickness stitching. Three different stitch patterns were used to explore the ability to arrest impact damage during and after impact. Simply supported boundary conditions were maintained during the impact testing. Non-destructive evaluations were performed using ultrasonic C-scans and X-ray computed tomography (CT) imaging to determine the shape, size, and location of delaminations. Results indicate that while the dynamic response during the impact event was almost the same for the unstitched specimens and all stitch patterns considered, the extent of delamination and the compression strength after impact varied greatly. For both the 60 J and 80 J impact energies, the delamination area was significantly less for the stitched specimens than for the unstitched specimens, but the range of delamination areas among the stitch patterns was much larger for the lower impact energy than for the greater impact energy. Similarly, while the presence of stitching influenced the compression after impact strength, the strength values for all stitch patterns were very similar. These results are a step toward quantifying the influence of through-the-thickness stitching.

Ranatunga, Vipul↗

Adaptive Data-based Predictive Control for Short Take-off and Landing (STOL) Aircraft

Data-based Predictive Control is an emerging control method that stems from Model Predictive Control (MPC). MPC computes current control action based on a prediction of the system output a number of time steps into the future and is generally derived from a known model of the system. Data-based predictive control has the advantage of deriving predictive models and controller gains from input-output data. Thus, a controller can be designed from the outputs of complex simulation code or a physical system where no explicit model exists. If the output data happens to be corrupted by periodic disturbances, the designed controller will also have the built-in ability to reject these disturbances without the need to know them. When data-based predictive control is implemented online, it becomes a version of adaptive control. The characteristics of adaptive data-based predictive control are particularly appropriate for the control of nonlinear and time-varying systems, such as Short Take-off and Landing (STOL) aircraft. STOL is a capability of interest to NASA because conceptual Cruise Efficient Short Take-off and Landing (CESTOL) transport aircraft offer the ability to reduce congestion in the terminal area by utilizing existing shorter runways at airports, as well as to lower community noise by flying steep approach and climb-out patterns that reduce the noise footprint of the aircraft. In this study, adaptive data-based predictive control is implemented as an integrated flight-propulsion controller for the outer-loop control of a CESTOL-type aircraft. Results show that the controller successfully tracks velocity while attempting to maintain a constant flight path angle, using longitudinal command, thrust and flap setting as the control inputs.

Barlow, Jonathan Spencer↗

Investigation of ion-beam machining methods for replicated x-ray optics

The final figuring step in the fabrication of an optical component involves imparting a specified contour onto the surface. This can be expensive and time consuming step. The recent development of ion beam figuring provides a method for performing the figuring process with advantages over standard mechanical methods. Ion figuring has proven effective in figuring large optical components. The process of ion beam figuring removes material by transferring kinetic energy from impinging neutral particles. The process utilizes a Kaufman type ion source, where a plasma is generated in a discharge chamber by controlled electric potentials. Charged grids extract and accelerate ions from the chamber. The accelerated ions form a directional beam. A neutralizer outside the accelerator grids supplies electrons to the positive ion beam. It is necessary to neutralize the beam to prevent charging workpieces and to avoid bending the beam with extraneous electro-magnetic fields. When the directed beam strikes the workpiece, material sputters in a predicable manner. The amount and distribution of material sputtered is a function of the energy of the beam, material of the component, distance from the workpiece, and angle of incidence of the beam. The figuring method described here assumes a constant beam removal, so that the process can be represented by a convolution operation. A fixed beam energy maintains a constant sputtering rate. This temporally and spatially stable beam is held perpendicular to the workpiece at a fixed distance. For non-constant removal, corrections would be required to model the process as a convolution operation. Specific figures (contours) are achieved by rastering the beam over the workpiece at varying velocities. A unique deconvolution is performed, using series-derivative solution developed for the system, to determine these velocities.

Drueding, Thomas W.↗

Growth of organized flow coherent motions within a single-stream shear layer: 4D-PTV measurements

Abstract This study investigates the evolution of a single-stream shear layer (SSSL) originating from a wall boundary layer past a backward-facing step. Utilizing a time-resolved 3D-Particle Tracking Velocimetry (4D-PTV) technique, we track the trajectories of fluorescent particles to gain insight into the flow characteristics of the SSSL. A compact water tunnel facility ( $$\textrm{Re}_\tau =1\,240$$ Re τ = 1 240 ) is fabricated to obtain an SSSL with a perpendicular slow entrainment stream past the separation edge. A hybrid interpolation approach that combines ensemble binning and Gaussian weighting is implemented to derive minimally filtered mean and instantaneous lower- and higher-order flow field parameters. Spanwise-dominant coherent motion accompanied by finer flow scales is observed to grow due to flow entrainment through “nibbling” actions of small-scale vortices, “engulfing” by large-scale vortices, and vortex pairing events. Furthermore, the non-zero-speed stream edge grows relatively faster than the zero-speed stream edge, showing a strong asymmetry in mixing composition across a mixing layer. The SSSL reaches self-similarity at a streamwise distance of $$\approx 55\,\theta _{0}$$ ≈ 55 θ 0 , where $$\theta _0$$ θ 0 is the initial momentum thickness from the separation edge, i.e., considerably shorter than reported in previous studies. A literature comparison of growth rate parameters raises intriguing questions regarding a potential inclusive growth scaling unifying the free shear layers. A turbulent kinetic energy (TKE) budget analysis reveals a negative production region immediately downstream of the separation edge attributed to a large positive streamwise gradient of streamwise velocity. In the self-similar region, the phase-averaged flow mapping demonstrates a larger concentration of turbulence production rate around the outer edges of spanwise vortices, specifically at the intersection of braids and vortices. Furthermore, a spatial separation exists in the regions of peak production and dissipation rates within the vortex core region favoring dissipation. The braids exhibit a larger concentration of turbulence diffusion rates, indicating their function as a conduit for exchanging turbulence between neighboring coherent motions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kepler Mission: A Discovery-Class Mission Designed to Determine the Frequency of Earth-Size and Larger Planets Around Solar-Like Stars

The first step in discovering the extent of life in our galaxy is to determine the number of terrestrial planets in the habitable zone (HZ). The Kepler Mission is designed around a 0.95 in aperture Schmidt-type telescope with an array of 42 CCDs designed to continuously monitor the brightness of 100,000 solar-like stars to detect the transits of Earth-size and larger planets. The photometer is scheduled to be launched into heliocentric orbit in 2007. Measurements of the depth and repetition time of transits provide the size of the planet relative to the star and its orbital period. When combined with ground-based spectroscopy of these stars to fix the stellar parameters, the true planet radius and orbit scale, hence the position relative to the HZ are determined. These spectra are also used to discover the relationships between the characteristics of planets and the stars they orbit. In particular, the association of planet size and occurrence frequency with stellar mass and metallicity will be investigated. At the end of the four year mission, hundreds of terrestrial planets should be discovered in and near the HZ of their stars if such planets are common. Extending the mission to six years doubles the expected number of Earth-size planets in the HZ. A null result would imply that terrestrial planets in the HZ occur in less than 1% of the stars and that life might be quite rare. Based on the results of the current Doppler-velocity discoveries, detection of a thousand giant planets is expected. Information on their albedos and densities of those giants showing transits will be obtained.

William Borucki↗

Design and Experimental Demonstration of an Additive Manufactured Smart Oxy-Methane Burner for High Pressure and Supercritical Combustion

Pressurized oxy combustion-based systems can improve efficiency by recovering latent heat of the steam in the Flue Gas and achieving 90% CO2 capture. In addition, the novel Directly Heated Supercritical Carbon Dioxide (DH-SCO2) power cycles can achieve high thermal efficiencies and provide nearly full carbon capture. Additionally, due to the reduction of flue gas at higher pressure, smaller system size and capital cost reductions are also possible. Recent thermodynamic analysis of the DH-SCO2 cycle performed by the UTEP research team shows that combustion conditions in the vicinity of 300 bar pressure and 1000-1400 K temperature allow for relatively high system efficiencies while operating within the limit of available combustor materials. However, the realization of a directly heated supercritical power cycle requires combustion systems to operate in supercritical conditions and at temperature far below the blowout limit of conventional flames (above 1500 K). The thermodynamic properties along with the combustion properties and kinetics are unexplored at such conditions. Additionally, the interaction of the supercritical environment with energy components is also unknown. High-pressure combustion tests are performed using a smart burner at some intermediate pressure ranges (<20 bar) to help minimize these knowledge gaps. The knowledge obtained from the high-pressure test will assist in understanding the combustion chamber pressurization mechanism, ignition and flame behavior at the elevated pressure. The obtained data will act as a systematic first step in testing at higher pressures of 100 and 300 bar pressures.The primary purpose of this Dissertation is to demonstrate the operability of a low NOx smart burner for oxy-methane combustion at high pressure (< 20 bar) and scalable up to supercritical conditions. A shear co-axial smart burner is designed with a real-time temperature monitoring capability to understand the burner face interaction at high-pressure conditions. In addition, the burner has four independent injection ports to allow the independent injection of fuel and dilution gases in the combustor. The maximum operating capability of the burner is 575 kWth. The burner was fabricated using a Laser Powder Bed Fusion process with Nickel Alloy 718. A powder removal technology was invented comprising ultrasonic vibration, liquid nitrogen exposure and media blasting to remove powders from internal channels of the burner. The burner operability tests are performed in the high-pressure combustor and the swirl combustor. The high-pressure combustor was used to investigate the burner operability and thermal soak back at different pressurized conditions. The experimental tests in high-pressure combustor up to 275 kWth input resulted in 16.5 bar chamber pressure and 198°C thermal soaks back to the burner. The burner was capable of providing the required thermal input within a 3% deviation range. In addition, soot formation occurred at high-pressure tests. The swirl combustor was used to observe the flame stability of the burner. Flame lift-off was observed for jet velocities above 450 m/s. Additionally, lift-off decreased for low co-flow velocities. CO2 dilution experiments showed increased flame instabilities for all conditions above 50% dilution ratios. At lower thermal inputs, partial flame blow-offs occurred for dilution ratios above 50%. All conditions significantly reduced the flame temperature and increased the flame lift-off height. Finally, a 2nd generation AM smart burner was designed using the knowledge from the 1st generation burner experiments. The 2nd generation burner incorporated two sets of swirlers with 0.9 swirl no. A cooling system was also designed for long-duration tests at higher pressures. The thermal input and division of the burner's power are kept the same as the 1st generation burner. The burner is to be fabricated using nickel-alloy 718 for high-pressure handling capability. The design can sustain at high-pressure conditions up to 100 bar.

Islam, Md Nawshad Arslan↗

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First Annual Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Sensor Performance Model of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. Feasibility of the sensor was initially investigated with analytical derivations and with COMSOL RF Module computer simulations of resonant frequency spectral shift due to uniform load. We also investigated the mechanical integrity of the flowmeter’s membrane through analytical modelling and COMSOL Structural Mechanics Module computer simulations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the material’s yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure. Using results from models, we have developed an initial design for a microwave K-band sensor. A cylindrical resonator prototype was fabricated from brass for the initial tests. The external dimensions of the cavity are matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed consisting of a piping Tee with a bulkhead WR-42 microwave waveguide installed in a leak-proof assembly. A microwave waveguide circulator was installed in the setup to suppress the effect of reflections at the cavity entrance by increasing the isolation between the input and the output port. Preliminary spectral characterization of cavity spectral response was performed with a portable PXIe chassis microwave VNA with a custom GUI. Preliminary dry tests of the transducer response were conducted with a set of calibrated weights. Transducer frequency shift was shown to be monotonically increasing with increasing pressure. The next steps will involve investigation of the transducer performance for water flow sensing.

42 ENGINEERING↗

Steps toward determination of the size and structure of the broad-line region in active galactic nuclei. 8: An intensive HST, IUE, and ground-based study of NGC 5548

We present the data and initial results from a combined HST/IUE/ground-based spectroscopic monitoring campaign on the Seyfert 1 galaxy NGC 5548 that was undertaken in order to address questions that require both higher temporal resolution and higher signal-to-noise ratios than were obtained in our previous multiwavelength monitoring of this galaxy in 1988-89. IUE spectra were obtained once every two days for a period of 74 days beginning on 14 March 1993. During the last 39 days of this campaign, spectroscopic observations were also made with the HST Faint Object Spectrograph (FOS) on a daily basis. Ground-based observations, consisting of 165 optical spectra and 77 photometric observations (both CCD imaging and aperture photometry), are reported for the period 1992 October to 1993 September, although much of the data are concentrated around the time of the satellite-based program. These data constitute a fifth year of intensive optical monitoring of this galaxy. In this contribution, we describe the acquisition and reduction of all of the satellite and ground-based data obtained in this program. We describe in detail various photometric problems with the FOS and explain how we identified and corrected for various anomalies. During the HST portion of the monitoring campaign, the 1350 A continuum flux is found to have varied by nearly a factor of two. In other wavebands, the continuum shows nearly identical behavior, except that the amplitude of variability is larger at shorter wavelengths, and the continuum light curves appear to show more short time-scale variability at shorter wavelengths. The broad emission lines also vary in flux, with amplitudes that are slightly smaller than the UV continuum variations and with a small time delay relative to the UV continuum. On the basis of simple time-series analysis of the UV and optical continuum and emission line light curves, we find (1) that the ultraviolet and optical continuum variations are virtually simultaneous, with any lag between the 1350 A continuum and the 5100 A continuum amounting to less than about one day; (2) that the variations in the highest ionization lines observed, He II lambda 1640 and N V lambda 1240, lag behind the continuum variations by somewhat less than 2 days, and (3) that the velocity field of the C IV-emitting region is not dominated by radial motion. The results on the C IV velocity field are preliminary and quite uncertain, but there are some weak indications that the emission-line (wings absolute value of Delta upsilon is greater than or equal to 3000 km/s) respond to continuum variations slightly more rapidly than does the core. The optical observations show that the variations in the broad H beta line flux follow the continuum variations with a time lag of around two weeks, about twice the lag for Ly alpha and C IV, as in our previous monitoring campaign on this same galaxy. However, the lags measured for Ly alpha, C IV, and H Beta are each slightly smaller than previously determined. We confirm two trends reported earlier, namely (1) that the UV/optical continuum becomes 'harder' as it gets brighter, and (2) that the highest ionization emission lines have the shortest lags, thus indicating radial ionization stratificatin of a broad-line region that spans over an order of magnitude range in radius.

Korista, K.↗

Atmospheric Dust Modeling from Meso to Global Scales with the Online NMMB/BSC-Dust Model Part 2: Experimental Campaigns in Northern Africa

The new NMMB/BSC-Dust model is intended to provide short to medium-range weather and dust forecasts from regional to global scales. It is an online model in which the dust aerosol dynamics and physics are solved at each model time step. The companion paper (Perez et al., 2011) develops the dust model parameterizations and provides daily to annual evaluations of the model for its global and regional configurations. Modeled aerosol optical depth (AOD) was evaluated against AERONET Sun photometers over Northern Africa, Middle East and Europe with correlations around 0.6-0.7 on average without dust data assimilation. In this paper we analyze in detail the behavior of the model using data from the Saharan Mineral dUst experiment (SAMUM-1) in 2006 and the Bodele Dust Experiment (BoDEx) in 2005. AOD from satellites and Sun photometers, vertically resolved extinction coefficients from lidars and particle size distributions at the ground and in the troposphere are used, complemented by wind profile data and surface meteorological measurements. All simulations were performed at the regional scale for the Northern African domain at the expected operational horizontal resolution of 25 km. Model results for SAMUM-1 generally show good agreement with satellite data over the most active Saharan dust sources. The model reproduces the AOD from Sun photometers close to sources and after long-range transport, and the dust size spectra at different height levels. At this resolution, the model is not able to reproduce a large haboob that occurred during the campaign. Some deficiencies are found concerning the vertical dust distribution related to the representation of the mixing height in the atmospheric part of the model. For the BoDEx episode, we found the diurnal temperature cycle to be strongly dependant on the soil moisture, which is underestimated in the NCEP analysis used for model initialization. The low level jet (LLJ) and the dust AOD over the Bodélé are well reproduced. The remaining negative AOD bias (due to underestimated surface wind speeds) can be substantially reduced by decreasing the threshold friction velocity in the model.

photometers↗

Uranium Oxide Elutriation at Higher Flow Rates and Longer Times

This work supports the recovery of high-assay low-enriched uranium (i.e., U enriched in 235 U to >5% and <20%) from used nuclear fuel. One step includes elutriating triuranium octoxide (U 3 O 8 ) from a fluidized bed of alumina (Al 2 O 3 ) that fills the main reaction chamber. The experiments outlined in this report assessed the effects of the time and flow rate on the elutriation of U 3 O 8 powder from the reaction vessel. Testing at Oak Ridge National Laboratory (ORNL) was performed in a one-fourth scale experimental apparatus based on a pilot-scale system, which is the Material Recovery Pilot Plant (MRPP) at Idaho National Laboratory (INL). The experimental system comprises a supply of pressurized nitrogen (N 2 ), a gas heater, the reaction vessel, a knock-out pot, and a water trap where elutriated particles are collected for sampling. An Al 2 O 3 fluidized bed is contained within the reaction vessel. A bed of Al 2 O 3 with a nominal size of 180 µm and U 3 O 8 with a nominal size of 7 µm was fluidized at 2 different flow rates so the U 3 O 8 physically separated from the bed. The apparatus was simple and had no valves or sharp bends downstream from the reactor. The flow exited the reactor through a side port at the top. The flow passed through a knock-out vessel before the flow was passed through a bubbler to capture the U 3 O 8 . Seven tests were completed with U 3 O 8 using an elutriating gas of N 2 heated to 200°C. Tests were performed with bed loadings that comprise 25.0 g of U 3 O 8 , or 3.6 wt% of the bed mass. Except for one test performed at 12 SLPM, the gas flow rate through the bed for the other experiments was set at 10 SLPM as measured by the mass flow meter. Linear velocities through the bed were calculated as 28 and 34 cm/s, or 0.9 and 1.1 ft·s -1 . The duration of the tests ranged from 20 to 180 min, and one test had a flow cycled between 0 and 5 SLPM every 5 min for 40 min. The testing successfully demonstrated elutriation of U 3 O 8 from the fluidized bed as up to a 100% of the U 3 O 8 was transferred out of the bed. The removal of U 3 O 8 had a logarithmic dependence on flow rate. The design changes made the apparatus easier to clean, and the mass balances were better than those used in previous studies. Some elutriated Al 2 O 3 collected in the tubing that connected the reactor to the knock-out pot, but the amount of Al 2 O 3 was minimal for 10 SLPM. The flow at 12 SLPM showed considerable transport of Al 2 O 3 bed material into the reaction-vessel headspace and the transfer line to the knock-out pot.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows

New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.

42 ENGINEERING↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗

Absolute Position of Targets Measured Through a Chamber Window Using Lidar Metrology Systems

Lidar is a useful tool for taking metrology measurements without the need for physical contact with the parts under test. Lidar instruments are aimed at a target using azimuth and elevation stages, then focus a beam of coherent, frequency modulated laser energy onto the target, such as the surface of a mechanical structure. Energy from the reflected beam is mixed with an optical reference signal that travels in a fiber path internal to the instrument, and the range to the target is calculated based on the difference in the frequency of the returned and reference signals. In cases when the parts are in extreme environments, additional steps need to be taken to separate the operator and lidar from that environment. A model has been developed that accurately reduces the lidar data to an absolute position and accounts for the three media in the testbed air, fused silica, and vacuum but the approach can be adapted for any environment or material. The accuracy of laser metrology measurements depends upon knowing the parameters of the media through which the measurement beam travels. Under normal conditions, this means knowledge of the temperature, pressure, and humidity of the air in the measurement volume. In the past, chamber windows have been used to separate the measuring device from the extreme environment within the chamber and still permit optical measurement, but, so far, only relative changes have been diagnosed. The ability to make accurate measurements through a window presents a challenge as there are a number of factors to consider. In the case of the lidar, the window will increase the time-of-flight of the laser beam causing a ranging error, and refract the direction of the beam causing angular positioning errors. In addition, differences in pressure, temperature, and humidity on each side of the window will cause slight atmospheric index changes and induce deformation and a refractive index gradient within the window. Also, since the window is a dispersive media, the effect of both phase and group indices have to be considered. Taking all these factors into account, a method was developed to measure targets through multiple regions of different materials and produce results that are absolute measurements of target position in three-dimensional space, rather than simply relative position. The environment in which the lidar measurements are taken must be broken down into separate regions of interest and each region solved for separately. In this case, there were three regions of interest: air, fused silica, and vacuum. The angular position of the target inside the chamber is solved using only phase index and phase velocity, while the ranging effects due to travel from air to glass to vacuum/air are solved with group index and group velocity. When all parameters are solved simultaneously, an absolute knowledge of the position of each target within an environmental chamber can be derived. Novel features of this innovation include measuring absolute position of targets through multiple dispersive and non-dispersive media, deconstruction of lidar raw data from a commercial off-the-shelf unit into reworkable parameters, and use of group velocities to reduce range data. Measurement of structures within a vacuum chamber or other harsh environment, such as a furnace, may now be measured as easily as if they were in an ambient laboratory. This analysis permits transformation of the raw data into absolute spatial units (e.g., mm). This technique has also been extended to laser tracker, theodolite, and cathetometer measurements through refractive media.

Kubalak, David↗

Tidal Turbine Test, Downeast Turbines, July 12, 2021

Downeast Turbines tested a tidal turbine prototype with novel rotor/channel system and lateral effluent discharge apparatus (LEDA), during five days of testing in the flume at Alden Lab. Three days of testing (July 12-14, 2021) focused on turbine power metrics of torque and rpm, which were low, and then two days of follow-up testing (July 27-28, 2021) focused on LEDA performance metrics of pressure differential and rates of volumetric flow, with encouraging results. Next step is to to characterize, and even optimize, configurations of the LEDA, using 3D-CFD as a helpful tool, to refine its shape and explore its limits of performance as a means of effluent discharge that augments performance of an instream turbine. An improved configuration of the LEDA will be re-combined with the rotor/channel system of the turbine prototype, and ultimately, the rotor will be re-sized (enlarged), to better match the LEDA's performance capabilities in drawing through a rate of volumetric flow. This submission is Downeast Turbines' Post Access Report for the test event. It includes the files described here (next below), and several reference links. "Downeast TEAMER-Post-Access-Report...docx" is a document file containing the report. "Appendices A, B, and C" are included in this file, and so are "Figures #1-7." "Appendix D - Test Data Workbooks.zip" is an archive file containing all post access data (raw data tables, calculating tables, and graphs), presented in fourteen Excel workbooks as described in the report. "Appendix E - Post Access Figures.zip" is an archive file containing "Figures #8-52," of the report.

16 TIDAL AND WAVE POWER↗

Novel and effective surface passivation for high efficiency n- and p-type Silicon solar cell

The project objective was to develop a novel Si surface passivation method using chalcogens, sulfur (S) and/or selenium (Se), as passivating elements, to withstand industry-standard high temperature contacting and metallization schemes for p-type Si based passivated emitter and rear contact (p-PERC) solar cells. The back surface passivation of PERC cells has been improved drastically with the invention and successful application of an Al 2 O 3 passivation layer. However, the front n + diffused junction surface is still poorly passivated by the standard amorphous silicon nitride (SiNx) anti-reflection coating (ARC) layer. This project sought to address the passivation challenges of both front n+ emitter and undiffused p-Si back surface. Improved p-PERC solar cell performance with open circuit voltage (V OC ) > 680 mV and efficiency of 22% were targeted to validate superior defect passivation properties as compared to standard SiO 2 / Al 2 O 3 passivation. During this project, we systematically investigated process-structure-properties-performance relationships of this novel advanced defect passivation approach. The S/Se passivation was carried out by reacting industrial Czochralski (Cz) Si wafers in H 2 S and H 2 Se gases in a chemical vapor deposition (CVD) reactor at temperatures up to 700°C. After an exhaustive optimization of the process parameters (temperature, time, and gas concentration), we established an optimized process and demonstrated extremely low surface recombination velocities (SRVs) of 1.5 cm/s and 8 cm/s on n-type and p-type Si, respectively, by S-passivation. In-depth surface and interface characterization were performed using soft x-ray and photoelectron spectroscopies (XPS, UPS, XES), combined with capacitance-voltage-frequency (C-V-f) measurements, to decipher the surface chemical/electronic structure and interface defect state densities. These measurements provided critical understanding of the defect passivation mechanism and elucidated the presence of surface S-Si bonds, a reduction of surface dipoles, and low interface state densities (D it ) < 10 11 cm -2 ev -1 . We also found that the Se-passivation is inferior to the S-passivation (by at least one order of magnitude in SRV). Application of the optimized S-passivation to the n+ diffused emitter surface led to a low surface recombination current density, J0 ≈ 40 fA/cm 2 (~ 1/4 of the industry-standard SiNx-passivation), and high implied V OC (686 mV) in p-PERC solar cell structures. The S-passivation process also was found to improve the bulk quality of the p-type Si, better than the SiO 2 or Al 2 O 3 passivation processes. After successful demonstration of efficient passivation of Si surface defects by S, we extensively studied the air, thermal, and illumination stability of the passivation structure. S-passivation itself degrades in air due to competing reactions with moisture and oxygen to form oxides, which can be eliminated by a SiNx capping layer (also acting as a anti-reflective coating). After SiNx process optimization, we demonstrated illumination and thermally stable S-passivation with SRV < 5 cm/s and J 0 < 80 fA/cm 2 . These enhancements in Si passivation, incorporated into p-PERC cells, achieved an efficiency of 19.93% with V OC = 649 mV, using manufacturing metallization and contacting schemes. The low cell performance (cell V OC is much less than the implied V OC = 686 mV, anticipated from surface passivation) was identified due to degradation of S-passivation during the metal firing step (out-diffusion of S from the Si interface to the SiNx surface). The S-passivation of Si surfaces shows significant promise with excellent passivation quality, essential for high performance (high V OC , high efficiency) solar cells. Integration of this innovative defect passivation into devices, however, demands further development of the capping layer, low temperature (<700°C) metallization process, and/or engineering of advanced device structures. Surface passivation-dominated advanced Si solar cells, such as tunnel oxide passivated contacts and Si heterojunctions, are increasingly of interest due to their high-performance potential and will have a growing photovoltaic market share in the near future.

14 SOLAR ENERGY↗