Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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 721 records · Page 40

Exploratory Environmental Tests of Several Heat Shields

Exploratory tests have been conducted with several conceptual radiative heat shields of composite construction. Measured transient temperature distributions were obtained for a graphite heat shield without insulation and with three types of insulating materials, and for a metal multipost heat shield, at surface temperatures of approximately 2,000 F and 1,450 F, respectively, by use of a radiant-heat facility. The graphite configurations suffered loss of surface material under repeated irradiation. Temperature distribution calculated for the metal heat shield by a numerical procedure was in good agreement with measured data. Environmental survival tests of the graphite heat shield without insulation, an insulated multipost heat shield, and a stainless-steel-tile heat shield were made at temperatures of 2,000 F and dynamic pressures of approximately 6,000 lb/sq ft, provided by an ethylene-heated jet operating at a Mach number of 2.0 and sea-level conditions. The graphite heat shield survived the simulated aerodynamic heating and pressure loading. A problem area exists in the design and materials for heat-resistant fasteners between the graphite shield and the base structure. The insulated multipost heat shield was found to be superior to the stainless-steel-tile heat shield in retarding heat flow. Over-lapped face-plate joints and surface smoothness of the insulated multi- post heat shield were not adversely affected by the test environment. The graphite heat shield without insulation survived tests made in the acoustic environment of a large air jet. This acoustic environment is random in frequency and has an overall noise level of 160 decibels.

Goodman, George P.↗

A new approach to the solution of large, full matrix equations: A two-dimensional potential flow feasibility study

An approach to the solution of matrix problems resulting from integral equations of mathematical physics is presented. Based on the inherent smoothness in such equations, the problem is reformulated using a set of orthogonal basis vectors, leading to an equivalent coefficient problem which can be of lower order without significantly impairing the accuracy of the solution. This approach was evaluated using a two-dimensional Neumann problem describing the inviscid, incompressible flow over an airfoil. Two different kinds of mode functions were investigated, namely eigenfunction series and Fourier series. The method using Fourier series was found preferable. It uses all of the coefficients from a Fast Fourier Transform algorithm in an approximate method which exploits the known structure of the transformed coefficient matrix and very promising results for the flow over a realistic airfoil are obtained. On the basis of the results presented here, an order of magnitude reduction in this computer time can be expected for such problems as compared with the time for a direct matrix solution.

James, R. M.↗

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗

Three Dimensionally Interlinked, Dense, Solid Form of Single-Walled CNT Ropes

A 3D networked, dense form of single-walled carbon nanotubes (SWNT) has been made through isotropic shrinking of a gel-like SWNT-water paste by very slow evaporation. Approximately 35 g of Raw HiPco nanotubes were cleaned by the method of soft baking (250 C for 15 hours in air saturated with water vapor) in a glass beaker followed by leaching with concentrated hydrochloric acid. Typically, one liter of concentrated hydrochloric acid was added to the soft-baked voluminous mass in the same large beaker, and allowed to digest at room temperature with stirring overnight. The acid-digested SWNT slurry was filtered through a large porcelain Buchner funnel under atmospheric pressure. The slurry was continuously flushed, while still in the funnel, with a very slow but steady stream of deionized water employing a peristaltic pump. This process, referred to as gwashing, h continued until the filtrate water dripping from the Buchner funnel was clear, colorless, and neutral to a pH paper. This took about 15 liters of water to flow through the slurry over a day. At this point, the water pump was stopped and the SWNT-water slurry was allowed to drain the excess water for about 10 hours. The resulting thick paste of SWNT-neutral water was transferred to a beaker. The beaker was covered with aluminum foil with few holes and allowed to dry very slowly in a hood at room temperature. In about eight weeks, the sample gradually dried isotropically to a cylindrical dense mass referred to as a carbon nanotube block (CNB). There was no carbonaceous matter sticking to any of the glass surface where the SWNT-water paste made contact. The approximate dimensions of the cylindrical SWNT block that weighed 28 g were 1.5 in. (.3.8 cm) in diameter and 1.25 in. (.3.2 cm) in height. The bottom portion of the cylinder that was in contact with the beaker surface was slightly wider, indicating some resistance to shrinking. The cylindrical mass also consisted of several pores. The cylindrical mass was very tough and could not be broken with a small hammer using considerable force. The mass of the solid could be polished over a fine grain emery paper or even a smooth, stainless steel surface indicative of alignment at finer levels.

Smalley, Richard E.↗

Meshfree Methods

