Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hybrid Improved Empirical Mode Decomposition and Artificial Neural Network Model for the Prediction of Critical Heat Flux (CHF)

Three Hybrid artificial neural network (ANN) models namely radial basis function (RBF), generalized regression neural networks (GRNN), and multi-layer perceptron (MLP) combined with empirical mode decomposition (EMD) are developed for CHF predictive modelling using CHF experimental databases.First, the original experimental inputs data series are decomposed into several intrinsic mode functions (IMFs) and one residual by EMD, whose components are divided into high, medium and low components. The performance parameters of the hybrid models indicates that the root mean square error (RMSE) are 0.8831, 0.6522, and 0.4149; the mean absolute error (MAE) are 0.6697, 0.4636, and 0.1935. The values of the R-square of the developed prediction approach utilizing EMD-RBF, EMD-GRNN, and EMD-MLP models are 0.8553, 0.9302, and 0.9818, and the index of agreement are 0.9464, 0.9700, and 0.9894., The value of the R-square and the index of agreement of the proposed models are much higher than those of the simple models .The Pearson's test results show that the association strength between the measured and the predicted values of the proposed model EMD-MLP is the strongest. These results show the following: (a) compared with other related, recent studies, the prediction accuracy of the hybrid model EMD- MLP proposed in this research is the best hybrid model; (b) the proposed hybrid model (EMD-MLP) attains superior performance compared with simple models.

Djeddou, Messaoud↗
DOE OSTI2021

Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering

Angularly resolved Thomson scattering is a novel extension of Thomson scattering, enabling the measurement of the electron velocity distribution function over many orders of magnitude. Here, details of the theoretical basis of the technique and the instrument designed for this measurement are described. Angularly resolved Thomson-scattering data from several experiments are shown with descriptions of the corresponding distribution functions. A reduced model describing the distribution function is given and used to perform a Monte Carlo analysis of the uncertainty in the measurements. The electron density and temperature were determined to a precision of 12% and 21%, respectively, on average, while all other parameters defining the distribution function were generally determined to better than 20%. It was found that these uncertainties were primarily due to limited signal to noise and instrumental effects. Furthermore, measurements with this level of precision were sufficient to distinguish between Maxwellian and non-Maxwellian distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

A minimum entropy principle in the gas dynamics equations

Let u(x bar,t) be a weak solution of the Euler equations, governing the inviscid polytropic gas dynamics; in addition, u(x bar, t) is assumed to respect the usual entropy conditions connected with the conservative Euler equations. We show that such entropy solutions of the gas dynamics equations satisfy a minimum entropy principle, namely, that the spatial minimum of their specific entropy, (Ess inf s(u(x,t)))/x, is an increasing function of time. This principle equally applies to discrete approximations of the Euler equations such as the Godunov-type and Lax-Friedrichs schemes. Our derivation of this minimum principle makes use of the fact that there is a family of generalized entrophy functions connected with the conservative Euler equations.

Tadmor, E.↗
NASA NTRS1991

Fuzzy modeling based on generalized neural networks and fuzzy clustering objective functions

An approach to the formulation of fuzzy if-then rules based on clustering objective functions is proposed. The membership functions are then calibrated with the generalized neural networks technique to achieve a desired input-output mapping. The learning procedure is basically a gradient-descent algorithm. A Kalman filter algorithm is used to improve the overall performance.

Sun, Chuen-Tsai↗
NASA NTRS2002

Adjoint-Based, Three-Dimensional Error Prediction and Grid Adaptation

Engineering computational fluid dynamics (CFD) analysis and design applications focus on output functions (e.g., lift, drag). Errors in these output functions are generally unknown and conservatively accurate solutions may be computed. Computable error estimates can offer the possibility to minimize computational work for a prescribed error tolerance. Such an estimate can be computed by solving the flow equations and the linear adjoint problem for the functional of interest. The computational mesh can be modified to minimize the uncertainty of a computed error estimate. This robust mesh-adaptation procedure automatically terminates when the simulation is within a user specified error tolerance. This procedure for estimating and adapting to error in a functional is demonstrated for three-dimensional Euler problems. An adaptive mesh procedure that links to a Computer Aided Design (CAD) surface representation is demonstrated for wing, wing-body, and extruded high lift airfoil configurations. The error estimation and adaptation procedure yielded corrected functions that are as accurate as functions calculated on uniformly refined grids with ten times as many grid points.

Park, Michael A.↗
NASA NTRS2008

Transport Aircraft System Identification from Wind Tunnel Data

