Search NASA⌕ Search

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.

At least 235 records · Page 13

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Han Byul Woo↗

Assessing Potential Geophysical and Environmental Impacts from Frequent Rocket Launch Missions at Kennedy Space Center

Kennedy Space Center (KSC) in Florida has been utilized for space missions over many years with a gradually increasing number of rocket launches. There have been multiple studies where air-coupled acoustic waves and infrasound originating from launched rockets were used for operational purposes, such as locating booster trajectories as a function of changing atmospheric conditions. However, no study has utilized the acoustoelastic waves as a signal source for subsurface seismic investigations. We conducted a dispersion analysis using the seismic energy recorded from the Artemis I rocket launch at KSC in November 2022, and from these results we generated depth-sensitivity kernels at different wave frequencies. The kernels were compared with sedimentary core data to verify boundaries of carbonate layers above the Floridan Aquifer System. Accumulative information of bedrock-sediment boundary across the sedimentary platform could especially provide geo-structural evidence that manifests the configuration of coastal features. Continuous dispersion analysis of seismic recordings from consecutive rocket launches also has potential to identify non-stationary environmental effects, such as reorientation of sedimentary structures and fluctuation of the saltwater/groundwater lens from gravitational tides, which may affect erosional susceptibility of coastal features.

Rocket Launch Seismicity↗

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary systems. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗

Approximating the linear quadratic optimal control law for hereditary systems with delays in the control

The fundamental control synthesis issue of establishing a priori convergence rates of approximation schemes for feedback controllers for a class of distributed parameter systems is addressed within the context of hereditary schemes. Specifically, a factorization approach is presented for deriving approximations to the optimal feedback gains for the linear regulator-quadratic cost problem associated with time-varying functional differential equations with control delays. The approach is based on a discretization of the state penalty which leads to a simple structure for the feedback control law. General properties of the Volterra factors of Hilbert-Schmidt operators are then used to obtain convergence results for the controls, trajectories and feedback kernels. Two algorithms are derived from the basic approximation scheme, including a fast algorithm, in the time-invariant case. A numerical example is also considered.

Milman, Mark H.↗

A dual reciprocal boundary element formulation for viscous flows

The advantages inherent in the boundary element method (BEM) for potential flows are exploited to solve viscous flow problems. The trick is the introduction of a so-called dual reciprocal technique in which the convective terms are represented by a global function whose unknown coefficients are determined by collocation. The approach, which is necessarily iterative, converts the governing partial differential equations into integral equations via the distribution of fictitious sources or dipoles of unknown strength on the boundary. These integral equations consist of two parts. The first is a boundary integral term, whose kernel is the unknown strength of the fictitious sources and the fundamental solution of a convection-free flow problem. The second part is a domain integral term whose kernel is the convective portion of the governing PDEs. The domain integration can be transformed to the boundary by using the dual reciprocal (DR) concept. The resulting formulation is a pure boundary integral computational process.

Lafe, Olu↗

Nonparametric maximum likelihood estimation of probability densities by penalty function methods

When it is known a priori exactly to which finite dimensional manifold the probability density function gives rise to a set of samples, the parametric maximum likelihood estimation procedure leads to poor estimates and is unstable; while the nonparametric maximum likelihood procedure is undefined. A very general theory of maximum penalized likelihood estimation which should avoid many of these difficulties is presented. It is demonstrated that each reproducing kernel Hilbert space leads, in a very natural way, to a maximum penalized likelihood estimator and that a well-known class of reproducing kernel Hilbert spaces gives polynomial splines as the nonparametric maximum penalized likelihood estimates.

Demontricher, G. F.↗

VLA observations of impulsive solar flares at 4.9 GHz

A series of snapshot synthesis maps constructed as a function of time provided a complete mapping of source structure in one dimension and absolute location in two dimensions of five impulsive solar microwave bursts observed during July 15-18, 1978. It was found that all five flares were dominated by a single 12-18 arcsec source located between the H-alpha kernels, near the magnetic neutral line. Spectral data from Sagamore Hill (1.4-35 GHz) show that the spectral turnover at low frequencies results from the absorption by the ambient active region; these results, together with the observed magnetic configuration, suggest a model which could explain the evolution in source structure at 4.9 GHz in terms of magnetic reconnection which starts below the optical depth of the ambient active region equal to 1.

Marsh, K. A.↗

Error analysis for the ground-based microwave ozone measurements during STOIC

