SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Quantum Probability Cancellation Due to a Single-Photon State

When an N-photon state enters a lossless symmetric beamsplitter from one input port, the photon distribution for the two output ports has the form of Bernouli Binormial, with highest probability at equal partition (N/2 at one outport and N/2 at the other). However, injection of a single photon state at the other input port can dramatically change the photon distribution at the outputs, resulting in zero probability at equal partition. Such a strong deviation from classical particle theory stems from quantum probability amplitude cancellation. The effect persists even if the N-photon state is replaced by an arbitrary state of light. A special case is the coherent state which corresponds to homodyne detection of a single photon state and can lead to the measurement of the wave function of a single photon state.

Ou, Z. Y.↗

DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

The connection between nonzero density and spontaneous symmetry breaking for interacting scalars

We consider U(1)-symmetric scalar quantum field theories at zero temperature. At nonzero charge densities, the ground state of these systems is usually assumed to be a superfluid phase, in which the global symmetry is spontaneously broken along with Lorentz boosts and time translations. We show that, in d > 2 spacetime dimensions, this expectation is always realized at one loop for arbitrary non-derivative interactions, confirming that the physically distinct phenomena of nonzero charge density and spontaneous symmetry breaking occur simultaneously in these systems. We quantify this result by deriving universal scaling relations for the symmetry breaking scale as a function of the charge density, at low and high density. Moreover, we show that the critical value of μ above which a nonzero density develops coincides with the pole mass in the unbroken, Poincaré invariant vacuum of the theory. The same conclusions hold non-perturbatively for an O(N) theory with quartic interactions in d = 3 and 4, at leading order in the 1/N expansion. We derive these results by computing analytically the zero-temperature, finite-μ one-loop effective potential, paying special attention to subtle points related to the iε terms. We check our results against the one-loop low-energy effective action for the superfluid phonons in λΦ 4 theory in d = 4 previously derived by Joyce and ourselves, which we further generalize to arbitrary potential interactions and arbitrary dimensions. As a byproduct, we find analytically the one-loop scaling dimension of the lightest charge-n operator for the λΦ 6 conformal superfluid in d = 3, at leading order in 1/n, reproducing a numerical result of Badel et al. For a λΦ 4 superfluid in d = 4, we also reproduce the Lee-Huang-Yang relation and compute relativistic corrections to it. Finally, we discuss possible extensions of our results beyond perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Imaging elliptically polarized infrared near-fields on nanoparticles by strong-field dissociation of functional surface groups

Abstract We investigate the strong-field ion emission from the surface of isolated silica nanoparticles aerosolized from an alcoholic solution, and demonstrate the applicability of the recently reported near-field imaging at 720 nm [Rupp et al., Nat. Comm., 10(1):4655, 2019] to longer wavelength (2 $$\mu $$ μ m) and polarizations with arbitrary ellipticity. Based on the experimental observations, we discuss the validity of a previously introduced semi-classical model, which is based on near-field driven charge generation by a Monte-Carlo approach and classical propagation. We furthermore clarify the role of the solvent in the surface composition of the nanoparticles in the interaction region. We find that upon injection of the nanoparticles into the vacuum, the alcoholic solvent evaporates on millisecond time scales, and that the generated ions originate predominantly from covalent bonds with the silica surface rather than from physisorbed solvent molecules. These findings have important implications for the development of future theoretical models of the strong-field ion emission from silica nanoparticles, and the application of near-field imaging and reaction dynamics of functional groups on isolated nanoparticles. Graphical abstract

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Unidirectional luminescence from InGaN/GaN quantum-well metasurfaces

III–nitride light-emitting diodes (LEDs) are the backbone of ubiquitous lighting and display applications. Imparting directional emission is an essential requirement for many LED implementations. Although optical packaging, nanopatterning and surface roughening4 techniques can enhance LED extraction, directing the emitted light requires bulky optical components. Optical metasurfaces provide precise control over transmitted and reflected waveforms, suggesting a new route for directing light emission. However, it is difficult to adapt metasurface concepts for incoherent light emission, due to the lack of a phase-locking incident wave. Here, we demonstrate a metasurface-based design of InGaN/GaN quantum-well structures that generate narrow, unidirectional transmission and emission lobes at arbitrary engineered angles. We further demonstrate 7-fold and 100-fold enhancements of total and air-coupled external quantum efficiencies, respectively. The results present a new strategy for exploiting metasurface functionality in light-emitting devices.

inorganic LEDs↗

DOE OSTI2021

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1981

Search strategy effects on PN acquisition performance

