SEARCH · Search NASA

Results for “Approximation methods”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2019

A posteriori error analysis of fluid–structure interactions: Time dependent error

A posteriori error analysis is a technique to quantify the error in particular simulations of a numerical approximation method. In this article, we use such an approach to analyze how various error components propagate in certain moving boundary problems. We study quasi-steady state simulations where slowly moving boundaries remain in mechanical equilibrium with a surrounding fluid. Such problems can be numerically approximated with the Method of Regularized Stokelets (MRS), a popular method used for studying viscous fluid–structure interactions, especially in biological applications. Our approach to monitoring the regularization error of the MRS is novel, along with the derivation of linearized adjoint equations to the governing equations of the MRS with a elastic elements. Our main numerical results provide a clear illustration of how the error evolves over time in several MRS simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A General Framework for Bounding Approximate Dynamic Programming Schemes

For years, there has been interest in approximation methods for solving dynamic programming problems, because of the inherent complexity in computing optimal solutions characterized by Bellman’s principle of optimality. A wide range of approximate dynamic programming (ADP) methods now exists. It is of great interest to guarantee that the performance of an ADP scheme be at least some known fraction, say ß , of optimal. This letter introduces a general approach to bounding the performance of ADP methods, in this sense, in the stochastic setting. The approach is based on new results for bounding greedy solutions in string optimization problems, where one has to choose a string (ordered set) of actions to maximize an objective function. This bounding technique is inspired by submodularity theory, but submodularity is not required for establishing bounds. Instead, the bounding is based on quantifying certain notions of curvature of string functions; the smaller the curvatures the better the bound. The key insight is that any ADP scheme is a greedy scheme for some surrogate string objective function that coincides in its optimal solution and value with those of the original optimal control problem. The ADP scheme then yields to the bounding technique mentioned above, and the curvatures of the surrogate objective determine the value ß of the bound. The surrogate objective and its curvatures depend on the specific ADP.

discrete event systems↗

DOE OSTI2016

Calculating Time-Integral Quantities in Depletion Calculations

A method referred to as tally nuclides is presented for accurately and efficiently calculating the time-step averages and integrals of any quantities that are weighted sums of atomic densities with constant weights during the step. The method allows all such quantities to be calculated simultaneously as a part of a single depletion solution with existing depletion algorithms. Some examples of the results that can be extracted include step-average atomic densities and macroscopic reaction rates, the total number of fissions during the step, and the amount of energy released during the step. Furthermore, the method should be applicable with several depletion algorithms, and the integrals or averages should be calculated with an accuracy comparable to that reached by the selected algorithm for end-of-step atomic densities. The accuracy of the method is demonstrated in depletion calculations using the Chebyshev rational approximation method. Here, we demonstrate how the ability to calculate energy release in depletion calculations can be used to determine the accuracy of the normalization in a constant-power burnup calculation during the calculation without a need for a reference solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1975

Nozzle Free Jet Flows Within the Strong Curved Shock Regime

A study based on inviscid analysis was conducted to examine the flow field produced from a convergent-divergent nozzle when a strong curved shock occurs. It was found that a certain constraint is imposed on the flow solution of the problem which is the unique feature of the flow within this flow regime, and provides the reason why the inverse method of calculation cannot be employed for these problems. An approximate method was developed to calculate the flow field, and results were obtained for two-dimensional flows. Analysis and calculations were performed for flows with axial symmetry. It is shown that under certain conditions, the vorticity generated at the jet boundary may become infinite and the viscous effect becomes important. Under other conditions, the asymptotic free jet height as well as the corresponding shock geometry were determined.

Shih, Tso-Shin↗

DOE OSTI2018

A new computational approach for the solutions of generalized pantograph-delay differential equations