We present a formal error analysis and characterization of the microwave measurements made during the Stratospheric Ozone Intercomparison Campaign (STOIC). The most important error sources are found to be determination of the tropospheric opacity, the pressure-broadening coefficient of the observed line, and systematic variations in instrument response as a function of frequency ('baseline'). Net precision is 4-6% between 55 and 0.2 mbar, while accuracy is 6-10%. Resolution is 8-10 km below 3 mbar and increases to 17km at 0.2 mbar. We show the 'blind' microwave measurements from STOIC and make limited comparisons to other measurements. We use the averaging kernels of the microwave measurement to eliminate resolution and a priori effects from a comparison to SAGE 2. The STOIC results and comparisons are broadly consistent with the formal analysis.

Connor, Brian J.↗

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated. A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗

General theory of spherically symmetric boundary-value problems of the linear transport theory.

A general theory of spherically symmetric boundary-value problems of the one-speed neutron transport theory is presented. The formulation is also applicable to the 'gray' problems of radiative transfer. The Green's function for the purely absorbing medium is utilized in obtaining the normal mode expansion of the angular densities for both interior and exterior problems. As the integral equations for unknown coefficients are regular, a general class of reduction operators is introduced to reduce such regular integral equations to singular ones with a Cauchy-type kernel. Such operators then permit one to solve the singular integral equations by the standard techniques due to Muskhelishvili. We discuss several spherically symmetric problems. However, the treatment is kept sufficiently general to deal with problems lacking azimuthal symmetry. In particular the procedure seems to work for regions whose boundary coincides with one of the coordinate surfaces for which the Helmholtz equation is separable.

Kanal, M.↗

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Ott, L.↗

Dynamic topography and gravity anomalies for fluid layers whose viscosity varies exponentially with depth

