Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient 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 181 records · Page 10

Lyman-alpha clouds as a relic of primordial density fluctuations

Primordial density fluctuations are studied using a CDM model and primordial clouds some of which are expanding, driven by pressure gradients created when the medium is photionized, and some of which are massive enough to continue collapsing in spite of the pressure. Normalization of CDM models to the clustering properties on large scales are used to predict the parameters of collapsing clouds of subgalactic mass at early epochs. It is shown that the abundance and dimensions of these clouds are comparable to those of the Lyman-alpha systems. The evolutionary history of the clouds is computed, utilizing a spherically symmetric hydrodynamics code with the dark matter treated as a collisionless fluid, and the H I column density distribution is evaluated as a function of N(H I) and redshift. The observed cloud parameters come out naturally in the CDM model and suggest that Lyman-alpha clouds are the missing link between primordial density fluctuations and the formation of galaxies.

Bond, J. R.↗

Geomagnetic Field Distortion by a Solar Stream as a Mechanism for the Production of Polar Aurora and Electrojets

This paper describes a mechanism for charge separation in the geomagnetically trapped radiation which may account for some observed phenomena associated with the polar aurora and the electrojet current systems. The following development is proposed: given that there exist eastward or westward longitudinal gradients in the geomagnetic field resulting from distortion of the geomagnetic field by solar streams, if the trapped radiation is adiabatic in character, radial drift separation of positive and negative charged particles must occur. It follows that, for bounded or irregular distributions of plasma number density in such an adiabatic - drift region, electric fields will arise. The origin of such electric fields will not arrest the drift separation of the charged particles, but will contribute to exponential growth of irregularities in the trapped plasma density. An adiabatic acceleration mechanism is described, which is based on incorporating the electrostatic energy of the particle in the energy function for the particle. Direct consequences of polarization of the geomagnetically trapped radiation will be the polar electrojet current systems and the polar aurora.

Kern, J. W.↗

Effect of Interface Shape and Magnetic Field on the Microstructure of Bulk Ge:Ga

Thermal and compositional gradients induced during the growth process contribute significantly to the development of defects in the solidified boule. Thermal gradients and the solid-liquid interface shape can be greatly effected by ampoule material. Compositional gradients are strongly influenced by interface curvature and convective flow in the liquid. Results of this investigation illustrate the combined influences of interface shape and convective fluid flow. An applied magnetic field was used to reduce the effects of convective fluid flow in the electrically conductive melt during directional solidification. Several 8 mm diameter boules of Ga-doped Ge were grown at different field strengths, up to 5 Tesla, in four different ampoule materials. Compositional profiles indicate mass transfer conditions ranged from completely mixed to diffusion controlled. The influence of convection in the melt on the developing crystal microstructure and defect density was investigated as a function of field strength and ampoule material. Chemical etching and electron backscattered electron diffraction were used to map the crystal structure of each boule along the center plane. Dislocation etch pit densities were measured for each boule. Results show the influence of magnetic field strength and ampoule material on overall crystal quality.

Cobb, S. D.↗

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗

Fast Aircraft Separation Calculations for Gradient Based Optimization of Airspace Simulations

Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.

Optimization↗

Stellar occultation spikes as probes of atmospheric structure and composition

The characteristics of spikes observed in occultation light curves of Beta Scorpii by Jupiter are discussed in terms of the gravity-gradient model. The occultation of Beta Sco by Jupiter on May 13, 1971, is reviewed, and the gravity-gradient model is defined as an isothermal atmosphere of constant composition in which the refractivity is a function only of the radial coordinate from the center of refraction, which is assumed to lie parallel to the local gravity gradient. The derivation of the occultation light curve in terms of the atmosphere, the angular diameter of the occulted star, and the occultation geometry is outlined. It is shown that analysis of the light-curve spikes can yield the He/H2 concentration ratio in a well-mixed atmosphere, information on fine-scale atmospheric structure, high-resolution images of the occulted star, and information on ray crossing. Observational limits are placed on the magnitude of horizontal refractivity gradients, and it is concluded that the spikes are the result of local atmospheric density variations: atmospheric layers, density waves, or turbulence.