Recent studies have been undertaken to investigate and develop aerodynamic models that predict aircraft response in nonlinear unsteady flight regimes for transport configurations. The models retain conventional static and rotary dynamic terms but replace conventional acceleration terms with more general indicial functions. In the Integrated Resilient Aircraft Controls project of the NASA Aviation Safety Program one aspect of the research is to apply these current developments to transport configurations to facilitate development of advanced controls technology. This paper describes initial application of a more general modeling methodology to the NASA Langley Generic Transport Model, a sub-scale flight test vehicle.

Murphy, Patrick C.↗
DOE OSTI2022

Probing for the Trace Estimation of a Permuted Matrix Inverse Corresponding to a Lattice Displacement

We report thatpProbing is a general technique that is used to reduce the variance of the Hutchinson stochastic estimator for the trace of the inverse of a large, sparse matrix A. The variance of the estimator is the sum of the squares of the off-diagonal elements of A -1 . Therefore, this technique computes probing vectors that when used in the estimator annihilate the largest off-diagonal elements. For matrices that display decay of the magnitude of |A$^{-1}_{ij}$| with the graph distance between nodes i and j, this is achieved through graph coloring of increasing powers A k . Equivalently, when a matrix stems from a lattice discretization, it is computationally beneficial to find a distance-k coloring of the lattice. Previously, a hierarchical coloring was proposed so that k can be increased at runtime as needed without discarding previous work. In this work, we study probing for the more general problem of computing the trace of a permutation of A -1 , say PA -1 . The motivation comes from lattice quantum chromodynamics (QCD), where we need to construct “disconnected diagrams” to extract flavor-separated generalized parton functions. In lattice QCD, where the matrix has a four-dimensional toroidal lattice structure, these nonlocal operators correspond to a PA -1 , where P is the permutation relating to some displacement $\vec{p}$ in one or more dimensions. We focus on a single dimension displacement (p), but our methods are general. We show that probing on A k or (PA) k does not annihilate the largest magnitude elements. To resolve this issue, our displacement-based probing works on PA k using a new coloring scheme that works directly on appropriately displaced neighborhoods on the lattice. We prove lower bounds on the number of colors needed and study the effect of this scheme on variance reduction, both theoretically and experimentally on a real-world lattice QCD calculation. We achieve orders of magnitude speedup over the unprobed or the naively probed methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Polarizabilities and Other Properties of the td Muons Molecular Ion

Wavefunctions of Hylleraas type were used earlier to calculate energy levels of muonic systems. Recently, we found in the case of the molecular ions H2+, D2+ and HD+ that it was necessary to include high powers of the internuclear distance in the Hylleraas functions to localize the nuclear motion when treating the ions as three-body systems without invoking the Born-Oppenheimer approximation. We try the same approach in a muonic system, td(mu-). Improved convergence is obtained for J = 0 and 1 states for shorter expansions when we use this type of generalized Hylleraas function, but as the expansion length increases the high powers are no longer useful. We obtain good energy values for the two lowest J = 0 and J = 1 states and compare them with the best earlier calculations. Expectation values are obtained for various operators, the Fermi contact parameters, and the permanent quadrupole moment. The cusp conditions are also calculated. The polarizability of the ground state is then calculated using second-order perturbation theory with intermediate J = 1 pseudostates. It should be possible to measure the polarizability by observing Rydberg states of atoms with td(mu-) acting as the nucleus.

Bhatia, A. K.↗
NASA NTRS1987

H2O line emission from shocked gas

The H2O emission expected from a hot astrophysical plasma containing water is computed to obtain (1) a general cooling function for water, and (2) the individual H2O line intensities in the specific case of the shocked gas region in Orion-KL. It is found that for a shocked molecular region, such as has been previously proposed to account for H2, CO, and O I observations of Orion-KL, there are several hundred H2O lines with fluxes that exceed 10 to the -18th W/sq cm into a 1 arcmin beam. Though the strongest of these generally correspond to strong terrestrial water absorption features, making their detection difficult, future balloon and space experiments should be capable of detecting a large number of water lines. An analytic fit to the total cooling due to water is obtained as a function of temperature, H2 density, and H2O column density. At large optical depth, the result exceeds significantly that obtained from the 'universal cooling function' of Hollenbach and McKee (1979).

Neufeld, David A.↗
NASA NTRS1989

Some elements of a theory of multidimensional complex variables. I - General theory. II - Expansions of analytic functions and application to fluid flows

The paper introduces a new theory of N-dimensional complex variables and analytic functions which, for N greater than 2, is both a direct generalization and a close analog of the theory of ordinary complex variables. The algebra in the present theory is a commutative ring, not a field. Functions of a three-dimensional variable were defined and the definition of the derivative then led to analytic functions.

Martin, E. Dale↗
DOE OSTI2023

