Results for “Source and propagation modeling” — Page 38

DOE OSTI2022

Standoff Infrared Measurements of Chemical Plume Dynamics in Complex Terrain Using a Combination of Active Swept-ECQCL Laser Spectroscopy with Passive Hyperspectral Imaging

Chemical plume detection and modeling in complex terrain present numerous challenges. We present experimental results from outdoor releases of two chemical tracers (sulfur hexafluoride and Freon-152a) from different locations in mountainous terrain. Chemical plumes were detected using two standoff instruments collocated at a distance of 1.5 km from the plume releases. A passive long-wave infrared hyperspectral imaging system was used to show time- and space-resolved plume transport in regions near the source. An active infrared swept-wavelength external cavity quantum cascade laser system was used in a standoff configuration to measure quantitative chemical column densities with high time resolution and high sensitivity along a single measurement path. Both instruments provided chemical-specific detection of the plumes and provided complementary information over different temporal and spatial scales. The results show highly variable plume propagation dynamics near the release points, strongly dependent on the local topography and winds. Effects of plume stagnation, plume splitting, and plume mixing were all observed and are explained based on local topographic and wind conditions. Measured plume column densities at distances ~100 m from the release point show temporal fluctuations over ~1 s time scales and spatial variations over ~1 m length scales. The results highlight the need for high-speed and spatially resolved measurement techniques to provide validation data at the relevant spatial and temporal scales required for high-fidelity terrain-aware microscale plume propagation models.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Insights Into Lightning K‐Leader Initiation and Development From Three Dimensional Broadband Interferometric Observations

Abstract We report detailed observations of K‐leaders and the activity between them with the three‐dimensional Broadband Interferometric Mapping and Polarization system (BIMAP‐3D) at Los Alamos National Laboratory. It is found that K‐leaders have a general propagation trend of initial acceleration and then gradual deceleration, and the corresponding very high frequency (VHF) radiation power is exponentially correlated with the leader speed. Based on the 3D development and simultaneous electric field change measurement, some simple K‐leaders can be modeled with time‐evolving point charges at the propagating leader tip and at the stationary origin. We found that the charge magnitude increases during the initial acceleration stage and stays relatively constant for the rest of the development. K‐leaders are observed to interact with other branches; the branches affect the leader's propagation speed, and may affect the charge transfer. After the occurrence of a K‐leader, VHF emissions are quenched for several milliseconds. VHF sources then reappear in an impulsive and scattered manner as “twinkling,” and these sources are found not uniquely on the so‐called needles, but also on the main channel. These twinkling sources start near the apparent positive leader tip, and extend back toward the direction of the flash origin at about 10 5 m/s, while the apparent positive tip continues to extend forward at about 10 4 m/s. The twinkling extending toward the direction of flash origin appears to initiate the following K‐leader, although it may be interrupted by a K‐leader along a different branch, or simply die out without more K‐leader activity.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1972

The plasma physics of the Jovian decameter radiation.

We have assumed that the decameter radiation from Jupiter is produced near the local electron gyrofrequency and is amplified as it propagates out of the Jovian magnetosphere. We have derived the growth rate for radiation that propagates almost perpendicular to the direction of the magnetic field. When the electrons are described by a loss-cone distribution function, the growth rate is large enough to lead to a large amplification factor over a source of 100-4000 km, depending on the choice of parameters. Because we expect low-energy electrons to be trapped in the Jovian dipole field regardless of the position of the satellite Io, we maintain that this model provides a plausible mechanism for the decametric radiation not associated with Io.

Goldstein, M. L.↗

DOE OSTI2025

CONFLUX: A standardized framework to calculate reactor antineutrino flux

