Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytic 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
DOE OSTI2020

Analytical Functions for 200 West Pump-and-Treat SCADA Sensor Data

Historical operations at the U.S. Department of Energy’s Hanford Site included disposal of waste fluids to the subsurface in the 200 West Area on the Hanford Central Plateau. Subsequent infiltration of fluids has resulted in groundwater contamination with carbon tetrachloride, nitrate, uranium, technetium-99, and other contaminants. A pump-and-treat (P&T) system, with an extraction/injection well network and an aboveground treatment plant, was implemented as part of interim and final remedies in the 200 West area. The HYPATIA single-page web application (part of the SOCRATES suite) is being developed to provide access to and analysis of chemistry and treatment facility sensor data for this 200 West P&T system. For the web application, analytical algorithms were developed to perform summing, differencing, smoothing, outlier detection, change-point detection, mass flow rate, and injectivity calculations on the data. Candidate algorithms were identified and tested, with the best-performing algorithms then assembled for implementation in HYPATIA. Because HYPATIA is hosted on the Amazon Web Services (AWS) cloud computing platform, algorithms were implemented in a back-end AWS Lambda function that can be called by the HYPATIA front end. The Lambda function applies the requested data processing to specified data via functions written in R, Python, and JavaScript. Development, testing, and review of the data analysis algorithms was completed under an NQA-1 quality program. This new HYPATIA functionality will provide information to support site decisions regarding P&T system performance and optimization.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models

Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

A simple analytical function for bidirectional reflectance

A simple bidirectional reflectance model based on physical scattering laws is developed to calculate the bidirectional reflectance of a wide variety of surfaces. Chandrasekhar's radiative transfer solution, obtained for the anisotropically scattering semi-infinite medium, is used to compute the multiple-scattered radiances. However, since the radiative transfer solution does not take architectural effects into consideration, we adopt Hapke's (1986) approach and add an empirical term to explain the hot-spot phenomenon and use Cox and Munk's (1954) formulation to take into account the specular reflection. The physical parameters of the model are retrieved from bidirectional reflectance measurements. For each surface, only one set of model parameters is needed for application to all illumination and viewing geometries. The validity of the model is established by comparing the computed and measured reflectances for sets of viewing and illumination angles that were not included in the inversion algorithm. Good agreement is shown between the model-computed and the observed reflectances for dense prairie vegetation canopies, a sparse desert scrub community, a plowed agricultural field, and an alkali flat.

Ahmad, Suraiya P.↗
DOE OSTI2026

Precise Modeling of a Complex Solenoidal Magnetic Field Using a Combination of Analytic Functions and a PINN

We demonstrate an iterative approach to modeling a sparsely measured magnetic field in a large-bore solenoid. This approach uses a hybrid of traditional and machine learning techniques. The traditional technique is a linear least-squares fit using a series solution to Laplace's equation, while the machine learning technique involves the training of a physics-informed neural network (PINN) on the least-squares fit residuals. We use a newly defined activation function "DELTAsnake," a modification to the snake activation function proposed by Ziyin et al. that allows for stronger curvature and non-monotonicity. The combined model approximately obeys Maxwell's equations to a level sufficient for producing high quality physics simulations and analysis. Our approach is applied to a highly realistic calculation of the expected magnetic field in the Mu2e experiment's Detector Solenoid which includes a simple model for the expected statistical measurement uncertainties. Using ten toy measurement simulations, we demonstrate the capabilities of our model in comparison to the least-squares method alone; the least-squares method alone results in a reduced chi-squared statistic of ${2.15 \pm 0.01}$, while our approach improves the reduced chi-square to ${1.034 \pm 0.005}$. Furthermore, for an average toy simulation, we show that the range of the RMS of the three field component residuals reduces from ${0.07-0.37}$ Gauss to ${0.05-0.07}$ Gauss. We find that this novel method is robust against a realistic systematic uncertainty deriving from Hall probe calibration bias and can be used to significantly reduce the number of measurements required to achieve an accurate model.

Kampa, Cole [Caltech] (ORCID:0000000192972920)↗
DOE OSTI2016

