SEARCH · Engineering Papers

Results for “velocity steps”

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

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on ~30 m range gates, stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, below range, ran out of signal, cloud-topped). Cloud Base Height (Haar-gradient detection): 15 min estimates of cloud-base height (m) with a cloud-detection quality flag (0–3: none, low, moderate, high). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (2.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution, with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗

DOE OSTI2022

Commensuration effects on skyrmion Hall angle and drag for manipulation of skyrmions on two-dimensional periodic substrates

Here, we examine the dynamics of an individually driven skyrmion moving through a background lattice of skyrmions coupled to a 2D periodic substrate as we vary the ratio of the number of skyrmions to the number of pinning sites across commensurate and incommensurate conditions. As the skyrmion density increases, the skyrmion Hall angle is nonmonotonic, dropping to low or zero values in commensurate states and rising to an enhanced value in incommensurate states. Under commensuration, the driven skyrmion is channeled by a symmetry direction of the pinning array and exhibits an increased velocity. At fillings for which the skyrmion Hall angle is zero, the velocity has a narrow band noise signature, while for incommensurate fillings, the skyrmion motion is disordered and the velocity noise is broad band. Under commensurate conditions, multi-step depinning transitions appear and the skyrmion Hall angle is zero at low drives but becomes finite at higher drives, while at incommensurate fillings there is only a single depinning transition. As the gyrotropic component of the skyrmion dynamics, called the Magnus force, increases, peaks in the velocity that appear in commensurate regimes cross over to dips, and new types of directional locking effects can arise in which the skyrmion travels along other symmetry directions of the background lattice. At large Magnus forces, and particularly at commensurate fillings, the driven skyrmion can experience a velocity boost in which the skyrmion moves faster than the applied drive due to the alignment of the Magnus-induced velocity with the driving direction. In some cases, an increase of the Magnus force can produce regimes of enhanced pinning when the skyrmion is forced to move along a nonsymmetry direction of the periodic pinning array. This is in contrast to systems with random pinning, where increasing the Magnus force generally reduces the pinning effect. We demonstrate these dynamics for both square and triangular substrates and map out the different regimes as a function of filling fraction, pinning force, and the strength of the Magnus force in a series of dynamic phase diagrams.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on 30 m range gates (and 3 m range gates), stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, ran out of signal, below range, cloud-topped). Cloud Base Height (Haar-gradient detection): 10 min estimates of cloud-base height (m). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (3.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution (and 3 m for the year of 2025), with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min (10 min for Cloud Heights) summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗

DOE OSTI2026

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on 30 m range gates (and 3 m range gates), stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, ran out of signal, below range, cloud-topped). Cloud Base Height (Haar-gradient detection): 10 min estimates of cloud-base height (m). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (3.0.1), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution (and 3 m for the year of 2025), with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min (10 min for Cloud Heights) summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗

DOE OSTI2021

Towards cost-competitive middle distillate fuels from ethanol within a market-flexible C2 platform-based biorefinery concept

Ethanol to middle distillates (ETMD) is a promising pathway to produce sustainable liquid fuels to decarbonize the hard-to-electrify transportation sectors due to (1) the abundant sugar/starch and lignocellulosic biomass, (2) the existing deployment scale of fuel ethanol production (similar to 29 billion gallons per year globally), and (3) emerging opportunities in C2+ alcohol synthesis from CO2. Here we report a conceptual market-responsive biorefinery centered around a new ETMD pathway based on one-step ethanol to butene-rich olefins (ETO) over a Cu-Zn-Y/Beta catalyst. Specifically, this ethanol conversion pathway comprises one-step ETO, oligomerization, and hydrotreating. This ETO is distinct from that in the conventional ethanol-to-jet process which is based on two-step ethanol to ethylene and ethylene oligomerization to butenes. Butene-rich olefins can be shifted to butadiene-rich products by simply changing the reaction atmosphere from hydrogen to inert gas over the same ETO catalyst. Leveraging the experimental results, baseline techno-economic analysis (TEA) and sensitivity analysis indicate that the ethanol conversion cost is $0.60 per gallon gasoline equivalent (GGE), with opportunities for further cost reduction via improving the liquid hydrocarbon yield and space velocities, and process optimization on balancing dewatering of ethanol feed prior to the ETO step. The minimum fuel selling price (MFSP) of liquid hydrocarbons derived from corn starch ethanol with butadiene as coproduct is $1.64 per GGE, in the range that is cost competitive with petroleum kerosene-type jet fuel. Projected MFSP for cellulosic ethanol (corn stover) derived hydrocarbons is below $3.00 per GGE and co-production of butadiene further reduces the MFSP to $1.70 per GGE. The Well-to-Wake life-cycle analysis indicates that 85% greenhouse gas emission reduction can be achieved when using corn stover compared to petroleum reference and the associated carbon credits will provide significant economic incentives to favor the cellulosic ethanol-derived hydrocarbon fuels. This study demonstrates a low-cost pathway to middle distillate fuels leveraging existing ethanol infrastructure, where catalysis innovation drives the reduction of process complexity and flexible coproduction of a value-added chemical product.