Nuclear fission reactors are abundant sources of antineutrinos for neutrino physics experiments. The flux and spectrum of antineutrinos emitted by a reactor can indicate its activity and composition, suggesting potential applications of neutrino measurements beyond fundamental scientific studies that may be valuable to society. The utility of reactor antineutrinos for applications and fundamental science is dependent on the availability of precise predictions of these emissions. For example, in the last decade, disagreements between reactor antineutrino measurements and models have inspired revision of reactor antineutrino calculations and standard nuclear databases as well as searches for new fundamental particles not predicted by the Standard Model of particle physics. Past predictions and descriptions of the methods used to generate them are documented to varying degrees in the literature, with different modeling teams incorporating a range of methods, input data, and assumptions. The resulting difficulty in accessing or reproducing past models and reconciling results from differing approaches complicates the future study and application of reactor antineutrinos. The CONFLUX (Calculation Of Neutrino FLUX) software framework is a neutrino prediction tool built with the goal of simplifying, standardizing, and democratizing the process of reactor antineutrino flux calculations. CONFLUX includes three primary methods for calculating the antineutrino emissions of nuclear reactors or individual beta decays that incorporate common nuclear data and beta decay theory. The software is prepackaged with the current nuclear databases, including ENDF.B/VIII, JEFF-3.3, and ENSDF, and it includes the capability to predict time-dependent reactor emissions, adjust nuclear database or beta decay inputs/assumptions, and propagate related sources of uncertainty. Here, this paper describes the CONFLUX software structure, details the methods used for flux and spectrum calculations, and provides examples of potential use cases.

