Search NASA⌕ Search

SEARCH · Search NASA

Results for “mass flux”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Crossflow membrane filtration system for operando fouling characterization using transmission x-ray scattering

Membrane-based separations are widely used for wastewater treatment due to their low cost and efficiency. However, membrane fouling, which is the unwanted deposition or attachment of contaminants on membrane surfaces and/or within membrane pores, remains a major challenge as it increases the mass transfer resistance and reduces membrane productivity. Membrane fouling is typically probed by macroscopic performance metrics, such as flux decline, and ex situ characterization. However, this does not capture the membrane and fouling layer evolution under operating conditions, potentially masking important mechanisms and nonequilibrium pathways that impact fouling. Here, we present a remotely controlled crossflow membrane system and a custom membrane cell for operando fouling characterization using transmission small/wide angle x-ray scattering (SAXS/WAXS). This approach allows direct observation of the nanoscale changes occurring at the membrane surface during pressurized water treatment processes, enabling a new way to understand the connections between dynamic fouling behaviors and membrane performance. Nanoparticle fouling of porous membranes during ultrafiltration was investigated using operando SAXS, and mineral scaling of reverse osmosis membranes was investigated using operando WAXS. Furthermore, this system allows for tracking membrane fouling in real time and under realistic conditions, providing fundamental physical insights into how water chemistry and operating conditions affect macroscopic membrane performance. Moreover, this system opens the door for future in situ and operando studies, and it serves as a testbed for evaluating novel materials/processes for membrane-based separations.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Constraining Primordial Black Hole Abundance with the Galactic 511 keV Line

Models in which dark matter consists entirely of primordial black holes (PBHs) with masses around 10 17 g are currently unconstrained. However, if PBHs are a component of the Galactic dark matter density, they will inject a large flux of energetic particles into the Galaxy as they radiate. Positrons produced by these black holes will subsequently propagate throughout the Galaxy and annihilate, contributing to the Galactic 511 keV line. Here, using measurements of this line by the INTEGRAL satellite as a constraint on PBH positron injection, we place new limits on PBH abundance in the mass range 10 16 –10 17 g, ruling out models in which these PBHs constitute the entirety of dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Burn parameters affect PAH emissions at conditions relevant for prescribed fires

Wildfire smoke is a health hazard as it contains carcinogenic volatile compounds and fine particulate matter. In particular, exposure to polycyclic aromatic hydrocarbons (PAHs) is a major concern, since these compounds have been recognized as important contributors to the overall carcinogenic risk. In this work, gas and particle-phase PAH emissions from combustion of Eastern White Pine (Pinus strobus) were quantified using time-of-flight mass spectrometry over a range of burn conditions representative of wildfires and prescribed fires, including fuel moisture, heat flux, and oxygen concentration. We found that changing the burn environment lead to a variability of up to 77% in phenanthrene/anthracene emissions. This could explain a large part of the variability in PAH emission factors from biomass combustion reported in the literature. Here, we found that optimal conditions for fuel moisture content of 20–30%, sample heat load of 60 - 70 kW m -2 , and oxygen concentrations of 5–15% can significantly reduce the emissions of heavy molar weight PAHs. Our analysis showed that the relative carcinogenic risk from PAH exposure can be reduced by more than 50% under optimal conditions. In light of the increasing use of prescribed fire for forest management, the relationship between emissions and burn conditions that we have established provides a guidance for assessing the expected health impact from prescription burns, and can inform strategies to reduce PAH emissions from prescribed fire activities.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2021

Beneath the Bark: Assessing Woody Stem Water and Carbon Fluxes and Its Prevalence Across Climates and the Woody Plant Phylogeny