Elliot, J. L.↗

The radial gradients and collisional properties of solar wind electrons

The plasma instrument on Mariner 10 carried out measurements of electron density and temperature in the interplanetary medium between heliocentric distances of 0.85 and 0.45 AU. Due to the stable coronal configuration and low solar activity during the period of observation, the radial variations of these quantities could be obtained. The power-law exponent of the core temperature was measured to be -0.3 + or - 0.04, and the halo temperature was found to be almost independent of heliocentric distance. The exponent of the power law for the density variation was 2.5 + or - 0.2 and the extrapolated value at 1 AU was consistent with measured values during the same period. Calculations of the core electron self-collision time, and the core-halo equipartition time were made as a function of radial distance. These measurements indicate a macroscale picture of a Coulomb-collisional core and a collisionless isothermal halo. Extrapolating back to the sun, core and halo temperatures become equal at a radial distance of approx. 2-15 radii.

Gilvie, K. W.↗

Diffusion Driven Combustion Waves in Porous Media

Filtration of gas containing oxidizer, to the reaction zone in a porous medium, due, e.g., to a buoyancy force or to an external pressure gradient, leads to the propagation of Filtration combustion (FC) waves. The exothermic reaction occurs between the fuel component of the solid matrix and the oxidizer. In this paper, we analyze the ability of a reaction wave to propagate in a porous medium without the aid of filtration. We find that one possible mechanism of propagation is that the wave is driven by diffusion of oxidizer from the environment. The solution of the combustion problem describing diffusion driven waves is similar to the solution of the Stefan problem describing the propagation of phase transition waves, in that the temperature on the interface between the burned and unburned regions is constant, the combustion wave is described by a similarity solution which is a function of the similarity variable x/square root of(t) and the wave velocity decays as 1/square root of(t). The difference between the two problems is that in the combustion problem the temperature is not prescribed, but rather, is determined as part of the solution. We will show that the length of samples in which such self-sustained combustion waves can occur, must exceed a critical value which strongly depends on the combustion temperature T(sub b). Smaller values of T(sub b) require longer sample lengths for diffusion driven combustion waves to exist. Because of their relatively small velocity, diffusion driven waves are considered to be relevant for the case of low heat losses, which occur for large diameter samples or in microgravity conditions, Another possible mechanism of porous medium combustion describes waves which propagate by consuming the oxidizer initially stored in the pores of the sample. This occurs for abnormally high pressure and gas density. In this case, uniformly propagating planar waves, which are kinetically controlled, can propagate, Diffusion of oxidizer decreases the wave velocity. In addition to the reaction and diffusion layers, the uniformly propagating wave structure includes a layer with a pressure gradient, where the gas motion is induced by the production or consumption of the gas in the reaction as well as by thermal expansion of the gas. The width of this zone determines the scale of the combustion wave in the porous medium.

Aldushin, A. P.↗

Ballistic Performance Model of Crater Formation in Monolithic, Porous Thermal Protection Systems

Porous monolithic ablative systems insulate atmospheric reentry vehicles from reentry plasmas generated by atmospheric braking from orbital and exo-orbital velocities. Due to the necessity that these materials create a temperature gradient up to several thousand Kelvin over their thickness, it is important that these materials are near their pristine state prior to reentry. These materials may also be on exposed surfaces to space environment threats like orbital debris and meteoroids leaving a probability that these exposed surfaces will be below their prescribed values. Owing to the typical small size of impact craters in these materials, the local flow fields over these craters and the ablative process afford some margin in thermal protection designs for these locally reduced performance values. In this work, tests to develop ballistic performance models for thermal protection materials typical of those being used on Orion are discussed. A density profile as a function of depth of a typical monolithic ablator and substructure system is shown in Figure 1a.

Miller, J. E.↗

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗

Dynamic upper atmospheric force model on stabilized vehicles for a high-precision trajectory computer program