Adopting the formalism of Parsons and Daly (1983), analytical integral equations (Green's function integrals) are derived which relate gravity anomalies and dynamic boundary topography with temperature as a function of wavenumber for a fluid layer whose viscosity varies exponentially with depth. In the earth, such a viscosity profile may be found in the asthenosphere, where the large thermal gradient leads to exponential decrease of viscosity with depth, the effects of a pressure increase being small in comparison. It is shown that, when viscosity varies rapidly, topography kernels for both the surface and bottom boundaries (and hence the gravity kernel) are strongly affected at all wavelengths.

Revenaugh, Justin↗

Ultraviolet converter transients induced by electrons

The output of ultraviolet converters typically used in satellite astronomy was monitored during irradiation with electrons from a sealed SR-90 source which approximated the peak flux in earth's outer electron belt. The signal induced by irradiation was attributed to two mechanisms: (1) photoelectrons resulting from photons created in the MgF2 window and (2) the direct impact of electrons on the phosphor. For irradiation at about 1 x 10 to the 7th e/sq cm sec, these two effects produced signals which were, in order of magnitude, the same as those produced by an incident UV flux (254 nm) of 10 to the 8th and 10 to the 7th photons/sq cm sec, respectively. In addition, the induced signal was investigated as a function of electron energy by irradiating another converter with 0.4-1.8-MeV electrons from a Van de Graaff. These results suggest that the dominant contribution to the electron-induced signal is Cerenkov photon production in the MgF2 window.

Kernell, R. L.↗

Virtual Sensors: Using Data Mining to Efficiently Estimate Spectra

Detecting clouds within a satellite image is essential for retrieving surface geophysical parameters, such as albedo and temperature, from optical and thermal imagery because the retrieval methods tend to be valid for clear skies only. Thus, routine satellite data processing requires reliable automated cloud detection algorithms that are applicable to many surface types. Unfortunately, cloud detection over snow and ice is difficult due to the lack of spectral contrast between clouds and snow. Snow and clouds are both highly reflective in the visible wavelen,ats and often show little contrast in the thermal Infrared. However, at 1.6 microns, the spectral signatures of snow and clouds differ enough to allow improved snow/ice/cloud discrimination. The recent Terra and Aqua Moderate Resolution Imaging Spectro-Radiometer (MODIS) sensors have a channel (channel 6) at 1.6 microns. Presently the most comprehensive, long-term information on surface albedo and temperature over snow- and ice-covered surfaces comes from the Advanced Very High Resolution Radiometer ( AVHRR) sensor that has been providing imagery since July 1981. The earlier AVHRR sensors (e.g. AVHRR/2) did not however have a channel designed for discriminating clouds from snow, such as the 1.6 micron channel available on the more recent AVHRR/3 or the MODIS sensors. In the absence of the 1.6 micron channel, the AVHRR Polar Pathfinder (APP) product performs cloud detection using a combination of time-series analysis and multispectral threshold tests based on the satellite's measuring channels to produce a cloud mask. The method has been found to work reasonably well over sea ice, but not so well over the ice sheets. Thus, improving the cloud mask in the APP dataset would be extremely helpful toward increasing the accuracy of the albedo and temperature retrievals, as well as extending the time-series of albedo and temperature retrievals from the more recent sensors to the historical ones. In this work, we use data mining methods to construct a model of MODIS channel 6 as a function of other channels that are common to both MODIS and AVHRR. The idea is to use the model to generate the equivalent of MODIS channel 6 for AVHRR as a function of the AVHRR equivalents to MODIS channels. We call this a Virtual Sensor because it predicts unmeasured spectra. The goal is to use this virtual channel 6. to yield a cloud mask superior to what is currently used in APP . Our results show that several data mining methods such as multilayer perceptrons (MLPs), ensemble methods (e.g., bagging), and kernel methods (e.g., support vector machines) generate channel 6 for unseen MODIS images with high accuracy. Because the true channel 6 is not available for AVHRR images, we qualitatively assess the virtual channel 6 for several AVHRR images.

Srivastava, Ashok↗

Advanced Free Flight Planner and Dispatcher's Workstation: Preliminary Design Specification

The National Aeronautics and Space Administration (NASA) has implemented the Advanced Air Transportation Technology (AATT) program to investigate future improvements to the national and international air traffic management systems. This research, as part of the AATT program, developed preliminary design requirements for an advanced Airline Operations Control (AOC) dispatcher's workstation, with emphasis on flight planning. This design will support the implementation of an experimental workstation in NASA laboratories that would emulate AOC dispatch operations. The work developed an airline flight plan data base and specified requirements for: a computer tool for generation and evaluation of free flight, user preferred trajectories (UPT); the kernel of an advanced flight planning system to be incorporated into the UPT-generation tool; and an AOC workstation to house the UPT-generation tool and to provide a real-time testing environment. A prototype for the advanced flight plan optimization kernel was developed and demonstrated. The flight planner uses dynamic programming to search a four-dimensional wind and temperature grid to identify the optimal route, altitude and speed for successive segments of a flight. An iterative process is employed in which a series of trajectories are successively refined until the LTPT is identified. The flight planner is designed to function in the current operational environment as well as in free flight. The free flight environment would enable greater flexibility in UPT selection based on alleviation of current procedural constraints. The prototype also takes advantage of advanced computer processing capabilities to implement more powerful optimization routines than would be possible with older computer systems.

Wilson, J.↗

Diffusion Flame Stabilization

Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm) (Fig. 1, left). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated (Fig. 1, right). A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.

Takahashi, Fumiaki↗

Activity-Centric Approach to Distributed Programming

The first phase of an effort to develop a NASA version of the Cybele software system has been completed. To give meaning to even a highly abbreviated summary of the modifications to be embodied in the NASA version, it is necessary to present the following background information on Cybele: Cybele is a proprietary software infrastructure for use by programmers in developing agent-based application programs [complex application programs that contain autonomous, interacting components (agents)]. Cybele provides support for event handling from multiple sources, multithreading, concurrency control, migration, and load balancing. A Cybele agent follows a programming paradigm, called activity-centric programming, that enables an abstraction over system-level thread mechanisms. Activity centric programming relieves application programmers of the complex tasks of thread management, concurrency control, and event management. In order to provide such functionality, activity-centric programming demands support of other layers of software. This concludes the background information. In the first phase of the present development, a new architecture for Cybele was defined. In this architecture, Cybele follows a modular service-based approach to coupling of the programming and service layers of software architecture. In a service-based approach, the functionalities supported by activity-centric programming are apportioned, according to their characteristics, among several groups called services. A well-defined interface among all such services serves as a path that facilitates the maintenance and enhancement of such services without adverse effect on the whole software framework. The activity-centric application-program interface (API) is part of a kernel. The kernel API calls the services by use of their published interface. This approach makes it possible for any application code written exclusively under the API to be portable for any configuration of Cybele.

Levy, Renato↗