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 217 records · Page 12

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗

Stormtime ring current and radiation belt ion transport: Simulations and interpretations

We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Lyons, Larry R.↗

Electron Scattering by Low-frequency Whistler Waves at Earth’s Bow Shock

Electrons are accelerated to nonthermal energies at shocks in space and astrophysical environments. While shock drift acceleration (SDA) has been considered a key process of electron acceleration at Earth’s bow shock, it has also been recognized that SDA needs to be combined with an additional stochastic process to explain the observed power-law energy spectra. Here, we show mildly energetic (∼0.5 keV) electrons are locally scattered (and accelerated while being confined) by magnetosonic-whistler waves within the shock transition layer, especially when the shock angle is large (θ(sub Bn) approximately equal or greater than 70°). When measured by the Magnetospheric Multiscale mission at a high cadence, ∼0.5 keV electron flux increased exponentially in the shock transition layer. However, the flux profile was not entirely smooth and the fluctuation showed temporal/spectral association with large-amplitude (δB/B ~ 0.3), low-frequency (approximately equal or less than 0.1 Ω(sub ce) where Ω(sub ce) is the cyclotron frequency), obliquely propagating (θ(sub kB) ~ 30°–60°, where θ(sub kB) is the angle between the wave vector and background magnetic field) whistler waves, indicating that the particles were interacting with the waves. Particle simulations demonstrate that, although linear cyclotron resonances with ∼0.5 keV electrons are unlikely due to the obliquity and low frequencies of the waves, the electrons are still scattered beyond 90° pitch angle by (1) resonant mirroring (transit-time damping), (2) non-resonant mirroring, and (3) subharmonic cyclotron resonances. Such coupled nonlinear scattering processes are likely to provide the stochasticity needed to explain the power-law formation.

Oka, M.↗

Analysis of Satellite-to-Satellite Tracking (SST) and altimetry data from GEOS-C

Radar altimetry and satellite-to-satellite (SST) range and range rate tracking measurements were used to infer the exterior gravitational field of the earth and the structure of the geoid from GEOS-C metric data. Under the SST analysis, a direct point-by-point estimate of gravity disturbance by means of a recursive filter with backward smoothing was attempted but had to be forsaken because of poor convergence. The adopted representation consists of a more or less uniform grid of discrete masses at a depth of approximately 400 km from the earth's surface. The layer is superimposed on a spherical harmonics model. The procedure for smoothing the altimetry and inferring the fine-structured gravity field over the Atlantic test area is described. The local disturbances are represented by means of a density layer. The altimeter height biases were first estimated by a least squares adjustment at orbital crossover points. After taking out the bias, long wavelength contributions from GEM-6 as well as a calibration correction were subtracted. The residual heights were then represented by a mass distribution beneath the earth's surface.

Wong, L.↗

Discrete Tchebycheff orthonormal polynomials and applications

Discrete Tchebycheff orthonormal polynomials offer a convenient way to make least squares polynomial fits of uniformly spaced discrete data. Computer programs to do so are simple and fast, and appear to be less affected by computer roundoff error, for the higher order fits, than conventional least squares programs. They are useful for any application of polynomial least squares fits: approximation of mathematical functions, noise analysis of radar data, and real time smoothing of noisy data, to name a few.

Lear, W. M.↗

A 5.57 hour modulation of the X-ray flux of 4U1822-37: A new model for CYG X-3

The X-ray flux of 4U1822-37 is shown to be modulated with the 5.57 hour period of its optical counterpart. The X-ray light curve is two component with a smooth sinusoidal like 25 percent semiamplitude modulation and a 30 minute dip approximately 0.2 in phase following the other minimum. The X-ray spectrum is a relatively flat power law up to 17 keV, above which its steepens. Iron emission is detected at 6.7 keV with a 4 keV FWHM and an equivalent width of 1100 eV. There is an excee below 2 keV that is consistent with either a 0.25 keV thermal component or 350 eV equivalent width iron L emission. A slight softening of the spectrum is seen during both X-ray minima. The dip is interpreted as the partial occultation of an extended cloud of optically thick highly ionized material surrounding the central X-ray source. Modeling the eclipse gives a system inclination of 70-79 deg and a spherical cloud radius of 0.2-0.3 solar radius. Models for the long term modulation are considered. The properties of this source are compared to those of Cyg X-3. It is concluded that they are similar systems.

White, N. E.↗

Wipe Melt for InP Seed Substrate

Seed substrate pulled under InGaAsP wipe melt at temperature of 650 degrees C for approximately 15 to 60 seconds. After pullthrough and meltback, surface has high smoothness and luster without meniscus lines, exaggerated erosion, pits, and pearls characteristic of other wipe melts. Layer-to-substrate interface structure is more planar and of better quality. New combination of elements leaves smooth high luster surface.

Hawrylo, F. Z.↗

Performance comparison of axisymmetric and three-dimensional hydrogen film coolant injection in a 110 N hydrogen/oxygen rocket