Meshfree methods have undergone substantial development and have received much attention in the last two decades. This new family of numerical methods is designed to inherit the main advantages of the finite element method such as compact supports of shape functions and good approximation properties while, at the same time, overcome the main disadvantages of the finite element method caused by the mesh dependence. The meshfree methods share a common feature that no mesh is needed and shape functions are constructed from sets of points, thus eliminating the need for time consuming mesh generation. The most significant advantage of meshfree methods is the flexibility in customizing approximation functions for desired regularity and for capturing essential physics and features of the particular problems of interest. Adaptivity formulation and multiple-scale solution strategies also can be implemented with relative ease. It has become clear that the meshfree methods provide considerable advantages over the conventional finite element methods in solving problems involving moving discontinuities, evolving material interfaces, multiple-scale phenomena, large material distortion and structural deformation, and fracture and damage processes. This Chapter gives an overview of many classes of meshfree methods, with more detailed discussions on Smoothed Particle Hydrodynamics (SPH), the Reproducing Kernel Particle Method (RKPM), Peridynamics (PD), the Material Point Method (MPM), as well as their applications in various challenging engineering problems.2

Chen, Jiun-Shyan↗

Self-Cleaning Highly Porous TiO 2 Coating Designed by Swelling-Assisted Sequential Infiltration Synthesis (SIS) of a Block Copolymer Template

Photocatalytic self-cleaning coatings with a high surface area are important for a wide range of applications, including optical coatings, solar panels, mirrors, etc. Here, we designed a highly porous TiO 2 coating with photoinduced self-cleaning characteristics and very high hydrophilicity. This was achieved using the swelling-assisted sequential infiltration synthesis (SIS) of a block copolymer (BCP) template, which was followed by polymer removal via oxidative thermal annealing. The quartz crystal microbalance (QCM) was employed to optimize the infiltration process by estimating the mass of material infiltrated into the polymer template as a function of the number of SIS cycles. This adopted swelling-assisted SIS approach resulted in a smooth uniform TiO 2 film with an interconnected network of pores. The synthesized film exhibited good crystallinity in the anatase phase. The resulting nanoporous TiO 2 coatings were tested for their functional characteristics. Exposure to UV irradiation for 1 h induced an improvement in the hydrophilicity of coatings with wetting angle reducing to unmeasurable values upon contact with water droplets. Furthermore, their self-cleaning characteristics were tested by measuring the photocatalytic degradation of methylene blue (MB). The synthesized porous TiO 2 nanostructures displayed promising photocatalytic activity, demonstrating the degradation of approximately 92% of MB after 180 min under ultraviolet (UV) light irradiation. Thus, the level of performance was comparable to the photoactivity of commercial anatase TiO 2 nanoparticles of the same quantity. Our results highlight a new robust approach for designing hydrophilic self-cleaning coatings with controlled porosity and composition.

36 MATERIALS SCIENCE↗

Comet Shoemaker-Levy 9: Impact on Jupiter and plume evolution

The impact of fragments of Comet Shoemaker-Levy 9 on Jupiter and the resulting vapor plume expansion are investigated by conducting three-dimensional numerical simulations using the smoothed particle hydrodynamics (SPH) method. An icy body, representing the cometary fragments, with a velocity of 60 km/sec and a diameter of 2 km can penetrate to 350 km below the 1-bar pressure level in the atmosphere. Most of the initial kinetic energy of the fragment is transferred to the atmosphere between 50 km and 300 km below the 1-bar pressure level. The shock-heated atmospheric gas in the wake is totally dissociated and partially ionized. Scaling our SPH results to other sizes indicates that fragments larger than approximately 100 m in diameter can penetrate to below the visible cloud decks. The energy deposited in the atmosphere is explosively released in the upward expansion of the resulting plume. The plume preferentially expands upward rather than horizontally due to the density gradient of the ambient atmosphere. It rises greater than or equal to 10(exp 2) km in approximately 10(exp 2) sec. Eventually the total atmospheric mass ejected to above 1 bar is greater than or equal to 40 times the initial mass of the impactor. The plume temperature at a radius approximately 10(exp 3) km is greater than 10(exp 3) K for 10(exp 3) sec for a 2-km fragment. We predict that impact-induced plumes will be observable with the remote sensing instruments of the Galileo spacecraft. As the impact site rotates into the view of Earth some 20 min after the impact, the plume expansion will be observable using the Hubble Space Telescope (HST) and from visible and infrared instruments on groundbased telescopes. The rising plume reaches approximately 3000 km altitude in approximately 10 min and will be visible from Earth.

Takata, Toshiko↗

A Model for Backscattering from Quasi Periodic Corn Canopies at L-Band

