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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Microwave background fluctuations due to the Sunyaev-Zel'dovich effects in pancakes

We calculate distortions in the microwave background radiation from the Sunyaev-Zel'dovich effect, produced by hot gas in large (approximately 100 Mpc) pancakes. The large-scale distribution of the pancakes is taken to be that of a Voronoi foam. Fluctuations for this scenario are estimated to be on the order of delta T/T is approximately 10(exp -5). Using computer simulations, we produce several 32 deg x 32 deg images with 0.25 deg resolution. These images show characteristic linear features produced when a pancake is viewed nearly edge-on. By calculating the two-point and the degenerate three-point correlation functions, we are able to statistically detect such non-Gaussian features even in the presence of a relatively large amount of Gaussian noise. The degenerate three-point correlation function is found to be particularly useful since it is insensitive to correlated Gaussian noise. We also smooth our data over a 7 deg Full Width at Half Maximum (FWHM) Gaussian window to simulate the Cosmic Background Explorer Satellite (COBE) observations. We find that under such low-resolution conditions, the features are highly suppressed.

Subbarao, M. U.↗
NASA NTRS1994

Hydrodynamics of binary coalescence. 1: Polytropes with stiff equations of state

We performed a series of three-dimensional hydrodynamic calculations of binary coalescence using the smoothed particle hydrodynamics (SPH) method. The initial conditions are exact polytropic equilibrium configurations on the verge of dynamical instability. We consider synchronized equilbria only and concentrate on stiff equations of state, with adiabatic Gamma greater than 5/3. We assume that the polytropic constants (K identically equal to P/(rho(exp Gamma)) are the same for both components. These conditions apply well to models of neutron star binaries. Accordingly, we discuss our results in the context of the Laser Interferometer Gravitational-Wave Observatory (LIGO) project, and we calculate the emission of gravitational radiation in the quadruple approximation. The fully nonlinear development of the instability is followed using SPH until a new equilibrium configuration is reached by the system. We find that the properties of this final configuration depend sensitively on both the compressibility and mass ratio. An axisymmetric merged configuration is always produced when the adiabatic exponent Gamma approximately less than 2.3. As a consequence, the emission of gravitational radiation shuts off abruptly right after the onset of dynamical instability. In contrast, triaxial merged configurations are obtained when Gamma approximately greater than 2.3, and the system continues to emit gravitational waves after the final coalescence. Systems with mass ratios q not equal to 1 typically become dynamically unstable before the onset of mass transfer. Stable mass transfer from one neutron star to another in a close binary is therefore probably ruled out. For a mass ratio q approximately less than 0.5, however, dynamical mass transfer can temporarily retard the coalescence by causing a rapid reexpansion of the binary into a new, slightly eccentric but dynamically stable orbit. The maximum amplitude h(sub max) and peak luminosity L(sub max) of the gravitational waves emitted during the final coalescence are nearly independent of Gamma, but depend sensitively on the mass ratio q. The approximate scalings we find are h(sub max) varies as q(exp 2) and L(sub max) varies as q(exp 6) for q close to unity.

Rasio, Frederic A.↗
DOE OSTI2022

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

King Crater and its environs, part M

Superb photographs of King Crater were returned by the Apollo 16 crew. The metric and panoramic cameras provided stereoscopic coverage with which high quality topographic maps may be made. King crater, approximately 75 km in diameter and, on the average, 3.8 km deep, is among the most significant features photographed and visually studied from lunar orbit on the Apollo 16 mission. It is the freshest crater on the far side in its size range and has several unique characteristics. (1) Although its ejecta blanket displays very sharp features, the crater does not appear to have extensive ray systems and secondary crater chains. (2) As depicted for the first time by Apollo 16 photographs, the ejecta blanket is strewn with numerous flow fronts and scarps. (3) A flat and smooth, dark, lavalike deposit 15 km across is on the northern rim of the crater. (4) The central peak of the crater is a Y-shaped structure that begins at the southern rim. (5) Bright units that are aligned with the central peak on the northern wall and rim suggest that a preexisting tabular body was excavated by the cratering process.

El-Baz, F.↗
NASA NTRS1994

Dynamics of massive black holes as a possible candidate of Galactic dark matter

If the dark halo of the Galaxy is comprised of massive black holes (MBHs), then those within approximately 1 kpc will spiral to the center, where they will interact with one another, forming binaries which contract, owing to further dynamical friction, and then possibly merge to become more massive objects by emission of gravitational radiation. If successive mergers would invariably lead, as has been proposed by various authors, to the formation of a very massive nucleus of 10(exp 8) solar mass, then the idea of MBHs as a dark matter candidate could be excluded on observational grounds, since the observed limit (or value) for a Galactic central black hole is approximately 10(exp 6.5) solar mass. But, if successive mergers are delayed or prevented by other processes, such as the gravitational slingshot or rocket effect of gravitational radiation, then a large mass accumulation will not occur. In order to resolve this issue, we perform detailed N-body simulations using a modfied Aarseth code to explore the dynamical behavior of the MBHs, and we find that for a 'best estimate' model of the Galaxy a runaway does not occur. The code treates the MBHs as subject to the primary gravitational forces of one another and to the smooth stellar distribution, as well as the secondary perturbations in their orbits due to another and to the smooth stellar distribution, as well as the secondary perturbations in their orbits due to dynamical friction and gravitational radiation. Instead of a runaway, three-body interactions between hard binaries and single MBHs eject massive objects before accumulation of more than a few units, so that typically the center will contain zero, one, or two MBHs. We study how the situation depends in detail on the mass per MBH, the rotation of the halo, the mass distribution within the Galaxy, and other parameters. A runaway will most sensitively depend on the ratio of initial (spheroid/halo) central mass densities and secondarily on the typical values for the mass per MBH, with the rough dividing line, using Galactic parameters, being M(sub BH) less than or = 10(exp 6.5) solar mass. Using parameters from Lacey & Ostriker (1985) and our most accurate model for Galaxy, no runaway occurs.

Xu, Guohong↗
NASA NTRS2011

Machine Vision for Relative Spacecraft Navigation During Approach to Docking

This paper describes a machine vision system for relative spacecraft navigation during the terminal phase of approach to docking that: 1) matches high contrast image features of the target vehicle, as seen by a camera that is bore-sighted to the docking adapter on the chase vehicle, to the corresponding features in a 3d model of the docking adapter on the target vehicle and 2) is robust to on-orbit lighting. An implementation is provided for the case of the Space Shuttle Orbiter docking to the International Space Station (ISS) with quantitative test results using a full scale, medium fidelity mock-up of the ISS docking adapter mounted on a 6-DOF motion platform at the NASA Marshall Spaceflight Center Flight Robotics Laboratory and qualitative test results using recorded video from the Orbiter Docking System Camera (ODSC) during multiple orbiter to ISS docking missions. The Natural Feature Image Registration (NFIR) system consists of two modules: 1) Tracking which tracks the target object from image to image and estimates the position and orientation (pose) of the docking camera relative to the target object and 2) Acquisition which recognizes the target object if it is in the docking camera Field-of-View and provides an approximate pose that is used to initialize tracking. Detected image edges are matched to the 3d model edges whose predicted location, based on the pose estimate and its first time derivative from the previous frame, is closest to the detected edge1 . Mismatches are eliminated using a rigid motion constraint. The remaining 2d image to 3d model matches are used to make a least squares estimate of the change in relative pose from the previous image to the current image. The changes in position and in attitude are used as data for two Kalman filters whose outputs are smoothed estimate of position and velocity plus attitude and attitude rate that are then used to predict the location of the 3d model features in the next image.

