SEARCH · Engineering Papers

Results for “injection height”

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.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Effects of oscillating shock impingement on the mixing of a gaseous jet in a Mach 3 airstream

Oscillating shocks of different frequencies were impinged on a 15 deg downstream angled, sonic helium jet injected into a Mach 3 airflow. The oscillating shocks were produced by the normal injection of a steady liquid jet from an opposing wall into the Mach 3 airflow. The freestream Reynolds number per centimeter was 500,000 producing a boundary layer 1.6 times the helium injector diameter in thickness. Vertical helium concentration profiles were taken at nine transverse stations at an axial station 90 helium jet diameters downstream. Helium concentration was measured using an aspirating hot film probe with an integral thermocouple. Oscillating shock frequencies of 60 kHz, 120 kHz, and 240 kHz were studied, compared with the estimated frequencies of 60 kHz of the largest eddies in the approach boundary layer. Visualization using nanoshadowgraph photographs, showed large turbulent structures in all cases. In addition, there were clear changes in eddy size with changing frequency visible on the nanoshadowgraphs. Impingement of the oscillating shock upstream of the helium jet had no measurable effect on the mixing. Impingement downstream of the helium injector produced measurable changes in the mixing achieved. Increasing the frequency of the oscillating shock decreased the maximum concentration and the penetration height of the maximum concentration. There were additional measurable effects on the area of a typical concentration contour with increased oscillation frequency.

Thomas, R. H.↗

DOE OSTI2021

High Temperature Stability of Aqueous Foams for Potential Application in Enhanced Geothermal System (EGS)

Enhanced Geothermal Systems (EGS) are typically constructed by injecting high-pressure water into the deep hot dry rocks (HDR) under carefully controlled conditions to create new or re-open existing fractures, which usually utilizes an immense quantity of water. Alternatively, aqueous foam-based fracturing fluids may offer potential advantages over conventional water-based fracturing fluid, including reduced consumption of water and environmental impact. Although foam-based fracturing has shown promising results in oil and gas industries, its applicability to EGS is unknown. Foams are complex mixtures of liquids with gaseous phases that are thermodynamically unstable. Their stability will decrease over time due to liquid drainage, bubble coarsening, and coalescence. Therefore, it is essential to stabilize foam fluids at high temperatures for EGS related applications such as fracking of HDRs.This paper is focused on investigating high-temperature stability of N2-in-water or N2 foams containing different surfactants and stabilizing agents. Foam stability was characterized by the half-life, which is defined as the time a foam decreases to 50% of its original height due to drainage. Surfactants with various concentrations dissolved in deionized water were investigated, including alfa olefin sulfonate (AOS), sodium dodecyl sulfonate (SDS), TergitolTM (NP – 40), and cetyltrimethylammonium chloride (CTAC). Also, the effects of other stabilizing agents, including guar gum, bentonite clay, crosslinker, graphene oxide (GO) dispersion on foam stability were tested. Initial results showed that foam stability decreased dramatically as temperature increased. On the other hand, foams became relatively more stable as pressure increased. Certain stabilizing agents, such as guar gum and crosslinker, can enhance foam stability at a wide range of temperatures. Our results indicate that it may be possible to obtain stable foams at high-temperature, high-pressure conditions with appropriate stabilizing agents.

Thakore, Virensinh↗

DOE OSTI2023

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗

DOE OSTI2015

Vertical equilibrium, energetics, and star formation rates in magnetized galactic disks regulated by momentum feedback from supernovae

Recent hydrodynamic (HD) simulations have shown that galactic disks evolve to reach well-defined statistical equilibrium states. The star formation rate (SFR) self-regulates until energy injection by star formation feedback balances dissipation and cooling in the interstellar medium (ISM), and provides vertical pressure support to balance gravity. In this paper, we extend our previous models to allow for a range of initial magnetic field strengths and configurations, utilizing three-dimensional, magnetohydrodynamic (MHD) simulations. We show that a quasi-steady equilibrium state is established as rapidly for MHD as for HD models unless the initial magnetic field is very strong or very weak, which requires more time to reach saturation. Remarkably, models with initial magnetic energy varying by two orders of magnitude approach the same asymptotic state. In the fully saturated state of the fiducial model, the integrated energy proportions E{sub turb}:E{sub th}:δE{sub mag}: E-bar {sub mag} are 0.35:0.39:0.15:0.11, while the proportions of midplane support P{sub turb}:P{sub th}:δΠ{sub mag}: Π-bar {sub mag} are 0.49:0.18:0.18:0.15. Vertical profiles of total effective pressure satisfy vertical dynamical equilibrium with the total gas weight at all heights. We measure the “feedback yields” η{sub c}≡P{sub c}/Σ{sub SFR} (in suitable units) for each pressure component, finding that η{sub turb}∼3.5−4 and η{sub th}∼1.1−1.4 are the same for MHD as in previous HD simulations, and δη{sub mag}∼1.3−1.5. These yields can be used to predict the equilibrium SFR for a local region in a galaxy based on its observed gas and stellar surface densities and velocity dispersions. As the ISM weight (or dynamical equilibrium pressure) is fixed, an increase in η from turbulent magnetic fields reduces the predicted Σ{sub SFR} by ∼20−30% relative to the HD case.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Effects of edge-localized electron cyclotron current drive on edge-localized mode suppression by resonant magnetic perturbations in DIII-D

