Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth Basis 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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Determination of discontinuities in the derivative of a piecewise-smooth function on the basis of the results of measurements

A piecewise-smooth function with discontinuity in the first derivative on a given interval is considered. The values of the function are measured at a sequence of points in the interval and a random error is included in the measurements. A method is proposed to estimate the position of the discontinuity in the derivative. Regression lines are associated with each measurement point and account for k - 1 points preceding or following the point. The estimate for the position of the discontinuity is the measurement point with the largest angle between the regression lines. The error in the estimate is analyzed and the results are verified.

Gurin, L. S.↗
NASA NTRS1975

Solutions to Kuessner's integral equation in unsteady flow using local basis functions

The computational procedure and numerical results are presented for a new method to solve Kuessner's integral equation in the case of subsonic compressible flow about harmonically oscillating planar surfaces with controls. Kuessner's equation is a linear transformation from pressure to normalwash. The unknown pressure is expanded in terms of prescribed basis functions and the unknown basis function coefficients are determined in the usual manner by satisfying the given normalwash distribution either collocationally or in the complex least squares sense. The present method of solution differs from previous ones in that the basis functions are defined in a continuous fashion over a relatively small portion of the aerodynamic surface and are zero elsewhere. This method, termed the local basis function method, combines the smoothness and accuracy of distribution methods with the simplicity and versatility of panel methods. Predictions by the local basis function method for unsteady flow are shown to be in excellent agreement with other methods. Also, potential improvements to the present method and extensions to more general classes of solutions are discussed.

Fromme, J. A.↗
NASA NTRS1995

Exponentially accurate approximations to piece-wise smooth periodic functions

A family of simple, periodic basis functions with 'built-in' discontinuities are introduced, and their properties are analyzed and discussed. Some of their potential usefulness is illustrated in conjunction with the Fourier series representations of functions with discontinuities. In particular, it is demonstrated how they can be used to construct a sequence of approximations which converges exponentially in the maximum norm to a piece-wise smooth function. The theory is illustrated with several examples and the results are discussed in the context of other sequences of functions which can be used to approximate discontinuous functions.

Greer, James↗
NASA NTRS1998

Basis Selection for Wavelet Regression

A wavelet basis selection procedure is presented for wavelet regression. Both the basis and the threshold are selected using cross-validation. The method includes the capability of incorporating prior knowledge on the smoothness (or shape of the basis functions) into the basis selection procedure. The results of the method are demonstrated on sampled functions widely used in the wavelet regression literature. The results of the method are contrasted with other published methods.

Wheeler, Kevin R.↗
NASA NTRS1995

Multiresolution Representation Using Biorthogonal Multiwavelets

We generalize Harten's multiresolution representation to biorthogonal multiwavelets. Several variants are considered. For example, a given array of discrete point values is transformed to point values and derivatives or point 'values and cell averages'. Compact Hermite interpolation is used in the decomposition and reconstruction algorithm. The resulting basis functions that are symmetric or skewsymmetric, compact, and smooth with optimal order accuracy. Harten's approach has several advantages: the multiresolution scheme is inherently discrete, non-periodic boundary conditions are easy to implement, and the representation can be extended to unstructured grids in bounded domains. We demonstrate the compression features of the new mutliwavelets by application to variable scale piecewise smooth functions with jump discontinuities typical of numerical solutions of nonlinear hyperbolic conservation laws.

Warming, Robert F.↗
NASA NTRS1997

Complex Volume Grid Generation Through the Use of Grid Reusability

This paper presents a set of surface and volume grid generation techniques which reuse existing surface and volume grids. These methods use combinations of data manipulations to reduce grid generation time, improve grid characteristics, and increase the capabilities of existing domain discretization software. The manipulation techniques utilize physical and computational domains to produce basis function on which to operate and modify grid character and smooth grids using Trans-Finite Interpolation, a vector interpolation method and parametric re-mapping technique. With these new techniques, inviscid grids can be converted to viscous grids, multiple zone grid adaption can be performed to improve CFD solver efficiency, and topological changes to improve modeling of flow fields can be done simply and quickly. Examples of these capabilities are illustrated as applied to various configurations.

Alter, Stephen J.↗
NASA NTRSDate not supplied

Spacesuit and Mobility Performance Changes

