Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximation & analysis”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

BLSTA: A boundary layer code for stability analysis

A computer program is developed to solve the compressible, laminar boundary-layer equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows including the flow along the plane of symmetry, flow along the leading-edge attachment line, and swept-wing flows with a conical flow approximation. The finite-difference numerical procedure used to solve the governing equations is second-order accurate. The flow over a wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be calculated. Various wall boundary conditions, such as wall suction or blowing and hot or cold walls, can be applied. The results indicate that this boundary-layer code gives velocity and temperature profiles which are accurate, smooth, and continuous through the first and second normal derivatives. The code presented herein can be coupled with a stability analysis code and used to predict the onset of the boundary-layer transition which enables the assessment of the laminar flow control techniques. A user's manual is also included.

Wie, Yong-Sun↗
NASA NTRS1980

Grid generation using differential systems techniques

The errors in approximating the derivatives of a function by traditional central differences at grid points of a curvilinear coordinate system were examined. The implications concerning the accuracy of the numerical solution of a partial differential equation are explained by considering several numerical examples. Although this study only considers the two dimensional case, the techniques and implications are equally valid for three dimensional grids. An interesting feature of the error analysis is its simplicity. Most of the results follow by merely working with the truncation terms of some power series expansion. These series expansions also give rise to higher order difference approximations which can significantly reduce error when the grid spacing changes rapidly, as might be the case in problems with shock waves or thin boundary layers.

Thompson, J. F.↗
DOE OSTI2020

An investigation of Newton-Sketch and subsampled Newton methods

Sketching, a dimensionality reduction technique, has received much attention in the statistics community. In this paper, we study sketching in the context of Newton's method for solving finite-sum optimization problems in which the number of variables and data points are both large. In this work, we study two forms of sketching that perform dimensionality reduction in data space: Hessian subsampling and randomized Hadamard transformations. Each has its own advantages, and their relative tradeoffs have not been investigated in the optimization literature. Additionally, our study focuses on practical versions of the two methods in which the resulting linear systems of equations are solved approximately, at every iteration, using an iterative solver. The advantages of using the conjugate gradient method vs. a stochastic gradient iteration are revealed through a set of numerical experiments, and a complexity analysis of the Hessian subsampling method is presented.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Shock waves and drag in the numerical calculation of isentropic transonic flow

Properties of the shock relations for steady, irrotational, transonic flow are discussed and compared for the full and approximate governing potential in common use. Results from numerical experiments are presented to show that the use of proper finite difference schemes provide realistic solutions and do not introduce spurious shock waves. Analysis also shows that realistic drags can be computed from shock waves that occur in isentropic flow. In analogy to the Oswatitsch drag equation, which relates the drag to entropy production in shock waves, a formula is derived for isentropic flow that relates drag to the momentum gain through an isentropic shock. A more accurate formula for drag, based on entropy production, is also derived, and examples of wave drag evaluation based on these formulas are given and comparisons are made with experimental results.

Steger, J. L.↗
DOE OSTI2023

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗
NASA NTRS1974

Analysis of the effects of baffles on combustion instability

An analytical model has been developed for predicting the effects of baffles on combustion instability. This model has been developed by coupling an acoustic analysis of the wave motion within baffled chambers with a model for the oscillatory combustion response of a propellant droplet developed by Heidmann. A computer program was developed for numerical solution of the resultant coupled equations. Diagnostic calculations were made to determine the reasons for the improper prediction. These calculations showed that the chosen method of representing the combustion response was a very poor approximation. At the end of the program, attempts were made to minimize this effect but the model still improperly predicts the stability trends. Therefore, it is recommended that additional analysis be done with an improved approximation.

Oberg, C. L.↗
DOE OSTI2020

Quantum approximate optimization of the long-range Ising model with a trapped-ion quantum simulator