Chien, Chiun-Hong↗
NASA NTRS1992

Effect of Index of Refraction on Radiation Characteristics in a Heated Absorbing, Emitting, and Scattering Layer

The index of refraction can considerably influence the temperature distribution and radiative heat flow in semitransparent materials such as some ceramics. For external radiant heating, the refractive index influences the amount of energy transmitted into the interior of the material. Emission within a material depends on the square of its refractive index, and hence this emission can be many times that for a biackbody radiating into a vacuum. Since radiation exiting through an interface into a vacuum cannot exceed that of a blackbody, there is extensive reflection at the internal surface of an interface, mostly by total internal reflection. This redistributes energy within the layer and tends to make its temperature distribution more uniform. The purpose of the present analysis is to show that, for radiative equilibrium in a gray layer with diffuse interfaces, the temperature distribution and radiative heat flux for any index of refraction can be obtained very simply from the results for an index of refraction of unity. For the situation studied here, the layer is subjected to external radiative heating incident on each of its surfaces. The material emits, absorbs, and isotropically scatters radiation. For simplicity the index of refraction is unity in the medium surrounding the layer. The surfaces of the layer are assumed diffuse. This is probably a reasonable approximation for a ceramic layer that has not been polished. When transmitted radiation or radiation emitted from the interior reaches the inner surface of an interface, the radiation is diffused and some of it thereby placed into angular directions for which there is total internal reflection. This provides a trapping effect for retaining energy within the layer and tends to equalize its temperature distribution. An analysis of temperature distributions in absorbing-emitting layers, including index of refraction effects, was developed by Gardon (1958) to predict cooling and heat treating of glass plates. The interfaces were optically smooth; the resulting specular reflections were computed from the Fresnel reflection laws. This provides a somewhat different behavior than for diffuse interfaces. A similar application was for heating that occurs in a window of a re-entry vehicle (Fowle et al., 1969). A number of recent papers (Rokhsaz and Dougherty, 1989; Ping and Lallemand, 1989; Crosbie and Shieh, 1990) further examined the effects of Fresnel boundary reflections and nonunity refractive index. Other examples of analyses of both steady and transient heat transfer to single or multiple plane layers (Amlin and Korpela, 1979; Tarshis et al., 1969) have used diffuse assumptions at the interfaces as in the present study

