Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational discretization”

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 145 records · Page 8

Composite propellant combustion modeling with a porous plate burner

A burner is designed to model on a large scale (millimeters) the complex vapor phase processes that take place on awkwardly small scales (about a hundred microns) in the combustion of practical AP/ composite propellants at conventional pressures. Binder vapor evolution is modeled with gaseous fuel (ethane in the experiments reported) flow through a porous plate and oxidizer vapor with the flow of a gaseous oxidizer (air and enriched air in the experiments reported) through discrete holes in the porous plate. Measured flame standoff distance and surface temperature variations are consistent with theoretical predictions at atmospheric pressure. Data obtained at several atmospheres are discussed in the light of the current theories of propellant burning that place varying emphasis on the roles of chemical kinetics and fluid dynamic diffusion/mixing in the vapor phase of a burning composite propellant. The potential and proposed future applications of the burner are indicated.

Kumar, R. N.↗

Application of spatial features to satellite land-use analysis

A Level I land-use analysis of selected training areas of the Colorado Front Range was carried out using digital ERTS-A satellite imagery. Level I land-use categories included urban, agriculture (irrigated and dryland farming), rangeland, and forests. The spatial variations in spectral response for these land-use classes were analyzed using discrete two-dimensional Fourier transforms to isolate and extract spatial features. Analysis was performed on ERTS frame 1352-17134 (July 10, 1973) and frame number 1388-17131 (August 15, 1973). On training sets, spatial features yielded 80 to 100 percent classification accuracies with commission errors ranging from 0 to 20 percent.

Smith, J.↗

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

Global distribution and migration of subsurface ice on Mars

An analytical model is developed for the factors which changed the initial distribution and state of H2O, mainly trapped as subsurface ice, on Mars. A uniform initial ice distribution is assumed, and changes in ice layer thickness are calculated at discrete latitudes of Mars. The thermal/diffusive model used for the calculations accounts for variations in the solar luminosity, the planetary obliquity and eccentricity, with obliquity starting 3.5 Gyr ago, and insolation at various latitudes. Attention is given to variations in the albedo due to seasonal CO2 variations, in the regolith temperature, and in atmospheric pressure. Migration of H2O into and out of the regolith is a function of the H2O vapor temperature at the subsurface ice boundary and the annual average H2O at the base of the atmosphere. Projections of the ice thickness over 100 obliquity periods indicate that the intensity of the obliquity effects are latitudinally dependent. H2O atmospheric concentrations are projected to vary up to six orders of magnitude over an obliquity cycle.

Fanale, F. P.↗

On the nonlinear development of the most unstable Goertler vortex mode

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

On the nonlinear development of the most unstable Goertler vortex mode TLSP - Final report

The nonlinear development of the most unstable Goertler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Goertler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Goertler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Goertler equations valid for O(1) Goertler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

High dimensional similarity search with quantum assisted variational autoencoder

Recent progress in quantum algorithms and hardware is indicator of the potential importance of quantum computing in the next future. However, finding suitable application areas remains an active area of research. Quantum machine learning [1] is touted as a potential approach to demonstrate quantum advantage within both the gate-model [2,3] and the adiabatic [4,5] schemes. For instance, the Quantum-assisted Variational Autoencoder (QVAE) [6] has been proposed as a quantum enhancement to the discrete VAE [7]. We extend on previous work and study the real-world applicability of a QVAE, specifically, for similarity search in large-scale high dimensional datasets. While similarity search algorithms are available for low dimensional datasets, scaling to billion-scale datasets with thousands of dimensions is non-trivial. We show how the latent-space representation of a QVAE can be used to construct a space-efficient search index. We back up our claims by experimental results which show a correlation between the Hamming distance in the embedded space and the Euclidean distance in the original space on the Moderate Resolution Imaging Spectroradiometer (MODIS) dataset. Further, we show real-world speedups compared to linear search and demonstrate memory efficient scaling to large-scale datasets.

Nicholas D Gao↗

Stress analysis and stress intensity factors for finite geometry solids containing rectangular surface cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement field in a finite geometry bar containing a variable depth rectangular surface crack under extensionally applied uniform loading. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. Using the obtained displacement field, normal stresses and the stress intensity factor variation along the crack periphery are calculated for different crack depth to bar thickness ratios. Crack opening displacements and stress intensity factors are also obtained for a through-thickness, center cracked bar with variable thickness. The reported results show a considerable potential for using this method in calculating stress intensity factors for commonly encountered surface crack geometries in finite solids.

Gyekenyesi, J. P.↗

Interim prediction method for fan and compressor source noise

A method is presented for interim use in assessing the noise generated by fans and compressors in turbojet and turbofan engines. One-third octave band sound pressure levels consisting of broadband, discrete tone, and combination-tone noise components are predicted. Spectral distributions and directivity variations are specified. The method is based on that developed by other investigators with modifications derived from an analysis of full-scale, single-stage fan data. Comparisons of predicted and measured noise performance are presented, and requirements for improving the method are discussed.