Zhang, Junyan↗

DOE OSTI2021

Type Ia supernova magnitude step from the local dark matter environment

ABSTRACT Residuals in the Hubble diagram at optical wavelengths and host galaxy stellar mass are observed to correlate in type Ia supernovae (SNe Ia) (‘magnitude step’). Among possible progenitor channels for the associated explosions, those based on dark matter (DM) have attracted significant attention, including our recent proposal that ‘normal’ SNe Ia from bare detonations in sub-Chandrasekhar white dwarf stars are triggered by the passage of asteroid-mass primordial black holes (PBHs): the magnitude step could then originate from a brightness dependence on stellar properties, on DM properties, or both. Here, we present a method to estimate the local DM density and velocity dispersion of the environment of SN Ia progenitors. We find a luminosity step of $0.52\pm 0.11\,$ mag corresponding to bins of high versus low DM density in a sample of 222 low-redshift events from The Open Supernova Catalog. We investigate whether the magnitude step can be attributed to local DM properties alone, assuming asteroid-mass PBHs. Given the inverse correlation between SN Ia brightness and PBH mass, an intriguing explanation is a spatially inhomogeneous PBH mass function. If so, a strong mass segregation in the DM density-dependent PBH mass scale is needed to explain the magnitude step. While mass segregation is observed in dense clusters, it is unlikely to be realized on galactic scales. Therefore, if DM consists of asteroid-mass PBHs, the magnitude step is more likely to exist, and dominantly to be attributed to local stellar properties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Precursory Off-Fault Deformation in Restraining and Releasing Step Overs: Insights From Discrete Element Method Models

Accelerating geophysical activity is detected preceding some, but not all, large earthquakes. This observation may indicate that no precursors occur before some earthquakes, or that the instrumentation lacks the required sensitivity. Here, to aid crustal monitoring efforts, we use discrete element method models to identify the locations and styles of deformation that may provide useful information about approaching fault reactivation. We model the reactivation of two healed rough faults in a variety of step over configurations, embedded in a host rock with varying amounts of damage subject to shear velocity loading parallel to the faults. Both the fault geometry and ratio of fault to host rock strength control the amount of off-fault deformation. Consistent with field observations, models with larger steps and more preexisting host rock damage produce higher amounts of off-fault deformation. We assess the size of the continuous regions of high velocities and strains to compare the value of the precursory information of each velocity and strain component. Comparing the three components of the velocity vector suggests that the fault-parallel velocity produces the largest and most temporally continuous regions of elevated velocity. The size of these regions increases toward failure, indicating the usefulness of tracking this component. Comparing the volumetric and shear components of the three-dimensional strain tensor suggests that during most of the interseismic period, the shear strain provides more information about approaching fault slip than the volumetric strain. However, in the days and months preceding fault reactivation, both the shear and volumetric strains provide similarly valuable information.

58 GEOSCIENCES↗

DOE OSTI2025

A kinetic line-driven radiation operator and its application to Gyrokinetics

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $\nabla_v \cdot(v\nu(v)f(v))$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $\nu(v)$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational cost of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.

kinetic↗

DOE OSTI2021

Individual Wave Detection and Tracking within a Rotating Detonation Engine through Computer Vision Object Detection applied to High-Speed Images