The upper atmosphere model draws heavily on the behavior of the earth's upper atmosphere which exhibits cyclic as well as irregular variations in density profile, temperature, pressure, and composition in unison with solar activities as deduced from the more recent land-based and satellite observations. The lift and drag model is designed specifically for inertially stabilized vehicles of the Mariner class, with possible extension to gravity gradient stabilized vehicles of the GEOS class. The model considers operation in the free molecular flow regimes with large Knudsen numbers. The vehicle is considered a composite structure with basic components having well-defined shapes, each with its own surface characteristics in terms of temperature, reflectivity, and accommodation of free stream molecules. The model takes into account both the calculation of precise aerodynamic force coefficients in terms of expansion of modified Bessel functions in speed ratios and angle of attack, and approximate force coefficients when the speed ratios approach infinity. Other considerations include specular and diffused reflectivity, shielding, and shadow effects.

Khatib, A. R.↗

Quantitative Infrared Image Analysis Of Thermally-Thin Cellulose Surface Temperatures During Upstream and Downstream Microgravity Flame Spread from A Central Ignition Line

Surface view calibrated infrared images of ignition and flame spread over a thin cellulose fuel were obtained at 30 Hz during microgravity flame spread tests in the 10 second Japan Microgravity Center (JAMIC). The tests also used a color video of the surface view and color images of the edge view using 35 millimeter 1600 Kodak Ektapress film at 2 Hz. The cellulose fuel samples (50% long fibers from lumi pine and 50% short fibers from birch) were made with an area density of 60 grams per square meters. The samples were mounted in the center of a 12 centimeter wide by 16 centimeter tall flow duct that uses a downstream fan to draw the air through the flow duct. Samples were ignited after the experiment package was released using a straight hot wire across the center of the 7.5 centimeter wide by 14 centimeter long samples. One case, at 1 atmosphere 35%O2 in N2, at a forced flow of 10 centimeters per second, is presented here. In this case, as the test progresses, the single flame begins to separate into simultaneous upstream and downstream flames. Surface temperature profiles are evaluated as a function of time, and temperature gradients for upstream and downstream flame spread are measured. Flame spread rates from IR image data are compared to visible image spread rate data. IR blackbody temperatures are compared to surface thermocouple readings to evaluate the effective emissivity of the pyrolyzing surface. Preheat lengths are evaluated both upstream and downstream of the central ignition point. A surface energy balance estimates the net heat flux from the flame to the fuel surface along the length of the fuel.

Olson, Sandra L.↗

Growth hormone aggregates in the rat adenohypophysis

Although it has been known for some time that GH aggregates are contained within the rat anterior pituitary gland, the role that they might play in pituitary function is unknown. The present study examines this issue using the technique of Western blotting, which permitted visualization of 11 GH variants with apparent mol wt ranging from 14-88K. Electroelution of the higher mol wt variants from gels followed by their chemical reduction with beta-mercaptoethanol increased GH immunoassayability by about 5-fold. With the blot procedure we found 1) that GH aggregates greater than 44K were associated with a 40,000 x g sedimentable fraction; 2) that GH aggregates were not present in glands from thyroidectomized rats, but were in glands from the thyroidectomized rats injected with T4; 3) that GH aggregates were uniquely associated with a heavily granulated somatotroph subpopulation isolated by density gradient centrifugation; and 4) that high mol wt GH forms were released from the dense somatotrophs in culture, since treatment of the culture medium with beta-mercaptoethanol increased GH immunoassayability by about 5-fold. Taken together, the results show that high mol wt GH aggregates are contained in secretory granules of certain somatotrophs and are also released in aggregate form from these cells in vitro.

NASA Discipline Cell Biology↗

A Chandra X-Ray Study of Cygnus A: The Cluster of Galaxies - 3