Abstract According to recent DIII-D experiments (Logan et al 2024 Nucl. Fusion 64 014003), injecting edge localized electron cyclotron current drive (ECCD) in the counter-plasma-current (counter- I p ) direction reduces the n = 3 resonant magnetic perturbation (RMP) current threshold for edge-localized mode (ELM) suppression, while co- I p ECCD during the suppressed ELM phase causes a back transition to ELMing. This paper presents nonlinear two-fluid simulations on the ECCD manipulation of edge magnetic islands induced by RMP using the TM1 code. In the presence of a magnetic island chain at the pedestal-top, co- I p ECCD is found to decrease the island width and restore the initially degraded pedestal pressure when its radial deposition location is close to the rational surface of the island. With a sufficiently strong co- I p ECCD current, the RMP-driven magnetic island can be healed, and the pedestal pressure fully recovers to its initial ELMing state. On the contrary, counter- I p ECCD is found to increase the island width and further reduce the pedestal pressure to levels significantly below the peeling-ballooning-mode limited height, leading to even stationary ELM suppression. These simulations align with the results from DIII-D experiments. However, when multiple magnetic island chains are present at the pedestal-top, the ECCD current experiences substantial broadening, and its effects on the island width and pedestal pressure become negligible. Further simulations reveal that counter- I p ECCD enhances RMP penetration by lowering the penetration threshold, with the degree of reduction proportional to the amplitude of ECCD current. For the ∼1 MW ECCD in DIII-D, the predicted decrease in the RMP penetration threshold for ELM suppression is approximately 20%, consistent with experimental observations. These simulations indicate that edge-localized ECCD can be used to either facilitate RMP-driven ELM suppression or optimize the confinement degradation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

NASA's Space Launch System Begins Integration, Stacking in Preparation for Artemis I Launch