SU-G-206-06: Analytic Dose Function for CT Scans in Infinite Cylinders as a Function of Scan Length and Cylinder Radius

Purpose: The radiation dose absorbed at a particular radius ρ within the central plane of a long cylinder following a CT scan is a function of the length of the scan L and the cylinder radius R along with kVp and cylinder composition. An analytic function was created that that not only expresses these dependencies but is integrable in closed form over the area of the central plane. This feature facilitates explicit calculation of the planar average dose. The “approach to equilibrium” h(L) discussed in the TG111 report is seamlessly included in this function. Methods: For a cylindrically symmetric radiation field, Monte Carlo calculations were performed to compute the dose distribution to long polyethylene cylinders for scans of varying L for cylinders ranging in radius from 5 to 20 cm. The function was developed from the resultant Monte Carlo data. In addition, the function was successfully fit to data taken from measurements on the 30 cm diameter ICRU/TG200 phantom using a real-time dosimeter. Results: Symmetry and continuity dictate a local extremum at the center which is a minimum for the larger sizes. There are competing effects as the beam penetrates the cylinder from the outside: attenuation, resulting in a decrease; scatter, abruptly increasing at the circumference. This competition may result in an absolute maximum between the center and outer edge leading to a “gull wing” shape for the radial dependence. For the smallest cylinders, scatter may dominate to the extent that there is an absolute maximum at the center. Conclusion: An integrable, analytic function has been developed that provides the radial dependency of dose for the central plane of a scan of length L for cylinders of varying diameter. Equivalently, we have developed h(L,R,ρ).

60 APPLIED LIFE SCIENCES↗
NASA NTRS1974

A three parameter analytic phase function for multiple scattering calculations

A simple procedure was developed to fit the first three moments of an actual phase function with a three parameter analytic phase function. The exact Legendre Polynomial decomposition of this function is, suitable for multiple scattering calculations. The use of this function is expected to yield excellent flux values at all depths within a medium. Since it is capable of reproducing the glory, it can be used in synthetic spectra computations from planetary atmospheres. Accurate asymptotic radiance values can also be achieved as long as the single scattering albedo omega sub 0 is greater than or equal to 0.9.

Kattawar, G. W.↗
NASA NTRS1975

A three-parameter analytic phase function for multiple scattering calculations

A very simple procedure has been developed to fit the first three moments of an actual phase function with a three-parameter analytic phase function. The exact Legendre polynomial decomposition of this function is known which makes it quite suitable for multiple scattering calculations. The use of this function can be expected to yield excellent flux values at all depths within a medium. Since it is capable of reproducing the glory, it can be used in synthetic spectra computations from planetary atmospheres. Accurate asymptotic radiance values can also be achieved as long as the single scattering albedo is not less than 0.9.

Kattawar, G. W.↗
DOE OSTI2019

High Order Anchoring and Reinitialization of Level Set Function for Simulating Interface Motion

A second order interface anchoring method has been developed and used with fast sweeping algorithm for reinitialization of a level set function. The algebraic anchoring formulation ensures that the location of the actual interface is preserved, leading to better mass conservation property. It also provides high order accurate algebraic constraint for solving the Eikonal equation on a finite difference grid. Geometric properties of the interface such as surface normal and curvature are subsequently computed from the reinitialized distance function. Various analytical functions for modeling distortion in level set field are considered and accuracy of reinitialization is evaluated using first and second order anchoring schemes. It is also shown that accurate computation of interface curvature requires a high order anchor in addition to a high order fast sweeping method. Mass conservation property of reinitialization is also analyzed by considering test problems from literature including the classic Rider–Kothe single vortex problem. The formulation is suitable for efficient parallelization for both distributed memory and on-node shared memory parallel systems. In conclusion, scalability and performance of the reinitialization scheme on multiple architectures are demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Generalized Bessel–Struve Operator and Its Properties

We define the generalized Bessel–Struve operator in the space of functions analytic in an arbitrary domain. The conditions of equivalence of the generalized Bessel–Struve operator to the operator of second derivative are investigated. We also describe the commutant of the generalized Bessel–Struve operator and establish its hypercyclicity and chaotic nature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Approximating accelerator impedances with resonator networks

