SEARCH · Engineering Papers

Results for “Gas property variation”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1999

Instrumentation Working Group Summary

The Instrumentation Working Group compiled a summary of measurement techniques applicable to gas turbine engine aerosol precursors and particulates. An assessment was made of the limits, accuracy, applicability, and technology readiness of the various techniques. Despite advances made in emissions characterization of aircraft engines, uncertainties still exist in the mechanisms by which aerosols and particulates are produced in the near-field engine exhaust. To adequately assess current understanding of the formation of sulfuric acid aerosols in the exhaust plumes of gas turbine engines, measurements are required to determine the degree and importance of sulfur oxidation in the turbine and at the engine exit. Ideally, concentrations of all sulfur species would be acquired, with emphasis on SO2 and SO3. Numerous options exist for extractive and non-extractive measurement of SO2 at the engine exit, most of which are well developed. SO2 measurements should be performed first to place an upper bound on the percentage of SO2 oxidation. If extractive and non-extractive techniques indicate that a large amount of the fuel sulfur is not detected as SO2, then efforts are needed to improve techniques for SO3 measurements. Additional work will be required to account for the fuel sulfur in the engine exhaust. Chemical Ionization Mass Spectrometry (CI-MS) measurements need to be pursued, although a careful assessment needs to be made of the sampling line impact on the extracted sample composition. Efforts should also be placed on implementing non-intrusive techniques and extending their capabilities by maximizing exhaust coverage for line-of-sight measurements, as well as development of 2-D techniques, where feasible. Recommendations were made to continue engine exit and combustor measurements of particulates. Particulate measurements should include particle size distribution, mass fraction, hydration properties, and volatile fraction. However, methods to ensure that unaltered samples are obtained need to be developed. Particulate speciation was also assigned a high priority for quantifying the fractions of carbon soot, PAH, refractory materials, metals, sulfates, and nitrates. High priority was also placed on performing a comparison of particle sizing instruments. Concern was expressed by the workshop attendees who routinely make particulate measurements about the variation in number density measured during in-flight tests by different instruments. In some cases, measurements performed by different groups of researchers during the same flight tests showed an order of magnitude variation. Second priority was assigned to measuring concentrations of odd hydrogen and oxidizing species. Since OH, HO2, H2O2, and O are extremely reactive, non-extractive measurements are recommended. A combination of absorption and fluorescence is anticipated to be effective for OH measurements in the combustor and at the engine exit. Extractive measurements of HO2 have been made in the stratosphere, where the ambient level of OH is relatively low. Use of techniques that convert HO2 to OH for combustor and engine exit measurements needs to be evaluated, since the ratio of HO2/OH may be 1% or less at both the combustor and engine exit. CI-MS might be a viable option for H2O2, subject to sampling line conversion issues. However, H2O2 is a low priority oxidizing species in the combustor and at the engine exit. Two candidates for atomic oxygen measurements are Resonance Enhanced Multi-Photon Ionization (REMPI) and Laser-Induced Fluorescence (LIF). Particulate measurement by simultaneous extractive and non-extractive techniques was given equal priority to the oxidizer measurements. Concern was expressed over the ability of typical ground test sampling lines to deliver an unaltered sample to a remotely located instrument. It was suggested that the sampling probe and line losses be checked out by attempting measurements using an optical or non-extractive technique immediately upstream of the sampling probe. This is a possible application for Laser Induced Incandescence (LII) as a check on the volume fraction of soot. Optical measurements of size distribution are not well developed for ultrafine particles less than about 20 nm in diameter, so a non-extractive technique for particulate size distribution cannot be recommended without further development. Carbon dioxide measurements need to be made to complement other extractive measurement techniques. CO2 measurements enable conversion of other species concentrations to emission indices. Carbon monoxide, which acts as a sink for oxidizing species, should be measured using non-extractive techniques. CO can be rapidly converted to CO2 in extractive probes, and a comparison between extractive and non-extractive measurements should be performed. Development of non-extractive techniques would help to assess the degree of CO conversion, and might be needed to improve the concentration measurement accuracy. Measurements of NO(x) will continue to be critical due to the role of NO and NO2 in atmospheric chemistry, and their influence on atmospheric ozone. Time-resolved measurements of temperature, velocity, and species concentrations were included on the list of desired measurement. Thermocouples are typically adequate for engine exit measurements. PIV and LDV are well established for obtaining velocity profiles. The techniques are listed in the accompanying table; are divided into extractive and non-extractive techniques. Efforts were made to include a measurement uncertainty for each technique. An assessment of the technology readiness was included.