The complex interactions between the human body and spacesuit lead to changes inmovement patterns and mobility performancesof the wearer. In general, factors including the geometric properties, such as shape and size, mechanical properties of the suit, and pressurization of the suit are known to be associated with altered movement patterns as compared to an unsuited human.However, their relative contributions have not been explicitly quantified from the mobility performance perspectives. The goal of this study wasthus to assessthe effects fromthe different types of mobility constraintconditions, namely by wearing either a 3D printed hard upper torso (HUT) assembly orfully pressurized spacesuit. The outcome was also compared against the unsuited motions. For this study, an xEMU (exploration Extravehicular Mobility Unit) suitwas considered, which is the next generation spacesuit developed by NASA. In each test condition, the subject was asked to move the arm and hand as prescribedfor different task types, and the corresponding body segment locations were recordedusing a 3D motion capture system. The following three tasks were performedand analyzed: 1) Outward one-handed reaches:the subject in a standing pose made sweeping motions with the extended right arm from the extreme end-to-end positions, including side-to-side at different elevations and top-to-bottom at different azimuths. 2) Outward two-handed reaches:similarto the previous task, howeverthe subject kept the hands together during the motions in order to assessthe areas that can be reached by both hands. 3) Inward one-handed reaches:the subject made right-hand reach motions to the surface of the HUT.The hand traces collected from each task were modeled by a template shape parametrically deformed with a radial basis function. This process enabled foran abstraction of the hand traces into a smooth surface envelope representing the maximally reachable area of the test subject, of which the shapes and sizes were compared across the different test conditions. The preliminary analysis has shown that theoverall size of reachenvelopes decreases in a pressurized suit compared to 3D printed mockup HUT and unsuited conditions. The specific shape of the envelopes, which were determined by the reachable and unreachable zones, alsovary with the testconditions, and the differences werepronounced with the inward reaches to the HUT surface. The latter observation ispotentially relatedto the increased demandfor shoulder and elbow flexions.Overcomingthe resistance from the pressurizedsoft goods and mechanical constraints of the shoulder assemblywas seen to be associated with the difference in motion patterns between the suited and unsuited conditions. Overall, the information quantified from this study is expected to provide structured metrics for spacesuit mobility, which can improve design optimization and human-system integration.

K Han Kim↗
NASA NTRS2003

Molecular and functional analyses of the contractile apparatus in lymphatic muscle

Lymphatics are necessary for the generation and regulation of lymph flow. Lymphatics use phasic contractions and extrinsic compressions to generate flow; tonic contractions alter resistance. Lymphatic muscle exhibits important differences from typical vascular smooth muscle. In this study, the thoracic duct exhibited significant functional differences from mesenteric lymphatics. To understand the molecular basis for these differences, we examined the profiles of contractile proteins and their messages in mesenteric lymphatics, thoracic duct, and arterioles. Results demonstrated that mesenteric lymphatics express only SMB smooth muscle myosin heavy chain (SM-MHC), whereas thoracic duct and arterioles expressed both SMA and SMB isoforms. Both SM1 and SM2 isoforms of SM-MHC were detected in arterioles and mesenteric and thoracic lymphatics. In addition, the fetal cardiac/skeletal slow-twitch muscle-specific beta-MHC message was detected only in mesenteric lymphatics. All four actin messages, cardiac alpha-actin, vascular alpha-actin, enteric gamma-actin, and skeletal alpha-actin, were present in both mesenteric lymphatics and arterioles. However, in thoracic duct, predominantly cardiac alpha-actin and vascular alpha-actin were found. Western blot and immunohistochemical analyses corroborated the mRNA studies. However, in arterioles only vascular alpha-actin protein was detected. These data indicate that lymphatics display genotypic and phenotypic characteristics of vascular, cardiac, and visceral myocytes, which are needed to fulfill the unique roles of the lymphatic system.

Non-NASA Center↗
NASA NTRS2001

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗
NASA NTRS1995

A 3D High-Order Unstructured Finite-Volume Algorithm for Solving Maxwell's Equations

A three-dimensional finite-volume algorithm based on arbitrary basis functions for time-dependent problems on general unstructured grids is developed. The method is applied to the time-domain Maxwell equations. Discrete unknowns are volume integrals or cell averages of the electric and magnetic field variables. Spatial terms are converted to surface integrals using the Gauss curl theorem. Polynomial basis functions are introduced in constructing local representations of the fields and evaluating the volume and surface integrals. Electric and magnetic fields are approximated by linear combinations of these basis functions. Unlike other unstructured formulations used in Computational Fluid Dynamics, the new formulation actually does not reconstruct the field variables at each time step. Instead, the spatial terms are calculated in terms of unknowns by precomputing weights at the beginning of the computation as functions of cell geometry and basis functions to retain efficiency. Since no assumption is made for cell geometry, this new formulation is suitable for arbitrarily defined grids, either smooth or unsmooth. However, to facilitate the volume and surface integrations, arbitrary polyhedral cells with polygonal faces are used in constructing grids. Both centered and upwind schemes are formulated. It is shown that conventional schemes (second order in Cartesian grids) are equivalent to the new schemes using first degree polynomials as the basis functions and the midpoint quadrature for the integrations. In the new formulation, higher orders of accuracy are achieved by using higher degree polynomial basis functions. Furthermore, all the surface and volume integrations are carried out exactly. Several model electromagnetic scattering problems are calculated and compared with analytical solutions. Examples are given for cases based on 0th to 3rd degree polynomial basis functions. In all calculations, a centered scheme is applied in the interior, while an upwind matching scheme is employed at material interfaces and the Engquist-Majda non-reflecting boundary condition is implemented at the numerical outer boundaries. The staggered leapfrog scheme and the Runge-Kutta methods are utilized for the time integration. Excellent agreements are found between the numerical and analytical solutions.

