Engineering Papers⌕ Search

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Orthogonality of spherical harmonic coefficients

Orthogonality relations are obtained for the spherical harmonic coefficients of functions defined on the surface of a sphere. Following a brief discussion of the orthogonality of Fourier series coefficients, consideration is given to the values averaged over all orientations of the coordinate system of the spherical harmonic coefficients of a function defined on the surface of a sphere that can be expressed in terms of Legendre polynomials for the special case where the function is the sum of two delta functions located at two different points on the sphere, and for the case of an essentially arbitrary function. It is noted that the orthogonality relations derived have found applications in statistical studies of the geomagnetic field.

Mcleod, M. G.↗
NASA NTRS2005

Stroboscopic Interferometer for Measuring Mirror Vibrations

Stroboscopic interferometry is a technique for measuring the modes of vibration of mirrors that are lightweight and, therefore, unavoidably flexible. The technique was conceived especially for modal characterization of lightweight focusing mirror segments to be deployed in outer space; however, the technique can be applied to lightweight mirrors designed for use on Earth as well as the modal investigation of other optical and mechanical structures. To determine the modal structure of vibration of a mirror, it is necessary to excite the mirror by applying a force that varies periodically with time at a controllable frequency. The excitation can utilize sinusoidal, square, triangular, or even asynchronous waveforms. Because vibrational modes occur at specific resonant frequencies, it is necessary to perform synchronous measurements and sweep the frequency to locate the significant resonant modes. For a given mode it is possible to step the phase of data acquisition in order to capture the modal behavior over a single cycle of the resonant frequency. In order to measure interferometrically the vibrational response of the mirror at a given frequency, an interferometer must be suitably aligned with the mirror and adjustably phase-locked with the excitation signal. As in conventional stroboscopic photography, the basic idea in stroboscopic interferometry is to capture an image of the shape of a moving object (in this case, the vibrating mirror) at a specified instant of time in the vibration cycle. Adjusting the phase difference over a full cycle causes the interference fringes to vary over the full range of motion for the mode at the excitation frequency. The interference-fringe pattern is recorded as a function of the phase difference, and, from the resulting data, the surface shape of the mirror for the given mode is extracted. In addition to the interferometer and the mirror to be tested, the equipment needed for stroboscopic interferometry includes an arbitrary-function generator (that is, a signal generator), an oscilloscope, a trigger filter, and an advanced charge-coupled-device (CCD) camera. The optical components are positioned to form a pupil image of the mirror under test on the CCD chip, so that the interference pattern representative of the instantaneous mirror shape is imaged on the CCD chip.

Stahl, H. Philip↗
DOE OSTI2023

Octofitter: Fast, Flexible, and Accurate Orbit Modeling to Detect Exoplanets

As next-generation imaging instruments and interferometers search for planets closer to their stars, they must contend with increasing orbital motion and longer integration times. These compounding effects make it difficult to detect faint planets but also present an opportunity. Increased orbital motion makes it possible to move the search for planets into the orbital domain, where direct images can be freely combined with the radial velocity and proper motion anomaly, even without a confirmed detection in any single epoch. In this paper, we present a fast and differentiable multimethod orbit-modeling and planet detection code called Octofitter. This code is designed to be highly modular and allows users to easily adjust priors, change parameterizations, and specify arbitrary function relations between the parameters of one or more planets. Octofitter further supplies tools for examining model outputs including prior and posterior predictive checks and simulation-based calibration. We demonstrate the capabilities of Octofitter on real and simulated data from different instruments and methods, including HD 91312, simulated JWST/NIRISS aperture masking interferometry observations, radial velocity curves, and grids of images from the Gemini Planet Imager. We show that Octofitter can reliably recover faint planets in long sequences of images with arbitrary orbital motion. This publicly available tool will enable the broad application of multiepoch and multimethod exoplanet detection, which could improve how future targeted ground- and space-based surveys are performed. Finally, its rapid convergence makes it a useful addition to the existing ecosystem of tools for modeling the orbits of directly imaged planets.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1963

