SEARCH · Search NASA

Results for “KERNEL FUNCTION”

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

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Recovering pointwise values of discontinuous data within spectral accuracy

The pointwise values of a function, f(x), can be accurately recovered either from its spectral or pseudospectral approximations, so that the accuracy solely depends on the local smoothness of f in the neighborhood of the point x. Most notably, given the equidistant function grid values, its intermediate point values are recovered within spectral accuracy, despite the possible presence of discontinuities scattered in the domain. (Recall that the usual spectral convergence rate decelerates otherwise to first order, throughout). To this end, a highly oscillatory smoothing kernel is employed in contrast to the more standard positive unit-mass mollifiers. In particular, post-processing of a stable Fourier method applied to hyperbolic equations with discontinuous data, recovers the exact solution modulo a spectrally small error. Numerical examples are presented.

Gottlieb, D.↗

DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗

NASA NTRS2001

Aeroelastic Response of Nonlinear Wing Section by Functional Series Technique

This paper addresses the problem of the determination of the subcritical aeroelastic response and flutter instability of nonlinear two-dimensional lifting surfaces in an incompressible flow-field via indicial functions and Volterra series approach. The related aeroelastic governing equations are based upon the inclusion of structural and damping nonlinearities in plunging and pitching, of the linear unsteady aerodynamics and consideration of an arbitrary time-dependent external pressure pulse. Unsteady aeroelastic nonlinear kernels are determined, and based on these, frequency and time histories of the subcritical aeroelastic response are obtained, and in this context the influence of the considered nonlinearities is emphasized. Conclusions and results displaying the implications of the considered effects are supplied.

Silva, Walter A.↗

NASA NTRS2000

Aeroelastic Response of Nonlinear Wing Section By Functional Series Technique

Marzocca, Piergiovanni↗

DOE OSTI2026

A Discrete Hankel Transform Approach to Nuclear Data Processing for Fusion Applications

This study introduces advancements to the numerical solutions employed in the processing of nuclear data for fusion applications. It leverages the convolution theorem and Fourier transform techniques to enhance computational efficiency and broaden applicability. Building upon a previously reported discrete Hankel transform approach for Doppler broadening, this work refines the solution of convolution integrals central to these applications. The methodology provides a general and unified framework for evaluating any convolution operation, regardless of whether the underlying problem involves temperature effects in nuclear reactions. The applicability to the nuclear data processing for fusion is demonstrated by deriving the convolution integrals for some of the fusion-related quantities. As before, the convolution operation utilizes a Gaussian-based kernel; however, the discrete Hankel transform of order $𝛼$ = $\frac{1}{2}$ is now applied to the forward Fourier transform of the nonkernel argument, rather than the inverse Fourier transform. This modification eliminates the need for the integration of the nonkernel, cross section–based function, which is a step that posed challenges for certain pointwise cross-section representations. It also removes the requirement for cross-section linearization. Optimized for graphics processing unit architectures, the approach significantly improves computational performance. These advancements are currently under evaluation as the foundation for the next-generation thermonuclear data file processing codes being developed at Lawrence Livermore National Laboratory.

Nuclear science and engineering↗

NASA NTRS2005

Development of Ionic Liquid Monopropellants for In-Space Propulsion

A family of new, low toxicity, high energy monopropellants is currently being evaluated at NASA Marshall Space Flight Center for in-space rocket engine applications such as reaction control engines. These ionic liquid monopropellants, developed in recent years by the Air Force Research Laboratory, could offer system simplification, less in-flight thermal management, and reduced handling precautions, while increasing propellant energy density as compared to traditional storable in-space propellants such as hydrazine and nitrogen tetroxide. However, challenges exist in identifying ignition schemes for these ionic liquid monopropellants, which are known to burn at much hotter combustion temperatures compared to traditional monopropellants such as hydrazine. The high temperature combustion of these new monopropellants make the use of typical ignition catalyst beds prohibitive since the catalyst cannot withstand the elevated temperatures. Current research efforts are focused on monopropellant ignition and burn rate characterization, parameters that are important in the fundamental understanding of the monopropellant behavior and the eventual design of a thruster. Laboratory studies will be conducted using alternative ignition techniques such as laser-induced spark ignition and hot wire ignition. Ignition delay, defined as the time between the introduction of the ignition source and the first sign of light emission from a developing flame kernel, will be measured using Schlieren visualization. An optically-accessible liquid monopropellant burner, shown schematically in Figure 1 and similar in design to apparatuses used by other researchers to study solid and liquid monopropellants, will be used to determine propellant burn rate as a function of pressure and initial propellant temperature. The burn rate will be measured via high speed imaging through the chamber s windows.

Blevins, John A.↗

NASA NTRS1980