Liu, Yen↗
NASA NTRS1977

In-flight gyro drift rate calibration on the Viking Orbiters

The drift rates of the attitude control gyros onboard the Viking Orbiters were calibrated several times during flight. The calibration was performed by engaging the gyro control of attitude for a six-hour period, recording the Viking Orbiter coordinates of the sun and the roll reference star as a function of time, and processing these data by a computer program. The computer program calculated drift rates for each data increment, and then smoothed the calculated rate versus the time function for all the data on the basis of a Gauss-Markov gyro-drift model. The effects of data noise were minimized by using a relatively long time constant in the Gauss-Markov model and by empirically determining a set of corrections for celestial inertial scale factor error differences. In this way the rms error of the data fit was lessened.

Breckenridge, W. G.↗
NASA NTRS2011

Development of the Orion Crew Module Static Aerodynamic Database: Hypersonic - Part 1

The Orion aerodynamic database provides force and moment coefficients given the velocity, attitude, configuration, etc. of the Crew Exploration Vehicle (CEV). The database is developed and maintained by the NASA CEV Aerosciences Project team from computational and experimental aerodynamic simulations. The database is used primarily by the Guidance, Navigation, and Control (GNC) team to design vehicle trajectories and assess flight performance. The initial hypersonic re-entry portion of the Crew Module (CM) database was developed in 2006. Updates incorporating additional data and improvements to the database formulation and uncertainty methodologies have been made since then. This paper details the process used to develop the CM database, including nominal values and uncertainties, for Mach numbers greater than 8 and angles of attack between 140deg and 180deg. The primary available data are more than 1000 viscous, reacting gas chemistry computational simulations using both the Laura and Dplr codes, over a range of Mach numbers from 2 to 37 and a range of angles of attack from 147deg to 172deg. Uncertainties were based on grid convergence, laminar-turbulent solution variations, combined altitude and code-to-code variations, and expected heatshield asymmetry. A radial basis function response surface tool, NEAR-RS, was used to fit the coefficient data smoothly in a velocity-angle-of-attack space. The resulting database is presented and includes some data comparisons and a discussion of the predicted variation of trim angle of attack and lift-to-drag ratio. The database provides a variation in trim angle of attack on the order of +/-2deg, and a range in lift-to-drag ratio of +/-0.035 for typical vehicle flight conditions.

Bibb, Karen L.↗
NASA NTRS1993

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗
NASA NTRS1972

Solution of the Boltzmann kinetic equation for the relaxation of a gas mixture

The temporal behavior is considered of a quiescent mixture of gases of different temperatures with spatially uniform distribution. The process of heating a cold gas by a hot gas is treated on the basis of the Boltzmann kinetic equation. The mixture is assumed to be composed of absolutely hard smooth spheres, and the initial distribution functions for each gas is taken to the Maxwellian. With such a choice of initial distribution functions, it is shown that the solution of the Boltzmann kinetic equation depends only on the velocity modulus and the time.

Rykov, V. A.↗
NASA NTRS1975

The use of local basis functions in unsteady aerodynamics

This paper reports a successful demonstration of using local basis functions to solve Kussner's integral equation for unsteady airloads. Pressure basis functions are defined over relatively small subsets of the aerodynamic surface and are zero elsewhere. This formulation combines the accuracy and smoothness of distribution methods with the simplicity and versatility of panel methods. Application is made to subsonic compressible flow about oscillating planar surfaces with controls. Numerical results are presented. Convergence with respect to panel size, extensions to more general classes of solutions, and areas needing additional investigation are discussed.

Fromme, J.↗
NASA NTRS2013

RUSHMAPS: Real-Time Uploadable Spherical Harmonic Moment Analysis for Particle Spectrometers