Siegel, R.↗
DOE OSTI2022

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Influence of aerodynamic roughness length on aeolian processes: Earth, Mars, Venus

The aerodynamic roughness length (z sub 0) is the height at which a wind profile assumes a zero velocity. The lower part of the atmospheric boundary layer will be impeded by friction with the surface. An increase in surface roughness will also increase the shear stress required to initiate particle entrainment by the wind. Bagnold (1941) estimated z sub 0 as being 1/30 of the mean particle size. In Nature, surface roughness is composed of nonerodible elements as well as sand-size erodible particles. To assess z sub 0 values as a function of terrain, field experiments were conducted to obtain wind profiles monitored over natural surfaces at 15 sites in the Mojave Desert, Death Valley, and Nye County, Nevada. These sites span a variety of arid-land terrains, including smooth playas, alluvial fans, and lava flows; z sub 0 values ranged from 0.0001 cm to 1 cm. These values were incorporated in a threshold model and a flux model to assess transport efficiency over such terrains in three planetary environments (Venus, Earth, and Mars), and for particle sizes ranging from 60-500 micron. Threshold and flux are a function of planetary environment, particle density and size (Dp), and z sub 0, and the shear velocity of 1.2 x U sub *t (for Dp = 250 micron and z sub 0 = 0.84). Results show that flux on Mars is approximately 14 g/(cm x s), on Earth it is approximately 3 g/(cm x s), and on Venus 0.5 g/(cm x s). Under all planetary environments, the results also show a dramatic decrease in the flux for particles greater than 200 microns when z sub 0 increases above 0.0085 cm (corresponding to sites consisting of alluvium). When z sub 0 approaches 0.03 cm (corresponding to a mantled pahoehoe lava), the flux diminishes.

Blumberg, Dan G.↗
NASA NTRS2003

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗
NASA NTRS1938

A preliminary investigation of boundary-layer transition along a flat plate with adverse pressure gradient

Boundary-layer surveys were made throughout the transition region along a smooth flat plate placed in an airstream of practically zero turbulence and with an adverse pressure gradient. The boundary-layer Reynolds number at the laminar separation point was varied from 1,800 to 2,600. The test data, when considered in the light of certain theoretical deductions, indicated that transition probably began with separation of the laminar boundary layer. The extent of the transition region, defined as the distance from a calculated laminar separation point to the position of the first fully developed turbulent boundary-layer profile, could be expressed as a constant Reynolds number run of approximately 70,000. Some speculations are presented concerning the application of the foregoing concepts, after certain assumptions have been made, to the problem of the connection between transition on the upper surface of an airfoil at high angles of attack and the maximum lift.

Von Doenhoff, Albert E↗
NASA NTRS1993

An early solar dynamo prediction: Cycle 23 is approximately cycle 22

In this paper, we briefly review the 'dynamo' and 'geomagnetic precursor' methods of long-term solar activity forecasting. These methods depend upon the most basic aspect of dynamo theory to predict future activity, future magnetic field arises directly from the magnification of pre-existing magnetic field. We then generalize the dynamo technique, allowing the method to be used at any phase of the solar cycle, through the development of the 'Solar Dynamo Amplitude' (SODA) index. This index is sensitive to the magnetic flux trapped within the Sun's convection zone but insensitive to the phase of the solar cycle. Since magnetic fields inside the Sun can become buoyant, one may think of the acronym SODA as describing the amount of buoyant flux. Using the present value of the SODA index, we estimate that the next cycle's smoothed peak activity will be about 210 +/- 30 solar flux units for the 10.7 cm radio flux and a sunspot number of 170 +/- 25. This suggests that solar cycle #23 will be large, comparable to cycle #22. The estimated peak is expected to occur near 1999.7 +/- 1 year. Since the current approach is novel (using data prior to solar minimum), these estimates may improve when the upcoming solar minimum is reached.

Schatten, Kenneth H.↗
DOE OSTI2023

Framework for phase transitions between the Maxwell and Gibbs constructions