Zaller, Michelle↗

DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2009

Inert Welding/Brazing Gas Filters and Dryers

The use of hybridized carbon/silicon carbide (C/SiC) fabric to reinforce ceramic matrix composite face sheets and the integration of such face sheets with a foam core creates a sandwich structure capable of withstanding high-heat-flux environments (150 W/sq cm) in which the core provides a temperature drop of 1,000 C between the surface and the back face without cracking or delamination of the structure. The composite face sheet exhibits a bilinear response, which results from the SiC matrix not being cracked on fabrication. In addition, the structure exhibits damage tolerance under impact with projectiles, showing no penetration to the back face sheet. These attributes make the composite ideal for leading-edge structures and control surfaces in aerospace vehicles, as well as for acreage thermal protection systems and in high-temperature, lightweight stiffened structures. By tailoring the coefficient of thermal expansion (CTE) of a carbon fiber containing ceramic matrix composite (CMC) face sheet to match that of a ceramic foam core, the face sheet and the core can be integrally fabricated without any delamination. Carbon and SiC are woven together in the reinforcing fabric. Integral densification of the CMC and the foam core is accomplished with chemical vapor deposition, eliminating the need for bond-line adhesive. This means there is no need to separately fabricate the core and the face sheet, or to bond the two elements together, risking edge delamination during use. Fibers of two or more types are woven together on a loom. The carbon and ceramic fibers are pulled into the same "pick" location during the weaving process. Tow spacing may be varied to accommodate the increased volume of the combined fiber tows while maintaining a target fiber volume fraction in the composite. Foam pore size, strut thickness, and ratio of face sheet to core thickness can be used to tailor thermal and mechanical properties. The anticipated CTE for the hybridized composite is managed by the choice of constituents, varying fiber tow sizes and constituent part ratios. This structural concept provides high strength and stiffness at low density 1.06 g/cu cm in panels tested. Varieties of face sheet constructions are possible, including variations in fiber type and weave geometry. The integrated structures possible with this composite could eliminate the need for non-load-bearing thermal protection systems on top of a structural component. The back sheet can readily be integrated to substructures through the incorporation of ribs. This would eliminate weight and cost for aerospace missions.

Goudy, Jerry↗

NASA NTRS1997

Description and Analysis of Core Samples: The Lunar Experience