In this study, a model for backscattering at L-band from a corn canopy is proposed. The canopy consists of a quasi-periodic distribution of stalks and a random distribution of leaves. The Distorted Born Approximation (DBA) is employed to calculate the single scattered return from the corn field. The new feature of the method is that the coherence of the stalks in the row direction is incorporated in the model in a systematic fashion. Since the wavelength is on the order of the distance between corn stalks in a row, grating lobe behavior is observed at certain azimuth angles of incidence. The results are compared with experimental values measured in Huntsville, Alabama in 1998. The mean field and the effective dielectric constant of the canopy are obtained by using the Foldy approximation. The stalks are placed in the effective medium in a two dimensional lattice to simulate the row structure of a corn field. In order to mimic a real corn field, a quasi-periodic stalk distribution is assumed where the stalks are given small random perturbations about their lattice locations. Corn leaves are also embedded in the effective medium and the backscattered field from the stalks and the leaves is computed. The backscattering coefficient is calculated and averaged over successive stalk position perturbations. It is assumed that soil erosion has smoothed the soil sufficiently so that it can be assumed flat. Corn field backscatter data was collected from cornfields during the Huntsville 98 experimental campaign held at Alabama A&M University Research Station, Huntsville, Alabama in 1998 using the NASA/GW truck mounted radar. Extensive ground truth data was collected. This included soil moisture measurements and corn plant architectural data to be used in the model. In particular, the distances between the stalks in a single row have been measured. The L-band radar backscatter data was collected for both H and V polarizations and for look angles of 15o and 45o over a two week period under varying soil moisture conditions. These measured backscattering values will be compared with the model backscattering values and a discussion of the results will be presented.

Lang, R.↗

Use of multiple gauges and microwave attenuation of precipitation for satellite verification

In this paper both a microwave attenuation measurement along a horizontal line and multiple point gauge measurements are analyzed as possible ground-truth designs to validate satellite precipitation retrieval algorithms at the field of view spatial level (typically about 20 km). The design consists of comparing a sequence of pairs of contemporaneous measurements taken from the ground and from space. The authors examine theoretically the variance of expected differences between the two systems. The line average measurement leads to a smaller mean-square error compared to the case of a single point gauge, since some of the small-scale variability of the rain field is smoothed away by the line integration. The multiple point gauge measurements also give smaller mean-square error than that of a single point gauge. The centroid of the line and point gauge configurations are considered to be located randomly inside the field of view for different overpasses. A space-time spectral formalism is used with a noise-forced diffusive rain field to find the mean-square error. By considering instantaneous ground and satellite measurement pairs over about 50 visits when raining, we can reduce the expected error to approximately 10% of the standard deviation of climatological variability. This is considered to be a useful level of tolerance for identifying biases in the retrieval algorithms. It is found that the multiple point gauges (especially two gauges) are the economical ground-truth design compared to the microwave attenuation based on the mean-square error comparison. The major finding of this study is that a significant improvement over the point gauge is obtained by adding a single additional piece of information; adding more gauges or extending the line of attenuation is not an important improvement.

Ha, Eunho↗

Development of the endolymphatic sac and duct in the Japanese red-bellied newt, Cynops pyrrhogaster

The development and maturation of the endolymphatic sac (ES) and duct (ED) were studied in the newt Cynops pyrrhogaster. The ES first appears as an oval capsule at the dorsal-medial tip of the otic vesicle at stage 39, about 11 days after oviposition. The ES consists of polymorphous epithelial cells with a minimum of cytoplasm. The intercellular space (IS) between the epithelial cells is narrow and has a smooth surface. At stage 44, the size of the ES increases as many vacuoles in the IS become filled. At stage 46, 18 days after oviposition, the ES elongates markedly and a slit-like lumen is found in the ES. The epithelium contains a few cell organelles which are scattered in the cytoplasm. The vacuoles in the IS are fused, which expands the IS. Two days later (stage 48), floccular material (endolymph) is present in the expanded lumen. The IS dilates and has a wide and irregular appearance. At stage 50, approximately 26 days after oviposition, the ES extends and expands significantly and crystals (otoconia) can now be seen in the widened lumen of the ES. The cytoplasm of the cuboidal epithelial cells contains an abundance of vesicles surrounded by ribosomes and Golgi complexes. Intercellular digitations are formed in the expanded IS. At stage 54, the ES forms a large bellow-like pouch. Numerous otoconia accumulate in the lumen. Free floating cells and cell debris can be seen in the lumen at this stage. The epithelial cells contain numerous cytoplasmic organelles which are evenly distributed in the cytoplasm. Granules are found in the apical and lateral cytoplasm. The IS is loose and displays a labyrinthine appearance. The primitive ED first appears as a connection between the ES and the saccule but no lumen is present inside at stage 39. At stage 46, a narrow lumen is formed in the ED, which corresponds to the formation of the ES lumen. At stage 50, as the ED extends, floccular material is seen in the lumen. At stage 54, the ED bears numerous microvilli on its luminal surface. Otoconia and endolymph are present in the ED. Tight junctions between the epithelial cells are formed at stage 46. A fully developed intercellular junctional complex is produced at stage 54. Based on the development of the ES and ED, the maturation of function of the ES and ED are discussed.