Reducing the error of geoid undulation computations by modifying Stokes' function

The truncation theory as it pertains to the calculation of geoid undulations based on Stokes' integral, but from limited gravity data, is reexamined. Specifically, the improved procedures of Molodenskii et al. are shown through numerical investigations to yield substantially smaller errors than the conventional method that is often applied in practice. In this improved method, as well as in a simpler alternative to the conventional approach, the Stokes' kernel is suitably modified in order to accelerate the rate of convergence of the error series. These modified methods, however, effect a reduction in the error only if a set of low-degree potential harmonic coefficients is utilized in the computation. Consider, for example, the situation in which gravity anomalies are given in a cap of radius 10 deg and the GEM 9 (20,20) potential field is used. Then, typically, the error in the computed undulation (aside from the spherical approximation and errors in the gravity anomaly data) according to the conventional truncation theory is 1.09 m; with Meissl's modification it reduces to 0.41m, while Molodenskii's improved method gives 0.45 m. A further alteration of Molodenskii's method is developed and yields an RMS error of 0.33 m. These values reflect the effect of the truncation, as well as the errors in the GEM 9 harmonic coefficients. The considerable improvement, suggested by these results, of the modified methods over the conventional procedure is verified with actual gravity anomaly data in two oceanic regions, where the GEOS-3 altimeter geoid serves as the basis for comparison. The optimal method of truncation, investigated by Colombo, is extremely ill-conditioned. It is shown that with no corresponding regularization, this procedure is inapplicable.

Jekeli, C.↗

DOE OSTI2026

Comment on “QCD factorization with multihadron fragmentation functions”

We make several comments on the recent work in Rogers et al. [Phys. Rev. D 111, 056001 (2025)] while also reaffirming and adding to the work in Pitonyak et al. [Phys. Rev. Lett. 132, 011902 (2024)]. We show that the factorization formula for 𝑒 + ⁢𝑒 − → (ℎ 1 ⋯ ℎ 𝑛 )⁢𝑋 in Rogers et al. is equivalent to a version one can derive using the definition of a 𝑛-hadron fragmentation function (FF) introduced in Pitonyak et al.. In addition, we scrutinize how to generalize the number density definition of a single-hadron FF to a 𝑛-hadron FF, arguing that the definition given in Pitonyak et al. should be considered the standard one while the definition in Rogers et al. has no clear interpretation. We also emphasize that the evolution equations for dihadron FFs (DiFFs) in Pitonyak et al. have the same splitting functions as those for single-hadron FFs. Therefore, the DiFF (and 𝑛-hadron FF) definitions in Pitonyak et al. have a natural number density interpretation and, contrary to what is stated in Rogers et al., are consistent with collinear factorization using the standard hard factors and evolution kernels.