Although no samples yet have been returned from a comet, extensive experience from sampling another solar system body, the Moon, does exist. While, in overall structure, composition, and physical properties the Moon bears little resemblance to what is expected for a comet, sampling the Moon has provided some basic lessons in how to do things which may be equally applicable to cometary samples. In particular, an extensive series of core samples has been taken on the Moon, and coring is the best way to sample a comet in three dimensions. Data from cores taken at 24 Apollo collection stations and 3 Luna sites have been used to provide insight into the evolution of the lunar regolith. It is now well understood that this regolith is very complex and reflects gardening (stirring of grains by micrometeorites), erosion (from impacts and solar wind sputtering), maturation (exposure on the bare lunar surface to solar winds ions and micrometeorite impacts) and comminution of coarse grains into finer grains, blanket deposition of coarse-grained layers, and other processes. All of these processes have been documented in cores. While a cometary regolith should not be expected to parallel in detail the lunar regolith, it is possible that the upper part of a cometary regolith may include textural, mineralogical, and chemical features which reflect the original accretion of the comet, including a form of gardening. Differences in relative velocities and gravitational attraction no doubt made this accretionary gardening qualitatively much different than the lunar version. Furthermore, at least some comets, depending on their orbits, have been subjected to impacts of the uppermost surface by small projectiles at some time in their history. Consequently, a more recent post-accretional gardening may have occurred. Finally, for comets which approach the sun, large scale erosion may have occurred driven by gas loss. The uppermost material of these comets may reflect some of the features of this erosional process, such as crust formation, and variations with depth might be expected. Overall, the upper few meters of a comet may be as complex in their own way as the upper few meters of the lunar regolith have proven to be, and by analogy, detailed studies of core samples containing this depth information will be needed to understand these processes and the details of the accretional history and the subsequent alteration history of comets.

McKay, David S.↗

NASA NTRS2001

Diffusivity Measurements Made Instant and Easy

A compact common path interferometer (CPI) system has been developed to measure the diffusivity of liquid pairs. The CPI is an optical technique that can be used to measure changes in the gradient of the refraction index of transparent materials. It uses a shearing interferometer that shares the same optical path from a laser light source to the final imaging plane. The molecular diffusion coefficient of liquids can be determined from the physical relations between changes in the optical path length and liquid phase properties. When the data obtained by using the CPI have been compared with similar results from other techniques, the instrument has been demonstrated to be far superior to other instruments for measuring the diffusivity of miscible liquids while staying very compact and robust (ref. 1). Because of its compactness and ease of use, the CPI has been adopted for use in studies of interface dynamics as well as other diffusion-controlled process applications (ref. 2). This progress will permit experiments in microgravity that can quantitatively answer basic science questions about mass and thermal diffusion and their effect in transport processes. This instrument is a spinoff of a diagnostic development for microgravity fluid physics experiments at the NASA Glenn Research Center that has used optics and electronics existing in the fluid physics laboratory for feasibility studies. Optical diagnostic techniques have become an integral part of many areas of measurement applications in industrial and research laboratories. Many types of interferometers and their phase-shifted versions have been used as instruments for measuring optical wave fronts for lens testing and combustion and fluid flow diagnostics. One of these, the point diffraction interferometer, is considered to be robust (see, for example, ref. 3) because it has a common-path design. The point diffraction interferometer is difficult to align and has a limited measurement range for liquid-phase applications. Interferometry and schlieren techniques have been widely used for many years for gas-flow-visualization. On the other hand, the CPI, which is basically a Wollaston prism Polarizer in combination with an analyser (crystalline quartz retardation plate) is used instead of a point diffraction interferometer unit. The advantage of using the CPI over other optical techniques is that it can make quantitative measurements in liquids with the large index of refraction variations that often occur in interface dynamics studies. This can be simply accomplished by using different prism polarizers for each particular experimental condition.

Rashidnia, Nasser↗

NASA NTRS2013

Hydrometeor Size Distribution Measurements by Imaging the Attenuation of a Laser Spot