Preliminary Monte Carlo and Thermal Hydraulic Analysis using a Hybrid ETF-Corrected-Diffusion Prediction Block

This paper builds upon previous work to accelerate the Picard iteration (PI) method typically applied for coupled Monte Carlo-Thermal hydraulic (MC-TH) solutions. Previously, the use of the generalized transfer functions (GTFs) to predict variation in macroscopic cross sections following a perturbation in TH properties was demonstrated for a subset of simple 3D problems. In addition, the reduced-order transport prediction block relied on the first order perturbation (FOP) method, which was shown to have computational overheads. Recent work replaced the FOP block with a 1-group nodal diffusion solver to eliminate these overheads. While the use of diffusion is desirable for large-scale problems, the new solver introduces significant homogenization error. This work aims to address this issue by using the Jacobian-Free Newton Krylov (JFNK) method to generate a set of super homogenization (SPH) factors to improve the accuracy of the diffusion solution. The SPH factors will be used in conjunction with an improved cross section prediction method – the expanded transfer function (ETF) method – to produce a highly accurate flux prediction for an axial 1D boiling water reactor (BWR) pincell following a large perturbation in moderator density. The ETF-corrected diffusion (ETF-CD) block is shown to be highly accurate for the 1D test case. Future work will investigate the accuracy of the method for a realistic 3D pressurized water reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Charge-transfer satellites and chemical bonding in photoemission and x-ray absorption of SrTi O 3 and rutile Ti O 2 : Experiment and first-principles theory with general application to spectroscopic analysis

First-principles, real-time-cumulant, and Bethe-Salpeter-equation calculations fully capture the detailed satellite structure that occurs in response to the sudden creation of the core hole in both photoemission and x-ray absorption spectra of the transition-metal compounds SrTi O 3 and rutile Ti O 2 . Analysis of the excited-state, real-space charge-density fluctuations betrays the physical nature of these many electron excitations that are shown to reflect the materials’ solid-state electronic structure and chemical bonding. This first-principles development of the cumulant-based core hole spectral function is generally applicable to other systems and should become a standard tool for all similar spectroscopic analysis going beyond the quasiparticle physics of the photoelectric effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Full waveform analysis for high pressure ultrasonic measurement

In this paper we report a data analysis protocol for ultrasonic velocity measurements carried out in a synchrotron x-ray facility based multi-anvil high-pressure apparatus. Synthetic ultrasonic signals for the time period between echoes from the two ends of the sample assembly are created with a few simple parameters as the echo waveform is modeled by the first echo from the anvil-assembly boundary. Each echo is modeled as a delta function with an amplitude and arrival time and convolved with this waveform. The final waveform, fit to the data, enables analysis of signals with overlapped echoes as they are common in small cell assemblies and high pressure.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Accelerated coupled Monte Carlo-Thermal hydraulic calculations using a hybrid GTF-diffusion-based prediction block: first results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections' variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint Eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Accelerated Coupled Monte Carlo-Thermal Hydraulic Calculations using a Hybrid GTF-Diffusion-based Prediction Block: First Results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work [1, 2] investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections’ variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Differentiating densities on smooth manifolds

Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. Furthermore, we highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

On Fast Post-Processing of Global Positioning System Simulator Truth Data and Receiver Measurements and Solutions Data

Post-processing of data, related to a GPS receiver test in a GPS simulator and test facility, is an important step towards qualifying a receiver for space flight. Although the GPS simulator provides all the parameters needed to analyze a simulation, as well as excellent analysis tools on the simulator workstation, post-processing is not a GPS simulator or receiver function alone, and it must be planned as a separate pre-flight test program requirement. A GPS simulator is a critical resource, and it is desirable to move off the pertinent test data from the simulator as soon as a test is completed. The receiver and simulator databases are used to extract the test data files for postprocessing. These files are then usually moved from the simulator and receiver systems to a personal computer (PC) platform, where post-processing is done typically using PC-based commercial software languages and tools. Because of commercial software systems generality their functions are notoriously slow and more than often are the bottleneck even for short duration simulator-based tests. There is a need to do post-processing faster and within an hour after test completion, including all required operations on the simulator and receiver to prepare and move off the post-processing files. This is especially significant in order to use the previous test feedback for the next simulation setup or to run near back-to-back simulation scenarios. Solving the post-processing timing problem is critical for a pre-flight test program success. Towards this goal an approach was developed that allows to speed-up post-processing by an order of a magnitude. It is based on improving the post-processing bottleneck function algorithm using a priory information that is specific to a GPS simulation application and using only the necessary volume of truth data. The presented postprocessing scheme was used in support of a few successful space flight missions carrying GPS receivers.

Kizhner, Semion↗