Pitonyak, D. [Lebanon Valley College, Annville, PA↗

DOE OSTI2026

AstraAI v1

AstraAI is an open-source, structure-aware AI coding agent designed for large scientific and DOE-HPC codebases such as AMReX-based applications. Unlike general-purpose coding assistants, AstraAI combines retrieval-augmented generation (RAG) with compiler-level Abstract Syntax Tree (AST) analysis to perform precise, scope-constrained code modifications. It identifies exact function spans, enforces locality of edits, and maintains cross-file invariants, enabling deterministic and build-safe transformations in complex C++/GPU environments. AstraAI is intended for developers working on large, evolving HPC frameworks where correctness, reproducibility, and structural integrity are critical. Typical use cases include modifying physics kernels, updating GPU device lambdas, and performing multi-file refactors without breaking compilation or runtime semantics. Compared to conventional LLM-based coding agents - even those with repository access - AstraAI provides structural guarantees rather than free-form text patches. It minimizes unintended diffs, prevents scope drift, preserves formatting and build stability, and reduces structural hallucinations. By integrating compiler tooling directly into the generation loop, AstraAI transforms AI-assisted coding from probabilistic text editing into deterministic, structure-preserving program transformation suitable for mission-critical scientific software.

Natarajan, Mahesh [Lawrence Berkeley National Labo↗

NASA NTRS1985

A new approach to approximating the linear quadratic optimal control law for hereditary systems with control delays

A factorization approach is presented for deriving approximations to the optimal feedback gain for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the feedback kernels.

Milman, M. H.↗

NASA NTRS1991

Mark 4A antenna control system data handling architecture study

A high-level review was conducted to provide an analysis of the existing architecture used to handle data and implement control algorithms for NASA's Deep Space Network (DSN) antennas and to make system-level recommendations for improving this architecture so that the DSN antennas can support the ever-tightening requirements of the next decade and beyond. It was found that the existing system is seriously overloaded, with processor utilization approaching 100 percent. A number of factors contribute to this overloading, including dated hardware, inefficient software, and a message-passing strategy that depends on serial connections between machines. At the same time, the system has shortcomings and idiosyncrasies that require extensive human intervention. A custom operating system kernel and an obscure programming language exacerbate the problems and should be modernized. A new architecture is presented that addresses these and other issues. Key features of the new architecture include a simplified message passing hierarchy that utilizes a high-speed local area network, redesign of particular processing function algorithms, consolidation of functions, and implementation of the architecture in modern hardware and software using mainstream computer languages and operating systems. The system would also allow incremental hardware improvements as better and faster hardware for such systems becomes available, and costs could potentially be low enough that redundancy would be provided economically. Such a system could support DSN requirements for the foreseeable future, though thorough consideration must be given to hard computational requirements, porting existing software functionality to the new system, and issues of fault tolerance and recovery.

Briggs, H. C.↗

NASA NTRS1961

Consistent Pl Analysis of Aqueous Uranium-235 Critical Assemblies

The lethargy-dependent equations of the consistent Pl approximation to the Boltzmann transport equation for slowing down neutrons have been used as the basis of an IBM 704 computer program. Some of the effects included are (1) linearly anisotropic center of mass elastic scattering, (2) heavy element inelastic scattering based on the evaporation model of the nucleus, and (3) optional variation of the buckling with lethargy. The microscopic cross-section data developed for this program covered 473 lethargy points from lethargy u = 0 (10 Mev) to u = 19.8 (0.025 ev). The value of the fission neutron age in water calculated here is 26.5 square centimeters; this value is to be compared with the recent experimental value given as 27.86 square centimeters. The Fourier transform of the slowing-down kernel for water to indium resonance energy calculated here compared well with the Fourier transform of the kernel for water as measured by Hill, Roberts, and Fitch. This method of calculation has been applied to uranyl fluoride - water solution critical assemblies. Theoretical results established for both unreflected and fully reflected critical assemblies have been compared with available experimental data. The theoretical buckling curve derived as a function of the hydrogen to uranium-235 atom concentration for an energy-independent extrapolation distance was successful in predicting the critical heights of various unreflected cylindrical assemblies. The critical dimensions of fully water-reflected cylindrical assemblies were reasonably well predicted using the theoretical buckling curve and reflector savings for equivalent spherical assemblies.

Fieno, Daniel↗

NASA NTRS1994

A spectral boundary integral equation method for the 2-D Helmholtz equation

In this paper, we present a new numerical formulation of solving the boundary integral equations reformulated from the Helmholtz equation. The boundaries of the problems are assumed to be smooth closed contours. The solution on the boundary is treated as a periodic function, which is in turn approximated by a truncated Fourier series. A Fourier collocation method is followed in which the boundary integral equation is transformed into a system of algebraic equations. It is shown that in order to achieve spectral accuracy for the numerical formulation, the nonsmoothness of the integral kernels, associated with the Helmholtz equation, must be carefully removed. The emphasis of the paper is on investigating the essential elements of removing the nonsmoothness of the integral kernels in the spectral implementation. The present method is robust for a general boundary contour. Aspects of efficient implementation of the method using FFT are also discussed. A numerical example of wave scattering is given in which the exponential accuracy of the present numerical method is demonstrated.

Hu, Fang Q.↗

NASA NTRSDate not supplied

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Han Byul Woo↗

NASA NTRSDate not supplied

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Rocket Launch Seismicity↗

DOE OSTI2024

Vacancy-Induced Tunable Kondo Effect in Twisted Bilayer Graphene

In single sheets of graphene, vacancy-induced states have been shown to host an effective spin-1/2 hole that can be Kondo screened at low temperatures. Here, we show how these vacancy-induced impurity states survive in twisted bilayer graphene (TBG), which thus provides a tunable system to probe the critical destruction of the Kondo effect in pseudogap hosts. Ab initio calculations and atomic-scale modeling are used to determine the nature of the vacancy states in the vicinity of the magic angle in TBG, demonstrating that the vacancy can be treated as a quantum impurity. Utilizing this insight, we construct an Anderson impurity model with a TBG host that we solve using the numerical renormalization group combined with the kernel polynomial method. We determine the phase diagram of the model and show how there is a strict dichotomy between vacancies in the AA/BB versus AB/BA tunneling regions. In AB/BA vacancies, the Kondo temperature at the magic angle develops a broad distribution with a tail to vanishing temperatures due to multifractal wave functions at the magic angle. Finally, we argue that scanning tunneling microscopy in the vicinity of the vacancy can act as a probe of both the critical single-particle states and the underlying many-body ground state in magic-angle TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1987

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary systems. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗

NASA NTRS1988

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary schemes. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