The present paper focusses on 'random' and 'expanding window' PN acquisition search strategies and analytically develops the PN acquisition time statistics as functions of salient system parameters such as prediction SNR, detection and false alarm probabilities and a priori information on epoch location. The significance of this analysis is its general applicability to arbitrary postdetection processing schemes. Computed performance results account for the above salient parameters, wherein sequential detection is employed in conjunction with random and selected expanding window search strategies.

Weinberg, A.↗

NASA NTRS1992

Flexure-torsion behavior of sheat-deformable beams with applications to aircraft wing sections

The flexure-torsion behavior of a tip-loaded cantilever beam with an arbitrary cross-section is studied using Saint-Venant's semi-inverse method along with a power series solution for the out-of-plane flexure and torsion warping functions. The power series coefficients are determined by solving a set of variationally derived linear algebraic equations. For complex cross-sections, the calculated coefficients represent a 'best-fit approximation' to the exact warping function. The resulting warping functions are used to determine the cross-section properties including: the torsion constant, shear deformation coefficients, shear correction factors, and the shear center location. A new linear relation is developed for locating the shear center using the Saint-Venant flexure and torsion solutions, where the twist rate is zero about the line of shear centers (not the centroidal axis). Numerical results are presented for a triangular cross-section and different NACA airfoils.

Kosmatka, J. B.↗

NASA NTRS1992

A finite element-boundary integral method for cavities in a circular cylinder

Conformal antenna arrays offer many cost and weight advantages over conventional antenna systems. However, due to a lack of rigorous mathematical models for conformal antenna arrays, antenna designers resort to measurement and planar antenna concepts for designing non-planar conformal antennas. Recently, we have found the finite element-boundary integral method to be very successful in modeling large planar arrays of arbitrary composition in a metallic plane. We extend this formulation to conformal arrays on large metallic cylinders. In this report, we develop the mathematical formulation. In particular, we discuss the shape functions, the resulting finite elements and the boundary integral equations, and the solution of the conformal finite element-boundary integral system. Some validation results are presented and we further show how this formulation can be applied with minimal computational and memory resources.

Kempel, Leo C.↗

NASA NTRS1993

Modeling the densification of metal matrix composite monotape

We present a first model that enables prediction of the density (and its time evolution) of a monotape lay-up subjected to a hot isostatic or vacuum hot pressing consolidation cycle. Our approach is to break down the complicated (and probabilistic) consolidation problem into simple, analyzable parts and to combine them in a way that correctly represents the statistical aspects of the problem, the change in the problem's interior geometry, and the evolving contributions of the different deformation mechanisms. The model gives two types of output. One is in the form of maps showing the relative density dependence upon pressure, temperature, and time for step function temperature and pressure cycles. They are useful for quickly determining the best place to begin developing an optimized process. The second gives the evolution of density over time for any (arbitrary) applied temperature and pressure cycle. This has promise for refining process cycles and possibly for process control. Examples of the models application are given for Ti3Al + Nb, gamma TiAl, Ti6Al4V, and pure aluminum.

Elzey, D. M.↗

NASA NTRS2017

Trajectory Specification for Terminal Air Traffic: Pairwise Conflict Detection and Resolution

Trajectory Specification is the explicit bounding and control of aircraft trajectories such that the position at any point in time is constrained to a precisely defined volume of space. The bounding space is defined by cross-track, along-track, and vertical tolerances relative to a reference trajectory that specifies position as a function of time. The tolerances are dynamic and will be based on the aircraft navigation capabilities and the current traffic situation. Assuming conformance, Trajectory Specification can guarantee safe separation for an arbitrary period of time even in the event of an air traffic control (ATC) system or datalink failure; hence it can help to achieve the high level of safety and reliability needed for ATC automation. It can also reduce the reliance on tactical backup systems during normal operation. This paper applies it to the terminal area around a major airport and presents algorithms and software for detecting and resolving conflicts. A representative set of pairwise conflicts was generated, and a fast-time simulation was run on them. All conflicts were successfully resolved in real time, demonstrating the computational feasibility of the concept.

terminal area↗

DOE OSTI2023

Full forward model of galaxy clustering statistics with AbacusSummit light cones

ABSTRACT Novel summary statistics beyond the standard 2-point correlation function (2PCF) are necessary to capture the full astrophysical and cosmological information from the small-scale (r < 30h−1Mpc) galaxy clustering. However, the analysis of beyond-2PCF statistics on small scales is challenging because we lack the appropriate treatment of observational systematics for arbitrary summary statistics of the galaxy field. In this paper, we develop a full forward modelling pipeline for a wide range of summary statistics using the large high-fidelity AbacusSummit light cones that account for many systematic effects as well as remain flexible and computationally efficient to enable posterior sampling. We apply our forward model approach to a fully realistic mock galaxy catalog and demonstrate that we can recover unbiased constraints on the underlying galaxy–halo connection model using two separate summary statistics: the standard 2PCF and the novel k-th nearest neighbour (kNN) statistics, which are sensitive to correlation functions of all orders. We will demonstrate its strong constraining power on extended galaxy–halo connection models and cosmology in follow up papers. We expect this to become a powerful approach when applying to upcoming surveys such as DESI where we can leverage a multitude of summary statistics across a wide redshift range to maximally extract information from the non-linear scales.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