The optical extinction of a laser due to scattering of particles is a well-known phenomenon. In a laboratory environment, this physical principle is known as the Beer-Lambert law, and is often used to measure the concentration of scattering particles in a fluid or gas. This method has been experimentally shown to be a usable means to measure the dust density from a rocket plume interaction with the lunar surface. Using the same principles and experimental arrangement, this technique can be applied to hydrometeor size distributions, and for launch-pad operations, specifically as a passive hail detection and measurement system. Calibration of a hail monitoring system is a difficult process. In the past, it has required comparison to another means of measuring hydrometeor size and density. Using a technique recently developed for estimating the density of surface dust dispersed during a rocket landing, measuring the extinction of a laser passing through hail (or dust in the rocket case) yields an estimate of the second moment of the particle cloud, and hydrometeor size distribution in the terrestrial meteorological case. With the exception of disdrometers, instruments that measure rain and hail fall make indirect measurements of the drop-size distribution. Instruments that scatter microwaves off of hydrometeors, such as the WSR-88D (Weather Surveillance Radar 88 Doppler), vertical wind profilers, and microwave disdrometers, measure the sixth moment of the drop size distribution (DSD). By projecting a laser onto a target, changes in brightness of the laser spot against the target background during rain and hail yield a measurement of the DSD's second moment by way of the Beer-Lambert law. In order to detect the laser attenuation within the 8-bit resolution of most camera image arrays, a minimum path length is required. Depending on the intensity of the hail fall rate for moderate to heavy rainfall, a laser path length of 100 m is sufficient to measure variations in optical extinction using a digital camera. For hail fall only, the laser path may be shorter because of greater scattering due to the properties of hailstones versus raindrops. A photodetector may replace the camera in automated installations. Laser-based rain and hail measurement systems are available, but they are based on measuring the interruption of a thin laser beam, thus counting individual hydrometeors. These systems are true disdrometers since they also measure size and velocity. The method reported here is a simple method, requiring far less processing, but it is not a disdrometer.

Lane, John↗

DOE OSTI2020

BWR Spent Nuclear Fuel Acquisition and Testing to Support DOE-NE High Burnup Spent Fuel Data Project

The Office of Spent Fuel and Waste Disposition (SFWD) within the US Department of Energy (DOE) Office of Nuclear Energy (NE) established the Spent Fuel and Waste Science and Technology (SFWST) campaign to conduct research and development (R&D) activities related to the storage, transportation, and disposal of spent nuclear fuel (SNF) and high-level radioactive waste. The SFWST program was created within SFWD to address issues of extended or long-term SNF storage and transportation. Some near-term objectives of SFWST are to use a science-based, engineering-driven approach to: Support the enhancement of the technical bases to support the continued safe and secure dry storage of SNF for extended periods; Support the enhancement of the technical bases for retrieving SNF after extended dry storage; Support the enhancement of the technical bases for transporting high burnup (HBU) fuel and transporting low burnup fuel and HBU fuel after dry storage DOE-NE, in partnership with the Electric Power Research Institute, developed the High Burnup Spent Fuel Data Project to perform a large-scale demonstration and laboratory-scale testing of HBU pressurized water reactor (PWR) fuels (exceeding 45 gigawatt-days per metric ton of uranium [GWd/MTU]). Under this project, 25 sister rods—which are rods that have the same design, power histories, and other characteristics—were removed from assemblies at the North Anna Nuclear Power Station and sent to Oak Ridge National Laboratory (ORNL) in January 2016. ORNL performed detailed nondestructive examination (NDE) on all 25 rods. The NDE consisted of visual examinations, gamma and neutron scanning, profilometry and rod length measurements, and eddy current examinations. After completing the NDE, 10 of the sister rods were delivered to Pacific Northwest National Laboratory (PNNL) in a NAC International, Inc. legal-weight truck cask in September 2018 for destructive examination (DE). To date, SFWD work has focused on the PWR fuel that is part of the Sister Rod Test program. No boiling water reactor (BWR) fuel has been tested in the program, and the data needs that were identified for the PWR fuel have not been collected for BWR fuel. The goal to obtain six to nine BWR rods and test them at ORNL will support closing this important data gap. BWR fuel comprises approximately 56% of the total fuel assemblies currently in storage at nuclear power plants in the United States. BWR nuclear fuel and cladding designs and manufacturing are significantly different from PWRs. Differences include the following: BWR fuel pellets are larger than PWR pellets; Variations of Zircaloy-2 (including liners) are used instead of the Zircaloy-4 cladding materials used in PWRs; Clad manufacturing and stress-relief processes are different between PWRs and BWRs; The fuel rod dimensions are different because larger rod diameters and thicker cladding are used in BWRs; BWR fuel typically has lower internal rod pressures and sees vastly different operating conditions than PWR fuel (i.e., two-phase flow); BWR assemblies are “canned,” meaning each assembly is surrounded by a metal fuel channel; BWR cladding is often composed of an inner pure Zr liner that has widely different mechanical properties than the Zircaloy-2 alloy and exhibits a stronger affinity for hydrogen; The BWR SNF generally has more total hydrogen in the cladding/liner than typical PWR fuel; The construction of the PWR and BWR assemblies is vastly different; BWR rods are solidly attached to the assembly nozzles and experience a much different vibration and shock load than PWR rods, which are “floating” within a grid system attached to guide tubes, and the rods sit loosely on the bottom end plates. These numerous differences will affect the way the BWR SNF responds under dry storage preparation processes (e.g., vacuum drying) and during transportation. The results collected in the PWR experimental program must be compared with a subset of similar data collected on BWR SNF to establish a technical basis for whether the larger PWR database is sufficient to bound the BWR SNF end-of-life conditions as is currently assumed for several fuel/clad properties. Changes that occur in both fuel types at HBU could exacerbate any mechanical property differences. As the fuel burnup increases, several changes occur that might affect the performance of the fuel, cladding, and assembly hardware in storage and transportation. These changes include increased cladding corrosion layer thickness, increased cladding hydrogen content, increased cladding creep strains, increased fission gas release, and the formation of the HBU structure at the surface of the fuel pellets. The Nuclear Regulatory Commission (NRC) limits the current maximum rod-averaged burnup to 62 GWd/MTU due to these changes and the lack of data at higher burnups.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

