SEARCH · Engineering Papers
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.
This dataset contains all well logs from Utah FORGE well 16A(78)-32. This includes the mud log, Sanvean Technologies logs, and Schlumberger logs. Please see the file descriptions below for information about each log.
Shale formations are attractive prospects due to their potential in oil and gas production. Some of the largest shale formations in the mainland US, such as the Tuscaloosa Marine Shale (TMS), have reserves estimated to be around 7 billion barrels. Despite their huge potential, shale formations present major concerns for drilling operators. These prospects have unique challenges because of all their alteration and incompatibility issues with drilling and completion fluids. Most shale formations undergo numerous chemical and physical alterations, making their interaction with the drilling and completion fluid systems very complex to understand. In this study, a high-pressure, high-temperature (HPHT) drilling simulator was used to mimic real time drilling operations to investigate the performance of inhibitive drilling fluid systems in two major shale formations (Eagle Ford Shale and Tuscaloosa Marine Shale). A series of drilling experiments using the drilling simulator and clay swelling tests were conducted to evaluate the drilling performance of the KCl drilling fluid and cesium formate brine systems and their effectiveness in minimizing drilling concerns. Cylindrical cores were used to mimic vertical wellbores. It was found that the inhibitive muds systems (KCl and cesium formate) provided improved drilling performance compared to conventional fluid systems. Among the inhibitive systems, the cesium formate brine showed the best drilling performances due to its low swelling rate and improved drilling performance.
Corneal penetration is a key rate limiting step in the bioavailability of topical ophthalmic formulations that incorporate poorly permeable drugs. Recent advances have greatly aided the ocular delivery of such drugs using colloidal drug delivery systems. Ribavirin, a poorly permeable BCS class-III drug, was incorporated in bioadhesive multiple W/O/W microemulsion (ME) to improve its corneal permeability. The drug-loaded ME was evaluated regarding its physical stability, droplet size, PDI, zeta potential, ultrastructure, viscosity, bioadhesion, in vitro release, transcorneal permeability, cytotoxicity, safety and ocular tolerance. Our ME possessed excellent physical stability, as it successfully passed several cycles of centrifugation and freeze–thaw tests. The formulation has a transparent appearance due to its tiny droplet size (10 nm). TEM confirmed ME droplet size and revealed its multilayered structure. In spite of the high aqueous solubility and the low permeability of ribavirin, this unique formulation was capable of sustaining its release for up to 24 h and improving its corneal permeability by 3-fold. The in vitro safety of our ME was proved by its high percentage cell viability, while its in vivo safety was confirmed by the absence of any sign of toxicity or irritation after either a single dose or 14 days of daily dosing. Our ME could serve as a vehicle for enhanced ocular delivery of drugs with different physicochemical properties, including those with low permeability.
During the Surveyor program spacecraft were successfully landed at five widely separated lunar locations. Recent computer simulations of each landing have provided more comprehensive data on the mechanical properties of the lunar surface than have been obtained previously by this method of analysis. Results show that the variations in surface bearing pressure observed at the various lunar sites are probably due to surface slope effects and do not necessarily indicate differences in soil properties at these sites. Estimates of cohesion at two sites give almost identical results and further support the conclusion that the soil properties at all sites are probably very similar. Surface pressures that resist horizontal (plowing) motion are largely due to cohesion, and density and gravitational contributions are small. It is concluded that the lunar surface bearing strength is essentially zero at the surface and, for zero surface slope, increases with penetration depth at a rate of 1.87 (plus or minus 0.33) N/cu cm. The cohesion of the lunar soil is estimated to be between 0.11 and 0.17 N/sq cm.
Geothermal is a promising source of renewable energy with almost zero emissions. However, there remains untapped geothermal potential around the world. This is largely because geothermal projects have high development costs and high resource-value uncertainty, and returns on initial investments can be slow to materialize. Increasing drilling efficiency of geothermal wells is one way to decrease development costs, as drilling accounts for up to 50% of the upfront costs of a geothermal power project. Currently, little quantitative information exists about how fast geothermal drilling occurs on a global scale. This paper was originally intended to be an extension of Frone and Boyd's 2018 report on geothermal drilling rates in California and Nevada. Our project includes global data and considers additional measures, such as flat time. Data for the project was limited to published papers and drilling reports that are publicly available online. The goal of our project was to establish a baseline geothermal drilling rate that is representative of current global geothermal drilling practices can be used to gauge the impact of future improvements in geothermal drilling technology. We recorded the number of days from beginning to end of a project, the number of days spent drilling, and the number of days considered to be flat time. We found the average global rate at which a drilling project is completed (including non-drilling activities) to be 160 ft/day, and the average drilling rate (including only time during which active drilling occurred) to be 360 ft/day. Our data shows no clear trend in rate change from 2000 to 2017 and significant variability in drilling rates between countries and within countries.
Thermal runaway of lithium-ion batteries can involve various types of failure mechanisms each with their own unique characteristics. Using fractional thermal runaway calorimetry and high-speed radiography, the response of three different geometries of cylindrical cell (18650, 21700, and D-cell) to different abuse mechanisms (thermal, internal short circuiting, and nail penetration) are quantified and statistically examined. Correlations between the geometry of cells and their thermal behavior are identified, such as increasing heat output per amp-hour (kJ Ah −1 ) of cells with increasing cell diameter during nail penetration. High-speed radiography reveals that the rate of thermal runaway propagation within cells is generally highest for nail penetration where there is a relative increase in rate of propagation with increasing diameter, compared to thermal or internal short-circuiting abuse. For a given cell model tested under the same conditions, a distribution of heat output is observed with a trend of increasing heat output with increased mass ejection. Finally, internal temperature measurements using thermocouples embedded in the penetrating nail are shown to be unreliable thus demonstrating the need for care when using thermocouples where the temperature is rapidly changing. All data used in this manuscript are open access through the NREL and NASA Battery Failure Databank.
Passive detection of special nuclear material (SNM) is challenging due to its inherently low rate of spontaneous emission of penetrating radiation, the relative ease of shielding, and the fluctuating and frequently overwhelming background. Active interrogation, the use of external radiation to increase the emission rate of characteristic radiation from SNM, has long been considered to be a promising method to overcome those challenges. Current active-interrogation systems that incorporate radiography tend to use bremsstrahlung beams, which can deliver high radiation doses. Low-energy ion-driven nuclear reactions that produce multiple monoenergetic photons may be used as an alternative. The 12 C ( p , p ' ) 12 C reaction is one such reaction that could produce large yields of highly penetrating 4.4- and 15.1-MeV gamma rays. This reaction does not directly produce neutrons below the approximately 19.7 MeV threshold, and the 15.1-MeV gamma-ray line is well matched to the photofission cross section of 235 U and 238 U . In this article, we report the measurements of thick-target gamma-ray yields at 4.4 and 15.1 MeV from the 12 C ( p , p ' ) 12 C reaction at proton energies of 19.5, 25, and 30 MeV. Measurements are made with two 3 -in. EJ-309 cylindrical liquid scintillation detectors and thermoluminescent dosimeters placed at 0 ° and 90 ° , with an additional 1.5 -in. Na I ( Tl ) cylindrical scintillation detector at 0 ° . We estimate the highest yields of the 4.4- and 15.1-MeV gamma rays of 1.65 × 10 10 and 4.47 × 10 8 sr - 1 μ C - 1 at a proton energy of 30 MeV, respectively. The yields in all experimental configurations are greater than in a comparable deuteron-driven reaction that produces the same gamma-ray energies— 11 B ( d , n γ ) 12 C . However, a significant increase of the neutron radiation dose accompanies the proton energy increase from 19.5 to 30 MeV.
The author has identified the following significant results. A new device, called a Spectral Attenuation Board, has been designed and tested, which enables ERTS-1 sea truth collection teams to monitor the attenuation depths of three colors continuously, as the board is being towed behind a boat. The device consists of a 1.2 x 1.2 meter flat board held below the surface of the water at a fixed angle to the surface of the water. A camera mounted above the water takes photographs of the board. The resulting film image is analyzed by a micro-densitometer trace along the descending portion of the board. This yields information on the rate of attenuation of light penetrating the water column and the Secchi depth. Red and green stripes were painted on the white board to approximate band 4 and band 5 of the ERTS MSS so that information on the rate of light absorption by the water column of light in these regions of the visible spectrum could be concurrently measured. It was found that information from a red, green, and white stripe may serve to fingerprint the composition of the water mass. A number of these devices, when automated, could also be distributed over a large region to provide a cheap method of obtaining valuable satellite ground truth data at present time intervals.
Slides for Oxidation Activities, oxidation rate, penetration/lathing, and strength after oxidation work. Also includes recent and current oxidation studies, background for rate behavior, strength behavior, density profile, and penetration measurements, oxidation resistant graphite coating development, and irradiated and unirradiated graphite oxidation rate behavior.
During average solar wind flow conditions at 1 AU, ionization rates of interstellar neutrals that penetrate into the inner heliosphere are dominated by charge exchange with solar wind protons for H atoms, and by photoionization for He atoms. During occurrences of strong, coronal mass ejection (CME)-driven interplanetary shock waves near 1 AU, electron impact ionization can make substantial, if not dominating, contributions to interstellar neutral ionization rates in the regions downstream of the shocks. However, electron impact ionization is expected to be relatively less important with increasing heliocentric distance because of the decrease in electron temperature. Ulysses encountered many CME-driven shocks during its journey to and beyond Jupiter, and in addition, encountered a number of strong corotating interaction region (CIR) shocks. These shocks generally occur only beyond approximately 2 AU. Many of the CIR shocks were very strong rivalling the Earth's bow shock in electron heating. We have compared electron impact ionization rates calculated from electron velocity distributions measured downstream from CIR shocks using the Ulysses SWOOPS experiment to charge-exchange rates calculated from measured proton number fluxes and the photoionization rate estimated from an assumed solar photon spectrum typical of solar maximum conditions. We find that, although normally the ratio of electron-impact ionization rates to charge-exchange (for H) and to photoionization (for He) rates amounts to only about one and a few tens of percent, respectively, downstream of some of the stronger CIR shocks they amount to more than 10% and greater than 100%, respectively.
The spray performance of a fuel injection system applicable for use in main combustion chamber of an oxidizer-rich staged combustion (ORSC) cycles is presented. The experimental data reported here include mean drop size and drop size distribution, spray cone half-angle, and momentum rate (directly related to spray penetration). The maximum entropy formalism, MEF, method to predict drop size distribution is applied and compared to the experimental data. Geometric variables considered include the radius of the injector inlet orifice plate through which oxidizer flows (&) and the exposed length from the fuel inlet to the injector exit plane (L2). Operating conditions that were varied include the liquid mass flow rate and air mass flow rate. For orifices B and C there is a significant dependence of D3Z on both the air and liquid mass flow rates, as well as on L2. For the A orifice, the momentum rate of the air flow appears to exceed a threshold value above which a constant D32 is obtained. Using the MEF method, a semi-analytical process was developed to model the spray distribution using two input parameters (q = 0.4 and Dso). The momentum rate of the spray is directly related to the air and liquid mass flow rates. The cone half angle of the spray ranges from 25 to 17 degrees. The data resulting from this project will eventually be used to develop advanced rocket systems.
One of the most important aspects that need to be addressed to increase solar energy penetration is the power ramp-rate control. In weak grids such as the one found in Puerto Rico, it is important to smooth power fluctuations caused by the intermittence of passing clouds. In this work, a novel power ramp-rate control strategy is proposed. Additionally, a comparison with some of the most common power ramp-rate control methods is performed using a proposed model and real solar radiation data from the Coto Laurel photovoltaic power plant located in Ponce, Puerto Rico. The proposed model was validated using one-year real data from Coto Laurel. The power ramp-rate control methods were compared in real-time simulations using the OP5700 from Opal-RT Technologies considering power ramp rate fluctuations, power ramp-rate violations, fluctuations in the state-of-charge, among other indicators. Moreover, the proposed power ramp-rate control strategy, called predictive dynamic smoothing was explained and compared. Results indicate that the predictive dynamic smoothing produced a considerably reduced Levelized Cost of Storage compared to other power ramp-rate control methods and provided a higher lifetime expectancy for lithium batteries.
There is considerable variation between subduction zones in the largest characteristic earthquake within each zone. Assuming that coupling between downgoing and upper plates is directly related to characteristic earthquake size, tests for correlations between variation in coupling and other physical features of subduction zones are conducted: the lateral extent and penetration depth of Benioff zones, age of subducting lithosphere, convergence rate, and back-arc spreading. Using linear multivariate regression, coupling is correlated with two variables: convergence rate and lithosphere age. Secondary correlations within the data set are penetration depth versus lithosphere age, and lateral extent versus convergence rate. Taken together, the observed correlations suggest a simple qualitative model where convergence rate and lithosphere age determine the horizontal and sinking rates, respectively, of slabs: these parameters influence the seismic coupling in the subduction zone. In the limit of a fast sinking rate and slow convergence rate, back-arc spreading occurs and thereby appears to be a passive process.
Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.
Abstract Many existing and future tokamaks with tungsten divertors operate, or will operate, with low- Z impurity seeding, but the direct effect of these seeded impurities on tungsten Scrape-off-Layer (SOL) transport has not been explored in detail. This paper reports on a DIII-D experiment designed to test how tungsten divertor leakage from the Small-Angle Slot V-Shaped, tungsten-coated divertor is impacted by neon seeding at a variety of injection rates and poloidal injection locations. Measurements from the experiment show an inverse relationship between the neon injection rate and the tungsten core penetration factor. Interpretive modeling is performed with a combination of the SOLPS-ITER and DIVIMP codes to assess the underlying tungsten behavior. The modeling results show that the reduction in tungsten divertor leakage is driven by both an increase in the divertor collisionality as well as a reduction in the ion temperature gradient near the divertor target. Collisions between low- Z impurities and tungsten impurities are found to have a significant impact on the tungsten SOL transport, such that ignoring the low- Z impurity collisional effects on the tungsten transport can result in an overestimate of the divertor leakage by an order-of-magnitude. Given the importance of these localized interactions, neon seeding from the closed, slot-like divertor has a clear advantage in being able to reduce tungsten divertor leakage without the high levels of neon core contamination that occur when seeding from other poloidal locations.
This report summarizes the progress in developing a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose of this effort is to introduce to a user community—involving researchers, regulators, and industry—a generic, reliable numerical model for the prediction of aerosol transport while accounting for potential deposition and plugging of the leak paths to model spent nuclear fuel (SNF) release from postulated Stress Corrosion Cracks in canisters under storage or transportation. Current work focuses on expanding the model to predict aerosol (and gas) flow through complex microchannel (nozzle) geometries in addition to the rectangular and cylindrical geometries, as we approach more realistic Stress Corrosion Crack conditions. In this regard, a divergent nozzle geometry, used at Sandia National Laboratories (SNL) for testing aerosol release and retention (Durbin et al. 2021), was added to the model. The model was then validated with blowdown data for the particular microchannel, from experiments conducted at SNL. The report also presents preliminary non-benchmarked aerosol penetration fraction and mass flow rate for a monodisperse 10-micron (AED) particle concentration of 1.7e-08 kg/m 3 being released through the divergent microchannel geometry from the SNL’s aerosol experimental tank setup.Future work will involve validating the aerosol model and updating the model’s Graphical User Interface (GUI) to include the divergent nozzle geometry. This is expected to help stakeholders perform quick and easy first principles calculations without needing to understand the underlying MATLAB script.
Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.
A model of the zodiacal cloud is used to evaluate the effects of rotational bursting on interplanetary dust particles caused by the interaction of solar radiation pressure with the irregular surface geometry of the particles. The effectiveness of this mechanism is evaluated using several values for the effective moment arm and tensile strength of the particles. Due to repeated bursting, it is shown that particles spiraling in from the asteroid belt by the Poynting-Robertson effect do not penetrate inside 1.2 AU. Spin rate stabilization under the effect of magnetic damping is found to depend critically on the initial heliocentric distance of the particle and on its radius.