Quantum computers and simulators may offer significant advantages over their classical counterparts, providing insights into quantum many-body systems and possibly improving performance for solving exponentially hard problems, such as optimization and satisfiability. Here, we report the implementation of a low-depth Quantum Approximate Optimization Algorithm (QAOA) using an analog quantum simulator. We estimate the ground-state energy of the Transverse Field Ising Model with long-range interactions with tunable range, and we optimize the corresponding combinatorial classical problem by sampling the QAOA output with high-fidelity, single-shot, individual qubit measurements. We execute the algorithm with both an exhaustive search and closed-loop optimization of the variational parameters, approximating the ground-state energy with up to 40 trapped-ion qubits. We benchmark the experiment with bootstrapping heuristic methods scaling polynomially with the system size. We observe, in agreement with numerics, that the QAOA performance does not degrade significantly as we scale up the system size and that the runtime is approximately independent from the number of qubits. We finally give a comprehensive analysis of the errors occurring in our system, a crucial step in the path forward toward the application of the QAOA to more general problem instances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Heat transfer distribution on a discontinuously convective cooled surface

An active cooling scheme in which discrete cooling passages are used was investigated. In this scheme, the coolant passages are placed perpendicular to the main stream direction. Streamwise temperature oscillations result from the highly cooled regions, near the coolant channels, being next to regions not directly cooled between the channels. The accuracy of the heat transfer rate predictions by these approximations was evaluated by means of both experimental and numerical investigations. Experiments were conducted at free stream Mach number of 6 with several cooling passages. Heat transfer measurements were made and corresponding theoretical analysis were performed for comparison in evaluating both local and overall deviations in the heat transfer rate.

Piva, R.↗
NASA NTRSDate not supplied

Sensitivity of Simulated Radiative Emission from Nitrogen Flows to Chemical Parameters

Numerical models of hypersonic flows rely on chemistry data inferred from experiments and quantum calculations. Experiments from NASA's Electric Arc Shock Tube (EAST) measure the spectral emission of strong shocks in pure nitrogen flows. These spectral profiles are only approximately captured by numerical simulations. To address these discrepancies, one wishes to calibrate some of the chemical parameters to better match experiment. Because there are hundreds of parameters in the numerical model, a sensitivity analysis was carried out on the entire chemistry database to determine which parameters are most important to analyze in the future. In order to perform this work it was necessary to generate a consistent chemical database that interfaced with both the CFD and line-by-line radiation solver. Using Monte Carlo methods, the total Sobol index was calculated for important parameters. Twenty parameters were identified as important across the nitrogen system. However, it was noted that in some spectral regions, the parametric uncertainty was unable to bound EAST measurements. Further improvements could be made by increasing the parameter uncertainties or addressing other sources of modeling error.

Kaelan B Hansson↗
NASA NTRSDate not supplied

Sensitivity of Simulated Radiative Emission from Nitrogen Flows to Chemical Parameters

Numerical models of hypersonic flows rely on chemistry data inferred from experiments and quantum calculations. Experiments from NASA's Electric Arc Shock Tube (EAST) measure the spectral emission of strong shocks in pure nitrogen flows. These spectral profiles are only approximately captured by numerical simulations. To address these discrepancies, one wishes to calibrate some of the chemical parameters to better match experiment. Because there are hundreds of parameters in the numerical model, a sensitivity analysis was carried out on the entire chemistry database to determine which parameters are most important to analyze in the future. In order to perform this work it was necessary to generate a consistent chemical database that interfaced with both the CFD and line-by-line radiation solver. Using Monte Carlo methods, the total Sobol index was calculated for important parameters. Twenty parameters were identified as important across the nitrogen system. However, it was noted that in some spectral regions, the parametric uncertainty was unable to bound EAST measurements. Further improvements could be made by increasing the parameter uncertainties or addressing other sources of modeling error.

Kaelan Hansson↗
DOE OSTI2020

Efficient and Robust Dynamic Simulation of Power Systems With Holomorphic Embedding

