SEARCH · Engineering Papers
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.
Average energies of ground and excited configurations in highly ionized atoms, using orthogonalized screened hydrogenic radial functions
It is shown that the contribution of higher harmonics to the lunar transfer functions for the tangential components of the surface magnetic field is significant at frequencies greater than 0.01 Hz. The inclusion of the higher harmonics shows that there are two distinct transfer functions corresponding to the components of the tangential surface magnetic field perpendicular and parallel to the direction of the wave vector of the external disturbance forcing the lunar induction. The dependences of these transfer functions on frequency and location are determined. The effects of the higher harmonics can: (1) account for a hitherto unexplained feature in the Apollo 12-Explorer 35 transfer functions, namely the rolloff at high frequencies; and (2) offer a possible explanation for the frequency dependence of the difference between the transfer functions for the two orthogonal components of the surface magnetic field. The harmonic response of a simple current layer model of the moon is derived.
It is shown that the contribution of higher harmonics to the lunar transfer functions for the tangential components of the surface magnetic field is significant at frequencies greater than 0.01 Hz. The inclusion of the higher harmonics shows that there are two distinct transfer functions corresponding to the components of the tangential surface magnetic field perpendicular and parallel to the direction of the wave vector of the external disturbance forcing the lunar induction. The dependences of these transfer functions on frequency and location are determined. The effects of the higher harmonics can (1) account for a hitherto unexplained feature in the Apollo 12 and Explorer 35 transfer functions, namely the rolloff at high frequencies, and (2) offer a possible explanation for the frequency dependence of the difference between the transfer functions for the two orthogonal components of the surface magnetic field.
The question of determining when a Bernstein-Bezier cubic curve in the plane can be represented as the graph of function in some fixed orthogonal coordinate system is considered. The notion of a curve being monotone in a given direction is introduced to aid in the analysis.
An interactive, graphics-based grid generation system for 3D turbomachinery geometries has been developed. The system consists of separate modules for geometry modeling and grid generation. The grid generation method generates a series of 2D grids in the blade-to-blade passage to build up the 3D grid. A Poisson equation with forcing functions selected to control orthogonality and spacing on all boundaries is solved to generate the grid. Interactive definition of Bezier curves and surfaces as internal boundaries is used to improve control of grid quality. A multi-block data structure simpifies the creation of structured H-grids about complex turbomachinery geometries. The multi-block strategy facilitates the creation of a grid in the tip flow region, which is an important contributor to losses within the passage. The grid generation system is applied to several complex geometries including a simple radial impeller, a conventional turbine rotor, the SSME LOX turbine rotor with and without a tip flow grid, a multi-passage inducer-impeller, and a tip-flow cascade. In addition, a calculation of the tip-flow cascade flow field is shown.
One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as linearity, of being stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposition data, the HHT allows spectrum analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a near orthogonal adaptive basis, a basis that is derived from the data. The IMFs can be further analyzed for spectrum interpretation by the classical Hilbert Transform. A new engineering spectrum analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications post additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs near orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the developments of new HHT processing options, such as real-time and 2-D processing using Field Programmable Array (FPGA) computational resources, enhanced HHT synthesis, and broaden the scope of HHT applications for signal processing.
One of the main traditional tools used in scientific and engineering data spectral analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as being linear and stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectral analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposed data, the HHT allows spectral analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real-value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a nearly orthogonal derived from the data (adaptive) basis. The IMFs can be further analyzed for spectrum content by using the classical Hilbert Transform. A new engineering spectral analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications pose additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs nearly orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the development of new HHT processing options, such as real-time and 2-D processing using Field Programmable Gate Array (FPGA) computational resources,
The B-spline R-matrix method is employed to investigate the low-energy elastic electron scattering by atomic oxygen. Flexible non-orthogonal sets of radial functions are used to construct the target description and to represent the scattering functions. A detailed investigation regarding the dependence of the predicted partial and total cross sections on the scattering model and the accuracy of the target description is presented. The predicted angle-integrated elastic cross sections are in good agreement with experiment, whereas significant discrepancies are found in the angle-differential elastic cross sections near the forward direction. .The near-threshold results are found to strongly depend on the treatment of inner-core short-range correlation effects in the target description, as well as on a proper account of the target polarizability. A sharp increase in the elastic cross sections below 1 eV found in some earlier calculations is judged to be an artifact of an unbalanced description of correlation in the N-electron target structure and the (N+l)-electron-collision problems.
A methodology for surface grid generation in three dimensional space is described. The method solves a Poisson equation for each coordinate on arbitrary surfaces using successive line over-relaxation. The complete surface curvature terms were discretized and retained within the nonhomogeneous term in order to preserve surface definition; there is no need for conventional surface splines. Control functions were formulated to permit control of grid orthogonality and spacing. A method for interpolation of control functions into the domain was devised which permits their specification not only at the surface boundaries but within the interior as well. An interactive surface generation code which makes use of this methodology is currently under development.
In-flight estimation of large structure model errors is usually required as a means of detecting inevitable deficiencies in large structure controller/estimator models. The present paper deals with a least-squares formulation which seeks to minimize a quadratic functional of the model errors. The properties of these error estimates are analyzed. It is shown that an arbitrary model error can be decomposed as the sum of two components that are orthogonal in a suitably defined function space. Relations between true and estimated errors are defined. The estimates are found to be approximations that retain many of the significant dynamics of the true model errors. Current efforts are directed toward application of the analytical results to a reference large structure model.
The probability, differential, and integral scattering cross sections of the 1{sup 1}S→2{sup 1}S and 1{sup 1}S→2{sup 1}P transitions of helium have been calculated in the first Born approximation. The projectile-target interactions depending on the temperature and the density of plasma are described by the Debye-Hückel model. Wave functions of the target before and after collision were modeled by non orthogonal Hartree-Fock orbitals. The wave functions parameters are calculated with the Ritz variational method. We improve our unscreened first Born approximation integral cross sections by using the BE-scaled (B stands for binding energy and E excitation energy) method. The second Born approximation has also been used to calculate the excitation cross sections in Debye plasma. Our calculations are compared to other theoretical and experimental results where applicable.
Fluid functions are expanded in a series of functions constituting a complete and orthogonal system of wave solutions of the linearized Navier-Stokes system. The motion of characteristic waves is described by Hamiltonian equations of ray acoustics. Amplitude of the waves is described with the aid of nonlinear interaction terms. The distribution functions (squares of amplitudes) satisfy Boltzmann-type equations and completely describe the mean turbulence properties.
Orthogonality relations are obtained for the spherical harmonic coefficients of functions defined on the surface of a sphere. Following a brief discussion of the orthogonality of Fourier series coefficients, consideration is given to the values averaged over all orientations of the coordinate system of the spherical harmonic coefficients of a function defined on the surface of a sphere that can be expressed in terms of Legendre polynomials for the special case where the function is the sum of two delta functions located at two different points on the sphere, and for the case of an essentially arbitrary function. It is noted that the orthogonality relations derived have found applications in statistical studies of the geomagnetic field.
Post-translational modifications (PTMs) can occur on almost all amino acids in eukaryotes as a key mechanism for regulating protein function. The ability to study the role of these modifications in various biological processes requires techniques to modify proteins site-specifically. One strategy for this is genetic code expansion (GCE) in bacteria. The low frequency of post-translational modifications in bacteria makes it a preferred host to study whether the presence of a post-translational modification influences a protein’s function. Genetic code expansion employs orthogonal translation systems engineered to incorporate a modified amino acid at a designated protein position. Selenoproteins, proteins containing selenocysteine, are also known to be post-translationally modified. Selenoproteins have essential roles in oxidative stress, immune response, cell maintenance, and skeletal muscle regeneration. Their complicated biosynthesis mechanism has been a hurdle in our understanding of selenoprotein functions. As technologies for selenocysteine insertion have recently improved, we wanted to create a genetic system that would allow the study of post-translational modifications in selenoproteins. By combining genetic code expansion techniques and selenocysteine insertion technologies, we were able to recode stop codons for insertion of N ε -acetyl- l -lysine and selenocysteine, respectively, into multiple proteins. The specificity of these amino acids for their assigned position and the simplicity of reverting the modified amino acid via mutagenesis of the codon sequence demonstrates the capacity of this method to study selenoproteins and the role of their post-translational modifications. Moreover, the evidence that Sec insertion technology can be combined with genetic code expansion tools further expands the chemical biology applications.
The main bottleneck in modeling transport in molecular devices is to develop the correct formulation of the problem and efficient algorithms for analyzing the electronic structure and dynamics using, for example, the time-dependent density functional theory. We have divided this task into several steps. The first step is to developing the right mathematical formulation and numerical algorithms for analyzing the electronic structure using density functional theory. The second step is to study time-dependent density functional theory, particularly the far-field boundary conditions. The third step is to study electronic transport in molecular devices. We are now at the end of the first step. Under DOE support, we have made subtantial progress in developing linear scaling and sub-linear scaling algorithms for electronic structure analysis. Although there has been a huge amount of effort in the past on developing linear scaling algorithms, most of the algorithms developed suffer from the lack of robustness and controllable accuracy. We have made the following progress: (1) We have analyzed thoroughly the localization properties of the wave-functions. We have developed a clear understanding of the physical as well as mathematical origin of the decay properties. One important conclusion is that even for metals, one can choose wavefunctions that decay faster than any algebraic power. (2) We have developed algorithms that make use of these localization properties. Our algorithms are based on non-orthogonal formulations of the density functional theory. Our key contribution is to add a localization step into the algorithm. The addition of this localization step makes the algorithm quite robust and much more accurate. Moreover, we can control the accuracy of these algorithms by changing the numerical parameters. (3) We have considerably improved the Fermi operator expansion (FOE) approach. Through pole expansion, we have developed the optimal scaling FOE algorithm.
An optimal Routh approximant of a single-input single-output dynamic system is a reduced-order transfer function of which the denominator is obtained by the Routh approximation method while the numerator is determined by minimizing a time-response integral-squared-error (ISE) criterion. In this paper, a new elegant approach is presented for obtaining the optimal Routh approximants for linear time-invariant continuous-time systems. The approach is based on the Routh canonical expansion, which is a finite-term orthogonal series of rational basis functions, and minimization of the ISE criterion. A procedure for combining the above approach with the bilinear transformation is also presented in order to obtain the optimal bilinear Routh approximants of linear time-invariant discrete-time systems. The proposed technique is simple in formulation and is amenable to practical implementation.
The non-orthogonal configuration interaction with fragments (NOCI-F) approach is extended opening the possibility to study intramolecular processes and materials with covalent or ionic lattices. So far, NOCI-F has been applied to study intermolecular energy and electron transfer employing ensembles of fragments that do not have atoms or bonds in common. The here presented approach divides the target system into two overlapping fragments that share one or more atoms and/or one or more bonds. After the construction of a collection of (multiconfigurational) fragment wave functions in a state specific optimization procedure, the fragment wave functions are combined to form many-electron basis functions for the non-orthogonal configuration interaction of the whole system. The orbitals in the overlapping fragment are defined by a corresponding orbital transformation of the fragment orbitals through a singular value decomposition. The overlapping fragments approach is first illustrated for a model system and then used to highlight some possible applications of NOCI with overlapping fragments. In conclusion, the results of excited state diffusion in transition metal oxide, intramolecular singlet fission and magnetic interactions in organic biradicals and ionic compounds are discussed.
An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.