Search NASA⌕ Search

SEARCH · Search NASA

Results for “FUNCTIONAL INTEGRATION”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Cosmology Large Angular Scale Surveyor (CLASS): 90 GHz Telescope Pointing, Beam Profile, Window Function, and Polarization Performance

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background (CMB) over ∼75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the large angular scale CMB polarization to constrain the tensor-to-scalar ratio and the optical depth to last scattering. This paper presents the optical characterization of the 90 GHz telescope. Observations of the Moon establish the pointing while dedicated observations of Jupiter are used for beam calibration. The standard deviations of the pointing error in azimuth, elevation, and boresight angle are 1$^{'}_.$3, 2$^{'}_.$1, and 2$^{'}_.$0, respectively, over the first 3 yr of observations. This corresponds to a pointing uncertainty ∼7% of the beam’s full width at half-maximum (FWHM). The effective azimuthally symmetrized instrument 1D beam estimated at 90 GHz has an FWHM of 0$^°_.$620 ± 0$^°_.$003 and a solid angle of 138.7 ± 0.6(stats.) ± 1.1(sys.) μsr integrated to a radius of 4°. The corresponding beam window function drops to $b^2_ℓ$ = 0.93, 0.71, 0.14 at ℓ = 30, 100, 300, respectively. Far-sidelobes are studied using detector-centered intensity maps of the Moon and measured to be at a level of 10 −3 or below relative to the peak. The polarization angle of Tau A estimated from preliminary survey maps is 149°.6 ± 0°.2(stats.) in equatorial coordinates. The instrumental temperature-to-polarization (T → P) leakage fraction, inferred from per-detector demodulated Jupiter scan data, has a monopole component at the level of 1.7 × 10 −3 , a dipole component with an amplitude of 4.3 × 10 −3 , with no evidence of quadrupolar leakage.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Selective Reduction of Carbon Dioxide in Water Using [M(bpy2+)(CO)3(I)]2+ (M = Mn, Re) Electrocatalysts with Pendent Cations

Manganese(I) carbonyl complexes are promising electrocatalysts for CO2 reduction, yet their application in homogeneous aqueous media remains limited by poor solubility and selectivity. Here, we report water-soluble Mn(I) and Re(I) complexes fac-[M(bpy2+)(CO)3X]2+ (X = I or Cl), featuring bipyridine ligands functionalized with -Ph-CH2-(NMe3)+ cationic ammonium groups that integrate water solubility with secondary-sphere stabilization. In bicarbonate buffer at pH 6.8, the Mn catalyst is completely selective for CO production at a low overpotential (η = 0.3 V), operating by a protonation-first mechanism with observed rates of ~10 s−1. Pulse radiolysis reveals that the one-electron reduced Mn species undergoes dimerization in the absence of CO2 but uniquely reacts competitively with CO2 through an initial pre-equilibrium followed by fast formation of a dinuclear CO2-bridged species (ΔGo = −12.4 kcal mol−1). At a higher 0.6 V over-potential, a faster reduction-first pathway (~100 s−1) is available upon reduction of the metallo-carboxylic acid intermediate, Mn-CO2H2+; however, this regime is functionally limited by the formation of a resistive, noncatalytic film on the electrode surface. Comparison to the analo-gous water-soluble Re catalyst (kobs = 440 s−1, η = 0.6 V) highlights the distinct mechanistic ad-vantages of earth-abundant Mn in low-potential catalysis. These results demonstrate how cati-onic second-sphere design enables selective, homogeneous CO2 reduction in water while reveal-ing competing radical and electrode-mediated processes that govern catalytic performance.

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

Functional stimuli-responsive polymers on micro- and nano-patterned interfaces

Micro- and nano-patterned surfaces offer precise control over morphology and chemical composition, enhancing the stability, durability, and functionality of coating materials. When combined with stimuli-responsive polymers, these surfaces gain dynamic adaptability, enabling reversible binding, reusable sensing, and selective molecular capture. Furthermore, while recent review articles have explored various aspects of stimuli-responsive materials, from hydrogel patterns for bioanalytical applications to shape-morphing hydrogels for soft robotics and sensors, a comprehensive review focused on the integration of smart polymers with micro- or nano-patterned interfaces remains absent. This review addresses key surface patterning techniques, including soft lithography, colloidal lithography, and polymer brush photolithography, as well as advances in surface-initiated polymerization methods, such as surface-initiated controlled radical polymerization (SI-CRP). In addition, we discuss recent progress in integrating stimuli-responsive polymers with patterned surfaces to create advanced, functional materials.

Colloidal lithography↗
DOE OSTI2025

Fast HARDI Uncertainty Quantification and Visualization with Spherical Sampling

In this paper, we study uncertainty quantification and visualization of orientation distribution functions (ODF), which corresponds to the diffusion profile of high angular resolution diffusion imaging (HARDI) data. The shape inclusion probability (SIP) function is the state‐of‐the‐art method for capturing the uncertainty of ODF ensembles. The current method of computing the SIP function with a volumetric basis exhibits high computational and memory costs, which can be a bottleneck to integrating uncertainty into HARDI visualization techniques and tools. We propose a novel spherical sampling framework for faster computation of the SIP function with lower memory usage and increased accuracy. In particular, we propose direct extraction of SIP isosurfaces, which represent confidence intervals indicating spatial uncertainty of HARDI glyphs, by performing spherical sampling of ODFs. Our spherical sampling approach requires much less sampling than the state‐of‐the‐art volume sampling method, thus providing significantly enhanced performance, scalability, and the ability to perform implicit ray tracing. Our experiments demonstrate that the SIP isosurfaces extracted with our spherical sampling approach can achieve up to 8164× speedup, 37282× memory reduction, and 50.2% less SIP isosurface error compared to the classical volume sampling approach. We demonstrate the efficacy of our methods through experiments on synthetic and human‐brain HARDI datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High‐Temperature Single‐Photon Emission From Covalently Functionalized van der Waals Heterostructures

Two-dimensional (2D) transition metal dichalcogenides (TMDs) such as tungsten diselenide (WSe 2 ) are attractive nanomaterials for quantum information applications due to single-photon emission (SPE) from intrinsic atomic defects. Defect and strain engineering techniques have been developed to produce high purity, deterministically placed SPE in WSe 2 . However, a major challenge in the application of these techniques is the low temperature required to observe defect-bound TMD exciton emission, typically limiting SPE to T < 30 K. SPE at higher temperatures either loses purity or requires integration into complex devices such as optical cavities. Here, 2D heterostructure engineering and molecular functionalization are combined to achieve high purity (>90%) SPE in strained WSe 2 persisting to over T = 90 K. Covalent diazonium functionalization of graphite in layered WSe 2 /graphite heterostructures maintains high purity up to T = 90 K and single-photon source integrity up to T = 115 K. This method preserves the best qualities of SPE from WSe 2 while increasing working temperature to more than three times the typical range. This work demonstrates the versatility of surface functionalization and heterostructure design to synergistically improve the properties of quantum emission and offers new insights into the phenomenon of SPE from 2D materials.

2D materials↗
DOE OSTI2024

Domain decomposition efforts in Chroma

Multilevel integration is a technique used in lattice quantum chromodynamics (LQCD) simulations, where the correlation functions are decomposed into factors that depend only on fields localized within specific lattice subdomains. The approach significantly improves the efficiency of the calculations because the factors can be independently integrated. The main challenge in implementing multilevel integration is that the standard formulation of the LQCD action is not local in the presence of fermions. However, past studies have shown that approximations of the quark propagator and the fermionic determinant, which depend on the gauge fields within specific subdomains, can still be effective. We will present the current progress in supporting domain decomposition within the Chroma software, which is necessary for implementing multilevel integration approaches. Most of the efforts have focused on extending the Chroma code base to efficiently manipulate subdomains for lattice fields and operators, as well as performing inversions and eigendecompositions within these domains.

Alcalde, Eloy Romero↗
DOE OSTI2024

Two-loop master integrals for leading-color $$ pp\to t\overline{t}H $$ amplitudes with a light-quark loop

Abstract We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in$$ t\overline{t}H $$ t t ¯ H production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations inϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Physics↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation

This Final Scientific and Technical Report summarizes work performed under the Phase IIA SBIR project “Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation” (DE-SC0020906) from August 2023 through August 2025. The objective of the Phase IIA effort was to mature and harden capabilities developed during Phase II, with the goal of enabling practical interoperability between Coreform’s isogeometric analysis (IGA) technologies and the Multiphysics Object-Oriented Simulation Environment (MOOSE), while improving robustness, performance, and scalability for complex, nuclear-relevant geometries. Over the course of Phase IIA, the project established and validated an extraction-based interoperability pathway between Coreform tools and MOOSE. A combined mesh and matrix format was defined collaboratively with MOOSE developers and integrated into the solver, enabling standard MOOSE workflows to operate on data exported from Coreform’s IGA and Flex Representation Method (FRM) pipelines. Early demonstrations validated architectural compatibility using linear solid mechanics problems, while later efforts focused on benchmark testing and external use. By the end of the project period, engineers at BWXT were able to independently set up and execute a simulation using the Coreform–MOOSE workflow and provide direct feedback that informed further refinement. In parallel, substantial effort was devoted to improving the robustness of trimmed U-spline construction for complex CAD geometries. A growing test suite of nuclear-relevant models was compiled through collaboration with multiple stakeholders and used to drive extensive bug fixing and reliability improvements. These efforts resulted in improved robustness and performance, including the addition of fallback capabilities that enhance reliability when the underlying commercial CAD kernel fails. Performance-oriented work progressed later in the project, with the development and demonstration of methods to decompose complex geometries into structured subregions and updated data representations to support more efficient solver processing. Additionally, extensive enhancements to threadsafe parallel data structures and trimming operations established a foundation for scalable processing of large assemblies. Collaboration with Sandia National Laboratories on the SGM geometric modeling kernel advanced to a functioning interface test case, positioning the workflow for future kernel integration. Overall, the Phase IIA effort successfully transitioned the project from architectural proof-of-concept to externally exercised, solver-integrated capability, while clarifying remaining technical challenges related to standardization, performance optimization, and kernel integration.

42 ENGINEERING↗
DOE OSTI2025

Non-Destructive, Three-Dimensional Imaging of Processes in the Rhizosphere Utilizing High Energy Photons

Soil structure, which can be described as the aggregation and distribution of pore spaces, regulates carbon, nutrient, and water cycling across the Earth system. Yet the inability to make quantitative, dynamic, in situ measurements of soil structure and rhizosphere carbon flow has prevented meaningful incorporation of soil structural processes into Earth System Models (ESMs). The key limitations are: (i) Lack of scale integration between micron-scale soil structure and ecosystem-scale models, (ii) Poor functional linkage between soil structural properties and biogeochemical processes, and (iii) Absence of dynamic 3D measurements of rhizosphere structural changes and carbon transformations. To address these challenges, we developed a new integrated positron emission tomography (PET) - microcomputed tomography (CT) imaging platform that enables the first 4D (spatial 3D and time), non-invasive, quantitative imaging of carbon allocation and rhizosphere structural dynamics in living plants and intact soils. The system combines: • Rhizo-PET (R-PET): a high-resolution positron emission tomography scanner optimized for 11 CO 2 tracing in plant roots. • a-Se micro-CT: a high-contrast, high-spatial resolution CT system based on amorphous selenium (a-Se) direct conversion technology, enabling micron-scale visualization of soil structure and root–soil interfaces. Together, these advances allow us to quantify how carbon exudates move, transform, and stabilize within the rhizosphere, directly informing missing processes in BER-relevant carbon cycle models.

42 ENGINEERING↗
DOE OSTI2026

Test of Harmonic Coil Accuracy Using a Multi Vertex Probe

Magnets are integral components of particle accelerators, with different multipoles used to perform various beam conditioning functions, like steering and focusing. Testing accelerator magnets to ensure compliance with specification is a crucial step of fabrication. One of the most popular techniques for verifying magnetic field strength and integrity is the harmonic coil probe. Harmonic coil probes are comprised of passive loops that are rotated around the aperture of a magnet and are used to measure magnetic flux. Oftentimes the measurements of the harmonic coil probe are affected by transverse vibrations along the axis of rotation and torsional vibrations of the position closure sensor, both causing spurious harmonics. To make matters challenging, it is difficult to distinguish between the real harmonics of the magnet and the spurious harmonics during analysis. This research explores the use of a novel configuration of harmonic coil probe, Multi-Vertex Probe (MVP), to address multiple factors that affect the accuracy and repeatability of harmonic coil probe measurement. Most of the results are based on simulation data, in which artificial torsional distortions are introduced and removed with the proposed processing method. In addition to implementing torsional vibration removal in simulation, laboratory analysis regarding the source of vibrations was performed using a pre-fabricated MVP. More specifically, various mechanical conditions were used to explore the standard deviation of main harmonic’s flux amplitude. The conditions included: coupling type (rigid or soft), rotation speed, probe starting position within magnetic field, and cabling location inside the magnetic field. Along with this, the MVP’s measured noise floor in and outside a magnetic field was quantified.

Lofquist, Claire [Northern Illinois U.]↗
DOE OSTI2025

PNNL-Predictive-Phenomics/ProteoMeter

ProteoMeter is a Python package that assists in the statistical analysis of global proteomics, protein post-translation modification (PTM), and limited proteolysis (LiP) data. It contains batch correction, normalization, and statistical testing methods, as well as functions that "roll up" peptide-level data to the single-site level. It has a robust user configuration system, allowing it to flexibly integrate different types of experiment designs. For basic usage, a simple configuration file provides the essential functionality. Advanced users have access to the entire statistical pipeline for fine-tuning analyses. Processed data is easily exported to many common spreadsheet and data-frame formats.

Rozum, Jordan [Pacific Northwest National Lab]↗
DOE OSTI2025

Advancing 3D surface imaging: single-axis structured light illumination plenoptic camera with machine learning integration

Structured light illumination (SLI) is a configurable 3D surface imaging modality that can function largely independently of surface texture. At the same time, machine learning (ML) approaches are providing new ways to capture relevant information from SLI patterns, avoiding the need to develop advanced computer vision algorithms. By projecting an optical pattern onto a surface and measuring the apparent distortion of that pattern, one can determine surface topography from a single image. Common realizations of SLI 3D imaging use off-axis SLI to allow for parallax-based determination of depth; however, in constrained geometries, the ability to make single-axis measurements can be of major benefit. While plenoptic imaging (PI) cameras have long been developed for the purpose of single-axis 3D imaging, they are generally reliant on the surface texture of the measured object, thus making them unreliable in certain experimental conditions. Therefore, we present a single-axis 3D SLI plenoptic camera, which combines the single-axis benefits of PI technology while using coaxial SLI to maintain indifference to surface conditions. We also present a study of the camera capabilities paired with the development of several algorithms, including traditional feature tracking methods as well as ML methods, which are found to enhance resolution and range. We report depth sensitivity down to 0.2% $\frac{dz}{z_0}$. The single-axis SLI 3D plenoptic camera demonstrates potential applicability for in-situ topographical measurements under a wide range of conditions including, but not limited to, objects without trackable surface texture, high temperatures, and constrained geometry environments.

Imaging systems↗
DOE OSTI2025

Immersive Scientific Visualization of Molten-Salt Reactor Waste Characteristics Using Virtual Reality

Immersive visualization is changing how we explore, communicate, and understand complex scientific systems. In nuclear energy, an area in which data are often multidimensional, time-dependent, and difficult to interpret, virtual reality (VR) represents a powerful and intuitive informational medium. This work introduces a VR-based platform that visualizes the post-shutdown behavior and waste management lifecycle of molten-salt reactors (MSRs), a next-generation reactor type with unique operational and safety characteristics. The platform, built in Unity, is streamed on the Meta Quest 3 headset. It transforms high-fidelity simulation data into an interactive, immersive experience. Users can explore time-dependent reactor characteristics such as nuclide decay, which is a key factor for evaluating reactor waste strategies. The datasets were generated using the MOOSE (Multiphysics Object-Oriented Simulation Environment) framework and then processed through ParaView scripting for smooth integration into Unity. From a visualization standpoint, the platform emphasizes spatial storytelling, temporal exploration, and user-centered interaction. Users can navigate 3D reactor geometries, slice through volumetric data, and manipulate time to observe how physical phenomena evolve. Real-scale rendering and embodied interaction make the experience feel tangible. The interface is designed to be accessible, even to those without nuclear or simulation expertise. This lowers the barrier for stakeholders, policymakers, and the general public, while still supporting expert analysis and collaborative decision-making. This work shows how immersive visualization can function as both a scientific tool and a communication interface. By integrating simulation, processing, and visualization into a cohesive workflow, we offer a scalable framework for immersive scientific storytelling. The modular design supports future extensions to other reactor types and lifecycle stages, from shutdown to long-term storage, making the platform adaptable for both research and outreach.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Surface water and groundwater FTICR-MS, NPOC, and TN from nine wetlands and three upland wells at the Tanglewood Biological Station, Alabama

This dataset supports a broader study examining wetland hydrobiogeochemical responses to flood disturbance and the subsequent impacts on watershed nutrient export. The study was designed following ICON (integrated, coordinated, open, and networked) principles. Samples were collected from nine wetlands and three upland wells at the Tanglewood Biological Station, Alabama in August 2024 and February 2025, during the dry and wet season, respectively. The contents include geochemistry (dissolved organic carbon measured as non-purgeable organic carbon; total dissolved nitrogen) and organic matter characterization (FTICR-MS). Related water level data from the same locations can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253. Additional geochemistry will be published in a separate data package. For details on how to navigate this data package, see this infographic from the River Corridor SFA https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions.This dataset is comprised of (1) a folder containing environmental context photos; (2) file-level metadata; (3) data dictionary; (4) field metadata; (5) readme; (6) international generic sample number (IGSN) mapping file; (7) the field protocol; and (8) a subfolder with sample data. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (2) total nitrogen data and averages; (3) methods codes; and (4) a subfolder of 12 Tesla (12T) Fourier-transform ion cyclotron resonance mass spectrometry (FTICR-MS) data. All files are .csv, .pdf, .jpeg, or .jpg.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

PAC in DESI. I. Galaxy stellar mass function into the 10 6 M ⊙ frontier

The Photometric objects Around Cosmic webs (PAC) method integrates cosmological photometric and spectroscopic surveys, offering valuable insights into galaxy formation. PAC measures the excess surface density of photometric objects, $\bar{n}_2w_{{\rm {p}}}$, with specific physical properties around spectroscopic tracers. In this study, we improve the PAC method to make it more rigorous and eliminate the need for redshift bins. We apply the enhanced PAC method to the DESI Y1 BGS Bright spectroscopic sample and the deep Dark Energy Camera Legacy Survey (DECaLS) photometric sample, obtaining $\bar{n}_2w_{{\mathrm {p}}}$ measurements across the complete stellar mass range, from $10^{5.3}$ to $10^{11.5}\,{\rm M}_{\odot }$ for blue galaxies, and from $10^{6.3}$ to $10^{11.9}\,{\rm M}_{\odot }$ for red galaxies. We combine $\bar{n}_2w_{{\rm {p}}}$ with $w_{{\rm {p}}}$ measurements from the BGS sample, which is not necessarily complete in stellar mass. Assuming that galaxy bias is primarily determined by stellar mass and colour, we derive the galaxy stellar mass functions (GSMFs) down to $10^{5.3}\,{\rm M}_{\odot }$ for blue galaxies and $10^{6.3}\,{\rm M}_{\odot }$ for red galaxies, while also setting lower limits for smaller masses. The blue and red GSMFs are well described by single and double Schechter functions, respectively, with low-mass end slopes of $\alpha _{\rm {blue}}=-1.54^{+0.02}_{-0.02}$ and $\alpha _{\rm {red}}=-2.50^{+0.08}_{-0.08}$, resulting in the dominance of red galaxies below $10^{7.6}\,{\rm M}_{\odot }$. Stage-IV cosmological photometric surveys, capable of reaching 2–3 mag deeper than DECaLS, present an opportunity to explore the entire galaxy population in the local universe with PAC. This advancement allows us to address critical questions regarding the nature of dark matter, the physics of reionization, and the formation of dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