Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function approximation”

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 1,081 records · Page 60

Atomic Data and Spectral Line Intensities for Ni XI

Electron impact collision strengths, energy levels, oscillator strengths and spontaneous radiative decay rates are calculated for Ni XI. We include in the calculations the 10 lowest configurations, corresponding to 164 fine structure levels: 3s(sup 2)3p(sup 6), 3s(sup 2)3p(sup 5)3d, 3s(sup 2)3p(sup 4)3d(sup 2), 3s3p(sup 6)3d, 3s(sup 2)3p(sup 5)4l and 3s3p6 4l with l =.s, p, d. Collision strengths are calculated at five incident energies for all transitions: 7.1, 16.8, 30.2, 48.7 and 74.1 Ry above the threshold of each transition. An additional energy, very close to the transition threshold, has been added, whose value is between 0.06 Ry and 0.25 Ry depending on the lower level. Calculations have been carried out using the Flexible Atomic Code and the distorted wave approximation. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates of the present work, combined with Close Coupling collision excitation rate coefficient available in the literature for the lowest 17 levels, statistical equilibrium equations for level populations are solved at electron densities covering the 10(exp 8)-10(exp 14) cu cm range and at an electron temperature of logT(sub c)(K)=6.1, corresponding to the maximum abundance of Ni XI. Spectral line intensities are calculated, and their diagnostic relevance is discussed. This dataset will be made available in the next version of the CHIANTI database.

Bhatia, A. K.↗

Position Sensor Integral with a Linear Actuator

A noncontact position sensor has been designed for use with a specific two-dimensional linear electromagnetic actuator. To minimize the bulk and weight added by the sensor, the sensor has been made an integral part of the actuator: that is to say, parts of the actuator structure and circuitry are used for sensing as well as for varying position. The actuator (see Figure 1) includes a C-shaped permanent magnet and an armature that is approximately centered in the magnet gap. The intended function of the actuator is to cause the permanent magnet to translate to, and/or remain at, commanded x and y coordinates, relative to the armature. In addition, some incidental relative motion along the z axis is tolerated but not controlled. The sensor is required to measure the x and y displacements from a nominal central position and to be relatively insensitive to z displacement. The armature contains two sets of electromagnet windings oriented perpendicularly to each other and electrically excited in such a manner as to generate forces in the x,y plane to produce the required motion. Small sensor excitation coils are mounted on the pole tips of the permanent magnet. These coils are excited with a sine wave at a frequency of 20 kHz. This excitation is transformer-coupled to the armature windings. The geometric arrangement of the excitation coils and armature windings is such that the amplitudes of the 20-kHz voltages induced in the armature windings vary nearly linearly with x and y displacements and do not vary significantly with small z displacements. Because the frequency of 20 kHz is much greater than the maximum frequency characteristic of the actuation signals applied to the armature windings, there is no appreciable interference between actuator and sensor functions of the armature windings.

Howard, David E.↗

Analysis and Sizing for Transient Thermal Heating of Insulated Aerospace Vehicle Structures

An analytical solution was derived for the transient response of an insulated structure subjected to a simplified heat pulse. The solution is solely a function of two nondimensional parameters. Simpler functions of these two parameters were developed to approximate the maximum structural temperature over a wide range of parameter values. Techniques were developed to choose constant, effective thermal properties to represent the relevant temperature and pressure-dependent properties for the insulator and structure. A technique was also developed to map a time-varying surface temperature history to an equivalent square heat pulse. Equations were also developed for the minimum mass required to maintain the inner, unheated surface below a specified temperature. In the course of the derivation, two figures of merit were identified. Required insulation masses calculated using the approximate equation were shown to typically agree with finite element results within 10%-20% over the relevant range of parameters studied.

Blosser, Max L.↗

Atomic Data and Spectral Line Intensities for Ca IX