Optimizations of Human Restraint Systems for Short-Period Acceleration

A restraint system's main function is to restrain its occupant when his vehicle is subjected to acceleration. If the restraint system is rigid and well-fitting (to eliminate slack) then it will transmit the vehicle acceleration to its occupant without modifying it in any way. Few present-day restraint systems are stiff enough to give this one-to-one transmission characteristic, and depending upon their dynamic characteristics and the nature of the vehicle's acceleration-time history, they will either magnify or attenuate the acceleration. Obviously an optimum restraint system will give maximum attenuation of an input acceleration. In the general case of an arbitrary acceleration input, a computer must be used to determine the optimum dynamic characteristics for the restraint system. Analytical solutions can be obtained for certain simple cases, however, and these cases are considered in this paper, after the concept of dynamic models of the human body is introduced. The paper concludes with a description of an analog computer specially developed for the Air Force to handle completely general mechanical restraint optimization programs of this type, where the acceleration input may be any arbitrary function of time.

Payne, P. R.↗
DOE OSTI2020

Hamiltonian formulations for perturbed dissipationless plasma equations

The Hamiltonian formulations for the perturbed Vlasov–Maxwell equations and the perturbed ideal magnetohydrodynamics (MHD) equations are expressed in terms of the perturbation derivative ∂ F / ∂ ϵ ≡ [ F , S ] of an arbitrary functional F [ ψ ] of the Vlasov–Maxwell fields ψ = ( f , E , B ) or the ideal MHD fields ψ = ( ρ , u , s , B ) , which are assumed to depend continuously on the (dimensionless) perturbation parameter ϵ . In this study, [ , ] denotes the functional Poisson bracket for each set of plasma equations and the perturbation action functional S is said to generate dynamically accessible perturbations of the plasma fields. The new Hamiltonian perturbation formulation introduces a framework for functional perturbation methods in plasma physics and highlights the crucial roles played by polarization and magnetization in Vlasov–Maxwell and ideal MHD perturbation theories. One application considered in this paper is a formulation of plasma stability that guarantees dynamical accessibility and leads to a natural generalization to higher-order perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Optimization techniques in self-similar compressible flow

We investigate the one-dimensional (1D) inviscid compressible flow equations for an ideal gas through the lens of optimization techniques. It is the case that, to our knowledge, optimization analysis applied to the so-called “linear velocity” solutions of the Euler compressible flow equations has not been previously conducted. Through both gradient-based and variational techniques, new variants of well-studied flow scenarios, i.e., self-similar, 1D, linear velocity solution class to idealized inviscid compressible flow equations, are determined, as encoded in both the kinematic and thermodynamic properties of this self-similar solution class. With the kinematics of the said solutions being driven by a self-similar “scale radius” and the thermodynamics being driven separately through the appearance of an arbitrary function, a myriad of new solution classes is possible. Acting as a guide to more realistic physical circumstances as well as discovery, it is the hope that the presented cases serve as the framework for future investigations into the intersection of self-similarity and optimization techniques. Fields of study that may find this work to be of interest include aerodynamic design, flow control, inertial confinement fusion, physics-informed neural networks, and other related areas of interest.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Energy density and spectrum of single-cycle and sub-cycle electromagnetic pulses

Based on the exact solution of the Maxwell equations in the form of a collapsing spherical vector wave specified by an arbitrary function of time, we have calculated the maximum energy density that can be achieved when focusing extremely short pulses of various shapes. It is shown that our earlier formula expressing the maximum energy density in terms of the spectrum parameters for Gaussian quasi-monochromatic pulses is approximately (with an accuracy of ±40 %) valid for the main types of extremely short pulses. (extremely short laser pulses)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Metric DBSCAN

SAND2025-11725O Metric DBSCAN is an implementation of the popular DBSCAN clustering algorithm that works in general metric spaces. DBSCAN is a clustering algorithm, a fundamental building block in machine learning. It takes a set of objects and, given some notion of distance, identifies coherent groups of objects. With Metric DBSCAN, users can provide an arbitrary function to compute distance. Nearly all existing implementations of DBSCAN restrict distance to one of a few formulations. Metric DBScan accomplishes this cleanly and efficiently. The Python source code is on Github. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Dalbey, Keith↗
DOE OSTI2021

