Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “outflows”

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

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

At least 1,117 records · Page 62

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗

A large-scale comparison of Artificial Intelligence and Data Mining (AI&DM) techniques in simulating reservoir releases over the Upper Colorado Region

In recent years, the Artificial Intelligence and Data Mining (AI&DM) models have become popular tools in assisting various aspects of reservoir operation. However, the practical uses are still rarely reported. Comparison experiment of many AI&DM models over a large number of reservoir cases is particularly valuable to help reservoir operators first examine the usefulness and transferability of different AI&DM models, and then identify the most stable and reliable AI&DM model in assist of various decision-making processes. In this study, a total of 12 AI&DM models with different parameterizations and simulation scenarios are comprehensively tested out and compared in simulating the controlled reservoir outflows of 33 reservoir cases over the Upper Colorado Region, United States. Results show that the Random Forecast and the Long-Short-Term-Memory model could consistently derive the best statistical performance than other models under the baseline simulation scenario. The employed AI&DM models could obtain satisfactory statistical interquartile ranges (25–75%) between [0.6–0.9], [0.3–0.8], and [0.2–0.8], for CORR, NSE, and KGE measurements, respectively, and [1.5–6.5], [–15 to 20], and [0.5–8.5] for the normalized RMSE, PBIAS and RSR measurements, respectively. Results also show Multi-Layer Perceptron model and Extreme Gradient Boosting Tree Algorithm produced more stable and superior performance than other models under more complex input scenarios. We also found that the performance of different AI&DM models are closely relevant to the reservoir elevations, sizes, and functionalities. Discussions were made about the sensitivity of AI&DM models’ parameterizations and the key advantages of AI&DM models over the rule-based reservoir models. We further identify that the main advantage of AI&DM models is the flexibility in designing input structures, whereas the rule-based simulation model is rather limited. Future studies were suggested regarding the best way reservoir operators and researchers could use, select, and apply different AI&DM models in simulating reservoir releases under different natural and modeling environments. Finally, this comparison study also serves as a reference and a piece of groundwork for further promoting the practical uses of AI&DM models in assisting reservoir operation.

54 ENVIRONMENTAL SCIENCES↗

Estimation of process steady state with autoregressive models and Bayesian inference

To improve efficiency, separations engineers will typically design process circuits containing recirculating streams, which mix one or more of the process outputs with the feed material. Doing so can improve efficiency, but will cause a delay in the system reaching steady state conditions until the recirculating load mass flows stabilize. In testing separation circuits, engineers will often test a variety of factors and complete an analysis from sample results. Knowledge of if a process is at steady state, as well as the steady state conditions of a process, is essential for a valid techno-economic analysis. However, the definition of process steady state is often poorly defined, or does not include uncertainty quantification. If the performance of a process operating under two different sets of conditions are compared, an engineer who does not test for steady state or quantify steady state conditions risks producing a faulty analysis. In this work, a Bayesian statistical method for testing if all streams are at steady state is further motivated and then derived. Then after testing for steady state, the same model is used with a prior distribution that enforces a steady state assumption to estimate steady state conditions. Further, these methods were validated in a solvent extraction pilot plant where steady state conditions for all outflows were inferred with uncertainty quantification. Analysis is completed with functions available to the reader as part of the BayesMassBal (V 1.1.0) software package written in R.

01 COAL, LIGNITE, AND PEAT↗

AstroPix: A pixelated HVCMOS sensor for space-based gamma-ray measurement

A next-generation medium-energy gamma-ray telescope targeting the MeV range would address open questions in astrophysics regarding how extreme conditions accelerate cosmic-ray particles, produce relativistic jet outflows, and more. One concept, AMEGO-X, relies upon the mission-enabling CMOS Monolithic Active Pixel Sensor silicon chip AstroPix. AstroPix is designed for space-based use, featuring low noise, low power consumption, and high scalability. Desired performance of the device include an energy resolution of 5 keV (or 10% FWHM) at 122 keV and a dynamic range per-pixel of 25–700 keV, enabled by the addition of a high-voltage bias to each pixel which supports a depletion depth of 500 μ m. This work reports on the status of the AstroPix development process with emphasis on the current version under test, version three (v3), and highlights of version two (v2). Version 3 achieves energy resolution of 10.4 ± 3.2% at 59.5 keV and 94 ± 6 μ m depletion in a low-resistivity test silicon substrate.

