Two-dimensional splines
Noteworthy properties of one-dimensional splines are presented. Stochastic processes, kernel functions, interpolation, and data smoothing were discussed. Theorems and concomitant proofs are provided.
Utreras, F. I.↗
SEARCH · Search NASA
Results for “Kernel function”
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.
Exploratory flutter test in a cryogenic wind tunnel
A model consisting of a rigid wing with an integral, flexible beam support that was cantilever mounted from the wall in the NASA LaRC 0.3-m transonic cryogenic tunnel was used in a flutter analysis study. The wing had a rectangular planform of aspect ratio 1.5 and a 64A010 airfoil. Various considerations and procedures for conducting flutter tests in a cryogenic wind tunnel were evaluated. Flutter onset conditions were established from extrapolated subcritical response measurements. A flutter boundary was determined at cryogenic temperatures over a Mach number M range from 0.5 to 0.9. Flutter was obtained at two different Reynolds numbers R at M = 0.5 (R = 4.4 and 18.4 x 10 to the 6th power) and at M = 0.8 (R = 5.0 and 10.4 x 10 to the 6th power). Flutter analyses using subsonic lifting surface (kernel function) aerodynamics were made over the range of test conditions. To evaluate the Reynolds number effects at M = 0.5 and 0.8, the experimental results were adjusted using analytical trends to account for differences in the model test temperatures and mass ratios. The adjusted experimental results indicate that increasing Reynolds number from 5.0 to 20.0 x 10 to the 6th power decreased the dynamic pressure by 4.0 to 6.5 percent at M = 0.5 and 0.8.
Cole, S. R.↗
Divergence study of a high-aspect ratio, forward-swept wing
An experimental wind-tunnel study to determine the divergence characteristics of a high-aspect ratio, forward-swept wing has been conducted in the NASA Langley Research Center (LaRC) Transonic Dynamics Tunnel (TDT). The rectangular wing used for this study had a panel aspect ratio of 9.16 (lambda = 0 deg.) and the sweep angle could be set at lambda = 0 deg., -15 deg., -30 deg., -45 deg., or -60 deg. A rectangular wing tip shape was tested at each of these sweep angles. In addition, a tip shape parallel to the freestream flow was tested for a wing sweep angle of lambda = -45 deg. The root of the wing was cantilever mounted to the wall of the wind tunnel. Divergence conditions were measured at M = 0.4 for each sweep angle and tip configuration tested. Subcritical response techniques were used to extrapolate to the divergence conditions during the wind-tunnel test. The primary objective of this test was to obtain data which could be used to verify for this configuration the divergence prediction capability of an aeroelastic analysis code. Subsonic lifting surface theory (kernel function) aerodynamics are utilized by this particular code. The analytical predictions of divergence were found to be significantly conservative at all forward sweep angles. At lambda = -45 deg., the analysis was 14 percent conservative. The effect of the two tip shapes on the divergence dynamic pressure was predicted accurately by the analysis. The divergence condition for the tip shape parallel to the flow occurred at a dynamic pressure 14 percent higher than the divergence condition with a rectangular tip shape.
Cole, S. R.↗
Divergence study of a high-aspect ratio, forward-swept wing
A study has been conducted in the NASA-Langley Transonic Dynamics wind tunnel to determine the divergence characteristics of a high aspect ratio forward swept wing of rectangular planform, for sweep angles of zero, -15, -30, -45, and -60 deg. In addition to a rectangular wingtip, a tip whose geometry lay parallel with the freestream flow was tested at a sweep angle of -45 deg. The primary objective was to obtain data that could be used to verify the divergence prediction capabilities of an aeroelastic analysis code employing kernel function aerodynamics. The analytical predictions of convergence were found to be conservative for all forward sweep angles, and the effect of the two tip shapes on divergence dynamic pressure were accurately predicted.
Cole, S. R.↗
Some numerical calculations by using linear classical sonic theories approached from sub- or supersonic speeds
The unsteady aerodynamics of a two-dimensional wing at sonic speed are studied by using so-called classical sonic theories (linear), approached from supersonic flow (M=1+0) or subsonic flow (M=1-0). In the former approach, the exact expressions of lift and lift distribution are obtained in terms of Fresnel integrals, while in the latter approach an integral equation must be solved, the kernel function of which is obtained from the subsonic Possio's equation and has a root singularity. The discrete analysis is adopted on the basis of the semicircle method (SCM) and the weighting function for subsonic-flow-Gauss-quadrature, as well as modified characteristics obtained from both approaches agree quite well with each other. The results obtained by the present computations are compared with those of DLM-C (subsonic 2D code) developed by ANDO et al, and are found to give a reasonable outer boundary for subsonic unsteady aerodynamics.
Yamamoto, Y.↗
Flutter analysis of highly swept delta wings by conventional methods
The flutter boundaries of six thin highly-swept delta-platform wings have been calculated. Comparisons are made between experimental data and results using several aerodynamic methods. The aerodynamic methods used include a subsonic and supersonic kernel function, second order piston theory, and a transonic small disturbance code. The dynamic equations of motion are solved using analytically calculated mode shapes and frequencies.
Gibbons, M. D.↗
Some low-speed flutter characteristics of simple low-aspect-ratio delta wing models
Some results from a combined experimental and analytical study of the low-speed flutter characteristics of low-aspect-ratio delta wings are presented. Data are presented which show the effects of sweep angle on the flutter characteristics of some simple plate models of constant planform area. The range of sweep angles studied was from 30 to 72 degrees. In addition, flutter results are presented for two 30 deg-sweep clipped-delta wing models. Further, results are presented that show the effects of root clamping (percentage length of the root chord that is cantilevered) for a 45 deg-sweep delta wing. The experimental data are compared with analytical results obtained by using kernel function and doublet lattice subsonic unsteady lifting surface theories.
Doggett, Robert V., Jr.↗
Flutter of a Low-Aspect-Ratio Rectangular Wing
A flutter test of a low-aspect-ratio rectangular wing was conducted in the Langley Transonic Dynamics Tunnel (TDT). The model used in this flutter test consisted of a rigid wing mounted to the wind-tunnel wall by a flexible, rectangular beam. The flexible support shaft was connected to the wing root and was cantilever mounted to the wind-tunnel wall. The wing had an aspect ratio of 1.5 based on the wing semispan and an NACA 64A010 airfoil shape. The flutter boundary of the model was determined for a Mach number range of 0.5 to 0.97. The shape of the transonic flutter boundary was determined. Actual flutter points were obtained on both the subsonic and supersonic sides of the flutter bucket. The model exhibited a deep transonic flutter bucket over a narrow range of Mach number. At some Mach numbers, the flutter conditions were extrapolated using a subcritical response technique. In addition to the basic configuration, modifications were made to the model structure such that the first bending frequency was changed without significantly affecting the first torsion frequency. The experiment showed that increasing the bending stiffness of the model support shaft through these modifications lowered the flutter dynamic pressure. Flutter analysis was conducted for the basic model as a comparison with the experimental results. This flutter analysis was conducted with subsonic lifting-surface (kernel function) aerodynamics using the k method for the flutter solution.
Cole, Stanley R.↗
Some low-speed flutter characteristics of simple low-aspect-ratio delta wing models
Some results from a combined experimental and analytical study of the low-speed flutter characteristics of low-aspect-ratio delta wings are presented. Data are presented which show the effects of sweep angle on the flutter characteristics of some simple plate models of constant planform area. The range of sweep angles studied was from 30 to 72 degrees. In addition, flutter results are presented for two 30 deg-sweep clipped-delta wing models. Further results are presented that show the effects of root clamping (percentage length of the root chord that is cantilevered) for a 45 deg-sweep delta wing. The experimental data are compared with analytical results obtained by using kernel function and doublet lattice subsonic unsteady lifting surface theories.
Doggett, Robert V., Jr.↗
Application of the CAP-TSD unsteady transonic small disturbance program to wing flutter
The application and assessment of a computer program called CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) for flutter predictions are described. Flutter calculations are presented for two thin swept-and-tapered wing planforms with well-defined modal properties. One planform is a series of 45-degree swept wings and the other planform is a clipped delta wing. Comparisons are made between the results of CAP-TSD using the linear equation and no airfoil thickness and the results obtained from a subsonic kernel function analysis. The calculations cover a Mach number range from low subsonic to low supersonic values, including the transonic range, and are compared with subsonic linear theory and experimental data. It is noted that since both wings have very thin airfoil sections, the effects of thickness are minimal.
Bennett, Robert M.↗
Boundary-element shape sensitivity analysis for thermal problems with nonlinear boundary conditions
Implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solid objects subjected to nonlinear boundary conditions is shown to generate an accurate and economical approach for the computation of shape sensitivities for this class of problems. This approach involves the employment of analytical derivatives of boundary-element kernel functions with respect to shape design variables. A formulation is presented that can consistently account for both temperature-dependent convection and radiation boundary conditions. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems are presented that have closed-form solutions.
Kane, James H.↗
On the nonlinear three dimensional instability of Stokes layers and other shear layers to pairs of oblique waves
The nonlinear evolution of a pair of initially oblique waves in a high Reynolds Number Stokes layer is studied. Attention is focused on times when disturbances of amplitude epsilon have O(epsilon(exp 1/3)R) growth rates, where R is the Reynolds number. The development of a pair of oblique waves is then controlled by nonlinear critical-layer effects. Viscous effects are included by studying the distinguished scaling epsilon = O(R(exp -1)). This leads to a complicated modification of the kernel function in the integro-differential amplitude equation. When viscosity is not too large, solutions to the amplitude equation develop a finite-time singularity, indicating that an explosive growth can be introduced by nonlinear effects; we suggest that such explosive growth can lead to the bursts observed in experiments. Increasing the importance of viscosity generally delays the occurrence of the finite-time singularity, and sufficiently large viscosity may lead to the disturbance decaying exponentially. For the special case when the streamwise and spanwise wavenumbers are equal, the solution can evolve into a periodic oscillation. A link between the unsteady critical-layer approach to high-Reynolds-number flow instability, and the wave vortex approach is identified.
Wu, Xuesong↗
FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons
Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.
Joseph, Rose↗
Sensitivity Analysis of Wing Aeroelastic Responses
Design for prevention of aeroelastic instability (that is, the critical speeds leading to aeroelastic instability lie outside the operating range) is an integral part of the wing design process. Availability of the sensitivity derivatives of the various critical speeds with respect to shape parameters of the wing could be very useful to a designer in the initial design phase, when several design changes are made and the shape of the final configuration is not yet frozen. These derivatives are also indispensable for a gradient-based optimization with aeroelastic constraints. In this study, flutter characteristic of a typical section in subsonic compressible flow is examined using a state-space unsteady aerodynamic representation. The sensitivity of the flutter speed of the typical section with respect to its mass and stiffness parameters, namely, mass ratio, static unbalance, radius of gyration, bending frequency, and torsional frequency is calculated analytically. A strip theory formulation is newly developed to represent the unsteady aerodynamic forces on a wing. This is coupled with an equivalent plate structural model and solved as an eigenvalue problem to determine the critical speed of the wing. Flutter analysis of the wing is also carried out using a lifting-surface subsonic kernel function aerodynamic theory (FAST) and an equivalent plate structural model. Finite element modeling of the wing is done using NASTRAN so that wing structures made of spars and ribs and top and bottom wing skins could be analyzed. The free vibration modes of the wing obtained from NASTRAN are input into FAST to compute the flutter speed. An equivalent plate model which incorporates first-order shear deformation theory is then examined so it can be used to model thick wings, where shear deformations are important. The sensitivity of natural frequencies to changes in shape parameters is obtained using ADIFOR. A simple optimization effort is made towards obtaining a minimum weight design of the wing, subject to flutter constraints, lift requirement constraints for level flight and side constraints on the planform parameters of the wing using the IMSL subroutine NCONG, which uses successive quadratic programming.
Issac, Jason Cherian↗
Free-Spinning-Tunnel Investigation of a 1/20-Scale Model of an Unswept-Wing Jet-Propelled Trainer Airplane
A flutter analysis employing the kernel function for three- dimensional, subsonic, compressible flow is applied to a flutter-tested tail surface which has an aspect ratio of 3.5, a taper ratio of 0.15, and a leading-edge sweep of 30 deg. Theoretical and experimental results are compared at Mach numbers from 0.75 to 0.98. Good agreement between theoretical and experimental flutter dynamic pressures and frequencies is achieved at Mach numbers to 0.92. At Mach numbers from 0.92 to 0.98, however, a second solution to the flutter determinant results in a spurious theoretical flutter boundary which is at a much lower dynamic pressure and at a much higher frequency than the experimental boundary.
Bowman, James S., Jr.↗
Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations That Arise from Nonlinear Two-Point Boundary Value Problems
In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + definite integral of g(x, t)F(t,y(t))dt with limits between 0 and 1,0 less than or equal to x les than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integral equations arise, e.g., when one applied Green's function techniques to nonlinear two-point boundary value problems of the form y "(x) =f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and y(l) = y(sub l), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trepezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal rule, thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.
Sidi, Avram↗
Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations that Arise from Nonlinear Two-Point Boundary Value Problems
In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + integral(0 to 1) g(x,t) F(t, y(t)) dt, 0 less than or equal to x less than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integrals equations arise, e.g., when one applies Green's function techniques to nonlinear two-point boundary value problems of the form U''(x) = f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and g(l) = y(sub 1), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trapezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.
Sidi, Avram↗
Density Estimation with Mercer Kernels
We present a new method for density estimation based on Mercer kernels. The density estimate can be understood as the density induced on a data manifold by a mixture of Gaussians fit in a feature space. As is usual, the feature space and data manifold are defined with any suitable positive-definite kernel function. We modify the standard EM algorithm for mixtures of Gaussians to infer the parameters of the density. One benefit of the approach is it's conceptual simplicity, and uniform applicability over many different types of data. Preliminary results are presented for a number of simple problems.
Macready, William G.↗