Electron impact collision strengths, energy levels, oscillator strengths and spontaneous radiative decay rates are calculated for Ca IX. We include in the calculations the 33 lowest configurations in the n = 3, 4, 5 complexes, corresponding to 283 fine structure levels in the 3l3l ', 3l4l'' and 3l4l''' configurations, where l,l' = s, p, d, l '' = s, p, d, f and l''' = s, p, d, f, g. Collision strengths are calculated at five incident energies for all transitions: 5.8, 13.6, 24.2, 38.6 and 57.9 Ry above the threshold of each transition. An additional energy, very close to the transition threshold, has been added, whose value is between 0.0055 Ry and 0.23 Ry depending on the levels involved. Calculations have been carried out using the Flexible Atomic Code and the distorted wave approximation. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates calculated in the present work, statistical equilibrium equations for level populations are solved at electron densities covering the 10(exp 8)-10(exp 14)/cubic cm range and at an electron temperature of log T(sub e)(K)=5.8, corresponding to the maximum abundance of Ca IX. Spectral line intensities are calculated, and their diagnostic relevance is discussed.

Landi, E.↗

Atomic Data and Spectral Line Intensities for Ni XV

Electron impact collision strengths, energy levels, oscillator strengths, and spontaneous radiative decay rates are calculated for Ni XV.Weinclude in the calculations the 9 lowest configurations, corresponding to 126 fine structure levels: 3s23p2, 3s3p3, 3s23p3d, 3p4, 3s3p23d, and 3s2 3p4l with l =, s, p, d, f. Collision strengths are calculated at five incident energies for all transitions: 7.8, 18.5, 33.5, 53.5, and 80.2 Ry above the threshold of each transition. An additional energy, very close to the transition threshold, has been added, whose value is between 0.004 and 0.28 Ry depending on the levels involved. Calculations have been carried out using the Flexible Atomic Code and the distorted-wave approximation. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. Using the excitation rate coefficients and the radiative transition rates calculated in the present work, statistical equilibrium equations for level populations are solved at electron densities covering the 10(exp 8)-10(exp 14)/cu cm range and at an electron temperature of log T(sub e)(K) = 6.4, corresponding to the maximum abundance of Ni XV. Spectral line intensities are calculated, and their diagnostic relevance is discussed. This dataset will be made available in the next version of the CHIANTI database.

Landi, E.↗

Bayesian Revision of Residual Detection Power

This paper addresses some issues with quality assessment and quality assurance in response surface modeling experiments executed in wind tunnels. The role of data volume on quality assurance for response surface models is reviewed. Specific wind tunnel response surface modeling experiments are considered for which apparent discrepancies exist between fit quality expectations based on implemented quality assurance tactics, and the actual fit quality achieved in those experiments. These discrepancies are resolved by using Bayesian inference to account for certain imperfections in the assessment methodology. Estimates of the fraction of out-of-tolerance model predictions based on traditional frequentist methods are revised to account for uncertainty in the residual assessment process. The number of sites in the design space for which residuals are out of tolerance is seen to exceed the number of sites where the model actually fails to fit the data. A method is presented to estimate how much of the design space in inadequately modeled by low-order polynomial approximations to the true but unknown underlying response function.

DeLoach, Richard↗

An Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. II: Approach Based on the Dyson and Bethe-Salpeter Equations

In this paper, the vector radiative transfer equation is derived by means of the vector integral Foldy equations describing the electromagnetic scattering by a group of particles. By assuming that in a discrete random medium the positions of the particles are statistically independent and by applying the Twersky approximation to the order-of-scattering expansion of the total field, we derive the Dyson equation for the coherent field and the ladder approximated Bethe–Salpeter equation for the dyadic correlation function. Then, under the far-field assumption for sparsely distributed particles, the Dyson equation is reduced to the Foldy integral equation for the coherent field, while the iterated solution of the Bethe–Salpeter equation ultimately yields the vector radiative transfer equation.

Electromagnetic scattering↗

Adaptive Optimization for System Performance and Combined Bernstein Polynomial, Optimal Reciprocal Collision Avoidance, Differential Dynamic Programming for Trajectory Replanning and Collision Avoidance for UAM Vehicles