Zhang, Xianyi [Lawrence Livermore National Laborat↗

NASA NTRSDate not supplied

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

NASA NTRSDate not supplied

Assimilation of SMAP Observations Over Land Improves the Simulation and Prediction of Tropical Cyclone Idai

This work is focused on the role of soil moisture in the prediction of tropical cyclones (TCs) approaching land and after landfall. Soil moisture conditions can impact the circulation and structure of an existing tropical cyclone (TC) when part or all of the circulation is over land. For example, dry land surface conditions may lead to faster dissipation of a TC over land (often associated with changes in precipitation structure), whereas very wet conditions may help sustain or in rare cases re-intensify a TC. Moreover, the presence of strong soil moisture gradients may affect the symmetry and development of the TC circulation leading to changes in its over-land track. While the link between soil moisture conditions and TC evolution in proximity to land is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tb observations significantly improves modeled land surface states. Thus, it is expected that SMAP can be used to constrain land surface initial conditions and potentially benefit TC forecasts. Here we present two sets of retrospective forecasts of TC Idai that are compared in an Observing System Experiment framework at ¼ degree resolution: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP brightness temperature observations over land using a weakly-coupled land analysis. We find that the assimilation of SMAP meaningfully improves the representation of TC Idai’s structure as well as the prediction of its intensity and track. The analyzed TC size, as measured by the wind speed radius, is improved by up to 18% in the analysis with SMAP assimilation relative to the control run. The forecast intensity error, measured against the observed intensity, is reduced by up to 23%. At the 1/4-degree resolution used here, GEOS unavoidably under-estimates TC intensity and over-estimates TC size. The SMAP assimilation therefore corrects the model in the right direction, leading to a storm that is more energetic and more compact. Furthermore, we find that the along-track forecast error is reduced by up to 34%, indicating a more accurate propagation speed, which is consistent with the fact that TC speed over land is strongly affected by surface processes. The impact of SMAP assimilation on the forecast cross-track error is neutral. Across the TC forecast skill metrics used here, the improvements from SMAP DA are largest at lead times of 36 to 72 hours, suggesting that the predictability of forecasts at shorter lead times may be dominated by short-term convective processes, while the land and its longer memory gains in importance as a source of predictability on a 2-3 day timescale. We further investigated the underlying mechanisms leading to the skill improvements from SMAP data assimilation by isolating the land areas that directly influence TC Idai using a back trajectory analysis. We find that the assimilation of SMAP leads to wetter soil moisture conditions that cause an increased latent heat flux, which ultimately results in TC analyzed representation that has higher column-integrated total moisture content and total energy compared to the analysis in the control run without SMAP assimilation. Overall, the results highlight that the assimilation of SMAP observations into a global numerical weather prediction model can lead to pronounced improvements of TC predictions. This is a crucial step towards a better mitigation of the socio-economic impact of landfalling TCs and thus safeguarding human lives. Finally, our study presents an event-based approach that assesses the impact of land data assimilation for a particular weather event rather than by globally averaging differences in skill. We argue that global skill assessments – while necessary – can mute the impact of land data assimilation, because the land’s influence on the atmosphere is constrained to certain locations and certain times. Instead, the event-based approach better highlights the true potential of land data assimilation in the context of NWP, especially for extreme events when accurate predictions are critical.

Jana Kolassa↗

DOE OSTI2022

Dominant controls of cold-season precipitation variability over the high mountains of Asia

A robust understanding of the sub-seasonal cold season (November–March) precipitation variability over the High Mountains of Asia (HMA) is lacking. Here, we identify dynamic and thermodynamic pathways through which natural modes of climate variability establish their teleconnections over the HMA. First, we identify evaporative sources that contribute to the cold season precipitation over the HMA and surrounding areas. The predominant moisture contribution comes from the mid-latitude regions, including the Mediterranean/Caspian Seas and Mediterranean land. Second, we establish that several tropical and extratropical forcings display a sub-seasonally fluctuating influence on precipitation distribution over the region during the cold season. Many of them varyingly interact, so their impacts cannot be explained independently or at seasonal timescales. Lastly, a single set of evaporative sources is not identifiable as the key determinant in propagating a remote teleconnection because the sources of moisture anomalies depend on the pattern of sub-seasonally varying dynamical forcing in the atmosphere.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Infrasonic directivity of monopole, dipole and bipole ground-surface reflected sources

Infrasound (acoustic waves below 20 Hz) can be used to detect, locate and quantify activity in the atmosphere such as volcanic eruptions and anthropogenic explosions. Attempts to quantify volcanic eruption parameters such as exit velocity, plume height and mass flow rate using infrasound data depend strongly on assumptions of the acoustic source type. Infrasonic sources may produce omnidirectional or directional wavefields, while propagation effects, such as interaction with topography, can induce further wavefield directivity that is measured by field instrumentation. Limited sampling of these wavefields can hinder our ability to infer the underlying source, and thus our understanding of the eruption characteristics. Equivalent sources are often used to represent acoustic source mechanisms and resultant wavefields. In this study, we review equivalent acoustic sources as they pertain to infrasonic scale and wavelengths commonly encountered in very local (⁠<5 km range) geophysical field deployments. We highlight the equivalent infrasonic bipole source that can be induced by ground-reflection of an elevated monopole; we are not aware of any prior infrasound studies that use the bipole source concept. We use analytical and numerical methods to explore source directivity of monopole, dipole and bipole ground-reflected sources at infrasonic frequencies as well as the additional directivity complications introduced by interactions with topography. We illustrate that for typical volcano-infrasound wavelengths, increasing height above the ground as well as increasing source frequency leads to increased wavefield directivity. Numerical modelling using a simple omnidirectional monopole source embedded in topography further illustrates that both horizontal and vertical infrasound directionality can be induced by topography at the distance scales appropriate for local volcano infrasound monitoring. Information summarized in this analytical and numerical exploration of infrasound directivity may be used to help guide future volcano-infrasound field deployments intended to estimate source parameters or quantify wavefield directivity. Analytic solutions for simple whole-space or half-space atmospheres provide useful formulations for planning or initially analysing geophysical field-scale experimental data; however, especially at very local distances from the source (⁠<5 km), 3-D simulations are necessary to account for complex topography commonly encountered in volcano-infrasound applications.

Infrasound↗

NASA NTRS1980

Sound attenuation by liners in a blown flap environment

Sound propagation through a hot wall-jet flow over an absorbing wall is studied. The radiated sound field subject to the influence of flow convection and refraction is evaluated, and the nature of acoustic attenuation attributable to a sound absorbing liner is determined. Using a two-dimensional model, the noise field under the aircraft is also determined, and a slug-flow model is used to describe the influence of flow, density, and temperature on acoustic sources in jets. Results show significant changes in the radiated source due to the interference phenomenon, and a good absorber has the potential of changing the sound pressure range of variation to unity. A liner is also found to increase or decrease sound pressure, depending on the frequency.

Parthasarathy, S. P.↗

NASA NTRS2012

Dust in a Type Ia Supernova Progenitor: Spitzer Spectroscopy of Kepler's Supernova Remnant

Characterization of the relatively poorly-understood progenitor systems of Type Ia supernovae is of great importance in astrophysics, particularly given the important cosmological role that these supernovae play. Kepler's Supernova Remnant, the result of a Type Ia supernova, shows evidence for an interaction with a dense circumstellar medium (CSM), suggesting a single-degenerate progenitor system. We present 7.5-38 micron IR spectra of the remnant, obtained with the Spitzer Space Telescope, dominated by emission from warm dust. Broad spectral features at 10 and 18 micron, consistent with various silicate particles, are seen throughout. These silicates were likely formed in the stellar outflow from the progenitor system during the AGB stage of evolution, and imply an oxygen-rich chemistry. In addition to silicate dust, a second component, possibly carbonaceous dust, is necessary to account for the short-wavelength IRS and IRAC data. This could imply a mixed chemistry in the atmosphere of the progenitor system. However, non-spherical metallic iron inclusions within silicate grains provide an alternative solution. Models of collisionally-heated dust emission from fast shocks (> 1000 km/s) propagating into the CSM can reproduce the majority of the emission associated with non-radiative filaments, where dust temperatures are approx 80-100 K, but fail to account for the highest temperatures detected, in excess of 150 K. We find that slower shocks (a few hundred km/s) into moderate density material (n(sub o) approx 50-100 / cubic cm) are the only viable source of heating for this hottest dust. We confirm the finding of an overall density gradient, with densities in the north being an order of magnitude greater than those in the south.

Williams, Brian J.↗

NASA NTRSDate not supplied

A Hinode-IRIS HOP for Searching for Possible Solar Jet Sources for Switchbacks Observed in the Near-Sun Solar Wind by the Parker Solar Probe (PSP)

The Parker Solar Probe (PSP) detects “switchback” magnetic structures in the near-Sun solar wind. A candidate for the source of these is solar coronal jets. Recent studies show that coronal jets are frequently made by eruptions of small-filament-carrying flux ropes, or mini filament flux-rope eruptions, that can inject a magnetic twist wave onto high-reaching ambient magnetic field(Sterling et al. 2015; Moore et al 2015).That is, many coronal jets evidently result from the eruption of a small-scale flux rope of twisted field that holds a mini filament of chromospheric-temperature plasma. The erupting flux rope reconnects with ambient far-reaching field, ejecting both reconnection-heated plasma and cool filament plasma (forming the multi-temperature jet spire), and transferring its twist onto the far-reaching field via a mechanism proposed by Shibata & Uchida(1986).If the far-reaching field opens into the heliosphere, it is plausible that the twist propagates out into the solar wind, possibly appearing as a switchback in the magnetic field in the solar wind (Sterling & Moore 2020).As part of the "Connecting Heliophysics Observatories and Models with PSP" (CHOMP)investigation led by Pete Riley and funded by NASA’s Heliophysics System Observatory Connect (HSOC) Program, we are currently We are currently running a joint Hinode-IRIS HOP, HOP 410, for coordinated observations of the Sun during times of PSP perihelia. An output from this study is a series of soft X-ray movies from Hinode/XRT. If PSP sees a strong switchback during any of its front-side or partly-front-side perihelia, we will examine the corresponding XRT movies to see whether exceptional jet activity occurred at the expected time and location of the potential source of the switchback on the Sun. This work is supported by the NASA HSOC and HGI programs, and by the MSFC Hinode project.

Alphonse C Sterling↗

NASA NTRS1995

Radial alignment simulation of solar wind streams observed by Pioneers 10 and 11 in 1974

A particularly favorable radial lineup between spacecraft in Earth orbit and Pioneers 10 and 11 (near the ecliptic plane at approximately 5.5 AU and approximately 4.5 AU, respectively) occurred in mid-1974, when the solar corona was in a stable and well-defined warped-dipole configuration. The radial alignment study reported here differs from previous applications of the technique in two respects: (1) It is the first time a two-dimensional (2-D) MHD model has been tested over such a lengthy propagation interval; the 2-D capability is crucial for treatment of the nonradial shearing motions occurring across the stream interface. (2) The three-dimensional (3-D) structure observed in the white light corona at that time is related to the systematic patterns of nonradial flow deflections appearing at the Pioneer corotating interaction region (CIR) fronts. Comparison of predicted and observed flows for pairs of streams in two successive rotations reveals that when the parent coronal hole projects far across the solar equator in a predominantly north-south orientation (i.e., nearest the 2-D idealization), the mapping is accurate down to details of the flow structures. But where the spacecraft tracks along a latitudinal boundary of a hole or the associated stream front is inclined at a shallow angle to the equator, the numerical projections deviate systematically from the observations. Among the sources of error are 3-D dynamical interactions neglected in the model, differential rotation effects, and slow temporal evolution of the coronal structures. A characteristic pattern of north-south and east-west deflections indicative of the 3-D geometry appears across the CIRs, but these patterns evidently reflect primarily the local, as opposed to global, orientation of the 3-D stream fronts. Such patterns appear common in CIRs observed by Pioneer during this period. These findings thus hold special relevance for the analysis of Ulysses observations, since the present coronal configuration is similar to that of 1974.

Pizzo, V. J.↗

DOE OSTI2025

Reference solutions for linear radiation transport: the Hohlraum and Lattice Benchmarks

Radiation transport describes the propagation of energetic particles through space as they interact with a surrounding material medium. In a kinetic description, radiation transport is modeled by a radiation transport equation (RTE) that prescribes the density of the radiation in position-momentum phase space. The purpose of this dataset is to provide highly resolved solutions to two benchmark problems. These two benchmarks do not possess exact solutions; moreover, the construction of a manufactured solution may require a non-physical source that is not desirable, especially if it spoils the physical nature of the solution. Thus the goal of this computational study is to provide a highly resolved reference solution for testing newer, more cost efficient methods that are currently being developed in the research community.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2020

The Heliospheric Current Sheet in the Inner Heliosphere Observed by the Parker Solar Probe

The Parker Solar Probe (PSP) completed its first solar encounter in 2018 November, bringing it closer to the Sun than any previous mission. This allowed in situ investigation of the heliospheric current sheet (HCS) inside the orbit of Venus. The Parker observations reveal a well defined magnetic sector structure placing the spacecraft in a negative polarity region for most of the encounter. The observed current sheet crossings are compared to the predictions of both potential field source surface and magnetohydrodynamic models. All the model predictions are in good qualitative agreement with the observed crossings of the HCS. The models also generally agree that the HCS was nearly parallel with the solar equator during the inbound leg of the encounter and more significantly inclined during the outbound portion. The current sheet crossings at PSP are also compared to similar measurements made by the Wind spacecraft near Earth at 1 au. After allowing for orbital geometry and propagation effects, a remarkable agreement has been found between the observations of these two spacecraft underlying the large-scale stability of the HCS. Finally, the detailed magnetic field and plasma structure of each crossing is analyzed. Marked differences were observed between PSP and Wind measurements in the type of structures found near the HCS. This suggests that significant evolution of these small solar wind structures takes place before they reach 1 au.

Szabo, Adam↗

NASA NTRS1986

The synchrotron-maser theory of type II solar radio emission processes - The physical model and generation mechanism

A theory is proposed to explain the generation mechanism of type II solar radio bursts. It is suggested that the shock wave formed at the leading edge of a coronal transient can accelerate electrons. Because of the nature of the acceleration process, the energized electrons can possess a 'hollow-beam' type distribution function. When the electron beam propagates along the ambient magnetic field to lower altitudes and attains larger pitch angles, a synchrotron-maser instability can set in. This instability leads to the amplification of unpolarized or weakly polarized radiation. The present discussion incorporates a model which describes the ambient magnetic field and background plasma by means of MHD simulation. The potential emission regions may be located approximately, according to the time-dependent MHD simulation. Since the average local plasma frequency in the source region can be evaluated from the MHD model, the frequent drift associated with the radiation may be estimated. The result seems to be in good agreement with that derived from observations.

Wu, C. S.↗

NASA NTRS1986

Cyclotron maser radiation from a source structure localized perpendicular to the ambient magnetic field

Particle simulations were used to investigate the properties of electromagnetic radiation produced by the cyclotron maser instability in a source region which is localized in the direction perpendicular to the uniform magnetic field. In the procedure used, the simulation model of Pritchett and Strangeway (1985) for the generation of auroral kilometric radiation (AKR) along an auroral field line was modified to include inhomogeneities perpendicular to the magnetic field. The results are compared with the previous homogeneous simulations with regard to efficiency and saturation mechanism and with expectations based on the feedback model. The implications for the generation and propagation of AKR are discussed.

Pritchett, P. L.↗

NASA NTRS1995

The 1991 October 24 flare: A challenge for standard models

The M9.8 solar flare of 1991 October 24 22:30 UT presents several interesting characteristics: (1) energy release starts high in the corona; (2) the primary chromospheric ribbons are initially well separated and do not move apart at an observable rate; (3) no evidence is found for an erupting filament or other driver. To explain this flare, we consider several canonical flare models, including a filament eruption, a confined filament eruption, current interruption, and interacting loops. We conclude that none of these scenarios unequivocally explains this flare. Two possibilities which cannot be ruled out are (1) the eruption of a filament unobservable in H-alpha which starts high in the corona and produces no ribbon motions smaller than our detection threshold and no perceptible expansion of the coronal X-ray source, and (2) energy release due to spontaneous, propagating reconnection which allows the system to essentially brighten in place.

Beaujardiere, J.-F. De LA↗

DOE OSTI2022

MODELING AND PARAMETRIC STUDY OF END-GAS AUTOIGNITION TO ALLOW THE REALIZATION OF ULTRA-LOW EMISSIONS, HIGH-EFFICIENCY HEAVY-DUTY SPARK-IGNITED NATURAL GAS ENGINES

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas (NG) engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by the low emission limits and the wide range of chemical reactivity in pipeline-quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its Methane Number (MN), which is analogous to the Octane Number for liquid SI fuels, can span from 65 to 95. In order to realize diesel-like efficiencies, SI NG engines must be designed to operate at high Brake Mean Effective Pressures (BMEP), near or beyond knock limits, over a wide range of fuel reactivity. This requires a deep understanding of the combustion-engine interactions pertaining to flame propagation and End-Gas Autoignition (EGAI), i.e., the autoignition of the unburned gas (end gas) ahead of the flame front. However, EGAI, if controlled, provides an opportunity to increase SI NG engine efficiency by increasing the combustion rate and the total fraction of burned fuel, mitigating the effects of the slow flame speeds characteristic of natural gas fuels, which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, to realize diesel-like efficiencies and ultra-low emissions on SI NG engines, this work proposes the study of the main parameters influencing the modeling and prediction of NG EGAI to allow for its control. In this work, a novel EGAI detection and onset determination method was developed to reliably quantify EGAI for data analysis and engine control. The new method allowed the prediction of EGAI on SI NG engines without the need to use engine- and operating-condition-dependent thresholds and reduced the error in quantifying the fraction of the total energy released by the EGAI event by up to 40%pts. One- and three-dimensional engine models were then developed to study the engine/fuel interactions that lead to NG EGAI and its performance benefits. These models, although having decent agreement with experimental data, showed the need to account for NOx chemistry when predicting NG EGAI due to a consistently later prediction of the EGAI onset (~1.65 crank-angle degrees) and thus, a new reduced chemical mechanism for real NG fuels was developed containing NOx chemistry. The new reduced mechanism improved the EGAI onset prediction agreement to within ±0.5 crank-angle degrees and decreased simulation time during combustion by nearly 50% when using the further reduced AREIS50NOx chemical mechanism. These models were then used to study the role of NG composition on EGAI, evaluate the engine/fuel interactions leading to NG EGAI, and perform engine optimization while leveraging EGAI to increase thermal efficiency. Piston design optimization combined with a Controlled EGAI (C-EGAI) combustion mode allowed a Heavy-Duty (HD) SI NG engine to operate at diesel-like efficiencies, i.e., Brake Thermal Efficiency (BTE) ≥44%. Experimental and modeling data analysis revealed that earlier and faster heat release increases combustion efficiency by an average of 1%pts, increases work transferred to the piston resulting in a decrease in exhaust losses by 50% depending on the engine operating condition while slightly increasing heat losses. Finally, the simulation results revealed an opportunity to further enhance the BTE (up to 50%) by enabling C-EGAI combustion at leaner conditions, λ=1.4-1.6.

Bestel, Diego Bernardi↗