CO 2 –Brine Contact Angle Measurements on Navajo, Nugget, Bentheimer, Bandera Brown, Berea, and Mt. Simon Sandstones

In this study, contact angles were measured for CO 2 bubbles on six different sandstones (Navajo, Nugget, Bentheimer, Bandera Brown, Berea, and Mt. Simon) that could potentially represent properties of CO 2 storage depositional environments. The impacts of pressure and temperature were studied by focusing on the CO 2 phase behavior in three different scenarios: gaseous, liquid, and supercritical conditions. Despite controlling the sample preparation and cleanliness, CO 2 -brine equilibration conditions, and pressure and temperature, there were inconsistencies in contact angle trends that could largely be attributed to natural sample heterogeneity resulting from localized variations in topography, surface roughness, and mineral composition across the surface. Despite these variations, the analysis of 298 measurements from this study showed that 81% of the contact angles were below 40° representing strongly water-wet to (moderately) water-wet behaviors. Also, 17.3% of the measurements were between 40° and 60° (weakly water-wet) and primarily belonged to small CO 2 bubbles (< 500 μm) that were heavily impacted by localized surface heterogeneity in natural sandstone samples. Additionally, only 1.7% of the measurements had angles greater than 60° and could be classified as CO 2 -wet behavior, but these measurements occurred on extremely small bubbles (100-200 μm) that were very dependent on localized surface heterogeneity. While some of the detailed physics of the CO 2 -brine-sandstone interface is still poorly understood, from an application standpoint, the sandstones of this study were best characterized as strongly water-wet to (moderately) water-wet.

03 NATURAL GAS↗

DOE OSTI2022

Representations and strategies for transferable machine learning improve model performance in chemical discovery