It is common in the accelerator community to use the impedance of accelerator components to describe wake interactions in the frequency domain. However, it is often desirable to understand such wake interactions in the time domain in a general manner for excitations that are not necessarily Gaussian in nature. Here, the conventional method for doing this involves taking the inverse Fourier Transform of the component impedance, obtaining the Green's Function, and then convolving it with the desired excitation distribution. This method can prove numerically cumbersome, for a convolution integral must be evaluated for each individual point in time when the wake function is desired. An alternative to this method would be to compute the wake function analytically, which would sidestep the need for repetitive integration. Only a handful of cases, however, are simple enough for this method to be tenable. One of these cases is the case where the component in question is an RLC resonator, which has a closed-form analytical wake function solution. This means that a component which can be represented in terms of resonators can leverage this solution. As it happens, common network synthesis techniques may be used to map arbitrary impedance profiles to RLC resonator networks in a manner the accelerator community has yet to take advantage of. In this work, we will use Foster Canonical Resonator Networks and partial derivative descent optimization to develop a technique for synthesizing resonator networks that well approximate the impedances of real-world accelerator components. We will link this synthesis to the closed-form resonator wake function solution, giving rise to a powerful workflow that may be used to streamline beam dynamics simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Inverse Thermal Analysis of Ti-6Al-4V Friction Stir Welds Using Numerical-Analytical Basis Functions with Pseudo-Advection

Inverse thermal analysis of Ti-6Al-4V friction stir welds is presented that demonstrates application of a methodology using numerical-analytical basis functions and temperature-field constraint conditions. This analysis provides parametric representation of friction-stir-weld temperature histories that can be adopted as input data to computational procedures for prediction of solid-state phase transformations and mechanical response. These parameterized temperature histories can be used for inverse thermal analysis of friction stir welds having process conditions similar those considered here. Case studies are presented for inverse thermal analysis of friction stir welds that use three-dimensional constraint conditions on calculated temperature fields, which are associated with experimentally measured transformation boundaries and weld-stir-zone cross sections.

36 MATERIALS SCIENCE↗
DOE OSTI2015

SU-E-I-16: Scan Length Dependency of the Radial Dose Distribution in a Long Polyethylene Cylinder

Purpose: The area-averaged dose in the central plane of a long cylinder following a CT scan depends upon the radial dose distribution and the length of the scan. The ICRU/TG200 phantom, a polyethylene cylinder 30 cm in diameter and 60 cm long, was the subject of this study. The purpose was to develop an analytic function that could determine the dose for a scan length L at any point in the central plane of this phantom. Methods: Monte Carlo calculations were performed on a simulated ICRU/TG200 phantom under conditions of cylindrically symmetric conditions of irradiation. Thus, the radial dose distribution function must be an even function that accounts for two competing effects: The direct beam makes its weakest contribution at the center while the scatter begins abruptly at the outer radius and grows as the center is approached. The scatter contribution also increases with scan length with the increase approaching its limiting value at the periphery faster than along the central axis. An analytic function was developed that fit the data and possessed these features. Results: Symmetry and continuity dictate a local extremum at the center which is a minimum for the ICRU/TG200 phantom. The relative depth of the minimum decreases as the scan length grows and an absolute maximum can occur between the center and outer edge of the cylinders. As the scan length grows, the relative dip in the center decreases so that for very long scan lengths, the dose profile is relatively flat. Conclusion: An analytic function characterizes the radial and scan length dependency of dose for long cylindrical phantoms. The function can be integrated with the results expressed in closed form. One use for this is to help determine average dose distribution over the central cylinder plane for any scan length.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2019

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗
NASA NTRS2019

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗
DOE OSTI2017

Analytical trial functions for modelling a two-dimensional Bose condensate

Approximate analytical solutions to the time-independent two-dimensional Gross–Pitaevsky equation are obtained by a variational method. The solutions found are compared with the results of direct numerical calculations. The accuracy of analytical solutions is determined under various conditions of the problem, primarily the degree of nonlinear interaction of atoms in a Bose condensate. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