Astrophysics instrumentation↗

Valued peaks: Sustainable water allocation for small hydropower plants in an era of explicit ecological needs

Optimizing hydropower operations to balance economic profitability and support functioning ecosystem services is integral to river management policy. In this article, we propose a dynamic, constrained optimization framework for small hydropower plants (SHPs) to evaluate trade-offs between economic profitability and socio-ecological requirements. Specifically, we examine the balance between short-term losses in hydropower generation and the potential for compensatory benefits in the form of revenue from recreational ecosystem services, irrespective of the direct beneficiary. Our framework integrates a fish habitat model, a hydropower optimization model, and a recreational ecosystem service estimate to evaluate different environmental flow scenarios. The optimization process gives three outflow release scenarios, informed by previous streamflow realisations (dam inflow), and designed environmental flow constraints. The framework is applied and tested for the river Kuusinkijoki in North-eastern Finland, which is a habitat for migratory brown trout and grayling populations. We show that the revenue loss due to the environmental flow constraints arises through a reduction in revenue per generated energy unit and through a reduction in turbine efficiency. Additionally, the simulation results reveal that all the designed environmental flow constraints cannot be met simultaneously. Under the environmental flow scenario with both minimum flow and flow ramping rate constraints, the annual hydropower revenue decreases by 16.5 %. An annual increase of 8 % in recreational fishing visits offsets the revenue loss. In conclusion, the developed framework provides knowledge of the costs and benefits of hydropower environmental flow constraints and guides the prioritizing process of environmental measures.

13 HYDRO ENERGY↗

An Alfvenic reconnecting plasmoid thruster

A new concept for the generation of thrust for space propulsion is introduced. Energetic thrust is generated in the form of plasmoids (confined plasma in closed magnetic loops) when magnetic helicity (linked magnetic field lines) is injected into an annular channel. Using a novel configuration of static electric and magnetic fields, the concept utilizes a current-sheet instability to spontaneously and continuously create plasmoids via magnetic reconnection. Here, the generated low-temperature plasma is simulated in a global annular geometry using the extended magnetohydrodynamic model. Because the system-size plasmoid is an Alfvenic outflow from the reconnection site, its thrust is proportional to the square of the magnetic field strength and does not ideally depend on the mass of the ion species of the plasma. Exhaust velocities in the range of 20 to $500\ \textrm {km}\ \textrm {s}^{-1}$ , controllable by the coil currents, are observed in the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Observation of subcritical shocks in a collisional laboratory plasma: scale dependence near the resistive length

We present a study of subcritical shocks in a highly collisional laboratory plasma with a dynamically significant magnetic field. Shocks were produced by placing cylindrical obstacles into the supermagnetosonic ( $M_{{\rm MS}} \sim 1.9$ ) outflow from an inverse wire array z-pinch at the MAGPIE pulsed power facility ( $n_{e} \sim 8.5 \times 10^{17}\,{\rm cm}^{-3}$ , $v \sim 45\,{\rm km}\,{\rm s}^{-1}$ ). We demonstrate the existence of subcritical shocks in this regime and find that secondary stagnation shocks form in the downstream which we infer from interferometry and optical Thomson scattering measurements are hydrodynamic in nature. The subcritical shock width is found to be approximately equal to the resistive diffusion length and we demonstrate the absence of a jump in hydrodynamic parameters. Temperature measurements by collective optical Thomson scattering showed little temperature change across the subcritical shock ( ${<}10\,\%$ of the ion kinetic energy) which is consistent with a balance between adiabatic and Ohmic heating and radiative cooling. We demonstrate the absence of subcritical shocks when the obstacle diameter is less than the resistive diffusion length due to decoupling of the magnetic field from the plasma. These findings are supported by magnetohydrodynamic simulations using the Gorgon and AstroBEAR codes and discrepancies between the simulations and experiment are discussed.