The Artemis era of human lunar exploration is nearing take-off as NASA’s new super heavy-lift launch vehicle, the Space Launch System (SLS), begins stack-ing and integration operations in mid-2020 at Kennedy Space Center (KSC) in Florida. With a planned upgrade path to progressively more powerful vehicles and availability in crew and cargo configurations, SLS provides a unique and flexible launch solution to send crew, large-scale infrastructure and robotic probes to deep space. The SLS Block 1 vehicle, the initial variant to fly, is optimized for lunar missions with a proven propulsion system consisting of four liquid hydrogen (LH2)/liquid oxygen (LOX)-fed RS-25 engines and twin five-segment solid rocket boosters (SRBs). The Block 1 vehicle can also be outfitted with an industry-standard 5 m-class payload fairing (the “cargo” configuration) and will launch at least 27 metric tons (t) of mass to trans-lunar injection (TLI). SLS is the backbone of NASA’s Artemis program, which will return the agency’s human spaceflight program to the Moon for the first time since 1972. For the Artemis I mission, SLS will send an uncrewed Orion spacecraft to TLI, where it will enter a distant retrograde lunar orbit and fly 38,000 nmi past the Moon – farther than any spacecraft built for humans has ever traveled. The SLS Block 1 vehicle for Artemis I completed manufacturing in 2019. Several elements, including the upper stage, have been delivered to the Exploration Ground Systems (EGS) program at KSC and are being prepped for integration and stack-ing. The five-segment solid rocket boosters – the largest and most powerful ever built for flight – are also complete. The booster motor segments for the Artemis I flight are scheduled to ship from prime contractor Northrop Grumman’s Utah facilities and begin stacking and integration at KSC in June 2020. The SLS core stage is the largest rocket stage NASA has ever built in terms of volume and height, and includes the avionics and the tanks that feed cryogenic propellant to the four RS-25s (formerly Space Shuttle Main Engines [SSMEs]). They have been modified with an updated controller and nozzle insulation to protect them from the hotter launch environment. The SLS core stage is currently being test-ed at NASA’s Stennis Space Center (SSC) in a series of “green run” tests to verify it meets design and performance requirements. Following the green run test series, which is scheduled to culminate with a full-duration hot-fire of the four RS-25 engines, the core stage will ship to KSC and be stacked between the sol-id rocket boosters in the Vehicle Assembly Building (VAB). Integration of the vehicle will continue with the upper stage, known as the Interim Cryogenic Propulsion Stage (ICPS) and the Launch Vehicle Stage Adapter (LVSA) on the core stage. Another adapter, the Orion Stage Adapter (OSA), connects SLS to Orion and provides housing for 13 6U CubeSat payloads manifested on Artemis I. The CubeSats will be released in deep space after Orion separates from the vehicle, and the flight marks the first ride share opportunity for independent small-sats to deep space. The second major SLS variant to come online, Block 1B, replaces the single-engine ICPS with a four-engine LH2/LOX Exploration Upper Stage (EUS). This more powerful upper stage, along with other vehicle up-grades, will enable the Block 1B vehicle to launch 38-42 t to TLI, depending on crew or cargo configuration. The final evolution of the vehicle, Block 2, will onramp evolved solid rocket boosters to increase mass to TLI to 43-46 t, de-pending on crew or cargo configuration. The Block 1B/Block 2 vehicles can be outfitted with an 8.4 m-diameter payload fairing in 19.1 m or 27.4 m lengths, to provide unprecedented volume for payloads. With the initial Block 1 vehicle completely manufactured and the core stage in final testing before shipping to KSC, the SLS Program and its industry partners have made significant progress manufacturing subsequent vehicles. For the second Block 1 vehicle, the solid rocket motor segments are complete, as are the RS-25 engines with controllers. All five major components of the Artemis II core stage – the forward skirt, LOX and LH2 tanks, intertank and engine section – are manufactured and technicians are installing subsystems at NASA’s rocket factory, Michoud Assembly Facility. The RL-10 engine for the Artemis II ICPS is complete and panels have been machined for its LH2 tank. In addition, panels are machined for the vehicle’s two adapters, with welding scheduled to begin in summer 2020. Flight hard-ware is also in production for the third SLS vehicle, with several booster motor segments cast. The pace of development on the EUS has increased, with the goal to complete Critical Design Review (CDR) in December 2020. Several EUS test rings have been machined at Michoud. The EUS is designed to exe-cute a variety of missions – human spaceflight, deployment of deep-space infra-structure, or high-C3 missions to the outer solar system – with crew and cargo configurations available beginning in the mid-2020s. The near-term goal for the nation’s powerful new space exploration asset, however, is to launch the Arte-mis program, and send the first woman and the next man to the lunar surface. At the Astrodynamics Specialist Conference, the SLS program will update the community on the progress of the initial Block 1 vehicle in final green run test-ing, integration and stacking. In addition, this paper will provide an update to the community on the manufacturing status of subsequent Block 1 and Block 1B vehicles.

Steve Creech↗

DOE OSTI2016

Origin and heterogeneity of pore sizes in the Mount Simon Sandstone and Eau Claire Formation: Implications for multiphase fluid flow

The Mount Simon Sandstone and Eau Claire Formation represent a principal reservoir - caprock system for wastewater disposal, geologic CO 2 storage, and compressed air energy storage (CAES) in the Midwestern United States. Of primary concern to site performance is heterogeneity in flow properties that could lead to non-ideal injectivity and distribution of injected fluids (e.g., poor sweep efficiency). Using core samples from the Dallas Center Structure, Iowa, we investigate pore structure that governs flow properties of major lithofacies of these formations. Methods include gas porosimetry and permeametry, mercury intrusion porosimetry, thin section petrography, and X-ray diffraction. The lithofacies exhibit highly variable intra- and inter-informational distributions of pore throat and body sizes. Based on pore-throat size, samples fall into four distinct groups. Micropore-throat dominated samples are from the Eau Claire Formation, whereas the macropore-, mesopore-, and uniform-dominated samples are from the Mount Simon Sandstone. Complex paragenesis governs the high degree of pore and pore-throat size heterogeneity, due to an interplay of precipitation, non-uniform compaction, and later dissolution of cements. Furthermore, the cement dissolution event probably accounts for much of the current porosity in the unit. The unusually heterogeneous nature of the pore networks in the Mount Simon Sandstone indicates that there is a greater-than-normal opportunity for reservoir capillary trapping of non-wetting fluids — as quantified by CO 2 and air column heights — which should be taken into account when assessing the potential of the reservoir-caprock system for CO 2 storage and CAES.