Strategies for machine-learning (ML)-accelerated discovery that are general across material composition spaces are essential, but demonstrations of ML have been primarily limited to narrow composition variations. By addressing the scarcity of data in promising regions of chemical space for challenging targets such as open-shell transition-metal complexes, general representations and transferable ML models that leverage known relationships in existing data will accelerate discovery. Over a large set (~1000) of isovalent transition-metal complexes, we quantify evident relationships for different properties (i.e., spin-splitting and ligand dissociation) between rows of the Periodic Table (i.e., 3d/4d metals and 2p/3p ligands). We demonstrate an extension to the graph-based revised autocorrelation (RAC) representation (i.e., eRAC) that incorporates the group number alongside the nuclear charge heuristic that otherwise overestimates dissimilarity of isovalent complexes. To address the common challenge of discovery in a new space where data are limited, we introduce a transfer learning approach in which we seed models trained on a large amount of data from one row of the Periodic Table with a small number of data points from the additional row. We demonstrate the synergistic value of the eRACs alongside this transfer learning strategy to consistently improve model performance. Analysis of these models highlights how the approach succeeds by reordering the distances between complexes to be more consistent with the Periodic Table, a property we expect to be broadly useful for other material domains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2008

Reliability and Confidence Interval Analysis of a CMC Turbine Stator Vane

High temperature ceramic matrix composites (CMC) are being explored as viable candidate materials for hot section gas turbine components. These advanced composites can potentially lead to reduced weight, enable higher operating temperatures requiring less cooling and thus leading to increased engine efficiencies. However, these materials are brittle and show degradation with time at high operating temperatures due to creep as well as cyclic mechanical and thermal loads. In addition, these materials are heterogeneous in their make-up and various factors affect their properties in a specific design environment. Most of these advanced composites involve two- and three-dimensional fiber architectures and require a complex multi-step high temperature processing. Since there are uncertainties associated with each of these in addition to the variability in the constituent material properties, the observed behavior of composite materials exhibits scatter. Traditional material failure analyses employing a deterministic approach, where failure is assumed to occur when some allowable stress level or equivalent stress is exceeded, are not adequate for brittle material component design. Such phenomenological failure theories are reasonably successful when applied to ductile materials such as metals. Analysis of failure in structural components is governed by the observed scatter in strength, stiffness and loading conditions. In such situations, statistical design approaches must be used. Accounting for these phenomena requires a change in philosophy on the design engineer s part that leads to a reduced focus on the use of safety factors in favor of reliability analyses. The reliability approach demands that the design engineer must tolerate a finite risk of unacceptable performance. This risk of unacceptable performance is identified as a component's probability of failure (or alternatively, component reliability). The primary concern of the engineer is minimizing this risk in an economical manner. The methods to accurately determine the service life of an engine component with associated variability have become increasingly difficult. This results, in part, from the complex missions which are now routinely considered during the design process. These missions include large variations of multi-axial stresses and temperatures experienced by critical engine parts. There is a need for a convenient design tool that can accommodate various loading conditions induced by engine operating environments, and material data with their associated uncertainties to estimate the minimum predicted life of a structural component. A probabilistic composite micromechanics technique in combination with woven composite micromechanics, structural analysis and Fast Probability Integration (FPI) techniques has been used to evaluate the maximum stress and its probabilistic distribution in a CMC turbine stator vane. Furthermore, input variables causing scatter are identified and ranked based upon their sensitivity magnitude. Since the measured data for the ceramic matrix composite properties is very limited, obtaining a probabilistic distribution with their corresponding parameters is difficult. In case of limited data, confidence bounds are essential to quantify the uncertainty associated with the distribution. Usually 90 and 95% confidence intervals are computed for material properties. Failure properties are then computed with the confidence bounds. Best estimates and the confidence bounds on the best estimate of the cumulative probability function for R-S (strength - stress) are plotted. The methodologies and the results from these analyses will be discussed in the presentation.

Murthy, Pappu L. N.↗

DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗

DOE OSTI2023

Electrostatic modulation of the lateral carrier density profile in field effect devices with nonlinear dielectrics