Physics↗

Reconstruction of Electron and Ion Distribution Functions in a Magnetotail Reconnection Diffusion Region

In the diffusion region of magnetotail reconnection, particle distributions are highly structured, exhibiting triangular shapes and multiple striations that deviate dramatically from the Maxwellian distribution. Fully kinetic simulations have been demonstrated to be capable of producing the essential structures of the observed distribution functions, yet are computationally not feasible for 3D global simulations. The fluid models used for large-scale simulations, on the other hand, do not have the kinetic physics necessary for describing reconnection accurately. Our study aims to bridge fully kinetic and fluid simulations by quantifying the information required to capture the non-Maxwellian features in the distributions underlying the closures used in the fluid code. Here, we compare the results of fully kinetic simulations with observed electron velocity distributions in a magnetotail reconnection diffusion region and use the maximum entropy model to reconstruct electron and ion distributions using various numbers of moments obtained from the simulation. Our results indicate that using only local moments, the maximum entropy model can reproduce many of the features of the distributions: (1) the electron outflow distribution with a tilted triangular structure is reproduced with 21 or more moments in agreement with Ng et al. and (2) counterstreaming distributions can be captured with the 35-moment model when the separation in velocity space between the populations is large.

58 GEOSCIENCES↗

Multipoint Observations of Dispersive Scale Alfvénic Field‐Line Resonances Associated With Substorm Auroral Beads

We present a case study of the field-aligned current (FAC) systems that transpire within the high-altitude auroral acceleration region of an “auroral bead” initiated double oval substorm observed on 23 February 2001 by the Cluster fleet. Conjunctive Cluster measurements and auroral images from IMAGE reveal that auroral bead current system formation and evolution is a multi-scale, injection-mediated process. The FACs at large scales vary on substorm evolution time scales (~minutes) in response to the injection and evolution of hotter denser magnetospheric plasma. Embedded within the large-scale FACs are intense short-scale (≲ few 10s of km) currents comprising dispersive scale Alfvén wave (DAW) fluctuations. The DAWs are a complex mixture of ingoing and reflected components that regularly interfere to form a broad spectrum of kinetic (dispersive) scale Alfvénic field-line resonances (KFLRs). The Alfvénic currents appear as a nested series of upward and downward FAC densities with amplitudes reaching a few 100 nA/m 2 . Energized field-aligned or counterstreaming electrons near keV energies and below are observed with parallel skews that vary in concert with variations in the DAW current sense. Positive correlations between DAW electric field energy densities and the energies of energized H + , He + , and O + outflow are observed, indicative of ion energization within the DAW fields. Due to their L-shell location (L~5.8–7.0) and associations with injections, the KFLRs are interpreted as the high-altitude auroral zone analog of KFLRs observed in the equatorial inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Reconnection Inside Solar Wind Rotational Discontinuity During Its Interaction With the Quasi-Perpendicular Bow Shock and Magnetosheath

Here, using a three-dimensional global hybrid simulation, we investigate the formation and evolution of ion-scale magnetic reconnection inside an interplanetary rotational discontinuity (RD) owing to its interaction with the quasi-perpendicular (Q-⟂) bow shock and the magnetosphere. The interplanetary magnetic field (IMF) is initially predominantly northward, while it changes to purely southward across the RD. A significantly thinned RD current layer with a width ~ ion skin depth and normal magnetic field B n ≈0 is formed by a shock compression process as the RD interacts with the Q-⟂ shock. Magnetic reconnection thereupon takes place inside the thin RD current layer, where Hall magnetic and electric fields, reconnection electric field, and high-speed ion outflow jets are identified. Simultaneously, flux ropes form with an extension of a few ion inertial lengths. As the RD is transmitted into the magnetosheath, multiple reconnection sites lead to the formation of longer flux ropes. Moreover, magnetosheath reconnection is also found at new reconnection sites inside the RD. No reconnection is found in the Q-⟂ shock alone outside the RD. The flux ropes propagate poleward and tailward in the magnetosheath. In addition, magnetopause reconnection takes place under the southward IMF on the sunward side of the RD after the RD passes through the magnetopause.