58 GEOSCIENCES↗

NASA NTRS1994

Simultaneous observations of VLF ground transmitter signals on the DE 1 and COSMOC 1809 satellites: Detection of a magnetospheric caustic and a duct

Khabarovsk transmitter signals (15.0 kHz, 48 deg N, 135 deg E) were observed on the high-altitude (approximately 15000 km) Dynamic Explorer 1 (DE 1) and the low-altitude (approximately 960) km COSMOS 1809 satellites during a 9-day period in August 1989. On 7 out of 9 days the linear wave receiver (LWR) on the DE 1 satellite also detected signals from the Alpha transmitter (11.9-15.6 kHz) in Russia and an Omega transmitter (10.2-13.6 kHz) in Australia, as well as natural VLF emissions such as hiss, chorus, whistlers, and wideband impulsive signals. On two days, August 23 and 27, 1989, observations of the Khabarovsk transmitter signals were simultaneously carried out at high altitude on thre DE 1 satellite and at low altitude on the COSMOS 1809 satellite. Analysis of data from these 2 days has led to several new results on the propagation of whistler mode signals in the Earth's magnetosphere. New evidence was found of previously reported propagation phenomena, such as (1) confinement of transmitter signals in the conjugate hemisphere at ionospheric heights (approximately 1000 km), (2) observation of direct multipath propagation on both DE 1 and COSMOS 1809, (3) detection of ionospheric irregularities of is less than or equal to 100 km scale size with a few percent enhancement in electron density. We also report the first direct detection of a magnetospheric dust at L = 2.94 which was believed to be responsible for the ducted propagation of Khabarovsk signals observed in the COSMOS 1809 satellite. Our study also indicates that duct end points can extend down to at least approximately 1000 km. The peak electric and magnetic fields of ducted Khabarovsk transmitter signals at approximately 1000 km were 520 micro V/m and 36 pT respectively. Estimated field strengths of these signals inside the dust at the geomagnetic equator were 57 micro V/m and 12 pT for electric and magnetic field respectively. The results of two-dimensional ray tracing simulations were consistent with the observations of the nonducted whistler-mode propagation of Khabarovsk (15 kHz) and Alpha (11.9 kHz) signals from the transmitter location to the DE 1 and COSMOS 1809 satellites. Our results have direct implications for the question of accessibility of waves injected from the ground to various regions of the ionosphere and the magnetosphere.

Chmyrev, V. M.↗

DOE OSTI2021

Characterization of oceanic post-cold frontal clouds and their model representation. Final report

Low-level clouds are ubiquitous over the oceans, and their impact on the Earth’ radiative budget is important. However, general circulation models (GCMs) that are used to represent our climate and its evolution still experience problems to represent their correct amount and radiative impact. In the extratropical latitudes (~30-60 N/S), cloud properties are strongly modulated by the occurrence of anti-cyclones (high pressure systems) and cyclones (low pressure systems). While deep clouds dominate in cyclones near the center and along the cold and warm fronts, low-level clouds populate the colder regions of the cyclones, and more specifically the region in the wake of the cold fronts, the post-cold frontal region (PCF). The GCMs mentioned above were found to underestimate the amount of these PCF clouds, and this caused errors in their representation of the amount of solar radiation at the surface of the southern hemisphere oceans, causing an excess in absorption and errors in long-term predictions. To better understand which processes are ill-represented in the GCMs to cause such issues, we used ground-based observations from the Eastern North Atlantic (ENA) ARM site and a high resolution regional model (called WRF) to examine the properties of PCF clouds. With the WRF model, we simulated the passage of a cold front at the ENA site using theoretically distinct representations of convection and the physics of the boundary layer. While these physical representations had little impact on the timing, structure and circulation of the cold front, we found a high sensitivity to how the large scale information was fed into the model. Errors rapidly developed if the incoming flow information was injected too far from the site, but 1000 km was found to be an optimal distance between the outer boundary and the ENA site for a realistic cold front passage. With this model, we then conducted a series of experiment where convection and boundary layer schemes were changed, and revealed that 1) the clouds were more sensitive to the convection than boundary layer schemes and 2) the interaction between the two has a large impact. Furthermore, changing the representation of the physics in the model had greater impacts on PCF clouds than changing the initial conditions. As large-scale models increase their spatial resolution, convection will eventually be resolved but it remains that their representation of the physics of the boundary layer will matter for the correct representation of shallow (low-level) clouds. Observations were used to find which large-scale information had a clear relationship with PCF cloud and precipitation properties, in view of helping with GCM parameterizations and evaluation. Both lower troposphere stability and surface forcing were found to be important for PCF clouds that experience more vigorous dynamics than their tropical or subtropical counterparts. For cloud base and top heights, a simple measure of the change in potential temperature between the surface and the 800 hPa was found to be a good predictor. This relationship was found to hold as well when using data from campaigns recently conducted in the southern ocean south of Tasmania. A recent GCM was found to succeed in representing similar relations. For cloud depth, water content, and different characteristics of precipitation both the temperature contrast between the surface and near-air and the near-surface wind speed had a strong relationship. To summarize, through observational constraints and regional model experiments, the project has provided guidance to evaluate GCMs and tools to improve their representation of low-level clouds in the extratropics. It has also highlighted differences in behavior between these clouds and those typically found in the subtropics, thus opening new avenues of research for shallow clouds in general.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