Dynamic simulation is vitally important in power system analysis, but traditional approaches based on numerical integration over small time steps are time-consuming. Also, the Newton-Raphson method suffers from difficulty in convergence when solving nonlinear algebraic equations. In this paper, we propose a novel dynamic simulation approach based on holomorphic embedding. By obtaining a high-order approximation of system dynamics, it achieves a much larger time step and thus enhances the computational efficiency significantly. In addition, the new approach avoids non-convergence issues in solving algebraic equations, which improves robustness. The approach includes flexible modeling of synchronous generators and controllers, and we propose a method for modeling generator coordinate transformations. The approach is tested on the IEEE 39-bus, 10-generator system and a Polish 2383-bus, 327-generator system. The results demonstrate promising computational efficiency and satisfactory numerical robustness for the analysis of large-scale power systems.

42 ENGINEERING↗
DOE OSTI2021

High-Accuracy Mesh-Free Quadrature for Trimmed Parametric Surfaces and Volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. Additionally, the scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The Dark Energy Survey supernova program: a reanalysis of cosmology results and evidence for evolving dark energy with an updated Type Ia supernova calibration

We present improved cosmological constraints from a re-analysis of the Dark Energy Survey (DES) 5-year sample of Type Ia supernovae (DES-SN5YR). This re-analysis includes an improved photometric cross-calibration, recent white dwarf observations to cross-calibrate between DES and low-redshift surveys, retraining the salt3 light-curve model and fixing a numerical approximation in the host-galaxy colour law. Our fully recalibrated sample, which we call DES-Dovekie, comprises ~1600 likely Type Ia SNe from DES and ~200 low-redshift SNe from other surveys. With DES-Dovekie, we obtain Ω m = 0.330 ± 0.015 in flat Lambda-cold dark matter (⁠ΛCDM) which changes Ω m by –0.022 compared to DES-SN5YR. Combining DES-Dovekie with cosmic microwave background data from Planck, Atacama Cosmology Telescope, and South Pole Telescope and the DESI DR2 measurements in a flat CDM cosmology, we find ω 0 = –0.803 ± 0.054 and ω a = –0.72 ± 0.21⁠. Our results hold a significance of 3.2σ, reduced from 4.2σ for DES-SN5YR, to reject the null hypothesis that the data are compatible with the cosmological constant. This significance is equivalent to a Bayesian model preference odds of approximately 5:1 in favour of the flat ω 0 ω a CDM model. Using generally accepted thresholds for model preference, our updated data exhibits only a weak preference for evolving dark energy.

dark energy↗
NASA NTRS1974

Analytical and experimental study of mean flow and turbulence characteristics inside the passages of an axial flow inducer

The effort conducted to gather additional understanding of the complex inviscid and viscid effects existing within the passages of a three-bladed axial flow inducer operating at a flow coefficient of 0.065 is summarized. The experimental investigations included determination of the blade static pressure and blade limiting streamline angle distributions, and measurement of the three components of mean velocity, turbulence intensities and turbulence stresses at locations inside the inducer blade passage utilizing a rotating three-sensor hotwire probe. Applicable equations were derived for the hotwire data reduction analysis and solved numerically to obtain the appropriate flow parameters. Analytical investigations were conducted to predict the three-dimensional inviscid flow in the inducer by numerically solving the exact equations of motion, and to approximately predict the three-dimensional viscid flow by incorporating the dominant viscous terms into the exact equations. The analytical results are compared with the experimental measurements and design values where appropriate.

Gorton, C. A.↗
NASA NTRS1995

Computational analysis of forebody tangential slot blowing on the high alpha research vehicle

