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 577 records · Page 32

A Gravity Analysis of the Subsurface Structure of the Utopia Impact Basin

The large, shallow, circular depression in Utopia Planitia has been identified as a huge impact basin, based on both geological evidence and detailed analysis of MOLA topography. Its diameter (approximately 3000 km) is equivalent to that of the Hellas basin, as is its inferred age (early Noachian). However, there the similarity ends. Their appearance, both surficially and geophysically, are virtually polar opposites. Whereas Hellas is extremely deep with rough terrain and large slopes, high-precision MOLA measurements were required to unambiguously define the smooth, shallow, almost imperceptible bowl of the Utopia basin. Conversely, Utopia displays one of the largest (non-Tharsis-related) positive geoid anomalies on Mars, in contrast to a more subdued negative anomaly over Hellas. As these two features presumably formed roughly contemporaneously by similar mechanisms, it is reasonable to assume that they were originally quite similar, and that their differences are due largely to different paths of subsequent modification. The obvious source for these differences is in their elevations: Hellas is located in the southern highlands at a rim elevation of about 3km, whereas Utopia is in the lowlying northern plains, at an average elevation of 4 km. Thus Utopia has been in an especially gravitationally favorable position to be subjected to infilling, for example, by lava flows, sedimentation, or water. In fact, its floor was almost certainly the lowest point on the planet at one time, and it would have been the termination point for down-slope drainage from over two-thirds of Mars. Thus the nature of the material filling this basin has strong connections to the sedimentary and/or volcanic processes acting on Mars in the Noachian and Early Hesperian periods. In particular, it may be able to shed some light on amount and persistence of water on early Mars in general and in the Utopia basin in particular. In this study I will use the inferred early correspondence between Hellas and Utopia to investigate Utopia's subsurface structure.

Barnerdt, W. B.↗

ARISE (Antarctic Remote Ice Sensing Experiment) in the East 2003: Validation of Satellite-derived Sea-ice Data Product