NASA Discipline Neuroscience↗

Modular Flooring System

The modular flooring system (MFS) was developed to provide a portable, modular, durable carpeting solution for NASA fs Robotics Alliance Project fs (RAP) outreach efforts. It was also designed to improve and replace a modular flooring system that was too heavy for safe use and transportation. The MFS was developed for use as the flooring for various robotics competitions that RAP utilizes to meet its mission goals. One of these competitions, the FIRST Robotics Competition (FRC), currently uses two massive rolls of broadloom carpet for the foundation of the arena in which the robots are contained during the competition. The area of the arena is approximately 30 by 72 ft (approximately 9 by 22 m). This carpet is very cumbersome and requires large-capacity vehicles, and handling equipment and personnel to transport and deploy. The broadloom carpet sustains severe abuse from the robots during a regular three-day competition, and as a result, the carpet is not used again for competition. Similarly, broadloom carpets used for trade shows at convention centers around the world are typically discarded after only one use. This innovation provides a green solution to this wasteful practice. Each of the flooring modules in the previous system weighed 44 lb (.20 kg). The improvements in the overall design of the system reduce the weight of each module by approximately 22 lb (.10 kg) (50 %), and utilize an improved "module-to-module" connection method that is superior to the previous system. The MFS comprises 4-by-4-ft (.1.2-by- 1.2-m) carpet module assemblies that utilize commercially available carpet tiles that are bonded to a lightweight substrate. The substrate surface opposite from the carpeted surface has a module-to-module connecting interface that allows for the modules to be connected, one to the other, as the modules are constructed. This connection is hidden underneath the modules, creating a smooth, co-planar flooring surface. The modules are stacked and strapped onto durable, commercially available drywall carts for storage and/or transportation. This method of storage and transportation makes it very convenient and safe when handling large quantities of modules.

Thate, Robert↗

Aerodynamic design of turbomachinery blading in three-dimensional flow - An application to radial inflow turbines

A computational method, based on a theory for turbomachinery blading design in three-dimensional inviscid flow, is applied to a parametric design study of a radial inflow turbine wheel. As the method requires the specification of swirl distribution, a technique for its smooth generation within the blade region is proposed. Excellent agreements have been obtained between the computed results from this design method and those from direct Euler computations, demonstrating the correspondence and consistency between the two. The computed results indicate the sensitivity of the pressure distribution to a lean in the stacking axis and a minor alteration in the hub/shroud profiles. Analysis based on Navier-Stokes solver shows no breakdown of flow within the designed blade passage and agreement with that from design calculation; thus the flow in the designed turbine rotor closely approximates that of an inviscid one. These calculations illustrate the use of a design method coupled to an analysis tool for establishing guidelines and criteria for designing turbomachinery blading.

Yang, Y. L.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth Observing System) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the "Independent Pixel Approximation" or IPA. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM (atmospheric radiation measurement) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth observing system) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of 'radiative smoothing' of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. 'Offbeam' Lidar returns are now being used to estimate the cloud 'radiative Green's function' (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the 'Independent Pixel Approximation' or IPA. This and other measurements of 3D (three dimensional) transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international 'Intercomparison of 3D Radiation Codes' (I3RC) program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

Characterization of a Ultra-high Temperature Ceramic Composite