By taking the nucleon-to-quark phase transition within a neutron star as an example, we present a thermodynamically consistent method to calculate the equation of state of ambient matter so that transitions that are intermediate to those of the familiar Maxwell and Gibbs constructions can be described. This method does not address the poorly known surface tension between the two phases microscopically (as, for example, in the calculation of the core pasta phases via the Wigner-Seitz approximation) but instead combines the local and global charge neutrality conditions characteristic of the Maxwell and Gibbs constructions, respectively. Overall charge neutrality is achieved by dividing the leptons to those that obey local charge neutrality (Maxwell) and those that maintain global charge neutrality (Gibbs). The equation of state is obtained by using equilibrium constraints derived from minimizing the total energy density. The results of this minimization are then used to calculate neutron star mass-radius curves, tidal deformabilities, equilibrium and adiabatic sound speeds, and nonradial g-mode oscillation frequencies for several intermediate constructions. Various quantities of interest transform smoothly from their Gibbs structures to those of Maxwell as the local-to-total electron ratio η, introduced to mimic the hadron-to-quark interface tension from 0 (Gibbs) to ∞ (Maxwell), is raised from 0 to 1. As a result, a notable exception is the g-mode frequency for the specific case of η = 1 for which a gap appears between the quark and hadronic branches.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Ultraviolet photometry of the eclipsing variable CW Cephei

A series of photometric observations was made of the eclipsing variable CW Cephei on the OAO 2. Approximate elements were derived from the eclipse depths and shape of the secondary. Persistent asymmetries and anomalous light variations, larger than the expected experimental error, were also found, subsequent ground-based observations show H alpha entirely in emission, indicating the presence of an extended gaseous system surrounding one or both components. A detailed comparison was made of the flux distribution of the binary relative to that for the nominally unreddened stars delta Pic, BlIII, and eta Aur, B3V, to investigate the effects of interstellar extinction. The resultant extinction curves, normalized at a wavelength of 3330 A, show a relatively smooth increase with decreasing wavelength.

Sobieski, S.↗
NASA NTRS1994

Comments on the statistical analysis of excess variance in the COBE differential microwave radiometer maps

Cosmic anisotrophy produces an excess variance sq sigma(sub sky) in the Delta maps produced by the Differential Microwave Radiometer (DMR) on cosmic background explorer (COBE) that is over and above the instrument noise. After smoothing to an effective resolution of 10 deg, this excess sigma(sub sky)(10 deg), provides an estimate for the amplitude of the primordial density perturbation power spectrum with a cosmic uncertainty of only 12%. We employ detailed Monte Carlo techniques to express the amplitude derived from this statistic in terms of the universal root mean square (rms) quadrupole amplitude, (Q sq/RMS)(exp 0.5). The effects of monopole and dipole subtraction and the non-Gaussian shape of the DMR beam cause the derived (Q sq/RMS)(exp 0.5) to be 5%-10% larger than would be derived using simplified analytic approximations. We also investigate the properties of two other map statistics: the actual quadrupole and the Boughn-Cottingham statistic. Both the sigma(sub sky)(10 deg) statistic and the Boughn-Cottingham statistic are consistent with the (Q sq/RMS)(exp 0.5) = 17 +/- 5 micro K reported by Smoot et al. (1992) and Wright et al. (1992).

Wright, E. L.↗
DOE OSTI2025

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗
DOE OSTI2023

Thickness dependence of piezo-bimorph adaptive mirror bending

A new generation of adaptive x-ray optics (AXO) is being installed on high-coherent-flux x-ray beamlines worldwide to correct and control the optical wavefront with sub-nm precision. These ultra-smooth mirrors achieve high reflectivities at glancing angles of incidence and can be hundreds of mm long. One type of adaptive x-ray mirror relies on piezoelectric ceramic strips which are segmented into channels and actuated to induce local, longitudinal bending, generating one-dimensional shape changes in the mirror substrate. A recently described mirror model uses a three-layer geometry with parallel actuators on the front and back surfaces of a thicker mirror substrate. By analogy to a solved problem in the thermal actuation of a tri-metal strip, we show that the achievable bending radius varies approximately as the square of the substrate thickness. We provide an analytic solution and simulate bending using a finite-element model.

47 OTHER INSTRUMENTATION↗
NASA NTRS1975

Vertical temperature and density patterns in the Arctic mesosphere analyzed as gravity waves

Rocket soundings conducted from high latitude sites in the Arctic mesosphere are described. Temperature and wind profiles and one density profile were observed independently to obtain the thermodynamic structure, the wind structure, and their interdependence in the mesosphere. Temperature profiles from all soundings were averaged, and a smooth curve (or series of smooth curves) drawn through the points. A hydrostatic atmosphere based on the average, measured temperature profile was computed, and deviations from the mean atmosphere were analyzed in terms of gravity wave theory. The vertical wavelengths of the deviations were 10-20 km, and the wave amplitudes slowly increased with height. The experimental data were matched by calculated gravity waves having a period of 15-20 minutes and a horizontal wavelength of 60-80 km. The wind measurements are consistent with the thermodynamic measurements. The results also suggest that gravity waves travel from East to West with a horizontal phase velocity of approximately 60 m sec-1.

Eberstein, I. J.↗