In this study, a new computational approach is presented to solve the generalized pantograph-delay differential equations (PDDEs). The solutions obtained by our scheme represent by a linear combination of a special kind of basis functions, and can be deduced in a straightforward manner. Firstly, using the least squares approximation method and the Lagrange-multiplier method, the given PDDE is converted to a linear system of algebraic equations, and those unknown coefficients of the solution of the problem are determined by solving this linear system. Secondly, a PDDE related to the error function of the approximate solution is constructed based on the residual error function technique, and error estimation is presented for the suggested method. The convergence of the approximate solution is proved. Several numerical examples are given to demonstrate the accuracy and efficiency. Comparisons are made between the proposed method and other existing methods.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1982

Solution of axisymmetric fluid structure interaction problems with NASTRAN

The solution of axisymmetric acoustic fluid structure interaction problems, employing the NASTRAN computer program is presented. A previously developed 3-D Cartesian Coordinates pressure element formulation is adapted especially for axisymmetric elements. Analogous to the 3-D Cartesian Coordinate predecessor, the fluid portion of the problem is modeled with finite elements wherein one of the displacement components serves as a dummy variable for the pressure unknowns. Two alternatives for implementation of the analogy are presented: (1) an approximate method by which dummy values of G, and nu are used to approximately invoke the analogy wherein the accuracy of the approximation is made as close as desired to the proper analogy within an arbitrary small parameter epsiton; (2) an exact method whereby the NASTRAN FORTRAN coding is slightly changed to invoke the analogy exactly. Comparison of the finite element solution to the exact solution to the same problem is given.

Kalinowski, A. J.↗

NASA NTRS1997

Thermal Radiation Effects Analyzed in Translucent Composite and Thermal Barrier Coating

Ceramic parts and coatings are needed to withstand high temperatures in advanced aircraft engines. In hot environments, such as in the combustion chambers of these engines, infrared and visible radiation can penetrate into some ceramics and heat them internally. The internal temperatures depend on radiative effects combined with heat conduction, and on convection and radiation at the material boundaries. Since engine temperatures are high, radiant emission can be large from within translucent parts and coatings, and this must be included in the analysis. Transient and steady-state behavior are both important. During a transient, radiant penetration provides more rapid internal heating than conduction alone, and the temperature distributions are usually considerably different than for steady-state conditions; this can produce transient thermal stresses. Analytical and numerical methods are being used at the NASA Lewis Research Center to predict transient temperatures and heat flows in translucent materials. A transient analysis was done for a composite of two translucent layers. The layer refractive indices were larger than 1, producing internal reflections at the boundaries and at the internal interface. In addition, steady-state results were computed for a two-layer composite with one layer opaque. The results were used to assess the importance of internal radiation in a zirconia thermal barrier coating on a cooled metal wall. Since the radiative transfer equations are rather complex, especially when internal scattering is included, approximate methods are being investigated that might be convenient for computer design programs. An approximate two-flux method was used and verified by comparisons with solutions that used exact radiative transfer relations. The two-flux equations include scattering without increasing solution difficulty.

Siegel, Robert↗

DOE OSTI2018

Towards the Robust Small-Signal Stability Region of Power Systems Under Perturbations Such as Uncertain and Volatile Wind Generation

This paper addresses how to extend the concept of small-signal stability region (SSSR) to that of robust small-signal stability region (RSSSR), within which the system can remain stable even under perturbations caused by uncertain and volatile nodal injections, such as renewable generation. We first employ the structured perturbation theory to formulate the perturbations of nodal injections in state space. Both the intensity and the locations of perturbations can be taken into account. Then, we leverage the stability radius theory and structured singular value theory to define the RSSSR in parameter subspace, enabling a systematic analysis of small-signal stability of power systems under perturbations in a region-wise manner. The hyperplane-approximation method can be employed to construct a linear closed-form approximation of RSSSR boundaries. Case studies on the modified two-area system and New England 39-node system illustrate the new concept of RSSSR and its potential applications.

power system stability analysis↗

DOE OSTI2021

A Data-Driven Global Sensitivity Analysis Framework for Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.

14 SOLAR ENERGY↗

DOE OSTI2021

Conditional Point Sampling: A Monte Carlo Method for Radiation Transport in Stochastic Media.