The properties of two-dimensional (2D) electronic systems are often effectively controlled using electrostatic gating. The geometry of such field effect devices influences the effectiveness of the gate and the carrier density profile in the 2D device. Here, in this work, we analyze the gate-induced spatial variations in the lateral carrier density in patterned LaAlO 3 /SrTiO 3 devices. We model the electrostatics of the 2D interface using the Thomas-Fermi approximation and compute the gate-induced charge distribution at the interface. We show that the electric field lines generated by the gate are focused at the edges of the device, causing an increased depletion near its edges. This effect is accentuated in LaAlO 3 /SrTiO 3 due to the large, nonlinear dielectric constant of the substrate, and the large distance between the gate electrode and the interface. We experimentally demonstrate one consequence of this effect by directly imaging current distributions in gated heterostructures, finding that insulating regions nucleate at the edges of the device due to the gate. Our results suggest that device geometry and choice of dielectric materials control the charge distribution in 2D systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2023

Three-Wavelength Approach for Aerosol-Cloud Discrimination in the Sage III/ISS Aerosol Extinction Dataset

The tropical upper troposphere and lower stratosphere (UTLS) region is dominated by aerosols and clouds affecting Earth’s radiation budget and climate. Thus, satellites’ continuous monitoring and identification of these layers is crucial for quantifying their radiative impact. However, distinguishing between aerosols and clouds is challenging, especially under the perturbed UTLS conditions during post-volcanic eruptions and wildfire events. Aerosol-cloud discrimination is primarily based on their disparate wavelength-dependent scattering and absorption properties. In this study, we use aerosol extinction observations in the tropical (15°N-15°S) UTLS from June 2017 to February 2021, available from the latest generation of the Stratospheric Aerosol and Gas Experiment (SAGE) instrument-SAGE III onboard the International Space Station (ISS) to study aerosols and clouds. During this period, the SAGE III/ISS provided better coverage over the tropics at additional wavelength channels (relative to previous SAGE missions) and witnessed several volcanic and wildfire events that perturbed the tropical UTLS. We explore the advantage of having an extinction coefficient at an additional wavelength channel (1550 nm) from the SAGE III/ISS in aerosol-cloud discrimination using a method based on thresholds of two extinction coefficient ratios, R1R1 (520 nm/1020 nm) and R2R2 (1020 nm/1550 nm). This method was proposed earlier by Kent et al. [Appl. Opt. 36, 8639 (1997) [CrossRef] ] for the SAGE III-Meteor-3M but was never tested for the tropical region under volcanically perturbed conditions. We call this method the Extinction Color Ratio (ECR) method. The ECR method is applied to the SAGE III/ISS aerosol extinction data to obtain cloud-filtered aerosol extinction coefficients, cloud-top altitude, and seasonal cloud occurrence frequency during the entire study period. Cloud-filtered aerosol extinction coefficient obtained using the ECR method revealed the presence of enhanced aerosols in the UTLS following volcanic eruptions and wildfire events consistent with the Ozone Mapping and Profiler Suite (OMPS) and space-borne lidar-Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). The cloud-top altitude obtained from the SAGE III/ISS is within 1 km of the nearly co-located observations from OMPS and CALIOP. In general, the seasonal mean cloud-top altitude from the SAGE III/ISS events peaks during the December, January, and February months, with sunset events showing higher cloud tops than the sunrise events, indicating the seasonal and diurnal variation of the tropical convection. The seasonal altitude distribution of cloud occurrence frequency obtained from the SAGE III/ISS also agrees well with CALIOP observations within 10%. We show that the ECR method is a simple approach that relies on thresholds independent of the sampling period, providing cloud-filtered aerosol extinction coefficients uniformly for climate studies irrespective of the UTLS conditions. However, since the predecessor of SAGE III did not include a 1550 nm channel, the usefulness of this approach is limited to short-term climate studies after 2017.

Clouds Aerosols↗

NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗

NASA NTRS1995

Extinction coefficient (1 micrometer) properties of high-altitude clouds from solar occultation measurements (1985-1990): Evidence of volcanic aerosol effect