On hybrid scenarios in KSTAR

We report the status of hybrid scenario experiments in Korea Superconducting Tokamak Advanced Research (KSTAR). The hybrid scenario is defined as stationary discharges with ${\beta _{\rm{N}}}$ $ \geqslant $ 2.4 and ${H_{89}}$ $ \geqslant $ 2.0 at ${q_{95}}$ < 6.5 without or with very mild sawtooth activities in KSTAR. It is being developed towards reactor-relevant conditions. High performance of ${\beta _{\rm{N}}} \lesssim 3.0$, ${H_{89}} \lesssim 2.4$ and G-factor ($ \equiv {\beta _{\rm{N}}}{H_{89}}/q_{95}^2$) $ \lesssim 0.4$6 has been achieved and sustained for $ \gtrsim 40{{\tau }_{\text{E}}}$ at ne/nGW ~ 0.7 with heating power of $ \lesssim $5 MW. Some KSTAR hybrid discharges exhibit a unique feature of a slow transition from conventional H-mode to hybrid mode after the third neutral beam injection. The reason for the confinement enhancement is extensively studied in this transition period of a representative discharge exhibiting a common feature of KSTAR hybrid scenarios. 0D performance analysis with magnetohydrodynamic activities, 1D kinetic profile dynamics, power balance analysis, linear gyro-kinetic analysis and edge pedestal stability analysis were conducted. The enhancement is thought to be from both the core and the pedestal. The improvement in the core region of the ion energy channel is observed from the linear gyro-kinetic analysis considering the electromagnetic, the fast ion, the Shafranov shift, ${\it{\omega }_{E \times B}}$, and the magnetic shear effect. The electromagnetic finite $\beta $ stabilisation plays a role in the inner core region at ${\rho _{\rm{tor}}}\sim0.35$ together with the fast ion effect. The alpha stabilisation effect is also found at ${\rho _{\rm{tor}}}\sim0.5$. ${\omega _{E \times B}}$, which could reduce the linear growth of the ion temperature gradient mode in the outer core region at ${\rho _{\rm{tor}}}\sim0.5 - 0.7$ with the highest contribution from the toroidal rotation. Furthermore, regarding the improvement in the pedestal, Shafranov shift broadens the stability boundary of the pedestal in support of the diamagnetic effect. The pedestal height and width could be reproduced by the EPED model, while a realistic current profile is used to calculate the internal inductance for Shafranov shift. Based on these findings, a comprehensive confinement enhancement mechanism has been proposed by considering the core-edge interplay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1996

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗

DOE OSTI2014

Emergence of granular-sized magnetic bubbles through the solar atmosphere. I. Spectropolarimetric observations and simulations