End-to-End Differentiable Modeling and Management of the Environment

Focal Area: (2) Data acquisition and assimilation enabled by machine learning, AI, and advanced methods including experimental/network design/optimization. We emphasize the importance of leveraging optimization techniques from AI/machine learning (ML) to solve challenging problems in Earth system modeling. Science Challenge: Automatic differentiation has had a transformative effect on ML by allowing the calculation of gradients of arbitrary functions in an incredibly large class of models. We can potentially realize similar improvements in parameter estimation and control for Earth system models (ESMs) by reimplementing them in computational frameworks from ML. Practitioners working with large (>10 7 parameters) models in ML and AI can obtain good predictive performance in a range of spatiotemporal tasks by making use of optimization via stochastic gradient descent and incorporating prior knowledge at multiple levels. We propose writing ESMs in open-source computational frameworks such as Torch, Tensorflow, and JAX to greatly expand the scope of environmental forecasting and management challenges, which can be addressed by leveraging automatic differentiation and gradient descent-like algorithms. We do not call for a wholesale replacement of physical models with data-driven surrogates, but rather advocate for interleaving physical and empirical equations in a manner that is most faithful to the extent of our scientific knowledge and observational data. Central to this topic is the merging of differentiable physical simulations with differentiable optimization layers, which are now both beginning to come to the forefront.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Stochastic Approaches for Calculating and Aggregating Detection Probabilities for Nuclear Material Diversion

The authors built and tested a stochastic simulation to estimate achieved detection probabilities (DPs) on a stratum basis, over a tailorable range of diverted amounts from 0 to 2 significant quantities (SQ), using typical International Atomic Energy Agency (IAEA) inspection data: i.e., SQ in stratum, number of items, number of gross/partial/bias defect measurements conducted, and realistic relative standard deviation (RSD) values for typical IAEA verification measurements. For bulk strata, the model calculates achieved DP at 0.01 SQ diversion increments; for item strata, the model calculates DP using the smallest realistic diversion increment (e.g., a plate, pin, or coupon). After successfully benchmarking against IAEA deterministic models, the simulation was used to test the sensitivity of DP to certain standard assumptions and selected input parameters. First, the equal defect assumption was tested; the results suggest significant complexity in the effectiveness of partial defect measurements. Next, the authors explored the sensitivity of DP to the assumed RSD of attribute tests. Then, the authors compared non-normal models for instrument performance (e.g., logistic, step, or arbitrary functions) to the typical results from a normal distribution (characterized by RSD). This last comparison was supplemented with experimentally derived performance data for an HM-5 gamma spectrometer. The HM-5 was used to make enrichment measurements on both LEU and HEU MTR fuel elements as plates were removed, and the results fit with logistic and step curves and applied in the simulation. These stochastic DP results were compared to DP estimates from a deterministic model assuming a normal curve and typical RSD, yielding insights that could improve effectiveness in the field. These early results illustrate the potential of stochastic models to better understand achieved DP and to improve safeguards effectiveness.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2025

Optimization Through Multi-Fidelity Modeling

We present a novel method for optimizing parameter selection for simulations with an evaluation budget. We start with an existing method for building a multi-fidelity model out of many low-fidelity simulations and few high-fidelity simulations. We propose a novel method to simplify parameter selection without sacrificing performance. We verify these results and compare with existing literature. Next, we propose a novel algorithm which uses this difference model to suggest new points in the parameter design space to simulate. We add each point we simulate to the model to improve its quality for the next iteration. The algorithm trades off reducing the uncertainty of the existing model with optimization of the objective. The first is more useful when a large fraction of the computation budget remains. The second is more useful when a small fraction of the computation budget remains. Our method converges to the optimum by using a high-fidelity evaluation for just 16 of the 427 points. Our method is general enough to work if there is no low-fidelity model. Furthermore, it is agnostic to the underlying physics of the problem. Therefore, both the low-fidelity and high-fidelity models can be generated by any arbitrary function, including simulations and physical experiments.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Computer program for solving compressible nonsimilar-boundary-layer equations for laminar, transitional, or turbulent flows of a perfect gas