The properties of the 1-micrometer volume extinction coefficient of two geographically different high-altitude cloud systems have been examined for the posteruption period (1985-1990) of the April 1982 El Chichon volcanic event with emphasis on the effect of volcanic aerosols on clouds. These two high-altitude cloud systems are the tropical clouds in the tropopause region observed by the Stratospheric Aerosol and Gas Experiment (SAGE) 2 and the polar stratospheric clouds (PSCs) sighted by the Stratospheric Aerosol Measurement (SAM) 2. The results indicate that volcanic aerosols alter the frequency distributions of these high-altitude clouds in such a manner that the occurrence of clouds having high extinction coefficients (6 x 10(exp -3) - 2 x 10(exp -2)/km) is suppressed, while that of clouds having low extinction coefficients (2 x 10(exp -3) - 6 x 10(exp -2)/km) is enhanced. This influence of the volcanic aerosols appears to be opposite to the increase in the extinction coefficient of optically thick clouds observed by the Earth Radiation Budget Experiment (ERBE) during the initial posteruption period of the June 1991 Pinatubo eruption. A plausible explanation of this difference, based on the Mie theory, is presented. As a consequence of the Mie theory, the effective radius of most, if not all, of the high-altitude clouds, measured by the SAGE series of satellite instruments must be less than about 0.8 micrometers. This mean cloud particle size implied by the satellite extinction-coefficient data at a single wavelength (1 micrometer) is further substantiated by the particle size analysis based on cloud extinction coefficient at two wavelengths (0.525 and 1.02 micrometers) obtained by the SAGE 2 observations. Most of the radiation measured by ERBE is reflected by cloud systems comprised of particles having effective radii much greater than 1 micrometer. A reduction in the effective radius of these clouds due to volcanic aerosols is expected to increase their extinction-coefficient values, opposite the effect observed by SAGE 2 and SAM 2. This work further illustrates the capability of the solar occultation satellite sensor to provide particulate extinction-coefficient measurements important to the study of the aerosol-cloud interactions. It is important to examine the variations of the extinction coefficient of these two high-altitude cloud systems for the posteruption years of the Pinatubo volcanic event for further evidence of the impact of volcanic aerosols on high-altitude clouds.

Wang, Pi-Huan↗

DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗

DOE OSTI2026

Single-shot ultrafast dynamics of nanosecond pulsed plasmas: Transition from ps–ns nonequilibrium to near-full ionization

The ultrafast multi-stage evolution of state-defining properties in atmospheric-pressure nanosecond pulsed plasmas is quantified using single-discharge, jitter-free, continuous streak-sweep spectroscopy of N 2 (C → B) molecular and N + /O + ionic species with a single-shot time resolution as short as ∼160 ps. An early period of extreme nonequilibrium is observed with vibrational temperatures [T V (C)] dropping from ∼8000 K at <200 ps after plasma breakdown to <4000 K within ∼1 ns, with near-ambient rotational temperatures [T R (C)] ∼ 300 K due to limited collisional energy transfer. This early-time T V (C) trend is representative of a direct and unquenched look into the high-energy tail of the electron energy distribution function; thus, it tracks real-time changes in the mean electron energy via the N2(C) emission signatures. This is followed by the rapid onset of N + /O + ionic emission after a distinct time delay of ∼14.9 ns due to multistep chemical kinetics. The ionic emission enables determination of electron densities (n e ) >2 × 10 19 cm −3 and electron temperatures (T e ) >36 000 K, indicating the transition to a nearly fully ionized regime. The ps–ns temporal dynamics are also compared between air and N 2 plasmas to assess the influence of collisional partners, as well as across the anode, cathode, and central gap regions to identify spatial variations in plasma behavior. Finally, this work demonstrates versatile continuously probing approach capable of spectro-temporal and spatially resolved single-shot measurements of key state variables in the ps–ns evolution of atmospheric-pressure nanosecond pulsed plasmas.

Electric discharges↗