Heidmann, M. F.↗

GCR Simulator Reference Field and a Spectral Approach for Laboratory Simulation

The galactic cosmic ray (GCR) simulator at the NASA Space Radiation Laboratory (NSRL) is intended to deliver the broad spectrum of particles and energies encountered in deep space to biological targets in a controlled laboratory setting. In this work, certain aspects of simulating the GCR environment in the laboratory are discussed. Reference field specification and beam selection strategies at NSRL are the main focus, but the analysis presented herein may be modified for other facilities. First, comparisons are made between direct simulation of the external, free space GCR field and simulation of the induced tissue field behind shielding. It is found that upper energy constraints at NSRL limit the ability to simulate the external, free space field directly (i.e. shielding placed in the beam line in front of a biological target and exposed to a free space spectrum). Second, variation in the induced tissue field associated with shielding configuration and solar activity is addressed. It is found that the observed variation is likely within the uncertainty associated with representing any GCR reference field with discrete ion beams in the laboratory, given current facility constraints. A single reference field for deep space missions is subsequently identified. Third, an approach for selecting beams at NSRL to simulate the designated reference field is presented. Drawbacks of the proposed methodology are discussed and weighed against alternative simulation strategies. The neutron component and track structure characteristics of the simulated field are discussed in this context.

Slaba, Tony C.↗

Property enchancement of polyimide films by way of the incorporation of lanthanide metal ions