79 ASTRONOMY AND ASTROPHYSICS↗

Lower‐Hybrid Wave Structures and Interactions With Electrons Observed in Magnetotail Reconnection Diffusion Regions

Abstract We investigate waves close to the lower‐hybrid frequency in 12 magnetotail reconnection electron diffusion region (EDR) events with guide field levels of near‐zero to 30%. In about half of the events, the wave vector has a small component along the current sheet normal, consistent with known lower‐hybrid drift wave properties, but the perpendicular magnetic field fluctuations can be comparable or greater than the parallel component, a feature unique to the waves inside and adjacent to EDRs. Another new wave property is that the wave vector has a significant component along the current sheet normal in some events and completely along the normal for one event. In 1/4 of the events, the term has a significant contribution to the wave electric field, possibly a feature of lower‐hybrid waves more likely to exist in the diffusion region than further away from the X‐line. Electron temperature variations are correlated with the wave potential, due to wave electric field acceleration and crossings at the corrugated separatrix region with different amounts of mixing between reconnection inflowing and outflowing populations. The latter also leads to the anti‐correlation between parallel and perpendicular temperature components. Using four‐spacecraft measurements, the magnetic field line twisting is demonstrated by the correlated fluctuations in and . The lower‐hybrid wave in the EDR of weak guide field reconnection may be generated near separatrices and penetrate to the mid‐plane or locally generated, and the latter possibility is beyond the prediction of previous reconnection simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaling of Electron Heating by Magnetization During Reconnection and Applications to Dipolarization Fronts and Super-Hot Solar Flares