Ultra-high temperature ceramics (UHTC) are of interest for hypersonic vehicle leading edge applications. Monolithic UHTCs are of concern because of their low fracture toughness and brittle behavior. UHTC composites (UHTCC) are being investigated as a possible approach to overcome these deficiencies. In this study a small sample of a UHTCC was evaluated by limited mechanical property tests, furnace oxidation exposures, and oxidation exposures in a flowing environment. The composite was prepared from a carbon fiber perform using ceramic particulates and a preceramic polymer. The as-received composite plate was non-uniform from front to back surface. Plate dimensions were 150 x 150 x 6 mm. The back surface had a fibrous, uniform appearance; XRD analysis revealed the presence of Sic and C. The front surface was smooth and non-uniform in appearance with evidence of a coarse grain structure produced by a liquid phase; XRD analysis revealed the presence of HfB2. Microcracks were present throughout the thickness as one might expect from a carbon fiber reinforced composite with attendant large thermal expansion mismatch between the matrix phases and the fibers. The HfB2 phase on the front surface was comparable in thickness to a fiber ply or about 0.6 mm, and surface microcracks were evident. Limited four point flexural tests were carried out at span to depth ratios of approximately 14 and 16 with markedly different results. Tests were run with the front or the back surface in tension. At the shorter span to depth failures occurred under a loading pin for both orientations. At a span to depth of 16 failures occurred in the center of the span with fracture clearly initiating from a tensile failure. Ultimate flexural strength, strain at ultimate stress, stress and strain at deviation from linear elastic behavior are reported. Strains at ultimate stress ranged from about 0.6 to 0.7 % for the back surface in tension, and 0.4 to 0.6 for the front surface in tension. At constant span to depth the strain at ultimate stress was about 0.2% greater for the back surface in tension and the ultimate strength was also higher. Strengths were in line with predictions from theory. Furnace oxidation studies were carried out at 1627 and 1927OC in a static furnace environment using ten minute cycles and one, five, and ten cycles. Limited oxidation studies were also carried out in a flowing oxyacetylene torch environment. Specimens were photographed, and weight and dimensional changes were determined. XRD and SEM characterizations were performed. Weight losses were attributed primarily to carbon fiber oxidation. The composite survived the torch test with little visible distress. Further details will be determined once metallographic studies are completed.

Levine, Stanley R.↗

Combined Uncertainty and A-Posteriori Error Bound Estimates for General CFD Calculations: Theory and Software Implementation

This workshop presentation discusses the design and implementation of numerical methods for the quantification of statistical uncertainty, including a-posteriori error bounds, for output quantities computed using CFD methods. Hydrodynamic realizations often contain numerical error arising from finite-dimensional approximation (e.g. numerical methods using grids, basis functions, particles) and statistical uncertainty arising from incomplete information and/or statistical characterization of model parameters and random fields. The first task at hand is to derive formal error bounds for statistics given realizations containing finite-dimensional numerical error [1]. The error in computed output statistics contains contributions from both realization error and the error resulting from the calculation of statistics integrals using a numerical method. A second task is to devise computable a-posteriori error bounds by numerically approximating all terms arising in the error bound estimates. For the same reason that CFD calculations including error bounds but omitting uncertainty modeling are only of limited value, CFD calculations including uncertainty modeling but omitting error bounds are only of limited value. To gain maximum value from CFD calculations, a general software package for uncertainty quantification with quantified error bounds has been developed at NASA. The package provides implementations for a suite of numerical methods used in uncertainty quantification: Dense tensorization basis methods [3] and a subscale recovery variant [1] for non-smooth data, Sparse tensorization methods[2] utilizing node-nested hierarchies, Sampling methods[4] for high-dimensional random variable spaces.

CFD↗

Nonlinear estimation theory applied to the interplanetary orbit determination problem.

Martingale theory and appropriate smoothing properties of Loeve (1953) have been used to develop a modified Gaussian second-order filter. The performance of the filter is evaluated through numerical simulation of a Jupiter flyby mission. The observations used in the simulation are on-board measurements of the angle between Jupiter and a fixed star taken at discrete time intervals. In the numerical study, the influence of each of the second-order terms is evaluated. Five filter algorithms are used in the simulations. Four of the filters are the modified Gaussian second-order filter and three approximations derived by neglecting one or more of the second-order terms in the equations. The fifth filter is the extended Kalman-Bucy filter which is obtained by neglecting all of the second-order terms.

Tapley, B. D.↗

Statistical theories of Langmuir turbulence. II - Subsonic to sonic transition

The subsonic limit of the quadratic direct interaction approximation (DIA) applied to the Zakharov equations is compared with the cubic DIA applied to the nonlinear Schroedinger equation, which is the subsonic limit of the Zakharov equations. Comparisons with Monte Carlo simulations of a truncated system show that the first theory more accurately describes the regime of stationary turbulence, while the second theory more accurately describes the subsonic evolution of the modulational instability. The weak turbulence limits of the two theories describe the sonic and subsonic regimes, respectively. The addition of vertex corrections to the DIA leads to a hybrid weak turbulence theory that smoothly interpolates between the sonic and subsonic regimes.

Dubois, D. F.↗