While woody stems are known to influence carbon and water dynamics, direct exchange with the atmosphere is seldom quantified, limiting our understanding of how these processes influence the exchange of mass and energy. The presence of woody stem chlorophyll in a diversity of climates and across a range of species suggests an evolutionary advantage to sustaining carbon assimilation and water relations through permeable stem tissue. However, no formal evaluation of this hypothesis has been performed. In this mini-review, we explore the interactions between woody stems and the atmosphere by examining woody stem photosynthesis and bark-atmosphere water exchange. Specifically, we address the following questions: (1) How do water and carbon move between the atmosphere and woody stems? (2) In what climate space is woody stem photosynthesis and bark water uptake advantageous? (3) How ubiquitous across plant families is woody stem photosynthesis and bark-atmosphere water exchange? In the literature, only seven species have been identified as exhibiting bark water uptake while over 300 species are thought to conduct woody stem photosynthesis. The carbon dioxide and water gained from these processes can offset respiration costs and improve plant water balance. These species span diverse biomes suggesting a broad prevalence of bark-atmosphere permeability. Finally, our results demonstrate that there may be an evolutionary component as demonstrated by a high Pagel’s lambda for the presence of stem photosynthesis. We end with recommendations for future research that explores how bark water and carbon interactions may impact plant function and mass flow in a changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tuning plasma parameters to control reactive species fluxes to substrates in the context of plasma catalysis

The key reactive plasma-produced gas phase species responsible for the enhanced conversion of chemicals in plasma catalysis compared to thermal catalysis have to date not been identified. This outstanding question is mainly due to the inherent large variety of plasma-produced species and the challenge of controlling and measuring the flux of each constituent of the cocktail of reactive species to a (catalytic) substrate. In this paper, we explore the possibility to control the dominant reactive species fluxes, relevant for plasma-catalysis, to a substrate in the effluent of an RF driven Ar-O 2 plasma jet. The absolute species densities of the major reactive species (O, O 2 (a 1 Δ g ), O 3 and ions) were quantified by molecular beam mass spectrometry (MBMS) to assess the possibility of using treatment distance, O 2 admixture concentration, plasma dissipated power, RF modulation frequency and duty cycle as well as the feed gas flow rate to alter the dominant species densities. Selected experimental results were also compared with a pseudo-1D plug flow model. The short-lived and long-lived species can be effectively separated by changing the treatment distance and the RF modulation frequency. Furthermore, adjusting the O 2 admixture concentration enables to change the ratio of the O 2 (a 1 Δ g ) and O 3 density. The changes in the trend of ion and O flux were found to be very similar for nearly all investigated parameters. Nonetheless the gas flow rate was able to significantly change the ratio of the O and ion density in the plasma jet effluent. Here, the impact of the surface-dependent loss probability and boundary layer reactions on the species flux to a substrate and how this qualitatively relates to the MBMS density measurements is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Film condensation with high heat fluxes and scaled experiments using pure steam for reactor containment cooling

Condensation tests were performed using a newly developed test facility for scaling the passive containment cooling system (PCCS) to a small modular reactor (SMR). The PCCS of the SMR plays a pivotal role in ensuring greater safety, reliability, and compactness than what is afforded by traditional reactors. Therefore, a well-designed PCCS is essential to SMRs. However, previous studies and test data were unsuitable for scaling, due to high variation in the test geometry and operating conditions. This study intends to close this research gap by using a novel designed scaled test facility consisting of vertical condensing test sections featuring 1-, 2-, and 4-inch-diameter condensing tubes with annular water cooling, and by applying superheated and saturated steam with different steam mass flow ranges of 5–25 g/s. Further, the primary test data, including axial temperatures, mass flow rates, and pressures, were used in conjunction with a standard data reduction method to estimate critical parameters such as heat fluxes, heat transfer coefficients, and condensation rates. These scaled test data would support improving empirical correlations and validating condensation models to identify scaling distortion for SMR PCCSs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

FLASH pilot survey: detections of associated 21 cm H i absorption in GAMA galaxies at 0.42 < z < 1.00