Current methods for stochastic media transport are either computationally expensive or, by nature, approximate. Moreover, none of the well-developed, benchmarked approximate methods can compute the variance caused by the stochastic mixing, a quantity especially important to safety calculations. Therefore, we derive and apply a new conditional probability function (CPF) for use in the recently developed stochastic media transport algorithm Conditional Point Sampling (CoPS), which 1) leverages the full intra-particle memory of CoPS to yield errorless computation of stochastic media outputs in 1D, binary, Markovian-mixed media, and 2) leverages the full inter-particle memory of CoPS and the recently developed Embedded Variance Deconvolution method to yield computation of the variance in transport outputs caused by stochastic material mixing. Numerical results demonstrate errorless stochastic media transport as compared to reference benchmark solutions with the new CPF for this class of stochastic mixing as well as the ability to compute the variance caused by the stochastic mixing via CoPS. Using previously derived, non-errorless CPFs, CoPS is further found to be more accurate than the atomic mix approximation, Chord Length Sampling (CLS), and most of memory-enhanced versions of CLS surveyed. In addition, we study the compounding behavior of CPF error as a function of cohort size (where a cohort is a group of histories that share intra-particle memory) and recommend that small cohorts be used when computing the variance in transport outputs caused by stochastic mixing.

61 RADIATION PROTECTION AND DOSIMETRY↗

NASA NTRS2016

Optimization of Turbine Engine Cycle Analysis with Analytic Derivatives

A new engine cycle analysis tool, called Pycycle, was built using the OpenMDAO framework. Pycycle provides analytic derivatives allowing for an efficient use of gradient-based optimization methods on engine cycle models, without requiring the use of finite difference derivative approximation methods. To demonstrate this, a gradient-based design optimization was performed on a turbofan engine model. Results demonstrate very favorable performance compared to an optimization of an identical model using finite-difference approximated derivatives.

OpenMDAO↗

DOE OSTI2016

Numeric Modified Adomian Decomposition Method for Power System Simulations

This paper investigates the applicability of numeric Wazwaz El Sayed modified Adomian Decomposition Method (WES-ADM) for time domain simulation of power systems. WESADM is a numerical method based on a modified Adomian decomposition (ADM) technique. WES-ADM is a numerical approximation method for the solution of nonlinear ordinary differential equations. The non-linear terms in the differential equations are approximated using Adomian polynomials. In this paper WES-ADM is applied to time domain simulations of multimachine power systems. WECC 3-generator, 9-bus system and IEEE 10-generator, 39-bus system have been used to test the applicability of the approach. Several fault scenarios have been tested. It has been found that the proposed approach is faster than the trapezoidal method with comparable accuracy.

Adomian Decomposition↗

NASA NTRS2016

Optimization of Turbine Engine Cycle Analysis with Analytic Derivatives

A new engine cycle analysis tool, called Pycycle, was recently built using the OpenMDAO framework. This tool uses equilibrium chemistry based thermodynamics, and provides analytic derivatives. This allows for stable and efficient use of gradient-based optimization and sensitivity analysis methods on engine cycle models, without requiring the use of finite difference derivative approximation methods. To demonstrate this, a gradient-based design optimization was performed on a multi-point turbofan engine model. Results demonstrate very favorable performance compared to an optimization of an identical model using finite-difference approximated derivatives.

Turbine engine↗

NASA NTRS1975

Research study on stabilization and control: Modern sampled-data control theory. Design of the large space telescope system

Conditions of self-sustained oscillations in a two-axis model of the nonlinear LST system are studied. The describing function of the CMG frictional nonlinearity of the LST system is used for the analysis, as well as continuous-data and discrete-data models of the simplified LST control system. A numerical-iterative method is described for the analysis of the two-axis system. Approximation methods and the direct plotting of the stability equation are implemented in the study. It is shown that although the dynamics of the two axes are identical, the amplitudes of self-sustained oscillations in the two axes may in principle be different. Analysis shows that the LST systems are of equal amplitudes but with 180-degree phase shift.