Known for their simplistic design and continuous detonation, rotating detonation engines (RDEs) constitute a majority of current pressure gain combustion (PGC) research efforts. Experimental RDE operation times have been continuously extended through the use of rig cooling techniques. As the window of observable behavior is expanded, and as the technology matures toward eventual integration within gas turbines, monitoring techniques must evolve to better match industrial diagnostics. High-speed image analysis techniques prove useful to capture and evaluate the unsteady detonation behavior within the RDE. Traditional image analysis techniques, however, require extensive processing times which prohibit simultaneous monitoring. To better address this problem, a computer vision object detection methodology is proposed to quickly detect individual detonation waves within a single down-axis image. Detonation waves are detected in individual images by the implemented computer vision method You Only Look Once (YOLO) object detection network. In order to detect detonation waves, the network must first be trained using RDE images of interest, for which each required phase of network development is outlined. Detection of waves is improved through proper treatment of the collected image set, variation of Intersection over Union (IoU) and confidence thresholding, and through a parametric study of annotation dimensions. Each detected wave is described by its location and rotational direction, and locations are tracked to calculate wave velocity across each frame, leading to a timestep resolution of 20 µs. Wave velocities are also calculated through a series of frames, leading to a suitable average velocity estimation using as few as 10 frames. Uncertainty analysis accounting for variation in camera framerate, pixel width and annotation centroid locations estimates a total uncertainty of ±4.3% for velocity calculations, using the smallest annotation boxes. This method offers great reductions in processing times, as a step toward real-time monitoring of detonation waves within an RDE. Improving on previous studies, this technique is impartial to wave modes not included in the original training set and calculates wave velocities independent of high-speed pressure data. The ability to isolate waves within predicted bounding boxes will likely facilitate analysis of pixel intensity variation as an estimation of wave strength in future work.

Johnson, Kristyn↗

DOE OSTI2022

Electrical manipulation of skyrmions in a chiral magnet

Writing, erasing and computing are three fundamental operations required by any working electronic device. Magnetic skyrmions could be essential bits in promising in emerging topological spintronic devices. In particular, skyrmions in chiral magnets have outstanding properties like compact texture, uniform size, and high mobility. However, creating, deleting, and driving isolated skyrmions, as prototypes of aforementioned basic operations, have been a grand challenge in chiral magnets ever since the discovery of skyrmions, and achieving all these three operations in a single device is even more challenging. Here, by engineering chiral magnet Co 8 Zn 10 Mn 2 into the customized micro-devices for in-situ Lorentz transmission electron microscopy observations, we implement these three operations of skyrmions using nanosecond current pulses with a low current density of about 10 10 A·m -2 at room temperature. A notched structure can create or delete magnetic skyrmions depending on the direction and magnitude of current pulses. We further show that the magnetic skyrmions can be deterministically shifted step-by-step by current pulses, allowing the establishment of the universal current-velocity relationship. These experimental results have immediate significance towards the skyrmion-based memory or logic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Deactivation of Cu/ZSM-5 Catalysts during the Conversion of 2,3-Butanediol to Butenes

This work determines the deactivation mechanisms of Cu/ZSM-5 catalysts used for the conversion of 2,3-butanediol to butene as part of an alcohol-to-jet route. The deactivation of the catalyst, reflected by a drop in the rate of the limiting hydrogenation step by over 90% in 24 h at a weight hourly space velocity of 5.92 h−1, proceeds via both the agglomeration of copper particles and the obstruction of copper surfaces due to carbonaceous deposits, although the former has less impact on the decrease in the hydrogenation rate. To reduce the detrimental effect of carbonaceous deposits on catalytic activity, ZMS-5 is modified through desilication of the HZSM-5 support with NaOH and CsOH solutions to generate a hierarchical structure with mesopores. The catalyst with the CsOH-treated support generates the highest overall yield of desired olefin products and experiences the slowest deactivation. This is a result of the lower Brønsted acidity and larger mesopores found in the CsOH-treated catalyst, leading to the slower formation of carbonaceous deposits and the faster diffusion of their precursors out of the pores.

Chemistry↗

DOE OSTI2024

Flashback Studies of High-Hydrogen Flames using High-Speed OH Planar Laser-Induced Fluorescence