We present an analysis of the Chandra Advanced CCD Imaging Spectrometer (ACIS) observation of the intracluster gas associated with the cluster of galaxies surrounding Cygnus A. The dominant gaseous structure is a roughly elliptical (presumably prolate spheroidal in three dimensions) feature with semimajor axis approximately equal to 1.(prime)1 (approximately equal to 100 kpc). This structure apparently represents intracluster gas that has been swept up and compressed by a cavity inflated in this gas by relativistic material that has passed through the ends of the radio jets. The X-ray-emitting gas shows this prolate spheroidal morphology to approximately equal to 1.(prime)2 (110 kpc) from the radio galaxy but is spherical on larger scales. The X-ray emission from the intracluster gas extends to at least 8 prime (approximately equal to 720 kpc) from the radio galaxy, and a second, extended source of X-ray emission (probably associated with a second cluster of galaxies) is seen some 12 prime (approximately equal to 1 Mpc) to the northwest of Cygnus A. The X-ray spectrum of the integrated intracluster gas imaged on the S3 chip (dimensions 8 prime x 8 prime = 720 x 720 kpc), excluding the contribution from the radio galaxy and other compact sources of X-ray emission, has a gas temperature, metallicity, and unabsorbed 2-10 keV rest-frame luminosity of 7.7 keV, 0.34 times solar, and 3.5 x 10(exp 44) ergs per second, respectively. We have projected the X-ray spectra taken from 12 elliptical and circular annuli in order to derive a run of temperature, metallicity, density, and pressure as a function of radius. The temperature of the X-ray-emitting gas drops from approximately equal to 8 keV more than 100 kpc from the center to approximately equal to 5 keV some 80 kpc from the center, with the coolest gas immediately adjacent to the radio galaxy. "Belts" of slightly cooler (approximately equal to 4 keV) X-ray-emitting gas run around the minor dimension of the cavity created by the radio source, while the limb-brightened edges of the cavity are slightly hotter (approximately equal to 6 keV), perhaps as a result of heating by a bow shock driven by the probably expanding cavity into the intracluster gas. There is a metallicity gradient in the X-ray-emitting gas, with the highest metallicities (approximately solar) found close to the center, decreasing to approximately 0.3 times solar in the outer parts. We have used the assumption of hydrostatic equilibrium to derive a total cluster mass within 500 kpc of 2.0 x 10(exp 14) solar mass and 2.8 x 10(exp 14) solar mass for constant and centrally decreasing temperature profiles, respectively. The total mass of X-ray-emitting gas within the same radius is 1.1 x 10(exp 13) solar mass. Thus, the gas fraction of the cluster within 500 kpc is 0.055 and 0.039 for the constant and centrally decreasing temperature profiles, respectively.

Smith, David A.↗

Gas phase abundances and conditions along the sight line to the low-halo, inner galaxy star HD 167756

We present high-resolution (3.5 km/s) Goddard High Resolution Spectrograph (GHRS) measurements of the Mg II, Si II, Cr II, Fe II, and Zn II lines toward HD 167756, a low-latitude halo star at a distance of 4 kpc in the direction l = 351.5 deg, b = -12.3 and at a Galactic altitude of z = -0.85 kpc. Supplemental Na I, Ca II, and H I data are also presented for comparison with the UV lines. Our analysis centers on converting the observed absoprtion-line data into measures of the apparent column density per unit velocity. N(sub a)(v), over the velocity range -25 less than or = v(sub lsr) less than 30 km/s for each species observed. We use these N(sub a)(v) profiles to construct logarithmic abundance ratios of Mg II, Si II, Cr II, Fe II, and Ca II relative to Zn II, normalized to solar abundances, as a function of velocity. Compared to Zn, these species show an underabundance relative to their solar values, with the largest underabundances occurring in the v(sub lsr) approximately equals 5 km/s component(s), for which we find logarithmic abundances A(sub Si/Zn) greater than -0.38, A(Mg/Zn) = -0.82, A(sub Cr/Zn) = -1.18, and A(sub Fe/Zn) greater than 1.40 dex. We show that ionization effects, abundance gradients, or intrinsic abundance variability cannot be significant sources for the underabundances observed. The most likely explanation is gas phase depletion of elements onto dust grains. Comparisons with the gas phase abundances along other diffuse, warm gas sight lines, like the halo sight line to HD 93521, support this interpretation as do the derived physical properties of the sight line.

Cardelli, Jason A.↗

SOHO/CDS Measurements of Coronal EUV Polarization above the Limb