An experimental performance comparison of two geometrically different fuel film coolant injection sleeves was conducted on a 110 N gaseous hydrogen/oxygen rocket. One sleeve had slots milled axially down the walls and the other had a smooth surface to give axisymmetric flow. The comparison was made to investigate a conclusion in an earlier study that attributed a performance underprediciton to a simplifying modeling assumption of axisymmetric fuel film flow. The smooth sleeve had higher overall performance at one film coolant percentage and approximately the same or slightly better at another. The study showed that the lack of modeling of three-dimensional effects was not the cause of the performance underprediciton as speculated in earlier analytical studies.

Arrington, Lynn A.↗

Direct numerical simulation of instabilities in parallel flow with spherical roughness elements

Results from a direct numerical simulation of laminar flow over a flat surface with spherical roughness elements using a spectral-element method are given. The numerical simulation approximates roughness as a cellular pattern of identical spheres protruding from a smooth wall. Periodic boundary conditions on the domain's horizontal faces simulate an infinite array of roughness elements extending in the streamwise and spanwise directions, which implies the parallel-flow assumption, and results in a closed domain. A body force, designed to yield the horizontal Blasius velocity in the absence of roughness, sustains the flow. Instabilities above a critical Reynolds number reveal negligible oscillations in the recirculation regions behind each sphere and in the free stream, high-amplitude oscillations in the layer directly above the spheres, and a mean profile with an inflection point near the sphere's crest. The inflection point yields an unstable layer above the roughness (where U''(y) is less than 0) and a stable region within the roughness (where U''(y) is greater than 0). Evidently, the instability begins when the low-momentum or wake region behind an element, being the region most affected by disturbances (purely numerical in this case), goes unstable and moves. In compressible flow with periodic boundaries, this motion sends disturbances to all regions of the domain. In the unstable layer just above the inflection point, the disturbances grow while being carried downstream with a propagation speed equal to the local mean velocity; they do not grow amid the low energy region near the roughness patch. The most amplified disturbance eventually arrives at the next roughness element downstream, perturbing its wake and inducing a global response at a frequency governed by the streamwise spacing between spheres and the mean velocity of the most amplified layer.

Deanna, R. G.↗

Performance Comparison of Axisymmetric and Three-dimensional Hydrogen Film Coolant Injection in a 110N Hydrogen/oxygen Rocket

An experimental performance comparison of two geometrically different fuel film coolant injection sleeves was conducted on a 110 N gaseous hydrogen/oxygen rocket. One sleeve had slots milled axially down the walls and the other had a smooth surface to give axisymmetric flow. The comparison was made to investigate a conclusion in an earlier study that attributed a performance underprediction to a symplifying modeling assumption of axisymmetric fuel film flow. The smooth sleeve had higher overall performance at one film coolant percentage and approximately the same or slightly better at another. The study showed that the lack of modeling of three-dimensional effects was not the cause of the performance underprediction as speculated in earlier analytical studies.

Arrington, Lynn A.↗

Localization of the human 64kD autoantigen D1 to myofibrils in a subset of extraocular muscle fibers

PURPOSE. To evaluate the tissue-specific expression pattern of the 64kD human autoantigen D1, a tropomodulin-related protein that may be involved in thyroid-associated ophthalmopathy. METHODS. Recombinant 64kD human autoantigen D1 was generated in a bacterial expression system and used to immunize rabbits. Specific antibodies were affinity-purified and used for Western blots on normal and hyperthyroid rat and rabbit tissue, and immunofluorescence localization on cryosections of rat tissue. RESULTS. Anti-64kD human autoantigen D1 antibodies recognize specifically a approximately 70kD polypeptide in western blots of extraocular muscle, sternothyroid muscle, and smooth muscle. Immunofluorescence staining demonstrates that the 64kD human autoantigen D1 localizes to myofibrils in slow fibers from rat extraocular and sternothyroid muscle. The level of this protein is not altered in extraocular muscles from hyperthyroid rabbits. CONCLUSIONS. The 64kD human autoantigen D1 is expressed in slow fibers of extraocular and sternothyroid muscles as a component of myofibrils, and is not upregulated in conditions of hyperthyroidism.

Proteins/genetics/immunology/metabolism↗

An SDO/AIA-Observed Filament Eruption Triggered by a Lid-Removal Onset Mechanism

An eruption of a solar filament often presages the onset of a more general solar eruption, often leading to a solar flare and coronal mass ejection (CME). Among the mechanisms suggested for triggering eruptions are flux cancellation, flux emergence, tether‐cutting reconnection, and breakout reconnection. Here we present an example of a filament eruption due to a different trigger mechanism, which we call ``lid removal,'' whereby a magnetic structure overlying the filament is removed by a preceding adjacent eruption, rendering MHD unstable the magnetic system containing the filament and resulting in the subsequent eruption of the filament. This filament eruption occurred on 23 Jan 2013, and was well‐seen in SDO/AIA 193 Ang images. Prior to its eruption the filament was at an approximately constant height above the solar surface for approx. 4 hours, before smoothly lifting off. Evidence for the overlying "lid'' field was difficult to discern in 193 Ang images, but was apparent in hotter coronal images, such as SDO/AIA 335. Removal of the lid field was due to an eruption of that field visible in the hotter‐corona images. In this way, the lid‐removal filament‐eruption mechanism is similar to recent observations of connected or cascading eruptions originating from magnetically‐linked locations.