ABSTRACT We present the results of a search for associated 21 cm H i absorption at redshift 0.42 < z < 1.00 in radio-loud galaxies from three Galaxy And Mass Assembly (GAMA) survey fields. These observations were carried out as part of a pilot survey for the Australian SKA Pathfinder (ASKAP) First Large Absorption Survey in H i (FLASH). From a sample of 326 radio sources with 855.5 MHz peak flux density above 10 mJy, we detected two associated H i absorption systems, in SDSS J090331+010847 at z = 0.522 and SDSS J113622+004852 at z = 0.563. Both galaxies are massive (stellar mass $\gt 10^{11}\, \mathrm{M}_\odot$) and have optical spectra characteristic of luminous red galaxies, though spectral energy distribution fitting implies that SDSS J113622+004852 contains a dust-obscured starburst with star formation rate ∼69 M⊙ yr−1. The H i absorption lines have a high optical depth, with τpk of 1.77 ± 0.16 for SDSS J090331+010847 (the highest value for any z > 0.1 associated system found to date) and 0.14 ± 0.01 for SDSS J113622+004852. In the redshift range probed by our ASKAP observations, the detection rate for associated H i absorption lines (with τpk > 0.1 and at least 3σ significance) is $2.9_{-2.6}^{+9.7}$ per cent. Although the current sample is small, this rate is consistent with a trend seen in other studies for a lower detection rate of associated 21 cm H i absorption systems at higher redshift. We also searched for OH absorption lines at 0.67 < z < 1.34, but no detection was made in the 145 radio sources searched.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Effect of Formation Heterogeneity on CO 2 Dissolution in Subsurface Porous Media

Dissolution trapping is one of the most dominant mechanisms for the secure storage of CO 2 injected in porous subsurface formations saturated with brine. This trapping mechanism is enhanced by convective mixing, which occurs due to the gravitational instability between the different fluid layers in the aquifer. The reservoir permeability also plays a crucial role in the dissolution rate and overall fluid flow dynamics during the density-driven convection in porous media. Herein this study investigates the role of complex heterogeneity, i.e., irregular permeability distribution in CO 2 dissolution, using a novel experimental approach to create medium permeability heterogeneity in Hele–Shaw cells. Complex subsurface transport phenomena such as a preferential dissolution path, CO 2 sweep efficiency, changes in finger morphology, and CO 2 concentration distribution are visualized by creating heterogeneous media. Experimental results showed that reservoir permeability heterogeneity causes significant channeling effects and poor sweep efficiency. A scaling relationship between average finger growth rate (Gr) and permeability (k) was obtained as Gr [m s –1 ] = 266.8k [m 2 ] + 1.20 × 10 –6 . Furthermore, the mass of CO 2 dissolved is calculated using the spectrophotometric method to characterize the convective instability. The convective flux was analyzed by comparing the experimental dissolution flux with the theoretical diffusion flux, calculating a maximum Sherwood number of 6.8. The study’s findings improve the current understanding of the CO 2 convection morphology in heterogeneous media, allowing better assessment of long-term CO 2 storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Electronics integration for the GE2/1 and ME0 GEM detector systems for the CMS phase-2 muon system upgrade

Abstract The Large Hadron Collider is currently undergoing its High Luminosity upgrade, which is set to increase the instantaneous luminosity by about a factor of five. Consequently, the Compact Muon Solenoid experiment is upgrading its muon spectrometer to cope with the increased muon flux in the forward region. The GE2/1 triple-gas electron multiplier detector, which has recently entered the mass production phase, and the ME0 triple-GEM detector system, which is in the late prototyping phase, are undergoing electronics integration. These proceedings briefly discuss the frontend electronics for the GE2/1 and ME0 detector systems, the electronics integration testing process, and the future plans for the frontend electronics of these two detector systems by the CMS GEM collaboration.

Instruments & Instrumentation↗
DOE OSTI2020

Galaxy and Mass Assembly: A Comparison between Galaxy–Galaxy Lens Searches in KiDS/GAMA