The emerging urban air mobility (UAM) sector in aerospace is driving development of unconventional multi-modal vehicle configurations and autonomous flight. The combination of multi-modal vehicle dynamics, complex environment, requirements to deal with flight contingencies in an efficient and safe manner, as well as necessity for precise trajectory following and performance, are the driving influence behind adaptive optimization for system performance. We are interested in trajectory optimization algorithm that would system parameter estimation and identifying the optimal switching time between modes of hybrid dynamical systems. This presentation discusses a parameterized optimal control trajectory optimization algorithm that is an extended and generalized version of Differential Dynamic Programming (DDP), titled Parameterized Differential Dynamic Programming (PDDP). DDP is an efficient trajectory optimization algorithm relying on second order approximations of a system’s dynamics and cost function and has recently been applied to optimize systems with time invariant parameters. Experiments are presented applying PDDP to solve model predictive control (MPC) and moving horizon estimation (MHE) tasks simultaneously. In particular, PDDP is used to determine the optimal transition point between flight regimes of a complex urban air mobility (UAM) class vehicle exhibiting multiple phases of flight and to identify and compensate for actuation faults.

optimization↗

Time-varying Weights in Multi-Objective Optimal Control for Flexible Wing Aircraft

A multi-objective optimal control technique is modified to accommodate changing cost function weights and is used to control a flexible wing aircraft model. Variation of the weights is used to adjust the relative importance of each objective according to either a prescribed function of time or of the state. Several techniques for obtaining a practical approximation to the optimal control solution are presented, and stability of a specific weight structure with the optimal controller is demonstrated. Functionality of the multi-objective control design with weight variation is demonstrated in simulation of a flexible wing transport aircraft and is shown to improve performance over the fixed weight version both at a constant flight condition and across changing flight conditions.

Hashemi, Kelley E.↗

An algorithm for the rapid location of an extreme of a function subject only to geometric restrictions

The requirements of symmetry and convexity in the application of algorithms for the minimization or maximization of a function are discussed. It is argued that if a function of a single variable is convex and symmetric in a neighborhood of an extremum, the extremum may be approximated to the precision that increases by at least a power of two per functional evaluation. The procedure may be used to drive a complex optimization procedure in the multivariate area estimation problem encountered in remote sensing.

Terrell, G. R.↗

Static shape determination and control for large space structures. I - The flexible beam. II - A large space antenna

A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.

Weeks, C. J.↗

Femtolensing: Beyond the semiclassical approximation

Femtolensoing is a gravitational lensing effect in which the magnification is a function not only of the position and sizes of the source and lens, but also of the wavelength of light. Femtolensing is the only known effect of 10(exp -13) - 10(exp -16) solar mass) dark-matter objects and may possibly be detectable in cosmological gamma-ray burst spectra. We present a new and efficient algorithm for femtolensing calculation in general potentials. The physical optics results presented here differ at low frequencies from the semiclassical approximation, in which the flux is attributed to a finite number of mutually coherent images. At higher frequencies, our results agree well with the semicalssical predictions. Applying our method to a point-mass lens with external shear, we find complex events that have structure at both large and small spectral resolution. In this way, we show that femtolensing may be observable for lenses up to 10(exp -11) solar mass, much larger than previously believed. Additionally, we discuss the possibility of a search femtolensing of white dwarfs in the Large Magellanic Cloud at optical wavelengths.

Ulmer, Andrew↗

Imaginary-frequency polarizability and van der Waals force constants of two-electron atoms, with rigorous bounds

Variational functionals of Braunn and Rebane (1972) for the imagery-frequency polarizability (IFP) have been generalized by the method of Gramian inequalities to give rigorous upper and lower bounds, valid even when the true (but unknown) unperturbed wavefunction must be represented by a variational approximation. Using these formulas in conjunction with flexible variational trial functions, tight error bounds are computed for the IFP and the associated two- and three-body van der Waals interaction constants of the ground 1(1S) and metastable 2(1,3S) states of He and Li(+). These bounds generally establish the ground-state properties to within a fraction of a per cent and metastable properties to within a few per cent, permitting a comparative assessment of competing theoretical methods at this level of accuracy. Unlike previous 'error bounds' for these properties, the present results have a completely a priori theoretical character, with no empirical input data.

Glover, R. M.↗

First and second order convex approximation strategies in structural optimization

In this paper, various methods based on convex approximation schemes are discussed that have demonstrated strong potential for efficient solution of structural optimization problems. First, the convex linearization method (Conlin) is briefly described, as well as one of its recent generalizations, the method of moving asymptotes (MMA). Both Conlin and MMA can be interpreted as first-order convex approximation methods that attempt to estimate the curvature of the problem functions on the basis of semiempirical rules. Attention is next directed toward methods that use diagonal second derivatives in order to provide a sound basis for building up high-quality explicit approximations of the behavior constraints. In particular, it is shown how second-order information can be effectively used without demanding a prohibitive computational cost. Various first-order and second-order approaches are compared by applying them to simple problems that have a closed form solution.