Sterling, Alphonse C.↗

An SDO/AIA-Observed Filament Eruption Triggered by a Lid-Removal Onset Mechanism

An eruption of a solar filament often presages the onset of a more general solar eruption, often leading to a solar flare and coronal mass ejection (CME). Among the mechanisms suggested for triggering eruptions are flux cancellation, flux emergence, tether-cutting reconnection, and breakout reconnection. Here we present an example of a filament eruption due to a different trigger mechanism, which we call "lid removal,'' whereby a magnetic structure overlying the filament is removed by a preceding adjacent eruption, rendering MHD unstable the magnetic system containing the filament and resulting in the subsequent eruption of the filament. This filament eruption occurred on 23 Jan 2013, and was well-seen in SDO/AIA 193 Ang images. Prior to its eruption the filament was at an approximately constant height above the solar surface for approx 4 hours, before smoothly lifting off. Evidence for the overlying "lid'' field was difficult to discern in 193 Ang images, but was apparent in hotter coronal images, such as SDO/AIA 335. Removal of the lid field was due to an eruption of that field visible in the hotter‐corona images. In this way, the lid‐removal filament‐eruption mechanism is similar to recent observations of connected or cascading eruptions originating from magnetically‐linked locations.

Sterling, Alphonse C.↗

First Results from NuSTAR Observations of Mkn 421

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad nonthermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some intial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant f ux variations are smooth on timescales of hours, with typical intra-hour variations of approximately less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

active galactic nuclei (AGN)↗

Landing-Impact Characteristics of Load-Alleviating Struts on a Model of a Winged Space Vehicle

A brief experimental investigation was made of the landing-impact characteristics of a 1/9-scale dynamic model of a winged space vehicle. The landing tests were made by catapulting a free model onto a hard; surface runway and onto water. The model had a conical fuselage and a flat - plate wing with a basic delta planform and 75 deg sweepback of the leading edge. The use of yielding-metal shock absorbers and various landing-gear arrangements was investigated during landing impact. The basic landing gear consisted of a dual rubber-tired nose wheel and twin main skids aft of the center of gravity near the wing tips. landing motion and acceleration data were obtained over a range of landing attitudes, gross weights, and initial sinking speeds. Brief tests were made with an alternate nose-wheel location. An all-skid configuration also was briefly evaluated for hard-surface and water landings. The landing gear employing yielding struts for impact-energy absorption during hard-surface landings resulted in accelerations of approximately 5 1/2 g near the nose gear over a range of landing parameters. Replacing the nose wheel and tire with a skid did not significantly change the accelerations. Landings in smooth water with rigid struts and adequate planing area at the nose skid resulted in a maximum landing acceleration of approximately 4g.

Blanchard, Ulysse J.↗

A multiresolution wavelet representation in two or more dimensions

In the multiresolution approximation, a signal is examined on a hierarchy of resolution scales by projection onto sets of smoothing functions. Wavelets are used to carry the detail information connecting adjacent sets in the resolution hierarchy. An algorithm has been implemented to perform a multiresolution decomposition in n greater than or equal to 2 dimensions based on wavelets generated from products of 1-D wavelets and smoothing functions. The functions are chosen so that an n-D wavelet may be associated with a single resolution scale and orientation. The algorithm enables complete reconstruction of a high resolution signal from decomposition coefficients. The signal may be oversampled to accommodate non-orthogonal wavelet systems, or to provide approximate translational invariance in the decomposition arrays.

Bromley, B. C.↗

Derivative-free stochastic optimization via adaptive sampling strategies

In this paper, we present a novel derivative-free framework for solving unconstrained stochastic optimization problems. Many problems in fields ranging from simulation optimization to reinforcement learning to quantum computing involve settings where only stochastic function values are obtained via a zeroth-order oracle, which has no available gradient information and necessitates the usage of derivative-free optimization methodologies. Our approach includes estimating gradients using stochastic function evaluations and integrating adaptive sampling techniques to control the accuracy in these stochastic approximations. Our framework encapsulates several gradient estimation techniques, including standard finite-difference, Gaussian smoothing, sphere smoothing, randomized coordinate finite-difference, and randomized subspace finite-difference methods. We provide theoretical convergence guarantees for our framework and analyze the worst-case iteration and sample complexities associated with each gradient estimation method. Finally, we demonstrate the empirical performance of the methods on logistic regression and nonlinear least squares problems.

Adaptive sampling↗