We study a granular-sized magnetic flux emergence event that occurred in NOAA 11024 in 2009 July. The observations were made with the CRISP spectropolarimeter at the Swedish 1 m Solar Telescope achieving a spatial resolution of 0.''14. Simultaneous full Stokes observations of the two photospheric Fe I lines at 630.2 nm and the chromospheric Ca II 854.2 nm line allow us to describe in detail the emergence process across the solar atmosphere. We report here on three-dimensional (3D) semi-spherical bubble events, where instead of simple magnetic footpoints, we observe complex semi-circular feet straddling a few granules. Several phenomena occur simultaneously, namely, abnormal granulation, separation of opposite-polarity legs, and brightenings at chromospheric heights. However, the most characteristic signature in these events is the observation of a dark bubble in filtergrams taken in the wings of the Ca II 854.2 nm line. There is a clear coincidence between the emergence of horizontal magnetic field patches and the formation of the dark bubble. We can infer how the bubble rises through the solar atmosphere as we see it progressing from the wings to the core of Ca II 854.2 nm. In the photosphere, the magnetic bubble shows mean upward Doppler velocities of 2 km s{sup –1} and expands at a horizontal speed of 4 km s{sup –1}. In about 3.5 minutes it travels some 1100 km to reach the mid chromosphere, implying an average ascent speed of 5.2 km s{sup –1}. The maximum separation attained by the magnetic legs is 6.''6. From an inversion of the observed Stokes spectra with the SIR code, we find maximum photospheric field strengths of 480 G and inclinations of nearly 90° in the magnetic bubble interior, along with temperature deficits of up to 250 K at log τ = –2 and above. To aid the interpretation of the observations, we carry out 3D numerical simulations of the evolution of a horizontal, untwisted magnetic flux sheet injected in the convection zone, using the Bifrost code. The computational domain spans from the upper convection zone to the lower corona. In the modeled chromosphere, the rising flux sheet produces a large, cool, magnetized bubble. We compare this bubble with the observed ones and find excellent agreement, including similar field strengths and velocity signals in the photosphere and chromosphere, temperature deficits, ascent speeds, expansion velocities, and lifetimes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

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↗

NASA NTRS2022

Radiative Losses in the Chromosphere During A C-Class Flare