Kuo, B. C.↗

NASA NTRS1986

On compactness of admissible parameter sets: Convergence and stability in inverse problems for distributed parameter systems

A series of numerical examples is reported and several algorithms compared for estimation of coefficients in differential equation models. Unconstrained, constrained and Tikhonov regularization methods are tested for their behavior with regard to both convergence (of approximation methods for the states and parameters) and stability (continuity of the estimates with respect to perturbations in the data or observed states).

Banks, H. T.↗

NASA NTRS1997

Optimizing an Actuator Array for the Control of Multi-Frequency Noise in Aircraft Interiors

Techniques developed for selecting an optimized actuator array for interior noise reduction at a single frequency are extended to the multi-frequency case. Transfer functions for 64 actuators were obtained at 5 frequencies from ground testing the rear section of a fully trimmed DC-9 fuselage. A single loudspeaker facing the left side of the aircraft was the primary source. A combinatorial search procedure (tabu search) was employed to find optimum actuator subsets of from 2 to 16 actuators. Noise reduction predictions derived from the transfer functions were used as a basis for evaluating actuator subsets during optimization. Results indicate that it is necessary to constrain actuator forces during optimization. Unconstrained optimizations selected actuators which require unrealistically large forces. Two methods of constraint are evaluated. It is shown that a fast, but approximate, method yields results equivalent to an accurate, but computationally expensive, method.

Palumbo, D. L.↗

NASA NTRS1958

Stability-Derivative Determination from Flight Data

A comprehensive discussion of the various factors affecting the determination of stability and control derivatives from flight data is presented based on the experience of the NASA High-Speed Flight Station. Factors relating to test techniques, determination of mass characteristics, instrumentation, and methods of analysis are discussed. For most longitudinal-stability-derivative analyses simple equations utilizing period and damping have been found to be as satisfactory as more comprehensive methods. The graphical time-vector method has been the basis of lateral-derivative analysis, although simple approximate methods can be useful If applied with caution. Control effectiveness has been generally obtained by relating the peak acceleration to the rapid control input, and consideration must be given to aerodynamic contributions if reasonable accuracy is to be realized.. Because of the many factors involved In the determination of stability derivatives, It is believed that the primary stability and control derivatives are probably accurate to within 10 to 25 percent, depending upon the specific derivative. Static-stability derivatives at low angle of attack show the greatest accuracy.

Holowicz, Chester H.↗

NASA NTRSDate not supplied

An Overview of the Patch Integral Method (PIM), a New Heat Transfer Analysis Tool for Hypersonic Wind Tunnel Facilities at NASA Langley

NASA Langley’s hypersonic wind tunnels are heavily leveraged for planetary missions. The data collection method in these tunnels is thermography, and surface temperature measurements of the model surface are collected and reduced to produce surface heating data, as seen in Fig 1. However, during model injection, no temperature data are collected, and thus conventional, integral heat transfer methods cannot be used to solve for surface heating. A method was developed in the 1990’s to reduce this data despite the data gap, known as the step approximation method. The method assumes that the film coefficient behaves as a step function, the model is semi-infinite, and thermal properties are constant. With these simplifying assumptions, a Laplace transform can be performed to result in an equation that takes the initial temperature of the model and a temperature at some point in time to back out the film coefficient at that time. This is the method that is used in the current thermographic data reduction software, IHEAT. While computationally light-weight, the step approximation has several issues associated with it. The time-history of temperature is not accounted for, which is vital as heat transfer is an integral process. Additionally, the required semi-infinite assumption is unnecessary and might be violated during runtime. Thermal variation of material properties can have a sizeable impact on heating results and are not modeled by the method. This method also takes multiple seconds to “collapse” to a steady state, which is undesirable from both a facility and data reduction standpoint. The method is also very sensitive to the “effective time” approximation, an approximation of when heating instantaneously starts (which is a nonphysical simplification), and a small variation in this value can result in an error in heating results.

J. S. Cheatwood↗