Attempts to measure polarization in coronal EUV emission above the solar limb have been made using the SOHO/CDS normal-incidence spectrometer which has a polarization sensitivity of about 50%, a property that causes variations in intensity response as a function of the spacecraft's roll angle for polarized light. Such observations were made on the disk and up to 0.22 solar radii above the solar limb in a number of EUV lines during two special roll-maneuvers of the SOHO spacecraft. Measurements of intensity gradients were made above a modestly active equatorial region in 1997 and above a relatively cool polar region in 2001. Observed emission lines include He I 584A, He II 304A, 0 IV 555+610A, 0 V 630A, Mg IX 368A, Mg X 610+625A, and Si XI 303A, formed at temperatures that evenly cover the range in logT from 4.1 to 6.2. Near the disk, measured intensities of all lines fall off exponentially at different rates that can be used to determine the density scale-heights of the emitting plasma, since this emission is dominated by collisional excitation with an Ne-squared dependence. Assuming hydrostatic equilibrium, the intensity gradient for each line can then be converted into a 'scale-height temperature', which is found to be closely related to the ionization temperature of each line over the wide range of lines and solar conditions observed. Beyond a certain distance, intensity gradients of the cooler lines switch over to a flatter exponential slope, suggesting that this radiation is dominated by resonance scattering which varies as Ne to the first power. Such radiation should also be linearly polarized in the plane containing the line-of-sight and the solar center, a signature that would strongly confirm this interpretation.

Thomas, Roger J.↗

SOHO/CDS Measurements of Coronal EUV Polarization above the Limb

Attempts to measure polarization in coronal extreme ultraviolet (EUV) emission above the solar limb have been made using the SOHO/CDS normal-incidence spectrometer which has a polarization sensitivity of about 50%, a property that causes variations in intensity response as a function of the spacecraft's roll angle for polarized light. Such observations were made on the disk and up to 0.22 solar radii above the solar limb in a number of EUV lines during two special roll-maneuvers of the SOHO spacecraft. Measurements of intensity gradients were made above a modestly active equatorial region in 1997 and above a relatively cool polar region in 2001. Observed emission lines include He I 584 A, He II 304 A, O IV 555+610 A, O V 630 A, Mg IX 368 A, Mg X 610+625 A, and Si XI 303 A, formed at temperatures that evenly cover the range in log T from 4.1 to 6.2. Near the disk, measured intensities of all lines fall off exponentially at different rates that can be used to determine the density scale-heights of the emitting plasma, since this emission is dominated by collisional excitation with an Ne-squared dependence. Assuming hydrostatic equilibrium, the intensity gradient for each line can then be converted into a 'scale-height temperature', which is found to be closely related to the ionization temperature of each line over the wide range of lines and solar conditions observed. Thus the large-scale corona is remarkably uniform, even though clearly displaying a great deal of structure and non-uniformity on smaller spatial scales. Beyond a certain distance, intensity gradients of the cooler lines switch over to a flatter exponential slope, suggesting that this radiation is dominated by resonance scattering which varies as Ne to the first power. Such radiation should also be linearly polarized in the plane containing the line-of-sight and the solar center, a signature that would strongly confirm this interpretation.

Thomas, Roger J.↗

An Ensemble Approach to Building Mercer Kernels with Prior Information

This paper presents a new methodology for automatic knowledge driven data mining based on the theory of Mercer Kernels, which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly dimensional feature space. we describe a new method called Mixture Density Mercer Kernels to learn kernel function directly from data, rather than using pre-defined kernels. These data adaptive kernels can encode prior knowledge in the kernel using a Bayesian formulation, thus allowing for physical information to be encoded in the model. Specifically, we demonstrate the use of the algorithm in situations with extremely small samples of data. We compare the results with existing algorithms on data from the Sloan Digital Sky Survey (SDSS) and demonstrate the method's superior performance against standard methods. The code for these experiments has been generated with the AUTOBAYES tool, which automatically generates efficient and documented C/C++ code from abstract statistical model specifications. The core of the system is a schema library which contains templates for learning and knowledge discovery algorithms like different versions of EM, or numeric optimization methods like conjugate gradient methods. The template instantiation is supported by symbolic-algebraic computations, which allows AUTOBAYES to find closed-form solutions and, where possible, to integrate them into the code.

Srivastava, Ashok N.↗