Electron ring velocity space distributions have previously been seen in numerical simulations of magnetic reconnection exhausts and have been suggested to be caused by the magnetization of the electron outflow jet by the compressed reconnected magnetic fields (Shuster et al., 2014, https://doi. org/10.1002/2014GL060608). Here, we present a theory of the dependence of the major and minor radii of the ring distributions solely in terms of upstream (lobe) plasma conditions, thereby allowing a prediction of the associated temperature and temperature anisotropy of the rings in terms of upstream parameters. We test the validity of the prediction using 2.5-dimensional particle-in-cell (PIC) simulations with varying upstream plasma density and temperature, finding excellent agreement between the predicted and simulated values. We confirm the Shuster et al. suggestion for the cause of the ring distributions, and also find that the ring distributions are located in a region marked by a plateau, or shoulder, in the reconnected magnetic field profile. The predictions of the temperature are consistent with observed electron temperatures in dipolarization fronts, and may provide an explanation for the generation of plasma with temperatures in the 10s of MK in super-hot solar flares. A possible extension of the model to dayside reconnection is discussed. Since ring distributions are known to excite whistler waves, the present results should be useful for quantifying the generation of whistler waves in reconnection exhausts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating the Transport and Rupture of Pollen in the Atmosphere

Pollen, one type of primary biological aerosol particle (PBAP), is emitted from the terrestrial biosphere and can undergo physical changes in the atmosphere via particle rupture. To examine the fate of pollen and its atmospheric processing, a pollen emission and transport scheme is coupled to the Weather Research and Forecasting Model with Chemistry (WRF-Chem). We simulate the emission of pollen and its impacts on the cloud properties and precipitation in the Southern Great Plains from 12 to 19 April 2013, a period with both high pollen emissions and convective activity. We conduct a suite of ensemble runs that simulate primary pollen and three different pollen rupture mechanisms that generate subpollen particles, including (a) high humidity-induced surface rupture, (b) high humidity-induced in-atmosphere plus surface rupture, and (c) lightning-induced rupture, where in-cloud and cloud-to-ground lightning strikes trigger pollen rupture events. When relative humidity is high (>80%), coarse primary pollen (~1 μg m –3 ) is converted into fine subpollen particles (~1.2e –4 μg m –3 ), which produces 80% more subpollen particles than lightning-induced rupture. The in-atmosphere humidity-driven rupture predominantly produces subpollen particles, which is further enhanced during a frontal thunderstorm. During strong convection, vertical updrafts lift primary pollen and subpollen particles (~0.5e –4 μg m –3 ) to the upper troposphere (~12 km) and laterally transports the ruptured pollen in the anvil top outflow. In regions of high pollen and strong convection, ruptured pollen can influence warm cloud formation by decreasing low cloud (<4 km) cloud water mixing ratios and increasing ice phase hydrometeors aloft (>10 km)

54 ENVIRONMENTAL SCIENCES↗

Congo Basin Water Balance and Terrestrial Fluxes Inferred From Satellite Observations of the Isotopic Composition of Water Vapor

Large spatio-temporal gradients in the Congo basin vegetation and rainfall are observed. However, its water-balance (evapotranspiration minus precipitation, or ET - P) is typically measured at basin-scales, limited primarily by river-discharge data, spatial resolution of terrestrial water storage measurements, and poorly constrained ET. We use observations of the isotopic composition of water vapor to quantify the spatio-temporal variability of net surface water fluxes across the Congo Basin between 2003 and 2018. These data are calibrated at basin scale using satellite gravity and total Congo river discharge measurements and then used to estimate time-varying ET - P over four quadrants representing the Congo Basin, providing first estimates of this kind for the region. We find that the multi-year record, seasonality, and interannual variability of ET - P from both the isotopes and the gravity/river discharge based estimates are consistent. Additionally, we use precipitation and gravity-based estimates with our water vapor isotope-based ET - P to calculate time and space averaged ET and net river discharge within the Congo Basin. These quadrant-scale moisture flux estimates indicate (a) substantial recycling of moisture in the Congo Basin (temporally and spatially averaged ET/P > 70%), consistent with models and visible light-based ET estimates, and (b) net river outflow is largest in the Western Congo where there are more rivers and higher flow rates. Our results confirm the importance of ET in modulating the Congo water cycle relative to other water sources.

54 ENVIRONMENTAL SCIENCES↗

Radiocarbon Fingerprinting Black Carbon Source History in the Himalayas

Abstract Black carbon (BC) is considered as an important contributor to the Himalayan glaciers melt in the past few decades. However, the long‐term source apportionment of BC remains unclear. Here we present the first radiocarbon ( 14 C)‐based annual variation of BC source apportionment in an ice core spanning the period of 1959–2012 drilled from the Southeastern Tibetan Plateau, a receptor site of South Asia outflow. We find fossil fuel combustion is a major contribution (73% ± 5%), yet the biomass burning fraction ( ƒ biomass ) has grown from 24% ± 4% to 30% ± 4% since 1990. Intriguingly, we further find the ƒ biomass demonstrating a robust correlation with South Asian wildfires linked to climate oscillations. Thus, for mitigating BC impacts on Himalayan glaciers, South Asia's transition from fossil fuels to clean energy is a more efficient and urgent strategy than reducing residential biomass burning.

Wang, Mo [State Key Laboratory of Tibetan Plateau ↗

Impact of the Out‐Of‐Plane Flow Shear on Magnetic Reconnection at the Flanks of Earth's Magnetopause

Abstract Magnetic reconnection changes the magnetic field topology and facilitates the energy and particle exchange at magnetospheric boundaries such as the Earth's magnetopause. The flow shear perpendicular to the reconnecting plane prevails at the flank magnetopause under southward interplanetary magnetic field conditions. However, the effect of the out‐of‐plane flow shear on asymmetric reconnection is an open question. In this study, we utilize kinetic simulations to investigate the impact of the out‐of‐plane flow shear on asymmetric reconnection. By systematically varying the flow shear strength, we analyze the flow shear effects on the reconnection rate, the diffusion region structure, and the energy conversion rate. We find that the reconnection rate increases with the upstream out‐of‐plane flow shear, and for the same upstream conditions, it is higher at the dusk side than at the dawn side. The diffusion region is squeezed in the outflow direction due to magnetic pressure which is proportional to the square of the Alfvén Mach number of the shear flow. The out‐of‐plane flow shear increases the energy conversion rate , and for the same upstream conditions, the magnitude of is larger at the dusk side than at the dawn side. This study reveals that out‐of‐plane flow shear not only enhances the reconnection rate but also significantly boosts energy conversion, with more pronounced effects on the dusk‐side flank than on the dawn‐side flank. These insights pave the way for better understanding the solar wind‐magnetosphere interactions.

Liang, Haoming↗

Multi‐Species Energy‐Banded Ions in the Ionosphere During the 21 January 2005 Magnetic Storm: Low‐Altitude Edge of the Warm Plasma Cloak

During the 21–22 January 2005 magnetic storm, the FAST satellite observed warm (< few keV) ions in discrete energy bands on the dayside at ∼3,000 km altitude for more than 6.5 hr. We suggest that the ionospheric energy-banded ions represent the low-altitude edge of the warm plasma cloak observed simultaneously by magnetospheric satellites. This is a clear example of the multi-species ion energy bands (10 eV to several keV) observed during strong magnetic storms by the FAST satellite, stretching from the diffuse auroral region to the plasmapause with lifetimes up to 12 hr. The close association of these energy-banded ions with magnetic storms, their broad latitudinal extent, and the presence of multiple ion species in the same energy band, rather than at the same velocity, indicate that this is a distinct phenomenon from other types of energy-banded ions. During the 21–22 January 2005 magnetic storm, the dayside ion energy band structures, centered at 10 eV (H + ), 40 eV (H + and He + ), and 160 eV (H + , He + , and O + ), were consistent with a “time-of-flight and velocity filter” formation process acting on a near-cusp, impulsive outflow of a <200 eV multi-species ion-source population, poleward and in the same hemisphere as FAST. Understanding the sources and dynamics of warm energy-banded ions and their linkage to the warm plasma cloak is important because during superstorms these ions are transported to L values as low as L ∼ 1.2 in the dawn sector, significantly altering the energetics of the mid-latitude ionosphere.

58 GEOSCIENCES↗

Quantifying UAS Observation Error Variance Used in Data Assimilation Systems and Its Impact on Predictive Skill

Observation error determines the weights of the observations and background state used in data assimilation to generate analyses. Quantifying observation error is critical for the optimal assimilation of observational data sets. Uncrewed Aircraft System (UAS) observations have shown potential benefits in filling observational gaps in the lower atmosphere; however, characterization of their error characteristics has been limited. To optimize the use of UAS observations in numerical weather prediction, UAS observation error is estimated based on the 3‐cornered hat diagnostic approach which uses three independent estimates of the atmospheric state. This approach is applied to data from the 2018 Lower Atmospheric Profiling Studies at Elevation‐a Remotely‐piloted Aircraft Team Experiment field campaign using collocated UAS and rawinsonde observations along with output from a set of convection‐permitting model simulations. The estimated observation error values for UAS temperature, wind, and relative humidity measurements were found to be only weakly dependent on height AGL with mean values equal to 0.5°C, 0.8 m s −1 , and 3%, respectively. Only the newly estimated observation error for temperature differed from that previously used to assimilate commercial aircraft observations into global models (1.0°C). However, using this reduced temperature observation error produced more accurate mesoscale analyses and forecasts of both terrain‐driven flows and convection initiation generated by colliding outflow boundaries within the San Luis Valley of Colorado.

54 ENVIRONMENTAL SCIENCES↗