Context. Solar flares release an enormous amount of energy (~10(exp 32) erg) into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. Aims. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. Methods. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semiempirical model atmospheres, which were inferred from a multiline inversion of observed Stokes profiles obtained with the CRISP and CHROMIS instruments on the Swedish 1-m Solar Telescope. The radiative losses were computed taking into account the effect of partial redistribution and non-local thermodynamic equilibrium. To estimate the integrated radiative losses in the chromosphere, the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. We also compared our time series of radiative losses with those from the RADYN flare simulations. Results. We obtained a high spatial-resolution map of integrated radiative losses around the flare peak time. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H and K and Mgii h and k lines are strong and comparable to the Ca IR triplet (~32kW m(exp -2)). Since our flare is a relatively weak event, the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution (~7kW m(exp -2)) in the selected limits of the chromosphere. The temporal evolution of total integrated radiative losses exhibits sharply rising losses (0.4kW m(exp -2) (s(exp -1)) and a relatively slow decay (0.23kW m(exp -2) s(exp -1)). The maximum value of total radiative losses is reached around the flare peak time and can go up to 175kWm􀀀2 for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of 5 × 10(exp 10) erg s(exp -1) cm(exp -2).

Chromosphere↗

NASA NTRS2018

Radiative Heating of Large Meteoroids During Atmospheric Entry

A high-fidelity approach for simulating the aerothermodynamic environments of meteor entries was developed, which allows the commonly assumed heat transfer coefficient of 0.1 to be assessed. This model uses chemically reacting computational fluid dynamics (CFD), coupled with radiation transport and surface ablation. Coupled radiation accounts for the impact of radiation on the flowfield energy equations, while coupled ablation explicitly models the injection of ablation products within the flowfield and radiation simulations. For a meteoroid with a velocity of 20 km/s, coupled radiation is shown to reduce the stagnation point radiative heating by over 60%. The impact of coupled ablation (with coupled radiation) is shown to provide at least a 70% reduction in the radiative heating relative to cases with only coupled radiation. This large reduction is partially the result of the low ionization energies of meteoric ablation products relative to air species. The low ionization energies of ablation products, such as Mg and Ca, provide strong photoionization and atomic line absorption in regions of the spectrum that air species do not. MgO and CaO are also shown to provide significant absorption. Turbulence is shown to impact the distribution of ablation products through the shock-layer, which results in up to a 100% increase in the radiative heating downstream of the stagnation point. To create a database of heat transfer coefficients, the developed model was applied to a range of cases. This database considered velocities ranging from 14 to 20 km/s, altitudes ranging from 20 to 50 km, and nose radii ranging from 1 to 100 m. The heat transfer coefficients from these simulations are below 0.045 for the range of cases, for both laminar and turbulent, which is significantly lower than the canonical value of 0:1. When the new heat transfer model is applied to a Tunguska-like 15 Mt entry, the effect of the new model is to lower the height of burst by up to 2 km, depending on assumed entry angle. This, in turn, results in a significantly larger ground damage footprint than when the canonical heating assumption is used.

Christopher O Johnston↗

NASA NTRS1990

Near-infrared mapping of spiral barred galaxies

In external galaxies, near-infrared emission originates from stellar populations, hot dust, free-free emission from H+ regions, gaseous emission, non-thermal nucleus if any. Because of the low extinction compared to the visible, infrared wavelengths are useful to probe regions obscured by dust such as central parts where starburst phenomena can occur because of the large quantity of matter. The results presented were obtained with a 32 x 32 InSb charge injection device (CID) array cooled at 4K, at the f/36 cassegrain focus of the 3m60 Canada-France-Hawaii telescope with a spatial resolution of 0.5 inches per pixel. The objects presented are spiral barred galaxies mapped at J(1.25 microns), H(1.65 microns) and K(2.2 microns). The non-axisymetric potential due to the presence of a bar induces dynamical processes leading to the confinement of matter and peculiar morphologies. Infrared imaging is used to study the link between various components. Correlations with other wavelengths ranges and 2-colors diagrams ((J-H), (H-K)) lead to the identification of star forming regions, nucleus. Maps show structures connected to the central core. The question is, are they flowing away or toward the nucleus. Observations of M83 lead to several conclusions. The star forming region, detected in the visible and the infrared cannot be very compact and must extend to the edge of the matter concentration. The general shape of the near-infrared emission and the location of radio and 10 micron peaks suggest the confinement of matter between the inner Linblad resonances localized from CO measurements about 100 and 400 pc. The distribution of color indices in the arc from southern part to the star forming region suggests an increasing amount of gas and a time evolution eventually triggered by supernova explosions. Close to the direction of the bar, a bridge-like structure connects the arc to the nucleus with peculiar color indices. Perhaps, this structure can be linked to a height velocity component seen in UV and we can attribute it to a jet and/or a matter flow along the bar toward the nucleus, fuelling it. NGC 1068 is the nearest Seyfert 2 galaxy. It has been a subject of many studies at all wavelengths. This object was mapped at J, H, K, L and M, and in polaro-imagery. Results are given.

Gallais, P.↗

DOE OSTI2019

MINERALOGICAL AND TEXTURAL CONTROLS ON SHEAR STRENGTH, SLIP STABILITY AND PERMEABILITY OF FAULTS

Induced seismicity resulting from fluid injection into the subsurface related to water and CO2 disposal, hydraulic fracturing and the stimulation of geothermal reservoirs present an important societal concern. These human activities involve the injection of large volumes of pressurized fluid into the subsurface, potentially at high rates, raising local pore pressures and disturbing the pristine local stress regime by lowering effective normal stress on pre-existing faults and fractures. The reduction of effective normal stress may trigger fault/fracture reactivation and in some cases result in hazardous seismic ruptures. Effective management and engineering of anthropogenic seismic events requires substantial understanding in the mechanisms, especially for the controlling factors on coupled rheological and transport response, including fault shear strength, slip stability, and permeability evolution during such events. In this study, we explore the coupled rheological and transport response of faults and fractures during reactivation as controlled by two fundamental controlling properties, viz., mineralogy and textural features. We approach this problem through shear experiments on analog faults and fractures via laboratory and numerical experiments. Specifically, we investigate: (1) the influence of frictionally weak minerals (talc) in mixtures of mineral analogs featuring contrasting frictional properties, (2) the influence of iron oxide grain coatings on quartz aggregates, and (3) the influence of fracture roughness in mated fractures on the ensemble shear strength, slip stability, and permeability evolution during reactivation events. We address the following questions in this study: (1) how much and what distribution of frictionally weak minerals is required to induce significant weakening in faults consisting of a matrix of frictionally strong minerals, (2) how does a pre-imposed weak mineral layer influence the rheological and transport behavior of faults, (3) what is the influence of a trace amount of grain coating materials introduce on the coupled behavior of faults, and finally (4) how do asperity height and wave length control the ensemble behavior of faults. These questions are explicitly answered in the following.

58 GEOSCIENCES↗

DOE OSTI2018

High-resolution time-resolved experiments on mixing and entrainment of buoyant jets in stratified environments

Fluid jets interacting with a stratified layer play an important role in the safety of several reactor designs. In the containment of nuclear power plants, fluid jets dominate the transport and mixing of gaseous species and consequent hydrogen distribution in case of a severe accident. The mixing phenomena in the containment are driven by buoyant high-momentum injections (jets) and low momentum injection plumes. Mixing near the postulated break is initially dominated by high flow velocities. Plumes with moderate flow velocities are instead relevant in the break compartment during the long-term pressurization phase, or in any of the apertures between two connected compartments if the mass flows are sufficiently high and the density differences between efflux and ambient are sufficiently low. Phenomena of interest include free plumes (as produced by the efflux from the break compartment in a larger room or directly from a break flow), wall plumes (such those produced by low mass flows through inter-compartment apertures), and propagating stratification fronts in the ambient (for any stably stratified conditions). These phenomena have been highly ranked about nuclear reactor design, especially regarding of safety protocols. During a Pressurized Thermal Shock (PTS) scenario, the interaction between the cold ECCS injection plume and the stratified fluid present in the cold (or hot) leg is important in order to determine the temperature at the time-dependent temperature at the inlet of the reactor pressure vessel (RPV) and the potential to cause a thermal shock on the RPV wall. In sodium-cooled fast reactors (SFRs), core channels are typically hydro-dynamically isolated so that there exists a considerable temperature variation at the exit of adjacent fuel assemblies. All the above phenomena are characterized by the interaction of buoyant jets with the stratified flow. In stratified layers baroclinic forces create significant redistribution of turbulent kinetic energy and scales, which leads to anisotropy. This important physical phenomenon is highly three dimensional and is challenging to capture even with high-fidelity CFD simulations, due in part to lack of sufficiently resolved validation data. Furthermore, the experimental data available in the open literature do not feature the level of fidelity needed for an extensive validation of turbulence models in lower order CFD. To shed new lights into the crucial phenomena object of the present research project, it was proposed to conduct coordinated experiments and simulations at the University of Michigan and the George Washington University. The project has resulted in an experimental database of high-resolution time-resolved measurements of jets in uniform and stratified environments. The novel experimental data will be used to validate computational fluid dynamic (CFD) codes, including both Large Eddy Simulations (LES) and unsteady Reynolds-averaged Navier-Stokes equations (URANS) methodologies. In the Experimental and Multiphase flow (ECMF) laboratory at Univerisity of Michigan, we built two experimental facilities to investigate also the effect of scaling. The first facility, DESTROJER (DEnsity Stratified Turbulent ROund free Jet ExpeRiment), featuring a contoured jet nozzle with a diameter of D=12.7mm and a 1m×1m×1m cubic tank, which is made of acrylic glass for optical access. The ratio between the tank width and the nozzle diameter is equal to 78, which ensures that there is no direct interaction between the jet and the side walls. A second, modular experimental facility, features three different tank sizes of size 10×10×30, 20×20×30, 30×30×30 cm3 respectively (all tanks have the same height of 30 cm), and a jet diameter of 2mm. For the modular facility, tank-to-nozzle ratio of 50, 100, and 150 are obtained respectively. Experiments with different-density jet impactions and sharp interface with a density difference of 3.16% In the Laboratory at George Washington University, complementary experiments have been carried out at a facility featuring the same dimensions as DESTROYER. Experiments with different jet diameters have been carried out as well. High-resolution time-resolved measurements have been performed at all facilities using high-speed particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF). In this way, simultaneous measurements of velocity and density fields have been carried out. At the University of Michigan, wire-mesh sensors have been designed and built to gather high-resolution, time-dependent data of the stratified front in the tank. These measurements are important to provide accurate initial conditions for the CFD simulations. An important outcome of the project has been the development of novel techniques to achieve refractive index matching (RIM) of acqueous solutions with high density differences. RIM is crucial if optical techniques such as PIV are used to measure velocity fields. Previous to this project, experiments up to 3% density difference had been reported in the literature. With the methodology developed in this project, we have been able to push the envelop to density differences up to 9%. The experiments have been used to validate Reynolds Averaged Navier-Stokes (RANS) turbulence models and Large Eddy Simulations (LES) models. Both NEK5000 and the commercial CFD code STAR-CCM+ have been used. The experimental data will be made available to the community through the NE-KAMS databasehave been successfully conducted in the scaled facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