Fleury, C.↗

Thermophysical Properties of Liquid Te: Density, Electrical Conductivity, and Viscosity

The thermophysical properties of liquid Te, namely, density, electrical conductivity, and viscosity, were determined using the pycnometric and transient torque methods from the melting point of Te (723 K) to approximately 1150 K. A maximum was observed in the density of liquid Te as the temperature was increased. The electrical conductivity of liquid Te increased to a constant value of 2.89 x 10(exp 5 OMEGA-1m-1) as the temperature was raised above 1000 K. The viscosity decreased rapidly upon heating the liquid to elevated temperatures. The anomalous behaviors of the measured properties are explained as caused by the structural transitions in the liquid and discussed in terms of Eyring's and Bachiskii's predicted behaviors for homogeneous liquids. The Properties were also measured as a function of time after the liquid was coded from approximately 1173 or 1123 to 823 K. No relaxation phenomena were observed in the properties after the temperature of liquid Te was decreased to 823 K, in contrast to the relaxation behavior observed for some of the Te compounds.

Li, C.↗

The hard X-ray emission spectra from accretion columns in intermediate polars

We consider the hard (greater than 2 keV) X-ray emission from accretion columns in an intermediate polar system, GK Per, using a simple settling solution. The rate of photon emission per logarithmic energy interval can be fitted with a power law, E(exp -gamma), with gamma approximately 2.0, in agreement with observations. This index is only weakly dependent on the mass accretion rate, dot-M, for dot-M in the range of a few times 10(exp 16-18) g/s. The peak energy of the photon spectra (after photoelectric absorption) is expected to be E(sub p) approximately (5 keV) gamma(exp -1/3) (N(sub H)/10(exp 23)/sq cm)(exp 1/3) where N(sub H) is the hydrogen column density along the line of sight. The observed spectra of GK Per and possibly of V1223 Sgr suggest N(sub H) approximately 10(exp 23)/sq cm. This large N(sub H) may be due to partially ionized preshock column material. Alternatively, we also consider absorption by the cool outer parts of an accretion disk. In this case the photoelectric absorption depth in the disk is a sensitive function of inclination. For GK Per the required inclination is approximately 83 deg. For mass accretion rates larger than a critical rate of approximately 10(exp 18) g/s, X-ray emission from the column accretion is significantly affected by radiation drag. Although the mass accretion rate increases dramatically during outbursts, the observed hard (greater than 2 keV) X-ray luminosity will not rise proportionately. The slope and peak energy of the outburst spectra are only weakly affected. We conclude that the observed X-ray spectra can be explained by this simple analytic solution and that the production of hard X-rays from the accretion shock at the magnetic poles in the intermediate polars is in general agreement with the observations. However, since the X-ray emission and absorption depend on the mass accretion rate in a complicated manner, observed hard X-ray luminosities (greater than 2 keV) are not a good indicator of the mass accretion rate in the X-ray emitting regions.

Yi, Insu↗

Design of structure/control systems with transient response constraints exhibiting relative minima

Structural optimization problems involving dynamic behavior constraints often exhibit nonconvex design spaces. The direct application of a global optimization algorithm requires a large number of function evaluations which in term require a large number of dynamic structural analyses. This work presents a strategy aimed at finding the global optimum for problems with transient dynamic behavior constraints based on approximation concepts. The method consists of generating and solving a sequence of approximate problems using a global optimizer. The approximations are explicit and capture the inherent nonconvexity of the exact functions. A simple example problem is presented.

Sepulveda, A. E.↗

Simplified theory of large-amplitude wave propagation

An orbit perturbation procedure was applied to the description of monochromatic, large-amplitude, electrostatic plasma wave propagation. In the lowest order approximation, untrapped electrons were assumed to follow constant-velocity orbits and trapped electrons were assumed to execute simple harmonic motion. The deviations of these orbits from the actual orbits were regarded as perturbations. The nonlinear damping rate and frequency shift were then obtained in terms of simple functions. The results are in good agreement with previous less approximate analyses.

Kim, H.↗