A computer program is described which solves the compressible laminar, transitional, or turbulent boundary-layer equations for planar or axisymmetric flows. Three-point implicit difference relations are used to reduce the momentum and energy equations to finite-difference form. These equations are solved simultaneously without iteration. Turbulent flow is treated by the inclusion of either a two-layer eddy-viscosity model or a mixing-length formulation. The eddy conductivity is related to the eddy viscosity through a static turbulent Prandtl number which may be an arbitrary function of the distance from the wall boundary. The transitional boundary layer is treated by the inclusion of an intermittency function which modifies the fully turbulent model. The laminar-boundary-layer equations are recovered when the intermittency is zero, and the fully turbulent equations are solved when the intermittency is unity.

Price, J. M.↗
NASA NTRS1971

Thermal conductivity, electrical resistivity, and thermopower of aerospace alloys from 4 to 300 K. 6: Fe-22Cr-13Ni-5Mn stainless steel

The equipment and techniques for determining the thermal conductivity, electrical resistivity Lorenz ratio, and thermopower characteristics of Fe-22Cr-13Ni-5Mn stainless steel are discussed. The dimensions of the specimen and its preparation are described. The experimental data are represented by arbitrary functions over the entire range and smooth tables are generated from these functions.

Hust, J. G.↗
NASA NTRS1972

Lunar electromagnetic scattering. 1: Propagation parallel to the diamagnetic cavity axis

An analytic theory is developed for the time dependent magnetic fields inside the Moon and the diamagnetic cavity when the interplanetary electromagnetic field fluctuation propagates parallel to the cavity axis. The Moon model has an electrical conductivity which is an arbitrary function of radius. The lunar cavity is modelled by a nonconducting cylinder extending infinitely far downstream. For frequencies less than about 50 Hz, the cavity is a cylindrical waveguide below cutoff. Thus, cavity field perturbations due to the Moon do not propagate down the cavity, but are instead attenuated with distance downstream from the Moon.

Schwartz, K.↗
NASA NTRS1972

Development of flow in the entrance region of a converging channel.

An approximate solution for the laminar flow of an incompressible viscous fluid in the entrance region of a converging channel is obtained. The radial velocity distribution at the entrance of the channel is taken to be a symmetric but otherwise arbitrary function of the angular position. Expressions for the velocity components and pressure are given. The case of the uniformly distributed entrance velocity is considered as an example.

Atabek, H. B.↗
NASA NTRS1976

Thermal stresses in a spherical pressure vessel having temperature-dependent, transversely isotropic, elastic properties

Rayleigh-Ritz and modified Rayleigh-Ritz procedures are used to construct approximate solutions for the response of a thick-walled sphere to uniform pressure loads and an arbitrary radial temperature distribution. The thermoelastic properties of the sphere are assumed to be transversely isotropic and nonhomogeneous; variations in the elastic stiffness and thermal expansion coefficients are taken to be an arbitrary function of the radial coordinate and temperature. Numerical examples are presented which illustrate the effect of the temperature-dependence upon the thermal stress field. A comparison of the approximate solutions with a finite element analysis indicates that Ritz methods offer a simple, efficient, and relatively accurate approach to the problem.

Tauchert, T. R.↗
NASA NTRS1977

Large-scale inhomogeneity of the Universe - Spherically symmetric models

Dust-filled Tolman-Bondi (1934, 1947) spherically symmetric cosmological models are studied with application to superclustering of galaxies. Asymptotic expressions are given for the density and curvature contrast, valid near the singularity and at long times. The density profile of infalling material is expressed in terms of a single arbitrary function (the distribution of total energy per unit mass). Measurement of the distribution of excess galaxy counts and deviations from the Hubble flow in the vicinity of a supercluster could provide a determination of the cosmological deceleration parameter.

Silk, J.↗