Preliminary results are presented from the first validation of geophysical data products (ice concentration, snow thickness on sea ice (h(sub s) and ice temperature (T(sub i))fr om the NASA EOS Aqua AMSR-E sensor, in East Antarctica (in September-October 2003). The challenge of collecting sufficient measurements with which to validate the coarse-resolution AMSR-E data products adequately was addressed by means of a hierarchical approach, using detailed in situ measurements, digital aerial photography and other satellite data. Initial results from a circumnavigation of the experimental site indicate that, at least under cold conditions with a dry snow cover, there is a reasonably close agreement between satellite- and aerial-photo-derived ice concentrations, i.e. 97.2+/-.6% for NT2 and 96.5+/-2.5% for BBA algorithms vs 94.3% for the aerial photos. In general, the AMSR-E concentration represents a slight overestimate of the actual concentration, with the largest discrepancies occurring in regions containing a relatively high proportion of thin ice. The AMSR-E concentrations from the NT2 and BBA algorithms are similar on average, although differences of up to 5% occur in places, again related to thin-ice distribution. The AMSR-E ice temperature (T(sub i)) product agrees with coincident surface measurements to approximately 0.5 C in the limited dataset analyzed. Regarding snow thickness, the AMSR h(sub s) retrieval is a significant underestimate compared to in situ measurements weighted by the percentage of thin ice (and open water) present. For the case study analyzed, the underestimate was 46% for the overall average, but 23% compared to smooth-ice measurements. The spatial distribution of the AMSR-E h(sub s) product follows an expected and consistent spatial pattern, suggesting that the observed difference may be an offset (at least under freezing conditions). Areas of discrepancy are identified, and the need for future work using the more extensive dataset is highlighted.

Massom, Robert A.↗

Enabling Long-range Exploration in Minimization of Multimodal Functions

We consider the problem of minimizing multi-modal loss functions with a large number of local optima. Since the local gradient points to the direction of the steepest slope in an infinitesimal neighborhood, an optimizer guided by the local gradient is often trapped in a local minimum. To address this issue, we develop a novel nonlocal gradient to skip small local minima by capturing major structures of the loss’s landscape in black-box optimization. The nonlocal gradient is defined by a directional Gaussian smoothing (DGS) approach. The key idea of DGS is to conducts 1D long-range exploration with a large smoothing radius along d orthogonal directions in Rd, each of which defines a nonlocal directional derivative as a 1D integral. Such long-range exploration enables the nonlocal gradient to skip small local minima. The d directional derivatives are then assembled to form the nonlocal gradient. We use the Gauss-Hermite quadrature rule to approximate the d 1D integrals to obtain an accurate estimator. The superior performance of our method is demonstrated in three sets of examples, including benchmark functions for global optimization, and two real-world scientific problems.

Zhang, Jiaxin↗

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface Merging Technique to Design GA-Optimized Heat Sinks

High-power density inverters require efficient and small heat dissipation systems. For liquid cooling systems, an effective heat sink design is important to enable higher heat transfer and keep the pressure drop within reasonable limits. A heat sink geometry generated from an extruded, 1D fast Fourier transform can achieve more even temperature distribution and improved heat transfer than a conventional heat sink. Building on this concept, this paper proposes a new heat sink geometry creation algorithm. Instead of one plane profile being extruded, multiple plane profiles, each created from the fast Fourier transform method, are subsequently merged by smooth surfaces along the flow direction of the heat sink. Thus, a 3D heat sink geometry is created that enables 3D coolant flow and efficient heat transfer. The volume of this design was reduced by 50% compared with pin fin heat sinks. Compared with extruded 1D fast Fourier transform heat sinks, the proposed design showed an approximately 5% device temperature reduction and a 20% pressure drop reduction.

Barua, Himel↗

Large $$N_c$$ QCD phase diagram at $$\mu _B=0$$

Abstract Lattice studies suggest that at zero baryon chemical potential and increasing temperature there are three characteristic regimes in QCD that are connected by smooth analytical crossovers: a hadron gas regime at$$T < T_{ch}\sim 155$$ T < T ch ∼ 155 MeV, an intermediate regime, called stringy fluid, at$$T_{ch}< T < \sim 3 T_{ch}$$ T ch < T < ∼ 3 T ch , and a quark-gluon plasma regime at higher temperatures. These regimes have been interpreted to reflect different approximate symmetries and effective degrees of freedom. In the hadron gas the effective degrees of freedom are hadrons and the approximate chiral symmetry of QCD is spontaneously broken. The intermediate regime has been interpreted as lacking spontaneous chiral symmetry breaking along with the emergence of new approximate symmetry, chiral spin symmetry, that is not a symmetry of the Dirac Lagrangian, but is a symmetry of the confining part of the QCD Lagrangian. While the high temperature regime is the usual quark-gluon plasma which is often considered to reflect “deconfinement” in some way. This paper explores the behavior of these regimes of QCD as the number of colors in the theory,$$N_c$$ N c , gets large. In the large$$N_c$$ N c limit the theory is center-symmetric and notions of confinement and deconfinement are unambiguous. The energy density is$$\mathcal{O}(N_c^0)$$ O ( N c 0 ) in the meson gas,$${{\mathcal {O}}}(N_c^1)$$ O ( N c 1 ) in the intermediate regime and$${{\mathcal {O}}}(N_c^2)$$ O ( N c 2 ) in the quark-gluon plasma regime. In the large$$N_c$$ N c limit these regimes may become distinct phases separated by first order phase transitions. The intermediate phase has the peculiar feature that glueballs should exist and have properties that are unchanged from what is seen in the vacuum (up to$$1/N_c $$ 1 / N c corrections), while the ordinary dilute gas of mesons with broken chiral symmetry disappears and approximate chiral spin symmetry should emerge.

Physics↗

Space Flight Handbooks. Volume 3- Planetary Flight Handbook: Supplementary Trajectory Data: Venus to Earth and Mars to Earth - Part 2

Parts 2 and 3 present tabulations of trajectory data for scheduling flights to and from Venus and Mars during the period 1960-2000. Part 2 contains information for outbound flights to these planets; Part 3 contains information for trajectories returning from the planets to Earth. Each Part contains data for single-plane transfers, as well as for broken-plane transfers which employ a midcourse plane-change to eliminate the high speed "ridges." The mathematical analyses employed for all calculations are described in Part 1 of this handbook. To facilitate the construction of round-trip trajectories, the date at the target planet is held fixed while the trip duration is varied in 10-day increments from zero days to the length of that planet's synodic period with Earth. Dates of arrival at the target planet are presented in the extreme right-hand column of Part 2, and dates of departure from the target planet are presented in the extreme left-hand column of Part 3. Thus, by holding Part 2 directly to the left of Part 3, the analyst may easily and rapidly scan all trip possibilities which involve any desired stopover time at the target planet. Within approximately 200 days of each conjunction or opposition, data are presented in 10-day increments at the target planet. Only those trips are listed for which the hyperbolic excess speeds at either or both ends of the trajectory do not exceed 0.6 EMOS (Earth Mean Orbital Speed). For the remaining mission regions, the requirements are so smoothly varying that a 50-day interval in dates at the target planet may be employed; the 10-day interval in trip times is, however, preserved here. In these regions, only those trips are listed for which either or both speeds do not exceed 0.3 EMOS.

INTERPLANETARY FLIGHT↗

Design and evaluation of a 3 million DN series-hybrid thrust bearing

The bearing, consisting of a 150-mm ball bearing and a centrifugally actuated, conical, fluid-film bearing, was fatigue tested. Test conditions were representative of a mainshaft ball bearing in a gas turbine engine operating at maximum thrust load to simulate aircraft takeoff conditions. Tests were conducted up to 16000 rpm and at this speed an axial load of 15568 newtons (3500 lb) was safely supported by the hybrid bearing system. Through the series-hybrid bearing principle, the effective ball bearing speed was reduced to approximately one-half of the shaft speed. It was concluded that a speed reduction of this magnitude results in a ten-fold increase in the ball bearing fatigue life. A successful evaluation of fluid-film bearing lubricant supply failure was performed repeatedly at an operating speed of 10,000 rpm. A complete and smooth changeover to full-scale ball bearing operation was effected when the oil supply to the fluid-film bearing was cut off. Reactivation of the fluid-film oil supply system resulted in a flawless return to the original mode of hybrid operation.

Scibbe, H. W.↗

Some 5-13 micrometer airborne observations of Comet Wilson 1986l: Preliminary results

Comet Wilson was observed from the Kuiper Airborne Observatory approximately 23.6 and 25.7 Apr. 1987, UT (approx. 3 to 5 days after perihelion) using the NASA-Ames Faint Object Grating Spectrometer. Spectrophotometric data were observed with a 21 inch aperture between 5 and 13 micrometer and with a spectral resolution of 50 to 100. Spectra of the inner coma and nucleus reveal a fairly smooth continuum with little evidence of silicate emission. The 5 to 8 micrometer color temperature of the comet was 300 + or - 15 K, approx. 15 percent higher than the equilibrium blackbody temperature. All three spectra of the nucleus show a new emission feature at approx. 12.25 micrometer approx. two channels (.22 micrometer) wide. Visual and photographic observations made during the time of these observations showed a broad faint, possible two component tail. No outburst activity was observed.

Lynch, David K.↗

5- to 13-micron airborne observations of Comet Wilson 1986l

Comet Wilson was observed from the Kuiper Airborne Observatory approximately 23.6 and 25.7 Apr. 1987, UT (approx. 3 to 5 days after perihelion) using the NASA-Ames Faint Object Grating Spectrometer. Spectrophotometric data were observed with a 21 inch aperture between 5 and 13 micrometer and with a spectral resolution of 50 to 100. Spectra of the inner coma and nucleus reveal a fairly smooth continuum with little evidence of silicate emission. The 5 to 8 micrometer color temperature of the comet was 300 + or - 15 K, approx. 15 percent higher than the equilibrium blackbody temperature. All three spectra of the nucleus show a new emission feature at approx. 12.25 micrometer approx. two channels (.22 micrometer) wide. Visual and photographic observations made during the time of these observations showed a broad faint, possible two component tail. No outburst activity was observed.

Lynch, D. K.↗

SEDS1 mission software verification using a signal simulator

The first flight of the Small Expendable Deployer System (SEDS1) is schedule to fly as the secondary payload of a Delta 2 in March, 1993. The objective of the SEDS1 mission is to collect data to validate the concept of tethered satellite systems and to verify computer simulations used to predict their behavior. SEDS1 will deploy a 50 lb. instrumented satellite as an end mass using a 20 km tether. Langley Research Center is providing the end mass instrumentation, while the Marshall Space Flight Center is designing and building the deployer. The objective of the experiment is to test the SEDS design concept by demonstrating that the system will satisfactorily deploy the full 20 km tether without stopping prematurely, come to a smooth stop on the application of a brake, and cut the tether at the proper time after it swings to the local vertical. Also, SEDS1 will collect data which will be used to test the accuracy of tether dynamics models used to stimulate this type of deployment. The experiment will last about 1.5 hours and complete approximately 1.5 orbits. Radar tracking of the Delta II and end mass is planned. In addition, the SEDS1 on-board computer will continuously record, store, and transmit mission data over the Delta II S-band telemetry system. The Data System will count tether windings as the tether unwinds, log the times of each turn and other mission events, monitor tether tension, and record the temperature of system components. A summary of the measurements taken during the SEDS1 are shown. The Data System will also control the tether brake and cutter mechanisms. Preliminary versions of two major sections of the flight software, the data telemetry modules and the data collection modules, were developed and tested under the 1990 NASA/ASEE Summer Faculty Fellowship Program. To facilitate the debugging of these software modules, a prototype SEDS Data System was programmed to simulate turn count signals. During the 1991 summer program, the concept of simulating signals produced by the SEDS electronics systems and circuits was expanded and more precisely defined. During the 1992 summer program, the SEDS signal simulator was programmed to test the requirements of the SEDS Mission software, and this simulator will be used in the formal verification of the SEDS Mission Software. The formal test procedures specification was written which incorporates the use of the signal simulator to test the SEDS Mission Software and which incorporates procedures for testing the other major component of the SEDS software, the Monitor Software.

Pierson, William E.↗

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes - Part 1

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling.

Diskin, Boris↗

Space Flight Handbooks. Volume 3- Planetary Flight Handbook: Supplementary Trajectory Data: Venus to Earth and Mars to Earth - Part 3

Parts 2 and 3 present tabulations of trajectory data for scheduling flights to and from Venus and Mars during the period 1960-2000. Part 2 contains information for outbound flights to these planets; Part 3 contains information for trajectories returning from the planets to Earth. Each Part contains data for single-plane transfers, as well as for broken-plane transfers which employ a midcourse plane-change to eliminate the high speed "ridges." The mathematical analyses employed for all calculations are described in Part 1 of this handbook. To facilitate the construction of round-trip trajectories, the date at the target planet is held fixed while the trip duration is varied in 10-day increments from zero days to the length of that planet's synodic period with Earth. Dates of arrival at the target planet are presented in the extreme right-hand column of Part 2, and dates of departure from the target planet are presented in the extreme left-hand column of Part 3. Thus, by holding Part 2 directly to the left of Part 3, the analyst may easily and rapidly scan all trip possibilities which involve any desired stopover time at the target planet. Within approximately 200 days of each conjunction or opposition, data are presented in 10-day increments at the target planet. Only those trips are listed for which the hyperbolic excess speeds at either or both ends of the trajectory do not exceed 0.6 EMOS (Earth Mean Orbital Speed). For the remaining mission regions, the requirements are so smoothly varying that a 50-day interval in dates at the target planet may be employed; the 10-day interval in trip times is, however, preserved here. In these regions, only those trips are listed for which either or both speeds do not exceed 0.3 EMOS.

Source record↗

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

A discontinuous Galerkin method for general relativistic hydrodynamics in thornado

Discontinuous Galerkin (DG) methods provide a means to obtain high-order accurate solutions in regions of smooth fluid flow while still resolving strong shocks. These and other properties make DG methods attractive for solving problems involving hydrodynamics; e.g., the core-collapse supernova problem. With that in mind we are developing a DG solver for the general relativistic, ideal hydrodynamics equations under a 3+1 decomposition of spacetime, assuming a conformally-flat approximation to general relativity. With the aid of limiters we verify the accuracy and robustness of our code with several difficult test-problems: a special relativistic Kelvin-Helmholtz instability problem, a two-dimensional special relativistic Riemann problem, and a one- and two-dimensional general relativistic standing accretion shock (SAS) problem. We find good agreement with published results, where available. We also establish sufficient resolution for the 1D SAS problem and find encouraging results regarding the standing accretion shock instability (SASI) in 2D.

Dunham, Samuel↗

The analysis of Gemini 7 EEG data by three different computer oriented methods

The problem in assessment of electroencephalograms has been the subject element in visual evaluations of such records as well as the need for extensive training (2 or 3 years) before suitable candidate is capable to perform this task. There have been variety of computer type analyses of EEG during the past ten years or more, and we are still searching for a reliable, inexpensive, online data processing procedure for this task. We are aware our Soviet colleagues are routinely using computer assessment of EEGs, using relatively inexpensive hardware. In our report we endeavor to show how four different assessment procedures (visual assessment, zero crossing technique, smoothing and peak counting technique, and Weibull Statistic) may be utilized, and indicate the advantages and disadvantages peculiar to each. We wish to point out that in the Gemini 7 EEG analog records supplied to the Henry Ford Hospital, Channel II (mid-line central to mid-line occipital electrodes) could not be assessed reliably after approximately 26 hours of flight. Channel I (left central to left occipital electrodes) provided the more assessable EEG up to approximately 54 hours of flight. These channels were compared during the first 26 hours of flight, and no significant difference was found by our various assessment procedures. We, therefore, have limited this report to the findings on the analyses of Channel I during its approximate 54 hours of flight.

Proctor, L. D.↗

Planetary Waves

The mean temperature field at various pressure levels for January and July for both hemispheres are given. In summer (July in the Northern Hemisphere, January in the Southern Hemisphere) the fields are nearly zonally symmetric, but in winter large longitudinal variations are evident. These are mainly of low wave number, i.e., they may be represented by Fourier analysis around the globe using just a few (one or two) waves. This is the basis for representing the climatology of longitudinal variation in terms of wave components, since the fields may be defined with fewer values than by using a grid in longitude that has a sufficiently small interval to adequately represent the smooth variations. The amplitude and phase of temperature and geopotential height for wave numbers one and two, with In(pressure) as the vertical coordinate are given. The fields were calculated and plotted at pressure intervals of 0.2 in In(pressure) and at latitude intervals of 4 deg, but were interpolated to intervals of 0.5 in In(pressure) (approximately 3.5 km) and 10 deg latitude for tabulation. Tables and figures give wave coefficients for the monthly mean temperature fields. Thus they represent the quasi-statinary planetary waves.

Barnett, J. J.↗

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