RUSHMAPS is a new onboard data reduction scheme that gives real-time access to key science parameters (e.g. moments) of a class of heliophysics science and/or solar system exploration investigation that includes plasma particle spectrometers (PPS), but requires moments reporting (density, bulk-velocity, temperature, pressure, etc.) of higher-level quality, and tolerates a lowpass (variable quality) spectral representation of the corresponding particle velocity distributions, such that telemetry use is minimized. The proposed methodology trades access to the full-resolution velocity distribution data, saving on telemetry, for real-time access to both the moments and an adjustable-quality (increasing quality increases volume) spectral representation of distribution functions. Traditional onboard data storage and downlink bandwidth constraints severely limit PPS system functionality and drive cost, which, as a consequence, drives a limited data collection and lower angular energy and time resolution. This prototypical system exploit, using high-performance processing technology at GSFC (Goddard Space Flight Center), uses a SpaceCube and/or Maestro-type platform for processing. These processing platforms are currently being used on the International Space Station as a technology demonstration, and work is currently ongoing in a new onboard computation system for the Earth Science missions, but they have never been implemented in heliospheric science or solar system exploration missions. Preliminary analysis confirms that the targeted processor platforms possess the processing resources required for realtime application of these algorithms to the spectrometer data. SpaceCube platforms demonstrate that the target architecture possesses the sort of compact, low-mass/power, radiation-tolerant characteristics needed for flight. These high-performing hybrid systems embed unprecedented amounts of onboard processing power in the CPU (central processing unit), FPGAs (field programmable gate arrays), and DSP (digital signal processing) elements. The fundamental computational algorithm de constructs 3D velocity distributions in terms of spherical harmonic spectral coefficients (which are analogous to a Fourier sine-cosine decomposition), but uses instead spherical harmonics Legendre polynomial orthogonal functions as a basis for the expansion, portraying each 2D angular distribution at every energy or, geometrically, spherical speed-shell swept by the particle spectrometer. Optionally, these spherical harmonic spectral coefficients may be telemetered to the ground. These will provide a smoothed description of the velocity distribution function whose quality will depend on the number of coefficients determined. Successfully implemented on the GSFC-developed processor, the capability to integrate the proposed methodology with both heritage and anticipated future plasma particle spectrometer designs is demonstrated (with sufficiently detailed design analysis to advance TRL) to show specific science relevancy with future HSD (Heliophysics Science Division) solar-interplanetary, planetary missions, sounding rockets and/or CubeSat missions.

Figueroa-Vinas, Adolfo↗
NASA NTRS1997

Exponential Approximations Using Fourier Series Partial Sums

The problem of accurately reconstructing a piece-wise smooth, 2(pi)-periodic function f and its first few derivatives, given only a truncated Fourier series representation of f, is studied and solved. The reconstruction process is divided into two steps. In the first step, the first 2N + 1 Fourier coefficients of f are used to approximate the locations and magnitudes of the discontinuities in f and its first M derivatives. This is accomplished by first finding initial estimates of these quantities based on certain properties of Gibbs phenomenon, and then refining these estimates by fitting the asymptotic form of the Fourier coefficients to the given coefficients using a least-squares approach. It is conjectured that the locations of the singularities are approximated to within O(N(sup -M-2), and the associated jump of the k(sup th) derivative of f is approximated to within O(N(sup -M-l+k), as N approaches infinity, and the method is robust. These estimates are then used with a class of singular basis functions, which have certain 'built-in' singularities, to construct a new sequence of approximations to f. Each of these new approximations is the sum of a piecewise smooth function and a new Fourier series partial sum. When N is proportional to M, it is shown that these new approximations, and their derivatives, converge exponentially in the maximum norm to f, and its corresponding derivatives, except in the union of a finite number of small open intervals containing the points of singularity of f. The total measure of these intervals decreases exponentially to zero as M approaches infinity. The technique is illustrated with several examples.

Banerjee, Nana S.↗
NASA NTRS2005

Determining the Sun's Deep Meridional Flow Speed Using Active Latitude Drift Rates Since 1874

Dynamo models that incorporate a deep meridional return flow indicate that this flow regulates both the period and the amplitude of the sunspot cycle. We recently examined the equatorward drift of the active latitudes (as given by the centroid of the sunspot areas in each hemisphere) and found evidence supporting this view. In those studies we fit the equatorward drift in each hemisphere for each sunspot cycle with a simple parabola - giving us a drift rate and its deceleration for each hemisphere/cycle. Here we analyze the same data (the Royal Greenwich Observatory/JSAF/NOAA daily active region summaries) to determine the drift rates in each hemisphere on a yearly basis (rotation-by- rotation measurements smoothed to remove high frequencies) and fit them with a simple model for the meridional flow that provides the meridional flow speed as a function of latitude and time from 1874 to 2005. These flow speeds can be used to test dynamo models - some of which have predictive capabilities.

Hathaway, David H.↗