Flame flashback is a crucial challenge during the application of high-hydrogen fuels for the development of next-generation lean premixed combustors. The fundamentals of stable flame configuration, flashback phenomena, and spatiotemporal information of flame dynamics in a premixed, swirl-stabilized burner are investigated experimentally at atmospheric pressure conditions using high-repetition-rate nanosecond (ns)-duration hydroxyl radical planar laser-induced fluorescence (OH-PLIF) measurements. The inlet conditions are varied systematically with respect to equivalence ratio (ϕ), hydrogen percentage in a hydrogen-methane mixture (50%–90%) and pre-mixer velocities and the experiments are performed on the burner modified to provide optical access to the premixing section. Signal interferences such as saturation effects, background chemiluminescence, and laser non-uniformity issues were minimized or corrected during data acquisition and subsequent data processing steps. A stable flame characterization was performed at 20 kHz-repetition-rate for varied pre-mixer velocities, hydrogen content and ϕ in the reactant mixture and lift-off length (L) investigation showed dependence on each of these parameters. Flame dynamics such as flame flashback events, local ignition, local extinction, flame curvature as well as the main reaction zone of the combustion mechanism are well characterized. Flashback ϕ showed an expected linearly increasing trend with increasing pre-mixer velocity and decreasing hydrogen content and the conclusions drawn aligned well with detailed L investigation. Such spatially and temporally resolved kHz-rate OH-PLIF is a promising technique to observe rapidly occurring flashback dynamics and can further validate turbulence-chemistry interaction models of swirling flames.

Parajuli, Pradeep↗

DOE OSTI2025

Heterogeneous Dynamics in Shear Thickening Colloids Revealed by Intrinsic Heterodyne Correlation Spectroscopy

Shear induced frictional networks have been proposed to be responsible for the emergence of discontinuous shear thickening (DST) in complex fluids. However, little experimental evidence exists to support this model directly. Here, in this study, using x-ray photon correlation spectroscopy (XPCS), we show the existence of an intrinsic heterodyne feature during shear cessation, which originates from the relative motion of mobile particles against an aggregated or jammed network induced by shear thickening. Upon removing the shear, the shear stress dissipates rather quickly in a two-step fashion, whereas the heterogeneous particle dynamics persist much longer with the relative velocity decaying slowly with time as 𝑡 −1 . More importantly, both continuous shear thickening (CST) and DST show similar heterodyne features, indicating the intrinsic mechanisms causing shear thickening are similar in nature

Complex fluid↗

DOE OSTI2022

Development, characterization, and modeling of a high-performance Ru/B2CA catalyst for ammonia synthesis

This paper documents the development and performance of a nano-phase Ru catalyst on a (BaO) x (CaO) y (Al 2 O 3 ) support. Extensive screening of the support’s ternary composition shows the best stoichiometry is (BaO) 2 (CaO)(Al 2 O 3 ), denoted B2CA. The paper first describes catalyst preparation and characterization. The paper reports a detailed 12-step reaction mechanism that represents ammonia synthesis over wide ranges of temperature, pressure, space velocity, and feed composition. Additionally, the mechanism is developed and validated using results of packed-bed experiments. The elementary reaction pathways consider surface adsorbates, including catalyst-poisoning behaviors. The rate expressions include important coverage-dependent activation barriers. Machine learning models assist interpretation of the catalyst-support interactions. The detailed chemistry is much more predictive than is possible with global representations (N 2 + 3H 2 ⇌ 2NH 3 ). The validated models can be applied to assist optimizing reactor design and operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Generation and Life Cycle of Solar Spicules

The physical mechanism for the creation of solar spicules is proposed with three stages of their life cycle. It is assumed that at stage I the density hump is formed locally in the x-y plane in the lower chromosphere in the presence of temperature gradients of electrons and ions along the z-axis (the vertical direction). In this region, the density structure of quasi-neutral (n i ≃ n e = n) plasma after taking birth is accelerated in the vertical direction owing to the thermal force F th ∝ ∇n(x, y, t) × (∇T e + ∇T i ). The exact time-dependent analytical solution of two-fluid plasma equations is presented assuming that density is maximum at the center of the density structure and decays away from it gradually. The 2D density structure is created as a step function H(t) in time at the bottom of the chromosphere, and consequently, the vertical plasma velocity turns out to be the ramp function of time R(t) = tH(t), whereas the source term S(x, y, t) for the density follows the delta function δ(t) form. The upward acceleration a = a(x, y) $\hat{z}$ produced in this density structure is greater than the downward constant solar acceleration −g ⊙ in the chromosphere. In the transition region, the temperature gradients are steeper; therefore, the upward acceleration increases in magnitude g ⊙ ≪ a and the density hump spends less time there. This is stage II of its life cycle. In stage III, the density structure enters into the corona, where the gradients of temperatures vanish and the structure decelerates to zero velocity under the action of the solar gravitational force.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