Strong gravitational lenses are a rare and instructive type of astronomical object. Identification has long relied on serendipity, but different strategies—such as mixed spectroscopy of multiple galaxies along the line of sight, machine-learning algorithms, and citizen science—have been employed to identify these objects as new imaging surveys become available. We report on the comparison between spectroscopic, machine-learning, and citizen-science identification of galaxy–galaxy lens candidates from independently constructed lens catalogs in the common survey area of the equatorial fields of the Galaxy and Mass Assembly survey. In these, we have the opportunity to compare high completeness spectroscopic identifications against high-fidelity imaging from the Kilo Degree Survey used for both machine-learning and citizen-science lens searches. We find that the three methods—spectroscopy, machine learning, and citizen science—identify 47, 47, and 13 candidates, respectively, in the 180 square degrees surveyed. These identifications barely overlap, with only two identified by both citizen science and machine learning. We have traced this discrepancy to inherent differences in the selection functions of each of the three methods, either within their parent samples (i.e., citizen science focuses on low redshift) or inherent to the method (i.e., machine learning is limited by its training sample and prefers well-separated features, while spectroscopy requires sufficient flux from lensed features to lie within the fiber). These differences manifest as separate samples in estimated Einstein radius, lens stellar mass, and lens redshift. The combined sample implies a lens candidate sky density of ∼0.59 deg{sup −2} and can inform the construction of a training set spanning a wider mass–redshift space. A combined approach and refinement of automated searches would result in a more complete sample of galaxy–galaxy lens candidates for future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Massive neutrinos and degeneracies in Lyman-alpha forest simulations

Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-α (Lyα) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppression is a key method for inferring neutrino masses from cosmological data, and is one of the main goals of ongoing and future surveys like eBOSS, DES, LSST, Euclid or DESI . The clustering in the Lyα forest traces the quasi-linear power at late times and on small scales. In combination with observations of the cosmic microwave background, the forest therefore provides some of the tightest constraints on the sum of the neutrino masses. However there is a well-known degeneracy between Σmν and the amplitude of perturbations in the linear matter power spectrum. We study the corresponding degeneracy in the 1-D flux power spectrum of the Lyα forest, and for the first time also study this degeneracy in the 3-D flux power spectrum. We show that the non-linear effects of massive neutrinos on the Lyα forest, beyond the effect of linear power amplitude suppression, are negligible, and this degeneracy persists in the Lyα forest observables to a high precision. Here, we discuss the implications of this degeneracy for choosing parametrisations of the Lyα forest for cosmological analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Plasma Dynamics and Nonthermal Particle Acceleration in 3D Nonrelativistic Magnetic Reconnection

Understanding plasma dynamics and nonthermal particle acceleration in 3D magnetic reconnection has been a long-standing challenge. In this paper, we explore these problems by performing large-scale fully kinetic simulations of multi-X-line plasmoid reconnection with various parameters in both the weak- and strong-guide-field regimes. In each regime, we have identified its unique 3D dynamics that lead to field-line chaos and efficient acceleration, and we have achieved nonthermal acceleration of both electrons and protons into power-law spectra. The spectral indices agree well with a simple Fermi acceleration theory that includes guide-field dependence. In the low-guide-field regime, the flux rope kink instability governs the 3D dynamics for efficient acceleration. The weak dependence of the spectra on the ion-to-electron mass ratio and β ($\ll$1) implies that the particles are sufficiently magnetized for Fermi acceleration in our simulations. While both electrons and protons are injected at reconnection exhausts, protons are primarily injected by perpendicular electric fields through Fermi reflections and electrons are injected by a combination of perpendicular and parallel electric fields. The magnetic power spectra agree with in situ magnetotail observations, and the spectral index may reflect a reconnection-driven size distribution of plasmoids instead of the Goldreich–Sridhar vortex cascade. As the guide field becomes stronger, the oblique flux ropes of large sizes capture the main 3D dynamics for efficient acceleration. Intriguingly, the oblique flux ropes can also experience flux rope kink instability, to drive extra 3D dynamics. This work has broad implications for 3D reconnection dynamics and particle acceleration in heliophysics and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: The TRAPPIST-1 Case Study

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in X-ray and ultraviolet (XUV) irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on an ultraviolet continuum flux measurement (centered around ∼1900 Å) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photodissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 × 10{sup −4} L {sub *}, while the full flux including non-flaring periods is 6.1 × 10{sup −5} L {sub *}, when L {sub *} is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass loss and tidal evolution and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio Q{sup ′}=3Q/2k{sub 2} for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator VPLanet to study ocean retention for these planets using our derived system parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Surface Water Mass Transformation in the Southern Ocean: The Role of Eddies Revisited