4D Printing of Stretchable Supercapacitors via Hybrid Composite Materials

Stretchable supercapacitors (SCs) have attracted significant attention in developing power-independent stretchable electronic systems due to their intrinsic energy storage function and unique mechanical properties. Most current SCs are generally limited by their low stretchability, complicated fabrication process, and insufficient performance and robustness. This study presents a facile method to fabricate arbitrary-shaped stretchable electrodes via 4D printing of conductive composite from reduced graphene oxide, carbon nanotube, and poly(3,4-ethylenedioxythiophene) polystyrene sulfonate. The electrode patterns of an arbitrary shape can be deposited onto prestretched substrates by aerosol-jet printing, then self-organized origami (ridge) patterns are generated after releasing the substrates from holding stretchers due to the mismatched strains. The stretchable electrodes demonstrate superior mechanical robustness and stretchability without sacrificing its outstanding electrochemical performance. The symmetric SC prototype possesses a gravimetric capacitance of ≈21.7 F g -1 at a current density of 0.5 A g -1 and a capacitance retention of ≈85.8% from 0.5 to 5 A g -1 . A SC array with arbitrary-shaped electrodes is also fabricated and connected in series to power light-emitting diode patterns for large-scale applications. Finally, the proposed method paves avenues for scalable manufacturing of future energy-storage devices with controlled extensibility and high electrochemical performance.

36 MATERIALS SCIENCE↗

DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗

DOE OSTI2022

ExaTN: Scalable GPU-Accelerated High-Performance Processing of General Tensor Networks at Exascale

We present ExaTN (Exascale Tensor Networks), a scalable GPU-accelerated C++ library which can express and process tensor networks on shared- as well as distributed-memory high-performance computing platforms, including those equipped with GPU accelerators. Specifically, ExaTN provides the ability to build, transform, and numerically evaluate tensor networks with arbitrary graph structures and complexity. It also provides algorithmic primitives for the optimization of tensor factors inside a given tensor network in order to find an extremum of a chosen tensor network functional, which is one of the key numerical procedures in quantum many-body theory and quantum-inspired machine learning. Numerical primitives exposed by ExaTN provide the foundation for composing rather complex tensor network algorithms. We enumerate multiple application domains which can benefit from the capabilities of our library, including condensed matter physics, quantum chemistry, quantum circuit simulations, as well as quantum and classical machine learning, for some of which we provide preliminary demonstrations and performance benchmarks just to emphasize a broad utility of our library.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

Tail modeling in a stretched magnetosphere. I - Methods and transformations

A new method is developed for representing the magnetospheric field B as a distorted dipole field. Because Delta-B = 0 must be maintained, such a distortion may be viewed as a transformation of the vector potential A. The simplest form is a one-dimensional 'stretch transformation' along the x axis, concisely represented by the 'stretch function' f(x), which is also a convenient tool for representing features of the substorm cycle. One-dimensional stretch transformations are extended to spherical, cylindrical, and parabolic coordinates and then to arbitrary coordinates. It is shown that distortion transformations can be viewed as mappings of field lines from one pattern to another; the final result only requires knowledge of the field and not of the potentials. General transformations in Cartesian and arbitrary coordinates are derived, and applications to field modeling, field line motion, MHD modeling, and incompressible fluid dynamics are considered.

Stern, David P.↗

NASA NTRS1985

Finite-time arbitrary-motion unsteady cascade airfoil theory for helicopter rotors in hover

A complete and detailed derivation of finite-time arbitrary-motion cascade theory is presented for both Laplace and frequency domains. This theory includes the effect of returning wakes for both single- and multibladed rotors. The generalized cascade lift-deficiency function (CLDF) is shown to be consistent with the generalized Theodorsen lift-deficiency function when the wake spacing approaches infinity or when the reduced frequency tends to infinity. This function predicts correct zero-reduced-frequency limit. Accurate and efficient numerical procedures are presented for the evaluation of the CLDF. Numerical examples comparing the CLDF with Loewy's lift deficiency function in frequency domain are presented. Accurate Pade approximants of the CLDF are constructed using a Bode-plot approach which allows for complex poles.

Dinyavari, M. A. H.↗