A numerical analysis of forebody tangential slot blowing as a means of generating side force and yawing moment is conducted using an aircraft geometry. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved using a partially flux-split, approximately-factored algorithm. An algebraic turbulence model is used to determine the turbulent eddy viscosity values. Solutions are obtained using both patched and overset grid systems. In the patched grid model, and actuator plane is used to introduce jet variables into the flow field. The overset grid model is used to model the physical slot geometry and facilitate modeling of the full aircraft configuration. A slot optimization study indicates that a short slot located close to the nose of the aircraft provided the most side force and yawing moment per unit blowing coefficient. Comparison of computed surface pressure with that obtained in full-scale wind tunnel tests produce good agreement, indicating the numerical method and grid system used in the study are valid. Full aircraft computations resolve the changes in vortex burst point due to blowing. A time-accurate full-aircraft solution shows the effect of blowing on the changes in the frequency of the aerodynamic loads over the vertical tails. A study of the effects of freestream Mach number and various jet parameters indicates blowing remains effective through the transonic Mach range. An investigation of the force onset time lag associated with forebody blowing shows the lag to be minimal. The knowledge obtained in this study may be applied to the design of a forebody tangential slot blowing system for use on flight aircraft.

Gee, Ken↗
DOE OSTI2023

A unified and constructive framework for the universality of neural networks

Abstract One of the reasons why many neural networks are capable of replicating complicated tasks or functions is their universal approximation property. Though the past few decades have seen tremendous advances in theories of neural networks, a single constructive and elementary framework for neural network universality remains unavailable. This paper is an effort to provide a unified and constructive framework for the universality of a large class of activation functions including most of the existing ones. At the heart of the framework is the concept of neural network approximate identity (nAI). The main result is as follows: any nAI activation function is universal in the space of continuous functions on compacta. It turns out that most of the existing activation functions are nAI, and thus universal. The framework induces several advantages over the contemporary counterparts. First, it is constructive with elementary means from functional analysis, probability theory, and numerical analysis. Second, it is one of the first unified and constructive attempts that is valid for most of the existing activation functions. Third, it provides new proofs for most activation functions. Fourth, for a given activation and error tolerance, the framework provides precisely the architecture of the corresponding one-hidden neural network with a predetermined number of neurons and the values of weights/biases. Fifth, the framework allows us to abstractly present the first universal approximation with a favorable non-asymptotic rate. Sixth, our framework also provides insights into the developments, and hence providing constructive derivations, of some of the existing approaches.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2023

Equivalent Area Targets for Inverse Design Optimization With Changing Low-Boom Cruise Conditions

This paper provides theoretical and numerical justifications for inverse design optimization of reversed equivalent area (A e,r ) of low-boom supersonic aircraft when the cruise condition changes during the optimization iterations. A modified linear theory for steady flow around a supersonic projectile is used to establish the accuracy of A e,r -based body-of-revolution approximation of a low-boom supersonic aircraft for undertrack sonic boom analysis using computational fluid dynamics off-body pressure. As a result, designing a low-boom shape of the off-body pressure at three body lengths below the aircraft is equivalent to designing a low-boom A e,r shape. A Bezier curve with eight control points is used to define an A e,r target for inverse design optimization of A e,r of a supersonic aircraft. To make an A e,r target matchable by a low-boom supersonic aircraft, the A e,r target must have the minimal perceived level of decibels (PLdB) for sonic boom and satisfy two constraints defined by the cruise condition, effective length, and two nondimensional parameters. Numerical results are used to verify that, for different cruise conditions, all optimized A e,r targets with PLdB below 70 can be approximately generated by scaling one of them if the targets have the same nondimensional parameters and approximately the same effective length.

multidisciplinary optimization↗
NASA NTRS1980

Numerical simulation of supersonic inlets using a three-dimensional viscous flow analysis

A three-dimensional fully viscous computer analysis, which retains the viscous nature of the Navier-Stokes equations, was evaluated to determine its usefulness in the design of supersonic inlets. This procedure takes advantage of physical approximations to limit the high computer time and storage associated with complete Navier-Stokes solutions. Computed results are presented for a Mach 3.0 supersonic inlet with bleed and a Mach 7.4 hypersonic inlet. Good agreement was obtained between theory and data for both inlets. Results of a mesh sensitivity study are also shown.

Anderson, B. H.↗