The water mass transformation (WMT) framework describes how water of one class, such as a discrete interval of density, is converted into another class via air–sea fluxes or interior mixing processes. Here, this paper investigates how this process is modified at the surface when mesoscale ocean eddies are present, using a state-of-the-art high-resolution climate model with reasonable fidelity in the Southern Ocean. The method employed is to coarse-grain the high-resolution model fields to remove eddy signatures, and compare the results with those from the full model fields. This method shows that eddies reduced the WMT by 2–4 Sv (10%–20%; 1 Sv ≡ 10 6 m 3 s −1 ) over a wide range of densities, from typical values of 20 Sv in the smoothed case. The corresponding water mass formation was reduced by 40% at one particular density increment, namely, between 1026.4 and 1026.5 kg m −3 , which corresponds to the lighter end of the range of Indian Ocean Mode Water in the model. The effect of eddies on surface WMT is decomposed into three terms: direct modulation of the density outcrops, then indirectly, by modifying the air–sea density flux, and the combined effect of the two, akin to a covariance. It is found that the first and third terms dominate, i.e., smoothing the outcrops alone has a significant effect, as does the combination of smoothing both outcrops and density flux distributions, but smoothing density flux fields alone has little effect. Results from the coarse-graining method are compared to an alternative approach of temporally averaging the data. Implications for climate model resolution are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Atacama Cosmology Telescope: A Search for Planet 9

We use Atacama Cosmology Telescope (ACT) observations at 98 GHz (2015–2019), 150 GHz (2013–2019), and 229 GHz (2017–2019) to perform a blind shift-and-stack search for Planet 9. The search explores distances from 300 au to 2000 au and velocities up to 6farcm3 per year, depending on the distance (r). For a 5 Earth-mass Planet 9 the detection limit varies from 325 au to 625 au, depending on the sky location. For a 10 Earth-mass planet the corresponding range is 425 au to 775 au. The predicted aphelion and most likely location of the planet corresponds to the shallower end of these ranges. The search covers the whole 18,000 square degrees of the ACT survey. No significant detections are found, which is used to place limits on the millimeter-wave flux density of Planet 9 over much of its orbit. Overall we eliminate roughly 17% and 9% of the parameter space for a 5 and 10 Earth-mass Planet 9, respectively. These bounds approach those of a recent INPOP19a ephemeris-based analysis, but do not exceed it. We also provide a list of the 10 strongest candidates from the search for possible follow-up. More generally, we exclude (at 95% confidence) the presence of an unknown solar system object within our survey area brighter than 4–12 mJy (depending on position) at 150 GHz with current distance 300 au < r < 600 au and heliocentric angular velocityy <?CDATA $1\buildrel{\,\prime}\over{.} 5\,{\mathrm{yr}}^{-1}\lt v\cdot \tfrac{500\,\mathrm{au}}{r}\lt 2\buildrel{\prime\prime}\over{.} 3\,{\mathrm{yr}}^{-1}$?> 1 .' 5 yr – 1 < v · 500 au r < 2 ." 3 yr – 1 , corresponding to low-to-moderate eccentricities. These limits worsen gradually beyond 600 au, reaching 5–15 mJy by 1500 au.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Optimization-based, property-preserving finite element methods for scalar advection equations and their connection to Algebraic Flux Correction

In this paper, we continue our efforts to exploit optimization and control ideas as a common foundation for the development of property-preserving numerical methods. Here we focus on a class of scalar advection equations whose solutions have fixed mass in a given Eulerian region and constant bounds in any Lagrangian volume. Our approach separates discretization of the equations from the preservation of their solution properties by treating the latter as optimization constraints. This relieves the discretization process from having to comply with additional restrictions and makes stability and accuracy the sole considerations in its design. A property-preserving solution is then sought as a state that minimizes the distance to an optimally accurate but not property-preserving target solution computed by the scheme, subject to constraints enforcing discrete proxies of the desired properties. Furthermore, we consider two such formulations in which the optimization variables are given by the nodal solution values and suitably defined nodal fluxes, respectively. A key result of the paper reveals that a standard Algebraic Flux Correction (AFC) scheme is a modified version of the second formulation obtained by shrinking its feasible set to a hypercube. In conclusion, we present numerical studies illustrating the optimization-based formulations and comparing them with AFC

97 MATHEMATICS AND COMPUTING↗