Lanthanide metal ions were incorporated into the polyimide derived from 2,2-bis(3,4-dicarboxyphenyl) hexafluoropropane dianhydride (6FDA) and 1,3-bis(aminophenoxy) benzene (APB) in an attempt to produce molecular level metal-polymer composites. The lanthanide series of metal ions (including aluminum, scandium, and yttrium) provide discrete and stable metal ions in the 3+ oxidation state. Throughout the series there is a uniform variation in ionic size ranging from 50 pm for aluminum to a maximum of 103.4 pm for cerium and gradually decreasing again to 84.8 pm for lutetium. The high charge-to-size ratio for these ions as well as the ability to obtain large coordination numbers makes them excellent candidates for interacting with the polymer substructure. The distinct lack of solubility of simple lanthanide salts such as the acetates and halides has made it difficult to obtain metal ions distributed in the polymer framework as discrete ions or metal complexes rather than microcomposites of metal clusters. (Lanthanum nitrates are quite soluble, but the presence of the strongly oxidizing nitrate ion leads to serious degradation of the polymer upon thermal curing. This work was successful at extending the range of soluble metals salts by using chelating agents derived from the beta-diketones dipivaloylmethane, dibenzoylmethane, trifluoroacetylacetone, and hexafluoroacetylacetone. Metal acetates which are insoluble in dimethylacetamide dissolve readily in the presence of the diketones. Addition of the polyimide yields a homogeneous resin which is then cast into a clear film. Upon curing clear films were obtained with the dibenzoylmethane and trifluoroacetylacetone ligands. The dipavaloylmethane precipitates the metal during the film casting process, and hexafluoroacetylacetone gives cured films which are deformed and brittle. These clear films are being evaluated for the effect of the metal ions on the coefficient of thermal expansion, resistance to atomic oxygen, and on selective gas permeability. Much more commonly than above, polyimide films are prepared by casting the film as the poly(amic acid) precursor which is then converted to the imidized form during the thermal cure cycle. Very limited success was achieved in the past in adding lanthanide metal ions to the amide precursors because of gellation and lack of solubility. With the use of the diketone ligands cited above, the solubility and gellation problems were overcome. However, the films after curing were clear but unacceptably brittle. Attempts to overcome this cure embrittlement problem are in progress.

Thompson, David W.↗

The Use of RGPS Kinematic Data to Estimate Nonlinear Sea Ice Motion

In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motion In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motio

Pruis, M.↗

Stress analysis and stress-intensity factors for finite geometry solids containing rectangular surface cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement field in a finite geometry bar containing a variable depth rectangular surface crack under extensionally applied uniform loading. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. Using the obtained displacement field, normal stresses, and the stress-intensity factor variation along the crack periphery are calculated for different crack depth to bar thickness ratios. Crack opening displacements and stress-intensity factors are also obtained for a through-thickness, center-cracked bar with variable thickness. The reported results show a considerable potential for using this method in calculating stress-intensity factors for commonly encountered surface crack geometries in finite solids

Gyekenyesi, J. P.↗

The photometric period of the recurrent Nova T Pyxidis

The paper presents 1713 photometric measurements of T Pyx from 1966 to 1990. The light curve of T Pyx shows roughly sinusoidal variations with a typical amplitude of 0.09 mag and a time scale of 2 hr. Discrete Fourier transforms of data from individual nights and runs reveal a highly significant periodic modulation. This modulation is proven to be coherent only on time scales shorter than from 6 d to under 1 d. The present period for the modulation is 0.07616 +/- 0.00017 d. Since other stars have variable photometric periods (possibly related to the superhump phenomenon) which are slightly different from the orbital period, an orbital period of about 0.073 d for T Pyx is suggested.

Schaefer, Bradley E.↗

Seasonal and Non-Seasonal Variations of Jupiter's Atmosphere from Observations of Thermal Emission, 1994-2011

We analyzed mid-infrared images of Jupiter's thermal emission, covering approx.1.5 Jovian years, acquired in discrete filters between 7.8 and 24.5 microns. The behavior of stratospheric (approx.10-mbar) and tropospheric (approx.100-400 mbar) temperatures is generally consistent with predictions of seasonal variability, with differences between 100-mbar temperatures +/-50-60deg from the equator on the order of +/-2. Removing this effect, there appear to be long-term quasi-periodic variability of tropospheric temperatures, whose amplitude, phase and period depend on latitude. The behavior of temperatures in the Equatorial Zone (EZ) suggests a approx.4-6-year period with amplitude of about +/-1-1.5 K in temperature. At mid-latitudes, the periodicity is more distinct with amplitudes around +/-1.5-2.5 K and 4-8 year periods. The 4.2-year variation of stratospheric temperatures known as the quasiquadrennial oscillation or "QQO" (Leovy et al. 1991, Nature 354, 380) continued during this period. There were no variations of zonal mean temperatures associated with any of the "global upheaval" events that have produced dramatic changes of jupiter's visible appearance and cloud cover, although there are colder discrete regions associated with updrafts, e.g. the early stages of the re-darkening ("revival") of the South Equatorial Belt (SEB) in late 2010. On the other hand increases in the visible albedos ("fades") of belts are accompanied by increases in the thickness of a 700-mbar cloud layer (most likely NH3 ice) and clouds at higher pressures, together with the mixing ratio of NH3 gas near 400 mbar (above its condensation level). These quantities decrease during re-darkening ("revival") episodes, during which we note discrete features that are exceptions to the general correlation between dark albedos and minimal cloudiness. In contrast to all these changes, the meridional distribution of the 240-mbar para-H2 fraction appears to be invariant in time.

Orton, G.↗

Studies in a transonic rotor aerodynamics and noise facility

The design, construction and testing of a transonic rotor aerodynamics and noise facility was undertaken, using a rotating arm blade element support technique. This approach provides a research capability intermediate between that of a stationary element in a moving flow and that of a complete rotating blade system, and permits the acoustic properties of blade tip elements to be studied in isolation. This approach is an inexpensive means of obtaining data at high subsonic and transonic tip speeds on the effect of variations in tip geometry. The facility may be suitable for research on broad band noise and discrete noise in addition to high-speed noise. Initial tests were conducted over the Mach number range 0.3 to 0.93 and confirmed the adequacy of the acoustic treatment used in the facility to avoid reflection from the enclosure.

Wright, S. E.↗

Instability Mechanisms of Thermally-Driven Interfacial Flows in Liquid-Encapsulated Crystal Growth

During the past year, a great deal of effort was focused on the enhancement and refinement of the computational tools developed as part of our previous NASA grant. In particular, the interface mollification algorithm developed earlier was extended to incorporate the effects of surface-rheological properties in order to allow the study of thermocapillary flows in the presence of surface contamination. These tools will be used in the computational component of the proposed research in the remaining years of this grant. A detailed description of the progress made in this area is provided elsewhere. Briefly, the method developed allows for the convection and diffusion of bulk-insoluble surfactants on a moving and deforming interface. The novelty of the method is its grid independence: there is no need for front tracking, surface reconstruction, body-fitted grid generation, or metric evaluations; these are all very expensive computational tasks in three dimensions. For small local radii of curvature there is a need for local grid adaption so that the smearing thickness remains a small fraction of the radius of curvature. A special Neumann boundary condition was devised and applied so that the calculated surfactant concentration has no variations normal to the interface, and it is hence truly a surface-defined quantity. The discretized governing equations are solved subsequently using a time-split integration scheme which updates the concentration and the shape successively. Results demonstrate excellent agreement between the computed and exact solutions.

Haj-Hariri, Hossein↗

A Variational Principle for Reconstruction of Elastic Deformations in Shear Deformable Plates and Shells

A variational principle is formulated for the inverse problem of full-field reconstruction of three-dimensional plate/shell deformations from experimentally measured surface strains. The formulation is based upon the minimization of a least squares functional that uses the complete set of strain measures consistent with linear, first-order shear-deformation theory. The formulation, which accommodates for transverse shear deformation, is applicable for the analysis of thin and moderately thick plate and shell structures. The main benefit of the variational principle is that it is well suited for C(sup 0)-continuous displacement finite element discretizations, thus enabling the development of robust algorithms for application to complex civil and aeronautical structures. The methodology is especially aimed at the next generation of aerospace vehicles for use in real-time structural health monitoring systems.

Tessler, Alexander↗